平流层飞艇凭借静升力驻空, 因其具有高空侦察、通信中继和空间探索等诸多领域的百富策略白菜网优势, 成为美国、欧盟、日本等世界主要军事大国研究的热点[1,2,3,4,5].轻质高强的层压类织物膜材广泛百富策略白菜网在飞艇的主 (副) 气囊、尾翼等关键部件中, 其力学响应及弹性参数是飞艇蒙皮结构计算分析、工程设计、确定制造工艺的基础[2,3,4].飞艇蒙皮材料为层压类织物膜材, 和建筑用涂层织物类膜材相比, 在成型工艺、功能层结构及纤维类型方面存在一定差异[5], 但二者在纱线编织及力学响应分析方面存在相通之处.织物膜材 (包括涂层类及层压类) 因经纬纱线交叉编织及多层复合的结构, 其力学性能复杂, 具有明显的材料非线性、几何非线性、非弹性及正交异性等特征, 其弹性参数的确定也因此具有不确定性和复杂性, 一直是国内外学者的研究热点[6,7,8,9,10,11,12,13,14].

飞艇蒙皮结构的膜面基本处于平面受力状态, 双向加载下的弹性参数对结构设计分析非常关键.限于试验条件, 在国内外的工程百富策略白菜网中对膜材平面应力下变形响应特征及弹性参数的研究及百富策略白菜网仍有待深入[6,7,8,9,10].

目前, 对膜材双轴拉伸弹性参数的计算, 基本以线弹性假设为条件, 求取有限的弹性参数[8,9,10,15].而双轴拉伸受力下膜材的力学性质, 是一项多因素、多角度的内容, 其弹性常数宜涵盖膜材多方位的力学性质[16].另外, 蒙皮材料独特的纱线编织方式及多功能层的复合结构, 使其弹性参数具有多变性和复杂性, 简单少量的弹性参数信息很难全面描述膜材的力学性质, 势必降低工程设计时的精确度.因此, 进一步剖析膜材在双轴受力下弹性常数响应特征及变化规律, 探讨膜材力学响应的力学本质, 对准确把握蒙皮材料在实际飞艇结构中的力学性质及指导结构设计分析具有重要意义, 而国内外对此研究尚不足.本文针对高性能蒙皮材料Uretek3216LV?开展一系列多比例双轴拉伸试验, 深入研究蒙皮材料弹性参数在应力空间上的响应特征及分布规律, 探讨分析正交异性及循环加载作用对弹性参数的影响.提出基于应力空间的弹性参数加权计算方法, 获得了材料弹性参数随应力水平及应力比例的变化规律, 并且结合数值模型分析验证了所得参数变化规律的正确性和适用性.

试验材料为热致液晶 (LCP) 芳族聚酯类膜材Uretek3216LV?, 主要由Vectran?基布和聚氟乙烯 (PVF) 面层及各功能膜层压合而成.Uretek3216LV?膜材是新型高性能蒙皮材料的代表, 具有高强、高比强、高模量、耐氧化、性态稳定、耐磨等优点, 尤因其具有轻质、气密性优、耐强辐射、耐高温的特点, 而在国内外大中型飞艇、航天领域得到广泛使用[17].

蒙皮材料Uretek3216LV厚度为0.21mm, 面密度为200g/m2, 其织物组成结构及材料外观如图1所示.织物为平纹组织结构, 织物经纬向密度为17×12根/cm, 经纬纱线细度均为200Denier.膜材的粘贴层为弹性性能良好的丙烯酸化合物.

双轴拉伸试件尺寸如图2所示, 采用十字形切缝试样, 可实现多应力比荷载的施加, 是目前普遍采用的形式[10,18].试样按膜材的经纬向对称取样, 核心区域为16.0cm×16.0cm, 悬臂长16.0cm, 夹具夹持范围为4.0cm.为有效传递应力使中心区域应力分布较均匀, 悬臂间隔约4.0cm预制切缝, Bridgens[19]和Chen[20]已证实均匀切缝的存在可有效扩大应力均匀区面积.裁剪直角将引入较大的应力集中[5], 本试件采用半径15.0 mm的过渡圆弧以减弱应力集中影响.

试验环境参考GB/T 6529—2008, 试验室相对湿度 (65±4.0) %, 温度 (20±2) ℃.

