DYNA MAT24材料本构模型详解
2026年9月13日 19:51
一、 MAT_024 的灵魂：四大核心机制
MAT_024 的本质可以总结为一句话：基于 von Mises 屈服准则，结合关联流动法则，支持任意形状硬化曲线、任意应变率硬化效应及塑性失效删除的通用弹塑性模型。
弹性段 (E, ν) ──> von Mises 屈服 (J2) ──> 关联流动 ──> 硬化/率效应 (LCSS/C-P) ──> 塑性失效 (FAIL/TDEL)
1. 屈服准则：经典 J2（von Mises）屈服面
MAT_024 采用典型的 $J_2$ 异性各向同性屈服准则，其屈服函数 $\Phi$ 表达为：
$$\Phi(\boldsymbol{\sigma}) = \sigma_{\text{eq}} - \sigma_y(\bar{\varepsilon}^p) = 0$$
其中：
- $\sigma_{\text{eq}} = \sqrt{\frac{3}{2}\boldsymbol{S}:\boldsymbol{S}}$ 为 von Mises 等效应力（$\boldsymbol{S}$ 为偏应力张量）；
- $\sigma_y(\bar{\varepsilon}^p)$ 为当前等效塑性应变 $\bar{\varepsilon}^p$ 下的屈服应力。
- $\Phi < 0$：处于弹性阶段，满足 Hooke 定律；
- $\Phi = 0$：材料发生塑性屈服；
- $\Phi > 0$：通过径向返回算法（Radial Return Algorithm）将应力拉回至屈服面上。
2. 流动法则：关联流动（Associative Flow Rule）
塑性应变增量沿屈服面的外法线方向流动：
$$\Delta\boldsymbol{\varepsilon}^p = \Delta\gamma \frac{\partial \Phi}{\partial \boldsymbol{\sigma}} = \frac{3}{2} \Delta\gamma \frac{\boldsymbol{S}}{\sigma_{\text{eq}}}$$
这保证了塑性流动纯偏量化且体积不可压缩（$\text{tr}(\Delta\boldsymbol{\varepsilon}^p) = 0$），符合绝大多数金属材料在塑性变形时的物理表现。
3. 硬化机制：为什么它能击败 MAT_003？
传统 MAT_003（双线性塑性）的硬化段只能用一条直线（$E_{\text{tan}}$）拟合，这对于存在屈服平台、非线性抛物线硬化的真实金属（如低碳钢、高强钢、铝合金）来说误差极大。
MAT_024 允许用户直接输入实测的真实应力-真塑性应变曲线（通过 LCSS 定义）：
应力 σy ▲ / 真实曲线 (MAT_024) │ /---\___ │ / ----\___ │ / ---- 直线近似 (MAT_003) │ / │ ____/ (屈服平台) └──────────────────────────────► 塑性应变 εp
每个塑性应变区间的斜率都能精准匹配，吸能量与压溃载荷大幅提升。
4. 应变率效应（Rate Dependency）
动态加载下（如碰撞、冲击），金属通常会“越快越硬”。MAT_024 提供了三种精度的应变率处理方式：
- Cowper-Symonds 模型（简单快手）：
- $$\sigma_y' = \sigma_y \cdot \left[1 + \left(\frac{\dot{\varepsilon}}{C}\right)^{1/P}\right]$$
- （$C, P$ 为拟合参数，低碳钢典型值如 $C=40\text{ s}^{-1}, P=5$）。
LCSR缩放曲线：定义“屈服应力放大倍数 vs 应变率”曲线。LCSS表格模式（多率曲线族）：挂载包含不同应变率（如 0.001/s, 1/s, 100/s, 1000/s）下的实测硬化曲线族，求解时自动按当前塑性应变率进行线性或对数插值。这是精度最高的形式。
二、 参数卡片结构与输入逻辑
MAT_024 的输入卡片逻辑非常清晰，工程运用中最常用的是 Card 1 和 Card 2：
*MAT_PIECEWISE_LINEAR_PLASTICITY
$ Card 1
$# MID RO E PR SIGY ETAN FAIL TDEL
1 7.85E-09 210000.0 0.3 0.0 0.0 0.0 0.0
$ Card 2
$# C P LCSS LCSR VP
40.0 5.0 100 0 1
关键字段解析及优先级
字段 |
含义 |
调参建议/注意事项 |
LCSS |
硬化曲线/表格 ID |
核心输入。横坐标必须是等效塑性应变，绝对不能错填为总应变！ |
|
屈服强度 / 切线模量 |
当定义了 |
VP |
应变率定义标志 |
|
FAIL |
塑性失效应变删除 |
|
TDEL |
最小时间步删除 |
仅作为压溃严重、网格畸变导致时间步塌缩时的兜底手段，勿当作正常失效使用。 |
输入优先级逻辑：
LCSS表格/曲线 >EPS1~EPS8 / ES1~ES8内联点 >SIGY + ETAN线性硬化
三、 典型应用场景与对比选型
虽然 MAT_024 被称为“万能接盘侠”，但在复杂的工程场景中，盲目套用也会撞墙。
场景 |
推荐材料 |
为什么不选/慎选 MAT_024？ |
整车碰撞 / 低速跌落 |
MAT_024 |
完美契合（各向同性、率敏感、大塑性吸能）。 |
薄板成形与回弹 |
MAT_033 / 036 |
MAT_024 是各向同性屈服，无法反映轧制各向异性（$r$ 值），会导致回弹预测偏小。 |
高速穿甲 / 绝热温升 |
|
MAT_024 不具备热软化耦合机制，高速热塑性剪切带无法表现。 |
复杂应力状态断裂 |
|
MAT_024 仅提供常数塑性失效应变，无法表达三轴度对损伤演化的影响。 |
四、 工程实战避坑指南（80% 新手踩过的坑）
坑 1：直接把材料手册的“断后伸长率”填入 FAIL
- 错因：断后伸长率是宏观标距测量值（包含均匀延伸与颈缩）；而
FAIL是单元级的局部失效应变，与网格尺寸强相关。 - 对策：直接填伸长率会导致结构提前撕裂。建议先设
FAIL=0跑基线，若需模拟断裂，需通过单胞拉伸试验结合网格敏感性分析进行标定。
坑 2：LCSS 输入了总应变（Total Strain）
- 错因：若把总应变曲线直接挂在
LCSS上，求解器会误将弹性段也当成塑性段，导致材料在极小载荷下就开始发生塑性流动，整体响应大幅偏软。 - 对策：输入前必须剥离弹性应变：$\varepsilon^p = \varepsilon^{\text{total}} - \frac{\sigma}{E}$。
坑 3：碰撞分析忘记开率效应，或者 VP 选错
- 错因：在 100/s 的应变率下，低碳钢的屈服强度可能会提升一倍以上。不开率效应会导致压溃力与吸能量严重低估。
- 对策：确保输入 $C/P$ 或
LCSR；强烈建议将VP设置为 1（忽略弹性应变率波动的干扰）。
五、 标准标定三板斧 workflow
想让你的 MAT_024 模型打进整车级/工程级仿真，推荐遵循以下四步走流程：
┌─────────────────────────┐
│ 1. 准静态单轴拉伸试验 │ ───► 换算真应力-真塑性应变 ───► 建立准静态 LCSS 曲线
└─────────────────────────┘
│
▼
┌─────────────────────────┐
│ 2. 高速拉伸 / SHPB 试验 │ ───► 标定 C/P 参数或构建多率表格
└─────────────────────────┘
│
▼
┌─────────────────────────┐
│ 3. 单单元 (Single Element)│ ───► 准静态/动态拉伸数值测试 ───► 确保输出与输入曲线完全重合
└─────────────────────────┘
│
▼
┌─────────────────────────┐
│ 4. 组件级验证 (落锤/三点弯)│ ───► FAIL=0 跑基线 ───► 引入标定后的失效应变与网格敏感性验证
└─────────────────────────┘
总结
MAT_024 的成功秘诀就在于“简单、能打、鲁棒性强”。
弄懂它的 J2 屈服骨架、关联流动与率效应插值机制，避开失效应变与应变输入的常见误区，你就能发挥出这个“国民本构”的最大威力，为你的仿真模型提供坚实可靠的物理保障！
技术专长
全部评论 (0)