第三章 塑性材料模型基础：屈服准则、流动法则与硬化法则

各位同学，大家好。今天我们聊一个硬核话题——塑性材料模型。说实话，搞有限元分析这么多年，我见过太多人栽在材料模型上。模型选对了，结果八九不离十；模型选错了，算出来的东西就是一堆漂亮的数据垃圾。

塑性，说白了就是材料「回不去了」。你拉一根钢筋，弹性阶段松手它还能缩回去，但一旦过了屈服点，它就永久变形了。嗯，我们今天要讲的，就是怎么用数学语言描述这个「回不去」的过程。

3.1 屈服准则——材料什么时候「扛不住」？

屈服准则，说白了就是判断材料什么时候开始永久变形。你想想看，一个零件受力，到底多大的力才算「过界」？这就是屈服准则要回答的问题。

我个人习惯把屈服准则分成两大类：一类适用于金属，另一类适用于岩土和混凝土。金属材料通常用Von Mises或Tresca，而岩土材料则常用Drucker-Prager。

3.1.1 Von Mises屈服准则

Von Mises准则，也叫畸变能准则。它认为当材料的畸变能达到某个临界值时，材料开始屈服。这个准则在工程中用得最多，ANSYS里默认的金属塑性模型基本都基于它。

数学表达式长这样：

σ_vm = √[ ( (σ₁-σ₂)² + (σ₂-σ₃)² + (σ₃-σ₁)² ) / 2 ] ≤ σ_y

其中σ_vm是Von Mises等效应力，σ_y是屈服应力。

3.1.2 Tresca屈服准则

Tresca准则，也叫最大剪应力准则。它认为当最大剪应力达到临界值时，材料屈服。公式更直观：

τ_max = (σ₁ - σ₃) / 2 ≤ k

其中k是纯剪切屈服应力。

我记得有一次做压力容器分析，客户指定要用Tresca。说实话，Tresca比Von Mises保守一些，算出来的安全裕度更大。但ANSYS里Tresca的收敛性不如Von Mises好，我建议你慎用。

3.1.3 Drucker-Prager屈服准则

Drucker-Prager准则，是Von Mises的「升级版」。它考虑了静水压力对屈服的影响，所以特别适合岩土、混凝土这类材料。

公式长这样：

f = √J₂ + α·I₁ - k = 0

其中J₂是偏应力第二不变量，I₁是应力第一不变量，α和k是材料参数。

你想想看，沙子埋得越深，抗剪强度越大——这就是静水压力的影响。Von Mises描述不了这个，但Drucker-Prager可以。

准则名称	适用材料	考虑静水压力	ANSYS中关键字
Von Mises	金属（钢、铝、铜等）	否	MISO, BKIN
Tresca	金属（偏保守）	否	TB,TRESCA
Drucker-Prager	岩土、混凝土、聚合物	是	TB,DP

3.2 流动法则——塑性应变往哪个方向走？

屈服准则告诉我们「什么时候屈服」，但没告诉我们「屈服后应变怎么发展」。流动法则就是干这个的。

流动法则的核心思想是：塑性应变增量的方向，垂直于屈服面。嗯，这个叫「正交法则」。

3.2.1 关联流动 vs 非关联流动

关联流动：塑性势函数等于屈服函数。说白了，应变方向完全由屈服面决定。金属材料通常用这个。

非关联流动：塑性势函数不等于屈服函数。岩土材料常用这个，因为它的剪胀行为需要单独描述。

流动法则的数学形式：

dεᵖ = dλ · ∂g/∂σ

其中dεᵖ是塑性应变增量，dλ是塑性乘子，g是塑性势函数，σ是应力张量。

我曾经遇到一个案例：某土坝分析，用了关联流动，结果算出来的体积膨胀大得离谱。后来换成非关联流动，结果就合理了。嗯，岩土问题一定要小心这个。

3.3 硬化法则——屈服面会「长大」吗？

材料屈服之后，继续加载会发生什么？有的材料越压越硬（比如冷拔钢筋），有的材料屈服应力会漂移（比如循环加载）。硬化法则就是描述这个过程的。

3.3.1 各向同性硬化

各向同性硬化：屈服面均匀膨胀。说白了，材料在各个方向上都变强了。这个模型简单，适合单调加载。

在ANSYS中，用TB,MISO（多线性各向同性硬化）或TB,BISO（双线性各向同性硬化）来定义。

! ANSYS命令流示例：定义双线性各向同性硬化
TB,BISO,1,1,2
TBTEMP,20
TBDATA,1,250e6,2000e6  ! 屈服应力250MPa，切线模量2000MPa

3.3.2 随动硬化

随动硬化：屈服面大小不变，但中心位置移动。这个模型适合循环加载，能描述包辛格效应。

你想想看，一根钢筋先拉后压，拉伸屈服了，再压缩时屈服应力会降低——这就是包辛格效应。各向同性硬化描述不了这个，但随动硬化可以。

ANSYS中用TB,KINH（多线性随动硬化）或TB,BKIN（双线性随动硬化）。

! ANSYS命令流示例：定义双线性随动硬化
TB,BKIN,1,1,2
TBTEMP,20
TBDATA,1,250e6,2000e6  ! 屈服应力250MPa，切线模量2000MPa

3.3.3 混合硬化

混合硬化：各向同性硬化和随动硬化的组合。屈服面既膨胀又移动。这个模型最接近真实材料行为，但参数也多，标定起来比较麻烦。

ANSYS中可以用TB,CHABOCHE（Chaboche模型）来实现混合硬化。这个模型有多个背应力分量，能很好地描述循环塑性行为。

硬化类型	屈服面变化	适用场景	ANSYS关键字
各向同性硬化	均匀膨胀	单调加载、大变形	BISO, MISO
随动硬化	中心平移	循环加载、疲劳	BKIN, KINH
混合硬化	膨胀+平移	复杂加载路径	CHABOCHE

3.4 实战建议——怎么选？

讲了这么多理论，你可能会问：到底怎么选？我根据自己的经验，给你几个「傻瓜式」建议：

1. 金属结构静力分析：Von Mises + 各向同性硬化。简单、稳定、够用。
2. 金属疲劳/循环加载：Von Mises + 随动硬化。别忘了包辛格效应。
3. 岩土/混凝土：Drucker-Prager + 非关联流动。静水压力效应必须考虑。
4. 聚合物/泡沫材料：Drucker-Prager或自定义模型。这类材料体积变形大，要小心。
5. 复杂加载路径：混合硬化。虽然参数多，但精度最高。

好了，这一章的内容就到这里。塑性材料模型是有限元分析的「地基」，地基没打好，上面盖的楼再漂亮也是危楼。希望你能把屈服准则、流动法则、硬化法则这三个概念吃透，后面做复杂分析时才能游刃有余。

---

公众号:蓝海资料掘金营，微信deep3321