廖俊 袁俊杰 蔣祎 楊澤川 李珺 盧智勇 吳春暉 王寧

?

高空零壓氣球上升過程的運動特性研究

廖俊 袁俊杰 蔣祎 楊澤川 李珺 盧智勇 吳春暉 王寧

（中南大學航空航天學院，長沙 410083）

飛行在20km~100km高度的臨近空間浮空器，在區(qū)域監(jiān)視、偵查、應急通信與科學實驗等領域具有獨特的優(yōu)勢。高空零壓氣球是常見的臨近空間浮空平臺之一。對其上升過程運動特性的研究有利于增加駐空時間、減少橫向漂移距離。文章通過建立零壓氣球充氣量模型和氣球上升過程中的動力學模型，對氣球上升過程中動力學模型進行仿真分析，結合高空飛行試驗數(shù)據(jù)，驗證了所建模型的準確性。在模型的基礎上對比了不同浮升氣體對氣球上升速度以及駐空高度的影響；研究了浮升氣體充氣量對載荷總量的影響，并分析了氣球在上升過程中的體積變化以及駐空高度與浮升氣體質量的關系，最后討論了影響氣球上升特性的關鍵因素。可以得到零壓氣球上升過程運動特性與浮升氣體種類、充入氣體量、氣球規(guī)格、載荷質量有關，假設氣球可無限膨脹不發(fā)生爆炸，浮升氣體密度越小，充入氣體量越多，載荷質量越輕，氣球的駐空高度就越高，達平衡所需時間越少；實際工程應用中，氣球受材料限制不可無限膨脹，面密度小的氣球，可能在上升到駐空高度前發(fā)生爆炸，要使氣球攜帶一定質量載荷在某一高度范圍駐空，其充氣量不能過多，需與一些物理量滿足一定關系。

臨近空間 零壓氣球 運動特性 飛行試驗

0 引言

高空氣球又稱高空科學氣球，是指在臨近空間飛行的無動力浮空器[1]。臨近空間（Near Space）是指距離地面20~100km的空域，開發(fā)利用該區(qū)域將對拓展戰(zhàn)略空間、維護國家安全、推動科學社會發(fā)展等產(chǎn)生深遠影響，其重大開發(fā)價值在國際上引起了廣泛關注[2-5]。

高空氣球作為一種可到達臨近空間的實用浮空器，對其駐空特性的研究在開發(fā)廉價浮空器應用價值方面具有重要意義。由于高空氣球一般不裝備動力裝置,因此其上升過程的運動特性受浮力、風場力和大氣阻力影響，不同高度下大氣溫度、壓力和密度有很大差異，從而影響氣球上升過程中的浮力。高空氣球的駐留特性受氣球規(guī)格、充氣量的影響。文獻[6]計算了標準大氣模型下超壓氣球的上升軌跡和蒙皮及氦氣的平均溫度；文獻[7]基于對高空氣球熱力學環(huán)境分析，建立了熱力學與動力學耦合的高空氣球動力學模型并對高空氣球的上升與駐留過程進行仿真分析；文獻[8-9]在大氣參數(shù)和熱力學環(huán)境變化的基礎上,建立了氣球上升過程中的數(shù)學模型,分析了放飛參數(shù)對氣球上升過程的影響；文獻[10]建立了氣球運動的動力學方程和熱力學模型，考慮了影響預測氣球運動軌跡的一些不確定性因素；應用文獻[11]對氣球運動軌跡進行了預測；文獻[12]在動力學基礎上運用MATLAB編寫了零壓氣球高度控制程序；文獻[13]將風場數(shù)據(jù)引入到氣球熱動力模型中，可以初步預測高空氣球在水平方向上的運動軌跡；文獻[14]將平流層零壓氣球幾何模型、熱力學模型和運動學模型耦合到一起對氣球上升、駐空和下降過程的飛行性能進行了數(shù)值仿真。

本文基于前人的分析方法和研究成果，對零壓氣球上升過程進行建模和仿真分析，與所做的飛行試驗進行對比，驗證了模型預測的準確性，并通過仿真結果對影響氣球上升特性的關鍵因素作了歸納和分析。

1 高空氣球建模

平流層高空氣球是一種無動力飛行器，其升空是依靠氣球囊體內(nèi)的浮升氣體產(chǎn)生的浮力，而浮力大小與高空大氣密度、氣球囊體體積和重力加速度均有關系。因此，建立載荷質量與大氣密度、浮升氣體密度、浮升氣體體積等參數(shù)之間的關系尤為重要[15]。

