简单初识 Material Point Method（MPM）

目录

1. Preliminaries

介绍Eulerian和Lagrangian

两者方法区别

2. MPM

Reference

---

1. Preliminaries

介绍Eulerian和Lagrangian

物理模拟有两种方式——Eulerian grids & Lagrangian particles.

(https://www.youtube.com/watch?v=zUaD-GMARrA 这个youtube视频也很不错！超级安利！！！从Fluid Flow开始引入，简单介绍了Lagrangian vs. Eulerian (In Simple Terms))

下面的定义摘自上面的video😊

1. Lagrangian：track an object to determine its properties.
2. Eulerian: It's just observing fluid properties as a function of time and space.

形象理解如下图手绘所示意：

两者方法区别

• Eulerian：只track网格格点的状态-->会被初始化的网格所限定，所以一些奇怪的，不规则的形状就很难被规整的网格所做到（要用大量的计算，把网格变得很小很小，才能去贴近物体的边界去计算）

• Lagrangian：初始化是初始化不同数量的粒子-->可以很容易track物体边界

1. The Lagrangian view follows an individual fluid parcel as it moves through space and time. The Eulerian view focuses on specific locations in the space through which fluid flows.
2. In Lagrangian specifications, the coordinates move with the fluid parcels so it is useful for following trajectories, deformation, and rotations. Eulerian coordinates are fixed in space so better for analyzing flows, velocities, accelerations.

---

 

2. Material Point Method (MPM)

MPM是一种混合的(Hybrid)，结合Eulerian和Lagrangian的方法。

定义：MPM is a hybrid Eulerian-Lagrangian approach, which uses moving material points and a fixed computational grid. This approach is particularly useful in the context of large deformations where mesh distortions often produce unrealistic or undesirable results.

• Lagrangian：认为粒子带有材质属性。
• Eulerian：使用网格来做数值积分。
• lnteraction：粒子将属性传递给网格，网格执行更新，然后插值回粒子。

 示意图注释：

1）左上图（由很多点组成了一个field，然后放在网格上，粒子将属性传递给网络）

2）右上图（motion equations是在nodes点上计算，然后更新点的kinematics）

3）右下图（之后把网格点的kinematics，插值回粒子）

4）左下图（之后更新回粒子本身）

---

 

Reference

1. 超级无敌安利的参考资料：https://www.math.ucla.edu/~cffjiang/research/mpmcourse/mpmcourse.pdf 

关于上面的pdf有个知乎也不错：模拟连续介质的物质点法 - 知乎

2.  一些详细的公式和笔记也可以参考这篇详细的文章：图形学笔记（二十）粒子、刚体、流体的模拟—— 欧拉方法、Errors 和 Instability、中点法、自适应步长、隐式欧拉方法、Runge-Kutta方法、刚体与流体模拟（质点法、网格法、MPM）_隐式mpm-CSDN博客