Expansive Motions and the Polytope of Pointed Pseudo-Triangulations

We introduce the polytope of pointed pseudo-triangulations of a point set in the plane, defined as the polytope of infinitesimal expansive motions of the points subject to certain constraints on the increase of their distances. Its 1-skeleton is the graph whose vertices are the pointed pseudo-triangulations of the point set and whose edges are flips of interior pseudo-triangulation edges. For points in convex position we obtain a new realization of the associahedron, i.e., a geometric representation of the set of triangulations of an n-gon, or of the set of binary trees on n vertices, or of many other combinatorial objects that are counted by the Catalan numbers. By considering the 1-dimensional version of the polytope of constrained expansive motions we obtain a second distinct realization of the associahedron as a perturbation of the positive cell in a Coxeter arrangement. Our methods produce as a by-product a new proof that every simple polygon or polygonal arc in the plane has expansive motions, a key step in the proofs of the Carpenter's Rule Theorem by Connelly, Demaine and Rote (2000) and by Streinu (2000).Comment: 40 pages, 7 figures. Changes from v1: added some comments (specially to the "Further remarks" in Section 5) + changed to final book format. This version is to appear in "Discrete and Computational Geometry -- The Goodman-Pollack Festschrift" (B. Aronov, S. Basu, J. Pach, M. Sharir, eds), series "Algorithms and Combinatorics", Springer Verlag, Berli

Geometric auxetics

We formulate a mathematical theory of auxetic behavior based on one-parameter deformations of periodic frameworks. Our approach is purely geometric, relies on the evolution of the periodicity lattice and works in any dimension. We demonstrate its usefulness by predicting or recognizing, without experiment, computer simulations or numerical approximations, the auxetic capabilities of several well-known structures available in the literature. We propose new principles of auxetic design and rely on the stronger notion of expansive behavior to provide an infinite supply of planar auxetic mechanisms and several new three-dimensional structures

Pseudo-Triangulations, Rigidity and Motion Planning

This paper proposes a combinatorial approach to planning non-colliding trajectories for a polygonal bar-and-joint framework with n vertices. It is based on a new class of simple motions induced by expansive one-degree-of-freedom mechanisms, which guarantee noncollisions by moving all points away from each other. Their combinatorial structure is captured by pointed pseudo-triangulations, a class of embedded planar graphs for which we give several equivalent characterizations and exhibit rich rigidity theoretic properties. The main application is an efficient algorithm for the Carpenter\u27s Rule Problem: convexify a simple bar-and-joint planar polygonal linkage using only non-self-intersecting planar motions. A step of the algorithm consists in moving a pseudo-triangulation-based mechanism along its unique trajectory in configuration space until two adjacent edges align. At the alignment event, a local alteration restores the pseudo-triangulation. The motion continues for O(n3) steps until all the points are in convex position. © 2005 Springer Science+Business Media, Inc

Enumerating Constrained Non-Crossing Minimally Rigid Frameworks

In this paper we present an algorithm for enumerating without repetitions all the non-crossing generically minimally rigid bar-and-joint frameworks under edge constraints, which we call constrained non-crossing Laman frameworks, on a given set of n points in the plane. Our algorithm is based on the reverse search paradigm of Avis and Fukuda. It generates each output graph in O(n4) time and O(n) space, or, with a slightly different implementation, in O(n3) time and O(n2) space. In particular, we obtain that the set of all the constrained non-crossing Laman frameworks on a given point set is connected by flips which preserve the Laman property

Euclidean distance geometry and applications

Euclidean distance geometry is the study of Euclidean geometry based on the concept of distance. This is useful in several applications where the input data consists of an incomplete set of distances, and the output is a set of points in Euclidean space that realizes the given distances. We survey some of the theory of Euclidean distance geometry and some of the most important applications: molecular conformation, localization of sensor networks and statics.Comment: 64 pages, 21 figure

Unified Lagrangian formulation for fluid and solid mechanics, fluid-structure interaction and coupled thermal problems using the PFEM