1.1 幾何模型

如圖1所示，氣球在空中呈水滴形，設氣球俯視狀態(tài)下直徑為b；體積（或充入氣體體積）為b；高度為b；表面積為s；錐線長度為g，氣球設計的最大體積為D，氣球幾何參數(shù)如表1所示[14]。

1.2 充氣量模型

（1）充氣量下限

氣球升空是依靠往氣球囊體內(nèi)充入比空氣密度小的浮升氣體而產(chǎn)生浮力。氣球所受到的浮力由氣球內(nèi)氣體體積和氣球內(nèi)外氣體的密度差決定，根據(jù)浮力公式

表1 氣球幾何參數(shù)

Tab.1 Geometric parameters of Balloons

圖1 氣球幾何形狀

式中為氣球受到的浮力；a為空氣密度；

要使氣球上升，系統(tǒng)受到的浮力必須大于系統(tǒng)總重力，即

式中為系統(tǒng)總質量；g為球內(nèi)浮升氣體的質量。

氣球上升所需要充入的浮升氣體體積b需滿足如下條件：

式中g為浮升氣體密度。

（2）充氣量上限

由理想氣體狀態(tài)方程可知，氣球體積隨高度的上升而增大，每種氣球都有最小面密度，因此，若充入浮升氣體量過多時，會導致氣球未達到駐空高度前因體積過大而爆炸。

設氣球質量為b，氣球面密度可表示為

設氣球的最小面密度為min，帶入式（4），則氣球最大體積max可由下式表示

設氣球的駐空高度為m，則氣球要上升至高度m所需充入的最大浮升氣體體積bmax可表示為

式中m,m分別為m高度對應的大氣壓強、溫度；g為浮升氣體的摩爾質量；g為標準大氣壓下浮升氣體的密度。

由上述分析可知，充入浮生氣體體積b滿足如下關系：

1.3 運動模型

由式（7）可確定充入浮升氣體的體積b，在忽略上升過程中氣體泄露的情況下，氣體物質的量保持不變。由于大氣溫度、壓強隨高度的變化發(fā)生顯著變化，所以氣球上升過程中任意高度時的體積可表示為

（8）

式中0、0為氣球放飛時當?shù)氐臍鈮汉蜌鉁?；a、a為氣球在高度時對應的大氣壓強和溫度。

氣球受到的升力可表示為

式中 Δ為氣球球內(nèi)氣體與外界大氣的壓差[17]；p為載荷質量；為理想氣體常數(shù)

不考慮氣球在水平方向上的運動，即氣球在豎直方向上受到自身重力、升力和大氣阻力作用，根據(jù)牛頓第二定律可得：

式中為氣球上升的高度；為時間；為氣球豎直方向上的速度；為空氣阻力；0為氣球附加質量，考慮了被氣球拖動，隨氣球一起運動的那一層大氣。

式中d為阻力系數(shù)，對于球體，d一般介于0.10~0.47之間；b為氣球垂直投影面積；bv表示虛擬質量系數(shù)，假定值介于0.25~0.50之間[7]。同理，氣球在水平方向上的運動微分方程可表示為

式中、分別為水平面內(nèi)和方向的位移；x、y分別為水平面內(nèi)和方向的速度；x、y分別為水平面內(nèi)氣球受到的和方向的合力。

2 仿真結果與試驗驗證

為了使高空零壓氣球在20km高空駐留并驗證上述數(shù)學模型的正確性，進行了一系列飛行試驗，試驗大致包括充氣、連接、放飛、觀察數(shù)據(jù)4個部分。載荷艙放置探空儀，配合探空雷達用于測量氣象數(shù)據(jù)和跟蹤艙體位置。

浮空平臺各組成部分質量如表2所示：

表2 系統(tǒng)各部分質量

Tab.2 Total amount of each part of the system

使用上述數(shù)學模型在MATLAB環(huán)境下編程，所得到氣球上升過程中垂直方向飛行軌跡和水平方向飛行軌跡預測曲線，如圖2所示，并與試驗結果進行對比分析：

實際飛行試驗由于要達到氣球在確定高度范圍駐留的目的，使用了某種高度控制技術使氣球到達一定高度之后損失部分浮力，使氣球下降至駐空高度并維持駐空狀態(tài)，該過程在本文建模仿真過程中不予體現(xiàn)。故圖中仿真飛行軌跡與實際上升過程飛行軌跡基本吻合，證明了氣球上升過程動力學模型的準確性。

