汽车材料的冲击力学行为表征和仿真

官方工学老师·15 页·深入(追求细节与边界)·0 次浏览·2 天前
冲击力学材料表征CAE仿真汽车工程

汽车材料的冲击力学行为表征和仿真

搞懂材料冲击响应的测试方法、关键表征参数与主流仿真建模校核策略

按 空格/→ 演示下一步

1 / 15 页

全部页面点击任意一页,跳回舞台从这页播放

冲击力学材料表征CAE仿真汽车工程

汽车材料的冲击力学行为表征和仿真

搞懂材料冲击响应的测试方法、关键表征参数与主流仿真建模校核策略

1第 1 页 · 汽车材料的冲击力学行为表征和仿真

碰撞响应的复杂性从何而来

碰撞看上去只是金属瞬间被压扁——但车身真实响应远不止这一眼。同样一块钢板在不同位置为何表现迥异?同样一根梁在不同速度下为何塌缩方式完全不同?根源在结构与材料两层非线性的互相放大。

结构: 多路径分叉
多腔体薄壁件屈曲、分叉、瞬态接触,载荷在多条路径间动态切换
材料: 多场耦合
应力同时取决于应变、应变率、温度、损伤度,四者实时反馈
结构-材料耦合
局部屈曲瞬时改变应变率分布,率硬化反过来决定屈曲是否触发
台风天的大树对应 →汽车碰撞响应

枝杈交错对应车身多路径传力,木材随风变脆对应材料率效应,两者互相影响

σ=f(ε,ε˙,T,D)\sigma = f(\varepsilon,\, \dot{\varepsilon},\, T,\, D)
2第 2 页 · 碰撞响应的复杂性从何而来

结构-材料响应的层级关系

纵向读尺度递降(米→亚纳米),每行标出层级、特征尺度、主流仿真方法。跨越一级,物理模型通常也要换。

图解渲染中…
L4本构 = 应力-应变关系,连接工程结构与材料L5晶粒、相、析出物,决定强度与失效模式L6晶格、位错、键合,变形与断裂的源头
3第 3 页 · 结构-材料响应的层级关系

应变率效应:材料的速度敏感性

上页我们从结构讲到材料层级——那材料本身会不会跟着速度变?答案是会,而且变化很剧烈。同样一块钢板,缓慢拉伸和被子弹击穿时的力学响应,几乎不是同一种材料。这一页解释为什么。

应变率
单位时间内的应变变化量,量纲 s⁻¹;汽车碰撞中材料点典型范围 10²~10³ s⁻¹
应变率硬化
流动应力随应变率升高而增大,正的率敏感性是金属普遍特征
高率脆化
变形太快来不及散热,绝热温升使韧性下降,断裂模式由韧转脆
准静态 vs 冲击曲线
同一材料两条 σ-ε 曲线明显分离——这是碰撞仿真不可用准静态数据的根本原因
慢慢拉太妃糖对应 →应变率对塑性的影响

拉的速度=应变率;慢拉延展成丝(低率韧),快拉直接脆断(高率脆)

σdσs=1+(ε˙C)1/p\frac{\sigma_d}{\sigma_s} = 1 + \left(\frac{\dot{\varepsilon}}{C}\right)^{1/p}
4第 4 页 · 应变率效应:材料的速度敏感性

温度敏感性:冷脆与热软化

上一节我们看到,材料对冲击速度很敏感:快慢不同,强弱可能差好几倍。但还有第二个隐藏的「开关」在暗中控制它的性格——温度。北极的钢板脆如薄冰,赤道的车身却软得像面团。

冷脆现象
温度下降时,BCC金属(铁素体钢)冲击功急剧下降,断裂从韧性转为脆性解理
韧脆转变温度DBTT
表征该转变的临界温度;车用钢常要求-40℃下冲击功仍达标
热软化
温度升高,流动应力下降;热冲压(700–900℃)正是利用此特性
BCC敏感 / FCC不敏感
体心立方钢温度效应显著;面心立方铝、奥氏体钢几乎无DBTT
冰箱里的巧克力棒对应 →冷脆与热软化

冷藏的巧克力一掰就脆断,室温的可弯成弧;更冷=更脆、更热=更软

σ=(A+Bεn)(1+Clnε˙)(1Tm)\sigma = (A+B\varepsilon^n)(1+C\ln\dot{\varepsilon}^*)(1-T^{*m})
5第 5 页 · 温度敏感性:冷脆与热软化

应力状态多样性

同一块钢板,拽它会越拉越细再断开;压它会鼓起来,承载力还更大。为什么同一种材料在不同受力方式下表现差这么多?

单轴拉伸
材料颈缩拉长,流动应力最低,真实断裂应变大
单轴压缩
出现鼓形,承载力通常比拉伸高 30%–100%
纯剪切路径
流动应力最低,剪切带易局部化、温升集中
静水压力效应
压应力让材料'变硬',是 Drucker-Prager 类本构的关键项
碰撞中的体现
前纵梁压溃吸能;B 柱弯曲时一侧受拉一侧受压
捏气球 vs 拽气球 vs 扭气球对应 →压、拉、剪三种应力状态

捏气球变扁难爆(压加固);拽气球越拉越薄最后爆(颈缩);扭气球表面起皱(局部化)

σy=σ0+ασH\sigma_y = \sigma_0 + \alpha \cdot \sigma_H
6第 6 页 · 应力状态多样性

多场耦合环境

前几页我们分别看了应变率、温度、应力状态各自对材料的影响——但在真实碰撞里,这三者从不孤立出现。漠河冬夜的追尾、时速过百的正面冲击、大塑性变形时材料自身还会发热,三种效应同时发生。

