Material point method enhanced by modified gradient of shape function

A numerical scheme of computing quantities involving gradients of shape functions is introduced for the material point method (MPM), so that the quantities are continuous as material points move across cell boundaries. The noise and instability caused by cell crossing of the material points are then eliminated. In this scheme, the formulas used to compute these quantities can be expressed in the same forms as in the original material point method, but with the gradient of the shape function modified. For one-dimensional cases, the gradient of the shape function used in the generalized interpolation material point (GIMP) method is a special case of the modified gradient if the characteristic function of a material point is introduced. The characteristic function of a material point is not otherwise needed in this scheme, therefore difficulties in tracking its evolution are avoided. Although the support of the modified gradient of a shape function is enlarged from the cell containing the material point to also include the immediate neighbor cells, all the non-local effects of a material point can be accounted for by two consecutive local operations. Therefore this scheme can be used in calculations with unstructured grids. This scheme is proved to satisfy mass and momentum conservations exactly. The error in energy conservation is shown to be second order on both spatial and temporal discretizations. Although the error in energy conservation is the same order as that in the original material point method, numerical examples show that this scheme has significantly better energy conservation properties than those of the original material point method.

[1]  Charles R. Johnson,et al.  Matrix analysis , 1985, Statistical Inference for Engineers and Data Scientists.

[2]  A. A. Amsden,et al.  FLUID DYNAMICS. A LASL Monograph. , 1971 .

[3]  Qisu Zou,et al.  Material point method applied to multiphase flows , 2008, J. Comput. Phys..

[4]  Deborah Sulsky,et al.  Implicit dynamics in the material-point method , 2004 .

[5]  James Guilkey,et al.  An evaluation of explicit time integration schemes for use with the generalized interpolation material point method , 2008, J. Comput. Phys..

[6]  Deborah Sulsky,et al.  Mass matrix formulation of the FLIP particle-in-cell method , 1992 .

[7]  H. Benaroya,et al.  DYNAMICS OF TRANSVERSELY VIBRATING BEAMS USING FOUR ENGINEERING THEORIES , 1999 .

[8]  M. Berzins,et al.  Analysis and reduction of quadrature errors in the material point method (MPM) , 2008 .

[9]  A. Love.  A treatise on the mathematical theory of elasticity , 1892 .

[10]  L. F..  A Treatise on the Mathematical Theory of Elasticity , Nature.

[11]  D. Sulsky.  Erratum: Application of a particle-in-cell method to solid mechanics , 1995 .

[12]  Qisu Zou,et al.  Pressure calculations in disperse and continuous multiphase flows , 2007 .

[13]  S. Bardenhagen,et al.  The Generalized Interpolation Material Point Method , 2004 .

[14]  Balaji Jayaraman,et al.  Distribution coefficient algorithm for small mass nodes in material point method , 2010, J. Comput. Phys..

[15]  A. Sadeghirad,et al.  A convected particle domain interpolation technique to extend applicability of the material point method for problems involving massive deformations , 2011 .