The Human Touch: Measuring Contact with Real Human Soft Tissues

D. K. Pai, A. Rothwell, P. Wyder-Hodge, A. Wick, Y. Fan, E. Larionov, D. Harrison, D. R. Neog, and C. Shing

Simulating how the human body deforms in contact with clothing, wearables, and other objects is of central importance to many fields. However, the tissue material properties needed to accurately simulate real human bodies had been sorely lacking. We showed that these mechanical properties can be directly measured using a novel hand-held device. We have developed a complete pipeline for measurement, modelling, parameter estimation, and simulation using the finite element method. Our unique data may be used to create personalized models of an individual human or of a population. Consequently, our methods may have many potential applications in apparel design, e-commerce, computer animation, and medicine.

The Human Touch: Measuring Contact with Real Human Soft Tissues

Cable Joints

Matthias Müller, Nuttapong Chentanez, Stefan Jeschke, Miles Macklin

Robustly and efficiently simulating cables and ropes that are part of a larger system such as cable driven machines, cable cars or tendons in a human or robot is a challenging task. To be able to adapt to the environment, cables are typically modeled as a large number of small segments that are connected via joints. The two main difficulties with this approach are to satisfy the inextensibility constraint and to handle the typically large mass ratio between the small segments and the larger objects they connect. In this paper we present a new approach which solves these problems in a simple and effective way. Our method is based on the idea to simulate the effect of the cables instead of the cables themselves. To this end we propose a new special type of distance constraint we call cable joint that changes both its attachment points and its rest length dynamically. A cable connecting a series of objects is then modeled as a sequence of cable joints which reduces the complexity of the simulation from the order of the number of segments to just the number of connected objects. This makes simulations both faster and more robust as we will demonstrate on a variety of examples.

Cable Joints

Distributing and Load Balancing Sparse Fluid Simulations

Chinmayee Shah, David Hyde, Hang Qu, and Philip Levis

This paper describes a general algorithm and a system for load balancing sparse fluid simulations. Automatically distributing sparse fluid simulations efficiently is challenging because the computational load varies across the simulation domain and time. A key challenge with load balancing is that optimal decision making requires knowing the fluid distribution across partitions for future time steps, but computing this state for an arbitrary simulation requires running the simulation itself. The key insight of this paper is that it is possible to predict future load by running a low resolution simulation in parallel. This paper describes a system design and techniques to automatically distribute and load balance sparse fluid simulations, and presents speculative load balancing, a general technique to effectively balance future load using information about future load distribution obtained via a low resolution simulation. We mathematically formulate the problem of load balancing over multiple time-steps and present a polynomial time algorithm to compute an approximate solution to it. Our experimental results show that distributing and speculatively load balancing sparse FLIP simulations over 8 nodes speeds them up by 5.5x to 7.2x, and that speculative load balancing generates assignments that perform within 20% of optimal.

Distributing and Load Balancing Sparse Fluid Simulations

Fast Corotated FEM using Operator Splitting

Tassilo Kugelstadt, Dan Koschier, Jan Bender

In this paper we present a novel operator splitting approach for corotated FEM simulations. The deformation energy of the corotated linear material model consists of two additive terms. The first term models stretching in the individual spatial directions and the second term describes resistance to volume changes. By formulating the backward Euler time integration scheme as an optimization problem, we show that the first term is invariant to rotations. This allows us to use an operator splitting approach and to solve both terms individually with different numerical methods. The stretching part is solved accurately with an optimization integrator, which can be done very efficiently because the system matrix is constant over time such that its Cholesky factorization can be precomputed. The volume term is solved approximately by using the compliant constraints method and Gauss-Seidel iterations. Further, we introduce the analytic polar decomposition which allows us to speed up the extraction of the rotational part of the deformation gradient and to recover inverted elements. Finally, this results in an extremely fast and robust simulation method with high visual quality that outperforms standard corotated FEMs by more than two orders of magnitude and even the fast but inaccurate PBD and shape matching methods by more than one order of magnitude without having their typical drawbacks. This enables a very efficient simulation of complex scenes containing more than a million elements.

Fast Corotated FEM using Operator Splitting

An Extended Partitioned Method for Conservative Solid-Fluid Coupling

Muzaffer Akbay, Nicholas Nobles, Victor Zoran, Tamar Shinar

We present a novel extended partitioned method for two-way solid-fluid coupling, where the fluid and solid solvers are treated as black boxes with limited exposed interfaces, facilitating modularity and code reusability. Our method achieves improved stability and extended range of applicability over standard partitioned approaches through three techniques. First, we couple the black-box solvers through a small, reduced-order monolithic system, which is constructed on the fly from input/output pairs generated by the solid and fluid solvers. Second, we use a conservative, impulse-based interaction term to couple the solid and fluid rather than typical pressure-based forces. We show that both of these techniques significantly improve stability and reduce the number of iterations needed for convergence. Finally, we propose a novel boundary pressure projection method that allows for the partitioned simulation of a fully enclosed fluid coupled to a dynamic solid, a scenario that has been problematic for partitioned methods. We demonstrate the benefits of our extended partitioned method by coupling Eulerian fluid solvers for smoke and water to Lagrangian solid solvers for volumetric and thin deformable and rigid objects in a variety of challenging scenarios. We further demonstrate our method by coupling a Lagrangian SPH fluid solver to a rigid body solver

An Extended Partitioned Method for Conservative Solid-Fluid Coupling