双轴循环试验采用自主研制的双轴拉伸试验机SJTU-I, 如图3所示.双轴拉伸试验机应变测量范围:-10%~20%;夹具标准拉伸速率:2~4mm/min;实时性控制:1~5ms.采用精密伺服液压油缸作为动力装置, 通过比例阀、溢流阀等实现流量精确控制, 采用力传感器闭环反馈和PID控制器进行实时控制, 可实现任意载荷谱的精确跟踪.采用2个Green Pot LP-20F位移引伸计测量位移并计算应变, 其量程为20mm.

采用在等预应力水平基础上的多比例加载路径, 依据膜材单向拉伸应力应变关系并参考文献[21], 设置加载谷值2.5kN/m, 峰值因比例而异, 比如:1∶1峰值均为12.7kN/m, 1∶2峰值为12.7kN/m和6.25kN/m, 具体如图4所示.每个比例3个循环, 循环周期T=10 min, 各比例间隔ΔT=4min;在每个比例前, 先做3个1∶1循环, 以消除前一比例影响.由于飞艇蒙皮膜材与传统建筑膜结构有所差异, 其应力比相对集中, 基本上在1∶2~2∶1间, 为凸显这个区域内的应力比, 本文选取9种比例:0∶1, 1∶3, 1∶2, 2∶3, 1∶1, 3∶2, 2∶1, 3∶1和1∶0, 以获取更详细、更真实的材料力学响应.

参考标准 (MSAJ) [21], 采用应变残差法计算弹性参数.经计算知膜材Uretek3216LV不满足正交互补定理, 结果见表1.Ex和Ey分别为经纬向弹性模量, νx和νy分别为相应泊松比参数.

据表1, 3循环下, 膜材Uretek3216LV的K系数均大于1.0, 基本上在1.2~1.7.作者曾对聚酯类膜材Uretek3216L进行了分析, 结果和Uretek3216LV膜材基本一致.2类膜材的K系数数值和Gosling和Bridgens[22,23]对聚酯类PVC膜材及玻璃纤维类PTFE膜材的研究结果基本相似.这表明由于受到材料及几何非线性、非弹性特征的影响, 织物膜材本构关系是否满足正交互补关系值得商榷.

应变响应曲面可全方位呈现膜材应力应变响应特征, 是求取弹性参数响应曲面的基础.本文基于应变残差平方和最小的拟合原理 (式 (1) ) , 采用MAT-LAB编程对各应力比试验数据进行三维曲面拟合.所得应变响应曲面如图5所示.

式中:y (x) 为连续函数f (x) 在区间[a, b]内的逼近函数.

据图5, 响应曲面存在明显的起伏和卷曲, 这是由膜材力学参数的改变引起的.曲面的起伏和倾斜是膜材的模量、泊松比等参数改变的外在表现.响应曲面的起伏多变, 反映出弹性常数在应力空间的分布及膜材自身力学性质上的复杂性.基于上述响应曲面, 在应力空间上一点 (σx0, σy0) 的邻域对称选N个数据点 (σxi, σyi) (i, j为序号) , 百富策略白菜网应变残差最小二乘法, 求得 (σx0, σy0) 点处的弹性常数代表值Eij (i, j=1~4) .分区完毕后依次沿双向遍历所有应力空间点, 进而拟合可得膜材弹性参数的响应曲面.

弹性模量在应力空间上的分布如图6所示.

据图6可知, 弹性模量在所选应力空间内, 分布并不均匀, 存在明显的起伏特征.在曲面上形成典型的峰域和缓坡, 陡缓之间的梯度差异很大, 在峰域的等值线分布密集, 缓坡区域等值线分布稀疏.依据弹性模量高低将空间划分为多个区域, 如图7所示.

由图7可知, 从趋势来看, 经纬向弹性模量均随所在向应力的增加而升高.从所处应力空间位置来看, 经纬峰域更接近于区域B, 该区域经向应力弱于纬向应力;这表明经纬向刚度最大值并非出现在相同应力水平的区域, 这是因为经纬向纱线在编织几何参数、卷曲程度、预应力水平等方面存在差异, 纬向纱线一般卷曲度高, 可认为刚度稍“弱”于经纱.双向受力下, 经纬向高弹性模量的出现, 应首先克服纬纱的“弱”点, 使纬纱受力稍高于经纱, 二者的弹性模量才会同步处于高水平.低弹性模量区 (区C) , 基本位于经纬向应力的低水平区域, 显然在经纬向并不对称, 1∶1对称轴右侧 (纬向应力较大侧) 的低弹性模量区较大, 经纬向弹性模量高低的分布特征根源于膜材内在构成的差异性.

