Material

Point Method
Theory and Applica ons

Phu Nguyen

Department of Civil Engineering

Monash University

[email protected]

http://nvinhphu.wixsite.com/mysite

ti

Collaborators and funding agencies

Alban de Vaucorbeil
Deakin University

Sina Sinaie
Melbourne University

Tushar Mandal
PhD candidate

The problem: deforma on/fracture of solids

ti
Deforma on/fracture of solids: approaches

ti
 What computa onal engineer/scien st do

D⇢
Dt
Balance equa
C ⇢r v D 0ons
Cons tu ve models
Dv
⇢ Dr C ⇢b
Dt
De
⇢ DDW
Dt
r
<latexit sha1_base64="oQbK15w0nL9DZXgnd1W9dHV4C3M=">AAAD/XicdVLLbtQwFE0aHiW8WliysZiCihiNkhZBQUKqykiwLIJpK43TynacmWj8iGxn6CiK+AO28APsEFu+hTU/gpPJ0KZTrhTl6N5zjq+vL85Yqk0Q/HZXvCtXr11fveHfvHX7zt219XsHWuaK0AGRTKojjDRlqaADkxpGjzJFEceMHuLJm6p+OKVKp1J8NLOMRhyNRJqkBBmbOll3XRgnCpGiD6Aay7LomxI8rTGAAmGGACSxNABOeTEtwePXIAAQDncyE/kVayGf1+f6itUWYz3jBdTpiKMze6vApb8wA0120Q4952WZ/fJV2wWcCdv+x83RVvjhxBw3adAvN62PtkbdFv3JyVon6AV1gGUQNqDjNLFvh/YHxpLknApDGNJ6GAa2kQIpkxJG7ZVyTTNEJmhEhxYKxKmOivqxSvDIZmKQSGU/YUCdPa8oENe2PWyZHJmxvlirkpfVhrlJdqIiFVluqCDzg5KcASNB9fIgThUlhs0sQESltldAxsgO29j98CGjBu5hPIGKMnTqt6x1jrWVZIZIVZnY1elqhoSh8VtF6aQrZHWZVCSyy002VnFU2L+SW6VvnQT9RCTnSMQFNFNeDsPIAnpqcFJ0wrJcYuiasgGxZHE1C8kq3gYEF6gLr2oYl3otrP7PqFZhTmkfVtZr8bKK5/+WYBkcbPXC7d72+2ed3b1mQVadB85DZ9MJnRfOrvPO2XcGDnHH7hf3q/vN++x99354P+fUFbfR3Hda4f36C4hMUis=</latexit>
D S t D .D; /
ti
ti
ti

ti
ti
 Mesh based and mesh-free methods

Basic concept: a mesh/grid Basic concept: points/par cles

FDM, FVM, FEM EFG, SPH, PFEM, MPM, …
E cient, accurate Mesh distor on free
FEM: prone to mesh distor on Inaccurate geometry
ffi
ti

ti
ti

Material Point Method (MPM): par cles over a xed grid

grid nodes

posi on
volume
par cles momenta
stress
temperature
…

Sulsky et al, CMAME, 1994

ti
ti

ti
fi

Material Point Method: algorithm (explicit, ULMPM)

1 P2G 2 Grid update

4 Grid reset 3 G2P

MPM really can simulate complex problems

MPM: Example 1

Image based analysis

Foamed materials: network of struts
Very large deforma on
Lots of self contacts

Sulsky et al, CMAME, 1994

Bardenhangen et al, JMPS, 2005

ti

MPM: example 2

Hyper-velocity impact
Steel plate: Johnson-Cook model
Massive deforma on
Karamelo (in-house C++ MPI)
ti

 no mesh genera on
simple basis func ons
large deforma on
automa c no-slip contact
mul ple contacts
simple BC treatment
simple implementa on
ti
ti
ti
ti
ti

ti

 MPM really can simulate complex problems

