文章快速检索  
  高级检索
基于副气囊耦合的飞艇纵向运动影响分析
游颖捷1,2 , 糜攀攀2 , 吕明云2 , 刘文定1     
1. 北京林业大学 工学院, 北京 100083;
2. 北京航空航天大学 航空科学与工程学院, 北京 100083
摘要: 高空大型飞艇副气囊体积占比很大,甚至达90%以上,副气囊动力学对飞艇运动产生较大影响。对巡航高度为6 000 m的某大型飞艇,建立基于圆柱容器的副气囊动力学模型,进一步构建考虑副气囊的飞艇纵向耦合动力学方程,并对飞艇纵向运动进行了数值仿真,将其与不考虑副气囊时的结果进行对比。结果表明:副气囊对大型飞艇俯仰运动有较大影响,容易引起运动发散,且耦合影响随飞艇高度增大而变小,同时副气囊体积不变时,其直径越小,动力学耦合效应影响越小;副气囊效应会降低飞艇的操纵性,在进行飞艇设计时应该予以考虑。
关键词: 飞艇设计     副气囊     纵向运动     飞行动力学     平流层飞艇    
Airship's longitudinal motion impact analysis based on coupling ballonet
YOU Yingjie1,2 , MI Panpan2 , LYU Mingyun2 , LIU Wending1     
1. School of Engineering, Beijing Forestry University, Beijing 100083, China ;
2. School of Aeronautic Science and Engineering, Beijing University of Aeronautics and Astronautics, Beijing 100083, China
Received: 2015-06-25; Accepted: 2015-09-18; Published online: 2016-06-25 12: 00
Foundation item: Industry-research Project of China Aviation Industry(Cxy2014BH04)
Corresponding author. Tel 010-82338116 E-mail:lv503@buaa.edu.cn.
Abstract: Large high-altitude airship has a high volume ratio which even exceeds 90%, so it will induce significant impact on airship's motion because of the ballonets dynamics. Considering one large-scale airship whose cruise altitude is 6 000 m, ballonets dynamics based on the cylindrical container was established and airship longitudinal motion equations containing ballonets were further formulated. Finally, numerical simulation of airship longitudinal motion was performed and results were compared with that of airship model without ballonet dynamics. Results indicate that ballonets dynamics will make a difference to airship's pitch motion, easily inducing motion divergence and ballonets' coupling impacts reduce when the airship rises. In the meantime, the smaller the diameter of ballonets is, the smaller the dynamic coupling effect is when ballonets' volume is unchanged. In addition, the controllability of airship will be weakened because of ballonets, which should be taken into account during airship design.
Key words: airship design     ballonet     longitudinal motion     flight dynamics     stratospheric airship    

中高空大型飞艇(飞行高度为6~0 km)及平流层飞艇是一种主要靠浮力升空飞行的飞行器,由于其驻空时间长、载重量大和覆盖范围广等特点,能实现高分辨率对地观测及通信中继等,在军事和民用等方面有着广泛的应用前景,近年来受到各国广泛关注,成为各国航空领域的研究热点[1]

一个飞艇通常有2个或3个副气囊。低空飞艇体积较小,副气囊体积约为艇体体积的20%~30%,而20 km左右的平流层飞艇体积巨大,副气囊体积占比会超过90%[2]。所以,大型飞艇的俯仰运动会引起艇体内副气囊的运动,反过来,副气囊运动也会影响飞艇的俯仰运动,这种耦合动力学会进而影响飞艇的纵向稳定性和操纵性,而且一般副气囊体积越大,影响越明显。NASA相关研究报告[3]对液体在运动容器中的动力学模型进行了研究,并分析了其在航空航天领域的应用,奠定了理论基础; Maekawa和Saito[4]将副气囊简化为圆柱容器,分析了其对海平面上的某小型飞艇纵向运动的影响; Delaurier[5]通过建立弹簧阻尼系统模拟了系留飞艇副气囊的运动,分析了其对飞艇运动的影响。Nakadate[6]和Li等[7]对SPF-2飞艇进行了飞行试验,发现飞艇升空过程中副气囊体积变化很大,其质心纵向位置也发生较大变化,证明艇体与副气囊之间存在较强的动力学耦合效应。

