可压碎泡沫材料模型
可压碎泡沫材料的理论基础
泡沫材料的力学
老师,泡沫材料（发泡体）的力学与金属有何区别？
泡沫材料（EPS、PU泡沫、金属泡沫）在压缩时体积大幅减少。金属为非压缩塑性（$\Delta V = 0$）,而泡沫为压缩性塑性（$\Delta V \neq 0$）。通过壁面座屈和压碎吸收能量。
压缩应力-应变曲线
泡沫的典型压缩曲线：
1. 弹性域 — 壁面的弹性变形
2. 平台域 — 壁面的座屈/压碎。应力基本恒定,大变形
3. 密实化域 — 壁面相互接触。应力急剧上升
FEM中的建模
总结
泡沫模型诞生的动机
剑桥大学的Deshpande & Fleck在2000年发表的"Isotropic constitutive models for metallic foams"是为了量化铝制发泡金属（Alporas、Cymat）的冲击吸收特性而开发的。在此之前,von Mises模型完全无法表现泡沫由等静压压缩引起的大变形,因此需要允许体积塑性应变的新屈服面。该论文成为Abaqus"Crushable Foam"模型的理论基础。
可压碎泡沫材料的数值计算方法
泡沫的FEM设置
```
*MAT_LOW_DENSITY_FOAM
$ 压缩应力-应变表
*DEFINE_CURVE
0., 0.
0.1, 0.5
0.5, 0.6
0.8, 5.0
```
Abaqus: *CRUSHABLE FOAM + 表格输入。
总结
应力比参数k₀的实测方法
Crushable Foam的屈服面形状参数k₀（初始静水压屈服应力比）通过两种试验确定：单轴压缩试验和静水压压缩试验。住友电工生产的Alporas铝泡沫（密度0.25 g/cm³）的情况下,单轴压缩屈服应力≈1.6 MPa、静水压屈服≈1.9 MPa,由此得出k₀≈1.19,这是Deshpande本人的实验报告。在k₀不明确的情况下,实务上的惯例是先尝试k₀=1.1〜1.3。
可压碎泡沫材料的实务应用
泡沫的实务
包装（电子设备缓冲材）、汽车保险杠泡沫、头盔EPS内衬。
实务检查清单
汽车座椅垫圈冲击吸收
聚丙烯（PP）发泡泡沫（密度30〜60 kg/m³）用于头枕、膝垫等乘员保护部件。丰田、本田、大众等厂商从2010年代开始使用Abaqus的Crushable Foam模型进行虚拟试验,以符合FMVSS201U（乘用车内装冲击规范）标准。只需输入树脂发泡材的压缩-密实化曲线（应力~100%应变）,就能将冲击加速度-时间历程重现精度控制在±10%以内,各厂商都有相关报告。
可压碎泡沫材料的软件比较
泡沫的工具
LS-DYNA与Abaqus的发泡材实现差异
LS-DYNA的典型发泡材模型是MAT_LOW_DENSITY_FOAM（#57）,是一种准物理模型,可直接输入应变速率相关的加载-卸载曲线。而Abaqus的Crushable Foam基于屈服面的数学公式化,需要通过逆解析从试验曲线同定k₀、α等参数。业界的实务共识是,对于落体冲击等高速现象使用LS-DYNA MAT57,对于蠕变变形或长期疲劳使用Abaqus CDP系。
可压碎泡沫材料的先端研究
泡沫的先端
向各向异性泡沫的扩展
Deshpande-Fleck模型以等向性为前提,但挤出成形泡沫和纤维复合泡沫在厚度方向和面方向的屈服应力可相差2〜5倍。2006年,Chen和Lu提出了"Orthotropic Crushable Foam",引入三个方向独立屈服应力的定义方式。该公式已被实现为LS-DYNA MAT_ORTHOTROPIC_CRUSHABLE_FOAM（Material #142）,用于蜂窝板的OOB（面外压缩）分析。
可压碎泡沫材料的故障排除
泡沫的故障
密实化后的陡峭硬化与收敛问题
泡沫压缩应变超过70〜80%时开始密实化（densification）,应力呈指数增长。这种陡峭硬化常常阻碍增量分析的收敛。对策是在表格输入的最大应变处截断密实化区间,之后切换到线性外推。在Abaqus的Crushable Foam输入中,应力-应变表的最终斜率自动外推,因此在密实化开始点适当截断表格是解析稳定的关键。
"解析不符"时的对策
- 先深呼吸——仓促地随意改变设置只会使问题更复杂
- 制作最小再现案例——用最简单的形式重现可压碎泡沫材料模型的问题。"减法调试"最有效率
- 一次只改一项再执行——同时做多个改变会不清楚什么产生了效果。遵循科学实验的"对照实验"原则
- 回到物理本质——如果计算结果呈现出"东西无故浮起"这样非物理的结果,要怀疑输入数据的根本性错误
价值
详细
错误