HairControl: A Tracking Solution for Directable Hair Simulation

Antoine Milliez, Bob Sumner, Markus Gross, Bernhard Thomaszewski

We present a method for adding artistic control to physics-based hair simulation. Taking as input an animation of a coarse set of guide hairs, we constrain a subsequent higher-resolution simulation of detail hairs to follow the input motion in a spatially-averaged sense. The resulting high-resolution motion adheres to the artistic intent but is enhanced with detailed deformations and dynamics generated by physics-based simulation. The technical core of our approach is formed by a set of tracking constraints, requiring the center of mass of a given subset of detail hair to maintain its position relative to a reference point on the corresponding guide hair. As a crucial element of our formulation, we introduce the concept of dynamically changing constraint targets that allow reference points to slide along the guide hairs to provide sufficient flexibility for natural deformations. We furthermore propose to regularize the null space of the tracking constraints based on variance minimization, effectively controlling the amount of spread in the hair. We demonstrate the ability of our tracking solver to generate directable yet natural hair motion on a set of targeted experiments and show its application to production-level animations.

HairControl: A Tracking Solution for Directable Hair Simulation

Projective peridynamics for modeling versatile elastoplastic materials

Xiaowei He, Huamin Wang, Enhua Wu

Unified simulation of versatile elastoplastic materials and different dimensions offers many advantages in animation production, contact handling, and hardware acceleration. The unstructured particle representation is particularly suitable for this task, thanks to its simplicity. However, previous meshless techniques either need too much computational cost for addressing stability issues, or lack physical meanings and fail to generate interesting deformation behaviors, such as the Poisson effect. In this paper, we study the development of an elastoplastic model under the state-based peridynamics framework, which uses integrals rather than partial derivatives in its formulation. To model elasticity, we propose a unique constitutive model and an efficient iterative simulator solved in a projective dynamics way. To handle plastic behaviors, we incorporate our simulator with the Drucker-Prager yield criterion and a reference position update scheme, both of which are implemented under peridynamics. Finally, we show how to strengthen the simulator by position-based constraints and spatially varying stiffness models, to achieve incompressibility, particle redistribution, cohesion, and friction effects in viscoelastic and granular flows. Our experiments demonstrate that our unified, meshless simulator is flexible, efficient, robust, and friendly with parallel computing.

Projective peridynamics for modeling versatile elastoplastic materials

Cosserat Rods with Projective Dynamics

Carlota Soler, Tobias Martin, Olga Sorkine-Hornung

We present a novel method to simulate Cosserat rods with Projective Dynamics (PD). The proposed method is both numerically robust and accurate with respect to the underlying physics, making it suitable for a variety of applications in computer graphics and related disciplines. Cosserat theory assigns an orientation frame to each point and is thus able to realistically simulate stretching and shearing effects, in addition to bending and twisting. Within the PD framework, it is possible to obtain accurate simulations given the implicit integration over time and its decoupling of the local-global solve. In the proposed method, we start from the continuous formulation of the Cosserat theory and derive the constraints for the PD solver. We extend the standard definition of PD and add body orientations as system variables. Thus, we include the preservation of angular momentum, so tha twisting and bending can be accurately simulated. Our formulation allows the simulation of different bending behaviors with respect to a user-defined Young’s modulus, the radius of the rod’s cross-section, and material density. We show how different material specifications in our simulations converge within a few iterations to a reference solution, generated with a high-precision finite element method. Furthermore, we demonstrate mesh independence of our formulation: Refining the simulation mesh still results in the same characteristic motion, which is in contrast to previous position based method.

A Temporally Adaptive Material Point Method with Regional Time Stepping

Yu Fang, Yuanming Hu, Shi-Min Hu, Chenfanfu Jiang

Spatially and temporally adaptive algorithms can substantially improve the computational efficiency of many numerical schemes in computational mechanics and physics-based animation. Recently, a crucial need for temporal adaptivity in the Material Point Method (MPM) is emerging due to the potentially substantial variation of material stiffness and velocities in multi-material scenes. In this work, we propose a novel temporally adaptive symplectic Euler scheme for MPM with regional time stepping (RTS), where different time steps are used in different regions. We design a time stepping scheduler operating at the granularity of small blocks to maintain a natural consistency with the hybrid particle/grid nature of MPM. Our method utilizes the Sparse Paged Grid (SPGrid) data structure and simultaneously offers high efficiency and notable ease of implementation with a practical multi-threaded particle-grid transfer strategy. We demonstrate the efficacy of our asynchronous MPM method on various examples including elastic objects, granular media, and fluid.

A Temporally Adaptive Material Point Method with Regional Time Stepping

Time-Domain Parallelization for Accelerating Cloth Simulation

Junbang Liang, Ming C. Lin

Cloth simulations, widely used in computer animation and apparel design, can be computationally expensive for real-time applications. Some parallelization techniques have been proposed for visual simulation of cloth using CPU or GPU clusters and often rely on parallelization using spatial domain decomposition techniques that have a large communication overhead. In this paper, we propose a novel time-domain parallelization technique that makes use of the two-level mesh representation to resolve the time-dependency issue and develop a practical algorithm to smooth the state transition from the corresponding coarse to fine meshes. A load estimation and a load balancing technique used in online partitioning are also proposed to maximize the performance acceleration. Our method achieves a nearly linear performance scaling on manycore clusters and outperforms spatial-domain parallelization on a diverse set of benchmarks.

Time-Domain Parallelization for Accelerating Cloth Simulation