Helical polynomial curves and double Pythagorean hodographs II. Enumeration of low-degree curves

AbstractA “double” Pythagorean-hodograph (DPH) curve r(t) is characterized by the property that |r′(t)| and |r′(t)×r″(t)| are both polynomials in the curve parameter t. Such curves possess rational Frenet frames and curvature/torsion functions, and encompass all helical polynomial curves as special cases. As noted by Beltran and Monterde, the Hopf map representation of spatial PH curves appears better suited to the analysis of DPH curves than the quaternion form. A categorization of all DPH curve types up to degree 7 is developed using the Hopf map form, together with algorithms for their construction, and a selection of computed examples of (both helical and non-helical) DPH curves is included, to highlight their attractive features. For helical curves, a separate constructive approach proposed by Monterde, based upon the inverse stereographic projection of rational line/circle descriptions in the complex plane, is used to classify all types up to degree 7. Criteria to distinguish between the helical and non-helical DPH curves, in the context of the general construction procedures, are also discussed

An Optimization Platform for High Speed Propellers

To improve the efficiency by which current power plants translate jet energy into useful thrust the use of turboprop and in particular open rotor aircraft are being revisited. One challenge in association with developing new powerplants for such aircraft is high speed propeller design in general and noise prediction in particular. The Boxprop was invented in 2009 by GKN Aerospace in order to mitigate the effects of the tip vortex on noise and to improve upon the aerodynamics of a conventional propeller blade. The Boxprop is composed of a double-bladed propeller joined at the tips, and the design has the potential to eliminate the tip vortex, and thereby decrease that particular noise source. The complex and highly three-dimensional shape of an advanced propeller blade is challenging to model with classical propeller design methods, requiring instead more sophisticated optimization methods. This paper presents an optimization platform developed for high speed propellers, and illustrates its use by performing a reduced aerodynamic optimization of the Boxprop. The optimization process starts by performing a Latin Hypercube Sampling of the design space, and analyzes the resulting geometries using CFD. A meta-model employing radial basis functions is then used to interpolate on the obtained CFD results, which the GA uses to find optimal candidates along the obtained Pareto front. These designs are then evaluated using CFD, and their data added to the meta-model. The process iterates until the meta-model converges. The results of this paper demonstrate the capability of the presented optimization platform, and applying it on the Boxprop has resulted in valuable design improvements and insights. The obtained designs show less blade interference, more efficiently loaded blades, and less produced swirl. The methodology for geometry generation, meshing and optimizing is fast, robust, and readily extendable to other types of optimization problems, and paves the way for future collaborative research in the area of turbomachinery

BEST : Bézier-Enhanced Shell Triangle : a new rotation-free thin shell finite element

