聚氨酯泡沫铝动力学性能实验及本构模型研究*

张 勇，陈 力，陈荣俊,谢卫红

(1.中国矿业大学力学与建筑工程学院，江苏 徐州 221008； 2.空军勤务学院机场工程系,江苏 徐州 221000； 3.中国人民解放军理工大学国防工程学院,江苏 南京 210007)

开孔泡沫铝具有三维连通的孔结构，通过向其孔洞中填入泡沫聚氨酯等黏弹性高分子材料，能够起到良好的缓冲、吸能和减震作用。利用聚氨酯泡沫铝复合材料的抗冲击性能和吸能特性，可以有效地提高抵抗爆炸等冲击作用所产生的危害，对制造抗冲击防护材料具有重要意义。然而，要成为性能优异的冲击吸能材料[1-7]，需要知道泡沫聚氨酯材料在填充到泡沫铝孔隙中后其压缩变形行为和吸能性。

在泡沫铝材料和聚氨酯材料的实验研究方面，V.S.Deshpande等[8]通过实验得出泡沫铝的力学性能与应变率无关的结论；K.A.Dannemann等[9]采用常规分离式霍普金森压杆技术实验,发现闭孔Alporas泡沫铝具有明显的应变率效应；王永刚等[10]利用改进了的SHPB技术对泡沫铝进行了高应变率下的动态压缩实验，王志华等[11]利用SHPB实验技术和MTS材料试验机对开孔泡沫铝进行了准静态和动态压缩实验研究，卢子兴等[12]对高密度硬质聚氨酯泡沫塑料的拉伸力学性能进行了研究。在本构模型的研究中，主要成果有可用于金属大变形、高应变率和高温情况下的本构模型----Johnson-Cook模型[13]，J.A.Sherwood 等[14]提出了考虑环境温度和相对密度及应变率等因素的经验型本构关系等，A.G.Hanssen等[15]提出了泡沫铝在单轴受压或静水压力下的应变强化模型，胡玲玲等[16]提出了泡沫铝材料的一维黏塑性本构模型。以上研究成果主要是对泡沫金属和聚氨酯材料单独进行的。本文中对聚氨酯泡沫铝复合材料的动态力学性能进行较系统的实验研究，并通过考虑聚氨酯填充量对复合材料动力学性能的影响，提出适用于该复合材料的本构模型。

1 聚氨酯泡沫铝复合材料制备

采用高压渗流法制造开孔泡沫铝板材，并通过线切割将板材加工成直径74.00 mm的48个圆柱形薄试件待用；聚氨酯由白料(包括聚醚三元醇、发泡剂、催化剂、泡沫稳定剂等)和黑料(甲苯二异氰酸酯：PM200)按其先后顺序进行1∶1的质量比例混合(所得泡沫聚氨酯的密度约为0.3 g/cm3)以后，迅速将其倒入装有泡沫铝试件的密闭容器(将废旧塑料瓶进行改装后作为反应容器)中加压发泡(加压方法采用活塞手动加压)，使发泡产生的聚氨酯泡沫能尽可能多地填充到泡沫铝试件的孔隙之中。

经过打磨以后就得到聚氨酯泡沫铝的48个复合试件。经过40 ℃恒温固化和自然风干后，再称其质量，通过泡沫铝的质量增量来确定聚氨酯的填充量。为便于分析，将所有试件按相对密度的升序变化排列。制备的48个试件分成4组，每组12个试件。

2 动态冲击实验

动态荷载利用直径74 mm的直锥变截面大直径SHPB装置施加冲击荷载，由于应变率只与子弹速度相关，而弹速取决于气压，因此设计了由0.3、0.4、0.6 MPa 3种气压大致得到的1 000、1 600、1 900 s-1等3种应变率。图1是聚氨酯泡沫铝试件在应变率为1 900 s-1下冲击前后的形状对比。

图1 聚氨酯泡沫铝试件在应变率1 900 s-1下冲击前后的形状对比Fig.1 Shape contrast of polyurethane foam aluminum test-pieces before and after impact at strain rate of 1 900 s-1

图2 不同应变率下聚氨酯泡沫铝试件的应力应变曲线Fig.2 Stress-strain curves of polyurethane foam aluminum test-pieces at different strain rate

