高应变率下航空透明聚氨酯的动态本构模型

张龙辉，张晓晴,姚小虎，臧曙光

(1.华南理工大学土木与交通学院，广东 广州 510640； 2.中国建筑材料检验认证中心有限公司，北京 100024)



高应变率下航空透明聚氨酯的动态本构模型

张龙辉1，张晓晴1,姚小虎1，臧曙光2

(1.华南理工大学土木与交通学院，广东 广州 510640； 2.中国建筑材料检验认证中心有限公司，北京 100024)

采用低阻抗分离式霍普金森压杆对航空透明聚氨酯进行了高应变率下的动态力学性能测试，得到的应力应变曲线表现出了显著的非线性黏弹性特征。基于本构理论和实验数据，构建了航空透明聚氨酯的松弛时间应变相关的超黏弹性本构形式。该本构模型由2部分组成 ：一部分表征准静态下的超弹性行为，另一部分描述非线性应变率的相关特性。利用超黏弹性本构模型对不同应变率下航空透明聚氨酯的动态应力应变曲线进行拟合，拟合曲线与实验曲线一致性良好。

固体力学；本构模型；SHPB；航空PU；高应变率；超黏弹性

近年来，随着科技的发展，一大批新型工程材料不断涌现出来，随着对安全与防护的日益重视，橡胶和泡沫等低阻抗材料作为吸能缓冲材料被广泛应用在交通工具、管道、钢结构和军用设备上。因此，研究这类材料在各种冲击载荷下的动态力学行为已成为设计部门日益关注的问题。聚氨酯(简称PU)作为橡胶的替代物，是用途最广的材料之一。弹性体动态力学性能是当今国际上高聚物弹性体力学性能研究的前沿问题，C.M.Roland等[1]设计了一个测试弹性体高达500 s-1应变率下的动态拉伸力学性能的滑轮冲击装置，描述了聚脲在应变率14～573 s-1下应变达到300%的动态拉伸行为。J.Yi等[2]研究了聚脲和三种聚氨酯在低应变率和高应变率下的动态压缩力学行为。S.S.Sarva等[3]研究了聚脲和聚氨酯从0.001 s-1到10 000 s-1宽泛应变率下的单轴压缩应力应变行为，依次在每一个应变率数量级下进行实验，聚脲和聚氨酯经过了低应变率下的橡胶态到高应变率下的玻璃态的转变。J.Shim等[4]结合波反褶积技术，采用改装SHPB装置进行了聚脲在中高应变率下的动态压缩实验，并且将改装的SHPB装置跟传统的SHPB装置以及万能液压器下的实验进行了对比，验证低应变率和高应变率下该实验的有效性。A.V.Amirkhizi等[5]研究了基于实验的考虑压力和温度效应的聚脲的线性黏弹性本构行为，该模型综合考虑了热力学能量耗散机理，WLF时温转换和压力相关性。C.Li等[6]建立了聚脲的超黏弹性本构关系，其中超弹性部分采用的是Ogden模型，非线性黏弹性部分采用非线性黏弹性模型，根据剪切模量的实验数据确定松弛时间。

航空聚氨酯胶片是用于风挡夹层玻璃的中间聚合物膜，服役温度广，基本满足航空透明件的实际服役温度；跟普通的聚氨酯相比透明性非常好，在紫外线照射下不会发黄；同时具有承受大变形和吸收冲击能的作用。透明聚氨酯胶片作为黏结层或中间层将被广泛应用于夹层透明材料的生产中，例如航空透明件(风挡玻璃、舷窗玻璃等)、高铁风挡玻璃等。上述结构和部件在使用的过程中有可能承受如高速撞击、爆炸等冲击载荷的作用，但是高性能透明聚氨酯胶片在冲击荷载下的动态力学行为和本构模型鲜有报道。因此对透明聚氨酯材料在高应变率动态加载下响应特性的研究十分重要，这将为飞机和高铁挡风玻璃的设计提供正确的材料模型和可靠的材料数据。

本文中将研究航空透明PU在高应变率作用下的动态力学性能。用低阻抗Hopkinson分离式压杆在应变率为2 000～7 000 s-1范围内对材料进行力学性能测试。同时结合已有的超黏弹性本构模型[7-8]构建航空透明聚氨酯的本构模型，并拟合实验结果、给出拟合参数。

1 实验装置

测试材料为美国PPG 公司提供的S-123 新型透明聚氨酯胶片，密度1 180 kg/m3。这种高性能航空级别透明聚氨酯胶片具有优异的动态力学性能、耐热性、耐化学品性、耐光降解性及耐候性。实验是在低阻抗Hopkinson分离式压杆装置上进行的。由于这种聚合物材料中的波阻抗低，实验中考虑了材料在SHPB实验中的试样尺寸设计[9]，将试样加工为长2 mm、直径4 mm的圆柱体。采用低阻抗的铝杆SHPB装置，其参数为：撞击杆的尺寸为∅14 mm×400 mm；入射杆的尺寸为∅14 mm×1 000 mm，透射杆的尺寸为∅14 mm×1 000 mm。总共做了18组实验，打击速度为6.28～17 m/s。在试样两端及压杆端面均匀地涂抹一层薄薄的凡士林，尽量消除端面摩擦力的影响。

