搞懂这5种屈服准则,有限元分析少走一半弯路
做有限元分析(FEA)时,屈服准则 的选择直接决定了塑性变形计算的准确性。简单来说,屈服准则就是一个“开关”——它告诉求解器:材料在什么复杂应力状态下开始进入塑性流动。
在主应力空间里,这个“开关”对应的数学表达式,就是我们常说的屈服方程。
今天我们谈谈工程中最常用的5种屈服准则:Tresca、Von-Mises、Mohr-Coulomb、Drucker-Prager 和 Zienkiewicz-Pande。前两种专攻金属,后三种则更“接地气”,专门对付混凝土、岩石和土体这类材料。
1. Tresca 准则——最大剪应力说了算
这是最古老的屈服准则之一,由 Tresca 在1864年提出。它的逻辑非常直白:只要材料内部某个方向上的最大剪应力达到极限值,材料就屈服了。
如果三个主应力满足 σ₁ ≥ σ₂ ≥ σ₃,那条件就简化为:
σ₁ – σ₃ = σₛ(单轴拉伸屈服应力)
如果主应力大小顺序不确定,则需要用这个通用形式:
max{|σ₁–σ₂|, |σ₂–σ₃|, |σ₃–σ₁|} = 2k(k 为纯剪屈服应力)
在 π 平面上,Tresca 准则画出来是一个正六边形;在主应力空间里,则是正六棱柱面。
✅ 优点:主应力排序已知时,手算很方便,公式简单。
❌ 缺点:完全无视中间主应力 σ₂ 的影响,也跟静水压力没关系;而且屈服面存在尖角,数值计算时容易“卡住”。
🎯 适用:金属材料的初步估算,尤其是主应力方向明确的问题。
2. von Mises 准则——形状改变能攒够了就屈服
这是工程中最常用的金属屈服准则,由 von Mises 在1913年提出。它的物理角度很巧妙:当单位体积内用于改变形状的那部分弹性能达到材料自身极限时,屈服就发生了。
数学表达式大家都很熟悉:
√{½[(σ₁–σ₂)² + (σ₂–σ₃)² + (σ₃–σ₁)²]} = σₛ
也可以写成偏应力第二不变量 J₂ 的形式:
√(3J₂) = σₛ
在主应力空间里,von Mises 屈服面是一个光滑的圆柱体,在 π 平面上投影是个标准圆形。
✅ 优点:考虑了中间主应力 σ₂ 的影响;屈服面光滑,数值计算非常顺畅;材料参数少,通过单轴拉伸实验就能确定。
❌ 缺点:仍然不涉及静水压力;对岩土类材料完全“不感冒”——既反映不了压力敏感性,也体现不了拉压强度不同的 S-D 效应。
🎯 适用:钢材、铝材等各向同性金属材料的塑性分析,CAE 中的绝对主力。
3. Mohr-Coulomb 准则——给摩擦材料量身定制的
前面两个准则放到岩石、土和混凝土上,偏差会很大。为什么?因为这些材料吃“压力”——静水压力越大,强度越高。
Mohr 换了个思路:他认为,材料某个平面上发生屈服的剪应力 τₙ 并不是常数,而是与这个平面上的正应力 σₙ 有关,线性关系如下:
τₙ = c – σₙ·tanφ
其中 c 为粘聚力,φ 为内摩擦角。当主应力排序为 σ₁ ≥ σ₂ ≥ σ₃ 时,可以改写成用主应力表达的形式:
(σ₁–σ₃)/2 = (σ₁+σ₃)/2·sinφ + c·cosφ
✅ 优点:有效反映了静水压力对屈服强度的提升;能表征岩土材料拉压强度不同的特性;参数 c 和 φ 物理意义明确,实验测定方便。
❌ 缺点:仍然忽略中间主应力 σ₂;屈服面在 π 平面是一个不等角六边形,存在尖点,有限元迭代时容易引起数值奇异。
🎯 适用:岩石、土体、混凝土等摩擦性材料的稳定性和承载力分析。
4. Drucker-Prager 准则——M-C 的光滑“改良版”
Drucker-Prager 准则(1952年)可以理解为对 von Mises 的扩展——在 J₂ 的基础上,加入了一项与静水压力相关的 I₁(应力第一不变量),从而让金属准则“学会”感受压力。
其屈服函数为:
f = √J₂ + α·I₁ – k = 0
其中 α 和 k 是由 c、φ 换算得来的材料常数(具体换算方式取决于与 M-C 外角点还是内角点匹配)。
在主应力空间中,D-P 屈服面是一个圆锥面,π 平面上是圆形。
✅ 优点:同时考虑了中间主应力和静水压力的影响;屈服面光滑,数值鲁棒性好;参数简单,可从 M-C 参数直接换算。
❌ 缺点:采用线性近似,无法反映岩土材料在高静水压力下的非线性屈服特性;在偏平面上也体现不了拉压强度的 S-D 差异。
🎯 适用:岩土工程和混凝土结构的简化弹塑性分析,计算效率较高。
5. Zienkiewicz-Pande 准则——M-C 的“无死角”升级
M-C 准则最大的痛点就是尖点奇异性,尤其在有限元迭代中,角点处无法确定唯一的塑性流动方向。为此,Zienkiewicz 和 Pande 在1977年提出了一个光滑化改进版本。
Z-P 准则在子午面(p–q 平面)上用双曲线、抛物线或椭圆型光滑函数替代了 M-C 的直线,从而在 π 平面上实现了圆角过渡,彻底消除了尖点。
其统一形式可写为:
f = F(p, q) – q_c = 0
其中 F(p, q) 是关于平均应力 p 和广义剪应力 q 的光滑凸函数。
✅ 优点:子午面和 π 平面均光滑,数值计算非常友好;一定程度上反映了静水压力与屈服强度的非线性关系;也能部分计入中间主应力效应。
❌ 缺点:参数标定相对复杂,工程普及度不如前几种高。
🎯 适用:对精度要求较高的岩石力学和混凝土非线性有限元分析。
一张表帮你理清怎么选
为了方便大家在项目中快速决策,我们把五个准则的核心特征和适用场景汇总成两张表。
表1 各准则核心数学特征对比
| 准则名称 | π平面形状 | 是否考虑静水压力 | 是否考虑中主应力σ₂ | 屈服面是否光滑 |
|---|---|---|---|---|
| Tresca | 正六边形 | 否 | 否 | 有尖角 |
| von Mises | 圆形 | 否 | 是 | 光滑 |
| Mohr-Coulomb | 不等角六边形 | 是 | 否 | 有尖角 |
| Drucker-Prager | 圆形 | 是 | 是 | 光滑 |
| Zienkiewicz-Pande | 圆角化形状 | 是(可非线性) | 部分考虑 | 光滑 |

