Browsing by Author "Kugelstadt, Tassilo"

Now showing 1 - 3 of 3
  • No Thumbnail Available
    Item
    Direct Position‐Based Solver for Stiff Rods
    (© 2018 The Eurographics Association and John Wiley & Sons Ltd., 2018) Deul, Crispin; Kugelstadt, Tassilo; Weiler, Marcel; Bender, Jan; Chen, Min and Benes, Bedrich
    In this paper, we present a novel direct solver for the efficient simulation of stiff, inextensible elastic rods within the position‐based dynamics (PBD) framework. It is based on the XPBD algorithm, which extends PBD to simulate elastic objects with physically meaningful material parameters. XPBD approximates an implicit Euler integration and solves the system of non‐linear equations using a non‐linear Gauss–Seidel solver. However, this solver requires many iterations to converge for complex models and if convergence is not reached, the material becomes too soft. In contrast, we use Newton iterations in combination with our direct solver to solve the non‐linear equations which significantly improves convergence by solving all constraints of an acyclic structure (tree), simultaneously. Our solver only requires a few Newton iterations to achieve high stiffness and inextensibility. We model inextensible rods and trees using rigid segments connected by constraints. Bending and twisting constraints are derived from the well‐established Cosserat model. The high performance of our solver is demonstrated in highly realistic simulations of rods consisting of multiple 10 000 segments. In summary, our method allows the efficient simulation of stiff rods in the PBD framework with a speedup of two orders of magnitude compared to the original XPBD approach.We present a novel direct solver for the efficient simulation of stiff, inextensible elastic rods. It is based on the XPBD algorithm, which extends Position‐Based Dynamics to simulate elastic objects with physically meaningful material parameters. However, the non‐linear Gauss‐Seidel solver of XPBD requires many iterations to converge for complex models and if convergence is not reached, the material becomes too soft. In contrast, we use Newton iterations in combination with our direct solver which significantly improves convergence by solving all constraints of an acyclic structure simultaneously. We model rods using rigid segments connected by constraints. Bending and twisting constraints are derived from the Cosserat model. The high performance of our solver allows the simulation of rods consisting of multiple 10 000 segments with a speedup of two orders of magnitude compared to the original XPBD approach.
  • Thumbnail Image
    Item
    Fast Corotated Elastic SPH Solids with Implicit Zero-Energy Mode Control
    (ACM, 2021) Kugelstadt, Tassilo; Bender, Jan; Fernández-Fernández, José Antonio; Jeske, Stefan Rhys; Löschner, Fabian; Longva, Andreas; Narain, Rahul and Neff, Michael and Zordan, Victor
    We develop a new operator splitting formulation for the simulation of corotated linearly elastic solids with Smoothed Particle Hydrodynamics (SPH). Based on the technique of Kugelstadt et al. [2018] originally developed for the Finite Element Method (FEM), we split the elastic energy into two separate terms corresponding to stretching and volume conservation, and based on this principle, we design a splitting scheme compatible with SPH. The operator splitting scheme enables us to treat the two terms separately, and because the stretching forces lead to a stiffness matrix that is constant in time, we are able to prefactor the system matrix for the implicit integration step. Solid-solid contact and fluid-solid interaction is achieved through a unified pressure solve. We demonstrate more than an order of magnitude improvement in computation time compared to a state-of-the-art SPH simulator for elastic solids. We further improve the stability and reliability of the simulation through several additional contributions. We introduce a new implicit penalty mechanism that suppresses zero-energy modes inherent in the SPH formulation for elastic solids, and present a new, physics-inspired sampling algorithm for generating highquality particle distributions for the rest shape of an elastic solid. We finally also devise an efficient method for interpolating vertex positions of a high-resolution surface mesh based on the SPH particle positions for use in high-fidelity visualization.
  • Thumbnail Image
    Item
    Higher-Order Time Integration for Deformable Solids
    (The Eurographics Association and John Wiley & Sons Ltd., 2020) Löschner, Fabian; Longva, Andreas; Jeske, Stefan; Kugelstadt, Tassilo; Bender, Jan; Bender, Jan and Popa, Tiberiu
    Visually appealing and vivid simulations of deformable solids represent an important aspect of physically based computer animation. For the temporal discretization, it is customary in computer animation to use first-order accurate integration methods, such as Backward Euler, due to their simplicity and robustness. Although there is notable research on second-order methods, their use is not widespread. Many of these well-known methods have significant drawbacks such as severe numerical damping or scene-dependent time step restrictions to ensure stability. In this paper, we discuss the most relevant requirements on such methods in computer animation and motivate the interest beyond first-order accuracy. Keeping these requirements in mind, we investigate several promising methods from the families of diagonally implicit Runge-Kutta (DIRK) and Rosenbrock methods which currently do not appear to have considerable popularity in this field. We show that the usage of such methods improves the visual quality of physical animations. In addition, we demonstrate that they allow distinctly more control over damping at lower computational cost than classical methods. As part of our theoretical contribution, we review aspects of simulations that are often considered more intricate with higher-order methods, such as contact handling. To this end, we derive an implicit linearized contact model based on a predictor-corrector approach that leads to consistent behavior with higher-order integrators as predictors. Our contact model is well suited for the simulation of stiff, nonlinear materials with the integration methods presented in this paper and more common methods such as Backward Euler alike.