本文针对某大型飞艇,建立其纵向小扰动运动方程,在此基础上,借鉴Maekawa和Saito[4]圆柱容器的副气囊模型,建立飞艇与副气囊动力学耦合方程,同时考虑在不同高度下飞艇的质量变化,分析了副气囊运动对飞艇纵向稳定性与操纵性的影响,并分析了副气囊不同设计尺寸对飞艇运动的影响,可为飞艇设计及副气囊参数确定提供参考。

1 飞艇运动建模 1.1 飞艇模型结构

本文研究的大型飞艇为一个基于YEZ-2A[8]飞艇改进设计的双椭圆常规布局的软式飞艇,尾翼采用“+”结构布局,吊舱位于艇囊下方,飞艇质心位于飞艇纵剖面内体积中心正下方,矢量推力螺旋桨在吊舱两侧对称布置,推力倾角为μ,如图 1所示。该飞艇巡航高度为6~8 km,部分参数如表 1所示。

图 1 大型飞艇示意图 Fig. 1 Large airship sketch
表 1 大型飞艇参数 Table 1 Large airship parameters
参数数值
飞艇长度/m129.5
飞艇最大直径/m32
飞艇囊体体积/ m369 413
起飞质量/t82.6
海平面副气囊体积占比/%74.6

1.2 飞艇动力学建模及其线性化

由于飞艇运动过程中高度变化和副气囊充放气的影响,飞艇重心位置发生变化,为便于描述和建模,选择飞艇艇体坐标系进行动力学建模,原点选为飞艇体积中心,同时引入以下几个假设:

1) 飞艇为刚体,忽略其弹性效应,即保持飞艇体积不变。

2) 飞艇具有纵向对称面,且重心位于纵向对称平面内,惯性积Ixy=Iyz=0。

3) 飞艇体积中心与浮心重合。

利用牛顿第二定律得到大型飞艇的六自由度动力学模型如下:

(1)

式中: u、v和w分别为艇体坐标系中轴向、侧向和垂直方向上的线速度;p、q和r分别为绕坐标轴的滚转、俯仰和偏航角速度;M为质量矩阵:

(2)

其中:mx=m+mairk1,m为飞艇质量,mair为飞艇所排开的空气质量;my=mz=m+mairk2;Jx=Ix;Jy=(1+k3)Iy;Jz=(+k3)Iz;xG和zG为飞艇重心坐标;k1、k2和k3为势流理论分析得到的椭圆惯性因子;,Ix、Iy和Iz分别为绕重心坐标Obxb、Obyb和Obzb的转动惯量,l和d分别为飞艇总长度和最大直径[9-12]

Fd为动力矢量,只与速度量有关:

(3)

FA为气动力矢量,计算采用Mueller和Paluszek[13]的半经验模型:

(4)

其中:α为迎角;β为偏滑角;δRUD为方向舵偏角;δELV为升降舵偏角;X1a、Y1a和Z1a为气动力;L1a、M1a和Na为气动力矩。

FB为浮力矢量,FG为重力矢量,FT为推力矢量,都表示为3个坐标轴的力和力矩分量,这里不再赘述。

将飞艇纵向运动分为基准运动和扰动运动,在一定的基准运动状态下,对飞艇动力学方程进行线性化,选择状态向量为xL=[u,w,q,θ]T(θ为俯仰角),控制输入向量为uL=[T u δe]T,得到飞艇纵向小扰动方程如下:

(5)
(6)
(7)
(8)

式中:dz为推力在艇体坐标系的z坐标;δe为纵向升降舵操纵量;为飞艇浮力;G为飞艇重力;下标e表示飞艇基准运动状态;Xua、Xwa、Zua和Zwa等为气动力导数,计算方法参考文献[14-17]的飞艇模型。

2 副气囊动力学模型

对于飞艇副气囊的动力学模型,其基本思想是借鉴流体在运动容器中的运动机理,将容器内的流体作为一个多体弹簧阻尼系统,同时将副气囊这个柔性结构近似为一个刚性容器[4]。事实上,由于薄膜的抗拉刚度、压力维形和张线的存在,副气囊的形变可能很小。因此,在本文中假设成一个刚性圆柱容器是合理的,考虑副气囊的飞艇如图 2所示。

CG—飞艇重心;—推力;lfx—前气囊轴向距离;lax—后气囊轴向距离;lfz—前气囊z向距离;laz—后气囊z向距离。 图 2 包含副气囊的飞艇 Fig. 2 An airship with ballonets

考虑到容器是一个直径为d的圆柱,且内部充满了高度为h流体,则该流体在运动中的自然频率为

