岩土数值模拟：从Mohr-Coulomb到复杂本构的选型逻辑

做岩土数值模拟的人，大概率在第一次选本构模型的时候被劝退过。Mohr-Coulomb、Drucker-Prager、Hardening-Soil、Soft-Soil……光FLAC3D里内置的就有十几种，每个都有一堆参数要填。参数没搞对，算出来的位移能差一个数量级。

这不是软件的问题，是岩土本身的复杂性——土和岩石不是均匀连续介质，它们的力学行为强烈依赖应力历史、排水条件和加载速率。

本构模型选型的根本判据

选型的核心问题是：你的工程问题里，哪些力学行为是不可忽略的？

Mohr-Coulomb：最简单，只有粘聚力和内摩擦角两个强度参数。适合初步估算、应力水平不高、变形精度要求不严格的场景。缺点是忽略了土的压缩硬化和剪胀性，算出来的位移通常在真实值的0.5-2倍之间。

Drucker-Prager：Mohr-Coulomb的光滑化版本，适合需要用塑性势理论的有限元程序。在PLAXIS里用得比较多，但对土的剪胀描述仍然粗糙。

Hardening-Soil（HS）：岩土数值模拟里用得最广的进阶模型。考虑了剪切硬化和压缩硬化，用三个模量参数（E50^ref, Eur^ref, Eoed^ref）描述土的刚度随应力水平的变化。PLAXIS和FLAC3D都有实现，对砂土和硬粘土的变形预测相当靠谱。

Soft-Soil（SS）：专门为正常固结粘土和高灵敏度粘土设计的。考虑了体积压缩（virgin compression）和次压缩（secondary compression），对软土地基的长期沉降预测比HS模型更准确。

FLAC3D vs PLAXIS：工具选型

FLAC3D（Itasca）：基于有限差分法，显式求解。最大优势是大变形和接触面问题的处理——边坡失稳的大位移、隧道开挖的渐进坍塌，FLAC3D的描述比有限元更自然。缺点是前处理弱，建模效率低，结果后处理需要自己写FISH脚本。

PLAXIS（Bentley）：基于有限元法，隐式求解。界面友好，建模效率高，输出结果的可视化效果好。对常规边坡、地基、挡土墙问题，PLAXIS的工作流更高效。缺点是对大变形和接触面分离的描述不如FLAC3D。

实际选型建议：如果是科研或者复杂接触问题（比如桩-土相互作用、锚杆加固的渐进破坏），用FLAC3D；如果是工程设计（边坡稳定分析、地基承载力），用PLAXIS效率更高。

网格密度的敏感度和收敛判断

岩土数值模拟里，网格密度对结果的影响比结构力学大得多。土的本构关系有强烈的非线性，单元尺寸不同，积分点的应力路径不一样，得到的荷载-位移曲线能差20-30%。

收敛判断的标准做法：做网格收敛性验证。用粗网格、中网格、细网格各算一遍，看关心的结果量（比如边坡安全系数、地基最大沉降）是否随网格加密趋于稳定。通常单元尺寸缩小一半，结果变化小于5%，可以认为网格收敛。

FLAC3D里有个实用技巧：在关心区域（比如边坡潜在滑动面、隧道周边）做局部网格加密，远离关心区域用粗网格。这样既能保证精度，又控制计算规模。

边界条件的设置陷阱

底面边界：通常设为固定边界（零位移）。但如果模型深度不够，固定底面会约束掉真实的沉降变形。经验法则：模型底面距离工程尺度的3-5倍以下，固定边界的影响可以忽略。

侧面边界：通常设水平向固定、竖向自由。这个设置在均质场里没问题，但如果有地形坡度或者非均质土层，侧面边界会引入不合理的约束。更合理的做法是设rollers（滚轴边界），允许侧面土体在竖向自由变形。

地下水头边界：这是最容易设置错误的部分。孔隙水压力的分布直接影响有效应力，进而影响土的强度和变形。PLAXIS里有专用的”地下水”模块，可以定义复杂的水头边界和排水条件。

接触面和界面单元

土和结构物的相互作用，大多通过接触面来描述。FLAC3D里的interface单元和PLAXIS里的interface单元，本质都是薄层单元，用剪切刚度和法向刚度两个参数描述接触行为。

这两个参数很难通过实验确定。实际做法是用反分析法——拿一个已知的工程案例（比如现场载荷试验），调整接触面刚度参数，让模拟结果和现场测量匹配，再用这套参数去做新的设计计算。

常见坑点

坑一：参数取值没有来源。岩土本构模型的参数，理论上都可以通过三轴试验、固结试验确定。但实际工程里，往往是拿经验值或者文献值填进去。这样算出来的结果，可信度要打折扣。建议至少做敏感性分析，把关键参数（粘聚力c、内摩擦角φ、压缩模量E）在合理范围内扫一遍，看结果的变化幅度。

坑二：忽略排水条件。不排水条件下（比如快速填筑的路堤），孔隙水压力会升高，有效应力降低，土的强度打折。如果本构模型里没正确设置排水条件，算出来的安全系数会偏于危险。

坑三：把数值结果当精确值。岩土数值模拟的输出，本质上是”在给定参数和模型假设下的力学响应”。参数的不确定性、本构模型的不完备性，都会让数值结果和真实情况有偏差。负责任的做法是给出结果的范围（比如安全系数的上下限），而不是一个精确到小数点后两位的数字。

岩土数值模拟的价值，不在于算出一个”精确”的位移或安全系数，而在于帮你看清楚力学机制——哪里是潜在的滑动面？哪些参数对结果影响最大？加固措施改变了哪一部分的应力路径？这些定性但深刻的洞察，往往比一个数字更有用。