极涡风场中平流层浮空器轨迹仿真研究

1、极涡风场中平流层浮空器轨迹仿真研究摘要：为实现平流层浮空器绕极地飞行，结合极地地区的环境特征，建立了极地风场模型与太阳辐照模型、平流层浮空器动 力学模型以及太阳电池与蓄电池相结合的循环能源模型。利用极地地区特有的极涡风场，设计了平流层浮空器轨迹控制策 略，以纬向风作为动力在东西方向上飞行，考虑风场紊流影响，通过螺旋桨对浮空器在南北方向上的偏移施加控制。仿真结 果表明，在推进功率等约束条件下，采用PID控制方法的浮空器在受到干扰产生偏移时仍然能够回到预定的轨迹，证明了利 用极涡风场和控制策略实现平流层浮空器绕极地飞行的可行性。关键词：轨迹仿真;极涡风场;平流层浮空器Trajectory Simu

2、lation of Stratosphere Aerostatsin Polar Vortex Wind FieldABSTRACT ：In order to achieve flying around the polar region, polar wind field model, the solar irradiation model, dynamic model, and energy model are established in this paper based on the environmental characteristics of the polar region. T

3、he trajectory control strategies of stratosphere aerostat were studied with the considering of polar vortex wind field and turbulence. The aerostats flew with latitude wind in the east and west direction and were controlled in the north and south direction by propellers. The simulation results revea

4、l that under the constraint of propeller power and other constraints, the aerostats can maintain the trajectory with the PID control method under the disturbance of turbulence ,which demonstrates the conceptual feasibility of flying around the polar region of the stratosphere aerostats using polar v

5、ortex wind field and control strategies.KEYWORDS: Trajectory simulation ； Polar vortex wind field ； Stratosphere aerostat1引言由于其独特的地缘战略价值和丰富的战略资源储备，极 地地区是当前世界各国争相开发和利用的前沿热点。现 有的极地勘探手段众多，主要包括极地科考船平台、卫星遥 感平台、极地科考站、无人机平台和浮标平台等。但是，它们 均存在不同程度上的局限性，如船舶破冰能力有限,受冰面 和天气限制严重;卫星过顶时间短，资源昂贵;科考站受天气 限制，需人员长期驻留，且覆盖范围

6、有限；飞机续航能力不 足，不能持久驻空等。为了巩固在极地地区的利益，世界各 国均在大力发展新型极地探测与保障手段。平流层浮空器 作为承载平台具备长航时大范围信息获取的能力，在极地地 区的持久大范围信息保障方面具有巨大潜力2-4。研究表明，极涡是一种存在于极地地区中上层对流层和 平流层大气的持续性大尺度气旋性环流5，利用极涡这种持 续性、稳定性、气旋性的特点，平流层浮空器可以实现大尺度 绕极飞行。与此同时，极地存在显著的极昼现象，这对于 利用太阳能飞行的平台而言，可大幅降低对储能电池的需 求，缩小系统运行规模，提高载荷能力和应用效费比。对于在极涡风场中工作的平流层浮空器，由于风场紊流 对绕极飞行

7、的干扰，其面临的一个重要问题就是在设计约束 下的大尺度轨迹保持能力。平流层浮空器在风场环境中的轨迹仿真与控制长期以来受到学者们关注。Kayhan等牝通 过控制气体压缩释放系统进行平流层浮空器的高度调控,分 析了平流层浮空器在准零风层的轨迹保持能力。Born等时 提出一种预测规划算法,使得平流层浮空器在只能垂直机动 的情况下，也能保持横向的轨迹。Ramesh等9比较了带有 控制装置的双气球系统和气球-帆系统在风场中的轨迹保 持能力。李魁等通过平流层风场预测建模，分析了预 测模型对浮空器轨迹控制的影响，研究了平流层浮空器区域 驻留策略。这些研究均为极涡环境下的平流层浮空器轨迹 控制提供了参考。本文

8、针对极地这一特定区域，对平流层浮空器在极地风 场环境下的飞行轨迹进行仿真分析。同时，针对极昼现象对 能源系统的影响,本文分析了在极昼条件下平流层浮空器能 源系统的功率变化情况,研究了功率限制条件下水平方向控 制策略，为平流层浮空器在极地地区的应用提供支撑。2极地风场模型平流层浮空器的充足可靠的动力是制约其应用的一个 瓶颈4。考虑到极地地区普遍存在极涡风场的规律性，利用 极地环流可以实现浮空器绕极地无动力飞行，降低对平台动 力系统的技术需求。下面建立极地风场的模型。参照国家军用标准GJB 366.2 -8712在20km处的大气 风场数据，选取南纬65度到90度之间的风速进行分析，得 到一月份不

