冲击载荷下径向密度排布对泡沫金属力学性能影响的研究*

王根伟，刘 冕，宋 辉，王 彬

（1. 太原理工大学机械与运载工程学院应用力学研究所，山西 太原 030024；2. 太原理工大学材料强度与结构冲击山西省重点实验室，山西 太原 030024；3. 伦敦布鲁内尔大学机械航空工程系，英国 伦敦 UB8 3PH）

材料与结构的轻量化是科技发展的重要趋势，其中泡沫金属材料由于轻质化、孔隙率高和可变形量大等，被广泛应用于能量吸收和抗冲击减振的保护结构。为了提高泡沫材料的力学性能，研究人员提出了具有密度梯度的泡沫材料[1-3]，以实现多目标优化设计。

由于材料制造技术的局限性，在实际应用中还不容易获得具有特定相对密度分布的梯度泡沫样本，也无法直接观察加载情况下胞元内部的变形情况，研究人员转而应用有限元方法研究梯度泡沫材料的力学响应[4-9]。研究结果显示，冲击载荷下梯度泡沫材料存在3 种变形模态[10-12]，且具有特定密度梯度的泡沫金属可以降低受保护结构的最大冲击应力，改善材料的能量吸收和抗冲击性能[13-19]。Zeng 等[20]测试了密度梯度空心球金属泡沫的动态响应，将最硬层作为第1 冲击层，将最弱层与受保护结构接触，可以减少传输力并获得更多的吸收能量。Ajdari 等[21]研究了具有规则和不规则排列的功能梯度多胞结构的动态压溃，表明沿加载方向降低相对密度，可以增强蜂窝在压溃初期的能量吸收。Zhang 等[22]建立了分层递变梯度蜂窝模型，建议把最弱的层放置在冲击端或输出端，最强的层被放置在中间阶段，以达到更高的能量吸收效率。Zheng 等[23]从理论和数值计算两方面发现，均匀泡沫的能量吸收效果可能最好，密度梯度较大的正、中高梯度泡沫可以在撞击端保护物体，密度梯度较小的负梯度泡沫可以在支撑端保护物体，而大梯度的中低、中高、正梯度泡沫可以提高吸能效率。

目前，对梯度泡沫金属材料的研究主要集中在密度单向变化，为了进一步研究梯度泡沫金属材料密度在非单向变化下的变形及能量吸收机理，我们设计了一种密度径向连续排布的梯度泡沫材料，研究它在冲击载荷作用下的动态力学性能。首先，利用3D-Voronoi 技术构建密度排布不同的层状、径向梯度泡沫结构；然后，导入有限元分析软件中生成梯度泡沫金属的冲击模型，观察他们在不同恒定冲击速度下的变形特点，比较他们的抗冲击性能和能量吸收能力，考察不同梯度泡沫材料在工程防护上的优劣；最后，分析密度梯度和相对密度参数对径向梯度泡沫材料力学性能的影响。

1 数值模拟

1.1 三维随机Voronoi 技术

3D-Voronoi 模型是将N 个成核点完全随机分布在体积为V0的立方体区域内生成，为了防止Voronoi模型产生过小胞元，立方体内任意两相邻核点之间的最小距离为：

式中：无量纲变量k 为不规则度，d0为正十四面体模型中任意两相邻成核点间的最小距离[24]。将成核点复制到周围区域内，编程生成Voronoi 构型，并通过切割算法得到原先指定区域中的Voronoi 结构作为最后的计算模型。

1.2 密度梯度

为了研究密度排布对泡沫金属力学性能的影响，基于密度沿着厚度方向上单向变化的层状梯度泡沫模型(见图1(a))，设计了一种径向梯度泡沫模型，即以正方形中心点出发密度沿着半径方向渐进变化，如图1(b)所示。

层状梯度泡沫的密度梯度γ 可描述为：

由此，推演径向梯度泡沫密度梯度γ 的定义为：

图1 密度梯度泡沫模型Fig. 1 Density gradient foam models

曹国剑等[25]介绍了利用离心沉积法制备梯度泡沫铝的方法，可直接用来生产本文中提出的径向负梯度泡沫。

1.3 有限元模型

表1 为层状正、负梯度与径向正、负梯度4 种泡沫金属模型的材料参数。4 种泡沫金属的冲击加载模型如图2 所示，垂直于所示平面方向泡沫密度均匀，冲击加载沿梯度方向进行。模型尺寸为50 mm×50 mm×25 mm，不规则度为0.3，有限元模型的核点个数为900，胞壁厚度由给定的相对密度确定。在有限元收敛分析后，泡沫模型使用S3R 和S4R 壳单元进行混合网格划分，网格尺寸设置为0.2 mm。在左端固定一个刚性板，右端的刚性板以给定的恒定速度撞击泡沫模型。基体材料设为双线性应变强化模型，密度ρ=2 700 kg/m3，弹性模量E=69 GPa，泊松比ν=0.3，屈服应力σy=76 MPa，切线模量为0.69 GPa。泡沫模型定义为自接触，两端刚性板和泡沫模型之间定义为面面接触，设定动摩擦因数为0.2。