誤差分析：

1）由于阻力系數(shù)與大氣密度、壓強、溫度等因素都有關系，是一個動態(tài)量，在本文所建立的模型中，阻力系數(shù)取的是一個定值，所以在氣球上升過程中高度隨時間變化曲線的斜率有所差異，但總體能反映出氣球上升過程中的運動特性。

2）由于本文建模所基于的大氣參數(shù)公式與飛行試驗時的大氣參數(shù)有誤差，因此建立的氣球動力學模型得到的仿真結果與實際飛行試驗的數(shù)據(jù)存在差異，但在誤差允許的范圍內(nèi)，具有一定參考價值。

圖2 飛行試驗數(shù)據(jù)與仿真結果

3 影響氣球上升軌跡的關鍵因素分析

由上述分析可知，載荷質量決定充氣量的下限值，充氣量、浮升氣體種類（氫氣與氦氣）決定了浮力的大小，從而影響了氣球的上升速度及氣球在上升過程中的體積，由于壓強隨高度的增加而下降，氣球體積將在上升過程中逐漸增大，由于氣球規(guī)格的不一樣，不同質量的氣球爆炸時膨脹到達的體積不同，因此所能到達的極限高度也不一樣。故充氣量、氣球規(guī)格（蒙皮質量）、浮升氣體種類、載荷質量都對氣球上升過程中的各運動參數(shù)如上升速度和駐空高度有極大的影響。

3.1 浮升氣體種類的影響

圖3為氣球上升過程的軌跡和上升速度變化仿真結果，從圖中可以看出，充氣量、載荷質量、氣球規(guī)格相同時，浮生氣體密度越小，氣球收到的凈舉力越大，上升過程中的速度和加速度越大，到達平衡高度的時間越短。若不考氣球最小面密度的影響，充入密度小的浮升氣體，氣球所能到達的平衡高度也越高。

3.2 充氣量的影響

圖4中的三條曲線分別表示充入1m3、2m3和3m3氫氣的情況下，氣球體積隨高度的變化。明顯可以發(fā)現(xiàn)充入浮升氣體體積越多，相同高度上氣球體積越大。圖5表示充1m3氫氣、2m3氫氣和3m3氫氣時780g氣球升空過程中的運動特性，從圖中可以看出載荷質量相同時，充入浮升氣體量越多，上升速度越快，到達平衡高度所需的時間越短。

圖6為充入不同浮升氣體時，氣球駐空高度與充入氣體質量之間的關系，圖中可看出，在最小充氣量能帶起載荷和氣球球皮可無限延伸的條件下，充氣量越多，氣球駐空高度越高，充氣量相同的情況下，氣球充入氫氣時最后的駐空高度高于充入氦氣時最后的駐空高度。

3.3 氣球規(guī)格（質量）的影響

圖7表示350g、750g和1 600g氣球上升過程中面密度的變化。很顯然，相同的高度上，氣球的面密度與質量成正比。圖8表示充入1.6m3氫氣時，350g氣球、750g氣球和1 600g氣球升空過程中的運動特性，從圖中可以發(fā)現(xiàn)充入浮升氣體量相同時，氣球蒙皮質量越大，平衡高度越低，到達平衡高度所需的時間越長。

圖3 浮升氣體種類對氣球上升特性的影響

圖4 氣球體積隨高度變化

圖5 不同充氣量對氣球上升特性的影響

圖6 不同質量的氫氣和氦氣所能達到的駐空高度

圖7 不同質量的氣球面密度隨高度變化

Fig.7 The surface density of balloons in different model varies with height

圖8 氣球規(guī)格對上升特性的影響

3.4 載荷質量的影響

圖9表示充氫氣量相同時載荷為400g、600g和1 000g時，氣球升空過程中的運動特性，從圖中可以看出，充氣量相同時，載荷質量越大，上升速度越慢，上升加速度越小，不考慮氣球最小面密度限制的情況下，氣球駐空高度越低。

圖9 載荷質量對上升特性的影響

4 結束語

