• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

YB-2航空有机玻璃的应变率和温度敏感性及其本构模型

史飞飞 索涛 侯兵 李玉龙

史飞飞, 索涛, 侯兵, 李玉龙. YB-2航空有机玻璃的应变率和温度敏感性及其本构模型[J]. 爆炸与冲击, 2015, 35(6): 769-776. doi: 10.11883/1001-1455(2015)06-0769-08
引用本文: 史飞飞, 索涛, 侯兵, 李玉龙. YB-2航空有机玻璃的应变率和温度敏感性及其本构模型[J]. 爆炸与冲击, 2015, 35(6): 769-776. doi: 10.11883/1001-1455(2015)06-0769-08
Shi Fei-fei, Suo Tao, Hou Bing, Li Yu-long. Strain rate and temperature sensitivity and constitutive model of YB-2 of aeronautical acrylic polymer[J]. Explosion And Shock Waves, 2015, 35(6): 769-776. doi: 10.11883/1001-1455(2015)06-0769-08
Citation: Shi Fei-fei, Suo Tao, Hou Bing, Li Yu-long. Strain rate and temperature sensitivity and constitutive model of YB-2 of aeronautical acrylic polymer[J]. Explosion And Shock Waves, 2015, 35(6): 769-776. doi: 10.11883/1001-1455(2015)06-0769-08

YB-2航空有机玻璃的应变率和温度敏感性及其本构模型

doi: 10.11883/1001-1455(2015)06-0769-08
基金项目: 国家自然科学基金项目(11372256,11272267,11228206,11202168);高等学校学科创新引智计划项目(B07050)
详细信息
    作者简介:

    史飞飞(1981—), 女, 博士研究生

    通讯作者:

    索涛, suotao@nwpu.edu.cn

  • 中图分类号: O347.1