图2给出了聚氨酯泡沫铝在应变率为1 000、1 600、1 900 s-1等3种应变率下动态压缩时的应力应变曲线。从图2(a)中可以看到，低应变率的聚氨酯泡沫铝基本上只表现出了2个变形阶段的特征，即弹性段和应力上升的塑性段，没有出现明显的密实段，这是由于撞击杆的速度较低，撞击能量不足以将聚氨酯泡沫铝压实，所以压缩曲线不完整。

应变率相对较高的图2(b)～(c)中的应力应变曲线却表现出了线弹性段、应力上升的塑性段和密实段3个变形阶段的特征，同种应变率下线弹性阶段基本都是重合的，说明复合材料的动态弹性模量与相对密度无关；聚氨酯泡沫铝的塑性段曲线有所抖动，而且随着相对密度和应变率的上升波动减弱，这说明抖动除了泡沫铝的非线性行为以外，还有聚氨酯的影响因素；另外由图2可以看出，聚氨酯泡沫铝的屈服强度和流变应力随应变率的增加而提高。

由图3可知，泡沫聚氨酯质量的增加量与屈服强度的提高呈近似线性关系。

图3 屈服强度与聚氨酯质量分数的关系Fig.3 The yield strength vs the mass of the polyurethane foam

3 聚氨酯泡沫铝动态本构模型的建立

文献[14]在建立聚氨酯泡沫本构关系时将环境温度和相对密度及应变率等因素考虑进去，提出了经验型本构关系：

(1)

胡时胜等[17]对本构关系(1)进行了修正，在修正的本构关系中不包含温度，关于应变率对应力的影响则采用基于热激活机制的Seeger模型，得到拟合的本构关系：

(2)

图4 屈服强度与聚氨酯质量分数的数据拟合Fig.4 Data fitting of yield intensity vs mass of polyurethane

图4是聚氨酯质量分数与屈服强度的关系曲线图。根据图4，本构关系的影响可近似取线性关系表示为：

P(m)=Hm+k

式中：m为聚氨酯含量，H、k为常数。

由本文实验数据进行分析可知，相对密度对应力的影响近似为线性关系，与本构关系(1)结论一致，可表示为：

G(ρ)=J(ρ-ρ0)+1

式中：J为常数。

关于应变率对应力的影响仍采用本构关系(2)中的关系。最终可得拟合的本构方程为：

(3)

由图4拟合直线可得：P(m)=63.2m+13.723。相关系数R1=0.929 96，拟合效果很好。

y=-18 860x4+21 117x3-6 818.4x2+913.79x

由相关系数R2=0.962 9可知，其拟合性比较好。

则得聚氨酯泡沫铝的本构方程为：

(-18 860ε4+21 117ε3-6 818.4ε2+913.79ε)

(4)

下面分析聚氨酯含量、应变率和相对密度不同时对拟合的本构模型应力应变曲线的影响并与实验数据对比。

由图5可知，聚氨酯泡沫铝的本构模型(4)在相对密度为0.274～0.293之间，应变率在1 000～1 900 s-1之间，其拟合结果与实际应力应变曲线吻合较好，特别是相对密度越接近0.293，应变率越接近1 900 s-1，其拟合效果越好。

由图5(e)可知，试件36的实验条件都处于相对密度和应变率适用范围的下限，所以拟合效果就比较差，但是在线弹性阶段和进入塑性阶段初期，其拟合效果比较好，塑性阶段进入密实阶段，拟合效果开始变差。

图5 本构模型曲线与实验数据曲线对比Fig.5 Constitutive model vs experimental data curves

4 结 论

通过对聚氨酯泡沫铝冲击实验和动态本构模型的分析，可以得到以下结论：

(1)聚氨酯泡沫铝的动态弹性模量与相对密度无关。

(2)屈服强度和流变应力与应变率和泡沫铝的相对密度成正比。

(3)拟合的本构模型在相对密度为0.274～0.293之间，应变率在1 000～1 900 s-1之间，拟合结果与实际应力应变曲线吻合较好，特别是相对密度越接近0.293，应变率越接近1 900 s-1，拟合效果越好；实验条件都处于相对密度和应变率适用范围的下限时，拟合效果开始变差，说明本构方程的使用有一定的局限性。

[1] Li Xu-hua, Yuan Qiao-long, Wang De-ning, et al.Cross-linked polyurethane with high thermal stability and low optical loss[J].Materials and Devices for Optical and Wireless Communications, 2006,49(5):405-408.

[2] Sun Ming-cheng, Liu Fu-chun, Shi Hong-wei, et al.A study on water absorption in freestanding polyurethane films filled with nano-TiO2pigments by capacitance measurements[J].Acta Metallurgica Sinica, 2009,22(1):27-34.

