[论文解读] A Material Point Method Formulation for Plasticity

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

Fluid Dynamics Simulations and Interactions参考文献 10被引用 3

一句话总结

本文提出了一种用于率无关塑性的材料点法（MPM）公式，推导出一种弱形式，指出对变形粒子的体积积分是主要误差来源。该方法采用基于粒子的积分与变分框架，表明通过李导数和变形梯度的乘法分解改进应力更新算法，可提升MPM在塑性材料模拟中的精度。

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.

研究动机与目标

• 开发一种适用于加法分解率无关塑性的通用MPM公式。
• 识别出在MPM中对变形粒子进行体积积分是主要误差来源。
• 推导一种弱形式，以实现使用MPM框架稳定且准确地求解塑性问题。
• 为基于李导数与变形梯度的乘法分解改进应力更新算法奠定基础。

提出的方法

• 采用变分方法，推导塑性背景下质量、动量与能量守恒定律的弱形式。
• 通过使用基于粒子的权函数与有限元形函数，将材料点法应用于弱形式的离散化。
• 采用对变形粒子的一致积分方案，假设每个粒子具有恒定材料属性，以简化弱形式。
• 引入热惯性和刚度矩阵公式，以模拟系统中的热传导与能量平衡。
• 采用变分框架推导半离散化运动方程，其中包含温度及其节点值的时间导数。
• 提出改进的应力更新算法——可能通过李导数与变形梯度的乘法分解——可减少积分误差。

实验结果

研究问题

• RQ1 如何通过弱形式系统性地构建适用于率无关塑性的材料点法？
• RQ2 为何对变形粒子的体积积分会在MPM模拟中引入显著误差？
• RQ3 变形梯度的加法或乘法分解在提升应力更新精度中起什么作用？
• RQ4 弱形式如何被离散化以在塑性模拟中保持一致性和稳定性？
• RQ5 MPM在塑性中的稳定半离散化方程组的关键组成部分是什么？

主要发现

• 对变形粒子的体积积分被识别为MPM在塑性中主要的误差来源。
• 弱形式推导表明，必须谨慎处理基于粒子的积分，以保持动量与能量平衡的准确性。
• 热惯性和刚度矩阵由粒子积分项导出，从而实现一致的能量与热传导建模。
• 半离散化方程组以节点温度及其时间导数表示，其贡献来自粒子源与边界通量。
• 该公式表明，改进的应力更新算法——特别是基于李导数与变形梯度乘法分解的算法——可显著减少积分误差。
• 本文为未来在大变形与复杂材料行为背景下提升MPM塑性算法的研究奠定了基础。