ICCM Conferences, The 8th International Conference on Computational Methods (ICCM2017)

Font Size: 
A smoothed MPM algorithm for elastodynamics
Lisha He, Zhen Chen, Yong Gan

Last modified: 2017-06-28


To accurately capture the dynamics in elastic materials by the computational methods, the artificial viscosity is necessary to damp out the numerical oscillation. The choice of the artificial viscosity is always problem-dependent, which is different equations for different types of problems [1]. The material point method (MPM) has been developed to model the multi-phase (solid-fluid-gas) interactions including large deformation and failure evolution [2, 3]. The SPH method has been also widely applied to simulate the fluid and solid mechanics [1]. However, the performance of both MPM and SPH methods in elastodynamics will be influenced by an unsuitable artificial viscosity. This work will focus on improving the capability of MPM in elastodynamics without using artificial viscosity, by combining the advantages of the SPH method and the MPM to develop a new smoothed MPM algorithm. The proposed method is more accurate and stable than both MPM and SPH for elastodynamics. It is different from the coupling of the MPM and SPH method in work [4]. The smoothed MPM is based on introducing the SPH reconstruction scheme into the MPM mapping and remapping procedure. This algorithm will include the advantages of both two methods, such as the high efficiency and convenient boundary conditions in the MPM method, the capability of modelling the free surface in the SPH method. Three elastodyamical examples will be presented to demonstrate the capability of the proposed algorithm. The first example is a one-dimensional wave propagation problem, which was simulated by classical MPM in previous work [5]. The compression problem on a two-dimensional elastic plane with circular hole is modeled by the SMPM and then compared with the FEM and SPH simulation in work [6] in the second example. In the third example, the impact problem between two-dimensional elastic bars is simulated by the SMPM. There is no artificial viscosity in all these three examples. Future work will be discussed based on the current work.


[1]     Liu, G.R., Liu, M.B. (2004) Smoothed Particles Hydrodynamics: A Meshfree Particle Method. World Scientific Publishing, Singapore.

[2]     Sulsky, D., Chen, Z. and Schreyer H.L. (1994) A particle method for history-dependent materials. Computer Methods in Applied Mechanics and Engineering, 118:179–196, (1994).

[3]     Zhang, X., Chen, Z. and Liu Y. (2016) Material Point Method – A Continuum-Based Particle Method for Extreme Loading Cases, Academic Press, Elsevier.

[4]     Raymond, S.J., Jones, B. and Williams, J.R. (2016) A strategy to couple the material point method (MPM) and smoothed particle hydrodynamics (SPH) computational techniques. Computational Particle Mechanics, 10.1007/s40571-016-0149-9.

[5]     Chen, Z. and Brannon, R. (2002) An evaluation of the material point method. Sandia Report SAND2002-0482, Sandia National Laboratories, Albuquerque, NM, USA.

[6]     He, L. and Seaid, M. (2016) A Runge–Kutta–Chebyshev SPH algorithm for elastodynamics. Acta Mechanica, 227:1813–1835.


MPM, SPH, large deformation, elastodynamics, Numerical simulations

An account with this site is required in order to view papers. Click here to create an account.