岩土工程博0一枚，也是数值算法开发的小白，被taichi简洁优雅的代码以及精美的效果吸引。于是我趁着毕业的空档研究了taichi的一些关于模拟的项目，比如SPH和MPM，发现图形学在本构计算中常用的一些变量和我们有较大的差别，比如我们在工程中会经常遇到到柯西应力张量，偏应力，偏应变等这些变量。
于是我冷饭新炒，对张雄老师一篇关于MPM的文章：Zhang X, Chen Z, Liu Y. The material point method: a continuum-based particle method for extreme loading cases[M]. Academic Press, 2016 进行复现，实现了非关联流动D-P砂土模型的MLS-MPM。
同时也参考了蒋老师文章中从粒子应力计算网格节点力的策略：Jiang C, Schroeder C, Teran J, et al. The material point method for simulating continuum materials[M]//Acm siggraph 2016 courses. 2016: 1-52. 另外边界加上了摩擦（感谢 g1n0st）。
二维的效果看起来还行，可是三维的就有点勉强（如果代码没错的话）。
代码托管在：GitHub - AUTUMN9452/MLS_MPM_sand_DPmodel: Application of the DP constitutive model to the MLS-MPM
总结：欢迎各个领域的大佬提提意见和交流 
，土木老铁代码水平有限，后面有空再写写砂-水的耦合。
