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 /

Geometric Data Analysis: Advancements of the Statistical Methodology and Applications

Data analysis has become fundamental to our society and comes in multiple facets and approaches. Nevertheless, in research and applications, the focus was primarily on data from Euclidean vector spaces. Consequently, the majority of methods that are applied today are not suited for more general data types. Driven by needs from fields like image processing, (medical) shape analysis, and network analysis, more and more attention has recently been given to data from non-Euclidean spaces–particularly (curved) manifolds. It has led to the field of geometric data analysis whose methods explicitly take the structure (for example, the topology and geometry) of the underlying space into account. This thesis contributes to the methodology of geometric data analysis by generalizing several fundamental notions from multivariate statistics to manifolds. We thereby focus on two different viewpoints. First, we use Riemannian structures to derive a novel regression scheme for general manifolds that relies on splines of generalized Bézier curves. It can accurately model non-geodesic relationships, for example, time-dependent trends with saturation effects or cyclic trends. Since Bézier curves can be evaluated with the constructive de Casteljau algorithm, working with data from manifolds of high dimensions (for example, a hundred thousand or more) is feasible. Relying on the regression, we further develop a hierarchical statistical model for an adequate analysis of longitudinal data in manifolds, and a method to control for confounding variables. We secondly focus on data that is not only manifold- but even Lie group-valued, which is frequently the case in applications. We can only achieve this by endowing the group with an affine connection structure that is generally not Riemannian. Utilizing it, we derive generalizations of several well-known dissimilarity measures between data distributions that can be used for various tasks, including hypothesis testing. Invariance under data translations is proven, and a connection to continuous distributions is given for one measure. A further central contribution of this thesis is that it shows use cases for all notions in real-world applications, particularly in problems from shape analysis in medical imaging and archaeology. We can replicate or further quantify several known findings for shape changes of the femur and the right hippocampus under osteoarthritis and Alzheimer's, respectively. Furthermore, in an archaeological application, we obtain new insights into the construction principles of ancient sundials. Last but not least, we use the geometric structure underlying human brain connectomes to predict cognitive scores. Utilizing a sample selection procedure, we obtain state-of-the-art results

Representing archaeological uncertainty in cultural informatics

This thesis sets out to explore, describe, quantify, and visualise uncertainty in a cultural informatics context, with a focus on archaeological reconstructions. For quite some time, archaeologists and heritage experts have been criticising the often toorealistic appearance of three-dimensional reconstructions. They have been highlighting one of the unique features of archaeology: the information we have on our heritage will always be incomplete. This incompleteness should be reflected in digitised reconstructions of the past. This criticism is the driving force behind this thesis. The research examines archaeological theory and inferential process and provides insight into computer visualisation. It describes how these two areas, of archaeology and computer graphics, have formed a useful, but often tumultuous, relationship through the years. By examining the uncertainty background of disciplines such as GIS, medicine, and law, the thesis postulates that archaeological visualisation, in order to mature, must move towards archaeological knowledge visualisation. Three sequential areas are proposed through this thesis for the initial exploration of archaeological uncertainty: identification, quantification and modelling. The main contributions of the thesis lie in those three areas. Firstly, through the innovative design, distribution, and analysis of a questionnaire, the thesis identifies the importance of uncertainty in archaeological interpretation and discovers potential preferences among different evidence types. Secondly, the thesis uniquely analyses and evaluates, in relation to archaeological uncertainty, three different belief quantification models. The varying ways that these mathematical models work, are also evaluated through simulated experiments. Comparison of results indicates significant convergence between the models. Thirdly, a novel approach to archaeological uncertainty and evidence conflict visualisation is presented, influenced by information visualisation schemes. Lastly, suggestions for future semantic extensions to this research are presented through the design and development of new plugins to a search engine

System- and Data-Driven Methods and Algorithms

An increasing complexity of models used to predict real-world systems leads to the need for algorithms to replace complex models with far simpler ones, while preserving the accuracy of the predictions. This two-volume handbook covers methods as well as applications. This first volume focuses on real-time control theory, data assimilation, real-time visualization, high-dimensional state spaces and interaction of different reduction techniques