表1 模型材料参数Table 1 Model material parameters

图2 层状和径向的梯度泡沫Fig. 2 Layered and radial graded foams

2 结果与讨论

2.1 变形特点和名义应力应变曲线

根据梯度泡沫在不同冲击速度下的准静态模态、过渡模态和冲击模态[26-27]，将冲击速度设置为30、80 和200 m/s。其中，层状正、负梯度泡沫的变形特点与以往文献描述的相似。图3 为径向梯度泡沫在不同冲击速度下的应变云图。

图3 应变云图Fig. 3 Strain distributions

当冲击速度较小 (v=30 m/s) 时，变形属于准静态模态。径向正梯度泡沫先在冲击端变形（低密度），随后支撑端（低密度）的胞元发生变形。这是由于，当速度相对较低时，泡沫内部处于均匀的应力状态，两端首先达到准静态屈服应力。正梯度泡沫中间区域的密度较高，且模型结构并不完全对称，在ε＞0.5 后出现沿支撑板两个切线方向滑移现象。径向负梯度泡沫因中间密度较小，变形先从中间区域处产生，两端在后期同样出现滑移现象。

当冲击速度中等(v=80 m/s)时，变形受准静态屈服应力和惯性效应两个因素的影响。径向正梯度泡沫的变形特点属于典型的过渡模态，而径向负梯度泡沫没有表现明显的过渡模态。负梯度泡沫先在冲击端发生变形，随后中间区域变形，并逐渐向支撑端传播，此现象与冲击模态类似。发现在冲击速度为60 m/s 时，径向负梯度泡沫处于过渡模态。冲击端先发生少量变形，接着中间区域持续发生变形，待中间区域压实后，变形继续向支撑端传播，直到支撑端完全压实。从变形特点可以看出，径向负梯度泡沫比径向正梯度泡沫在更低速度下进入冲击模态。

当冲击速度较高(v=200 m/s)时，惯性作用成为了控制变形的主导因素，其影响效果远超过准静态屈服应力，泡沫进入冲击模态。径向正梯度泡沫在冲击端发生大量变形后，仅在支撑端发生少量变形；而径向负梯度泡沫只在冲击端发生变形，并逐渐变形至压实。

图4 为不同冲击速度下4 种泡沫冲击端与支撑端的名义应力应变曲线。随着冲击速度的增大，4 种泡沫冲击端应力与压实应变逐渐增大，支撑端比冲击端越来越晚达到初始屈服应力，且应力平台区越来越长，应力值差异逐渐明显。层状正梯度泡沫的冲击端初始峰值应力较小，到达峰值应力后曲线逐渐上升，但支撑端曲线应力值高于层状负梯度泡沫，初始屈服应力较大。层状负梯度泡沫的冲击端初始峰值应力较大，到达峰值应力后曲线逐渐下降。层状梯度泡沫的数值模拟结果与其他文献描述相似，佐证了有限元模型的合理性。

图4 4 种梯度泡沫冲击端和支撑端的名义应力应变曲线Fig. 4 Nominal stress-strain curves of impact end and support end of four graded foams

与其他3 种泡沫相比，径向正梯度泡沫由于两端密度小、相对较“软”，冲击端初始峰值应力和支撑端初始屈服应力均最低。径向正梯度泡沫中间区域密度大、强度高，冲击端曲线在应变0.2～0.6 区间时先上升后下降。径向负梯度泡沫因两端密度大、中间区域密度小，冲击端曲线在应变0.2～0.6 区间时先下降后上升，而支撑端曲线始终保持基本稳定。

图5 为4 种梯度泡沫冲击端与支撑端的最大应力和标准差。最大应力表示泡沫达到密实化应变前冲击端或支撑端所承受的最高载荷；标准差表示数据之间的离散程度，离散程度越小，曲线波动幅度越小。通过比较发现：冲击速度越高，4 种泡沫冲击端的最大应力与波动幅度呈整体增加趋势；径向正梯度泡沫冲击端的最大应力在4 种泡沫中最小；径向负梯度泡沫的两端应力波动幅度明显小于其他3 种，支撑端应力尤其稳定。

图5 4 种梯度泡沫冲击端和支撑端的最大应力和应力标准差Fig. 5 Maximum stresses and stress standard deviations of impact end and support end of four graded foams

2.2 能量吸收能力