BUT

(UL)MPM: issues

cell-crossing instability
(UL)MPM: issues

numerical fracture
ULMPM: issues and par al solu ons

ti
ti

We present two solu ons to solve issues of ULMPM

ti

Total Lagrangian MPM (TLMPM): ini al con gura on

use the ini al undeformed con gura on

no cell-crossing instability
no numerical fracture
most accurate quadrature
no inherent no-slip contact
small me step for compression

A. de Vaucorbeil, V. P. Nguyen, and C. R. Hutchinson. A total-lagrangian

material point method for solid mechanics problems involving large
deforma ons. CMAME, 2020.
A. de Vaucorbeil and V.P. Nguyen. Modeling contacts with a total lagrangian
mater al point method. CMAME 2021.
ti
ti
ti
ti



fi
ti

ti
fi
ti

 Generalised Par cle in Cell (GPIC): mesh again

solids are meshed
ini al con gura on: forces
current con gura on: contact

no cell crossing issue

no numerical fracture
accurate quadrature
accurate boundary
easy boundary condi ons
easy material interfaces

fragmenta on?
extreme deforma on?

1. A. Sadeghirad, R. M. Brannon, and J. Burghardt. A convected par cle domain interpola on technique to
extend applicability of the material point method for problems involving massive deforma ons. Interna onal
Journal for Numerical Methods in Engineering, 86(12):1435–1456, 2011.
2. V. P. Nguyen, C. T. Nguyen, T. Rabczuk, and S. Natarajan. On a family of convected par cle domain
interpola ons in the material point method. Finite Elements in Analysis and Design, 126:50–64, 2017.
3. V.P. Nguyen, A. de Vaucorbeil, C. Nguyen-Thanh, and T. Mandal. A generalized par cle in cell method for
explicit solid dynamics. Computer Methods in Applied Mechanics and Engineering, 360:112783, 2020.
ti
ti
fi
ti
fi

ti
ti
ti

ti

ti

ti

ti
ti
ti
ti

ti
Generalised Par cle in Cell (GPIC): example 1

ti

Generalised Par cle in Cell (GPIC): example 1

ABAQUS GPIC

ti

Generalised Par cle in Cell (GPIC): example 2

ti

Scratch test of copper: TLMPM, Johnson-Cook plas c model

de Vaucorbeil et al, CMAME, 2019

de Vaucorbeil et al, AAM, 2019

ti

Scratch test of copper: groove topology

 Modelling fracture

Damage mechanics approach: Fracture mechanics approach:

+ easy to code (2D, 3D) + capture discon nui es

- less e cient (re ned meshes) + e cient
- hard to capture discon nui es - hard to implement
- required non locality - tracking 3D non-planar crack surfaces
ffi
ffi

ti
fi

ti

ti

ti

TLMPM: modeling duc le fracture

GRID
Damage mechanics approach
Local models
Nonlocal models

mesh dependent
mesh biased

Local model

ti

 Nonlocal gradient enhanced JC damage: length scale

R. H. J. Peerlings, R. de Borst, W. A. M. Brekelmans, and J. de Wree. Gradient enhanced damage for quasi bri le materials.
Interna onal Journal for Numerical Methods in Engineering, 39:3391–3403, 1996
ti
tt

TLMPM—modelling duc le fracture: result 1

ti

TLMPM—modelling duc le fracture: result 2

ti

TLMPM—modelling duc le fracture: result 2

ti

TLMPM—modelling duc le fracture: result 2

ti

Summary

no mesh genera on
simple basis func ons
large deforma on
mul ple contacts
simple BC treatment
simple implementa on

approximate geometry
blur material interfaces
numerical fracture: ULMPM
ti

ti
ti
ti

ti

 The end
Phu Nguyen

Department of Civil Engineering

Monash University

[email protected]

http://nvinhphu.wixsite.com/mysite

Tri-Uni Workshop, November 2020