[Paper Review] A Material Point Method Formulation for Plasticity

B. Banerjee|arXiv (Cornell University)|Jan 13, 2012

Fluid Dynamics Simulations and Interactions10 references3 citations

TL;DR

This paper presents a material point method (MPM) formulation for rate-independent plasticity, deriving a weak form that identifies volume integration over deforming particles as a key source of error. The approach uses a variational framework with particle-based integration and demonstrates that improved stress update algorithms—potentially via Lie derivatives and multiplicative decomposition—can enhance accuracy in MPM simulations of plastic materials.

ABSTRACT

This paper discusses a general formulation of the material point method in the context of additive decomposition rate-independent plasticity. The process of generating the weak form shows that volume integration over deforming particles can be major source of error in the MPM method. Several useful identities and other results are derived in the appendix.

Motivation & Objective

• To develop a general MPM formulation for additive decomposition rate-independent plasticity.
• To identify volume integration over deforming particles as a major source of error in MPM.
• To derive a weak form that enables stable and accurate solution of plasticity problems using the MPM framework.
• To lay the groundwork for improved stress update algorithms using Lie derivatives and multiplicative decomposition of the deformation gradient.

Proposed method

• Derives the weak form of the balance laws for mass, momentum, and energy in the context of plasticity using a variational approach.
• Applies the material point method by discretizing the weak form using particle-based weighting functions and finite element shape functions.
• Uses a consistent integration scheme over deforming particles, assuming constant material properties per particle to simplify the weak form.
• Introduces a thermal inertia and stiffness matrix formulation to model heat transfer and energy balance in the system.
• Employs a variational framework to derive semidiscrete equations of motion, with time derivatives of temperature and nodal values.
• Proposes that improved stress update algorithms—potentially using Lie derivatives and multiplicative decomposition of the deformation gradient—can reduce integration errors.

Experimental results

Research questions

• RQ1 How can the material point method be systematically formulated for rate-independent plasticity using a weak form?
• RQ2 Why does volume integration over deforming particles introduce significant errors in MPM simulations?
• RQ3 What is the role of the deformation gradient's additive or multiplicative decomposition in improving stress update accuracy?
• RQ4 How can the weak form be discretized to maintain consistency and stability in plasticity simulations?
• RQ5 What are the key components of a stable semidiscrete system of equations for MPM in plasticity?

Key findings

• Volume integration over deforming particles is identified as a primary source of error in the material point method for plasticity.
• The weak form derivation reveals that particle-based integration must be carefully handled to preserve accuracy in momentum and energy balance.
• The thermal inertia and stiffness matrices are derived from particle-integrated terms, enabling consistent energy and heat transfer modeling.
• The semidiscrete system of equations is expressed in terms of nodal temperature and its time derivative, with contributions from particle sources and boundary fluxes.
• The formulation suggests that improved stress update algorithms—particularly those based on Lie derivatives and multiplicative decomposition of the deformation gradient—can significantly reduce integration errors.
• The paper establishes a foundation for future work on enhanced MPM algorithms for plasticity, especially in the context of large deformations and complex material behavior.