A new thin shell finite element is presented. This new element doesn’ t have rotational degrees of freedom. Instead, in order to overcome the C1 continuity requirement across elements, the author resorts to enhance the geometric description of the flat triangles of a mesh made out of linear triangles, by means of Bernstein polynomials and triangular Bernstein-Bézier patches. The author estimates the surface normals at the nodes of a mesh of triangles, in order to use them to define the Bernstein-Bézier patches. Ubach, Estruch and García-Espinosa performed a comprehensive statistical comparison of different weighting factors. The conclusion of that work is that the inverse of the area of the circumscribed circle to the triangle and the internal angle of the triangle at the node considered, should be used as weighting factor. Using this new weighting factor, we reduce by about 10% the root mean square error in the estimation of normals of randomly generated surfaces with respect to the previous best weighting factor found in the literature. The author uses the information of the normal vectors at the nodes and the triangular Bernstein-Bézier patches to build cubic Bézier triangles. These cubic Bézier triangles are surface interpolants; C1 continuous at the nodes and C0 continuous across the edges. Owing to this approach, the new element is called Bézier-enhanced shell triangle (BEST). The BEST element takes advantage of all the nodes’ connectivities in each triangle of the mesh. The computation of the normal vectors at the nodes doesn’ t depend on the number of triangles surrounding each node of the mesh. The BEST element is independent from the mesh topology. A new paradigm is presented consisting on the reconstruction of the geometry of a cubic triangular element. This geometric reconstruction exploits the properties of cubic B-spline functions (cubic Bézier triangle). This way, the author builds a conforming continuum-based shell finite element. A cubic Bézier triangle has 30 parameters (3 coordinates for each of the 10 control points). Therefore it needs to apply 30 independent conditions. 15 of these conditions are given directly by the positions of the 3 vertices of the triangle and the orientations of the normal vectors at the 3 vertices. 8 of the remaining conditions are imposed introducing energy minimization considerations. These energy minimization considerations serve also to define a well-posed element. The author defines 3 different reduced problems for the 3 different shell deformation modes: bending deformation, membrane (in-plane extension) deformation and in-plane shear (drilling rotation) deformation. The only degrees of freedom of the BEST element are the vertices’ coordinates (9 variables). The remaining 21 parameters are solved internally. In order to fix the values of these 21 internal parameters, each BEST element solves 9 systems of linear equations of rank 3. The BEST element is successfully applied to the analysis of thin shells in linear and geometrically non-linear regimes using an implicit method. The non-linearity is solved using a Total Lagrangian formulation. The author succeeds at pre-integrating through-the-thickness efficiently and accurately. The through-the-thickness integrals are evaluated just once: at the reference configuration. There are just 14 through-the-thickness scalar integrals to perform for each Gauss point. The numerical examples results show that the BEST element has the potential to achieve cubic convergence. Although they also cast doubts on the possibility of reproducing this result for a wide range of problems. For in-plane shear dominated problems, the formulation used in this thesis only achieves linear convergence. For membrane oriented tests with curvature, the convergence is quadratic. The BEST element exhibits membrane locking behavior. The author suggests exploiting further the drilling rotations kinematics in order to solve membrane locking.Se presenta un nuevo elemento finito de lámina delgada. Este nuevo elemento no usa rotaciones como grados de libertad. En su lugar, para sortear el requisito de mantener continuidad C1 entre elementos, el autor mejora la descripción geométrica de los triángulos planos de una malla de triángulos lineales, por medio de polinomios de Bernstein y particiones triangulares de Bernstein-Bézier. Para definir las particiones de Bernstein-Bézier, el autor estima las normales a la superficie en los nodos de una malla de triángulos. Ubach, Estruch y García-Espinosa hicieron una comparación estadística exhaustiva entre distintos factores de ponderación. La conclusión de dicho trabajo conduce a usar como factor de ponderación: el inverso del área de la circunferencia circunscrita al triángulo y el ángulo interno del triángulo en el nodo considerado. Con este nuevo factor de ponderación, se reduce en aproximadamente un 10% el error medio cuadrático cometido en la estimación de las normales de superficies generadas aleatoriamente, respecto del mejor factor usado previamente en la literatura. Con la información de los vectores normales en los nodos, el autor construye triángulos cúbicos de Bézier. Estos triángulos cúbicos de Bézier interpolan la superficie; con continuidad C1 en los nodos y C0 en las aristas. En virtud a este planteamiento, el nuevo elemento recibe el nombre de BEST. El elemento BEST aprovecha todas las conectividades nodales de cada triángulo de la malla. El número de triángulos que rodean cada nodo de la malla no afecta al cálculo de los vectores normales. El elemento BEST es independiente de la topología de la malla. Se propone un nuevo paradigma que consiste en reconstruir la geometría de un elemento triangular cúbico. Esta reconstrucción geométrica aprovecha las propiedades de las funciones cúbicas B-spline (triángulo cúbico de Bézier). Así, el autor crea un elemento de lámina conforme basado en el continuo. Un triángulo cúbico de Bézier tiene 30 parámetros (3 coordenadas para cada uno de los 10 puntos de control). Es necesario aplicar 30 condiciones independientes. 15 de estas condiciones se deducen de la posición de los 3 vértices del triángulo y de los vectores normales en los 3 vértices. De las otras 15 condiciones, 8 se obtienen a partir de criterios de minimización de la energía. Estos criterios de minimización de la energía sirven para definir un elemento bien planteado. El autor desarrolla 3 problemas reducidos para los 3 modos de deformación de la lámina: deformación de flexión, de membrana (extensión en el plano) y de cortante en el plano (rotación de taladro). Los únicos grados de libertad del elemento BEST son las posiciones de los vértices (9 variables). Los otros 21 parámetros se resuelven internamente. Para obtener estos 21 parámetros internos, hay que resolver 9 sistemas de ecuaciones lineales de rango 3 para cada elemento BEST. Se ha aplicado el elemento BEST con éxito al cálculo de láminas delgadas en régimen lineal y geométricamente no-lineal con un método implícito. La no-linealidad se plantea con una formulación Lagrangiana total. Se demuestra cómo pre-integrar en el espesor de manera eficiente y precisa. Solo es preciso evaluar las integrales en el espesor una vez: en la configuración de referencia. Solo hay 14 integrales escalares en el espesor para cada punto de Gauss. Los ejemplos numéricos muestran que el elemento BEST tiene potencial para converger cúbicamente. Pero también existen dudas sobre la capacidad de reproducir de manera consistente este resultado en un amplio rango de problemas. En problemas dominados por la deformación de cortante en el plano, la formulación utilizada en esta tesis solo alcanza convergencia lineal. En ejemplos orientados a la deformación de membrana que incluyen curvatura, la convergencia es cuadrática. El elemento BEST sufre de bloqueo por membrana. El autor sugiere desarrollar más profundamente la cinemática de las rotaciones de taladro para resolver el bloqueo por membrana.Es presenta un nou element finit de làmina prima. Aquest nou element no fa servir rotacions com a graus de llibertat. Enlloc d'això, per esquivar el requisit de mantenir continuïtat C1 entre els elements, l'autor millora la descripció geomètrica dels triangles plans d'una malla de triangles lineals, mitjançant polinomis de Bernstein i particions triangulars de Bernstein-Bézier.Per definir les particions de Bernstein-Bézier, l'autor estima les normals a la superfície en els nodes d'una malla de triangles. Ubach, Estruch i García-Espinosa varen fer una comparació estadística exhaustiva entre diferents factors de ponderació. La conclusió d'aquest treball condueix a fer servir com a factor de ponderació: l'invers de l'àrea de la circumferència circumscrita al triangle i l'angle intern del triangle en el node considerat. Amb aquest nou factor de ponderació, es redueix aproximadament en un 10% l'error quadràtic mig comès en l'estimació de les normals de superfícies generades aleatòriament, respecte del millor factor usat prèviament a la literatura.Amb la informació dels vectors normals en els nodes, l'autor construeix triangles cúbics de Bézier. Aquests triangles cúbics de Bézier interpolen la superfície; amb continuïtat C1 als nodes i C0 a les arestes. En virtut d'aquest plantejament, el nou element rep el nom de BEST (Bézier-enhanced shell triangle).L'element BEST aprofita totes les connectivitats nodals de cada triangle de la malla. El nombre de triangles que envolten cada node de la malla no afecta al càlcul dels vectors normals. L'element BEST és independent de la topologia de la malla.Es proposa un nou paradigma que consisteix en reconstruir la geometria d'un element triangular cúbic. Aquesta reconstrucció geomètrica aprofita les propietats de les funcions cúbiques B-spline (triangle cúbic de Bézier). D'aquesta manera l'autor crea un element de làmina que és conforme i basat en el continu.Un triangle cúbic de Bézier té 30 paràmetres (3 coordenades per cadascun dels 10 punts de control). Cal aplicar 30 condicions independents. 15 d'aquestes condicions es dedueixen de la posició dels 3 vèrtexs del triangle i dels vectors normals en els 3 vèrtexs.De les 15 condicions restants, 8 s'obtenen a partir de criteris de minimització de l'energia. Aquests criteris de minimització de l'energia serveixen per definir un element ben plantejat. L'autor desenvolupa 3 problemes reduïts per als 3 modes de deformació de la làmina: deformació de flexió, de membrana (extensió en el pla) i de tallant en el pla (rotació de barrina).Els únics graus de llibertat de l'element BEST són les posicions dels vèrtexs (9 variables). Els altres 21 paràmetres es resolen internament. Per obtenir aquests 21 paràmetres interns, cal resoldre 9 sistemes d'equacions lineals de rang 3 per cada element BEST.S'ha aplicat l'element BEST amb èxit al càlcul de làmines primes en règim lineal i geomètricament no-lineal fent servir un mètode implícit. La no-linealitat es planteja amb una formulació Lagrangiana total. Es demostra com es pot pre-integrar a través del gruix de manera eficient i precisa. Només cal avaluar les integrals a través del gruix un cop: a la configuració de referència. Només hi ha 14 integrals escalars a través del gruix per a cada punt de Gauss. Els exemples numèrics mostren que l'element BEST té potencial per convergir cúbicament. Però també hi ha dubtes de que aquest resultat es pugui reproduir de manera consistent per un ventall ampli de problemes. En problemes dominats per la deformació de tallant en el pla, la formulació emprada en aquesta tesi només assoleix convergència lineal. En exemples orientats a la deformació de membrana que incloguin curvatura, la convergència és quadràtica. L'element BEST pateix de bloqueig per membrana. L'autor suggereix desenvolupar en més profunditat la cinemàtica de les rotacions de barrina per resoldre el bloqueig per membrana

