Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
skip to main content
research-article
Public Access

Decomposed optimization time integrator for large-step elastodynamics

Published: 12 July 2019 Publication History

Abstract

Simulation methods are rapidly advancing the accuracy, consistency and controllability of elastodynamic modeling and animation. Critical to these advances, we require efficient time step solvers that reliably solve all implicit time integration problems for elastica. While available time step solvers succeed admirably in some regimes, they become impractically slow, inaccurate, unstable, or even divergent in others --- as we show here. Towards addressing these needs we present the Decomposed Optimization Time Integrator (DOT), a new domain-decomposed optimization method for solving the per time step, nonlinear problems of implicit numerical time integration. DOT is especially suitable for large time step simulations of deformable bodies with nonlinear materials and high-speed dynamics. It is efficient, automated, and robust at large, fixed-size time steps, thus ensuring stable, continued progress of high-quality simulation output. Across a broad range of extreme and mild deformation dynamics, using frame-rate size time steps with widely varying object shapes and mesh resolutions, we show that DOT always converges to user-set tolerances, generally well-exceeding and always close to the best wall-clock times across all previous nonlinear time step solvers, irrespective of the deformation applied.

Supplementary Material

ZIP File (a70-li.zip)
Supplemental material
MP4 File (papers_172.mp4)

References

