We use cookies to ensure that we give you the best experience on our website. By continuing to browse this repository, you give consent for essential cookies to be used. You can read more about our Privacy and Cookie Policy.

Durham Research Online
You are in:

iGIMP : an implicit generalised interpolation material point method for large deformations.

Charlton, T.J. and Coombs, W.M. and Augarde, C.E. (2017) 'iGIMP : an implicit generalised interpolation material point method for large deformations.', Computers and structures., 190 . pp. 108-125.


The Material Point Method (MPM) uses a combined Eulerian-Lagrangian approach to solve problems involving large deformations. A problem domain is discretised as material points which are advected on a background grid. Problems are encountered with the original MPM when material points cross between grid cells, and this has been tackled by the development of the Generalised Interpolation MPM, where material points’ domains of influence extend beyond the currently occupied grid cell. In this paper, the Generalised Interpolation Material Point (GIMP) Method has been implemented implicitly in a manner that allows a global stiffness matrix to be constructed similar to that in the Finite Element Method (FEM) by combining contributions from individual elements on the background grid. An updated Lagrangian finite deformation framework has been used to ensure non-linear behaviour within each of the loadsteps. The weighting functions used for this which make the GIMP method different to standard MPM are presented and the implementation is explained. Specific details on computing the deformation gradient to be consistent with the updated Lagrangian framework and the updating of the material point influence domains are outlined, both of which are currently unclear in the published literature. It is then shown through numerical examples that for both small and large deformation elastic and elasto-plastic problems, the implicit GIMP method agrees well with analytical solutions and exhibits convergence properties between that of linear and quadratic FEM.

Item Type:Article
Full text:(AM) Accepted Manuscript
Available under License - Creative Commons Attribution.
Download PDF
Full text:(VoR) Version of Record
Available under License - Creative Commons Attribution.
Download PDF
Publisher Web site:
Publisher statement:© 2017 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license (http://
Date accepted:13 May 2017
Date deposited:15 May 2017
Date of first online publication:31 May 2017
Date first made open access:No date available

Save or Share this output

Look up in GoogleScholar