4. 材料模型与应变率效应
碰撞仿真与准静态分析最大的差别，在于材料在毫秒级加载下表现出的应变率相关性。当应变率从 10⁻³ s⁻¹ 提升到 10²～10³ s⁻¹ 时，低碳钢屈服强度可提升 30%～80%，铝合金提升 10%～40%，而碳纤维复合材料则表现为模量基本不变、强度略升、失效应变下降。如果材料卡片仍沿用准静态参数，仿真会低估吸能量、高估侵入量，轻量化方案就会被迫"加厚保安全"，从而丧失减重空间。本节围绕 OptiStruct 中金属与复合材料碰撞材料卡片的设置方法，以及应变率参数对轻量化结果的影响评估展开。
一、金属材料：弹塑性 + 应变率 + 损伤三件套
OptiStruct 中金属碰撞常用 MATS1（弹塑性）、MATX42（应变率相关塑性）或 MAT1 + TABLES1 组合。推荐做法是：用 MAT1 定义弹性模量 E、泊松比 ν、密度 ρ，用 MATS1 定义硬化曲线，用 TABLES1 定义屈服应力随塑性应变率变化的缩放因子，再用 MATADD 或 MATT1 叠加失效准则。
以某 590DP 高强钢为例，准静态屈服 420 MPa，抗拉 780 MPa，密度 7.85e-9 t/mm³，弹性模量 210000 MPa。硬化曲线可用 Swift 模型拟合：σ = K(ε₀ + εₚ)ⁿ，取 K = 1200 MPa、ε₀ = 0.003、n = 0.16。应变率缩放因子按 Cowper-Symonds 形式给出：
σ_dyn / σ_static = 1 + (ε̇ / D)^(1/p)
590DP 常用 D = 8000 s⁻¹, p = 5.5
ε̇ = 1 s⁻¹ → 因子 ≈ 1.20
ε̇ = 100 s⁻¹ → 因子 ≈ 1.42
ε̇ = 1000 s⁻¹→ 因子 ≈ 1.55在 OptiStruct 中，这些因子以 TABLES1 的 (ε̇, 因子) 对形式输入，通过 MATX42 的 RATE 字段引用。损伤方面，金属常用等效塑性应变失效，设置 MATS1 的 FAIL 或 MATX42 的 EPSF，590DP 可取 0.35～0.55；若需模拟裂纹扩展，可叠加 MATDAMG 或 GISSMO 类损伤卡片。
二、复合材料：正交各向异性 + 多失效准则 + 应变率修正
碳纤维复合材料在 OptiStruct 中通常用 MAT8（二维正交各向异性）或 MAT9（三维）定义弹性常数，用 MAT8 的 XT/XC/YT/YC/S 定义拉伸、压缩、剪切强度，再通过 MATADD 叠加 MATUSR 或 MATDAMG 实现 Hashin、Puck 或 LaRC 失效判据。应变率效应在复合材料中相对温和，但冲击下基体剪切强度可提升 20%～50%，失效应变下降 10%～30%。
工程上常用做法是：对 MAT8 的 G12、S、YT 按应变率分档给出，例如 1 s⁻¹、100 s⁻¹、1000 s⁻¹ 三档，用 TABLES1 插值。若软件版本支持 MAT8 的率相关字段，可直接引用；否则用多工况材料卡片切换，在显式分析中按单元应变率自动选择。需要注意：复合材料层合板还要配合 PCOMP 定义铺层角度与厚度，碰撞中常见的分层可用 COHESIVE 或 MATDAMG 界面单元模拟。
三、应变率参数对轻量化结果的影响评估方法
要评估应变率效应到底带来多少减重空间，建议做一组对照仿真：同一白车身或电池包下箱体模型，分别用"准静态材料卡片"和"含应变率材料卡片"跑同一碰撞工况（如 50 km/h 正面 100% 重叠刚性墙），对比三个指标——最大侵入量、峰值加速度、总吸能量。若含率效应模型在侵入量满足目标的前提下吸能更高，说明原设计有减薄余量；反之则说明原设计偏薄，需要补强。
具体操作上，可在 OptiStruct 中建立两个 MAT 集合，用 INCLUDE 分别引用，通过 LOADSTEP 切换。评估时重点关注：
- 吸能量分配：前纵梁、门槛梁、B 柱等主要吸能件的内能占比是否合理，一般前纵梁应承担 40%～60%。
- 应变率分布：用
H3D后处理查看单元应变率云图，确认高应变率区域是否落在材料卡片覆盖范围内（如 100～1000 s⁻¹）。 - 失效单元比例：删除单元过多说明失效应变设置偏小，会低估吸能；删除过少则可能高估。
- 质量-性能灵敏度：以板厚为变量做 DOE，观察含率效应与不含率效应两条灵敏度曲线的斜率差异，斜率差越大，说明应变率对轻量化的贡献越显著。
四、一个可复用的材料卡片检查清单
| 检查项 | 金属（MATS1/MATX42） | 复合材料（MAT8/PCOMP） |
|---|---|---|
| 弹性常数 | E、ν、ρ 与准静态一致 | E1、E2、G12、ν12 完整 |
| 塑性/强度 | 硬化曲线覆盖 0～0.5 塑性应变 | XT/XC/YT/YC/S 四项齐全 |
| 应变率 | TABLES1 至少 3 档（1/100/1000 s⁻¹） | G12、S、YT 分档修正 |
| 失效 | EPSF 按实测标定，0.3～0.6 | Hashin/Puck 判据 + 分层界面 |
| 单位制 | t-mm-s-N-MPa 统一 | 同左，注意密度量级 |
| 验证 | 单胞拉伸/压缩仿真对标试验 | 层合板三点弯对标 |
完成材料卡片设置后，建议先用一个简单帽形梁或圆管做轴向压溃仿真，与试验的力-位移曲线对比，确认峰值力、平均压溃力和吸能量误差在 10% 以内，再投入整车或电池包级别的轻量化优化。只有材料卡片本身可信，后续的厚度优化、拓扑优化和形貌优化结果才具有工程意义。