Increasing the critical time step: micro-inertia, inertia penalties and mass scaling

Explicit time integration is a popular method to simulate the dynamical behaviour of a system. Unfortunately, explicit time integration is only conditionally stable: the time step must be chosen not larger than the so-called “critical time step”, otherwise the numerical solution may become unstable. To reduce the CPU time needed to carry out simulations, it is desirable to explore methods that increase the critical time step, which is the main objective of our paper. To do this, first we discuss and compare three approaches to increase the critical time step: micro-inertia formulations from continuum mechanics, inertia penalties which are used in computational mechanics, and mass scaling techniques that are mainly used in structural dynamics. As it turns out, the similarities between these methods are significant, and in fact they are identical in 1D if linear finite elements are used. This facilitates interpretation of the additional parameters in the various methods. Next, we derive, for a few simple finite element types, closed-form expressions for the critical time step with micro-structural magnification factors. Finally, we discuss computational overheads and some implementational details.

[1]  J. Huetink,et al.  Simulation of aluminium sheet forming at elevated temperatures , 2006 .

[2]  A. Eringen.  On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves , 1983 .

[3]  J. Awrejcewicz,et al.  Asymptotic approaches in mechanics: New parameters and procedures , 2003 .

[4]  Antonio Rodríguez-Ferran,et al.  Bipenalty method for time domain computational dynamics , 2010, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[5]  H. Askes,et al.  Element size and time step selection procedures for the numerical analysis of elasticity with higher-order inertia , 2008 .

[6]  R. D. Mindlin.  Micro-structure in linear elasticity , 1964 .

[7]  Jacob Fish,et al.  Non‐local dispersive model for wave propagation in heterogeneous media: one‐dimensional case , 2002 .

[8]  Harm Askes,et al.  The Representative Volume Size in Static and Dynamic Micro-Macro Transitions , 2005 .

[9]  Julius Kaplunov,et al.  On Timoshenko-Reissner type theories of plates and shells , 1993 .

[10]  Mattias Unosson,et al.  Selective mass scaling for explicit finite element analyses , 2005 .

[11]  Eiji Nakamachi,et al.  Dynamic‐explicit elastic plastic finite‐element simulation of hemispherical punch‐drawing of sheet metal , 1996 .

[12]  Lars Olovsson,et al.  Iterative solution technique in selective mass scaling , 2005 .

[13]  M. B. Rubin,et al.  Continuum model of dispersion caused by an inherent material characteristic length , 1995 .

[14]  Jeong Kim,et al.  A comparative study of implicit and explicit FEM for the wrinkling prediction in the hydroforming process , 2003 .

[15]  B. H. Aubert,et al.  A mass penalty technique to control the critical time increment in explicit dynamic finite element analyses , 1995 .

[16]  Shaoqun Zeng,et al.  The key technology and realization of virtual ring rolling , 2007 .

[17]  A. Tyas,et al.  Asymptotic equivalence of homogenisation procedures and fine-tuning of continuum theories , 2008 .

[18]  J. Fish,et al.  A Dispersive Model for Wave Propagation in Periodic Heterogeneous Media Based on Homogenization With Multiple Spatial and Temporal Scales , 2001 .

[19]  C. Sun,et al.  Modeling micro-inertia in heterogeneous materials under dynamic loading , 2002 .

[20]  Graeme W Milton,et al.  On modifications of Newton's second law and linear continuum elastodynamics , 2007, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[21]  J. Awrejcewicz,et al.  Continuous models for 2D discrete media valid for higher-frequency domain , 2008 .

[22]  Petr Plechác,et al.  Implicit Mass-matrix Penalization of Hamiltonian Dynamics with Application to Exact Sampling of Stiff Systems , 2009, Multiscale Model. Simul..

[23]  Andrei V. Metrikine,et al.  An isotropic dynamically consistent gradient elasticity model derived from a 2D lattice , 2006 .

[24]  J. Awrejcewicz,et al.  Continuous models for 1D discrete media valid for higher-frequency domain , 2005 .

[25]  H. Askes,et al.  Penalty methods for time domain computational dynamics based on positive and negative inertia , 2009 .

[26]  Mattias Unosson,et al.  Selective mass scaling for thin walled structures modeled with tri-linear solid elements , 2004 .

[27]  J. Engelbrecht,et al.  Waves in microstructured materials and dispersion , 2005 .

[28]  H. Askes,et al.  One-dimensional dynamically consistent gradient elasticity models derived from a discrete microstructure: Part 2: Static and dynamic response , 2002 .

[29]  Andrei V. Metrikine,et al.  One-dimensional dynamically consistent gradient elasticity models derived from a discrete microstructure: Part 1: Generic formulation , 2002 .

[30]  Ted Belytschko,et al.  On the dynamic effects of explicit FEM in sheet metal forming analysis , 1998 .

[31]  E. Aifantis,et al.  On the role of microstructure in the behavior of soils: Effects of higher order gradients and internal inertia , 1994 .