New Techniques for the Modeling, Processing and Visualization of Surfaces and Volumes

With the advent of powerful 3D acquisition technology, there is a growing demand for the modeling, processing, and visualization of surfaces and volumes. The proposed methods must be efficient and robust, and they must be able to extract the essential structure of the data and to easily and quickly convey the most significant information to a human observer. Independent of the specific nature of the data, the following fundamental problems can be identified: shape reconstruction from discrete samples, data analysis, and data compression. This thesis presents several novel solutions to these problems for surfaces (Part I) and volumes (Part II). For surfaces, we adopt the well-known triangle mesh representation and develop new algorithms for discrete curvature estimation,detection of feature lines, and line-art rendering (Chapter 3), for connectivity encoding (Chapter 4), and for topology preserving compression of 2D vector fields (Chapter 5). For volumes, that are often given as discrete samples, we base our approach for reconstruction and visualization on the use of new trivariate spline spaces on a certain tetrahedral partition. We study the properties of the new spline spaces (Chapter 7) and present efficient algorithms for reconstruction and visualization by iso-surface rendering for both, regularly (Chapter 8) and irregularly (Chapter 9) distributed data samples

Generative Mesh Modeling

Generative Modeling is an alternative approach for the description of three-dimensional shape. The basic idea is to represent a model not as usual by an agglomeration of geometric primitives (triangles, point clouds, NURBS patches), but by functions. The paradigm change from objects to operations allows for a procedural representation of procedural shapes, such as most man-made objects. Instead of storing only the result of a 3D construction, the construction process itself is stored in a model file. The generative approach opens truly new perspectives in many ways, among others also for 3D knowledge management. It permits for instance to resort to a repository of already solved modeling problems, in order to re-use this knowledge also in different, slightly varied situations. The construction knowledge can be collected in digital libraries containing domain-specific parametric modeling tools. A concrete realization of this approach is a new general description language for 3D models, the "Generative Modeling Language" GML. As a Turing-complete "shape programming language" it is a basis of existing, primitv based 3D model formats. Together with its Runtime engine the GML permits - to store highly complex 3D models in a compact form, - to evaluate the description within fractions of a second, - to adaptively tesselate and to interactively display the model, - and even to change the models high-level parameters at runtime.Die generative Modellierung ist ein alternativer Ansatz zur Beschreibung von dreidimensionaler Form. Zugrunde liegt die Idee, ein Modell nicht wie üblich durch eine Ansammlung geometrischer Primitive (Dreiecke, Punkte, NURBS-Patches) zu beschreiben, sondern durch Funktionen. Der Paradigmenwechsel von Objekten zu Geometrie-erzeugenden Operationen ermöglicht es, prozedurale Modelle auch prozedural zu repräsentieren. Statt das Resultat eines 3D-Konstruktionsprozesses zu speichern, kann so der Konstruktionsprozess selber repräsentiert werden. Der generative Ansatz eröffnet unter anderem gänzlich neue Perspektiven für das Wissensmanagement im 3D-Bereich. Er ermöglicht etwa, auf einen Fundus bereits gelöster Konstruktions-Aufgaben zurückzugreifen, um sie in ähnlichen, aber leicht variierten Situationen wiederverwenden zu können. Das Konstruktions-Wissen kann dazu in Form von Bibliotheken parametrisierter, Domänen-spezifischer Modellier-Werkzeuge gesammelt werden. Konkret wird dazu eine neue allgemeine Modell-Beschreibungs-Sprache vorgeschlagen, die "Generative Modeling Language" GML. Als Turing-mächtige "Programmiersprache für Form" stellt sie eine echte Verallgemeinerung existierender Primitiv-basierter 3D-Modellformate dar. Zusammen mit ihrer Runtime-Engine erlaubt die GML, - hochkomplexe 3D-Objekte extrem kompakt zu beschreiben, - die Beschreibung innerhalb von Sekundenbruchteilen auszuwerten, - das Modell adaptiv darzustellen und interaktiv zu betrachten, - und die Modell-Parameter interaktiv zu verändern