(9)

式中:g为重力加速度;ξn为1阶第一类贝塞尔函数的解。

(10)

图 3所示,用该等效模型来表示流体的动态运动。有

c.g.—副气囊的质心位置;I0—固定块的转动惯量。 图 3 副气囊动力学模型 Fig. 3 Dynamic model of ballonet
(11)

式中:m1为容器中流体的总质量;mn为与第n个自然频率一致的质量;kn为弹簧刚度;hn为mn的位置;m0为固定块的质量;h0为固定块的高度;ρ为密度。

使用这些参数来表示各质量块的运动方程为

(12)

式中:y0为重心的位移;yn为第n个质量的附加位移;φ0为重心的转动角度;γn为第n个系统的临界阻尼比。通过每一个质量块的运动来表示对容器的反作用力:

(13)

因此,作用在容器上重心的反作用力为

(14)

式中:

3 副气囊与飞艇方程耦合

飞艇通常有一个副气囊在前面,一个在后面,如图 3所示。从第2节可以得出,无限质量系统可以当作副气囊的力学模型。然而,高于二阶的质量比第一阶系统小很多,甚至可以忽略不计,故本文计算时忽略二阶以上的部分。表 2给出了某飞艇副气囊的各阶等效参数[16],其中副气囊几何参数如下:d=3.24 m,h=4 m,ρ=1.225 kg/m3,mT=40.4 kg。可以发现,忽略二阶以上质量块是合理的。

表 2 副气囊的各阶等效参数 Table 2 Each order equivalent parameters of ballonet
n ξnmn/kgωn/(rad·s-1)
11.8417.4383.338
25.3310.2245.681
38.5360.0537.188
411.7060.0218.418
514.8640.019.486

副气囊坐标系和飞艇运动坐标系之间的关系为

(15)

式中:lx和lz分别为副气囊重心与艇体的中心在x和z方向上的距离。因此,

(16)

令y1为飞艇前部副气囊的一阶质量m1位移,y2为它的速度(y2=y1);y3为飞艇尾部副气囊的一阶质量,y4为它的速度(y4=y3)。然后把式(2)转换为以下形式:

(17)

式中:

(18)
(19)
(20)
(21)

其中:下标fa分别表示飞艇前部和尾部的副气囊。

由副气囊晃动引起的作用在飞艇上的力可以转换到艇体坐标系表示为

(22)
(23)

把式(23)、式(22)用矩阵表示:

(24)

式中:P、G、H及J定义如下:

(25)
(26)
(27)
(28)

把副气囊作用耦合到飞艇运动方程中,并定义以下状态变量:

(29)

可得加上副气囊作用的飞艇耦合运动方程为

(30)

式中:

4 数值仿真 4.1 飞艇稳定性分析

采用前面的大型飞艇模型,分别仿真飞艇在海平面(0 m)、3 000 m及巡航高度(6 000 m)的俯仰稳定性,结果见图 4。其中,假设飞艇前后副气囊形状相同,副气囊相关参数见表 3

图 4 不同高度下飞艇的稳定性特性 Fig. 4 Airship stability characteristics at different heights
表 3 飞艇副气囊参数(海平面) Table 3 Ballonet parameters of airship (sea level)
参数dh lfx>laxlfzlaz
数值4022.9525403.223.52

图 4给出了飞艇在不同高度下受到0.05 rad俯仰角扰动情况下的动态响应。

可以看出,在H=0 m处,当不考虑副气囊动力学效应时,飞艇参数经过小幅振荡后很快收敛;而考虑副气囊效应后,飞艇俯仰稳定性受到很大影响,纵向参数u、w、q和θ在50 s后均出现振荡加剧而发散的结果。这是由于海平面飞艇副气囊体积占比最大(达74.6%),副气囊与飞艇耦合运动明显,极易导致飞艇俯仰运动发散。

当飞艇在3 000 m和6 000 m时,副气囊体积逐渐减小,导致副气囊对飞艇运动的影响减小。从图 4(c)发现,考虑副气囊时,飞艇除参数u收敛较慢外,w、q和θ等均经过小幅振荡,在50 s后达到收敛。当飞艇飞至6 000 m高空时,副气囊体积占比减小至52.8%,副气囊运动耦合效应明显减小,各项参数在40 s后很快收敛,具有很好的俯仰稳定性。

