An energy-preserving Discrete Element Method for elastodynamics

Monasse, Laurent; Mariotti, Christian

ESAIM: Mathematical Modelling and Numerical Analysis , Tome 46 (2012) no. 6, pp. 1527-1553.

We develop a Discrete Element Method (DEM) for elastodynamics using polyhedral elements. We show that for a given choice of forces and torques, we recover the equations of linear elastodynamics in small deformations. Furthermore, the torques and forces derive from a potential energy, and thus the global equation is an Hamiltonian dynamics. The use of an explicit symplectic time integration scheme allows us to recover conservation of energy, and thus stability over long time simulations. These theoretical results are illustrated by numerical simulations of test cases involving large displacements.

[1] H.C. Andersen, RATTLE : A “velocity” version of the SHAKE algorithm for molecular dynamics calculations. J. Comput. Phys. 52 (1983) 24-34. | Zbl

[2] C. Antoci, M. Gallati and S. Sibilla, Numerical simulation of fluid-structure interaction by SPH. 4th MIT Conference on Computational Fluid and Solid Mechanics. Comput. Struct. 85 (2007) 879-890.

[3] J. Bonet and T.S.L. Lok, Variational and momentum preservation aspects of Smooth particle hydrodynamic formulations. Comput. Meth. Appl. Mech. Eng. 180 (1999) 97-115. | MR | Zbl

[4] P.A. Cundall and O.D.L. Strack, A discrete numerical model for granular assemblies. Geotech. 29 (1979) 47-65.

[5] G.A. D'Addetta, F. Kun and E. Ramm, On the application of a discrete model to the fracture process of cohesive granular materials. Granul. Matter 4 (2002) 77-90. | Zbl

[6] A.T. De Hoop, A modification of Cagniard's method for solving seismic pulse problem. Appl. Sci. Res. B 8 (1960) 349-356. | Zbl

[7] A.C. Eringen, Theory of micropolar elasticity, in Fracture, edited by H. Liebowitz. Academic Press, New York 2 (1968) 621-729. | Zbl

[8] E.P. Fahrenthold and B.A. Horban, An improved hybrid particle-element method for hypervelocity impact simulation. Symposium on Hypervelocity Impact, Galveston. Texas (2000). Int. J. Impact Eng. 26 (2001) 169-178.

[9] E.P. Fahrenthold and R. Shivarama, Extension and validation of a hybrid particle-finite element method for hypervelocity impact simulation. Hypervelocity Impact Symposium. Int. J. Impact Eng. 29 (2003) 237-246.

[10] Y.T. Feng, K. Han, C.F. Li and D.R.J. Owen, Discrete thermal element modelling of heat conduction in particle systems : Basic formulations. J. Comput. Phys. 227 (2008) 5072-5089. | Zbl

[11] S. Forest, F. Pradel and K. Sab, Asymptotic analysis of heterogeneous Cosserat media. Int. J. Solids Struct. 38 (2001) 4585-4608. | MR | Zbl

[12] R.A. Gingold and J.J. Monaghan, smoothed particle hydrodynamics : Theory and application to nonspherical stars. Mon. Not. R. Astron. Soc. 181 (1977) 375-389. | Zbl

[13] O. Gonzalez, Exact energy and momentum conserving algorithms for general models in nonlinear elasticity. Comput. Meth. Appl. Mech. Eng. 190 (2000) 1763-1783. | MR | Zbl

[14] E. Hairer and G. Vilmart, Preprocessed discrete Moser-Veselov algorithm for the full dynamics of a rigid body. J. Phys. A 39 (2006) 13225-13235. | MR | Zbl

[15] E. Hairer, C. Lubich and G. Wanner, Geometric Numerical Integration : Structure-Preserving Algorithms for Ordinary Differential Equations, 2nd edition. Springer Series in Comput. Math. 31 (2006). | MR | Zbl

[16] K. Han, Y.T. Feng and D.R.J. Owen, Coupled lattice Boltzmann and discrete element modelling of fluid-particle interaction problems, in 4th MIT Conference on Computational Fluid and Solid Mechanics. Comput. Struct. 85 (2007) 1080-1088.

[17] P. Hauret and P. Le Tallec, Energy-controlling time integration methods for nonlinear elastodynamics and low-velocity impact. Comput. Meth. Appl. Mech. Eng. 195 (2006) 4890-4916. | MR | Zbl

[18] D.L. Hicks, J.W. Swegle and S.W. Attaway, Conservative smoothing stabilizes discrete-numerical instabilities in SPH material dynamics computations. Appl. Math. Comput. 85 (1997) 209-226. | MR | Zbl

