Kugelstadt, TassiloBender, JanFernández-Fernández, José AntonioJeske, Stefan RhysLöschner, FabianLongva, AndreasNarain, Rahul and Neff, Michael and Zordan, Victor2022-02-072022-02-0720212577-6193https://doi.org/10.1145/3480142https://diglib.eg.org:443/handle/10.1145/3480142We 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.Computing methodologiesPhysical simulationSmoothed Particle Hydrodynamicsdeformable solidsfluid simulationsolidfluid coupling10.1145/3480142