4.2 飞艇操纵性分析

针对该大型飞艇模型,数值仿真飞艇在海平面(0 m)、3 000 m及巡航高度(6 000 m)的飞艇俯仰操纵响应,结果如图 5所示。其中,飞艇基准状态为速度V=11.0 m/s,迎角为.962°,升降舵输入为0.1 rad的阶跃响应。

图 5 飞艇升降舵0.1 rad响应特性 Fig. 5 Airship response of elevator at 0.1 rad

对比图 5(a)~图 5(c)可以看出,不考虑副气囊效应时,飞艇在输入响应下50 s后就很快达到稳定状态。考虑副气囊效应时,与稳定性类似,在海平面(H=0 m),副气囊体积占比较大,耦合效应明显,其运动基频为0.15 Hz,相应当量质量为13.6 t,故飞艇在0. rad阶跃输入下纵向各项参数均发生振荡,并逐渐发散;当飞艇飞行高度增加时,如到达3 000 m和6 000 m时,随着副气囊占比下降,其振荡基频分别为0.48 Hz和0.143 Hz,相应当量质量为9.9 t和6.8 t,其耦合效应减小,u、w、q和θ飞艇各项参数均能够经过小幅振荡后达到收敛,但与不考虑副气囊情况相比,二者u、w和θ的稳态值并不相等,例如,6 000 m处飞艇在副气囊耦合影响下,其俯仰角θ稳态响应值与不考虑耦合相比减小近0.03 rad(1.8),这说明副气囊效应降低了飞艇的操纵性,在进行飞艇设计时应该予以考虑。

4.3 飞艇副气囊形状影响分析

为了研究副气囊形状对其耦合效应的影响,选取该飞艇在3 000 m处基准状态下,保持其副气囊体积不变,直径增大为,高度缩小为20.24×0.5 m(即),使其形状更加“矮胖”,给定升降舵0. rad阶跃输入,得到飞艇响应如图 6所示;同时给定0.05 rad的俯仰角扰动得到飞艇的动态响应如图 7所示。

图 6 飞艇响应 Fig. 6 Airship response
图 7 飞艇的稳定性特性 Fig. 7 Airship stability characteristics

图 6图 7可以看出,副气囊体积不变的情况下,当副气囊直径扩大后,飞艇耦合动力学效应的影响有所变强,振荡周期增大。特别是飞艇轴向速度u、发散较快,其他参数也都在50 s后逐渐发散,这是因为副气囊形状改变后,其基频变为0.097 Hz(变形前为0.48 Hz),很接近飞艇纵向钟摆频率0. Hz,同时相应当量质量变为6.9 t(变形前为9.9 t),导致飞艇产生较大的动力学耦合现象,引起飞艇俯仰运动发散。

相反地,当保持仿真条件不变,让副气囊直径缩小为原来的,高度就增加为原来的2倍(即),通过仿真得到该情况下飞艇的基频为0.8 Hz,相应当量质量为3.67 t,这使得飞艇基频远离飞艇钟摆频率,同时当量质量也减小了很多,耦合效应明显减小,飞艇的输入响应几乎与不考虑副气囊效应时的情况相同,如图 8图 9所示。

图 8 飞艇响应 Fig. 8 Airship response
图 9 飞艇的稳定性特性 Fig. 9 Airship stability characteristics
5 结 论

本文对某大型飞艇,建立其纵向小扰动方程,借鉴Maekawa和Saito圆柱容器的副气囊模型,建立了基于副气囊的飞艇耦合动力学模型,着重分析了副气囊效应对大型飞艇俯仰稳定性及操纵性的影响,得到如下结论:

1) 对于体积巨大的高空大型飞艇,副气囊占比很大,其动力学耦合效应不可忽略,对飞艇俯仰运动产生很大影响,随着飞艇升空,副气囊体积减小,耦合效应也逐渐减小。

2) 大型飞艇在副气囊耦合效应下的操纵稳态响应值小于忽略耦合效应的响应值,即耦合效应将降低飞艇的操纵性,且随着飞艇高度增加,这种影响会逐渐减小。所以,在进行飞艇设计的时候需要予以考虑。

3) 保持飞艇副气囊体积不变,副气囊形状越“瘦高”,其基频越远离飞艇钟摆频率,动力学耦合效应的影响越小;形状越“矮胖”,副气囊运动基频越接近飞艇钟摆特征频率,耦合现象越明显。

