Physically-based forehead animation including wrinkles

Physically-based animation techniques enable more realistic and accurate animation to be created. We present a fully physically-based approach for efficiently producing realistic-looking animations of facial movement, including animation of expressive wrinkles. This involves simulation of detailed voxel-based models using a graphics processing unit-based total Lagrangian explicit dynamic finite element solver with an anatomical muscle contraction model, and advanced boundary conditions that can model the sliding of soft tissue over the skull. The flexibility of our approach enables detailed animations of gross and fine-scale soft-tissue movement to be easily produced with different muscle structures and material parameters, for example, to animate different aged skins. Although we focus on the forehead, our approach can be used to animate any multi-layered soft body

Comparing small visual differences between conforming meshes

This paper gives a method of quantifying small visual differences between 3D mesh models with conforming topology, based on the theory of strain fields. Our experiments show that our difference estimates are well correlated with human perception of differences. This work has applications in the evaluation of 3D mesh watermarking, 3D mesh compression reconstruction, and 3D mesh filtering

Discontinuities without discontinuity: The Weakly-enforced Slip Method

Tectonic faults are commonly modelled as Volterra or Somigliana dislocations in an elastic medium. Various solution methods exist for this problem. However, the methods used in practice are often limiting, motivated by reasons of computational efficiency rather than geophysical accuracy. A typical geophysical application involves inverse problems for which many different fault configurations need to be examined, each adding to the computational load. In practice, this precludes conventional finite-element methods, which suffer a large computational overhead on account of geometric changes. This paper presents a new non-conforming finite-element method based on weak imposition of the displacement discontinuity. The weak imposition of the discontinuity enables the application of approximation spaces that are independent of the dislocation geometry, thus enabling optimal reuse of computational components. Such reuse of computational components renders finite-element modeling a viable option for inverse problems in geophysical applications. A detailed analysis of the approximation properties of the new formulation is provided. The analysis is supported by numerical experiments in 2D and 3D.Comment: Submitted for publication in CMAM

Error estimation and adaptive mesh refinement for parallel analysis of shell structures

The formulation and application of element-level, element-independent error indicators is investigated. This research culminates in the development of an error indicator formulation which is derived based on the projection of element deformation onto the intrinsic element displacement modes. The qualifier 'element-level' means that no information from adjacent elements is used for error estimation. This property is ideally suited for obtaining error values and driving adaptive mesh refinements on parallel computers where access to neighboring elements residing on different processors may incur significant overhead. In addition such estimators are insensitive to the presence of physical interfaces and junctures. An error indicator qualifies as 'element-independent' when only visible quantities such as element stiffness and nodal displacements are used to quantify error. Error evaluation at the element level and element independence for the error indicator are highly desired properties for computing error in production-level finite element codes. Four element-level error indicators have been constructed. Two of the indicators are based on variational formulation of the element stiffness and are element-dependent. Their derivations are retained for developmental purposes. The second two indicators mimic and exceed the first two in performance but require no special formulation of the element stiffness mesh refinement which we demonstrate for two dimensional plane stress problems. The parallelizing of substructures and adaptive mesh refinement is discussed and the final error indicator using two-dimensional plane-stress and three-dimensional shell problems is demonstrated

Direct medical image-based Finite Element modelling for patient-specific simulation of future implants

[EN] In patient specific biomedical simulation, the numerical model is usually created after cumbersome, time consuming procedures which often require highly specialized human work and a great amount of man-hours to be carried out. In order to make numerical simulation available for medical practice, it is of primary importance to reduce the cost associated to these procedures by making them automatic. In this paper a method for the automatic creation of Finite Element (FE) models from medical images is presented. This method is based on the use of a hierarchical structure of nested Cartesian grids in which the medical image is immersed. An efficient h-adaptive procedure conforms the FE model to the image characteristics by refining the mesh on the basis of the distribution of elastic properties associated to the pixel values. As a result, a problem with a reasonable number of degrees of freedom is obtained, skipping the geometry creation stage. All the image information is taken into account during the calculation of the element stiffness matrix, therefore it is straightforward to include the material heterogeneity in the simulation. The proposed method is an adapted version of the Cartesian grid Finite Element Method (cgFEM) for the FE analysis of objects defined by images. cgFEM is an immersed boundary method that uses h-adaptive Cartesian meshes non-conforming to the boundary of the object to be analysed. The proposed methodology, used together with the original geometry-based cgFEM, allows prosthesis geometries to be easily introduced in the model providing a useful tool for evaluating the effect of future implants in a preoperative framework. The potential of this kind of technology is presented by mean of an initial implementation in 2D and 3D for linear elasticity problems.With the support of the European Union Framework Programme (FP7) under grant agreement No. 289361 'Integrating Numerical Simulation and Geometric Design Technology (INSIST)', the Ministerio de Economia y Competitividad of Spain (DPI2010-20542) and the Generalitat Valenciana (PROMETEO/2016/007).Giovannelli, L.; Ródenas, J.; Navarro-Jiménez, J.; Tur Valiente, M. (2017). Direct medical image-based Finite Element modelling for patient-specific simulation of future implants. Finite Elements in Analysis and Design. 136:37-57. https://doi.org/10.1016/j.finel.2017.07.010S375713

Local enrichment of NURBS patches using a non-intrusive coupling strategy: Geometric details, local refinement, inclusion, fracture

International audienceIn this work, we apply a non-intrusive global/local coupling strategy for the modelling of local phenomena in a NURBS patch. The idea is to consider the NURBS patch to be enriched as the global model. This results in a simple, flexible strategy: first, the global NURBS patch remains unchanged, which completely eliminates the need for costly re-parametrization procedures (even if the local domain is expected to evolve); then, easy merging of a linear NURBS code with any other existing robust codes suitable for the modelling of complex local behaviour is possible. The price to pay is the number of iterations of the non-intrusive solver but we show that this can be strongly reduced by means of acceleration techniques. The main development for NURBS is to be able to handle non-conforming geometries. Only slight changes in the implementation process, including the setting up of suitable quadrature rules for the evaluation of the interface reaction forces, are made in response to this issue. A range of numerical examples in two-dimensional linear elasticity are given to demonstrate the performance of the proposed methodology and its significant potential to treat any case of local enrichment in a NURBS patch simply

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