本文將大氣模型、零壓氣球幾何模型和動力學模型耦合到一起，得到氣球充氣量的取值范圍，對氣球升空過程中的飛行特性進行了數(shù)值仿真，并以飛行試驗驗證了模型的正確性，最后通過模型分析了影響零壓氣球飛行特性的關鍵因素，主要得出以下結論：1）零壓氣球上升過程中的運動特性與浮升氣體種類、充入浮升氣體量、載荷質量有關，浮升氣體密度越小，充入浮升氣體量越多，載荷質量越輕，氣球上升過程中的速度和加速度就越大。2）在不考慮氣球最小面密度，及認為氣球可無限膨脹的情況下，氣球的平衡高度與浮升氣體種類、充入浮升氣體量、載荷質量有關，浮升氣體密度越小，充入浮升氣體量越多，載荷質量越輕，氣球的駐空高度就越高。3）由于在實際應用中，氣球受蒙皮材料的限制，存在一個最小面密度，不可能一直上升，若充入浮升氣體的量過多，可能會在到達駐空高度前因體積過大而爆炸，因此，要使懸掛一定載荷的零壓氣球在某一高度范圍駐空，其充氣量需滿足一定關系。

[1] 李小建. 臨近空間浮空器熱—結構耦合數(shù)值模擬研究[D]. 南京: 南京航空航天大學, 2013. LI Xiaojian. Numerical Simulation of Thermal-structure Coupling for Near Space Airship[D]. Nanjing: Nanjing University of Aeronautics & Astronautics, 2013. (in Chinese)

[2] CARICHNER G E, NICOLAI L M. Fundamentals of Aircraft & Airship Design, Volume 2-Airship Design & Case Studies[M]. Washington DC: American Institute of Aeronautics & Astronautics, 2013. (in Chinese)

[3] 田莉莉, 方賢德. NASA高空氣球的研究及其進展[J]. 航天返回與遙感, 2012, 33(1): 81-87. TIAN Lili, FANG Xiande. Research & Progress of NASA’s Balloon[J]. Spacecraft Recovery & Remote Sensing, 2012, 33(1): 81-87. (in Chinese)