非线性叠加
三场共同决定材料响应,效应不是简单相加而是相互放大
温-率耦合
低温放大应变率敏感效应;高温下率效应反而弱化
变形生热耦合
大塑性变形约 90% 转热,绝热温升可让局部瞬时升 200°C
需耦合本构
必须建立描述三场交互的方程,孤立拟合任一曲线都会失真
寒冬雨夜高速驾驶对应 →三场耦合环境

低温(冷脆)、高速(高应变率)、复杂路况(多应力)同时作用,驾驶难度非线性放大,非各条件简单相加

σ=(A+Bεn)(1+Clnε˙)(1Tm)\sigma = \left(A + B\varepsilon^n\right)\left(1 + C\ln\dot{\varepsilon}^*\right)\left(1 - T^{*m}\right)
7第 7 页 · 多场耦合环境

分离式霍普金森杆试验原理

横看是物理过程:撞击→入射→反射/透射→采集→反演;右回流是波的回弹。

图解渲染中…
a2入射杆贴应变片,记录入射波εᵢ和反射波εᵣa4透射杆贴应变片,记录透射波εₜa7三波法反演:εᵣ算应变率与应变,εₜ算应力
8第 8 页 · 分离式霍普金森杆试验原理

冲击试验方法对比

霍普金森杆与落锤都叫"冲击",但应变率、试样、目的完全不同——选错方法,结论就废了。

分离式霍普金森杆
  • 应变率 10²–10⁴ s⁻¹,中高应变率区段
  • 毫米级小试样,须满足均匀应力假设
  • 输出完整动态应力-应变曲线
  • 适用:材料本构与应变率效应标定
落锤冲击试验
  • 应变率 10¹–10² s⁻¹,中低应变率
  • 大试样乃至真实结构件,几何保真
  • 输出冲击力-位移与吸收能量曲线
  • 适用:部件抗撞性与整车碰撞复现
材料本构选霍普金森杆,部件/整车抗撞性选落锤——两者不互替。
9第 9 页 · 冲击试验方法对比

试验设计优化策略

前页说明冲击试验手段很多,但“测得到”不等于“能代表实物”。优化设计应先回答:哪些物理机制必须相似,又怎样用最少的试验覆盖主要因素?

先定目标响应
先明确比较峰值力、吸能、失效模式或变形场;目标不同,相似条件也会改变。
无量纲匹配
用Buckingham Π定理组合变量,再按目标筛选几何、材料、惯性、速率、接触和热学量组。
正交组合
如L₉(3⁴)以9次试验估计4个三水平因素的主效应,而非穷举3⁴=81次。
估计主效应
用极差、方差或信噪比识别敏感因素;交互作用须预留列,并保证重复和误差估计。
缩小模型预测真实碰撞对应 →动态相似准则

模型与实物须成组匹配关键无量纲比值;只保持几何比例,不能保证响应等效。

Πdyn=ρv2E,Πε˙=ε˙Lv,Πm=mρL3\Pi_{\mathrm{dyn}}=\frac{\rho v^2}{E},\quad \Pi_{\dot{\varepsilon}}=\frac{\dot{\varepsilon}L}{v},\quad \Pi_m=\frac{m}{\rho L^3}
10第 10 页 · 试验设计优化策略

碰撞仿真有限元流程

从CAD到结果对标,五步构成碰撞仿真完整闭环。

1
几何与离散
CAD清理、中面抽取、网格划分,单元类型与尺寸决定精度与成本
2
材料本构
率相关硬化、失效准则如GISSMO,参数依赖试验标定
3
接触与载荷
自接触算法、摩擦模型、初始速度与约束边界
4
显式求解
时间步控制、沙漏抑制、能量监控与稳定性保障
5
验证迭代
能量平衡检查、与试验对标、参数敏感性迭代优化
11第 11 页 · 碰撞仿真有限元流程

显式与隐式算法对比

显式和隐式都能算冲击,但适用边界截然不同,选错代价很大。

显式算法
  • 时间步极小,约微秒级
  • 无需迭代,显式推进时间
  • 高速冲击、大变形首选
  • 条件稳定,受CFL限制
隐式算法
  • 时间步可大,毫秒级以上
  • 需迭代求解平衡方程
  • 准静态、屈曲等问题适用
  • 理论上无条件稳定
冲击碰撞必须用显式;准静态或屈曲等低速问题才考虑隐式。
12第 12 页 · 显式与隐式算法对比

典型材料卡片参数示例

LS-DYNA中常见材料模型的参数配置

典型材料卡片参数示例
LS-DYNA中常见材料模型的参数配置
13第 13 页 · 典型材料卡片参数示例

核心概念自测

点击作答

车身结构件碰撞工况应变率约 50–500 s⁻¹,若仅以 SHPB 试验数据建立材料卡进行仿真,下列判断最准确?

14第 14 页 · 核心概念自测

课后思考

先独立思考三分钟再看参考答案,重点体会思考路径而非结论本身。

1为什么准静态试验得到的材料参数不能直接用于碰撞仿真?

参考答案因为材料具有应变率敏感性——准静态约 0.001/s,碰撞可达 100–1000/s,屈服应力与断裂行为差异显著,必须通过 SHPB 等高速试验重新标定。

2设计吸能盒冲击试验时,应变率与三轴度应如何同时复现?

参考答案试件需逼近真实应变率并复现实际应力三轴度,常用带缺口 SHPB 试件调整应力状态,但要权衡端部摩擦、应力均匀性等次生效应。

3若材料局部出现相变或损伤累积,连续介质假设下的仿真还成立吗?

参考答案此时均匀连续假设失效,需引入内聚力模型(CZM)或相场断裂等非连续方法,结合微观–宏观跨尺度耦合,是当前研究前沿。

15第 15 页 · 课后思考