9、同经纬度在20km处的风速，由于极地环流的特 殊性，在南极圈内同一纬度,不同经度处的风速相差不大，为 了便于进行仿真,本文在同一纬度选取了四个不同经度处的 风速取平均值，以此作为该纬度的平均风速。见表1。表1 一月份20km高度处的平均风速经度65N70N75N80N85N90N80W31241714.5141420E292722. 519181780E32经度65N70N75N80N85N90N80W31241714.5141420E292722. 519181780E322721171413.5140E3031.526231817均值30.527.421. 618.41615.4不同的纬度

10、虽然极涡风场在极地区域普遍存在，但该风场并不是理 想的绕极环流，而是存在不同尺度的风场扰动等现象。例 如,NASA在南极开展的高空科学气球试验发现，气球通常在 绕南极飞行数圈后回到放飞点附近13-15,但是由于大气紊 流的存在，气球飞行轨迹存在较大不确定性，气球在没有控 制装置的情况下，其实际飞行轨迹难以得到保证，存在飞出 极地区域的可能。因此,对于绕极持续工作的平流层浮空器 38 图1南纬65度到90图1南纬65度到90度之间一月份海拔20km处的平均风速本文利用MATLAB产生符合标准正态分布的随机数，用 于模拟风场紊流,可以较好地模拟平流层浮空器在大气紊流 中受到的侧向干扰力。3平流层浮

11、空器动力学模型平流层浮空器在飞行过程中，水平方向主要受到风场的 作用力以及推进系统提供的控制力，而竖直方向上主要受到 浮升气体的浮力以及浮空器自身的重力。平流层浮空器气 囊形状的确定需要综合考虑空气动力学，结构强度等因 素2，为了便于仿真，本文采用球形气囊进行计算。以中国 南极中山站为坐标原点，建立东北天坐标系，可得到动力学 方程10mallX = Fdl + F峥maliy = Fd, + F岫 + Fw心=B + Fdz + 膈-maUg(1)式中为系统总质量;B为浮力;g为当地重力加速度; Fd A 后为气动力兀为附加惯性力;F斗为螺 旋桨动力系统的径向控制力。(1)Fdx = - Cd