[3] Li Hai-bin, Tian Zhen, Zhang Ai-ying, et al.Synthesis and characterization of novel triblock copolymers comprising poly (tetrahydrofuran) as a central block and poly(γ-benzyl L-glutamate)s as outer blocks[J].Frontiers of Materials Science in China, 2008,2(1):84-90.

[4] You Cao, Yu Jiang, Zhao Shu-lu, et al.Effect of montmorillonite on kinetics of polyurethane preparation reaction[J].Chinese Chemical Letters, 2008,19(1):115-118.

[5] 穆建春,习会峰,龙志勤.不同孔隙率及孔径泡沫铝的力学与吸能特性研究[J].实验力学,2009(3):223-227.Mu Jian-chun, Xi Hui-feng, Long Zhi-qin.Study on the mechanics and energy absorb property of aluminum foam in different void rate and different aperture[J].Experiment Mechanics, 2009(3):223-227.

[6] Zhang Guo-qi, Wang Bing, Ma Li, et al.Energy absorption and low velocity impact response of polyurethane foam filled pyramidal lattice core sandwich panels[J].Composite Structures, 2014,108:304-310.

[7] 赵军,侯兴勃,张丹.泡沫铝冲击吸能器工程设计研究[J].噪声与振动控制,2009,29(3):130-131,151.Zhao Jun, Hou Xing-bo, Zhang Dan.Engineering design and study on the impact and energy absorb implement of aluminum foam[J].Yawp and Vibraten Control, 2009,29(3):130-131,151.

[8] Deshpande V S, Fleck N A.High strain rate compressive behaviour of aluminium alloy foams[J].International Journal of Impact Engineering, 2000,24(3):277-298.

[9] Dannemann K A, Lankford J J.High strain rate compression of closed-cell aluminium foams[J].Material Science and Engineering: A, 2000,293(1/2):157-164.

[10] 王永刚,胡时胜,王礼立.爆炸荷载下泡沫铝材料中冲击波衰减特性的实验和数值模拟研究[J].爆炸与冲击,2003,23(6):516-522.Wang Yong-gang, Hu Shi-sheng, Wang Li-li.Shock attenuation in aluminum foams under explosion loading[J].Explosion and Shock Waves, 2003,23(6):516-522.

[11] 王志华,曹晓卿,马宏伟,等.泡沫铝合金动态力学性能实验研究[J].爆炸与冲击,2006,26(1):46-52.Wang Zhi-hua, Cao Xiao-qing, Ma Hong-wei, et al.Experimental studies on the dynamic compressive properties of open celled aluminum alloy foams[J].Explosion and Shock Waves, 2006,26(1):46-52.

[12] 卢子兴,高镇同,朱汪鲲,等.聚氨酯泡沫塑料的强度与断裂韧性[J].力学学报,2000,32(5):625-632.Lu Zi-xing, Gao Zhen-tong, Zhu Wang-kun, et al.Investigation on the fracture properties of polyurethane rigid foam plastics[J].Mechanics Transaction, 2000,32(5):625-632.

[13] Johnson G R，Cook W H.A constitutive model and data for metals subjected to large strain, high strain rates and high temperatures∥Proceedings of the 7th International Symposium on Ballistics[C].The Netherlands: The Hague, 1983:541-547.

[14] Sherwood J A, Frost C C.Constitutive modeling and simulation of energy absorbing polyurethane foam under impact loading[J].Polymer Engineering & Science, 1992,32(16):1138-1146.

[15] Hanssen A G, Hopperstad O S, Langseth M, et al.Validation of constitutive models applicable to aluminium foams[J].International Journal of Mechanical Sciences, 2002,44(2):359-406.

[16] 胡玲玲,黄小清,张红,等.泡沫铝材料一维粘塑性本构关系[J].华南理工大学学报,2004,32(4):87-90.Hu Ling-ling, Huang Xiao-qing, Zhang Hong, et al.One dimension viscoplasticity constitutive relationship of aluminum foam material[J].South China Science and Engineering Transaction, 2004,32(4):87-90.

[17] 胡时胜,刘剑飞,王悟.硬质聚氨酯泡沫塑料本构关系的研究[J].力学学报,1998,30(2)：151-156.Hu Shi-sheng, Liu Jian-fei, Wang Wu.Study on the constitutive relationship of horniness polyurethane foam[J].Mechanics Transaction, 1998,30(2):151-156.