Material point method enhanced by modified gradient of shape function

  • Authors:

  • Duan Z. Zhang;Xia Ma;Paul T. Giguere

  • Affiliations:

  • Theoretical Division, Fluid Dynamics and Solid Mechanics Group (T-3, B216), Los Alamos National Laboratory, Los Alamos, NM 87545, USA;Theoretical Division, Fluid Dynamics and Solid Mechanics Group (T-3, B216), Los Alamos National Laboratory, Los Alamos, NM 87545, USA;Theoretical Division, Fluid Dynamics and Solid Mechanics Group (T-3, B216), Los Alamos National Laboratory, Los Alamos, NM 87545, USA

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2011

Quantified Score

Hi-index 31.45

Visualization

Abstract

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.