2、土力学基础回顾：有效应力原理、土的抗剪强度、土的本构模型概述
各位同学，大家好。欢迎来到《单桩基础桩土相互作用模拟教程》的第二讲。
在正式开始用数值软件“折腾”桩土相互作用之前，我觉得有必要先带大家把土力学里最核心的几个概念捋一遍。为什么？因为数值模拟说白了就是“用数学公式去模仿物理现象”。如果物理概念本身是模糊的，那算出来的结果再漂亮，也只是“垃圾进，垃圾出”。
我个人习惯是，不管做多复杂的模拟，第一步永远是先问自己三个问题：土里的力是怎么传递的？土什么时候会坏？土坏的时候是什么样子的？ 这三个问题，正好对应我们今天要讲的三个核心内容：有效应力原理、抗剪强度、以及本构模型。
2.1 有效应力原理——土力学的“牛顿定律”
有效应力原理，是太沙基老爷子在1923年提出来的。你想想看，如果没有这个原理，我们根本没法解释为什么饱和砂土在震动下会液化，也没法解释为什么打桩时桩端阻力会随着深度增加。
原理本身很简单，就一个公式：
σ' = σ - u其中：
- σ：总应力（土体自重+外荷载）
- u：孔隙水压力（静水压力或超静水压力）
- σ'：有效应力（真正作用在土骨架上的力）
说白了，土体不是一块实心的铁疙瘩。它是由土颗粒、孔隙水和孔隙气组成的“三相体”。你施加一个荷载，这个荷载一部分被水“扛”了（孔隙水压力），另一部分才真正压在土颗粒骨架上（有效应力）。而只有有效应力，才会引起土体的变形和强度变化。
在数值模拟中，有效应力原理的应用体现在哪里？
- 初始地应力平衡： 必须先计算自重作用下的有效应力场和孔压场。
- 排水/不排水分析： 排水分析时，孔压会消散，有效应力增加；不排水分析时，孔压会累积，有效应力不变。
- 液化判断： 当孔压升高到等于总应力时，有效应力为零，土体失去抗剪强度——这就是液化。
2.2 土的抗剪强度——土什么时候会“坏”？
土体的破坏，本质上都是剪切破坏。不管是边坡滑坡、地基失稳，还是桩基的刺入破坏，背后都是剪应力超过了土的抗剪强度。
抗剪强度的经典表达式是摩尔-库仑准则：
τ_f = c' + σ' · tan(φ')其中：
- τ_f：抗剪强度
- c'：有效粘聚力（颗粒间的胶结力）
- σ'：破坏面上的有效正应力
- φ'：有效内摩擦角（颗粒间的咬合摩擦）
这个公式告诉我们两件事：
- 土越压越强： 有效正应力越大，抗剪强度越高。这就是为什么深基础比浅基础承载力高——桩端处的土被压得更紧实。
- 粘聚力是“底牌”： 对于砂土（c'≈0），抗剪强度完全来自摩擦；对于粘土，粘聚力提供了初始强度。
在数值模拟中，抗剪强度参数（c', φ'）是定义材料破坏准则的核心输入。无论是摩尔-库仑模型还是Drucker-Prager模型，最终都要用到这两个参数。
2.3 土的本构模型概述——从“线性”到“非线性”
本构模型，说白了就是描述土体“应力-应变”关系的数学公式。你给它一个应力增量，它告诉你土体会产生多少变形。
土体的应力-应变关系非常复杂：非线性、弹塑性、剪胀性、硬化/软化…… 没有一个模型能完美描述所有情况。我们只能根据工程问题的特点，选择合适的模型。
2.3.1 摩尔-库仑模型（M-C模型）
这是岩土工程中最经典、最常用的模型。它的核心思想是：
- 弹性阶段： 应力低于抗剪强度时，土体表现为线弹性（用弹性模量E和泊松比ν描述）。
- 塑性阶段： 应力达到抗剪强度时，土体发生屈服（破坏），屈服面由摩尔-库仑准则定义。
摩尔-库仑模型的屈服面在主应力空间中是一个六棱锥。它的优点是参数物理意义明确（c', φ'），工程经验丰富。但缺点也很明显：
- 在三维应力空间中，棱锥的角点会导致数值计算不稳定（应力奇异）。
- 不能很好地模拟土体的剪胀性和硬化/软化行为。
2.3.2 Drucker-Prager模型（D-P模型）
Drucker-Prager模型可以看作是摩尔-库仑模型的“光滑版本”。它的屈服面在主应力空间中是一个圆锥，没有角点，数值稳定性更好。
D-P模型的屈服准则表达式为：
f = √J₂ + α·I₁ - k = 0其中：
- J₂：偏应力第二不变量（反映剪切作用）
- I₁：应力第一不变量（反映静水压力作用）
- α, k：材料参数，可以通过与摩尔-库仑参数（c', φ'）换算得到
D-P模型在有限元软件（如Abaqus、ANSYS）中非常常见。它的优点是：
- 计算效率高，收敛性好。
- 可以考虑静水压力对屈服的影响（土体的压硬性）。
但缺点也很明显：
- 参数α和k的物理意义不如c'和φ'直观。
- 对于纯摩擦型材料（c'=0），D-P模型的预测结果可能与摩尔-库仑模型有偏差。
2.4 三种模型对比总结
| 模型名称 | 屈服面形状 | 核心参数 | 优点 | 缺点 | 适用场景 |
|---|---|---|---|---|---|
| 线弹性模型 | 无屈服面 | E, ν | 简单、计算快 | 不能模拟破坏 | 小变形、低荷载 |
| 摩尔-库仑 | 六棱锥 | c', φ', ψ | 参数直观、经验丰富 | 角点奇异、无硬化 | 极限承载力、边坡稳定 |
| Drucker-Prager | 圆锥 | α, k (由c', φ'换算) | 数值稳定、收敛好 | 参数不直观 | 通用有限元分析 |
好了，这一讲的内容就到这里。我们回顾了有效应力原理、抗剪强度公式，以及两个最常用的本构模型。这些概念是后续进行桩土相互作用模拟的“地基”。地基不牢，楼盖得再高也是危楼。
下一讲，我们会正式进入数值模拟的世界，聊聊如何用Python和有限元软件搭建一个简单的单桩模型。到时候，今天讲的这些概念都会派上用场。