Numerical Simulation of Frictional Contact Problems using Nagata Patches in Surface Smoothing

Tese de doutoramento em Engenharia Mecânica, na especialidade de Tecnologias de Produção, apresentada ao Departamento de Engenharia Mecânica da Faculdade de Ciências e Tecnologia da Universidade de CoimbraAll movements in the world involve contact and friction, from walking to car driving. The contact mechanics has application in many engineering problems, including the connection of structural members by bolts or screws, design of gears and bearings, sheet metal or bulk forming, rolling contact of car tyres, crash analysis of structures, as well as prosthetics in biomedical engineering. Due to the nonlinear and non-smooth nature of contact mechanics (contact area is not known a priori), such problems are currently solved using the finite element method within the field of computational contact mechanics. However, most of the commercial finite element software packages presently available are not entirely capable to solve frictional contact problems, demanding for efficient and robust methods. Therefore, the main objective of this study is the development of algorithms and numerical methods to apply in the numerical simulation of 3D frictional contact problems between bodies undergoing large deformations. The presented original developments are implemented in the in-house finite element code DD3IMP. The formulation of quasi-static frictional contact problems between bodies undergoing large deformations is firstly presented in the framework of the continuum mechanics, following the classical scheme used in solid mechanics. The kinematic description of the deformable bodies is presented adopting an updated Lagrangian formulation. The mechanical behaviour of the bodies is described by an elastoplastic constitutive law in conjunction with an associated flow rule, allowing to model a wide variety of contact problems arising in industrial applications. The frictional contact between the bodies is established by means of two conditions: the principle of impenetrability and the Coulomb’s friction law, both imposed to the contact interface. The augmented Lagrangian method is applied for solving the constrained minimization incremental problem resulting from the frictional contact inequalities, yielding a mixed functional involving both displacements and contact forces. The spatial discretization of the bodies is performed with isoparametric solid finite elements, while the discretization of the contact interface is carried out using the classical Node-to-Segment technique, preventing the slave nodes from penetrating on the master surface. The geometrical part of the contact elements, defined by a slave node and the closest master segment, is created by the contact search algorithm based on the selection of the closest point on the master surface, defined by the normal projection of the slave node. In the particular case of contact between a deformable body and a rigid obstacle, the master rigid surface can be described by smooth parameterizations typically used in CAD models. However, in the general case of contact between deformable bodies, the spatial discretization of both bodies with low order finite elements yields a piecewise bilinear representation of the master surface. This is the central source of problems in solving contact problems involving large sliding, since it leads to the discontinuity of the surface normal vector field. Thus, a surface smoothing procedure based on the Nagata patch interpolation is proposed to describe the master contact surfaces, which led to the development of the Node-to-Nagata contact element. The accuracy of the surface smoothing method using Nagata patches is evaluated by means of simple geometries. The nodal normal vectors required for the Nagata interpolation are evaluated from the CAD geometry in case of rigid master surfaces, while in case of deformable bodies they are approximated using the weighted average of the normal vectors of the neighbouring facets. The residual vectors and tangent matrices of the contact elements are derived coherently with the surface smoothing approach, allowing to obtain quadratic convergence rate in the generalized Newton method used for solving the nonlinear system of equations. The developed surface smoothing method and corresponding contact elements are validated through standard numerical examples with known analytical or semi-analytical solutions. More advanced frictional contact problems are studied, covering the contact of a deformable body with rigid obstacles and the contact between deformable bodies, including self-contact phenomena. The smoothing of the master surface improves the robustness of the computational methods and reduces strongly the non-physical oscillations in the contact force introduced by the traditional faceted description of the contact surface. The presented results are compared with numerical solutions obtained by other authors and experimental results, demonstrating the accuracy and performance of the implemented algorithms for highly nonlinear problems.Todos os movimentos no mundo envolvem contato e atrito, desde andar até conduzir um carro. A mecânica do contacto tem aplicação em muitos problemas de engenharia, incluindo a ligação de elementos estruturais com parafusos, projeto de engrenagens e rolamentos, estampagem ou forjamento, contato entre os pneus e a estrada, colisão de estruturas, bem como o desenvolvimento de próteses em engenharia biomédica. Devido à natureza não-linear e não-suave da mecânica do contato (área de contato desconhecida a priori), tais problemas são atualmente resolvidos usando o método dos elementos finitos no domínio da mecânica do contato computacional. No entanto, a maioria dos programas comerciais de elementos finitos atualmente disponíveis não é totalmente capaz de resolver problemas de contato com atrito, exigindo métodos numéricos mais eficientes e robustos. Portanto, o principal objetivo deste estudo é o desenvolvimento de algoritmos e métodos numéricos para aplicar na simulação numérica de problemas de contato com atrito entre corpos envolvendo grandes deformações. Os desenvolvimentos apresentados são implementados no programa de elementos finitos DD3IMP. A formulação quasi-estática de problemas de contato com atrito entre corpos deformáveis envolvendo grandes deformações é primeiramente apresentada no âmbito da mecânica dos meios contínuos, seguindo o método clássico usado em mecânica dos sólidos. A descrição cinemática dos corpos deformáveis é apresentada adotando uma formulação Lagrangeana reatualizada. O comportamento mecânico dos corpos é descrito por uma lei constitutiva elastoplástica em conjunto com uma lei de plasticidade associada, permitindo modelar uma grande variedade de problemas de contacto envolvidos em aplicações industriais. O contacto com atrito entre os corpos é definido por duas condições: o princípio da impenetrabilidade e a lei de atrito de Coulomb, ambas impostas na interface de contato. O método do Lagrangeano aumentado é aplicado na resolução do problema de minimização com restrições resultantes das condições de contato e atrito, produzindo uma formulação mista envolvendo deslocamentos e forças de contato. A discretização espacial dos corpos é realizada com elementos finitos sólidos isoparamétricos, enquanto a discretização da interface de contacto é realizado utilizando a técnica Node-to-Segment, impedindo os nós slave de penetrar na superfície master. A parte geométrica do elemento de contacto, definida por um nó slave e o segmento master mais próximo, é criada pelo algoritmo de deteção de contacto com base na seleção do ponto mais próximo na superfície master, obtido pela projeção normal do nó slave. No caso particular de contato entre um corpo deformável e um obstáculo rígido, a superfície rígida master pode ser descrita por parametrizações normalmente utilizadas em modelos CAD. No entanto, no caso geral de contato entre corpos deformáveis, a discretização espacial dos corpos com elementos finitos lineares origina uma representação da superfície master por facetas. Esta é a principal fonte de problemas na resolução de problemas de contato envolvendo grandes escorregamentos, uma vez que a distribuição dos vetor normais à superfície é descontínua. Assim, é proposto um método de suavização para descrever as superfícies de contacto master baseado na interpolação Nagata, que conduziu ao desenvolvimento do elemento de contacto Node-to-Nagata. A precisão do método de suavização das superfícies é avaliada através de geometrias simples. Os vetores normais nodais necessários para a interpolação Nagata são avaliados a partir da geometria CAD no caso de superfícies rígidas, enquanto no caso de corpos deformáveis são aproximados utilizando a média ponderada dos vetores normais das facetas vizinhas. Tanto os vetores de segundo membro como as matrizes residuais tangentes dos elementos de contacto são obtidas de forma coerente com o método de suavização da superfície, permitindo obter convergência quadrática no método de Newton generalizado, o qual é utilizado para resolver o sistema de equações não lineares. O método de suavização das superfícies e os elementos de contacto desenvolvidos são validados através de exemplos com soluções analíticas ou semi-analíticas conhecidas. Também são estudados outros problemas de contato mais complexos, incluindo o contato de um corpo deformável com obstáculos rígidos e o contato entre corpos deformáveis, contemplando fenómenos de auto-contato. A suavização da superfície master melhora a robustez dos métodos computacionais e reduz fortemente as oscilações na força de contato, associadas à descrição facetada da superfície de contato. Os resultados são comparados com soluções numéricas de outros autores e com resultados experimentais, demonstrando a precisão e o desempenho dos algoritmos implementados para problemas fortemente não-lineares.Fundação para a Ciência e Tecnologia - SFRH/BD/69140/201