经纬向弹性模量在变化趋势、分布特征上也存在相似之处, 表现出一定的同步特征.由图8可知, 1st循环的弹性模量图, 经纬向弹性模量从 (0, 0) 点出发, 沿1∶1 (近纬向侧) 方向爬升, 当然同时经纬向沿各自应力增大向爬升, 形成了相似的峰域和缓坡区域.因经向弹性模量大于纬向弹性模量, 二者形成上下交叉叠合的形状.

和弹性模量相比, 经纬向泊松比的差异性较显著 (如图8所示) .比如变化幅度, 以第1循环为例, νs变化范围为0.17~0.75, νy范围稍小为0.20~0.50.此处大泊松比 (>0.5) 的出现反映了膜材经纬向相互作用及变形机理的复杂性.

当然, 2个泊松比在形状方面也存在一些相似之处:二者均从 (0, 0) 出发近似沿应力比1∶1方向下降, 形成“V”形, 如图8 (c) 所示.在近 (0, 0) 应力点处较大, 应力比1∶1方向最大应力点 (σx-max, σy-max) 处几乎最小.但整体上泊松比规律性不明显, 变化趋势不规律.

循环荷载的作用使弹性模量明显提高, 第3循环时平缓区域的弹性模量, 经向在1 500kN/m左右, 纬向在1 200kN/m左右, 显著大于第1循环的平缓区域数值:经向1 100kN/m左右, 纬向900kN/m左右.峰域的分布也大有变化, 第1循环时的单峰, 变成了第3循环的周边多峰特征 (如图8 (b) 所示) .第3循环时, 经纬向在3个区域出现峰域, 高弹性模量出现在经纬向应力均较大或仅一向大的区域.泊松比数值有所变化, 虽然个别区域的数值较大, 但大部分区域的数值变得更均匀, 中间区域的等值线变得稀疏, 泊松比随循环增加有趋稳特征.

基于应力空间内的分析结果, 提出弹性常数加权计算法, 对弹性常数响应曲面在应力空间内积分, 积分的加权均值作为相应范围的弹性常数.弹性常数Emn (m, n=1, 2) 的表达式为:

式中:Eij (σx, σy) 为应力场 (σx, σy) 处的弹性常数;Di为第i组应力场 (共N组) ;WDi≤1为应力场Di的弹性常数权函数, 对应Di应力场面积占总面积的比重 (见式 (3) ) , 其和满足∑WDi=1.在工程结构分析时, 权函数可依据应力场在膜面中的分布比例来定.比如, 飞艇的主气囊结构, 膜面应力比基本处于1∶2~2∶1, 在此范围内权值建议取0.6~0.8.

依据所提方法, 将应力场进行分区处理, 依次选取Ⅰ~Ⅳ及全域Ⅴ (分区如图9 (a) 所示) 进行计算, 结果见表2.此处, 所用权函数均为1.

由表2, Ⅴ (图9 (a) ) 对应整个应力场的弹性参数数值, 是整个试验加载范围内膜材的整体参数响应, 从数值大小看, 和前面章节中计算结果 (见表1) 略有差异, 但基本相当.但在沿着1∶1线的4组分区中, 弹性常数差异显著.从趋势看, Ⅰ~Ⅳ4个分区, 应力水平逐渐增加, 膜材弹性模量也依次增加.可见, 应力水平的改变, 可对膜材的弹性参数产生明显影响, 且应力水平增加可引起膜材刚度的增大.

依据应力比分区 (如图9 (b) 所示) , 弹性常数计算结果如图10所示.据图10, 无论第1还是第3循环, 经纬弹性模量均表现出“X”型交叉特征.随名义应力的增加, 经向弹性模量以增大趋势为主, 纬向弹性模量以减小趋势为主.飞艇主气囊结构, 应力比以1∶2~2∶1为主, 因此建议优先选取该比例范围的

