106 research outputs found

    Curves with rational chord-length parametrization

    Get PDF
    It has been recently proved that rational quadratic circles in standard Bezier form are parameterized by chord-length. If we consider that standard circles coincide with the isoparametric curves in a system of bipolar coordinates, this property comes as a straightforward consequence. General curves with chord-length parametrization are simply the analogue in bipolar coordinates of nonparametric curves. This interpretation furnishes a compact explicit expression for all planar curves with rational chord-length parametrization. In addition to straight lines and circles in standard form, they include remarkable curves, such as the equilateral hyperbola, Lemniscate of Bernoulli and Limacon of Pascal. The extension to 3D rational curves is also tackled

    A novel parallel algorithm for surface editing and its FPGA implementation

    Get PDF
    A thesis submitted to the University of Bedfordshire in partial fulfilment of the requirements for the degree of Doctor of PhilosophySurface modelling and editing is one of important subjects in computer graphics. Decades of research in computer graphics has been carried out on both low-level, hardware-related algorithms and high-level, abstract software. Success of computer graphics has been seen in many application areas, such as multimedia, visualisation, virtual reality and the Internet. However, the hardware realisation of OpenGL architecture based on FPGA (field programmable gate array) is beyond the scope of most of computer graphics researches. It is an uncultivated research area where the OpenGL pipeline, from hardware through the whole embedded system (ES) up to applications, is implemented in an FPGA chip. This research proposes a hybrid approach to investigating both software and hardware methods. It aims at bridging the gap between methods of software and hardware, and enhancing the overall performance for computer graphics. It consists of four parts, the construction of an FPGA-based ES, Mesa-OpenGL implementation for FPGA-based ESs, parallel processing, and a novel algorithm for surface modelling and editing. The FPGA-based ES is built up. In addition to the Nios II soft processor and DDR SDRAM memory, it consists of the LCD display device, frame buffers, video pipeline, and algorithm-specified module to support the graphics processing. Since there is no implementation of OpenGL ES available for FPGA-based ESs, a specific OpenGL implementation based on Mesa is carried out. Because of the limited FPGA resources, the implementation adopts the fixed-point arithmetic, which can offer faster computing and lower storage than the floating point arithmetic, and the accuracy satisfying the needs of 3D rendering. Moreover, the implementation includes Bézier-spline curve and surface algorithms to support surface modelling and editing. The pipelined parallelism and co-processors are used to accelerate graphics processing in this research. These two parallelism methods extend the traditional computation parallelism in fine-grained parallel tasks in the FPGA-base ESs. The novel algorithm for surface modelling and editing, called Progressive and Mixing Algorithm (PAMA), is proposed and implemented on FPGA-based ES’s. Compared with two main surface editing methods, subdivision and deformation, the PAMA can eliminate the large storage requirement and computing cost of intermediated processes. With four independent shape parameters, the PAMA can be used to model and edit freely the shape of an open or closed surface that keeps globally the zero-order geometric continuity. The PAMA can be applied independently not only FPGA-based ESs but also other platforms. With the parallel processing, small size, and low costs of computing, storage and power, the FPGA-based ES provides an effective hybrid solution to surface modelling and editing

    Gn blending multiple surfaces in polar coordinates

    Get PDF
    International audienceThis paper proposes a method of Gn blending multiple parametric surfaces in polar coordinates. It models the geometric continuity conditions of parametric surfaces in polar coordinates and presents a mechanism of converting a Cartesian parametric surface into its polar coordinate form. The basic idea is first to reparameterize the parametric blendees into the form of polar coordinates. Then they are blended simultaneously by a basis function in the complex domain. To extend its compatibility, we also propose a method of converting polar coordinate blending surface into N NURBS patches. One application of this technique is to fill N-sided holes. Examples are presented to show its feasibility and practicability

    Robust multigrid methods for Isogeometric discretizations applied to poroelasticity problems

    Get PDF
    El análisis isogeométrico (IGA) elimina la barrera existente entre elementos finitos (FEA) y el diseño geométrico asistido por ordenador (CAD). Debido a esto, IGA es un método novedoso que está recibiendo una creciente atención en la literatura y recientemente se ha convertido en tendencia. Muchos esfuerzos están siendo puestos en el diseño de solvers eficientes y robustos para este tipo de discretizaciones. Dada la optimalidad de los métodos multimalla para elementos finitos, la aplicación de estosmétodos a discretizaciones isogeométricas no ha pasado desapercibida. Nosotros pensamos firmemente que los métodos multimalla son unos candidatos muy prometedores a ser solvers eficientes y robustos para IGA y por lo tanto en esta tesis apostamos por su aplicación. Para contar con un análisis teórico para el diseño de nuestros métodos multimalla, el análisis local de Fourier es propuesto como principal análisis cuantitativo. En esta tesis, a parte de considerar varios problemas escalares, prestamos especial atención al problema de poroelasticidad, concretamente al modelo cuasiestático de Biot para el proceso de consolidación del suelo. Actualmente, el diseño de métodos multimalla robustos para problemas poroelásticos respecto a parámetros físicos o el tamaño de la malla es un gran reto. Por ello, la principal contribución de esta tesis es la propuesta de métodos multimalla robustos para discretizaciones isogeométricas aplicadas al problema de poroelasticidad.La primera parte de esta tesis se centra en la construcción paramétrica de curvas y superficies dado que estas técnicas son la base de IGA. Así, la definición de los polinomios de Bernstein y curvas de Bézier se presenta como punto de partida. Después, introducimos los llamados B-splines y B-splines racionales no uniformes (NURBS) puesto que éstas serán las funciones base consideradas en nuestro estudio.La segunda parte trata sobre el análisis isogeométrico propiamente dicho. En esta parte, el método isoparamétrico es explicado al lector y se presenta el análisis isogeométrico de algunos problemas. Además, introducimos la formulación fuerte y débil de los problemas anteriores mediante el método de Galerkin y los espacios de aproximación isogeométricos. El siguiente punto de esta tesis se centra en los métodos multimalla. Se tratan las bases de los métodos multimalla y, además de introducir algunos métodos iterativos clásicos como suavizadores, también se introducen suavizadores por bloques como los métodos de Schwarz multiplicativos y aditivos. Llegados a esta parte, nos centramos en el LFA para el diseño de métodos multimalla robustos y eficientes. Además, se explican en detalle el análisis estándar y el análisis basado en ventanas junto al análisis de suavizadores por bloques y el análisis para sistemas de ecuaciones en derivadas parciales.Tras introducir las discretizaciones isogeométricas, los métodos multimalla y el LFA como análisis teórico, nuestro propósito es diseñar métodos multimalla eficientes y robustos respecto al grado polinomial de los splines para discretizaciones isogeométricas de algunos problemas escalares. Así, mostramos que el uso de métodos multimalla basados en suavizadores de tipo Schwarz multiplicativo o aditivo produce buenos resultados y factores de convergencia asintóticos robustos. La última parte de esta tesis está dedicada al análisis isogeométrico del problema de poroelasticidad. Para esta tarea, se introducen el modelo de Biot y su discretización isogeométrica. Además, presentamos una novedosa estabilización de masa para la formulación de dos campos de las ecuaciones de Biot que elimina todas las oscilaciones no físicas en la aproximación numérica de la presión. Después, nos centramos en dos tipos de solvers para estas ecuaciones poroelásticas: Solvers desacoplados y solvers monolíticos. En el primer grupo, le dedicamos una especial atención al método fixed-stress y a un método iterativo propuesto por nosotros que puede ser aplicado de forma automática a partir de la estabilización de masa ya mencionada.Por otro lado, realizamos un análisis de von Neumann para este método iterativo aplicado al problema de Terzaghi y demostramos su estabilidad y convergencia para los pares de elementos Q1 Q1, Q2 Q1 y Q3 Q2 (con suavidad global C1). Respecto al grupo de solvers monolíticos, nosotros proponemos métodos multimalla basados en suavizadores acoplados y desacoplados. En esta parte, métodosIsogeometric analysis (IGA) eliminates the gap between finite element analysis (FEA) and computer aided design (CAD). Due to this, IGA is an innovative approach that is receiving an increasing attention in the literature and it has recently become a trending topic. Many research efforts are being devoted to the design of efficient and robust solvers for this type of discretization. Given the optimality of multigrid methods for FEA, the application of these methods to IGA discretizations has not been unnoticed. We firmly think that they are a very promising approach as efficient and robust solvers for IGA and therefore in this thesis we are concerned about their application. In order to give a theoretical support to the design of multigrid solvers, local Fourier analysis (LFA) is proposed as the main quantitative analysis. Although different scalar problems are also considered along this thesis, we make a special focus on poroelasticity problems. More concretely, we focus on the quasi-static Biot's equations for the soil consolidation process. Nowadays, it is a very challenging task to achieve robust multigrid solvers for poroelasticity problems with respect physical parameters and/or the mesh size. Thus, the main contribution of this thesis is to propose robust multigrid methods for isogeometric discretizations applied to poroelasticity problems. The first part of this thesis is devoted to the introduction of the parametric construction of curves and surfaces since these techniques are the basis of IGA. Hence, with the definition of Bernstein polynomials and B\'ezier curves as a starting point, we introduce B-splines and non-uniform rational B-splines (NURBS) since these will be the basis functions considered for our numerical experiments. The second part deals with the isogeometric analysis. In this part, the isoparametric approach is explained to the reader and the isogeometric analysis of some scalar problems is presented. Hence, the strong and weak formulations by means of Galerkin's method are introduced and the isogeometric approximation spaces as well. The next point of this thesis consists of multigrid methods. The basics of multigrid methods are explained and, besides the presentation of some classical iterative methods as smoothers, block-wise smoothers such as multiplicative and additive Schwarz methods are also introduced. At this point, we introduce LFA for the design of efficient and robust multigrid methods. Furthermore, both standard and infinite subgrids local Fourier analysis are explained in detail together with the analysis for block-wise smoothers and the analysis for systems of partial differential equations. After the introduction of isogeometric discretizations, multigrid methods as our choice of solvers and LFA as theoretical analysis, our goal is to design efficient and robust multigrid methods with respect to the spline degree for IGA discretizations of some scalar problems. Hence, we show that the use of multigrid methods based on multiplicative or additive Schwarz methods provide a good performance and robust asymptotic convergence rates. The last part of this thesis is devoted to the isogeometric analysis of poroelasticity. For this task, Biot's model and its isogeometric discretization are introduced. Moreover, we present an innovative mass stabilization of the two-field formulation of Biot's equations that eliminates all the spurious oscillations in the numerical approximation of the pressure. Then, we deal with two types of solvers for these poroelastic equations: Decoupled and monolithic solvers. In the first group we devote special attention to the fixed-stress split method and a mass stabilized iterative scheme proposed by us that can be automatically applied from the mass stabilization formulation mentioned before. In addition, we perform a von Neumann analysis for this iterative decoupled solver applied to Terzaghi's problem and demonstrate that it is stable and convergent for pairs Q1-Q1, Q2-Q1 and Q3-Q2 (with global smoothness C1). Regarding the group of monolithic solvers, we propose multigrid methods based on coupled and decoupled smoothers. Coupled additive Schwarz methods are proposed as coupled smoothers for isogeometric Taylor-Hood elements. More concretely, we propose a 51-point additive Schwarz method for the pair Q2-Q1. In the last part, we also propose to use an inexact version of the fixed-stress split algorithm as decoupled smoother by applying iterations of different additive Schwarz methods for each variable. For the latter approach, we consider the pairs of elements Q2-Q1 and Q3-Q2 (with global smoothness C1). Finally, thanks to LFA we manage to design efficient and robust multigrid solvers for the Biot's equations and some numerical results are shown.<br /

    Desenvolvimento de ferramentas computacionais para Análise Isogeométrica (IGA)

    Get PDF
    Mestrado em Engenharia MecânicaA few years ago drawings were made in the drawing boards and using pencils on vellum. There were no computers helping the designers in the parts modeling. After designing the object, the design was passed to the analysts. The designers and analysts were in constant communication. Nowadays, the designers used Computer Aided Design (CAD) tools in the parts modeling. For application the analysis at the geometries, initially a mesh to approximate the geometries is generated. After this, on the mesh the Finite Element Method (FEM) was applied. In complex engineering design, the generation and manipulation of meshes in FEA was estimated to take over 80% of the overall analysis time. The form to break down the barriers between engineering design and the analysis is with reconstruction the entire process, but at the same time maintaining compatibility with existing practices. Create only one geometric model is the focus of reconstruction process. This geometric model is used in the representation of the geometry, as well as in the analysis, and this concept is designated by Isogeometric Analysis (IGA). In this present work the development of the tools for generate the CAD and calculate the basis function for representation the object are proposed. Initially, the mathematical formulations for Bézier, B-Spline and NURBS, for curves and surfaces are presented. The algorithms developed to generate the curves and surfaces are demonstrated. The IGA and FEM formulation for tridimensional and bidimensional spaces are introduced. In this work, a development of a tools for application this method are proposed. The convergence of the results for FEM and IGA programs are studied and compared to the theoretical values and Abaqus comercial program. The results obtained with IGA formulation converge to the reference values.Há alguns anos atrás, os objectos eram feitos pelos designers e a criação do desenho era feita com lápis e papel vegetal. Não existiam computadores nos gabinetes de desenho para ajudar na modelação dos objectos. Após o desenho estar concluído este era entregue aos analistas para calcularem a resistência do mesmos quando solicitados por cargas externas. Assim, o gabinete de design e o gabinete de análise estavam em constante comunicação. Nos tempos de hoje os designers utilizam as ferramentas de Computer-Aided Design (CAD) para gerar os objectos, representando assim a geometria original. Por outro lado, os analistas fazem a análise baseada no Método dos Elementos Finitos (MEF). Neste método, inicialmente, gera-se uma malha para fazer a aproximação do objecto e utiliza-se esta malha gerada na análise. A forma de combater esta barreira é a construção de um novo processo de análise, mas ao mesmo tempo manter a compatibilidade com a análise do Método de Elementos Finitos. Este novo método foca-se na geração de um modelo geométrico, sendo este modelo utilizado tanto para a representação da geometria como para a análise. A principal sustentação deste novo método é a utilização das funções de base da criação e representação dos objectos, posteriormente, utilizadas na análise dos mesmos. Este novo conceito é designado por Análise Isogeométrica. Neste trabalho é exposto o desenvolvimento de ferramentas para gerar curvas e superfícies utilizando as formulações de Bézier, B-spline e NURBS. Assim, desenvolveram-se sub-rotinas para calcular as funções de base. Inicialmente apresentaram-se as formulações matemáticas e posteriormente os algoritmos desenvolvidos para a representação das curvas e superfícies. O desenvolvimento de ferramentas de análise para problemas no espaço bidimensional e tridimensional utilizando o Método de Elementos Finitos e a Análise Isogeométrica também é abordado neste trabalho. Para ser mais fácil a sua aplicação, foi desenvolvida um interface. Por fim utilizaram-se problemas e estudaram-se as curvas de convergência dos resultados e compararando-os com as referência analíticas e com o programa Abaqus. Em termos de conclusão, os resultados obtidos com a Análise Isogeométrica convergem mais rapidamente para os valores de referência do que o Abaqus e o programa desenvolvido com base no método de elementos finitos

    Görbék és felületek a geometriai modellezésben = Curves and surfaces in geometric modelling

    Get PDF
    B-spline görbék/felületek pontjai által, az alakzat két csomóértékének szimmetrikus változtatásakor leírt pályagörbéket vizsgáltuk, és olyan alakmódosítási eljárást adtunk, amivel a felület adott pontját/paramétervonalát előre megadott helyre mozgathatjuk a csomóértékek változtatásával. A C-Bézier, C-B-spline és F-B-spline görbék pályagörbéinek geometriai tulajdonságait írtuk le, és erre alapozva geometriai kényszereket kielégítő alakmódosításokat vizsgáltuk. Olyan általános leírási módot (linear blending) adtunk, mely egységesen kezeli az alakparaméterekkel rendelkező görbék széles osztályát, továbbá konkrét esetekben e paraméterek geometriai hatását írtuk le és kényszeres alakmódosításokra adtunk megoldást. A csomóértékeknek az interpoláló görbére gyakorolt hatását vizsgáltuk, mely alapján a harmadfokú interpoláció esetére interaktív alakmódosító eljárást dolgoztunk ki. Kontrollpontokkal adott görbék szingularitásainak detektálására a kontrollpontok helyzetén alapuló megoldást adtunk. Kontrollpont alapú szükséges és elégséges feltételt adtunk arra, hogy a Bézier-felület paramétervonalai egyenesek legyenek. Olyan Monte Carlo módszert dolgoztunk ki, amely rendezetlen ponthalmaz felülettel való interpolálásához négyszöghálót hoz létre a pontfelhő (mely elágazásokat és hurkokat is tartalmazhat) és annak topológikus gráfja ismeretében. A csonkolt Fourier-sorok terében olyan ciklikus bázist adtunk meg, amellyel végtelen simaságú zárt görbéket és felületeket írhatunk le. | We studied paths of points of B-spline curves/surfaces obtained by the symmetric alteration of two knot values and provided a constrained shape modification method that is capable of moving a point/isoparametric line of the surface to a user specified position. We described the geometric properties of paths of C-Bézier, C-B-spline and F-B-spline curves and on this basis we studied shape modifications subject to geometric constraints. We developed the general linear blending method that treats a wide class of curves with shape parameters in a uniform way; in special cases we described the geometric effects of shape parameters and provided constrained shape modification methods. We examined the impact of knots on the shape of interpolating curves, based on which we developed an interactive shape modification method for cubic interpolation. We proposed a control point based solution to the problem of singularity detection of curves described by control points. We provided control point based necessary and sufficient conditions for Bézier surfaces to have linear isoparametric lines. We developed a Monte Carlo method to generate a quadrilateral mesh (for surface interpolation) from point clouds (with possible junctions and loops) and their topological graph. We specified a cyclic basis in the space of truncated Fourier series by means of which we can describe closed curves and surfaces with C^infinity

    Comparison between specialized quadrature rules for method of moments with NURBS modelling applied to periodic multilayer structures

    Get PDF
    A comparison between Ma-Rokhlin-Wandzura (MRW) and double exponential (DE) quadrature rules for numerical integration of method of moments (MoM) matrix entries with singular behavior is presented for multilayer periodic structures. Non Uniform Rational B-Splines (NURBS) modelling of the layout surfaces is implemented to provide high-order description of the geometry. The comparison is carried out in order to show that quadrature rule is more suitable for MoM matrix computation in terms of sampling, accuracy of computation of MoM matrix, and CPU time consumption. The comparison of CPU time consumption shows that the numerical integration with MRW samples is roughly 15 times faster than that numerical integration using DE samples for results with similar accuracies. These promising results encourage to carry out a comparison with results obtained in previous works where a specialized approach for the specific analysis of split rings geometries was carried out. This previous approach uses spectral MoM version with specific entire domain basis function with edge singularities defined on split ring geometry. Thus, the previous approach provides accurate results with low CPU time consumption to be compared. The comparison shows that CPU time consumption obtained by MRW samples is similar to the CPU time consumption required by the previous work of specific analysis of split rings geometries. The fact that similar CPU time consumptions are obtained by MRW quadrature rules for modelling of general planar geometries and by the specialized approach for split ring geometry provides an assessment for the usage of the MRW quadrature rules and NURBS modelling. This fact provides an efficient tool for analysis of reflectarray elements with general planar layout geometries, which is suitable for reflectarray designs under local periodicity assumption where a huge number of periodic multilayer structures have to be analyzed.Agencia Estatal de InvestigaciónJunta de Comunidades de Castilla-La Manch
    corecore