An improved positivity preserving odd degree-n Said-Ball boundary curves on rectangular grid using partial differential equation

This paper discusses the sufficient conditions for positivity preserving odd degree-n Said-Ball boundary curves defined on a rectangular grid.We derive a sufficient condition on boundary curves of rectangular Said-Ball patches where the lower bound ordinates are adjusted independently.To construct the boundary curves for each rectangular patch, the Said-Ball polynomial solution of fourth order PDE will be considered where its coefficients can be calculated using edge Said-Ball ordinates which fulfill the positivity preserving conditions.Graphical examples are presented using well-known test functions

G1 Range Restricted Data Interpolation Using Bézier Triangular Patch

Pembinaan permukaan berparameter G 1 julat terhad kepada data yang semua terletak di sebelah satu satah kekangan dipertimbangkan. Permukaan interpolasi dibina secara cebis demi cebis daripada gabungan cembung tiga tampalan segi tiga Bézier kuartik. Syarat cukup untuk keselanjaran satah tangen sepanjang sempadan dua tampalan Bézier kuartik dibentangkan. The construction of range restricted G 1 parametric surface to data that all lie on one side of a constraint plane is considered. The interpolating surface is developed piecewise as the convex combination of three quartic Bézier triangular patches. Sufficient tangent plane continuity conditions along the common boundary of two adjacent quartic Bézier triangular patches are presented

Algorithms for curve design and accurate computations with totally positive matrices

Esta tesis doctoral se enmarca dentro de la teoría de la Positividad Total. Las matrices totalmente positivas han aparecido en aplicaciones de campos tan diversos como la Teoría de la Aproximación, la Biología, la Economía, la Combinatoria, la Estadística, las Ecuaciones Diferenciales, la Mecánica, el Diseño Geométrico Asistido por Ordenador o el Álgebra Numérica Lineal. En esta tesis nos centraremos en dos de los campos que están relacionados con matrices totalmente positivas.This doctoral thesis is framed within the theory of Total Positivity. Totally positive matrices have appeared in applications from fields as diverse as Approximation Theory, Biology, Economics, Combinatorics, Statistics, Differential Equations, Mechanics, Computer Aided Geometric Design or Linear Numerical Algebra. In this thesis, we will focus on two of the fields that are related to totally positive matrices.<br /

2D and 3D surface image processing algorithms and their applications

This doctoral dissertation work aims to develop algorithms for 2D image segmentation application of solar filament disappearance detection, 3D mesh simplification, and 3D image warping in pre-surgery simulation. Filament area detection in solar images is an image segmentation problem. A thresholding and region growing combined method is proposed and applied in this application. Based on the filament area detection results, filament disappearances are reported in real time. The solar images in 1999 are processed with this proposed system and three statistical results of filaments are presented. 3D images can be obtained by passive and active range sensing. An image registration process finds the transformation between each pair of range views. To model an object, a common reference frame in which all views can be transformed must be defined. After the registration, the range views should be integrated into a non-redundant model. Optimization is necessary to obtain a complete 3D model. One single surface representation can better fit to the data. It may be further simplified for rendering, storing and transmitting efficiently, or the representation can be converted to some other formats. This work proposes an efficient algorithm for solving the mesh simplification problem, approximating an arbitrary mesh by a simplified mesh. The algorithm uses Root Mean Square distance error metric to decide the facet curvature. Two vertices of one edge and the surrounding vertices decide the average plane. The simplification results are excellent and the computation speed is fast. The algorithm is compared with six other major simplification algorithms. Image morphing is used for all methods that gradually and continuously deform a source image into a target image, while producing the in-between models. Image warping is a continuous deformation of a: graphical object. A morphing process is usually composed of warping and interpolation. This work develops a direct-manipulation-of-free-form-deformation-based method and application for pre-surgical planning. The developed user interface provides a friendly interactive tool in the plastic surgery. Nose augmentation surgery is presented as an example. Displacement vector and lattices resulting in different resolution are used to obtain various deformation results. During the deformation, the volume change of the model is also considered based on a simplified skin-muscle model

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

Modelado jerárquico de objetos 3D con superficies de subdivisión

Las SSs (Superficies de Subdivisión) son un potente paradigma de modelado de objetos 3D (tridimensionales) que establece un puente entre los dos enfoques tradicionales a la aproximación de superficies, basados en mallas poligonales y de parches alabeados, que conllevan problemas uno y otro. Los esquemas de subdivisión permiten definir una superficie suave (a tramos), como las más frecuentes en la práctica, como el límite de un proceso recursivo de refinamiento de una malla de control burda, que puede ser descrita muy compactamente. Además, la recursividad inherente a las SSs establece naturalmente una relación de anidamiento piramidal entre las mallas / NDs (Niveles de Detalle) generadas/os sucesivamente, por lo que las SSs se prestan extraordinariamente al AMRO (Análisis Multiresolución mediante Ondículas) de superficies, que tiene aplicaciones prácticas inmediatas e interesantísimas, como la codificación y la edición jerárquicas de modelos 3D. Empezamos describiendo los vínculos entre las tres áreas que han servido de base a nuestro trabajo (SSs, extracción automática de NDs y AMRO) para explicar como encajan estas tres piezas del puzzle del modelado jerárquico de objetos de 3D con SSs. El AMRO consiste en descomponer una función en una versión burda suya y un conjunto de refinamientos aditivos anidados jerárquicamente llamados "coeficientes ondiculares". La teoría clásica de ondículas estudia las señales clásicas nD: las definidas sobre dominios paramétricos homeomorfos a R" o (0,1)n como el audio (n=1), las imágenes (n=2) o el vídeo (n=3). En topologías menos triviales, como las variedades 2D) (superficies en el espacio 3D), el AMRO no es tan obvio, pero sigue siendo posible si se enfoca desde la perspectiva de las SSs. Basta con partir de una malla burda que aproxime a un bajo ND la superficie considerada, subdividirla recursivamente y, al hacerlo, ir añadiendo los coeficientes ondiculares, que son los detalles 3D necesarios para obtener aproximaciones más y más finas a la superficie original. Pasamos después a las aplicaciones prácticas que constituyen nuestros principal desarrollo original y, en particular, presentamos una técnica de codificación jerárquica de modelos 3D basada en SSs, que actúa sobre los detalles 3D mencionados: los expresa en un referencial normal loscal; los organiza según una estructura jerárquica basada en facetas; los cuantifica dedicando menos bits a sus componentes tangenciales, menos energéticas, y los "escalariza"; y los codifica dinalmente gracias a una técnica similar al SPIHT (Set Partitioning In Hierarchical Tress) de Said y Pearlman. El resultado es un código completamente embebido y al menos dos veces más compacto, para superficies mayormente suaves, que los obtenidos con técnicas de codificación progresiva de mallas 3D publicadas previamente, en las que además los NDs no están anidados piramidalmente. Finalmente, describimos varios métodos auxiliares que hemos desarrollado, mejorando técnicas previas y creando otras propias, ya que una solución completa al modelado de objetos 3D con SSs requiere resolver otros dos problemas. El primero es la extracción de una malla base (triangular, en nuestro caso) de la superficie original, habitualmente dada por una malla triangular fina con conectividad arbitraria. El segundo es la generación de un remallado recursivo con conectividad de subdivisión de la malla original/objetivo mediante un refinamiento recursivo de la malla base, calculando así los detalles 3D necesarios para corregir las posiciones predichas por la subdivisión para nuevos vértices

Global multi-disciplinary design optimisation of super-hypersonic components and vehicles

Executive summary available: pages i-ii