[4] 趙達, 劉東旭, 孫康文, 等. 平流層飛艇研制現(xiàn)狀、技術難點及發(fā)展趨勢[J]. 航空學報, 2016, 37(1): 45-56. ZHAO Da, LIU Dongxu, SUN Kangwen, et al. Research Status, Technical Difficulties & Development Trend of Stratospheric Airship[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(1): 45-56. (in Chinese)

[5] 陳旭. 長時高空氣球的研制發(fā)展[J]. 航天器工程, 2007, 16(4): 83-88. CHEN Xu. Development of Long Duration High Altitude Balloon[J]. Spacecraft Engineering, 2007, 16(4): 83-88. (in Chinese)

[6] 田莉莉, 方賢德, 戴秋敏, 等. 高空氣球運動特性研究[C]//中國浮空器大會. 荊門, 2011. TIAN Lili, FANG Xiande, DAI Qiumin, et al. Study on the Motion Characteristics of High-altitude Balloons[C]// China Aerostat Conference. Jingmen, 2011. (in Chinese)

[7] 呂明云, 巫資春. 高空氣球熱力學模型與上升過程仿真分析[J]. 北京航空航天大學學報, 2011, 37(5): 505-509. LYU Mingyun, WU Zichun. Thermodynamic Model & Numerical Simulation of High Altitude Balloon Ascending Process[J]. Journal of Beijing University of Aeronautics & Astronautics, 2011, 37(5): 505-509. (in Chinese)

[8] SALEH S. Floating Performance Analysis & Extended Lifetime for High Altitude Zero Pressure Balloon[J]. International Journal of Modeling & Optimization, 2016, 6(4): 6-4.

[9] SALEH S, HE W. Ascending Performance Analysis for High Altitude Zero Pressure Balloon[J]. Advances in Space Research, 2017, 59(8): 2158-2172.

[10] LEE Y, YEE K. Numerical Prediction of Scientific Balloon Trajectories while Considering Various Uncertainties[J]. Journal of Aircraft, 2016, 54(2): 1-15.

[11] SóBESTER A, CZERSKI H, ZAPPONI N, et al. High-altitude Gas Balloon Trajectory Prediction: A Monte Carlo Model[J]. AIAA Journal, 2015, 52(4): 832-842.

[12] KAYHAN O, YUCEL O, HASTAOGLU M A. Control & Simulation of High Altitude Zero Pressure Balloons[C]// International Conference on Engineering & Natural Sciences. Bali, Indonesia, 2017.

[13] PALUMBO R, RUSSO M, FILIPPONE E, et al. ACHAB: Analysis Code for High-altitude Balloons[C]//AIAA. Atmospheric Flight Mechanics Conference and Exhibit, Hilton Head Island, SC, AIAA Paper 6642, 2007. DOI:10.2514/6.2007-6642.

[14] 蘇潤, 李小建. 平流層零壓氣球飛行控制仿真研究[J]. 合肥工業(yè)大學學報, 2018(3): 325-332.SU Run, LI Xiaojian. Numerical Simulation Study on Flight Performance of Stratospheric Zero-pressure Balloon[J], Journal of Hefei University of Techology, 2018(3): 325-332. (in Chinese)

[15] 戴秋敏, 方賢德, 王昊, 等. 大氣模型對高空氣球運動特性和熱特性的影響[J]. 計算機仿真, 2013, 30(9): 79-82.DAI Qiumin, FANG Xiande , WANG Hao, et al. Simulation of Influence of Different Atmospheric Models on Dynamic & Thermal Properties of High Altitude Balloons[J]. Computer Simulation, 2013, 30(9): 79-82. (in Chinese)

[16] 楊炳尉. 標準大氣參數(shù)的公式表示[J]. 宇航學報, 1983, 4(1): 86-89.YANG Bingwei. Formulization of Standard Atmospheric Parameters[J]. Journal of Astronautics, 1983, 4(1): 86-89. (in Chinese)

[17] 劉東旭, 樊彥斌, 馬云鵬, 等. 氦氣滲透對高空長航時浮空器駐空能力影響[J]. 宇航學報, 2010, 31(11): 2477-2482.LIU Dongxu, FAN Yanbin, MA Yunpeng, et al. Effect of Helium Permeability on Working Endurance High Altitude Long Duration LTA Vehicle[J]. Journal of Astronautics, 2010, 31(11): 2477-2482. (in Chinese)

Motion Characteristics of Zero-pressure Balloon in Ascending Process

LIAO Jun YUAN Junjie JIANG Yi YANG Zechuan LI Jun LU Zhiyong WU Chunhui WANG Ning

（School of Aeronautics & Astronautics, Central South University, Changsha 410083, China）

Near-space aerostats, which flying at an altitude from 20 kilometers to 100 kilometers, have advantages in the fields of local monitoring, reconnaissance, emergency communication and scientific experiment, etc. high-altitude zero-pressure balloon a familiar kind of near space aerostat the study of its movement characteristics is beneficial to increase its resident time and decrease lateral drift in space. In the present study, a dynamic model of zero-pressure balloon's ascending process is built and simulation of the ascending process is conducted. Accuracy of the model is confirmed with the data from actual flight test. This study focuses on how ascending rate and resident altitude are effected when using different lifting gas. This study also focuses on how weight of payload is effected by different volume of lifting gas and analyzes the relationship between the variation of balloon's volume in ascending process and the resident altitude or the mass of lifting gas. The critical factors that effect the character of balloon's ascension are discussed. As a result, The motion characteristics of the zero-pressure balloon in the rising process are related to the type of lifting gas, the amount of gas filled, the balloon specifications and the mass of payload. With the assumption that the balloon can expand indefinitely without exploding, the smaller the density of lifting gas is, the more the amount of lifting gas is filled, the lighter the payload mass is, the higher the resident altitude of the balloon, and the less time it takes to reach the equilibrium height. In practical engineering applications, a balloon cannot expand indefinitely because of its material, a balloon with a low surface density may explode before rising to the resident altitude. In order to make the balloon carry a certain weight payload floating in a certain height range in the air, the volume of lifting gas cannot be too much, requires to meet a certain relationship with some quantities.

near space; zero-pressure balloon; motion characteristics; flight test

V211.5

A

1009-8518(2019)01-0011-09

10.3969/j.issn.1009-8518.2019.01.002

廖俊，男，1982年生，博士，副教授。研究方向為飛行器總體設計，飛行器結構動力學分析，新概念浮空型飛行器設計。E-mail：liaojun@csu.edu.cn。

袁俊杰，男，1995年生，碩士生。研究方向為飛行器總體設計。E-mail：yuanjunjie@csu.edu.cn。

李珺，男，1989年生，博士，副教授。研究方向為浮空器設計。E-mail：lijun215@csu.edu.cn。

2018-09-12

（編輯：劉穎）