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 /

Isogeometric treatment of large deformation contact and debonding problems with T-splines: a review

AbstractWithin a setting where the isogeometric analysis (IGA) has been successful at bringing two different research fields together, i.e. Computer Aided Design (CAD) and numerical analysis, T-spline IGA is applied in this work to frictionless contact and mode-I debonding problems between deformable bodies in the context of large deformations. Based on the concept of IGA, the smooth basis functions are adopted to describe surface geometries and approximate the numerical solutions, leading to higher accuracy in the contact integral evaluation. The isogeometric discretizations are here incorporated into an existing finite element framework by using Bézier extraction, i.e. a linear operator which maps the Bernstein polynomial basis on Bézier elements to the global isogeometric basis. A recently released commercial T-spline plugin for Rhino is herein used to build the analysis models adopted in this study.In such context, the continuum is discretized with cubic T-splines, as well as with Non Uniform Rational B-Splines (NURBS) and Lagrange polynomial elements for comparison purposes, and a Gauss-point-to-surface (GPTS) formulation is combined with the penalty method to treat the contact constraints. The purely geometric enforcement of the non-penetration condition in compression is generalized to encompass both contact and mode-I debonding of interfaces which is approached by means of cohesive zone (CZ) modeling, as commonly done by the scientific community to analyse the progressive damage of materials and interfaces. Based on these models, non-linear relationships between tractions and relative displacements are assumed. These relationships dictate both the work of separation per unit fracture surface and the peak stress that has to be reached for the crack formation. In the generalized GPTS formulation an automatic switching procedure is used to choose between cohesive and contact models, depending on the contact status. Some numerical results are first presented and compared in 2D for varying resolutions of the contact and/or cohesive zone, including frictionless sliding and cohesive debonding, all featuring the competitive accuracy and performance of T-spline IGA. The superior accuracy of T-spline interpolations with respect to NURBS and Lagrange interpolations for a given number of degrees of freedom (Dofs) is always verified. The isogeometric formulation is also extended to 3D bodies, where some examples in large deformations based on T-spline discretizations show an high smoothness of the reaction history curves

Subdivision and manifold techniques for isogeometric design and analysis of surfaces

Design of surfaces and analysis of partial differential equations defined on them are of great importance in engineering applications, e.g., structural engineering, automotive and aerospace. This thesis focuses on isogeometric design and analysis of surfaces, which aims to integrate engineering design and analysis by using the same representation for both. The unresolved challenge is to develop a desirable surface representation that simultaneously satisfies certain favourable properties on meshes of arbitrary topology around the extraordinary vertices (EVs), i.e., vertices not shared by four quadrilaterals or three triangles. These properties include high continuity (geometric or parametric), optimal convergence in finite element analysis as well as simplicity in terms of implementation. To overcome the challenge, we further develop subdivision and manifold surface modelling techniques, and explore a possible scheme to combine the distinct appealing properties of the two. The unique advantages of the developed techniques have been confirmed with numerical experiments. Subdivision surfaces generate smooth surfaces from coarse control meshes of arbitrary topology by recursive refinement. Around EVs the optimal refinement weights are application-dependent. We first review subdivision-based finite elements. We then proceed to derive the optimal subdivision weights that minimise finite element errors and can be easily incorporated into existing implementations of subdivision schemes to achieve the same accuracy with much coarser meshes in engineering computations. To this end, the eigenstructure of the subdivision matrix is extensively used and a novel local shape decomposition approach is proposed to choose the optimal weights for each EV independently. Manifold-based basis functions are derived by combining differential-geometric manifold techniques with conformal parametrisations and the partition of unity method. This thesis derives novel manifold-based basis functions with arbitrary prescribed smoothness using quasi-conformal maps, enabling us to model and analyse surfaces with sharp features, such as creases and corners. Their practical utility in finite element simulation of hinged or rigidly joined structures is demonstrated with Kirchhoff-Love thin shell examples. We also propose a particular manifold basis reproducing subdivision surfaces away from EVs, i.e., B-splines, providing a way to combine the appealing properties of subdivision (available in industrial software) for design and manifold basis (relatively new) for analysis.Cambridge International Scholarship Scheme (CISS) by Cambridge Trus