浮空器膜材循环拉伸力学性能及弹性常数研究
发布时间:2021年12月25日 点击数:1558
飞艇等浮空器自20世纪30年代, 因飞机的发展及自身安全性问题而衰落, 70年代后, 飞艇再度活跃, 航空大国一致看好飞艇的发展前景[1].蒙皮材料作为飞艇的主体结构材料, 其性能的高低直接影响飞艇的应用效能, 如浮空高度、持续飞行时间、有效载荷和服役寿命等[2,3].国内外对飞艇蒙皮膜材力学性能研究[3,4,5,6,7,8,9,10,11], 主要针对单轴试验, 对双轴试验研究很少, 缺乏系统性的研究.Nakadate和Maekawa等[5,6,7]在飞艇膜材的抗撕裂性能、耐候性和开孔加强方式等方面进行了较详细研究.黄赛帅等[9]对飞艇膜材进行了双向拉伸测试, 得到了单一比例下的近似膜材弹性常数.高海健等[10,11]对3种浮空器囊体膜材进行了单轴拉伸、往复加载和徐变试验, 得到了膜材单向加载下的力学性能参数.
在建筑膜结构领域, 国内外对建筑膜材力学性能进行了较广泛研究, 并取得了许多科研成果[12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29].在国外, Blum等[13]应用双轴拉伸试验, 建立了线性增量应力应变模型的弹性常数测试方法、计算方法.Minami等[14]通过不同比例双轴拉伸试验, 由应变残差平方和最小二乘法求弹性常数.Bridgens等[15]指出膜材细观结构中存在经纬纱的卷曲及相互摩擦等作用, 膜材不满足正交异性互补定律, 不能用单轴拉伸试验预测膜材料的双轴向拉伸弹性常数.在中国, 卫东等[18]测试了PVDF/PES膜强度、徐变和弹性模量.易红雷等[19,20]以膜材满足麦克斯韦定理为前提, 以单轴拉伸试验为基础, 研究了弹性本构关系、强度准则和弹性模量估算方法.罗仁安等[23]基于双轴试验, 提出了广义泊松比, 研究了双轴受力下应力-应变、残余应变、滞回曲线、双模量等.李阳等[24]研究了膜材抗拉强度、撕裂强度、徐变等, 采用双轴拉伸试验和应变残差和最小法求膜材弹性常数.张营营等[28]对PTFE/GF和PVC/PES膜材的单双轴拉伸破坏、循环加载变形及力学性能进行了研究.
Uretek3216L为聚酯纤维平纹织物基布, 涂层为聚碳酸脂, 表层为耐候涂层PVF, 厚度0.40 mm, 面质量300 g/m2, 可用于中型尺度飞艇 (约50.0 m下) , 在国内外飞艇结构设计中得到广泛应用.本文对Uretek3216L薄膜进行一系列单双轴循环拉力试验, 得到了滞回环曲线、双轴应力应变响应面和单双轴弹性常数等材料力学性能, 为复杂飞艇构型的设计、结构分析和裁剪提供参考.
1 循环试验概况
1.1 试验仪器及加载制度
单轴循环采用ZWICK/R011/Z100试验机, 按正弦波循环加载15次, 加载速度为10 mm/min.试验环境:温度为 (20±2) ℃, 湿度为65±2%.最大拉力为1/4单轴拉伸强度, 最小拉力取1/40单轴拉伸强度.试验膜材经纬向单轴拉伸强度分别为50.74和35.89 kN/m.
双轴循环试验采用自主研制的双轴拉伸试验机, 见图1.试验机应变测量范围0%~20%;夹具标准拉伸速率2~4 mm/min;实时性控制1~5 μs.采用精密伺服液压油缸作为动力装置、通过比例阀和溢流阀等实现流量精确控制, 采用力传感器闭环反馈, 通过PID控制器进行实时控制, 中心引伸位移计采集应变, 非常速拉伸, 从而实现任意载荷谱精确跟踪.
依据膜材单向拉伸应力应变关系及文献[25], 确定加载曲线峰值14.0 kN/m (7.0 kN/m) , 谷值2.5 kN/m.双轴拉伸加载比例 (纬∶经) 包括1∶1, 1∶2, 2∶1, 1∶0和0∶1, 以纬经向应力比1∶2为例, 加载谱见图2.
1.2 试件尺寸
单轴循环试验经纬向各5个试件, 试件总长300 mm, 宽50 mm, 两端夹具夹持50 mm, 有效长度200 mm.
双轴循环试验采用十字形试件如图3所示, 膜材伸臂纵横向交叉, 每个伸臂长宽16.0 cm, 夹持范围为4.0 cm, 伸臂开3道缝.
2 应力应变关系
2.1 单轴循环
图4为试件循环拉伸试验应力应变曲线.由图4可知, 经纬向初次加载均有较大残余应变, 经向约0.6%, 纬向约1.0%, 且在此载荷阶段膜材具有非线性.经向应力应变曲线具有3阶段特征:初始变形平缓, 然后应力变化率加大, 10 kN/m后趋缓;而纬向为2段曲线特征, 第3段平缓段不明显, 这与纬向加载值偏小有关.
选取加载上升段, 去除残余应变后应力应变曲线见图5.
由图5可知, 随循环次数增加, 曲线斜率逐渐变大, 第1次循环曲线与第2次曲线差别最大.以经向为例, 当第1次循环时, 曲线非线性特征明显, 具有小三段特征;而第2次循环, 在8.0 kN/m之前曲线几乎线性;后续循环之间差异性逐渐减小, 第15次应力10.0 kN/m之前几乎线性.纬向特征与经向基本相似, 随次数增加同应力水平下应变逐渐变小, 其第1和第2次循环的曲线差异比经向更显著.
2.2 双轴循环
以1∶2 (纬∶经) 为例, 应力应变曲线见图6, 各试件曲线相对偏移经向1.0%, 纬向0.5%.
由图6可知, 初次循环后残余应变较大, 随着循环次数的增加, 残余应变增量减小, 滞回曲线形状趋于稳定;卸载时, 应变不是直接减小, 而是先增大后减小.去除残余应变后应力应变曲线见图7.
由图7可知, 第1次循环应力应变表现为显著的非线性, 在较高应力下表现出斜率变小, 应变快速增大的非线性特征.第2次和第3次线性特征增强, 应力应变基本线性变化, 第2次与第3次循环应力应变的差异明显小于第1次与第2次循环之间的差异.这也证明了选取3次循环加载制度的可行性.第1次与第2次循环之间的差异, 纬向比经向更明显, 原因是经向纤维有预拉紧, 双向力学特征受循环的影响小于盘绕状的纬向纤维.
取每个循环周期的试验数据, 进行多项式插值拟合可得应力应变响应面, 以第3次循环为例 (图8) .由图8可知, 经纬向应力对应变影响较复杂, 呈三维非线性关系, 响应面有显著曲率变化域和相对平缓域, 当膜材双向受力在显著变化区域以下的平滑区域内, 膜材处于较有利双轴受力态.
响应面可近似为三维曲面, 采用拟合公式:
ε=aσx+bσy+c.
式中:ε为应变, %;σx和σy分别为纬向和经向应力, kN/m;a, b和c为参数, 见表1.
表1 公式参数 导出到EXCEL
Tab.1 Parameters of formula
循环次数 |
应变 方向 |
参数 |
||
a |
b | c | ||
1 |
经向 | -0.045 8 | 0.245 3 | -0.715 2 |
|
纬向 | 0.325 7 | -0.009 9 | -0.827 7 |
2 |
经向 | -0.085 7 | 0.364 7 | -0.802 0 |
|
纬向 | 0.439 4 | -0.071 2 | -0.347 9 |
3 |
经向 | -0.198 0 | 0.511 3 | -0.963 5 |
|
纬向 | 0.592 0 | -0.213 0 | -0.144 4 |
3 弹性常数
3.1 单轴循环
随循环次数增大, 膜材弹性模量逐渐增大, 采用最小二乘法拟合弹性模量, 结果见表2.
表2 循环加载弹性模量 导出到EXCEL
Tab.2 Elastic modulus with cycle test
循环 次数 |
经向 弹性模量 / (kN·m-1) |
幅度f / % |
纬向 弹性模量 / (kN·m-1) |
幅度f / % |
1 |
504.25 | 0.00 | 286.45 | 0.00 |
2 |
584.37 | 60.98 | 397.24 | 60.36 |
3 |
592.61 | 6.27 | 410.75 | 7.36 |
4 |
601.35 | 6.65 | 421.47 | 5.84 |
5 |
602.15 | 0.61 | 430.63 | 4.99 |
6 |
607.38 | 3.98 | 436.66 | 3.29 |
7 |
612.98 | 4.26 | 441.68 | 2.73 |
8 |
615.39 | 1.83 | 447.79 | 3.33 |
9 |
618.89 | 2.66 | 452.35 | 2.48 |
10 |
622.07 | 2.42 | 456.16 | 2.08 |
11 |
624.75 | 2.04 | 459.69 | 1.92 |
12 |
626.76 | 1.53 | 462.63 | 1.60 |
13 |
629.27 | 1.91 | 466.04 | 1.86 |
14 |
633.24 | 3.02 | 468.59 | 1.39 |
15 |
635.64 | 1.83 | 470.01 | 0.77 |
注:f=[Et (i+1) -Et (i) ]/[Et (15) -Et (1) ], Et (i) 为第i次循环弹性模量, t为膜材厚度.
随循环次数增加, 经纬向弹性模量均增大, 第15次和初始相比, 经向增大131.39 kN/m, 纬向增大183.56 kN/m, 分别占初始值的20.7%和39.1%.可见循环加载对纬向影响更显著.这与经纬纤维编织工艺有关, 纬向纤维以波浪卷曲方式绕经向纤维编织;纬向纤维受力被拉直其弹性模量明显大于初始卷曲状态, 而经向纤维有预应力拉紧, 变化会稍小.相邻循环间增大幅度f中第2次循环最大, 经向为60.98%, 纬向为60.36%.后续相邻循环间增大的幅度变化经纬向也存在差异, 纬向是逐次减小, 而经向整体趋势减小, 个别循环时出现起伏, 起伏随循环次数增加逐渐变弱, 具体见图9.该起伏特征与经向纤维在膜材中的伸展状态有关, 反映了纤维纱线在循环工况下力学特征存在阶段性.
弹性模量增大幅度可以拟合为循环次数的负幂次函数关系, 经纬向可统一公式为:
f=9.8 (n-1.9) -0.75. (1)
式中:n≥2为循环次数.
3.2 双轴循环
首先推导耦合弹性模量与纬经向应力比的关系
公式.纬经向应力应变关系方程为:
[σxσy]=[E11E21E12E22][εxεy][σxσy]=[E11E21E12E22][εxεy]
. (2)
式中:σx和σy分别为纬向和经向应力;εx和εy分别为纬向和经向应变.
E′xt=νE11,E′yt=νE22. (3)Ex′t=νE11,Ey′t=νE22.(3)
式中:ν=1-νxνy, νx=E21E22νx=E21E22, νy=E12E11νy=E12E11;t为膜材厚度;E′xt和E′yt分别为纬向和经向 (弹性主向) 弹性模量.
本膜材经验算知不满足麦克斯韦定理, 即
νxE′ytνyE′xt=k≠1. (4)νxEy′tνyEx′t=k≠1.(4)
双轴拉伸纬经向应力比σx/σy=Rx/Ry =R (Rx和Ry分别为纬向和经向应力约简值, 不同时为0) .
由式 (2) ~ (4) 可得:
σx=E′xtεx1−νx/Rk,Rx,Ry>0.σy=E′ytεy1−Rkνy. (5)σx=Ex′tεx1-νx/Rk,Rx,Ry>0.σy=Ey′tεy1-Rkνy.(5)
当Rx=0或Ry=0时, 双轴向拉伸试验转化为单轴拉伸, 其应力应变关系为:
σy=E′ytεy=-E′yt/νyεx, Rx=0,
σx=E′xtεx=-E′xt/νxεy, Ry=0. (6)
对式 (5) 和式 (6) 进行求导, 得双轴耦合弹性模量:
Ext=E′xt1−νx/Rk,Eyt=E′yt1−Rkνy,Rx,Ry>0. (7)Ext=Ex′t1-νx/Rk,Eyt=Ey′t1-Rkνy,Rx,Ry>0.(7)
Eyt=E′yt=−E′ytνy,Rx=0,Ext=E′xt=E′xνx,Ry=0. (8)Eyt=Ey′t=-Ey′tνy,Rx=0,Ext=Ex′t=Ex′νx,Ry=0.(8)
双轴耦合弹性模量试验分析结果见表3.
表3 膜材弹性常数 导出到EXCEL
Tab.3 Elastic constants of fabric
比例 |
循环次数 |
耦合弹性模量 / (kN·m-1) |
泊松比 | ||
Ext |
Eyt | νx | νy | ||
1∶2 |
1 | 344.82 | 505.33 | 0.31 | 0.37 |
|
2 | 533.72 | 596.27 | 0.28 | 0.34 |
|
3 | 569.81 | 651.65 | 0.22 | 0.28 |
|
均值 | 482.78 | 584.42 | 0.27 | 0.33 |
1∶1 |
1 | 302.44 | 519.76 | 0.32 | 0.42 |
|
2 | 568.62 | 618.39 | 0.21 | 0.42 |
|
3 | 613.03 | 665.76 | 0.19 | 0.29 |
|
均值 | 494.69 | 601.30 | 0.24 | 0.38 |
2∶1 |
1 | 303.58 | 884.01 | 0.17 | 0.66 |
|
2 | 527.24 | 1 163.4 | 0.27 | 0.29 |
|
3 | 567.93 | 1 212.3 | 0.21 | 0.25 |
|
均值 | 466.25 | 1 086.52 | 0.22 | 0.40 |
均值 |
1 | 316.947 | 636.32 | 0.27 | 0.48 |
|
2 | 543.193 | 792.61 | 0.25 | 0.35 |
|
3 | 583.587 | 843.26 | 0.21 | 0.27 |
由表3可知, 随循环次数增加, 耦合弹性模量增大, 第1次和第2次之间差距明显大于第2次和第3次之间.第1次和第2次弹性模量平均值相差41.6% (纬) 和19.7% (经) , 第2次和第3次相差6.9% (纬) 和5.9% (经) .双轴耦合弹性模量明显大于单轴循环值, 单轴循环纬、经向值分别和双轴纬向第2次、经向第1次循环均值相当.该规律可由公式 (7) 解释, 由1-νx/Rk<1, 1-Rkνy<1知耦合弹性模量大于单轴弹性模量;此外, 由公式 (7) 知, 随Rk数值增大, 经向弹性模量Eyt将增大, 纬向弹性模量Ext将减小.由表3数据, 经向耦合弹性模量Eyt, 纬向Ext (R≥1) 符合上述规律;在R<1时纬向Ext没有表现此特征.推断是由纬向的非线性、非弹性等特性所致, 正交各向异性的平面应力假设在解释膜材某些力学行为时存在不足.
4 结 论
本文对聚酯纤维平纹织物基布浮空器膜材进行了单双轴循环拉伸力学性能研究, 得出如下结论:
1) Uretek3216L膜材的非线性及正交各向异性的力学性能在应力应变关系、弹性常数中得到体现.
2) 单轴循环时, 随循环次数增加, 经纬向弹性模量均增大, 第15次和初始相比, 经向和纬向分别增大初始值的20.7%和39.1%.增大幅度第2次循环最大, 经纬向均占总增量的60%以上.后续循环纬向增大幅度逐次减小, 而经向个别循环时出现起伏.这与经向纤维在膜材中的伸展状态有关, 反映了纤维纱在循环工况下力学特征存在阶段性.拟合得到了弹性模量增大幅度和循环次数的负幂次函数关系公式.
3) 双轴拉伸试验, 随循环次数增加, 耦合弹性模量增大, 第1次和第2次之间差距明显大于第2次和第3次.推导了双轴耦合弹性模量和纬经向应力比R关系式;公式及试验结果表明耦合弹性模量大于单轴弹性模量;随R增大, 经向耦合弹性模量增大, 纬向耦合模量 (在R≥1时) 减小.本文得到了不同比例加载下的应力应变响应曲面, 并得到了曲面拟合公式.所得结果可为浮空器膜结构的设计分析提供参考.