A Survey of Spatial Deformation from a User-Centered Perspective

The spatial deformation methods are a family of modeling and animation techniques for indirectly reshaping an object by warping the surrounding space, with results that are similar to molding a highly malleable substance. They have the virtue of being computationally efficient (and hence interactive) and applicable to a variety of object representations. In this paper we survey the state of the art in spatial deformation. Since manipulating ambient space directly is infeasible, deformations are controlled by tools of varying dimension - points, curves, surfaces and volumes - and it is on this basis that we classify them. Unlike previous surveys that concentrate on providing a single underlying mathematical formalism, we use the user-centered criteria of versatility, ease of use, efficiency and correctness to compare techniques

Isogeometric Approximation of Variational Problems for Shells

The interaction of applied geometry and numerical simulation is a growing field in the interplay of com- puter graphics, computational mechanics and applied mathematics known as isogeometric analysis. In this thesis we apply and analyze Loop subdivision surfaces as isogeometric tool because they provide great flexibility in handling surfaces of arbitrary topology combined with higher order smoothness. Compared with finite element methods, isogeometric methods are known to require far less degrees of freedom for the modeling of complex surfaces but at the same time the assembly of the isogeo- metric matrices is much more time-consuming. Therefore, we implement the isogeometric subdivision method and analyze the experimental convergence behavior for different quadrature schemes. The mid-edge quadrature combines robustness and efficiency, where efficiency is additionally increased via lookup tables. For the first time, the lookup tables allow the simulation with control meshes of arbitrary closed connectivity without an initial subdivision step, i.e. triangles can have more than one vertex with valence different from six. Geometric evolution problems have many applications in material sciences, surface processing and modeling, bio-mechanics, elasticity and physical simulations. These evolution problems are often based on the gradient flow of a geometric energy depending on first and second fundamental forms of the surface. The isogeometric approach allows a conforming higher order spatial discretization of these geometric evolutions. To overcome a time-error dominated scheme, we combine higher order space and time discretizations, where the time discretization based on implicit Runge-Kutta methods. We prove that the energy diminishes in every time-step in the fully discrete setting under mild time-step restrictions which is the crucial characteristic of a gradient flow. The overall setup allows for a general type of fourth-order energies. Among others, we perform experiments for Willmore flow with respect to different metrics. In the last chapter of this thesis we apply the time-discrete geodesic calculus in shape space to the space of subdivision shells. By approximating the squared Riemannian distance by a suitable energy, this approach defines a discrete path energy for a consistent computation of geodesics, logarithm and exponential maps and parallel transport. As approximation we pick up an elastic shell energy, which measures the deformation of a shell by membrane and bending contributions of its mid-surface. Bézier curves are a fundamental tool in computer-aided geometric design. We extend these to the subdivision shell space by generalizing the de Casteljau algorithm. The evaluation of Bézier curves depends on all input data. To solve this problem, we introduce B-splines and cardinal splines in shape space by gluing together piecewise Bézier curves in a smooth way. We show examples of quadratic and cubic Bézier curves, quadratic and cubic B-splines as well as cardinal splines in subdivision shell space

A geometric framework for immersogeometric analysis

The purpose of this dissertation is to develop a geometric framework for immersogeometric analysis that directly uses the boundary representations (B-reps) of a complex computer-aided design (CAD) model and immerses it into a locally refined, non-boundary-fitted discretization of the fluid domain. Using the non-boundary-fitted mesh which does not need to conform to the shape of the object can alleviate the challenge of mesh generation for complex geometries. This also reduces the labor-intensive and time-consuming work of geometry cleanup for the purpose of obtaining watertight CAD models in order to perform boundary-fitted mesh generation. The Dirichlet boundary conditions in the fluid domain are enforced weakly over the immersed object surface in the intersected elements. The surface quadrature points for the immersed object are generated on the parametric and analytic surfaces of the B-rep models. In the case of trimmed surfaces, adaptive quadrature rule is considered to improve the accuracy of the surface integral. For the non-boundary-fitted mesh, a sub-cell-based adaptive quadrature rule based on the recursive splitting of quadrature elements is used to faithfully capture the geometry in intersected elements. The point membership classification for identifying quadrature points in the fluid domain is based on a voxel-based approach implemented on GPUs. A variety of computational fluid dynamics (CFD) simulations are performed using the proposed method to assess its accuracy and efficiency. Finally, a fluid--structure interaction (FSI) simulation of a deforming left ventricle coupled with the heart valves shows the potential advantages of the developed geometric framework for the immersogeomtric analysis with complex moving domains

Isogeometric boundary element analysis using unstructured T-splines

We couple collocated isogeometric boundary element methods and unstructured analysis-suitable T-spline surfaces for linear elastostatic problems. We extend the definition of analysis-suitable T-splines to encompass unstructured control grids (unstructured meshes) and develop basis functions which are smooth (rational) polynomials defined in terms of the Bézier extraction framework and which pass standard patch tests. We then develop a collocation procedure which correctly accounts for sharp edges and corners, extraordinary points, and T-junctions. This approach is applied to several three-dimensional problems, including a real-world T-spline model of a propeller. We believe this work clearly illustrates the power of combining new analysis-suitable computer aided design technologies with established analysis methodologies, in this case, the boundary element method. © 2012 Elsevier B.V

Doctor of Philosophy

dissertationVolumetric parameterization is an emerging field in computer graphics, where volumetric representations that have a semi-regular tensor-product structure are desired in applications such as three-dimensional (3D) texture mapping and physically-based simulation. At the same time, volumetric parameterization is also needed in the Isogeometric Analysis (IA) paradigm, which uses the same parametric space for representing geometry, simulation attributes and solutions. One of the main advantages of the IA framework is that the user gets feedback directly as attributes of the NURBS model representation, which can represent geometry exactly, avoiding both the need to generate a finite element mesh and the need to reverse engineer the simulation results from the finite element mesh back into the model. Research in this area has largely been concerned with issues of the quality of the analysis and simulation results assuming the existence of a high quality volumetric NURBS model that is appropriate for simulation. However, there are currently no generally applicable approaches to generating such a model or visualizing the higher order smooth isosurfaces of the simulation attributes, either as a part of current Computer Aided Design or Reverse Engineering systems and methodologies. Furthermore, even though the mesh generation pipeline is circumvented in the concept of IA, the quality of the model still significantly influences the analysis result. This work presents a pipeline to create, analyze and visualize NURBS geometries. Based on the concept of analysis-aware modeling, this work focusses in particular on methodologies to decompose a volumetric domain into simpler pieces based on appropriate midstructures by respecting other relevant interior material attributes. The domain is decomposed such that a tensor-product style parameterization can be established on the subvolumes, where the parameterization matches along subvolume boundaries. The volumetric parameterization is optimized using gradient-based nonlinear optimization algorithms and datafitting methods are introduced to fit trivariate B-splines to the parameterized subvolumes with guaranteed order of accuracy. Then, a visualization method is proposed allowing to directly inspect isosurfaces of attributes, such as the results of analysis, embedded in the NURBS geometry. Finally, the various methodologies proposed in this work are demonstrated on complex representations arising in practice and research