[19] W.G. Hoover, Smooth Particle Applied Mechanics : The State of the Art (World Scientific). Adv. Ser. Nonlinear Dyn. 25 (2006). | MR | Zbl

[20] W.G. Hoover, W.T. Arhurst and R.J. Olness, Two-dimensional studies of crystal stability and fluid viscosity. J. Chem. Phys. 60 (1974) 4043-4047.

[21] A. Ibrahimbegovic and A. Delaplace, Microscale and mesoscale discrete models for dynamic fracture of structures built of brittle material. Comput. Struct. 81 (2003) 1255-1265.

[22] J.C. Koo and E.P. Fahrenthold, Discrete Hamilton's equations for arbitrary Lagrangian-Eulerian dynamics of viscous compressible flow. Comput. Meth. Appl. Mech. Eng. 189 (2000) 875-900. | Zbl

[23] S. Koshizuka and Y. Oka, Moving-particle semi-implicit method for fragmentation of incompressible fluid. Nucl. Sci. Eng. 123 (1996) 421-434.

[24] S. Koshizuka, A. Nobe and Y. Oka, Numerical analysis of breaking waves using the moving particle semi-implicit method. Int. J. Numer. Meth. Fluids 26 (1998) 751-769. | Zbl

[25] S. Koshizuka, M.S. Song and Y. Oka, A particle method for three-dimensional elastic analysis, in Proc. of 6th World Cong. Computational Mechanics (WCCM VI). Beijing (2004).

[26] F. Kun and H. Herrmann, A study of fragmentation processes using a discrete element method. Comput. Meth. Appl. Mech. Eng. 138 (1996) 3-18. | Zbl

[27] H. Lamb, On the propagation of tremors over the surface of an elastic solid. Philos. Trans. R. Soc. Lond. A 203 (1904) 1-42. | JFM

[28] T.A. Laursen and X.N. Meng, A new solution procedure for application of energy-conserving algorithms to general constitutive models in nonlinear elastodynamics. Comput. Meth. Appl. Mech. Eng. 190 (2001) 6309-6322. | MR | Zbl

[29] C.J.K. Lee, H. Noguchi and S. Koshizuka, Fluid-shell structure interaction analysis by coupled particle and finite element method, in 4th MIT Conference on Computational Fluid and Solid Mechanics. Comput. Struct. 85 (2007) 688-697.

[30] B.J. Leimkuhler and R.D. Skeel, Symplectic numerical integrators in constrained Hamiltonian systems. J. Comput. Phys. 112 (1994) 117-125. | MR | Zbl

[31] A. Lew, J.E. Marsden, M. Ortiz and M. West, Variational time integrators. Int. J. Numer. Meth. Eng. 60 (2004) 153-212. | MR | Zbl

[32] L.D. Libersky, A.G. Petschek, T.C. Carney, J.R. Hipp and F.A. Allahdadi, High strain Lagrangian hydrodynamics : A three-dimensional SPH code for dynamic material response. J. Comput. Phys. 109 (1993) 76-83. | Zbl

[33] L.B. Lucy, A numerical approach to the testing of the fission hypothesis. Astron. J. 82 (1977) 1013-1024.

[34] C. Mariotti, Lamb's problem with the lattice model Mka3D. Geophys. J. Int. 171 (2007) 857-864.

[35] J.J. Monaghan, Simulating free surface flows with SPH. J. Comput. Phys. 110 (1994) 399-406. | Zbl

[36] D.O. Potyondy and P.A. Cundall, A bonded-particle model for rock. Int. J. Rock Mech. Min. Sci. 41 (2004) 1329-1364.

[37] A. Ries, D.E. Wolf and T. Unger, Shear zones in granular media : Three-dimensional contact dynamics simulation. Phys. Rev. E 76 (2007) 051301.

[38] J.C. Simo, N. Tarnow and K.K. Wong, Exact energy-momentum conserving algorithms and symplectic schemes for nonlinear dynamics. Comput. Meth. Appl. Mech. Eng. 100 (1992) 63-116. | MR | Zbl

[39] Y. Suzuki and S. Koshizuka, A Hamiltonian particle method for non-linear elastodynamics. Int. J. Numer. Meth. Eng. 74 (2008) 1344-1373. | MR | Zbl

[40] J.W. Swegle, D.L. Hicks and S.W. Attaway, smoothed particle hydrodynamics stability analysis. J. Comput. Phys. 116 (1995) 123-134. | MR | Zbl

[41] K.Y. Sze, X.H. Liu and S.H. Lo, Popular benchmark problems for geometric nonlinear analysis of shells. Finite Elem. Anal. Des. 40 (2004) 1551-1569.

[42] H. Yserentant, A new class of particle methods. Numer. Math. 76 (1997) 87-109. | MR | Zbl