Strain rate and temperature sensitivity and constitutive model of YB-2 of aeronautical acrylic polymer

  • 摘要: 为了理解和评价YB-2航空有机玻璃在极端环境下的动态力学性能,采用电子万能试验机和分离式Hopkinson压杆对YB-2航空有机玻璃在218~373 K温度范围、10-3~3 000 s-1应变率范围内的压缩力学行为进行了研究,得到了材料的应力应变曲线。结果表明:随着温度的升高,材料的流动应力逐渐减小而破坏应变呈现增大的趋势;温度相同时,材料的流动应力随应变率的增加而增大,破坏应变随应变率的增加而减小。随着应变率的提高,材料的应变软化效应更加剧烈。基于朱-王-唐(ZWT)本构模型,得到了考虑温度效应的本构参数。结果显示,在8%应变范围内,改进的考虑温度效应的本构模型可以较为理想地表征该材料的应力应变响应。
  • 飞机座舱风挡、座舱盖和弦舱材料为高聚物材料, 在飞机飞行时需要承受一定的温度和载荷作用。在起飞和着陆时, 还有可能会遇到飞鸟的撞击, 军用飞机可能在空战中被非爆弹丸或高爆碎片击中, 这相当于对这些位置的材料施加一定的高速率载荷。同时, 由于飞机在低、高空飞行时, 环境温度相差很大, 高聚物材料的力学性能受外界温度的影响很大。因此, 开展对高聚物材料在高温和高应变率下的力学性能研究具有重要意义。

    C.G’Sell等[1]对聚合物在常应变率下的变形行为进行了研究, 发现应变率影响材料的应变硬化现象。陈建桥等[2]、吴衡毅等[3]通过恒应变率的实验对有机玻璃进行了研究, 结果表明, 在较低应变率下, 有机玻璃的力学性能具有明显的应变率相关性。E.M.Arruda等[4]、O.S.Lee等[5]分别研究了准静态及中、高应变率下温度和应变率对PMMA材料力学性能的影响, 发现随着温度的升高, 材料的应力应变行为发生了显著变化, 同时伴随有热软化现象发生。G.W.Adams等[6]、Machida等[7]通过研究温度和应变率对PC材料力学性能的影响后发现:温度在293~338K之间时, 材料变形过程中储存的能量对实验温度的依赖性并不明显。S.M.Walley等[8]在一定的应变率范围内对一些高聚物材料进行了力学性能测试, 解释了在高应变率下由于绝热效应的原因, 使得材料出现应变软化行为。索涛等[9-10]对几种航空有机玻璃进行了研究, 认为在某一应变率和某一温度下, 由于试样内部的微裂纹使得有机玻璃表现出明显的应变软化行为; 在相同温度下, 随着应变率的升高, 有机玻璃的流动应力增大, 而破坏应变减小, 并建立了考虑温度效应的本构模型。

    本文中对YB-2航空有机玻璃在不同温度下的准静态和动态力学行为进行系统的研究。分别采用Hopkinson分离式压杆和CSS44100型电子万能试验机对材料进行初始温度分别为218、253、293、343、373 K, 应变率分别为0.001、0.1、700、3 000 s-1下的压缩实验。同时, 结合实验数据, 基于朱-王-唐模型, 拟合考虑应变率效应以及温度效应的本构模型, 并给出模型预测和实验所得的应力应变曲线。

    有机玻璃是一种高分子材料, 对温度等环境比较敏感。由于聚合体材料的导热系数很低, 约是金属的1%, 即聚合体的导热系数为0.2 mm/min, 试样在加热或冷却的过程中, 如果加温或冷却的时间不充分, 试样表面温度和心部温度就会不一致, 导致实验结果不准确。因此, 在涉及温度的实验时, 为了确保实验中试样受热均匀和热平衡, 从而保证实验数据的准确性, 首先测试了YB-2有机玻璃在不同温度下试样表面和试样内部达到相同温度时所需要的时间。实验中, 采用正方体试样, 边长为10 mm。在试样中心沿着厚度方向钻孔, 孔的深度为5 mm。通过孔洞, 在试样心部引入一对热电偶。将另外一对热电偶缠绕在试样表面。用密闭的电热炉, 对试样进行加热或在密闭的容器中喷洒低温氮气对试样进行冷却。当经过标定的两组热电偶显示的温度相同时, 认为试样内部和外部的温度达到了平衡。

    低应变率压缩实验在CSS44100型电子万能试验机上进行, 压缩试样为直径6 mm、长4 mm的圆柱体, 如图 1所示。两种低应变率分别为0.1、0.001 s-1, 温度分别为218、253、293、343、373 K。实验前将压头降温或加热至相对应的实验温度后对试验机进行空压, 得到位移随着压力变化关系, 以此标定试验机刚度, 用于消除试样位移测量中试验机位移的影响。低温通过液氮循环获得, 高温实验采用电阻丝加热法对试样进行加热, 并用一个与试样相连的热电偶测量和显示温度, 通过电磁继电器的通断控制温度, 使之保持在实验值的±5 ℃误差范围内。

    图  1  试样尺寸图和试样的不同取向
    Figure  1.  Sample dimensions and its cutting directions

    700、3 000 s-1时的高应变率压缩实验是在一套直径12.7 mm的分离式Hopkinson压杆上进行的, 压杆所用材料为马氏体时效钢(18Ni(C250)), 实验装置见图 2。采取同低应变率下相同的试样。利用撞击杆撞击入射杆产生的应力波对放置于入射杆和透射杆之间的试样进行加载; 通过粘贴在两根杆子上与惠斯通电桥相连接的应变片采集输出的电压信号; 利用超动态应变仪和高频数据采集系统记录到实测波形, 并根据一维应力波原理可以得到材料的应力应变关系。高低应变下低温和高温的获得采用相同的方法。

    图  2  分离式Hopkinson压杆示意图
    Figure  2.  Schematic drawing of split Hopkinson pressure bar

    图 3(a)和(b)为加热或者冷却温度与保温时间的曲线以及不同的实验温度所需最短时间值。可以看出图 3(a)中, 当经过一定时间后, 无论是高温加热还是低温冷却, 试样内部和外部的温度达到一致, 记录加热或冷却的时间, 就得到了该温度下试样内外温度达到平衡所需的时间。从图 3(b)可以看出, 低温冷却试样内外温度平衡所需时间比高温时所需时间要短, 随着温度的逐渐升高, 平衡所需时间呈现逐渐变短的趋势。当温度高于有机玻璃的玻璃化转变温度后, 平衡时间变长。但如果温度继续升高, 平衡所需的时间几乎相同。

    图  3  不同温度下保温时间曲线和保温所需最小时间
    Figure  3.  Time-temperature curves and minimum times at different temperatures

    YB-2航空有机玻璃有别于定向有机玻璃, 是未经过取向的高聚物, 分子链段的排列是随机的, 因此理论上材料体现出各向同性。在开始实验之前, 为了测试该材料的各向性能, 首先对其进行了考虑方向性的准静态压缩实验。试样分别取自于垂直于板材和平行于板材两个方向, 见图 1。实验结果(见图 4)表明:沿2个方向取样所测出的应力应变曲线基本吻合, 这说明材料力学各向同性。在后续的实验中采用的试样沿V方向取样。

    图  4  不同取向试样的真实应力应变曲线
    Figure  4.  True strain-stress curves in two directions

    图 5给出了在不同温度下试样的真实应力应变曲线。可以看出, 高分子材料的力学性能强烈依赖于温度和应变率。图 5(a)中, 当应变率为0.001 s-1时, 随着温度的升高, 材料的弹性模量呈现降低的趋势, 屈服强度减小, 塑性变形能力增加。在其余的3种应变率下(见图 5(b)~(d)), 材料呈现相似的趋势。在室温及以下的高应变率时, 材料的破坏应变随着温度的升高而增大, 此时, 材料为玻璃体呈现出典型的脆性破坏形式, 这是因为在低温高应变率下, 材料处于玻璃态, 且加载时间非常短暂, 材料的分子结构没有足够的时间发生弛豫, 因而断裂特性表现为脆性; 而当实验温度在材料的玻璃化转变温度(90 ℃)附近时, 材料处于高弹态, 其承载能力急剧下降, 变形能力增强, 不再呈现脆性破坏, 这是由于分子结构的活化能降低、自由体积增大, 因此受力后分子结构中的链段运动更加容易[11]

    图  5  不同温度下的真实应力应变曲线
    Figure  5.  True strain-stress curves at different temperatures

    图 6给出了同一温度下不同应变率时的应力应变曲线。在218 K时(见图 6(a)), 应变率对材料的弹性模量和屈服应力的影响并不明显, 但屈服应力仍呈现随应变率升高而升高的趋势, 破坏应变则反之。当温度升高时(见图 6(b)~(d)), 应变率对材料的弹性模量的影响逐渐变大, 屈服应力以及破坏应变呈现同218K时相似的趋势。4种应变率下, 材料均出现不同程度的应变软化。低应变率时, 由于加载脉冲宽度有限, 材料在实验中并未发生破坏。已有研究结果[12]表明:低温时, 聚合物处于玻璃态, 分子链的链段被冻结发生松弛的时间是无穷大, 此时, 加载速率(即应变率)的改变对聚合物材料变形的影响作用非常微弱, 表现为弹性模量对应变率的不敏感; 而在高温时, 聚合物处于高弹态, 其变形实际上是一种强迫高弹形变的过程, 分子链的链段已经被激活, 并且分子链的链段运动受到松弛时间的影响非常大。由于应变率的变化对松弛时间的影响很大, 因此在高温时弹性模量对应变率较为敏感。

    图  6  不同应变率下的真实应力应变曲线
    Figure  6.  True strain-stress curves at different strain rates

    图 7为不同应变率时温度-应力峰值曲线。当实验温度为218 K时, 不同应变率下试样的应力峰值差值较小。随着温度的升高, 应力峰值的差值变大, 在室温时应力峰值的差值达到最大。随后, 该差值又逐渐变小。这表明:在室温时材料的应变率敏感性达到最强。在高应变率加载下, 应变率对材料应力峰值的影响逐渐弱化。

    图  7  温度-应力峰值曲线
    Figure  7.  Temperature-peak stress curves

    朱兆祥等[13]对典型工程塑料进行的一系列研究发现:在变形小于8%, 应变率在10-4~103 s-1范围内, 朱-王-唐模型(ZWT)可以较好地描述高聚物材料的本构模型。该模型描述如下:

    $$ \sigma(\varepsilon, \dot{\varepsilon})=E_{0} \varepsilon+\alpha \varepsilon^{2}+\beta \varepsilon^{2}+E_{1} \int_{0}^{t} \dot{\varepsilon} \exp \left(-\frac{t-\tau}{\theta_{1}}\right) \mathrm{d} t+E_{2} \int_{0}^{t} \dot{\varepsilon} \exp \left(-\frac{t-\tau}{\theta_{2}}\right) \mathrm{d} t $$ (1)

    式中:E0αβ为材料的弹性常数。2个积分项分别用来描述低应变率和高应变率下的黏弹性响应, E1θ1E2θ2分别是所对应的Maxwell单元的弹性常数和松弛时间。

    为描述温度对材料的影响, 索涛等[10]对朱-王-唐模型进行修正, 并分别对4种航空有机玻璃实验所得的真实应力应变曲线进行拟合, 实验与模型预测的曲线吻合较好。修正的模型中加入的温度效应项, 其形式为/f0(ε, T0), 它随温度的升高而减小, 随应变的增大而增大。改进后的模型描述为:

    $$ \sigma(\varepsilon, \dot{\varepsilon}, T)=\sigma(\varepsilon, \dot{\varepsilon}) f(\varepsilon, T) / f\left(\varepsilon, T_{0}\right) $$ (2)
    $$ f(\varepsilon, T)=a(\varepsilon) \exp [b(\varepsilon) T] $$ (3)

    式中: )为式(1)对应的温度和应变率下得到的应力值, T0为选定的参考温度。

    根据公式(1)~(3)可得到温度影响项的表达式: , 由实验数据拟合以多项式给出。

    基于该模型, 以293 K作为参考温度, 针对YB-2试样, 经过拟合实验曲线得到常温时的6个参数:E0=3.26 GPa、α=-20.12 GPa、E1=2.4 GPa、θ1=0.205 s、E2=2.622 GPa、θ2=6 μs。

    在此基础上, 采用最小二乘法拟合可以得到YB-2航空有机玻璃温度效应影响项中b(ε)的估计表达式为:

    $$ b(\varepsilon ) = 0.6429{\varepsilon ^2} - 0.1336\varepsilon - 0.004$$ (4)

    对于拟合参数优度的评估, 我们将对回归参数的方差和置信区间进行分析。根据实验结果, 对公式(4)中各系数的方差和置信区间的分析如下:假定置信系数为0.95, 计算得出拟合函数的方差和置信区间见表 1

    表  1  温度影响项中各系数的方差和置信区间
    Table  1.  Variances and confidence intervals of the coefficients for thermo-item
    参数方差置信区间
    a0.381 131(0.527 425, 1.273 1)
    b0.009 826(-0.235 72, -0.103 88)
    c0.000 295(-0.004 56, -0.003 54)
    下载: 导出CSV 
    | 显示表格

    图 8给出了拟合结果与实验结果的比较。不难看出, 图 8(a)~(c)中拟合曲线同实验曲线吻合较好; 对于图 8(d), 拟合曲线和实验曲线有较大差异, 这是因为在PMMA的热变形温度(70 ℃)左右时, 材料将开始发生玻璃态与高弹态之间的相态转变。而ZWT模型主要用来描述材料的黏弹性行为, 因此在拟合343 K下材料的变形行为时表现出较大的误差。总体来说, 改进后的ZWT模型可以较好的描述一定温度范围内材料的非线性黏弹性力学行为。

    图  8  修正的ZWT模型拟合结果与实验结果的比较
    Figure  8.  Comparison of prediction using modified ZWT model with experimental results

    通过对YB-2航空有机玻璃的研究得出以下结论:

    (1) YB-2有机玻璃的力学性能强烈依赖于应变率和温度:在低应变率下, 随着温度的升高, 材料的弹性模量和流动应力逐渐减小, 弹性变形能力降低, 塑性变形能力增加; 在高应变率下, 材料的流动应力随着温度的升高逐渐减小, 破坏应变随着温度的升高而增大。在室温及以下时, 材料为玻璃体, 呈现出典型的脆性破坏, 弹性模量对应变率不敏感; 当实验温度在材料的玻璃化转变温度附近时, 高应变率下材料破坏不再呈现脆性破坏, 弹性模量对应变率较为敏感。

    (2) 随着应变率的升高, 材料的流动应力显著升高, 但破坏应变减小。在室温时, 材料的应变率敏感性达到最强。在高应变率下, 材料的屈服强度受应变率的影响逐渐弱化。

    (3) 加入温度效应的ZWT模型能够较好地拟合应变范围在8%以内的实验曲线, 但由于模型中未考虑损伤以及高分子高温下自由体积和复杂的分子链运动的影响, 因而不能较准确地描述应变软化以及大应变时的应变硬化现象。

  • 图  1  试样尺寸图和试样的不同取向

    Figure  1.  Sample dimensions and its cutting directions

    图  2  分离式Hopkinson压杆示意图

    Figure  2.  Schematic drawing of split Hopkinson pressure bar

    图  3  不同温度下保温时间曲线和保温所需最小时间

    Figure  3.  Time-temperature curves and minimum times at different temperatures

    图  4  不同取向试样的真实应力应变曲线

    Figure  4.  True strain-stress curves in two directions

    图  5  不同温度下的真实应力应变曲线

    Figure  5.  True strain-stress curves at different temperatures

    图  6  不同应变率下的真实应力应变曲线

    Figure  6.  True strain-stress curves at different strain rates

    图  7  温度-应力峰值曲线

    Figure  7.  Temperature-peak stress curves

    图  8  修正的ZWT模型拟合结果与实验结果的比较

    Figure  8.  Comparison of prediction using modified ZWT model with experimental results

    表  1  温度影响项中各系数的方差和置信区间

    Table  1.   Variances and confidence intervals of the coefficients for thermo-item

    参数方差置信区间
    a0.381 131(0.527 425, 1.273 1)
    b0.009 826(-0.235 72, -0.103 88)
    c0.000 295(-0.004 56, -0.003 54)
    下载: 导出CSV
  • [1] G'Sell C, Jonas J J. Determination of the plastic behavior of solid polymer at constant true strain rate[J]. Journal of Materials Science, 1979, 14(3): 583-591. doi: 10.1007/BF00772717
    [2] 陈建桥, 李铁萍, 李之达, 等.有机玻璃本构关系研究[J].机械科学与技术, 2006, 25(3): 371-374.

    Chen Jian-qiao, Li Tie-ping, Li Zhi-da, et al. Experimental research on the constitutive relationship of PMMA[J]. Mechanical Science and Technology, 2006, 25(3): 371-374.
    [3] 吴衡毅, 马钢, 夏源明. PMMA低、中应变率单向拉伸力学性能的实验研究[J].实验力学, 2005, 20(2): 193-197. http://www.oalib.com/paper/1413732

    Wu Heng-yi, Ma Gang, Xia Yuan-ming. Experimental study on mechanical properties of PMMA under unidirectional tensile at low and intermediate strain rates[J]. Journal of Experimental Mechanics, 2005, 20(2): 193-197. http://www.oalib.com/paper/1413732
    [4] Arruda E M, Boyce M C, Jayachandran R. Effects of strain rate, temperature and thermomechanical coupling on the finite strain deformation of glassy polymers[J]. Mechanics of Materials, 1995, 19: 193-212. doi: 10.1016/0167-6636(94)00034-E
    [5] Lee O S, Kim M S. Dynamic material property characterization by using split Hopkinson pressure bar(SHPB)technique[J]. Nuclear Engineering and Design, 2003, 226(2): 119-125. doi: 10.1016/S0029-5493(03)00189-4
    [6] Adams G W, Farris R J. Latent energy of deformation of bisphenolapolycarbonate[J]. Journal of Polymer Science, 1988, 26(2): 433-445. doi: 10.1002/polb.1988.090260216/pdf
    [7] Machida T, Lee D. Deep drawing of polypropylene sheets under different heating conditions[J]. Polymer Engeering & Science, 1988, 28(7): 405-412. doi: 10.1002/pen.760280702
    [8] Walley S M, Field J E, Pope P H, et al. A study of the rapid deformation behavior of a range of polymers[J]. Philosphical Transactions of the Royal Society of London: A, 1989, 328: 783-811.
    [9] 索涛, 李玉龙, 刘元镛.温度、应变率对航空PMMA压缩力学性能的影响研究[J].材料科学与工程学报, 2006, 24(4): 547-550. http://www.cnki.com.cn/Article/CJFDTotal-CLKX200604018.htm

    Suo Tao, Li Yu-long, Liu Yuan-yong. Study on temperature and strain rate effects on mechanical behavior of aeronautical PMMA[J]. Journal of Materials Science & Engineering, 2006, 24(4): 547-550. http://www.cnki.com.cn/Article/CJFDTotal-CLKX200604018.htm
    [10] Suo T, Li Y L, Yu H, et al. Temperature effect on the mechanical behavior of acrylic polymers under quasi-static and dynamic loading[J]. Jounal of Shanghai Jiaotong University, 2004, 38: 89-95. http://d.wanfangdata.com.cn/Periodical/shjtdxxb2004z2017
    [11] 焦剑, 雷渭媛.高聚物结构、性能与测试[M].北京: 化学工业出版社, 2003: 438-444.
    [12] 何曼君, 张红东, 陈维孝, 等.高分子物理[M].上海: 复旦大学出版社, 2007: 105-198.
    [13] 朱兆祥, 徐大本, 王礼立.环氧树脂在高应变率下的热粘弹性本构方程和时温度等效性[J].宁波大学学报:理工版, 1988, 1(1): 58-68. http://www.cnki.com.cn/Article/CJFDTotal-LBJY198801009.htm

    Zhu Zhao-xiang, Xu Da-ben, Wang Li-li. Thermoviscoelastic constitutive equation and time-temperature equivalence of epoxy resin at high strain rates[J]. Journal of Ningbo University: Natural Science & Engineering Edition, 1988, 1(1): 58-68. http://www.cnki.com.cn/Article/CJFDTotal-LBJY198801009.htm
  • 期刊类型引用(1)

    1. 张世文,龙建华,贾宏志,刘仓理. 柱壳约束对散心冲击波空间分布的影响. 爆炸与冲击. 2018(02): 345-352 . 本站查看

    其他类型引用(3)

  • 加载中
图(8) / 表(1)
计量
  • 文章访问数:  4670
  • HTML全文浏览量:  707
  • PDF下载量:  745
  • 被引次数: 4
出版历程
  • 收稿日期:  2014-11-11
  • 修回日期:  2015-02-04
  • 刊出日期:  2015-12-10

目录

/

返回文章
返回