Grasp plannind under task-specific contact constraints

Several aspects have to be addressed before realizing the dream of a robotic hand-arm system with human-like capabilities, ranging from the consolidation of a proper mechatronic design, to the development of precise, lightweight sensors and actuators, to the efficient planning and control of the articular forces and motions required for interaction with the environment. This thesis provides solution algorithms for a main problem within the latter aspect, known as the {\em grasp planning} problem: Given a robotic system formed by a multifinger hand attached to an arm, and an object to be grasped, both with a known geometry and location in 3-space, determine how the hand-arm system should be moved without colliding with itself or with the environment, in order to firmly grasp the object in a suitable way. Central to our algorithms is the explicit consideration of a given set of hand-object contact constraints to be satisfied in the final grasp configuration, imposed by the particular manipulation task to be performed with the object. This is a distinguishing feature from other grasp planning algorithms given in the literature, where a means of ensuring precise hand-object contact locations in the resulting grasp is usually not provided. These conventional algorithms are fast, and nicely suited for planning grasps for pick-an-place operations with the object, but not for planning grasps required for a specific manipulation of the object, like those necessary for holding a pen, a pair of scissors, or a jeweler's screwdriver, for instance, when writing, cutting a paper, or turning a screw, respectively. To be able to generate such highly-selective grasps, we assume that a number of surface regions on the hand are to be placed in contact with a number of corresponding regions on the object, and enforce the fulfilment of such constraints on the obtained solutions from the very beginning, in addition to the usual constraints of grasp restrainability, manipulability and collision avoidance. The proposed algorithms can be applied to robotic hands of arbitrary structure, possibly considering compliance in the joints and the contacts if desired, and they can accommodate general patch-patch contact constraints, instead of more restrictive contact types occasionally considered in the literature. It is worth noting, also, that while common force-closure or manipulability indices are used to asses the quality of grasps, no particular assumption is made on the mathematical properties of the quality index to be used, so that any quality criterion can be accommodated in principle. The algorithms have been tested and validated on numerous situations involving real mechanical hands and typical objects, and find applications in classical or emerging contexts like service robotics, telemedicine, space exploration, prosthetics, manipulation in hazardous environments, or human-robot interaction in general

