We propose a novel method for obtaining more accurate tangential velocities for solid fluid coupling. Our method works for both rigid and deformable objects as well as both volumetric objects and thin shells. The fluid can be either one phase such as smoke or two phase such as water with a free surface. The coupling between the solid and the fluid can either be one-way with kinematic solids or fully two-way coupled. The only previous scheme that was general enough to handle both two-way coupling and thin shells required a mass lumping strategy that did not allow for freely flowing tangential velocities. Similar to that previous work, our method prevents leaking of fluid across a thin shell, however unlike that work our method does not couple the tangential velocities in any fashion, allowing for the proper slip independently on each side of the body. Moreover, since it accurately and directly treats the tangential velocity, it does not rely on grid refinement to obtain a reasonable solution. Therefore, it gives a highly improved result on coarse meshes.
Category: Deformables
Real-Time Deformation and Fracture in a Game Environment
This paper describes a simulation system that has been developed to model the deformation and fracture of solid objects in a real-time gaming context. Based around a corotational tetrahedral finite element method, this system has been constructed from components published in the graphics and computational physics literatures. The goal of this paper is to describe how these components can be combined to produce an engine that is robust to unpredictable user interactions, fast enough to model reasonable scenarios at real-time speeds, suitable for use in the design of a game level, and with appropriate controls allowing content creators to match artistic direction. Details concerning parallel implementation, solver design, rendering method, and other aspects of the simulation are elucidated with the intent of providing a guide to others wishing to implement similar systems. Examples from in-game scenes captured on the Xbox 360, PS3, and PC platforms are included.
Real-Time Deformation and Fracture in a Game Environment
A Point-based Method for Animating Elastoplastic Solids
In this paper we describe a point-based approach for animating elastoplastic materials. Our primary contribution is a simple method for computing the deformation gradient for each particle in the simulation. The deformation gradient is computed for each particle by finding the affine transformation that best approximates the motion of neighboring particles over a single timestep. These transformations are then composed to compute the total deformation gradient that describes the deformation around a particle over the course of the simulation. Given the deformation gradient we can apply arbitrary constitutive models and compute the resulting elastic forces. Our method has two primary advantages: we do not store or compare to an initial rest configuration and we work directly with the deformation gradient. The first advantage avoids poor numerical conditioning and the second naturally leads to a multiplicative model of deformation appropriate for finite deformations. We demonstrate our approach on a number of examples that exhibit a wide range of material behaviors.
Simple, yet accurate tensile stiffness
Recent Particle System models have evolved toward accurate representation of elastic stiffness based on continuum mechanics, converging to formulations that make them quite analogous to fast Finite Element methods. These formulations usually involve the linearization of tensors that help their formulation in the context of linear elasticity. Toward our objective of simulating the nonlinear properties of cloth accurately, we show through this work that this linearization can indeed be suppressed and replaced by adapted strain-stress laws relating precisely the nonlinear behavior of the material. This leads to very streamlined computations that are particularly efficient for simulating the nonlinear anisotropic tensile elasticity of deformable surfaces. Through a simple and efficient implementation using the Particle System formalism, we demonstrate the efficiency of this method with examples related to garment simulation.
An Edge-Based Computationally Efficient Formulation of Saint-Venant Kirchhoff Tetrahedral Finite Elements
This article describes a computationally efficient formulation and an algorithm for tetrahedral finite-element simulation of elastic objects subject to Saint Venant-Kirchhoff (StVK) material law. The number of floating point operations required by the algorithm is in the range of 15% to 27% for computing the vertex forces from a given set of vertex positions, and 27% to 38% for the tangent stiffness matrix, in comparison to a well-optimized algorithm directly derived from the conventional Total Lagrangian formulation. In the new algorithm, the data is associated with edges and tetrahedron-sharing edge-pairs (TSEPs), as opposed to tetrahedra, to avoid redundant computation. Another characteristic of the presented formulation is that it reduces to that of a spring-network model by simply ignoring all the TSEPs. The technique is demonstrated through an interactive application involving haptic interaction, being combined with a linearized implicit integration technique employing a preconditioned conjugate gradient method.
Implicit Contact Handling for Deformable Objects
We present an algorithm for robust and efficient contact handling of deformable objects. By being aware of the internal dynamics of the colliding objects, our algorithm provides smooth rolling and sliding, stable stacking, robust impact handling, and seamless coupling of heterogeneous objects, all in a unified manner. We achieve dynamicsawareness through a constrained dynamics formulation with implicit complementarity constraints, and we present two major contributions that enable an efficient solution of the constrained dynamics problem: a time stepping algorithm that robustly ensures non-penetration and progressively refines the formulation of constrained dynamics, and a new solver for large mixed linear complementarity problems, based on iterative constraint anticipation. We show the application of our algorithm in challenging scenarios such as multi-layered cloth moving at high velocities, or colliding deformable solids simulated with large time steps.
A Geometric Deformation Model for Stable Cloth Simulation
We propose an adapted shape-matching approach for the efficient and robust simulation of clothing. A combination of two different cluster types is employed to account for high stretching and shearing, and low bending resistance. Due to the inherent handling of overshooting issues, the proposed shape-matching deformation model is robust. The proposed cluster types allow for a computationally efficient handling of bending. The geometric deformation model is combined with a novel collision handling approach. The technique employs spatial subdivision to detect collisions and self-collisions. The response scheme is derived from an existing approach for elastic rods. To illustrate the physically plausible dynamics of our approach, it is compared to a traditional physically-based deformation model. Experiments indicate that similar cloth properties can be reproduced with both models. The computational efficiency of the proposed scheme enables the interactive animation of clothing and shells.
Seams and Bending in Cloth Simulation
Accurate modeling of bending behavior is one of the most important tasks in the field of cloth simulation. Bending stiffness is probably the most significant material parameter describing a given textile. Much work has been done in recent years to allow a fast and authentic reproduction of the effect of bending in cloth simulation systems. However, these approaches usually treat the textiles as consisting of a single, homogeneous material. The effects of seams, interlining and multilayer materials have not been considered so far. Recent work showed that the bending stiffness of a textile is greatly influenced by the presence of seams and that a good cloth simulation system needs to consider these effects.
In this work we show how accurate modeling of bending and seams can be achieved in a state-of-the-art cloth simulation system. Our system can make use of measured bending stiffness data, but also allows intuitive user control, if desired. We verify our approach using virtual draping tests and garments in the simulation and comparing the results to their real-world counterparts. Furthermore, we provide heuristics derived from measurements that can be used to approximate the influence of several common types of seams.
Corotated Finite Elements Made Fast and Stable
Multigrid finite-element solvers using the corotational formulation of finite elements provide an attractive means for the simulation of deformable bodies exhibiting linear elastic response. The separation of rigid body motions from the total element motions using purely geometric methods or polar decomposition of the deformation gradient, however, can introduce instabilities for large element rotations and deformations. Furthermore, the integration of the corotational formulation into dynamic multigrid elasticity simulations requires to continually rebuild consistent system matrices at different resolution levels. The computational load imposed by these updates prohibits the use of large numbers of finite elements at rates comparable to the small-strain finite element formulation. To overcome the first problem, we present a new method to extract the rigid body motion from total finite element displacements based on energy minimization. This results in a very stable corotational formulation that only slightly increases the computational overhead. We address the second problem by introducing a novel algorithm for computing sparse products of the form RKRT, as they have to be evaluated to update the multigrid hierarchy. By reformulating the problem into the simultaneous processing of a sequential data and control stream, cache miss penalties are significantly reduced. Even though the algorithm increases memory requirements, it accelerates the multigrid FE simulation by a factor of up to 4 compared to previous multigrid approaches. Due to the proposed improvements, finite element deformable body simulations using the corotational formulation can be performed at rates of 17 tps for up to 12k elements.
Hierarchical Position Based Dynamics
The Position Based Dynamics approach (PBD) recently introduced allows robust simulations of dynamic systems in real time. The simplicity of the method is due to the fact, that the solver processes the constraints one by one in a Gauss-Seidel type manner. In contrast to global Newton-Raphson solvers, the local solver can easily handle non-linear constraints as well as constraints based on inequalities. Unfortunately, this advantage comes at the price of much slower convergence. In this paper we propose a multi-grid based process to speed up the convergence of PBD significantly while keeping the power of the method to process general non-linear constraints. Several examples show that the new approach is significantly faster than the original one. This makes real time simulation possible at a higher level of detail in interactive applications such as computer games.