[1]
A. Abdelfattah, A. Haidar, S. Tomov, and J. Dongarra. 2017. Fast Cholesky factorization on GPUs for batch and native modes in MAGMA. J of Comp Sci 20 (2017).
[2]
U. M Ascher. 2008. Numerical methods for evolutionary differential equations.
[3]
S. Bouaziz, S. Martin, T. Liu, L. Kavan, and M. Pauly. 2014. Projective Dynamics: Fusing Constraint Projections for Fast Simulation. ACM Trans Graph 33, 4 (2014).
[4]
S. Boyd, N. Parikh, E. Chu, B. Peleato, J. Eckstein, et al. 2011. Distributed optimization and statistical learning via the alternating direction method of multipliers. Foundations and Trends® in Machine learning 3, 1 (2011).
[5]
J. Brown and P. Brune. 2013. Low-rank quasi-Newton updates for robust Jacobian lagging in Newton-type methods. In Int Conf Math Comp Meth App Nucl Sci Eng.
[6]
J. C. Butcher. 2016. Numerical methods for ordinary differential equations.
[7]
I. Chao, U. Pinkall, P. Sanan, and P. Schröder. 2010. A simple geometric model for elastic deformations. ACM Trans Graph (SIGGRAPH) 29, 4 (2010).
[8]
Y. Chen, T. A. Davis, W. W. Hager, and S. Rajamanickam. 2008. Algorithm 887: CHOLMOD, Supernodal Sparse Cholesky Factorization and Update/Downdate. ACM Trans on Mathematical Software (TOMS) 35, 3 (2008).
[9]
P. Deufihard. 2011. Newton Methods for Nonlinear Problems: Affine Invariance and Adaptive Algorithms.
[10]
V. Dolean, P. Jolivet, and F. Nataf. 2015. An introduction to domain decomposition methods: algorithms, theory and parallel implementation.
[11]
T. Gast, C. Schroeder, A. Stomakhin, C. Jiang, and J. M Teran. 2015. Optimization integrator for large time steps. IEEE Trans Vis Comp Graph 21, 10 (2015).
[12]
E. Hairer, C. Lubich, and G. Wanner. 2006. Geometric Numerical Integration.
[13]
E. Hairer, S. P Nørsett, and G. Wanner. 2008. Solving Ordinary Differential Equations I.
[14]
E. Hairer and G. Wanner. 1996. Solving Ordinary Differential Equations II.
[15]
F. Hecht, Y.J. Lee, J. R. Shewchuk, and J. F. O'Brien. 2012. Updated Sparse Cholesky Factors for Corotational Elastodynamics. ACM Trans Graph 31, 5 (2012).
[16]
J. Huang, X. Liu, H. Bao, B. Guo, and H. Shum. 2006. An efficient large deformation method using domain decomposition. Comp & Graph 30, 6 (2006).
[17]
C. Kane, J. E Marsden, M. Ortiz, and M. West. 2000. Variational integrators and the Newmark algorithm for conservative and dissipative mechanical systems. Int J for Numer Meth in Eng 49, 10 (2000).
[18]
G. Karypis and V. Kumar. 1998. A fast and high quality multilevel scheme for partitioning irregular graphs. SIAM J on Sci Comp 20 (1998).
[19]
L. Kharevych, W. Yang, Y. Tong, E. Kanso, J. E Marsden, P. Schröder, and M. Desbrun. 2006. Geometric, variational integrators for computer animation. In Symp Comp Anim.
[20]
T. Kim and D. L James. 2012. Physics-based character skinning using multidomain subspace deformations. IEEE Trans on visualization and Comp Graph 18, 8 (2012).
[21]
H. Liu, N. Mitchell, M. Aanjaneya, and E. Sifakis. 2016. A scalable schur-complement fluids solver for heterogeneous compute platforms. ACM Trans Graph 35, 6 (2016).
[22]
T. Liu, A. W. Bargteil, J. F. O'Brien, and L. Kavan. 2013. Fast Simulation of Mass-Spring Systems. ACM Trans Graph 32, 6 (2013). Proc of ACM SIGGRAPH Asia.
[23]
T. Liu, S. Bouaziz, and L. Kavan. 2017. Quasi-Newton Methods for Real-Time Simulation of Hyperelastic Materials. ACM Trans Graph 36, 4 (2017).
[24]
M. Macklin, M. Müller, and N. Chentanez. 2016. XPBD: position-based simulation of compliant constrained dynamics. In Proc of the 9th Int Conf on Motion in Games.
[25]
S. Martin, B. Thomaszewski, E. Grinspun, and M. Gross. 2011. Example-based elastic materials. ACM Trans Graph (SIGGRAPH) 30, 4 (2011).
[26]
A. McAdams, A. Selle, R. Tamstorf, J. Teran, and E. Sifakis. 2011. Computing the singular value decomposition of 3X 3 matrices with minimal branching and elementary floating point operations. University of Wisconsin Madison (2011).
[27]
M. Müller, B. Heidelberger, M. Hennix, and J. Ratcliff. 2007. Position based dynamics. J. Vis. Commun. Imag Represent. 18, 2 (2007).
[28]
R. Narain, M. Overby, and G. E Brown. 2016. ADMM ⊇ projective dynamics: fast simulation of general constitutive models. In Symp on Comp Anim.
[29]
JW Neuberger. 1985. Steepest descent and differential equations. J of the Mathematical Society of Japan 37, 2 (1985).
[30]
J. Nocedal and S. Wright. 2006. Numerical Optimization.
[31]
M. Ortiz and L. Stainier. 1999. The variational formulation of viscoplastic constitutive updates. Comp Meth in App Mech and Eng 171, 3--4 (1999).
[32]
M. Overby, G. E Brown, J. Li, and R. Narain. 2017. ADMM ⊇ Projective Dynamics: Fast Simulation of Hyperelastic Models with Dynamic Constraints. IEEE Trans Vis Comp Graph 23, 10 (2017).
[33]
N. Parikh and S. Boyd. 2012. Block splitting for distributed optimization.
[34]
A. Quarteroni, A. Valli, and P.M.A. Valli. 1999. Domain Decomposition Methods for Partial Differential Equations.
[35]
T. Schneider, Y. Hu, J. Dumas, X. Gao, D. Panozzo, and D. Zorin. 2018. Decoupling simulation accuracy from mesh quality. ACM Trans Graph (2018).
[36]
S. Sellán, H. Y. Cheng, Y. Ma, M. Dembowski, and A. Jacobson. 2018. Solid Geometry Processing on Deconstructed Domains. CoRR (2018).
[37]
VE Shamanskii. 1967. A modification of Newton's method. Ukrainian Mathematical J 19, 1 (1967).
[38]
A. Shtengel, R. Poranne, O. Sorkine-Hornung, S. Z. Kovalsky, and Y. Lipman. 2017. Geometric Optimization via Composite Majorization. ACM Trans Graph 36, 4 (2017).
[39]
B. Smith, F. De Goes, and T. Kim. 2018. Stable Neo-Hookean Flesh Simulation. ACM Trans Graph 37, 2 (2018).
[40]
A. Stomakhin, R. Howes, C. Schroeder, and J. M Teran. 2012. Energetically consistent invertible elasticity. In Symp Comp Anim.
[41]
J. Teran, E. Sifakis, G. Irving, and R. Fedkiw. 2005. Robust Quasistatic Finite Elements and Flesh Simulation. In Symp Comp Anim.
[42]
C. Xiao-Chuan and D. Maksymilian. 1994. Domain Decomposition Methods for Monotone Nonlinear Elliptic Problems. In Contemporary Math.
[43]
Y. Zhu, R. Bridson, and D. M. Kaufman. 2018. Blended Cured Quasi-Newton for Distortion Optimization. ACM Trans. on Graph (SIGGRAPH) (2018).