12、*p财S*以 1Fd, = CdP/y% 1 % I(2)Fdz = jCdPefS&z 11(3)式中:P矿为当地空气密度；&，土为浮空器有效横截面积; 七二 - Mwind , ry ”,ywmd，rz z，Mwind(2)Fdz = jCdPefS&z 11(3)y轴方向风速;虬为气体总体积;当浮空器形状为球形时，阻力系数cdi &计算公式如下Cdl 1.456/ ( 2财 Er 广q, - 1.456/ (2气 I r广(4)q, - 1.456/ (2气 I r广(4)式中:r为浮空器半径计算公式如下T四 1.714 x 10 _5 x (- )4-256 f*273.15，(5)(

13、6)式中：孔为当地大气温度。(5)(6)浮力大小B与当地空气密度、浮空器内浮升气体体积和 当地重力加速度有关：B - PrefTg式中：当地重力加速度g计算公式如下R2(7)g (z + R)2g。(7)g (z + R)2g。式中:R为地球半径,go为地球表面重力加速度。能源分系统是浮空器系统持续运行的关键，针对采用的 太阳电池与锂离子电池的能源组合，分析太阳电池发电功率 的变化情况。太阳电池发电功率与太阳辐射强度有关，因 此，首先建立太阳辐照模型：真太阳时与平太阳时之间的差值为时差,时差的计算公 式如下e 二-0.0002786409 +0. 1227715cos( W + 1.49831

14、1)-0. 1654575cos(2V -1. 261546)-0. 005353830cos(3 W -1. 1571)式中：e表示时差；W = 2nn/365；n为日期序数。 当地真太阳时计算公式为L _ Lbe峋 + 15 + 60(8)(9)式中b表示该地区标准时间;L表示当地经度;Lb表示当地 标准时间位置的经度。(8)(9)真太阳时用角度来表示为当地太阳时角，公式如下W (t _ 12) x 15太阳赤纬可由下式确定(10)8 23 45 x sin(360 x 284 + n x n )(10)-x Ain(xx )365180(11)(12)对于太阳高度角h有(11)(12)s

15、in h sin sin 8 + cos ?cos8coss 式中淬为当地纬度。设Y为太阳方位角、0为斜面倾角、。为太阳入射角, 则有sin sin h 一 sin 8cos y cos cos h(13)cos 0 sin 8( sincoS 一 cossincosy)(13)+ cos8cosw( cos cos S 一 sinsiScosy)+ cos 8sin jSsinysinw( 14)设Io为大气层外太阳辐射强度，则有10 _ 1353 x 1 + 0. 033 x cos( 2n x n/365) x cos。(15)大气层内太阳辐射由三部分构成，分别为太阳直射辐 射、天空散射

16、辐射和地面反射辐射，由于可供太阳电池发电 的辐射主要为太阳直射辐射，下面主要对太阳直射辐射进行 建模。太阳直射辐射强度IDN是指在垂直于太阳光射线的表面 上，单位时间内投射到单位面积上的太阳直射辐射能量， 且有IDN I0T( 16 )式中:T为大气透过率，计算公式如下-0. 65仇，-0. 095mt 0.5 (e+ e )(17)式中:m为大气质量,计算公式如下m ( P ) J1229 + (614sin h)2 一 614sinh (18) P。式中:P为浮空器所在高度处的大气压强,P0为海平面大气 压强,h为太阳高度角。太阳电池吸收的太阳直射辐射能量Qd* 16为QdN nS IDN

17、dt(19)式中:n为太阳辐射吸收率;s为辐射面积。由于极地存在高纬度地区的极昼现象，与中低纬度地区 相比，在极昼期间有持续的太阳辐照，可为太阳电池提供持 续的能源。本文比较了南极中山站和长沙市在极昼期间同 一天内20km处太阳电池发电量的情况。计算得到,在长沙 地区,一天内太阳电池发电总量为7896. 5W - h,而南极中山 站发电总量为14537W - h,是长沙的1. 84倍。考虑通信、测 控、球载计算机等功耗，设置总功率需求为400W。太阳电池 功率变化曲线如图2所示，在中低纬度地区,太阳电池一天 内只能在一定的时间段工作，以长沙为例，太阳电池在上午7 点到下午18点这一区间内可以提

18、供400W以上的功率，其它 13个小时的时间则需要蓄电池为系统工作提供能源。而在 南极地区，由于存在极昼现象，太阳电池在一天内基本能满 足400W的功率需求，只在上午1点到5点的4个小时内需 要蓄电池辅助其提供能源，对蓄电池的工作时间需求仅为长 沙地区的30%。因此对南极地区的平流层浮空器而言，在设 计阶段便可降低储能电池的比重，减小系统规模与质量。5仿真结果极地地区由于存在环流，东西方向上风速较大，本文利 用极地环流作为浮空器在东西方向上飞行的主要动力。然 而,环流并不是规则地存在于东西方向，其在南北方向上也 存在分量,极地环流在宏观上表现为不规则的圆或者椭圆， 因此，在利用环流飞行的过程中

19、，需要排除环流在南北方向图2功率随时间的变化情况上的分量对平流层浮空器轨迹的影响,通过螺旋桨提供的控 制力，将浮空器保持在预定的轨迹上飞行。环流在东西方向上的分量是浮空器的主要动力来源，因 此在东西方向上不对浮空器施加控制力，只在南北方向上施 加控制力。根据实际轨迹与预定轨迹的偏差量，采用PID控 制使得浮空器回到预定轨迹。其中，螺旋桨控制力的表达 式为t七=郊（y。-y） + 农 Z（ -y） +&（ - %）（20） 0在simulink上建立如图3的浮空器仿真模型。5.1小时空尺度下的无控飞行图3极地平流层浮空器系统simulink模型时间人（C）图4无控制的浮空器纬向位移随时间的变化情

20、况仿真过程假定在南北方向上受到一股持续100秒的大 小为10m/s的干扰风作用，整个仿真过程持续1000秒，当平 流层浮空器没有轨迹控制时，其轨迹仿真图如图4所示，在 东西方向上，由于仿真过程设置的风速不变，因此，纬向位移 与时间呈线性关系;在南北方向上，受到短暂的干扰风作用, 经向位移先是随时间快速增大，当干扰风消失之后，由于阻 力的作用，经向速度逐渐减小，图3极地平流层浮空器系统simulink模型时间人（C）图4无控制的浮空器纬向位移随时间的变化情况由仿真结果可以看出，不施加控制的时候，平流层浮空 器受到侧向干扰风力作用，会产生不可逆的偏移。5.2小时空尺度下的有控飞行平流层浮空器在南北方向上受到持续100秒的大小为 10m/s的干扰风作用，仿真过程持续1000秒，当浮空器在南 北方向上有螺旋桨提供的控制力时，其轨迹仿真图如图5所 示，在东西方向上，由于风速不变，平流层浮空器纬向位移与 无控制的浮空器一样，都是与时间呈线性关系的；在南北方 向上，受到短暂的干扰风作用后，经向位移先是随时间快速 增大,但同时，螺旋桨的控制力也开始产生作用，当干扰风消 失之后，在螺旋桨控制力的作用下，平流层浮空器能够在一 定时间内回到预定轨迹上。图5有控制的浮空器纬向位移随时间的变化情况由仿真结果