# Soft body dynamics

Soft body dynamics is a field of computer graphics that focuses on visually realistic physical simulations of the motion and properties of deformable objects (or soft bodies).[1] The applications are mostly in video games and film. Unlike in simulation of rigid bodies, the shape of soft bodies can change, meaning that the relative distance of two points on the object is not fixed. While the relative distances of points are not fixed, the body is expected to retain its shape to some degree (unlike a fluid). The scope of soft body dynamics is quite broad, including simulation of soft organic materials such as muscle, fat, hair and vegetation, as well as other deformable materials such as clothing and fabric. Generally, these methods only provide visually plausible emulations rather than accurate scientific/engineering simulations, though there is some crossover with scientific methods, particularly in the case of finite element simulations. Several physics engines currently provide software for soft-body simulation.[2][3][4][5][6][7]

Softbody objects react to forces and are able to collide with other 3D objects. This example has been created with Blender.

## Deformable solids

The simulation of volumetric solid soft bodies can be realised by using a variety of approaches.

### Spring/mass models

Two nodes as mass points connected by a parallel circuit of a spring and a damper.

In this approach, the body is modeled as a set of point masses (nodes) connected by ideal weightless elastic springs obeying some variant of Hooke's law. The nodes may either derive from the edges of a two-dimensional polygonal mesh representation of the surface of the object, or from a three-dimensional network of nodes and edges modeling the internal structure of the object (or even a one-dimensional system of links, if for example a rope or hair strand is being simulated). Additional springs between nodes can be added, or the force law of the springs modified, to achieve desired effects. Applying Newton's second law to the point masses including the forces applied by the springs and any external forces (due to contact, gravity, air resistance, wind, and so on) gives a system of differential equations for the motion of the nodes, which is solved by standard numerical schemes for solving ODEs.[8] Rendering of a three-dimensional mass-spring lattice is often done using free-form deformation,[9] in which the rendered mesh is embedded in the lattice and distorted to conform to the shape of the lattice as it evolves. Assuming all point masses equal to zero one can obtain the Stretched grid method aimed at several engineering problems solution relative to the elastic grid behavior.

### Finite element simulation

This is a more physically accurate approach, which uses the widely used finite element method to solve the partial differential equations which govern the dynamics of an elastic material. The body is modeled as a three-dimensional elastic continuum by breaking it into a large number of solid elements which fit together, and solving for the stresses and strains in each element using a model of the material.[10] The elements are typically tetrahedral, the nodes being the vertices of the tetrahedra (relatively simple methods exist[11][12] to tetrahedralize a three dimensional region bounded by a polygon mesh into tetrahedra, similarly to how a two-dimensional polygon may be triangulated into triangles). The strain (which measures the local deformation of the points of the material from their rest state) is quantified by the strain tensor $\boldsymbol{\epsilon}$. The stress (which measures the local forces per-unit area in all directions acting on the material) is quantified by the Cauchy stress tensor $\boldsymbol{\sigma}$. Given the current local strain, the local stress can be computed via the generalized form of Hooke's law: $\boldsymbol{\sigma} = \mathsf{C} \boldsymbol{\varepsilon} \ ,$ where $\mathsf{C}$ is the "elasticity tensor" which encodes the material properties (parametrized in linear elasticity for an isotropic material by the Poisson ratio and Young's modulus).

The equation of motion of the element nodes is obtained by integrating the stress field over each element and relating this, via Newton's second law, to the node accelerations.

Pixelux (developers of the Digital Molecular Matter system) use a finite-element-based approach for their soft bodies, using a tetrahedral mesh and converting the stress tensor directly into node forces.[13] Rendering is done via a form of free-form deformation.[9]

### Energy minimization methods

This approach is motivated by variational principles and the physics of surfaces, which dictate that a constrained surface will assume the shape which minimizes the total energy of deformation (analogous to a soap bubble). Expressing the energy of a surface in terms of its local deformation (the energy is due to a combination of stretching and bending), the local force on the surface is given by differentiating the energy with respect to position, yielding an equation of motion which can be solved in the standard ways.[14][15]

### Shape matching

In this scheme, penalty forces or constraints are applied to the model to drive it towards its original shape[16] (i.e. the material behaves as if it has shape memory). To conserve momentum the rotation of the body must be estimated properly, for example via polar decomposition. To approximate finite element simulation, shape matching can be applied to three dimensional lattices and multiple shape matching constraints blended.[17]

