《Computational Particle Mechanics》：Assessment of the hydromechanical higher-order MPM for the simulation of geotechnical problems
编辑推荐：
本研究聚焦于利用高阶B样条材料点法（BS-MPM）框架，评估其在模拟岩土工程中流体与固体耦合大变形问题（如边坡失稳、地基贯入）的能力。研究通过验证该框架在固结、边坡破坏和地基贯入等真实尺度问题上的稳定性和准确性，证明BS-MPM能有效缓解数值计算中常见的“单元穿越”误差和体积自锁问题，为复杂岩土工程问题提供了更可靠的模拟工具。
岩土工程大变形模拟的数值挑战
在岩土工程领域，例如预测山体滑坡或桩基贯入过程，我们需要精确模拟土壤这类多孔颗粒介质在固体和流体（孔隙水）耦合作用下发生的大变形。这是一个至关重要的课题，因为它直接影响工程安全性和设计的可靠性。传统的数值方法，如有限元法，在处理材料发生巨大形变时常会遇到困难。近几十年来，材料点法作为一种基于粒子的数值方法脱颖而出，它能够有效模拟大变形。然而，常规使用线性插值函数的MPM自身也存在着固有的数值误差，这阻碍了其在真实世界复杂岩土问题中的应用。
其中，最突出的问题之一是“单元穿越”误差。当材料点（MP）在背景网格的单元间移动时，线性插值函数导数的突变会导致节点内力的振荡，产生类似冲击波的伪影，最终扭曲应力解并导致位移解的漂移。此外，当采用流体力学（Hydromechanical）公式模拟多孔介质（涉及两相、各自不同的速度场）时，这些误差会与“体积自锁”误差叠加。体积自锁是所有使用数值网格进行空间积分的技术都会遇到的典型问题，当模拟近乎不可压缩材料（如泊松比接近0.5的材料、在塑性流动下等体积变形的材料或初始不排水的多孔介质）的行为时，如果不进行稳定化处理，会表现为“棋盘式”应力解和受限的变形轮廓。
本研究的目标与框架
为了应对这些挑战，本文旨在评估、验证并深入讨论一种高阶B样条材料点法框架。该框架的核心是利用高阶（如三次、四次）B样条函数代替传统的线性函数进行插值，从根本上消除“单元穿越”误差，并结合一个简单的应变平滑算法来缓解体积自锁问题。这种方法在结构化的背景网格下，能够以比更复杂的等几何MPM更低的计算成本，获得平滑且稳定的数值解。
主要技术方法概览
研究采用了显式更新-应力后置的两相单点MPM公式来描述饱和多孔介质的固-液耦合行为。空间积分则通过B样条函数完成，采用Cox-de Boor递推公式生成高阶插值函数，确保了插值在整个定义域内的正定性和导数光滑性。为了对抗体积自锁，研究采用了一种应变平滑算法，该算法将材料点的体积应变映射到节点，再将平滑后的节点值映射回材料点，从而在不显著改变计算流程的前提下稳定应变场。此外，研究还利用了BS-MPM的多片能力，通过连接位于相同坐标的节点自由度，来模拟不同材料（如地基与土体）间的相互作用。
研究结果与分析
- 1.
验证大应变固结仪固结问题
首先，研究通过大应变固结仪固结问题验证了BS-MPM框架。通过与解析解对比发现：在较小荷载（1 kPa）下，线性MPM结果与解析解高度吻合。但随着荷载增大（10 kPa），材料点发生单元穿越，线性MPM的孔隙水压力解出现失真，沉降预测被高估。而采用更高阶的插值函数（二次、三次、四次）能显著改善这一状况。当使用三次（p=3）或四次（p=4）B样条函数时，即使在高达200 kPa的荷载下，数值解依然稳定且与解析解高度一致。
- 2.
验证真实尺度边坡破坏实验
其次，研究模拟了由孔隙水压力补给触发的Selborne真实尺度边坡破坏实验。通过对比现场测量数据，评估了高阶BS-MPM捕捉滑坡特征（如孔隙水压力、变形）的能力。结果表明，在模拟这种涉及大变形和渐进破坏的复杂岩土问题时，高阶BS-MPM框架能够稳定地捕捉滑坡发生和发展的过程，其预测结果与实测数据吻合良好，凸显了该方法处理真实世界问题的潜力。
- 3.
分析地基贯入问题的多片能力
最后，研究在多片BS-MPM框架下模拟了条形基础贯入饱和均质土体的过程。通过变化贯入速率和水力传导系数，考察了固结效应对基础反应力的影响。研究验证了BS-MPM在模拟不同材料界面相互作用（如地基与土体）时的稳定性，表明尽管插值函数导数在片间交界处不连续，但框架仍能提供稳定的数值解。
结论与意义
本研究表明，高阶B样条材料点法为模拟岩土工程中的大应变流体力学耦合问题提供了一种强大而稳定的数值工具。通过采用三次及以上的插值函数，该框架能够有效消除“单元穿越”误差，并结合应变平滑算法缓解体积自锁，从而在真实尺度的固结、边坡稳定性和地基贯入问题中，获得收敛且合理的数值解。相比于线性MPM和其他变体，高阶BS-MPM在保证计算精度的同时，提供了更好的数值稳定性，特别适用于涉及复杂本构关系和显著材料变形的工程应用场景。这项工作为MPM在岩土工程中的进一步应用和推广奠定了坚实的基础。
