A High Order Material Point Method

A High Order Material Point Method

Tielen, R.P.W.M. (TU Delft Numerical Analysis)
Wobbes, Elizaveta (TU Delft Numerical Analysis)
Möller, M. (TU Delft Numerical Analysis)
Beuth, Lars (Deltares)

2017

The classical material point method (MPM) developed in the 90s is known for drawbacks which affect the quality of results. The movement of material points from one element to another leads to non-physical oscillations known as ‘grid crossing errors’. Furthermore, the use of material points as integration points renders a numerical quadrature rule of limited quality. Different solutions have been proposed in recent years to overcome these drawbacks. In this paper the approach of combining quadratic B-spline basis functions with a reconstruction based quadrature rule is pursued to solve these numerical problems. High-order B-spline basis functions solve the problem of grid crossing completely, whereas the considered reconstruction based quadrature rule reduces the quadrature error observed with MPM. In addition, the use of quadratic B-splines leads to a more accurate piecewise linear approximation of the stress field compared to the piecewise constant one obtained with linear Lagrangian basis functions commonly used with MPM. Two 1D benchmarks are considered involving large deformations, a vibrating bar and a column under self-weight. They render excellent results when adopting this high-order MPM.

material point method
B-splines
cubic splines

http://resolver.tudelft.nl/uuid:40e1b7ee-07b2-41c1-9389-fd12cbf7190e

https://doi.org/10.1016/j.proeng.2017.01.022

1877-7058

Procedia Engineering, 175, 265-272

1st International Conference on the Material Point Method, 2017-01-10 → 2017-01-13, Deltares, Delft, Netherlands

Institutional Repository

journal article

© 2017 R.P.W.M. Tielen, Elizaveta Wobbes, M. Möller, Lars Beuth