Collision Detection and Merging of Deformable B-Spline Surfaces in Virtual Reality Environment

This thesis presents a computational framework for representing, manipulating and merging rigid and deformable freeform objects in virtual reality (VR) environment. The core algorithms for collision detection, merging, and physics-based modeling used within this framework assume that all 3D deformable objects are B-spline surfaces. The interactive design tool can be represented as a B-spline surface, an implicit surface or a point, to allow the user a variety of rigid or deformable tools. The collision detection system utilizes the fact that the blending matrices used to discretize the B-spline surface are independent of the position of the control points and, therefore, can be pre-calculated. Complex B-spline surfaces can be generated by merging various B-spline surface patches using the B-spline surface patches merging algorithm presented in this thesis. Finally, the physics-based modeling system uses the mass-spring representation to determine the deformation and the reaction force values provided to the user. This helps to simulate realistic material behaviour of the model and assist the user in validating the design before performing extensive product detailing or finite element analysis using commercially available CAD software. The novelty of the proposed method stems from the pre-calculated blending matrices used to generate the points for graphical rendering, collision detection, merging of B-spline patches, and nodes for the mass spring system. This approach reduces computational time by avoiding the need to solve complex equations for blending functions of B-splines and perform the inversion of large matrices. This alternative approach to the mechanical concept design will also help to do away with the need to build prototypes for conceptualization and preliminary validation of the idea thereby reducing the time and cost of concept design phase and the wastage of resources