弹性模量作为分析参数.在该应力比例范围时经纬弹性模量第1循环分别为:1 279和987kN/m, 第3循环分别为:1 627和1 430kN/m, 计算结果和表1相比, 纬向稍大、经向稍小.当应力比扩大到1∶3~3∶1时, 弹性模量有所改变, 2个应力比范围的结果均列于图10的椭圆圈内.与应力比范围1∶2~2∶1相比, 1∶3~3∶1范围内弹性模量约存在5%的变化.

泊松比随名义应力比的变化结果如图11所示, 基本上呈现出中间凹两端高的特征.经纬向应力相差较大时, 泊松比出现大于0.5的数值, 表明经纬向间存在显著的相互作用, 同时说明蒙皮材料有别于均质材料, 泊松比变化范围较大.

应力空间加权均值法, 可实现对不同应力水平和应力比范围内的双轴拉伸弹性参数的计算和分析;对具体结构或部件, 依据其应力分布特征, 针对性地确定相应应力水平和应力比范围的弹性参数, 可提高结构计算分析的精确度.对于飞艇结构, 可将主 (副) 气囊、尾翼、吊屏等结构部件分别处理, 依据各自的应力分布特征, 确定相应的弹性参数, 并可按应力比或应力水平所占比重, 设定合理权值, 进一步提高所计算弹性参数的精确性和适用性.

基于3.4节弹性参数变化规律, 进行膜材双轴拉伸非线性分析, 验证所提出的结论.模型基本假设如下:

1) 材料性质满足正交异性的平面应力假设, 在特定应力比下, 应力应变关系满足线弹性:

2) 经纬向弹性模量Ex和Ey随名义应力比γy发生非线性变化, 存在如下关系:

3) 泊松比νxy随应力比的变化趋势不明显, 认为泊松比数值不变, 以减小模型复杂性.

模型含有7个参数:ax, bx, ay和by为系数参数;Ex (0.1~1.3) 和Ey (0.1~1.3) 为经纬向0∶1~1∶3比例范围内, 采用应力空间加权均值法计算出的弹性模量;泊松比νx∈ (0, 0.5) .上述模型参数可采用模量残差平方和最小二乘法求取, 模量残差平方和最小值Θ表达式:

式中:m为应力比个数;Exijmod () 和Exitest () 分别为模型和试验弹性模量结果, 括号内为其参数;ax, bx, ay, by为系数参数;νx为泊松比.依据式 (6) , 编制相应程序求取参数, 结果见表3.

运用有限元软件Abaqus的UMAT二次开发技术, 将上述本构模型百富策略白菜网于有限元数值分析.依据材料试样尺寸, 建立了相应的数值模型, 该数值模型采用四边形壳单元S4R, 模型单元数目4 816.荷载独立地加载于伸肢端部, 以模拟所有应力比工况, 图12所示为1∶1和1∶3加载时纬向应力应变关系.

图12中, NL模型为本文非线性模型, S模型的弹性参数为一组固定值.据图12, 在应力比1∶1时, NL材料模型和S材料模型差异较小, 均和试验曲线接近.但是在两肢应力差异较大时, NL比S模型具有更好的精确度.这表明随名义应力比改变的本构关系式在膜材应力应变预测上具有适用性, 同时证明了应力空间内弹性参数分析结论的科学性.

基于蒙皮材料的多比例双轴拉伸试验, 获得了材料在平面应力空间内的应变及弹性参数响应特征.主要结论如下:

1) 芳纶聚酯类膜材正交互补系数均大于1.0, 基本上在1.2~1.7范围内.表明, 正交互补性质百富策略白菜网于该膜材的弹性参数计算分析时适用性不足.

2) 在应力空间上, 经纬向弹性模量变化趋势的相似度很大, 经纬向均沿各自应力增大向爬升;主要差异是弹性模量的数值和峰值域的位置.而随名义应力比的增加, 经向弹性模量以增大趋势为主, 纬向弹性模量以减小趋势为主, 经纬弹性模量呈现出“X”型交叉特征.

3) 经纬向泊松比差异性较显著, 2个泊松比在应力空间上, 形成类似“V”形形状.泊松比存在大于0.5的数值, 表明蒙皮材料异于均质材料, 反映了其经纬向相互作用及变形机理的复杂性.

4) 应力空间加权均值法, 可实现对不同应力水平和应力比范围弹性参数的计算和分析;可针对具体结构或部件, 确定相应应力水平和应力比范围的弹性参数, 提高结构设计分析的精确度.