第32卷第4期 2011年4月 宇 航 学 报 Journal of Astronautics Vo1.32 No4 .April 2011 ~~~一一一~一一~一一~~一撇一一 一~~c雪 一 一一 ~ 一 ~脚一 一平流层验证飞艇结构体系比较研究 高海健,陈务军,付功义 (上海交通大学空间结构研究中心,上海200030) 摘要:为研究平流层飞艇的受力、变形及合理的结构体系布局,根据平流层飞艇结构特点,提出了十种结构 体系布局方案。基于浮力与重力作用下平衡形态的原理,建立平流层飞艇结构分析方法和各体系分析模型。以一 m一 喊~Ⅺd mc~ s 1 Ⅱ甜 .一 s 川 r e. 一 {.m~=. M 曼 个25m验证飞艇作为对象进行了系统的数值分析,与工程弹性理论进行比较,论证了各体系的结构特点,最后提出 合理的结构体系。本文对平流层飞艇结构设计具有参考价值。 关键词:平流层飞艇;结构体系;平衡形态 d 。 k 一 中图分类号:V274 文献标识码:A 文章编号:1000—1328(2011)04-0713-08 DoI:10.3873/j.issn.1000一l328.2011.04.002 血 Comparison Investigation for Architecture of Stratospheric Demonstration Airship 一 一 . ~一 一 GAO Hai—jian,CHEN Wu—jun,Fu Gong—yi (Space Structures Research Centre,Shanghai Jiaotong University,Shanghai 200030,ChiⅪ na) m ∞ 0 引 言 大型飞艇的受力及适用性;文献[7—9]研究了飞艇 ~一~~一一 囊体裁剪设计、张力估算和控制动力学模型等;文献 平流层飞艇重量轻、体积大,一般以软式飞艇为 [10]利用CFD软件对浮空器进行了热数值模拟,并 主要型式。就飞艇的结构而言,传统的硬式飞艇以 初步进行了实验研究;文献[1 1]基于工程弹性理论 刚性构架作为承力体系,结构体系简单明确¨ 。软 对软式飞艇应力、变形进行了深入分析,给出了多个 式飞艇具有大变形几何非线性特征,悬挂体系复杂, 结构设计参数;文献[12]采用非线性动态流体~结 采用工程弹性理论分析具有较大误差,且不能评价 构交错积分耦合法,对某平流层飞艇的突风干扰响 形态变化 。目前,对软式飞艇结构的研究还不 应特性进行了分析。 多。文献[4]研究了大型飞艇非线性气弹分析方 飞艇基于“轻于空气(LTA)”设计原理,浮重平 法;文献[5]基于平衡形态原理分析了245m软式飞 衡、浮重分布及有效传递是决定飞艇结构体系及承 艇形态;文献[6]基于工程弹性理论分析了200m级 载能力的核心,浮重对低速飞艇结构体系的制约大 收稿日期:2010 ̄1—15; 修回日期:2011432—23 基金项目:国家自然科学基金(50808122、50878128) 714 宇航学报 第32卷 于气动力,浮重与推力比一般大于50。因此,本文 半径,r ( )、r ( )分别为回转曲线半径函数一、二 基于浮重“平衡形态”原理,根据平流层飞艇的结构 阶导数 为囊体环向、纵向张力,t为囊体膜材厚 特点,提出了十种结构体系布局及分析模型,通过系 度,r( )为囊体截面半径,P为飞艇内外压差(无梯 统的计算分析论证结构体系,并与工程弹性理论比 度变化均匀压差), 。( )、or ( )分别为囊体环向、 较,验证其有效性,同时深入分析各体系的受力与变 纵向应力。 形特征。 在飞艇中间部分,囊体接近圆柱面,则式(1)、 1工程弹性理论 (2)可简化为: =P・r (5) 飞艇结构初步设计时常应用工程弹性理论¨ , 1 =÷p・r (6) 飞艇囊体假设为理想回转体,在内压力作用下,囊体 产生双向张力,环向和纵向张力分别为: 2结构体系设计方案 = -( )。t=P‘r ( )(1一 ro ,)(1) 平流层飞艇结构系统主要结构组件包括外气 =or:( )・ =÷p・ ( ) (2) 囊、副气囊、头锥、尾锥、尾翼、悬挂屏、悬挂索、分舱 膜、吊舱等,其中悬挂系统的布置方式是决定结构整 I I (3) 个体系受力的主要因素。根据飞艇大小、用途,吊舱 可在主囊体内、外吊、悬吊,悬挂系统可分为悬挂屏 r ( )=r( )√l+r ( ) (4) 结合悬索,以及采用分舱膜等,因此,可将结构体系 式中: ( )、 ( )分别为囊体纵向和环向的曲率 综合概括为十种型式,如表1。 表1平流层飞艇结构体系模型分类 Table 1 Classification of the structure system model for stratospheric airship 第4期 高海健等:平流层验证飞艇结构体系比较研究 715 3计算分析方法 束反力为零(尽量小),在吊舱的一定位置布置可移 动荷载,此荷载同时可以改变飞艇的重心,可作为飞 飞艇在浮力、重力和竖向荷载作用下的平衡形 艇姿态控制的一种手段。线约束为模拟底部龙骨整 态是决定飞艇结构体系的关键,而基于平衡形态建 体刚度及约束,四点约束模拟吊舱与囊体的点连接 立结构有限元分析模型的难点在于约束的处理。本 形式。 文考虑浮力和自重,以囊体底腹部约束点竖向反力 作为净浮力的平衡力,保证其余约束点无 、l,方向 4结构体系分析 反力,使模型的平衡形态与实际静态悬浮飞艇一致。 以三段线型25m验证飞艇为模型,采用ANSYS 对于提出的十种结构体系模型简化为三类约束条 进行分析,外气囊膜单元采用三角形shell41单元, 件,一类模型约束条件为:浮心正下方节点 、l,、z 悬吊索采用linklO单元,考虑几何非线性与应力刚 一瞬燃麟峨口口四- 三方向约束;正上方节点 、z两方向约束;头锥、尾 化效应。膜材参数取:厚0.482mm,弹性模量 锥处z方向约束。二类模型约束条件为:浮心正下方 706.39MPa,密度597.82Kg/m ,泊松比0.38;索参 主吊舱长度方向 、y、z三方向约束。第三类模型约 数取:直径2.5mm,弹性模量1.65×10 MPa,密度 束条件为:浮心正下方吊舱与囊体四个连接点 、】,、 7.8×10 Kg/m ,泊松比0.3。 z三方向约束。三类模型所承受载荷有:浮力、外气 4.1气压影响分析 囊重力作为面荷载,头锥、尾锥、副气囊、尾翼的重力 首先对同一模型在不同气压下的囊体应力 作为集中力加载。其中浮力是以内外压差的梯度作 ( )和位移(d)变化进行分析。对外气囊的顶线 用于外气囊上,压差为(p 一P )h,P。为空气密度,P (11)和浮心位置处飞艇横截面囊体周线(r1)上应 浮升气体氦气密度,h为距飞艇最低点高度。对于第 力和位移进行了比较。s表示距起始点(一类模型正 一类模型,若正上方节点存在 方向约束反力,则模 下方约束点)弧长。 型与实际受力条件不一致,为了使该点处 方向约 仃_圈礅嘲蕊∞口 _ 帆 ”””"磊 躬 ”拍 图1模型F不同气压下的应力云图 Fig.1 The stress of model F with various pressure 7l6 宇航学报 第32卷 P--400 Pa 蛊器 高一 Ⅱ Mmil2 2 2 3 向0站鲒"”l∞.蛊嗣髫吕_ 3 5 8攒嚣 l 3 ”越5 9 2 扑6 9二 搿 P 600Pa P=tO00 Pa 图2模型F不同气压下的位移云图 Fig.2 The displacement of model F with various pressure 表2各气压下模型F—rl处竖向间距与横向间距 Table 2 Vertical and horizon distance at rl of model F 。500O _一黼_一圈口团一 警 with various pressure 爱 lm m "n" '=0 S 81il 盯 nn拈 幅”"辣n 2 2 Z :.0皇衄宣圈霉 40O0 3000 m ∞ 曲乱札n眙 越 2000 目1000 0 l000 图1和图2分别是模型F在不同气压下的应力 —2000 和位移云图。计算结果表明:模型F外气囊的应力 .3000 40O0 -随气压的增大而增大,最大位移随气压的增大反而 减小。其原因是气压增大使外气囊膜材所受张力增 大,从而使应力增大;而气压增大又导致飞艇整体的 刚度变大,反而减小了外气囊的最大位移。最大应力 位置在约束处,最大位移在囊体顶部。图3可以看到 囊体在r1处的位移变化。由于浮力实际产生于外气 -40O0.3000.2000.1000 0 1000 2000 3000 4OOO z/mm 图3模型F—r1在不同气压下的位移示意图 Fig.3 The displacement at rl of model F with various pressure 囊内外压差沿y轴方向的梯度变化,在梯度压差作 图5一图8为不同气压下B、c、D、E四种体系 布局的飞艇囊体在r1处的位移,结果表明:B、C、E 用下,艇身的圆形截面在y轴方向被拉长(如表2), 呈“梨型”变化趋势。由式(1)、(2)可得rl处 = : =三种结构体系布局r1处囊体位移随着气压的增大 而减小,因为气压较小时,压力梯度的作用效果比较 0.5073。而数值计算主应力比如图4所示。 第4期 高海健等:平流层验证飞艇结构体系比较研究 717 0.7 0.6 0.5 越0.4 0_3 0.2 0.1 0.O /arm 图4模型F—rl处解析值与数值解的主应力比 Fig.4 The principal tension ratio at rl of model F 35O 3O0 25O 2 200 飞150 l0O 50 O 0 5000 l0000 15000 20000 25000 /mm 图5模型B—rl在不同气压下的位移 Fig.5 Displacement with various pressure at model B—r1 25。 2。。 重 。 飞1。。 5。 0 0 5000 l0000 l5000 20000 25000 /mm 图6模型C—rl在不同气压下的位移 Fig.6 Displacement with various pressure at model C—rl 明显,从而使囊体位移较大。随着气压增大,囊体刚 度逐渐增大,囊体在压力梯度等荷载作用下的变形 反而减小;D模型由于横截面内有悬挂索张拉增强 了横截面方向内的刚度,使飞艇囊体的横向膨胀不 如前述三种模型容易,导致囊体在r1处的位移先随 气压增大稍许增大,然后再减小。由此可见,结构体 400 350 300 g 250 20o 15O l0O 5O 0 0 5000 l0000 15000 20000 25000 s/mm 图7模型D—rl在不同气压下的位移 Fig.7 Displacement with various pressure at model D—rl 22O 2OO l8O l60 l40 目120 100 80 60 4O 2O 0 0 5000 l0000 l5000 20000 25000 /mm 图8模型E—rl在不同气压下的位移 Fig.8 Displacement with various pressure at model E—rl 系布局不同对飞艇囊体的形态影响较大。 4.2平衡约束影响分析 根据结构体系布局不同,假设不同的平衡约束条 件使计算模型更合理,并进行了在相同气压下飞艇不 同结构体系布局囊体的应力和位移的比较分析。 50 4O 一3O 芝 b 20 l0 O 0 5000 10000 15000 20000 25000 /mm 图9 400Pa气压作用下模型A—Ell应力变化 Fig.9 Stress change at 11 of model A—E under 400Pa 718 宇航学报 第32卷 H表3 400Pa气压作用下各模型rl处竖向与横向间距及r1和11处应力 Table 3 Distances at rl and stresses at rl/11 for the model under 400Pa 400 400 300 量 量 00 200 200 100 l0O 0 0 0 5000 10o00 l5000 200O0 25000 O 5000 100oo l5000 20000 25000 z/mm /nun 图10 400Pa气压作用下模型A—Ell位移变化 图12 400Pa气压作用下模型A—Erl位移变化 Fig.10 Displacement change at 11 of model A—E 11 Fig.12 Displacement change at rl of model A—E under 400Pa under 400Pa 1O 8 日6 室 b 4 2 0 0 5000 l0oOO l5000 20000 25000 0 5000 10000 15000 20000 25000 /nun /nnn 图13 400Pa气压作用下模型F—Jll应力变化 图11 400Pa气压作用下模型A—Erl应力变化 Fig.1 3 Stress change at 11 of model F—J under 400Pa Fig.1 1 Stress change at rl of model A—E under 400Pa 囊体应力和位移都比较小,且比较均匀,而r1处的 图9一图16为400Pa气压作用下分别在I1、r1 位移较小但应力相对较大,且在悬挂索与囊体连接 处的各模型囊体应力和位移。结果表明,对于第一 处有应力集中。对于第二类约束条件的四种结构体 类约束条件的五种结构体系布局,模型E在l1处的 第4期 高海健等:平流层验证飞艇结构体系比较研究 719 35O 300 250 目200 目 150 1O0 50 0 0 5000 10000 15000 20000 25o00 /mm 图14 400Pa气压作用下模型F—Jll位移变化 Fig.14 Displacement change at 11 of model F—J under 400Pa 4.0 3.5 3.0 b 2.5 2.O 0 5000 10000 15000 20000 25000 /ram 图15 400Pa气压作用下模型F—Jrl应力变化 Fig.1 5 Stress change at rl of model F~J under 400Pa 350 300 250 吕200 目 150 lO0 50 O 0 5000 10000 l5000 20000 25000 /mm 图16 400Pa气压作用下模型F—Jr1位移变化 Fig.1 6 Displacement change at rl of model F—J under 400Pa 系布局,模型H在l1和rl处的囊体应力和位移都 是最小的。第三类约束条件模型J的受力和变形与 模型H比较接近。 对各模型r1处竖向与横向间距及rl和l1处应 力进行分析(如表三),可以看到,模型B、C、E、H的 竖向间距变形较小,分别是1.8%、0.9%、2.3%、 1.7%,模型D、G、H、I的横向间距变形较小,分别是 0.2%、一0.2%、~0.5%、一0.1%,模型B、C、H、I 在r1处应力较小,分布较均匀,其应力均方差分别 为0.91、0.69、0.18、0.26。同样,考察l1处应力,模 型A、E、F、H较合理。所以,通过综合比较可以得 到模型H是较好的一种结构体系,可实现较好的结 构作用,同时满足高空平台特殊副气囊布局要求。 5 结 论 本文提出了十种平流层飞艇结构体系布局,并 基于平衡形态原理建立了分析模型,以一个25m验 证飞艇作为对象进行了系统的数值计算和分析,且 与工程弹性理论计算结果进行了比较。得到以下结 论:(1)不同的结构布局对平流层飞艇应力和位移 的影响较大,其囊体变形对不同气压的响应也不同, 悬挂索的合理布设能改变囊体的受力形态,改变飞 艇的刚度。(2)由于平流层飞艇的实际结构很复 杂,与简单的回旋体差异较大,所以工程弹性理论计 算结果与数值结果相差较大,只适用于飞艇总体评 估时作参考。(3)合理的体系设计能充分发挥材料 的受力性能,使囊体受力、变形较小,且比较均匀。参 考 文 献 []一l 一 720 craft Research,2006,1(1):41—46.] 宇航学报 第32卷 lation and experiment for an airship under low altitude environ— [6]Chen W J,Xiao W W,Kroplin B H,et al Structural perform- ance evaluation procedures for large flexible airship of HALE ment[J].Jourual of Astronautics,2010,31(10):2417— 2421.] stratospheric platform conception[J].Journa1.of Shanghai J'iao- tong University,2007,E一12(2):293—298. [11] 高海健,陈务军,付功义,等.平流层平台柔性飞艇结构弹性 分析理论[J].上海交通大学学报,2010,44(11):1583— 1588.[Gao Hai—jian,Chen Wu-jun,Fu Gong—yi,et a1.Strut— tural analysis oriented engineering elastic theory for flexible air- [7] 王文隽,李勇,姚伟,等.飞艇气囊压力与蒙皮张力的估算 [J].宇航学报,2007,28(5):1109—1112.[Wang Wen. jun,Li Yong,Yao Wei,et a1.Estimation of the relationship be— ship of stratospheric platform[J].Journal of Shanghai Jiaotong tween the pressure in airship ballonet and the tension in its envel- op[J].Journal of Astronautics,2007,28(5):1109—1112.] [8] 陈务军,肖薇薇,任小强,等.大型柔性飞艇外囊体裁剪设 计分析方法[J].浮空器研究,2007,1(37):35—41.[Chen Wu-jun,Xiao Wei-wei,Ren Xiao-qiang,et a1.Cutting design procedure for large non—rigid airship envelope[J].Floating Air- craft Research,2007,1(37):35—41.] [9] 高明伟,单雪雄.改变飞艇重心位置控制纵向运动[J].力学 季刊,2006,27(4):714—718.[Gao Ming.wei,Shan Xue. xiong.By changing position research of airship’s center of gravity to control longitudinal motion[J].Chinese Quarterly of Mechan- ies,2006,27(4):714—718.] [10] 程雪涛,徐向华,梁新刚.低空环境中浮空器的热数值模拟 与实验研究[J].宇航学报,2010,31(10):2417—2421. [Cheng Xue—tao,Xu Xiang—hua,Liang Xin—gang.Thermal simu一 University,2010,44(11):1583—1588.] [12] 王晓亮,单雪雄,陈丽.平流层飞艇流固耦合分析方法研究 [J].宇航学报,2011,32(1):22—29.[Wang Xiao—liang, Shah Xue-xiong,Chen Li.Study on fluid-sturcture coupled tom- putational method for stratosphere aisrhip[J].Journal of Astro— nautics,2011,32(1):22—28.] [13] 铁木辛柯S,奥诺斯基S.板壳理论[M].北京:科学出版 社,1977. 作者简介:高海健(1976一),男,博士研究生,从事空间结构 研究。 通信地址:上海市华山路1954号上海交通大学徐汇校区浩 然高科2304室(200030) 电话:(021)62932191 E—mail:gaohj@sjtu.edu.cn (编辑:沃云峰)