参考文献
[1] CHU A,BLACKMORE M,CHOLENDT R,et al.A novel concept for stratospheric communications and surveillance:Star light:AIAA-2007-2601[R].Reston:AIAA,2007.
Click to display the text
[2] YOKOMAKU Y.Overview of stratospheric platform airship R&D program in Japan[C]//The 2nd Stratospheric Platform Systems Workshop,2000:21-22.
Click to display the text
[3] ABRAMSON H N. The dynamic behavior of liquids in moving containers[M]. 版本 Washington,D.C.: NASA, 1996 : 1 -21.
Click to display the text
[4] MAEKAWA S, SAITO K. The effect of ballonet slosh on an airship's longitudinal motion[J]. Transactions of the Japan Society for Aeronautical and Space Sciences,2004, 47 (155) : 44 –50.
Click to display the text
[5] DELAURIER J. Influence of ballonet motions on the longitudinal stability of tethered aerostats[J]. Journal of Aircraft,1980, 17 (5) : 305 –312.
Click to display the text
[6] NAKADATE M.Development and flight test of SPF-2 low altitude stationary flight test vehicle[C]//AIAA 5th ATIO and AIAA 16th Lighter-than-Air Systems Technology Conference and Balloon Systems.Reston:AIAA,2005:999-1011.
Click to display the text
[7] LI Y W, NAHON M, SHARF I. Airship dynamics modeling:A literature review[J]. Progress in Aerospace Sciences,2011, 47 (3) : 217 –239.
Click to display the text
[8] ASHRAF M Z, CHOUDHRY M A. Dynamic modeling of the airship with Matlab using geometrical aerodynamic parameters[J]. Aerospace Science and Technology,2013, 25 (1) : 56 –64.
Click to display the text
[9] ATKINSON C J,URSO R G.Modeling of apparent mass effects for the real-time simulation of a hybrid airship[C]//AIAA Modeling and Simulation Technologies Conference.Reston:AIAA,2006:821-832.
Click to display the text
[10] TUCKERMAN L B.Inertia factors of ellipsoids for use in airship design:NASA-TR-210[R].Washington,D.CNASA,1926.
Click to display the text
[11] ZHAO X W, YANG B C, HUANG S X. Calculation of the rotary inertia of ellipsoid[J]. Physics and Engineering,2007, 17 (2) : 28 –29.
Click to display the text
[12] KHOURY G A, GILLETT J D. Airship technology[M]. Cambridge: Cambridge University Press, 2012 : 34 -56.
Click to display the text
[13] MUELLER J B,PALUSZEK M A.Development of an aerodynamic model and control law design for a high altitude airship[C]//AIAA Unmanned Unlimited Technical Conference,Workshop and Exhibit.Reston:AIAA,2004,1:415-431.
Click to display the text
[14] KULCZYCKI E A,JOHNSON J R,BAYARD D S,et al.On the development of parameterized linear analytical longitudinal airship models:AIAA-2008-7260[R].Reston:AIAA,2008.
Click to display the text
[15] SEMERARO N, GRIMALDI T, PESA M L, et al. New approach and results on the stability and control of airship[J]. Pathologica,2014, 62 (911) : 344 –354.
Click to display the text
[16] SASA S, HARADA K. Flight control testing of an unmanned airship model of 25 m long[J]. Aeronautical & Space Science Japan,2004, 52 : 30 –36.
Click to display the text
[17] SEMERARO N, GRIMALDI T, PESA M L, et al. New approach and results on the stability and control of airship[J]. Pathologica,2014, 62 (911) : 344 –354.
Click to display the text
http://dx.doi.org/10.13700/j.bh.1001-5965.2015.0427
北京航空航天大学主办。
0

文章信息

游颖捷, 糜攀攀, 吕明云, 刘文定
YOU Yingjie, MI Panpan, LYU Mingyun, LIU Wending
基于副气囊耦合的飞艇纵向运动影响分析
Airship's longitudinal motion impact analysis based on coupling ballonet
北京航空航天大学学报, 2016, 42(6): 1303-1310
Journal of Beijing University of Aeronautics and Astronsutics, 2016, 42(6): 1303-1310
http://dx.doi.org/10.13700/j.bh.1001-5965.2015.0427

文章历史

收稿日期: 2015-06-25
录用日期: 2015-09-18
网络出版时间: 2016-06-25 12: 00

相关文章

工作空间