### Rigid-body based deformation

Deformation can also be handled by a traditional rigid-body physics engine, modeling the soft-body motion using a network of multiple rigid bodies connected by constraints, and using (for example) matrix-palette skinning to generate a surface mesh for rendering. This is the approach used for deformable objects in Havok Destruction.[18]

## Cloth simulation

In the context of computer graphics, cloth simulation refers to the simulation of soft bodies in the form of two dimensional continuum elastic membranes, that is, for this purpose, the actual structure of real cloth on the yarn level can be ignored (though modeling cloth on the yarn level has been tried).[19] Via rendering effects, this can produce a visually plausible emulation of textiles and clothing, used in a variety of contexts in video games, animation, and film. It can also be used to simulate two dimensional sheets of materials other than textiles, such as deformable metal panels or vegetation. In video games it is often used to enhance the realism of clothed characters, which otherwise would be entirely animated.

Cloth simulators are generally based on mass-spring models, but a distinction must be made between force-based and position-based solvers.

### Force-based cloth

The mass-spring model (obtained from a polygonal mesh representation of the cloth) determines the internal spring forces acting on the nodes at each timestep (in combination with gravity and applied forces). Newton's second law gives equations of motion which can be solved via standard ODE solvers. To create high resolution cloth with a realistic stiffness is not possible however with simple explicit solvers (such as forward Euler integration), unless the timestep is made too small for interactive applications (since as is well known, explicit integrators are numerically unstable for sufficiently stiff systems). Therefore implicit solvers must be used,[20] requiring solution of a large sparse matrix system (via e.g. the conjugate gradient method), which itself may also be difficult to achieve at interactive frame rates. An alternative[21][22] is to use an explicit method with low stiffness, with ad hoc methods to avoid instability and excessive stretching (e.g. strain limiting corrections).

### Position-based dynamics

To avoid needing to do an expensive implicit solution of a system of ODEs, many real-time cloth simulators (notably PhysX, Havok Cloth, and Maya nCloth) use position based dynamics (PBD),[23] an approach based on constraint relaxation. The mass-spring model is converted into a system of constraints, which demands that the distance between the connected nodes be equal to the initial distance. This system is solved sequentially and iteratively, by directly moving nodes to satisfy each constraint, until sufficiently stiff cloth is obtained. This is similar to a Gauss-Seidel solution of the implicit matrix system for the mass-spring model. Care must be taken though to solve the constraints in the same sequence each timestep, to avoid spurious oscillations, and to make sure that the constraints do not violate linear and angular momentum conservation. Additional position constraints can be applied, for example to keep the nodes within desired regions of space (sufficiently close to an animated model for example), or to maintain the body's overall shape via shape matching.

## Collision detection for deformable objects

Realistic interaction of simulated soft objects with their environment may be important for obtaining visually realistic results. Cloth self-intersection is important in some applications for acceptably realistic simulated garments. This is challenging to achieve at interactive frame rates, particularly in the case of detecting and resolving self collisions and mutual collisions between two or more deformable objects.

Collision detection may be discrete/a posteriori (meaning objects are advanced in time through a pre-determined interval, and then any penetrations detected and resolved), or continuous/a priori (objects are advanced only until a collision occurs, and the collision is handled before proceeding). The former is easier to implement and faster, but leads to failure to detect collisions (or detection of spurious collisions) if objects move fast enough. Real-time systems generally have to use discrete collision detection, with other ad hoc ways to avoid failing to detect collisions.

Detection of collisions between cloth and environmental objects with a well defined "inside" is straightforward since the system can detect unambiguously whether the cloth mesh vertices and faces are intersecting the body and resolve them accordingly. If a well defined "inside" does not exist (e.g. in the case of collision with a mesh which does not form a closed boundary), an "inside" may be constructed via extrusion. Mutual- or self-collisions of soft bodies defined by tetrahedra is straightforward, since it reduces to detection of collisions between solid tetrahedra.

However, detection of collisions between two polygonal cloths (or collision of a cloth with itself) via discrete collision detection is much more difficult, since there is no unambiguous way to locally detect after a timestep whether a cloth node which has penetrated is on the "wrong" side or not. Solutions involve either using the history of the cloth motion to determine if an intersection event has occurred, or doing a global analysis of the cloth state to detect and resolve self-intersections. Pixar has presented a method which uses a global topological analysis of mesh intersections in configuration space to detect and resolve self-interpenetration of cloth.[24] Currently, this is generally too computationally expensive for real-time cloth systems.