The objective of this thesis is the derivation and implementation of a unified Finite Element formulation for the solution of uid and solid mechanics, Fluid-Structure Interaction (FSI) and coupled thermal problems. The unified procedure is based on a stabilized velocity-pressure Lagrangian formulation. Each time step increment is solved using a two-step Gauss-Seidel scheme: first the linear momentum equations are solved for the velocity increments, next the continuity equation is solved for the pressure in the updated configuration. The Particle Finite Element Method (PFEM) is used for the fluid domains, while the Finite Element Method (FEM) is employed for the solid ones. As a consequence, the domain is remeshed only in the parts occupied by the fluid. Linear shape functions are used for both the velocity and the pressure fields. In order to deal with the incompressibility of the materials, the formulation has been stabilized using an updated version of the Finite Calculus (FIC) method. The procedure has been derived for quasi-incompressible Newtonian fluids. In this work, the FIC stabilization procedure has been extended also to the analysis of quasi-incompressible hypoelastic solids. Specific attention has been given to the study of free surface flow problems. In particular, the mass preservation feature of the PFEM-FIC stabilized procedure has been deeply studied with the help of several numerical examples. Furthermore, the conditioning of the problem has been analyzed in detail describing the effect of the bulk modulus on the numerical scheme. A strategy based on the use of a pseudo bulk modulus for improving the conditioning of the linear system is also presented. The unified formulation has been validated by comparing its numerical results to experimental tests and other numerical solutions for fluid and solid mechanics, and FSI problems. The convergence of the scheme has been also analyzed for most of the problems presented. The unified formulation has been coupled with the heat tranfer problem using a staggered scheme. A simple algorithm for simulating phase change problems is also described. The numerical solution of several FSI problems involving the temperature is given. The thermal coupled scheme has been used successfully for the solution of an industrial problem. The objective of study was to analyze the damage of a nuclear power plant pressure vessel induced by a high viscous fluid at high temperature, the corium. The numerical study of this industrial problem has been included in the thesis.El objectivo de la presente tesis es la derivación e implementación de una formulación unificada con elementos finitos para la solución de problemas de mecánica de fluidos y de sólidos, interacción fluido-estructura (Fluid-Structure Interaction (FSI)) y con acoplamiento térmico. El método unificado està basado en una formulación Lagrangiana estabilizada y las variables incognitas son las velocidades y la presión. Cada paso de tiempo se soluciona a través de un esquema de dos pasos de tipo Gauss-Seidel. Primero se resuelven las ecuaciones de momento lineal por los incrementos de velocidad, luego se calculan las presiones en la configuración actualizada usando la ecuación de continuidad. Para los dominios fluidos se utiliza el método de elementos finitos de partículas (Particle Finite Element Method (PFEM)) mientras que los sólidos se solucionan con el método de elementos finitos (Finite Element Method (FEM)). Por lo tanto, se ramalla sólo las partes del dominio ocupadas por el fluido. Los campos de velocidad y presión se interpolan con funciones de forma lineales. Para poder analizar materiales incompresibles, la formulación ha sido estabilizada con una nueva versión del método Finite Calculus (FIC). La técnica de estabilización ha sido derivada para fluidos Newtonianos casi-incompresibles. En este trabajo, la estabilización con FIC se usa también para el análisis de sólidos hipoelásticos casi-incompresibles. En la tesis se dedica particular atención al estudio de flujo con superficie libre. En particular, se analiza en profundidad el tema de las pérdidas de masa y se muestra con varios ejemplos numéricos la capacidad del método de garantizar la conservación de masa en problemas de flujos en supeficie libre. Además se estudia con detalle el condicionamiento del esquema numérico analizando particularmente el efecto del módulo de compresibilidad. Se presenta también una estrategia basada en el uso de un pseudo módulo de compresibilidad para mejorar el condicionamiento del problema. La formulación unificada ha sido validada comparando sus resultados numéricos con pruebas de laboratorio y resultados numéricos de otras formulaciones. En la mayoría de los ejemplos también se ha estudiado la convergencia del método. En la tesis también se describe una estrategia segregada para el acoplamiento de la formulación unificada con el problema de transmisión de calor. Además se presenta una simple estrategia para simular el cambio de fase. El esquema acoplado ha sido utilizado para resolver varios problemas de FSI donde se incluye la temperatura y su efecto. El esquema acoplado con el problema térmico ha sido utilizado con éxito para resolver un problema industrial. El objetivo del estudio era la simulación del daño y la fusión de la vasija de un reactor nuclear provocados por el contacto con un fluido altamente viscoso y a gran temperatura. En la tesis se describe con detalle el estudio numérico realizado para esta aplicación industrialPostprint (published version

Intraoperative Endoscopic Augmented Reality in Third Ventriculostomy

In neurosurgery, as a result of the brain-shift, the preoperative patient models used as a intraoperative reference change. A meaningful use of the preoperative virtual models during the operation requires for a model update. The NEAR project, Neuroendoscopy towards Augmented Reality, describes a new camera calibration model for high distorted lenses and introduces the concept of active endoscopes endowed with with navigation, camera calibration, augmented reality and triangulation modules