能量吸收能力是判断泡沫材料抗冲击性能的一个典型指标。可以用材料压缩至某个变形量时所吸收的能量来衡量吸收能力的强弱，定义为：

图6 和表2 为4 种梯度泡沫能量吸收能力，可以看出，几种泡沫的能量吸收能力在变形过程中呈现交替变化的现象。

当速度为30 m/s 时：在压缩初期，径向负梯度泡沫表现出优异的能量吸收能力；在压缩后期，径向负梯度泡沫与层状正、负梯度泡沫几乎相同，而径向正梯度泡沫能量吸收值相对较小。当速度提高到80 和200 m/s 时：在压缩初期，层状、径向负梯度泡沫的能量吸收能力大于层状、径向正梯度泡沫；在压缩后期，径向负梯度泡沫的能量吸收值仍较高，层状正梯度泡沫的能量吸收能力增长加快。

图6 不同冲击速度下4 种梯度泡沫的能量吸收能力Fig. 6 Energy absorption capacities of four graded foams at different impact velocities

表2 4 种梯度泡沫的能量吸收能力Table 2 Energy absorption capacities of four graded foams

比较4 种梯度泡沫的力学响应发现，径向正梯度泡沫不仅拥有层状正梯度泡沫冲击端应力较小、保护冲击端冲击物体的特点，而且继承了层状负梯度泡沫支撑端应力较小、保护支撑端受冲击物体的特点。径向负梯度泡沫的两端应力波动幅度较小，曲线变化趋势较为平缓，能够保证被保护物体受力稳定，同时它的能量吸收能力始终处于较高水平，在低速冲击或大变形工况下是最优的吸能结构。因此，径向梯度泡沫是缓冲减压吸能装置的良好选择。

2.3 径向密度梯度的影响

为了分析径向正、负梯度泡沫对密度梯度的敏感程度，图7～8 分别给出了0.4、0.8 时的应力应变曲线和能量吸收曲线。

从图7 可以看出：对于冲击端，当密度梯度改变时，泡沫初始峰值应力和密实化时的应力变化明显；对于支撑端，径向正梯度泡沫在两种密度梯度下的曲线出现分离，而负梯度泡沫曲线几乎相同。对于径向正梯度泡沫，γ 越大，冲击端初始峰值应力与支撑端初始屈服应力越小；对于径向负梯度泡沫，当增加，冲击端初始峰值应力与支撑端初始屈服应力增大。

图7 在不同密度梯度下径向梯度泡沫的名义应力应变曲线Fig. 7 Nominal stress-strain curves of radial graded foams with different density gradients

从图8 可以看出，泡沫的能量吸收能力对密度梯度大小不敏感、对梯度方向敏感。当密度梯度改变时，两种梯度下的泡沫所吸收的能量几乎相同；在达到完全密实化前，径向负梯度泡沫的能量吸收能力始终大于径向正梯度泡沫。由于试样具有相同的平均相对密度，不同梯度泡沫完全密实化后所吸收的总能量基本相同。

图8 在不同密度梯度下径向梯度泡沫的能量吸收能力Fig. 8 Energy absorption capacities of radial graded foams with different density gradients

2.4 平均相对密度的影响

图9 在不同相对密度下径向梯度泡沫的名义应力应变曲线Fig. 9 Nominal stress-strain curves of radial graded foams with different relative densities

图10 在不同相对密度下径向梯度泡沫的能量吸收特性Fig. 10 Energy absorption per mass of radial graded foams with different relative densities

3 结 论

利用3D-Voronoi 技术设计了层状正、负梯度与径向正、负梯度4 种密度排布不同的梯度泡沫模型，通过有限元分析软件分别从变形特点、载荷响应历程及能量吸收机理比较了4 种梯度泡沫金属在恒速冲击载荷下的动态响应结果，同时探讨了径向梯度泡沫中密度梯度、相对密度参数对力学性能的影响，得到如下结论。

（1）径向正梯度泡沫与层状正、负梯度泡沫相比冲击端与支撑端的应力值均处于较低水平，不仅保护冲击端冲击物体，而且保护支撑端受冲击物体；径向负梯度泡沫两端应力应变曲线平缓，波动幅度较小，能够保证被保护物体受力稳定。

（2）在不同冲击速度下，4 种梯度泡沫金属的能量吸收能力交替变化，其中径向负梯度泡沫能量吸收能力始终处于较高水平，在低速冲击或大变形工况下是最优的吸能结构。

（3）能量吸收能力对密度梯度大小不敏感、对梯度方向敏感，在达到完全密实化前，径向负梯度泡沫的能量吸收能力始终大于径向正梯度泡沫；平均相对密度越大，径向正、负梯度泡沫两端应力越大、吸能效果越好。