表2 工程应用建议速查
| 准则名称 | 主要优势 | 主要短板 | 推荐应用场景 |
|---|---|---|---|
| Tresca | 主应力排序明确时计算简单 | 忽略σ₂和静水压力,面有尖角 | 金属材料初步估算 |
| von Mises | 公式简洁,数值友好,应用最广 | 不适用于压力敏感材料 | 钢材/铝材等塑性分析 |
| Mohr-Coulomb | 参数物理意义明确,反映压力敏感 | 忽略σ₂,尖点导致数值困难 | 岩土/混凝土稳定性分析 |
| Drucker-Prager | 兼顾压力与σ₂,计算高效 | 线性近似,无S-D效应 | 岩土工程简化弹塑性分析 |
| Zienkiewicz-Pande | 全光滑曲面,精度较高 | 参数标定略复杂 | 高精度岩石/混凝土分析 |

评论
等天黑

干货推荐

lammps中文教程12
春风得意马蹄疾
86 2

二氧化碳捕集及资源化利用:从"废气"到"碳资源"的产业革命
春风得意马蹄疾
76 2

lammps中文教程5
春风得意马蹄疾
47 2

【代码分享】费米面嵌套(Fermi surface nesting)函数计算
张阳小号1
650 2

小胶质细胞
Mr弘🔬
62 8

分子动力学模拟:从原始轨迹到科学结论的推理路径解析
吃鱼不
123 3