Cited By

View all

Index Terms

  1. Decomposed optimization time integrator for large-step elastodynamics

    Recommendations

    Comments

    Information & Contributors

    Information

    Published In

    cover image ACM Transactions on Graphics
    ACM Transactions on Graphics  Volume 38, Issue 4
    August 2019
    1480 pages
    ISSN:0730-0301
    EISSN:1557-7368
    DOI:10.1145/3306346
    Issue’s Table of Contents
    Permission to make digital or hard copies of all or part of this work for personal or classroom use is granted without fee provided that copies are not made or distributed for profit or commercial advantage and that copies bear this notice and the full citation on the first page. Copyrights for components of this work owned by others than ACM must be honored. Abstracting with credit is permitted. To copy otherwise, or republish, to post on servers or to redistribute to lists, requires prior specific permission and/or a fee. Request permissions from [email protected]

    Publisher

    Association for Computing Machinery

    New York, NY, United States

    Publication History

    Published: 12 July 2019
    Published in TOG Volume 38, Issue 4

    Permissions

    Request permissions for this article.

    Check for updates

    Author Tags

    1. computational optimization
    2. domain decomposition

    Qualifiers

    • Research-article

    Funding Sources

    Contributors

    Other Metrics

    Bibliometrics & Citations

    Bibliometrics

    Article Metrics

    • Downloads (Last 12 months)243
    • Downloads (Last 6 weeks)67
    Reflects downloads up to 15 Oct 2024

    Other Metrics

    Citations

    Cited By

    View all
    • (2024)Vertex Block DescentACM Transactions on Graphics10.1145/365817943:4(1-16)Online publication date: 19-Jul-2024
    • (2024)Position-Based Nonlinear Gauss-Seidel for Quasistatic HyperelasticityACM Transactions on Graphics10.1145/365815443:4(1-15)Online publication date: 19-Jul-2024
    • (2024)Super-Resolution Cloth Animation with Spatial and Temporal CoherenceACM Transactions on Graphics10.1145/365814343:4(1-14)Online publication date: 19-Jul-2024
    • (2024)Bistable Aerial Transformer: A Quadrotor Fixed-Wing Hybrid That Morphs Dynamically Via Passive Soft MechanismJournal of Mechanisms and Robotics10.1115/1.406515916:7Online publication date: 24-Apr-2024
    • (2024)Augmented Incremental Potential Contact for Sticky InteractionsIEEE Transactions on Visualization and Computer Graphics10.1109/TVCG.2023.329565630:8(5596-5608)Online publication date: Aug-2024
    • (2024)A Reconfigurable Data Glove for Reconstructing Physical and Virtual GraspsEngineering10.1016/j.eng.2023.01.00932(202-216)Online publication date: Jan-2024
    • (2024)Primal residual reduction with extended position based dynamics and hyperelasticityComputers and Graphics10.1016/j.cag.2024.103902119:COnline publication date: 1-Apr-2024
    • (2023)Primal Extended Position Based Dynamics for HyperelasticityProceedings of the 16th ACM SIGGRAPH Conference on Motion, Interaction and Games10.1145/3623264.3624437(1-10)Online publication date: 15-Nov-2023
    • (2023)Power Plastics: A Hybrid Lagrangian/Eulerian Solver for Mesoscale Inelastic FlowsACM Transactions on Graphics10.1145/361834442:6(1-11)Online publication date: 5-Dec-2023
    • (2023)Second-order Stencil Descent for Interior-point HyperelasticityACM Transactions on Graphics10.1145/359210442:4(1-16)Online publication date: 26-Jul-2023
    • Show More Cited By

    View Options

    View options

    PDF

    View or Download as a PDF file.

    PDF

    eReader

    View online with eReader.

    eReader

    Get Access

    Login options

    Full Access

    Media

    Figures

    Other

    Tables

    Share

    Share

    Share this Publication link

    Share on social media