2 聚氨酯动态实验结果分析

实验中实测到的入射波、反射波及透射波原始波形的典型曲线如图1所示。图2显示了多组航空透明聚氨酯材料高应变率下的真实应力应变曲线，该材料表现出了明显的应变率效应。由不同打击速度(6.28～17.00 m/s)下的实验发现,随着打击速度的提高,航空聚氨酯材料的应变率范围为2 500～6 500 s-1。对于航空聚氨酯这种高聚物材料，很难给出试样的弹性模量、屈服强度等参数。定义应变为0.15处对应点到原点连线的斜率来表示杨氏模量Eg，称之为割线模量[10]。为了比较航空聚氨酯的应变率效应,取每条曲线真实应变0.15 处对应的平台应力值σs进行比较。随着加载应变率的提高，应力值σs从22.73MPa增加到35.51MPa。另外，该材料的割线模量随着应变率的增加而提高，从2 500s-1下的151.53MPa增加到6 500s-1时的236.67MPa。

图1 原始波形Fig.1 Typical original waves

图2 航空聚氨酯SHPB压缩实验应力应变曲线Fig.2 True stress-strain curves of the aviation polyurethane

3 超黏弹性本构模型

典型高聚物的非线性黏弹性本构行为可以用ZWT本构模型[11]来描述，它由一个非线性弹簧，一个低频Maxwell体和高频Maxwell体三者并联所组成。在冲击载荷下，ZWT本构模型的低频Maxwell单元没有足够时间来松弛，所以由一个松弛时间可以描述材料在高应变率下的行为。航空透明聚氨酯的实验中，即使压缩应变达到55%以上，卸载后试件也基本恢复了原状而没有残余变形。借鉴ZWT本构模型的思路和方法，以及C.Li等[6]建立的聚脲的超黏弹性本构，和L.M.Yang[7-8]提出的高应变率下的超黏弹性本构模型，本文构建了描述航空透明聚氨酯在高应变率下的超黏弹性本构模型，如图3所示的并联模型。

图3 A和B单元的并联本构模型Fig.3 Parallel mechanical elements A and B

3.1 超弹性模型

根据R.S.Rivlin[12]的分析，各向同性不可压缩超弹性材料的本构关系可以用左Cauchy-Green变形张量的不变量形式表示为：

(1)

材料在单轴压缩载荷下，在加载方向上的伸长比计为λ，则3个主伸长比λ1=λ，λ2=λ3=λ1/2。变形梯度F和左Cauchy-Green变形张量B可由下式表示：

(2)

(3)

基于航空透明聚氨酯准静态实验的结果，应变势能的形式取为3项Mooney-Rivlin模型的应变能函数去拟合实验数据：

图4 准静态压缩实验曲线与模型拟合曲线的比较Fig.4 Comparison between quasi-static curves of experimental data and proposed model

(4)

所以，由式(3)～(4)得到

(5)

静态实验设备为INSTRON5567电子万能材料试验机。试样长为2 mm、直径为4 mm。实验中采用位移加载控制方式，应变率为0.001 s-1。采用最小二乘法数据拟合技术，方程(4)中的材料参量C10、C01、C11由其准静态实验应力应变数据拟合确定，图4表明:模型的响应和实验曲线是一致的。

3.2 黏弹性本构模型

黏弹性材料最大的特点就是力学性能受之前变形的影响，其应力依赖于应变和应变率历史，如聚合物，该类材料的响应特征是衰减惯性。对于不可压缩各向同性材料的黏弹性本构模型通常用下式来示：

(6)

式中：σv是Cauchy真应力张量，pv是静水压力，反映材料的不可压缩性。Ω为本构泛函，反映应变历史对应力的影响，为了用较少的参数描述材料有限变形下的黏弹性力学行为，将Ω近似取为：

(7)

其中应变率为：

(8)

这种积分形式中的松弛函数通常被定义为一个随时间t衰减的松弛函数，这里假设为

(9)

式中：θi为松弛时间。基于高应变率下ZWT模型，本文中取N=1，用来描述高应变率时的情况，而不考虑低应变率范围的应变率敏感性。

聚合物的应力松弛中，不同分子链条的重新调整对应着不同的松弛时间。在材料的局部区域，松弛包括短链的快速恢复，而长链的重新调整需要更多的松弛时间[13]。对应不同的应变，弹性体包含了宽泛的链条长度，从而对应不同的松弛时间。作用较小的应变时对应较短的链条的快速恢复从而需要较短的松弛时间，而作用较大的应变则包括了不同范围的链条长度的调整从而对应较长的松弛时间。

因为松弛时间是随着应变增加而增加的，这里假设松弛时间的形式为

(10)

本文中松弛时间取为最简单的形式

(11)

(12)

(13)

(14)

将准静态本构方程与应变率相关项相加得到航空聚氨酯的超黏弹性本构方程：

(15)

(16)

因此高应变率单轴加载下的应力可以表示为

(17)

由横向应力σ22=σ33=0，得到

(18)

