干货有限元 干货详情

搞懂这5种屈服准则,有限元分析少走一半弯路

等天黑99
有限元分析有限元仿真弹塑性分析数值模拟非线性有限元塑性力学屈服准则本构模型强度理论静水压力中间主应力金属材料岩土材料混凝土岩石力学土体岩土工程结构稳定性承载力分析

做有限元分析(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全光滑曲面,精度较高参数标定略复杂高精度岩石/混凝土分析


评论

0/1000发布评论
全部评论

等天黑

关注
TA的主页

干货推荐