von Mises塑性理论
von Mises塑性理论的理论基础
von Mises塑性理论是什么
老师,von Mises塑性理论是FEM材料非线性的基础吧。
von Mises屈服条件是描述金属塑性变形的最基本理论。「等价应力(von Mises应力)达到屈服应力时,塑性变形开始」。
von Mises等价应力
或分量表示:
屈服条件
$f < 0$ 时处于弹性范围。$f = 0$ 时在屈服面上(发生塑性变形)。$f > 0$ 不允许(不能超过屈服面)。
应力空间中的「球面」是屈服面吧。
在偏差应力空间中,von Mises屈服面是一个圆柱。不依赖于静水压力(体积应力)是von Mises的特点。金属的塑性变形不伴随体积变化(不可压缩塑性流动),因此在物理上是合理的。
硬化法则
屈服后的应力-应变关系(硬化法则):
| 硬化类型 | 屈服面变化 | 应用 |
|---|---|---|
| 完全弹塑性(完全塑性) | 屈服面固定 | 塌陷荷载评估 |
| 各向同性硬化 | 屈服面膨胀 | 单调加载 |
| 运动硬化(运动学) | 屈服面平移 | 循环加载(疲劳) |
| 混合硬化 | 膨胀+平移 | 最通用 |
Abaqus
```
*MATERIAL, NAME=steel
*ELASTIC
200000., 0.3
*PLASTIC
250., 0.0 $ 屈服应力250 MPa, 塑性应变0
400., 0.1 $ 400 MPa, 塑性应变10%
500., 0.3 $ 500 MPa, 塑性应变30%
```
Nastran
```
MAT1, 1, 200000., , 0.3
MATS1, 1, , PLASTIC, , , 1, 1
TABLES1, 1, , ,
, 0.0, 250., 0.1, 400., 0.3, 500., ENDT
```
用应力-塑性应变表格来定义硬化曲线啊。
将拉伸试验的工程应力-工程应变曲线转换为真应力-真应变后输入FEM。大变形分析中真应力-真应变是必需的。
总结
要点:
- $\sigma_{vm} = \sigma_Y$ 时屈服 — 金属塑性的基础
- 不依赖静水压力 — 无体积变化(金属特性)
- 硬化法则 — 完全塑性/各向同性硬化/运动硬化/混合硬化
- 真应力-真应变输入FEM — 需要从工程值进行转换
- 所有FEM求解器标配 — 最基本的材料非线性模型
von Mises的1913年论文
Richard von Mises在1913年于哥廷根科学协会杂志上提出了用第二偏差应力不变量J₂=k²表示的屈服准则。Hencky在1924年添加了剪切应变能达到临界值时屈服的物理解释。该准则在主应力空间中呈现为圆柱面,是目前最广泛使用的屈服基准。
von Mises塑性理论的数值计算方法
返回映射算法
老师,塑性的数值处理怎么做啊?
返回映射(Return Mapping)算法是标准方法:
1. 弹性预报(elastic predictor) — 假设应变增量全为弹性,计算临时应力
2. 屈服判定 — 临时应力是否在屈服面外?
3. 塑性修正(plastic corrector) — 若在屈服面外,将应力「返回」到屈服面上
「先用弹性计算→再返回屈服面」的两步法啊。
von Mises塑性的返回可以用径向返回(radial return)严格计算。计算效率高,数值稳定。所有商用求解器都实现了这个方法。
切线刚度矩阵(CTO)
塑性状态下的切线刚度(Consistent Tangent Operator, CTO):
$H$ 是硬化系数。$\{n\}$ 是屈服面的法向量。CTO保证牛顿-拉夫逊法的二阶收敛。
总结
径向返回映射
von Mises+各向同性硬化的FEM实现使用「径向返回映射」。包含弹性预报→超过屈服面的确认→沿切线方向返回三个步骤。由Simo & Taylor(1985年)证明了线性收敛性。通过使用切线弹性模量可以在一次迭代内获得严格解,计算成本低,因此几乎所有通用求解器都采用了这个方法。
von Mises塑性理论的工程应用
塑性分析的工程应用
von Mises塑性用在什么场景呢?
金属结构的全部非线性强度评估:
| 应用 | 目标 |
|---|---|
| 耐压试验的弹塑性分析 | 压力容器的ASME Div.2 Part 5 |
| 塑性塌陷荷载评估 | 极限荷载法(是否在2倍设计荷载处收敛) |
| 板料成形(冲压) | 变形后的形状和回弹 |
| 焊接残余应力 | 焊接→冷却的热弹塑性分析 |
| 地震弹塑性时程 | 塑性铰的形成 |
真应力-真应变转换
从拉伸试验数据(工程应力-工程应变)转换为FEM输入:
1. 均匀变形段以内 — $\sigma_{true} = \sigma_{eng}(1+\varepsilon_{eng})$, $\varepsilon_{true} = \ln(1+\varepsilon_{eng})$
2. 均匀变形段以后(颈缩后) — 上述变换不准确。需要逆分析或修正式
3. 塑性应变 — $\varepsilon_{pl} = \varepsilon_{true} - \sigma_{true}/E$
颈缩后的转换很困难啊。
颈缩开始后,应力状态不再是单轴(变成三轴应力)。简单的变换式不准确。需要Bridgman修正或逆FEM法(用FEM结果和试验的力-位移曲线匹配拟合)。
工程检查清单
汽车碰撞分析的主角
汽车前碰分析(NCAP符合)中,前侧成员通常采用SPCC(冷轧钢板)的von Mises+各向同性硬化模型。使用LS-DYNA的MAT_024+MATSUMOTO回弹模型的组合,碰撞行程200mm的最大荷载预测误差控制在±8%以内是2000年代以后的行业标准。
von Mises塑性理论的软件对比
von Mises塑性的工具
所有FEM求解器都标配支持。没有区别。
| 求解器 | 设置 |
|---|---|
| Abaqus | *PLASTIC 表 |
| Nastran | MATS1 + TABLES1 |
| Ansys | TB, BISO or TB, MISO |
| LS-DYNA | *MAT_24(弹塑性) |
LS-DYNA的MAT_24在碰撞安全中应用最广吧。
MAT_24是von Mises弹塑性 + 各向同性硬化 + 应变速率相关性(Cowper-Symonds)。汽车钢板碰撞分析中几乎是唯一的选择。
选择指南
所有主要求解器都支持
von Mises屈服准则在所有商用CAE求解器中都是标准实现。Abaqus(*PLASTIC)、LS-DYNA(MAT_024)、MSC Nastran(SOL 400 MATS1)、ANSYS Mechanical(Bilinear/Multilinear Isotropic)、Marc(Yield Criterion=VonMises)。这是从1950年代FEM黎明期开始就实现的历史最悠久的非线性材料模型。
von Mises塑性理论的前沿研究
晶体塑性(Crystal Plasticity)
von Mises是宏观各向同性塑性,但实际金属是多晶体。按各个晶粒的滑移系分别计算塑性的晶体塑性(Crystal Plasticity FEM, CPFEM)是研究热点。能预测织构(晶体取向)发展和各向异性演化。
延性破坏的耦合
von Mises塑性 + 延性破坏准则(Johnson-Cook、Gurson等)的耦合。大塑性应变时材料破坏。对碰撞和成形的破坏预测至关重要。
机器学习构成模型
用神经网络学习应力-应变关系的「数据驱动构成模型」。超越von Mises假设,能表现任意材料响应。
总结
GTN多孔塑性损伤模型
Gurson(1977年)的多孔塑性模型在von Mises准则中加入孔隙体积率f,表现延性破坏。Tvergaard和Needleman(1984年)引入了参数q₁~q₃并完成为GTN模型。现在在Abaqus UMAT中广泛用于厚板深拉伸分析和高强度钢(HSS)的孔扩张分析。
von Mises塑性理论故障排除
von Mises塑性的故障
塑性分析常见的故障有哪些?
应力超过屈服应力
von Mises应力超过屈服应力…硬化的影响。各向同性硬化中屈服面膨胀,所以$\sigma_{vm} > \sigma_{Y,initial}$ 是正常的。应该和对应塑性应变的硬化后屈服应力对比。
工程应力-工程应变直接输入
将工程值直接输入FEM的后果:
- 小应变(< 5%)时基本准确
- 大应变(> 10%)时应力被低估、应变被高估
对策:转换为真应力-真塑性应变后输入。
牛顿-拉夫逊法不收敛
塑性变形很大时:
- 减小荷载增量
- 启用自动时间步
- 确认NLGEOM=YES已设置(大塑性应变伴随大变形)
体积锁定
塑性变形是不可压缩的(ΔV = 0)。1阶单元的完全积分会产生体积锁定。
对策:
- C3D8R(低减积分)或C3D8RH(混合法)
- C3D10M(改进的TET10)
- B-bar法(LS-DYNA ELFORM=2)
总结
体积锁定的对应
完全不可压缩塑性(ν≈0.5)使用通常的8节点六面体单元时会产生体积锁定,导致位移被低估。对策是选择性低减积分(Abaqus C3D8R)或处理不可压缩条件的混合单元(C3D8H)。一般在无硬化的完全塑性分析中,使用低减积分单元加沙漏控制是最安全的选择。
相关主题
更详细
报告