(19)

图5 高应变率压缩实验曲线与模型拟合曲线的比较Fig.5 Comparison between high strain rates of experimental data and proposed model

表1 由实验数据拟合确定的模型参量

4 结 论

利用低阻抗的SHPB装置对航空透明聚氨酯进行了动态压缩测试，结果表明，测试的航空聚氨酯材料的动态压缩行为具有明显的应变率相关性，且表现出一定的黏弹性特征。应用相关本构理论，构建了描述航空透明聚氨酯动态压缩力学性能的率相关本构模型，该模型由描述准静态响应的Mooney-Rivlin超弹性模型和描述高应变率响应的非线性黏弹性模型组成，理论曲线与实验得到的应力应变曲线有很好的一致性,说明提出的模型能够很好地描述航空PU的力学性能，为航空挡风玻璃的高性能透明聚氨酯胶片材料的研发和设计奠定了基础。

[1] Roland C M, Twigg J N, Vu Y, et al. High strain rate mechanical behavior of polyurea[J]. Polymer, 2007,48(2):574-578.

[2] Yi J, Boyce M C, Lee G F, et al. Large deformation rate-dependent stress-strain behavior of polyurea and polyurethanes[J]. Polymer, 2006,47(1):319-329.

[3] Sarva S S, Deschanel S, Boyce M C, et al. Stress-strain behavior of a polyurea and a polyurethane from low to high strain rates[J]. Polymer, 2007,48(8):2208-2213.

[4] Shim J, Mohr D. Using split Hopkinson pressure bars to perform large strain compression tests on polyurea at low, intermediate and high strain rates[J]. International Journal of Impact Engineering, 2009,36(9):1116-1127.

[5] Amirkhizi A V, Isaacs J, Mcgee J, et al. An experimentally-based viscoelastic constitutive model for polyurea, including pressure and temperature effects[J]. Philosophical Magazine, 2006,86(36):5847-5866.

[6] Li C, Lua J. A hyper-viscoelastic constitutive model for polyurea[J]. Materials Letters, 2009,63(11):877-880.

[7] Yang L M, Shim V, Lim C T. A visco-hyperelastic approach to modelling the constitutive behaviour of rubber[J]. International Journal of Impact Engineering, 2000,24(6):545-560.

[8] Pouriayevali H, Guo Y B, Shim V. A constitutive description of Elastomer behaviour at high strain rates-A strain-dependent relaxation time approach[J]. International Journal of Impact Engineering, 2012,47:71-78.

[9] Chen W, Zhang B, Forrestal M J. A split Hopkinson bar technique for low-impedance materials[J]. Experimental Mechanics, 1999,39(2):81-85.

[10] 林玉亮,卢芳云,卢力.高应变率下硅橡胶的本构行为研究[J].高压物理学报,2007,21(3):289-294. Lin Yu-liang, Lu Fang-yun, Lu Li.Constitutive behaviors of a silicone rubber at high strain rates[J]. Chinese Journal of High Pressure Physics, 2007,21(3):289-294.

[11] 王礼立.冲击动力学进展[M].合肥：中国科学技术大学出版社，1992.

[12] Rivlin R S. Collected papers of R. S. Rivlin[M]. Berlin： Springer, 1997.

[13] Hoo Fatt M S, Xin O. Integral-based constitutive equation for rubber at high strain rates[J]. International Journal of Solids and Structures, 2007,44(20):6491-6506.

(责任编辑 曾月蓉)

Constitutive model of transparent aviation polyurethane at high strain rates

Zhang Long-hui1, Zhang Xiao-qing1, Yao Xiao-hu1, Zang Shu-guang2

(1.SchoolofCivilEngineeringandTransportation,SouthChinaUniversityofTechnology,Guangzhou510640,Guangdong,China; 2.ChinaBuildingMaterialTestandCertificationCenter,Beijing100024,China)

The uniaxial compressive properties of aviation polyurethane were investigated experimentally by using a modified aluminum split Hopkinson pressure bar apparatus. The obtained stress-strain curves presented distinct non-linear viscoelastic characteristic.Based on the constitutive theory and the experimental data, a hyper-viscoelastic constitutive model that incorporated a strain-dependent relaxation time was proposed to describe the large compressive deformation response of incompressible aviation polyurethane at high strain rates. The proposed model was made up of two parallel mechanical elements-one component to characterize quasi-static hyperelastic behavior, and the other to define rate-sensitivity and strain history dependence. The predictions of the mechanical behavior using a hyper-viscoelastic constitutive model based on strain energy functions and hereditary approach had a good agreement with experimental results.

solid mechanics; constitutive model; SHPB; aviation polyurethane; high strain rate; hyper-viscoelastic

10.11883/1001-1455(2015)01-0051-06

2013-05-23；

2013-10-17

国家自然科学基金项目(11372113,11472110)；国家国际科技合作项目(2011DFA53080)； 爆炸科学与技术国家重点实验室基金项目(KFJJ14-2M)

张龙辉(1991— )，男，硕士研究生;通讯作者： 张晓晴,tcqzhang@scut.edu.cn。

O347.3 国标学科代码： 13015

A