CN108241769A - 高空气球航迹的预测方法及装置 - Google Patents
高空气球航迹的预测方法及装置 Download PDFInfo
- Publication number
- CN108241769A CN108241769A CN201611229475.3A CN201611229475A CN108241769A CN 108241769 A CN108241769 A CN 108241769A CN 201611229475 A CN201611229475 A CN 201611229475A CN 108241769 A CN108241769 A CN 108241769A
- Authority
- CN
- China
- Prior art keywords
- aerostat
- model
- current state
- trajectory prediction
- field data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Development Economics (AREA)
- Geometry (AREA)
- Game Theory and Decision Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明公开了一种高空气球航迹的预测方法及装置。其中,该方法包括:建立高空气球的航迹预测模型,其中,航迹预测模型至少包括:动力学模型和热力学模型;获取高空气球的当前状态和风场数据,其中,当前状态至少包括:位置信息、速度信息、氦气量、温度和压强,风场数据至少包括:风向数据和风速数据;根据当前状态和风场数据,通过航迹预测模型,得到高空气球的航迹。本发明解决了现有技术中的高空气球航迹的预测方法无法同时满足简便预测和高精确度的要求的技术问题。
Description
技术领域
本发明涉及高空气球飞行力学与仿真分析领域,具体而言,涉及一种高空气球航迹的预测方法及装置。
背景技术
平流层高空气球的研究越来越受到重视,考虑到高空气球放飞、平飞以及回收过程中,需要对其航迹进行控制,特别是落点位置十分重要,必须进行航迹实时预测研究工作。
目前,航迹预测模拟方法主要包括:记录空中真实飞行器航迹和基于地球球体模型进行的航迹计算,第一种方法精确但是不简便,第二种方法简便却不够精确。
针对现有技术中的高空气球航迹的预测方法无法同时满足简便预测和高精确度的要求的问题,目前尚未提出有效的解决方案。
发明内容
本发明实施例提供了一种高空气球航迹的预测方法及装置,以至少解决现有技术中的高空气球航迹的预测方法无法同时满足简便预测和高精确度的要求的技术问题。
根据本发明实施例的一个方面,提供了一种高空气球航迹的预测方法,包括:建立高空气球的航迹预测模型,其中,航迹预测模型至少包括:动力学模型和热力学模型;获取高空气球的当前状态和风场数据,其中,当前状态至少包括:位置信息、速度信息、氦气量、温度和压强,风场数据至少包括:风向数据和风速数据;根据当前状态和风场数据,通过航迹预测模型,得到高空气球的航迹。
根据本发明实施例的另一方面,还提供了一种高空气球航迹的预测装置,包括:建立单元,用于建立高空气球的航迹预测模型,其中,航迹预测模型至少包括:动力学模型和热力学模型;获取单元,用于获取高空气球的当前状态和风场数据,其中,当前状态至少包括:位置信息、速度信息、氦气量、温度和压强,风场数据至少包括:风向数据和风速数据;处理单元,用于根据当前状态和风场数据,通过航迹预测模型,得到高空气球的航迹。
在本发明实施例中,建立高空气球的航迹预测模型,获取高空气球的当前状态和风场数据,根据当前状态和风场数据,通过航迹预测模型,得到高空气球的航迹。容易注意到的是,通过航迹预测模型预测高空气球的航迹,可以简化预测过程,将根据当前状态和风场数据作为航迹预测模型的输入,通过高空气球的真实飞行航迹更新实时预测后的航迹,可以减小航迹预测模型预测的误差,解决了现有技术中的高空气球航迹的预测方法无法同时满足简便预测和高精确度的要求的技术问题。因此,通过本发明上述实施例提供的方案,可以达到简便易用,精确预测的效果。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是根据本发明实施例的一种高空气球航迹的预测方法的流程图;
图2是根据本发明实施例的一种可选的高空气球航迹的预测方法的示意图;以及
图3是根据本发明实施例的一种高空气球航迹的预测装置的示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
实施例1
根据本发明实施例,提供了一种高空气球航迹的预测方法实施例,需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
图1是根据本发明实施例的一种高空气球航迹的预测方法的流程图,如图1所示,该方法包括如下步骤:
步骤S102,建立高空气球的航迹预测模型,其中,航迹预测模型至少包括:动力学模型和热力学模型。
具体的,高空气球可以是平流层高空气球,由一个浮升囊体和挂接在囊体下端的吊舱组成,囊体内充入氦气,以高空气球所需要的浮力。
在一种可选的方案中,动力学模型是航迹预测的基础,高空气球的飞行航迹需要由动力学模型获得。高空气球是一个气体浮空器,其浮力大小受热环境影响,氦气在球体中满足理想气体状态方程为:PV=nRT,其中,P为氦气的压强,V为氦气的体积,即球体体积,n为氦气的物质的量,R气体常数,T为氦气的热力学温度。在上升过程中,囊体体积会随温度变化而变化,在囊体到达最大体积后,温度继续增大就会产生超压量,所以需要建立高空飞行过程中的热力学模型,得到囊体内的氦气的温度变化。因此,在对高空气球进行航迹预测之前,需要建立高空气球的动力学模型和热力学模型。
步骤S104,获取高空气球的当前状态和风场数据,其中,当前状态至少包括:位置信息、速度信息、氦气量、温度和压强,风场数据至少包括:风向数据和风速数据。
步骤S106,根据当前状态和风场数据,通过航迹预测模型,得到高空气球的航迹。
在一种可选的方案中,高空气球在飞行过程中,由于其体型巨大,受到风的作用影响巨大,航迹预测需要有准确的风速风向信息,在建立高空气球的动力学模型和热力学模型之后,可以获取高空气球的当前时刻的真实状态和风场数据,通过动力学模型与热力学模型预测之后时间的飞行轨迹即航迹。
根据本发明上述实施例,建立高空气球的航迹预测模型,获取高空气球的当前状态和风场数据,根据当前状态和风场数据,通过航迹预测模型,得到高空气球的航迹。容易注意到的是,通过航迹预测模型预测高空气球的航迹,可以简化预测过程,将根据当前状态和风场数据作为航迹预测模型的输入,通过高空气球的真实飞行航迹更新实时预测后的航迹,可以减小航迹预测模型预测的误差,解决了现有技术中的高空气球航迹的预测方法无法同时满足简便预测和高精确度的要求的技术问题。因此,通过本发明上述实施例提供的方案,可以达到简便易用,精确预测的效果。
可选地,在本发明上述实施例中,步骤S102,建立高空气球的航迹预测模型至少包括:
步骤S1022,根据高空气球的动力学信息,得到动力学模型,其中,动力学信息至少包括:高空气球的质量、高空气球的参考面积和高空气球的速度。
具体的,上述的高空气球的质量可以包括整个高空气球的结构材料的质量,囊体内充入的氦气质量和附加空气质量,上述的高空气球的速度可以包括高空气球在垂直方向速度和水平方向速度。
在一种可选的方案中,可以采用单刚体动力学建模方法,可以将高空气球考虑成三自由度,不考虑姿态,也可以考虑六自由度刚体建模,受力分为重力,浮力,阻力,附加质量力,控制力。可以将东、北、天三个方向分别定义为地面惯性坐标系的x,y,z正向,即向东为x正向,向北为y正向,向上为z正向。对于高空气球,飞行过程中,受到的作用力主要有浮力,重力,附加质量力,气动力,由于高空气球的姿态变化对于受力影响较小,可以忽略球的姿态变化而将其看成质点,故竖直方向其质点动力学模型如下:
其中,M为整个系统的质量,包括结构材料质量,氦气质量,附加空气质量,充放气过程中h为高空气球的高度,Fhe为高空气球所受的浮力,Fdz为高空气球垂直方向所受的空气阻力,G为高空气球所受的重力。
其竖直速度与高度的关系如下式所示:
空气阻力大小与空气密度ρ,垂直方向的阻力系数Cdz,囊体参考面积S,高空气球的垂直速度vz(忽略垂直风)有关,Fdz如下式所示:
重力与高空气球重量M,当地重力加速度g有关,重力加速度随高度变化,G如下式所示:
G=Mg。
浮力与空气密度ρ,囊体体积V,重力加速度g有关,Fhe如下式所示:
Fhe=ρVg。
水平面受到的力是风引起的空气阻力,动力学方程如下所示:
其中,Cdx和Cdy分别为x和y方向的阻力系数,vx和vy分别为高空气球在x和y方向的速度,vxwind和vywind分别为风在x和y方向的速度。
水平方向的运动学方程如下所示:
步骤S1024,根据高空气球的热环境信息,得到热力学模型,其中,热环境信息至少包括:囊体光学特性参数和热环境参数。
具体的,上述的囊体光学特性参数可以包括:囊皮材料对太阳光的吸收率和透射率,高空气球的照射投影面积,囊体的外露面积,囊体蒙皮材料对红外辐射的吸收率和透射率,蒙皮的红外发射率,蒙皮平均温度,外部和内部对流换热系数,有效对流面积,氦气的平均温度,蒙皮比热容,蒙皮质量,氦气的定压比热容和定体积比热容,囊体体积;上述的热环境参数可以包括:太阳辐射强度,大气对太阳光的透射率,太阳高度角,地面平均反照率,高空气球表面对地球表面的角系数,地面平均发射率,大气对地面红外辐射的透射率,地面温度,大气红外辐射强度,空气的平均温度。
在一种可选的方案中,热力学模型包含太阳直射辐射、太阳反射、大气散射、地面红外辐射、蒙皮红外辐射、对流换热等热源,可以根据实际状况获取囊体光学特性参数以及热环境参数,并根据获取到的热环境信息得到热力学模型。
可选地,在本发明上述实施例中,步骤S1024,根据高空气球的热环境信息,得到热力学模型包括:
步骤S10242,根据高空气球的热环境信息,得到高空气球的热源模型,其中,热源至少包括:太阳直射辐射、太阳反射、大气散射、地面红外辐射、蒙皮红外辐射、大气红外辐射和对流换热。
在一种可选的方案中,热源主要有以下几种:
(1)直接太阳辐射,对直接太阳辐射的吸收包括外表面的吸收和透过蒙皮由内表面的吸收,而内表面由于多次反射的作用会提高内部的有效吸收率。假设囊皮材料对太阳光的吸收率与透射率分别为α和τ,反射率为r=1-α-τ,有效反射率reffective=r+r2+r3+…,容易得到对太阳直接辐射的吸收为:
其中,Aproject为对高空气球的照射投影面积,由几何模型获得;ISun为太阳辐射强度;τatm是大气对太阳光的透射率。
(2)大气散射,对于大气散射,可以由下式得到:
其中,κ为大气散射经验系数,通常取0.1;Asurf为囊体的外露面积。
(3)地球反照辐射,对地面反照辐射吸收为:
其中,qAlbedo=AlbedoISun sin(ELV);ELV为太阳高度角,特别地,太阳高度角小于0时,qAlbedo=0;Albedo为地面平均反照率,与经度和当地气象条件有关,一般取0.3;ViewFactor为高空气球表面对地球表面的角系数,与高度有关。
(4)地面红外辐射,地面红外辐射强度可以由斯特潘-波尔茨曼定理计算:
其中,αIR和τIR分别为囊体蒙皮材料对红外辐射的吸收率和透射率;εground为地面平均发射率,取0.95,σ为斯特潘-波尔茨曼常数,τatmIR为大气对地面红外辐射的透射率,Tground为地面温度,蒙皮材料对红外辐射的反射率rIR=1-αIR-τIR。
(5)蒙皮的红外辐射,囊体蒙皮的内外表面均向外辐射热量,内表面辐射的热量会有部分被蒙皮吸收,故蒙皮的红外辐射包括蒙皮对外的红外辐射以及蒙皮吸收的红外辐射,可以由下式得到蒙皮对外的红外辐射QIRFilm和蒙皮吸收的红外辐射QIRout:
其中,ε为蒙皮的红外发射率,假设满足基尔霍夫定理,红外发射率等于红外吸收率;Tfilm为蒙皮平均温度。
(6)大气红外辐射,对于大气红外辐射,可以由下式得到:
其中,大气等效温度IIRsky为大气红外辐射强度。
(7)对流换热,对流换热包括外表面与大气环境的对流换热以及内表面与氦气间的对流换热,蒙皮的对流换热为:
QconvExt=HCexternalAeffective(Tair-Tfilm),
QconvInt=HCinternalAeffective(Tfilm-Tgas),
其中,HCexternal和HCinternal分别为外部和内部对流换热系数,可由经验公式获得;Aeffective为有效对流面积,由几何模型获得;Tair和Tgas分别为空气和氦气的平均温度。
步骤S10244,根据高空气球的热源模型,得到热力学模型。
在一种可选的方案中,高空气球温度微分方程,即热力学微分方程为:
其中,cf为蒙皮比热容,Mfilm为蒙皮质量;γ=cp/cv,cp,cv分别为氦气的定压比热容和定体积比热容,Mgas为氦气质量,Volume为囊体体积,Q表示各种热源。
可选地,在本发明上述实施例中,步骤S104,获取高空气球的当前状态和风场数据包括:
步骤S1042,通过惯性导航设备和测控设备得到当前状态。
在一种可选的方案中,可以由高空气球上安装的惯性导航设备和测控设备得到位置、速度信息,以及高空气球的氦气量,温度,压强等当前状态。
步骤S1044,通过探空气球得到风场数据。
在一种可选的方案中,风速风向等风场数据可以由放飞探空气球获得或者由气象部门获取。
可选地,在本发明上述实施例中,航迹预测模型还包括:放气模型,其中,步骤S102,建立高空气球的航迹预测模型包括:
步骤S1026,根据高空气球的压力信息,得到放气模型,其中,压力信息至少包括:高空气球的内部压力和高空气球的外部压力,放气模型至少包括:放气速度和放气质量。
在一种可选的方案中,在高空气球放气时,内部氦气的气压要大于外部空气的气压,根据伯努利原理可知,有下式成立:
其中,ρhe为氦气的密度,Phe为球内氦气的压强,P为球外空气的压强,v为放气速度,由此可得到的表达式如下所示:
其中,△P为压强差,即,△P=Phe-P,则氦气质量的微分方程如下式所示:
其中,S′为放气阀门面积,m为氦气的放气质量。
此处需要说明的是,高空气球的囊体由于加工等原因造成的囊体泄露也可以用等效泄露面积来计算氦气泄露量。
可选地,在本发明上述实施例中,在根步骤S106,据当前状态和风场数据,通过航迹预测模型,得到高空气球的航迹之后,该方法还包括:
步骤S108,根据当前状态,通过放气模型,控制高空气球进行放气。
在一种可选的方案中,高空气球在飞行过程中会产生超压量,超压量超过囊体的耐压性能后会破坏囊体,为了保证飞行安全性,需要安装放气阀门进行超压放气控制。
图2是根据本发明实施例的一种可选的高空气球航迹的预测方法的示意图,下面结合图2对本发明一种优选的实施例进行详细说明,如图2所示,动力学模型:动力学模型是航迹预测的基础,高空气球之后的飞行航迹需要由动力学模型获得。动力学模型可以采用单刚体动力学建模方法,可以将高空气球考虑成三自由度,不考虑姿态,也可以考虑六自由度刚体建模,受力分为重力,浮力,阻力,附加质量力,控制力。热力学模型:热力学模型包含太阳直射辐射、太阳反射、大气散射、地面红外辐射、蒙皮红外辐射、对流换热等热源,可以根据实际状况获取囊体光学特性参数以及热环境参数。当前状态:包括位置、速度信息,以及高空气球的氦气量,温度,压强等信息,可以由惯性导航设备和测控设备得到。预测风场数据:高空气球在飞行过程中,由于其体型巨大,受到风的作用影响巨大,航迹预测需要有准确的风速风向信息,风速风向等风场数据可以由放飞探空气球获得或者由气象部门获取。
此处需要说明的是,动力学模型和预测风场数据可以由不同的方式得到,其组成的任意方式都是本发明上述实施例的具体实施方式。
通过本发明上述实施例,提供一种高空气球航迹实时预测方法,结合真实飞行航迹与模型预测航迹的优点,利用当前真实飞行航迹作为模型预测的输入,基于高空气球动力学与热力学模型预测之后的航迹,真实航迹更新可以实时预测之后的航迹,同时,真实航迹的输入可以减小模型预测的误差,使得本技术方案同时拥有简便易用与精确预测的优点。上述方法不仅适用于无动力高空气球飞行过程中的航迹预测,同时适用于无动力高空气球在地面的落点预测与飞行过程的轨迹预测,一般平流层飞行器动力小,飞行主要受风场影响,也可应用于小动力高空气球航迹预测。
实施例2
根据本发明实施例,提供了一种高空气球航迹的预测装置实施例。
图3是根据本发明实施例的一种高空气球航迹的预测装置的示意图,如图3所示,该装置包括:
建立单元31,用于建立高空气球的航迹预测模型,其中,航迹预测模型至少包括:动力学模型和热力学模型。
具体的,高空气球可以是平流层高空气球,由一个浮升囊体和挂接在囊体下端的吊舱组成,囊体内充入氦气,以高空气球所需要的浮力。
在一种可选的方案中,动力学模型是航迹预测的基础,高空气球的飞行航迹需要由动力学模型获得。高空气球是一个气体浮空器,其浮力大小受热环境影响,氦气在球体中满足理想气体状态方程为:PV=nRT,其中,P为氦气的压强,V为氦气的体积,即球体体积,n为氦气的物质的量,R气体常数,T为氦气的热力学温度。在上升过程中,囊体体积会随温度变化而变化,在囊体到达最大体积后,温度继续增大就会产生超压量,所以需要建立高空飞行过程中的热力学模型,得到囊体内的氦气的温度变化。因此,在对高空气球进行航迹预测之前,需要建立高空气球的动力学模型和热力学模型。
获取单元33,用于获取高空气球的当前状态和风场数据,其中,当前状态至少包括:位置信息、速度信息、氦气量、温度和压强,风场数据至少包括:风向数据和风速数据。
处理单元35,用于根据当前状态和风场数据,通过航迹预测模型,得到高空气球的航迹。
在一种可选的方案中,高空气球在飞行过程中,由于其体型巨大,受到风的作用影响巨大,航迹预测需要有准确的风速风向信息,在建立高空气球的动力学模型和热力学模型之后,可以获取高空气球的当前时刻的真实状态和风场数据,通过动力学模型与热力学模型预测之后时间的飞行轨迹即航迹。
根据本发明上述实施例,建立单元建立高空气球的航迹预测模型,获取单元获取高空气球的当前状态和风场数据,处理单元根据当前状态和风场数据,通过航迹预测模型,得到高空气球的航迹。容易注意到的是,通过航迹预测模型预测高空气球的航迹,可以简化预测过程,将根据当前状态和风场数据作为航迹预测模型的输入,通过高空气球的真实飞行航迹更新实时预测后的航迹,可以减小航迹预测模型预测的误差,解决了现有技术中的高空气球航迹的预测方法无法同时满足简便预测和高精确度的要求的技术问题。因此,通过本发明上述实施例提供的方案,可以达到简便易用,精确预测的效果。
可选地,在本发明上述实施例中,建立单元至少包括:
第一处理模块,用于根据高空气球的动力学信息,得到动力学模型,其中,动力学信息至少包括:高空气球的质量、高空气球的参考面积和高空气球的速度。
具体的,上述的高空气球的质量可以包括整个高空气球的结构材料的质量,囊体内充入的氦气质量和附加空气质量,上述的高空气球的速度可以包括高空气球在垂直方向速度和水平方向速度。
在一种可选的方案中,可以采用单刚体动力学建模方法,可以将气球考虑成三自由度,不考虑姿态,也可以考虑六自由度刚体建模,受力分为重力,浮力,阻力,附加质量力,控制力。可以将东、北、天三个方向分别定义为地面惯性坐标系的x,y,z正向,即向东为x正向,向北为y正向,向上为z正向。对于高空气球,飞行过程中,受到的作用力主要有浮力,重力,附加质量力,气动力,由于高空气球的姿态变化对于受力影响较小,可以忽略球的姿态变化而将其看成质点,故竖直方向其质点动力学模型如下:
其中,M为整个系统的质量,包括结构材料质量,氦气质量,附加空气质量,充放气过程中h为高空气球的高度,Fhe为高空气球所受的浮力,Fdz为高空气球垂直方向所受的空气阻力,G为高空气球所受的重力。
其竖直速度与高度的关系如下式所示:
空气阻力大小与空气密度ρ,垂直方向的阻力系数Cdz,囊体参考面积S,高空气球垂直速度vz(忽略垂直风)有关,Fdz如下式所示:
重力与气球重量M,当地重力加速度g有关,重力加速度随高度变化,G如下式所示:
G=Mg。
浮力与空气密度ρ,囊体体积V,重力加速度g有关,Fhe如下式所示:
Fhe=ρVg。
水平面受到的力是风引起的空气阻力,动力学方程如下所示:
其中,Cdx和Cdy分别为x和y方向的阻力系数,vx和vy分别为高空气球在x和y方向的速度,vxwind和vywind分别为风在x和y方向的速度。
水平方向的运动学方程如下所示:
第二处理模块,用于根据高空气球的热环境信息,得到热力学模型,其中,热环境信息至少包括:囊体光学特性参数和热环境参数。
具体的,上述的囊体光学特性参数可以包括:囊皮材料对太阳光的吸收率和透射率,高空气球的照射投影面积,囊体的外露面积,囊体蒙皮材料对红外辐射的吸收率和透射率,蒙皮的红外发射率,蒙皮平均温度,外部和内部对流换热系数,有效对流面积,氦气的平均温度,蒙皮比热容,蒙皮质量,氦气的定压比热容和定体积比热容,囊体体积;上述的热环境参数可以包括:太阳辐射强度,大气对太阳光的透射率,太阳高度角,地面平均反照率,高空气球表面对地球表面的角系数,地面平均发射率,大气对地面红外辐射的透射率,地面温度,大气红外辐射强度,空气的平均温度。
在一种可选的方案中,热力学模型包含太阳直射辐射、太阳反射、大气散射、地面红外辐射、蒙皮红外辐射、对流换热等热源,可以根据实际状况获取囊体光学特性参数以及热环境参数,并根据获取到的热环境信息得到热力学模型。
可选地,在本发明上述实施例中,第二处理模块包括:
第一处理子模块,用于根据高空气球的热环境信息,得到高空气球的热源模型,其中,热源至少包括:太阳直射辐射、太阳反射、大气散射、地面红外辐射、蒙皮红外辐射、大气红外辐射和对流换热。
在一种可选的方案中,热源主要有以下几种:
(1)直接太阳辐射,对直接太阳辐射的吸收包括外表面的吸收和透过蒙皮由内表面的吸收,而内表面由于多次反射的作用会提高内部的有效吸收率。假设囊皮材料对太阳光的吸收率与透射率分别为α和τ,反射率为r=1-α-τ,有效反射率reffective=r+r2+r3+…,容易得到对太阳直接辐射的吸收为:
其中,Aproject为对高空气球的照射投影面积,由几何模型获得;ISun为太阳辐射强度;τatm是大气对太阳光的透射率。
(2)大气散射,对于大气散射,可以由下式得到:
其中,κ为大气散射经验系数,通常取0.1;Asurf为囊体的外露面积。
(3)地球反照辐射,对地面反照辐射吸收为:
其中,qAlbedo=AlbedoISun sin(ELV);ELV为太阳高度角,特别地,太阳高度角小于0时,qAlbedo=0;Albedo为地面平均反照率,与经度和当地气象条件有关,一般取0.3;ViewFactor为高空气球表面对地球表面的角系数,与高度有关。
(4)地面红外辐射,地面红外辐射强度可以由斯特潘-波尔茨曼定理计算:
其中,αIR和τIR分别为囊体蒙皮材料对红外辐射的吸收率和透射率;εground为地面平均发射率,取0.95,σ为斯特潘-波尔茨曼常数,τatmIR为大气对地面红外辐射的透射率,Tground为地面温度,蒙皮材料对红外辐射的反射率rIR=1-αIR-τIR。
(5)蒙皮的红外辐射,囊体蒙皮的内外表面均向外辐射热量,内表面辐射的热量会有部分被蒙皮吸收,故蒙皮的红外辐射包括蒙皮对外的红外辐射以及蒙皮吸收的红外辐射,可以由下式得到蒙皮对外的红外辐射QIRFilm和蒙皮吸收的红外辐射QIRout:
其中,ε为蒙皮的红外发射率,假设满足基尔霍夫定理,红外发射率等于红外吸收率;Tfilm为蒙皮平均温度。
(6)大气红外辐射,对于大气红外辐射,可以由下式得到:
其中,大气等效温度IIRsky为大气红外辐射强度。
(7)对流换热,对流换热包括外表面与大气环境的对流换热以及内表面与氦气间的对流换热,蒙皮的对流换热为:
QconvExt=HCexternalAeffective(Tair-Tfilm),
QconvInt=HCinternalAeffective(Tfilm-Tgas),
其中,HCexternal和HCinternal分别为外部和内部对流换热系数,可由经验公式获得;Aeffective为有效对流面积,由几何模型获得;Tair和Tgas分别为空气和氦气的平均温度。
第二处理子模块,用于根据高空气球的热源模型,得到热力学模型。
在一种可选的方案中,高空气球温度微分方程,即热力学微分方程为:
其中,cf为蒙皮比热容,Mfilm为蒙皮质量;γ=cp/cv,cp,cv分别为氦气的定压比热容和定体积比热容,Mgas为氦气质量,Volume为囊体体积,Q表示各种热源。
可选地,在本发明上述实施例中,获取单元包括:
第三处理模块,用于通过惯性导航设备和测控设备得到当前状态。
在一种可选的方案中,可以由高空气球上安装的惯性导航设备和测控设备得到位置、速度信息,以及高空气球的氦气量,温度,压强等当前状态。
第四处理模块,用于通过探空气球得到风场数据。
在一种可选的方案中,风速风向等风场数据可以由放飞探空气球获得或者由气象部门获取。
可选地,在本发明上述实施例中,航迹预测模型还包括:放气模型,其中,建立单元包括:
建立模块,用于根据高空气球的压力信息,得到放气模型,其中,压力信息至少包括:高空气球的内部压力和高空气球的外部压力,放气模型至少包括:放气速度和放气质量。
在一种可选的方案中,在高空气球放气时,内部氦气的气压要大于外部空气的气压,根据伯努利原理可知,有下式成立:
其中,ρhe为氦气的密度,Phe为球内氦气的压强,P为球外空气的压强,v为放气速度,由此可得到的表达式如下所示:
其中,△P为压强差,即,△P=Phe-P,则氦气质量的微分方程如下式所示:
其中,S′为放气阀门面积,m为氦气的放气质量。
此处需要说明的是,高空气球的囊体由于加工等原因造成的囊体泄露也可以用等效泄露面积来计算氦气泄露量。
可选地,在本发明上述实施例中,该装置还包括:
控制单元,用于根据当前状态,通过放气模型,控制高空气球进行放气。
在一种可选的方案中,高空气球在飞行过程中会产生超压量,超压量超过囊体的耐压性能后会破坏囊体,为了保证飞行安全性,需要安装放气阀门进行超压放气控制。
上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
在本发明的上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其他实施例的相关描述。
在本申请所提供的几个实施例中,应该理解到,所揭露的技术内容,可通过其它的方式实现。其中,以上所描述的装置实施例仅仅是示意性的,例如所述单元的划分,可以为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,单元或模块的间接耦合或通信连接,可以是电性或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (12)
1.一种高空气球航迹的预测方法,其特征在于,包括:
建立高空气球的航迹预测模型,其中,所述航迹预测模型至少包括:动力学模型和热力学模型;
获取所述高空气球的当前状态和风场数据,其中,所述当前状态至少包括:位置信息、速度信息、氦气量、温度和压强,所述风场数据至少包括:风向数据和风速数据;
根据所述当前状态和风场数据,通过所述航迹预测模型,得到所述高空气球的航迹。
2.根据权利要求1所述的方法,其特征在于,建立高空气球的航迹预测模型至少包括:
根据所述高空气球的动力学信息,得到所述动力学模型,其中,所述动力学信息至少包括:所述高空气球的质量、所述高空气球的参考面积和所述高空气球的速度;
根据所述高空气球的热环境信息,得到所述热力学模型,其中,所述热环境信息至少包括:囊体光学特性参数和热环境参数。
3.根据权利要求2所述的方法,其特征在于,根据所述高空气球的热环境信息,得到所述热力学模型包括:
根据所述高空气球的热环境信息,得到所述高空气球的热源模型,其中,所述热源至少包括:太阳直射辐射、太阳反射、大气散射、地面红外辐射、蒙皮红外辐射、大气红外辐射和对流换热;
根据所述高空气球的热源模型,得到所述热力学模型。
4.根据权利要求1所述的方法,其特征在于,获取所述高空气球的当前状态和风场数据包括:
通过惯性导航设备和测控设备得到所述当前状态;
通过探空气球得到所述风场数据。
5.根据权利要求1至4中任意一项所述的方法,其特征在于,所述航迹预测模型还包括:放气模型,其中,建立高空气球的航迹预测模型包括:
根据所述高空气球的压力信息,得到所述放气模型,其中,所述压力信息至少包括:所述高空气球的内部压力和所述高空气球的外部压力,所述放气模型至少包括:放气速度和放气质量。
6.根据权利要求5所述的方法,其特征在于,在根据所述当前状态和风场数据,通过所述航迹预测模型,得到所述高空气球的航迹之后,所述方法还包括:
根据所述当前状态,通过所述放气模型,控制所述高空气球进行放气。
7.一种高空气球航迹的预测装置,其特征在于,包括:
建立单元,用于建立高空气球的航迹预测模型,其中,所述航迹预测模型至少包括:动力学模型和热力学模型;
获取单元,用于获取所述高空气球的当前状态和风场数据,其中,所述当前状态至少包括:位置信息、速度信息、氦气量、温度和压强,所述风场数据至少包括:风向数据和风速数据;
处理单元,用于根据所述当前状态和风场数据,通过所述航迹预测模型,得到所述高空气球的航迹。
8.根据权利要求7所述的装置,其特征在于,所述建立单元至少包括:
第一处理模块,用于采用单刚体动力学建模装置,得到所述动力学模型;
第二处理模块,用于根据所述高空气球的热环境信息,得到所述热力学模型,其中,所述热环境信息至少包括:囊体光学特性参数和热环境参数。
9.根据权利要求8所述的装置,其特征在于,所述第二处理模块包括:
第一处理子模块,用于根据所述高空气球的热环境信息,得到所述高空气球的热源模型,其中,所述热源至少包括:太阳直射辐射、太阳反射、大气散射、地面红外辐射、蒙皮红外辐射、大气红外辐射和对流换热;
第二处理子模块,用于根据所述高空气球的热源模型,得到所述热力学模型。
10.根据权利要求7所述的装置,其特征在于,所述获取单元包括:
第三处理模块,用于通过惯性导航设备和测控设备得到所述当前状态;
第四处理模块,用于通过探空气球得到所述风场数据。
11.根据权利要求7至10中任意一项所述的装置,其特征在于,所述航迹预测模型还包括:放气模型,其中,所述建立单元包括:
建立模块,用于根据伯努利原理,建立所述放气模型,其中,所述放气模型至少包括:放气速度和放气质量。
12.根据权利要求11所述的装置,其特征在于,所述装置还包括:
控制单元,用于根据所述当前状态,通过所述放气模型,控制所述高空气球进行放气。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611229475.3A CN108241769A (zh) | 2016-12-27 | 2016-12-27 | 高空气球航迹的预测方法及装置 |
PCT/CN2017/092011 WO2018120736A1 (zh) | 2016-12-27 | 2017-07-06 | 高空气球航迹的预测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611229475.3A CN108241769A (zh) | 2016-12-27 | 2016-12-27 | 高空气球航迹的预测方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108241769A true CN108241769A (zh) | 2018-07-03 |
Family
ID=62702797
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611229475.3A Pending CN108241769A (zh) | 2016-12-27 | 2016-12-27 | 高空气球航迹的预测方法及装置 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN108241769A (zh) |
WO (1) | WO2018120736A1 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110889256A (zh) * | 2019-12-04 | 2020-03-17 | 中国特种飞行器研究所 | 一种中低空气球飞行航迹预测方法 |
CN111189656A (zh) * | 2019-11-01 | 2020-05-22 | 中国科学院光电研究院 | 一种浮空器搭载行星采样返回舱回收试验系统 |
CN112395686A (zh) * | 2020-11-16 | 2021-02-23 | 北京临近空间飞艇技术开发有限公司 | 基于气象探测的平流层浮空器飞行轨迹预测系统及方法 |
CN112461125A (zh) * | 2020-10-29 | 2021-03-09 | 北京空间机电研究所 | 一种密闭浮空气球位置姿态光学测量系统及测量方法 |
CN117421984A (zh) * | 2023-11-02 | 2024-01-19 | 中国人民解放军海军航空大学 | 基于多模态学习的短期航迹预测方法 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109977609B (zh) * | 2019-04-16 | 2022-08-23 | 哈尔滨工业大学 | 一种基于真实遥感数据的地面高温热源红外图像仿真方法 |
CN113592127B (zh) * | 2020-04-30 | 2024-02-02 | 国家气象中心 | 一种往返式探空系统轨迹预测方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201870A (zh) * | 2007-12-20 | 2008-06-18 | 四川川大智胜软件股份有限公司 | 空中交通飞行态势动态仿真的方法 |
CN104950908A (zh) * | 2015-07-02 | 2015-09-30 | 上海交通大学 | 平流层飞艇水平位置控制系统及实现方法 |
CN105836097A (zh) * | 2016-03-30 | 2016-08-10 | 北京航空航天大学 | 一种高度与轨迹可控的高空气球 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102902204B (zh) * | 2012-10-24 | 2016-06-01 | 北京控制工程研究所 | 一种飞艇遥控起飞/着陆预警方法 |
CN103472504B (zh) * | 2013-09-26 | 2015-11-11 | 中国人民解放军理工大学 | 船用高空气象探测放球装置及其制作方法 |
CN104007665A (zh) * | 2014-05-30 | 2014-08-27 | 北京航空航天大学 | 一种固液动力飞行器飞行仿真测试系统 |
-
2016
- 2016-12-27 CN CN201611229475.3A patent/CN108241769A/zh active Pending
-
2017
- 2017-07-06 WO PCT/CN2017/092011 patent/WO2018120736A1/zh active Application Filing
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201870A (zh) * | 2007-12-20 | 2008-06-18 | 四川川大智胜软件股份有限公司 | 空中交通飞行态势动态仿真的方法 |
CN104950908A (zh) * | 2015-07-02 | 2015-09-30 | 上海交通大学 | 平流层飞艇水平位置控制系统及实现方法 |
CN105836097A (zh) * | 2016-03-30 | 2016-08-10 | 北京航空航天大学 | 一种高度与轨迹可控的高空气球 |
Non-Patent Citations (2)
Title |
---|
FARLEY RODGER E.: "BalloonAscent: 3-D Simulation Tool for the Ascent and Float of High-Altitude Balloons", 《DOI:10.2514/6.2005-7412》 * |
吕明云,巫资春: "高空气球热力学模型与上升过程仿真分析", 《北京航空航天大学学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111189656A (zh) * | 2019-11-01 | 2020-05-22 | 中国科学院光电研究院 | 一种浮空器搭载行星采样返回舱回收试验系统 |
CN111189656B (zh) * | 2019-11-01 | 2022-03-22 | 中国科学院空天信息创新研究院 | 一种浮空器搭载行星采样返回舱回收试验系统 |
CN110889256A (zh) * | 2019-12-04 | 2020-03-17 | 中国特种飞行器研究所 | 一种中低空气球飞行航迹预测方法 |
CN110889256B (zh) * | 2019-12-04 | 2023-12-22 | 中国特种飞行器研究所 | 一种中低空气球飞行航迹预测方法 |
CN112461125A (zh) * | 2020-10-29 | 2021-03-09 | 北京空间机电研究所 | 一种密闭浮空气球位置姿态光学测量系统及测量方法 |
CN112461125B (zh) * | 2020-10-29 | 2023-02-28 | 北京空间机电研究所 | 一种密闭浮空气球位置姿态光学测量系统及测量方法 |
CN112395686A (zh) * | 2020-11-16 | 2021-02-23 | 北京临近空间飞艇技术开发有限公司 | 基于气象探测的平流层浮空器飞行轨迹预测系统及方法 |
CN117421984A (zh) * | 2023-11-02 | 2024-01-19 | 中国人民解放军海军航空大学 | 基于多模态学习的短期航迹预测方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2018120736A1 (zh) | 2018-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108241769A (zh) | 高空气球航迹的预测方法及装置 | |
Dai et al. | Performance simulation of high altitude scientific balloons | |
Nebbali et al. | Dynamic simulation of the distributed radiative and convective climate within a cropped greenhouse | |
Farley | BalloonAscent: 3-D simulation tool for the ascent and float of high-altitude balloons | |
SA519402008B1 (ar) | مركبة هوائية غير مزودة بأفراد بدرجة تحمل عالية | |
Allen | Updraft model for development of autonomous soaring uninhabited air vehicles | |
Wu et al. | Thermal modeling of stratospheric airships | |
Kaushik | Theoretical and experimental aerodynamics | |
Liu et al. | A comprehensive numerical model investigating the thermal-dynamic performance of scientific balloon | |
Zhang et al. | Thermal performance analysis of a high-altitude solar-powered hybrid airship | |
Zhang et al. | Influences of initial launch conditions on flight performance of high altitude balloon ascending process | |
CN110325442A (zh) | 防风空中投放方法和系统 | |
Schuler et al. | Long duration flights in Venus’ atmosphere using passive solar hot air balloons | |
Hall et al. | Prototype development of a variable altitude venus aerobot | |
Alam et al. | High altitude airship: A review of thermal analyses and design approaches | |
Alam et al. | A multi-node model for transient heat transfer analysis of stratospheric airships | |
Palumbo et al. | ACHAB: Analysis code for high-altitude balloons | |
Jiang et al. | Performance evaluation for scientific balloon station-keeping strategies considering energy management strategy | |
CN116576862A (zh) | 一种临近空间飞艇航迹规划方法 | |
Schuler et al. | Altitude control of a solar balloon for mars exploration | |
Jiang et al. | An optimization approach for improving the solar array output power of stratospheric aerostat | |
Zhang et al. | Effect of vapor condensation on ascending performance of stratospheric airship | |
Gai et al. | A scientific mission based on a high altitude stratospheric balloon | |
Yoder et al. | Aerodynamic Tethered Sails for Scientific Balloon Trajectory Control: Small-Scale Experimental Demonstration | |
Yang et al. | Thermal Analysis of a Stratospheric Airship in Working Process Thermal Analysis of a Stratospheric Airship in Working Process |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180703 |