To do collision detection efficiently, primitives which are certainly not colliding must be identified as soon as possible and discarded from consideration to avoid wasting time. To do this, some form of spatial subdivision scheme is essential, to avoid a brute force test of $O[n^2]$ primitive collisions. Approaches used include:

• Bounding volume hierarchies (AABB trees,[25] OBB trees, sphere trees)
• Grids, either uniform[26] (using hashing for memory efficiency) or hierarchical (e.g. Octree, kd-tree)
• Coherence-exploiting schemes, such as sweep and prune with insertion sort, or tree-tree collisions with front tracking.
• Hybrid methods involving a combination of various of these schemes, e.g. a coarse AABB tree plus sweep-and-prune with coherence between colliding leaves.

## Other applications

Other effects which may be simulated via the methods of soft-body dynamics are:

Simulating fluids in the context of computer graphics would not normally be considered soft-body dynamics, which is usually restricted to mean simulation of materials which have a tendency to retain their shape and form. In contrast, a fluid assumes the shape of whatever vessel contains it, as the particles are bound together by relatively weak forces.

## References

1. ^ Nealen, Müller, Keiser, Boxerman & Carlson (2005). "Physically Based Deformable Models in Computer Graphics". CiteSeerX: 10.1.1.124.4664.
2. ^
3. ^
4. ^
5. ^
6. ^
7. ^
8. ^ Desbrun, Schroder & Barr (1999). "Interactive animation of structured deformable objects".
9. ^ a b Sederberg & Parry (1986). "Free-Form Deformation of Solid Geometric Models".
10. ^ Kaufmann, Martin, Botsch & Gross (2008). "Flexible Simulation of Deformable Models Using Discontinuous Galerkin FEM".
11. ^ Spillmann, Wagner & Teschner (2006). "Robust Tetrahedral Meshing of Triangle Soups".
12. ^ Bridson, Teran, Molino & Fedkiw (2003). "Adaptive Physics Based Tetrahedral Mesh Generation Using Level Sets".
13. ^ a b Parker & O'Brien (2009). "Real-Time Deformation and Fracture in a Game Environment".
14. ^ Terzopoulos, Platt, Barr & Fleischer (1987). "Elastically Deformable Models".
15. ^ Wardetzky, Bergou, Harmon, Zorin & Grinspun (2007). "Discrete Quadratic Curvature Energies".
16. ^ a b Müller, Heidelberger, Teschner & Gross (2005). "Meshless Deformations Based on Shape Matching".
17. ^ Steinemann, Otaduy & Gross (2008). "Fast Adaptive Shape Matching Deformations".
18. ^
19. ^ Kaldor, James & Marschner (2008). "Simulating Knitted Cloth at the Yarn Level".
20. ^ Baraff & Witkin (1998). "Large Steps in Cloth Simulation".
21. ^ Provot (1997). "Collision and self-collision handling in cloth model dedicated to design garments". CiteSeerX: 10.1.1.89.9232.
22. ^ Bridson, Fedkiw & Anderson (2002). "Robust Treatment of Collisions, Contact and Friction for Cloth Animation".
23. ^ Müller, Heidelberger, Hennix & Ratcliff (2006). "Position Based Dynamics".
24. ^ Baraff, Witkin & Kass (2003). "Untangling Cloth".
25. ^ Van Den Bergen (1998). "Efficient Collision Detection of Complex Deformable Models using AABB Trees".
26. ^ Teschner, Heidelberger, Müller, Pomeranets & Gross (2003). "Optimized Spatial Hashing for Collision Detection of Deformable Objects".
27. ^ Steinemann, Otaduy & Gross (2006). "Fast Arbitrary Splitting of Deforming Objects".
28. ^ Müller, Keiser, Nealen, Pauly, Gross & Alexa (2004). "Point Based Animation of Elastic, Plastic and Melting Objects".
29. ^ Selle, Lentine & Fedkiw (2008). "A Mass Spring Model for Hair Simulation".
30. ^
31. ^ "Numerion Software Ltd". Retrieved 15 July 2011.