CN105426606B - 带太阳能电池的平流层飞艇平飞过程平均温度计算方法 - Google Patents
带太阳能电池的平流层飞艇平飞过程平均温度计算方法 Download PDFInfo
- Publication number
- CN105426606B CN105426606B CN201510772490.1A CN201510772490A CN105426606B CN 105426606 B CN105426606 B CN 105426606B CN 201510772490 A CN201510772490 A CN 201510772490A CN 105426606 B CN105426606 B CN 105426606B
- Authority
- CN
- China
- Prior art keywords
- airship
- heat
- solar cell
- half part
- radiation
- 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.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 20
- 238000012546 transfer Methods 0.000 claims abstract description 29
- 239000000463 material Substances 0.000 claims abstract description 27
- 238000013461 design Methods 0.000 claims abstract description 21
- 210000004027 cell Anatomy 0.000 claims description 106
- 230000005855 radiation Effects 0.000 claims description 77
- 239000001307 helium Substances 0.000 claims description 54
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 claims description 54
- 229910052734 helium Inorganic materials 0.000 claims description 53
- 238000000034 method Methods 0.000 claims description 32
- 230000008569 process Effects 0.000 claims description 18
- 239000012774 insulation material Substances 0.000 claims description 12
- 238000009413 insulation Methods 0.000 claims description 9
- 239000004576 sand Substances 0.000 claims description 2
- 238000010521 absorption reaction Methods 0.000 claims 1
- 230000007613 environmental effect Effects 0.000 abstract description 4
- 238000012360 testing method Methods 0.000 abstract description 4
- 230000001052 transient effect Effects 0.000 description 13
- 230000008859 change Effects 0.000 description 3
- 238000004971 IR microspectroscopy Methods 0.000 description 2
- 238000012060 immune response imaging Methods 0.000 description 2
- 230000002146 bilateral effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 239000011810 insulating material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000035882 stress Effects 0.000 description 1
- 230000008646 thermal stress Effects 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Photovoltaic Devices (AREA)
Abstract
本发明提供了一种带太阳能电池的平流层飞艇平飞过程平均温度计算方法,其根据飞艇飞行参数、飞艇设计参数、艇体材料特性参数、太阳能电池特性参数及电池隔热材料特性参数,计算大气环境参数及飞艇热环境参数,并基于飞艇几何特征及传热模式,将飞艇划分为多个节点,建立各节点的能量微分方程,通过求解飞艇多节点的能量微分方程组,计算飞艇平飞过程各节点平均温度数据。本发明在带太阳能电池的平流层飞艇设计、材料选择、飞行试验规划、规避潜在危险等方面具有指导意义,可以提高带太阳能电池的平流层飞艇设计一次成功率,缩短带太阳能电池的平流层飞艇设计周期,降低带太阳能电池的平流层飞艇设计成本。
Description
技术领域
本发明属于飞艇热控制技术领域,尤其涉及一种带太阳能电池的平流层飞艇平飞过程平均温度计算方法。
背景技术
平流层飞艇具有可定点飞行、滞空时间长和分辨率高等优点,在空中预警、监视监测、民用通信等领域具有广泛应用前景,受到世界各主要强国的高度重视。
平流层飞艇在平飞过程中,环境温度、密度、压力、风速、太阳辐射、大气辐射和地面辐射等因素会对飞艇温度特性产生影响。温度过高将提高飞艇内部氦气压力,对飞艇产生重要影响:1、温度过高将改变飞艇艇体材料承力特性、增大飞艇艇体热应力、增大飞艇艇体张力,对飞艇艇体的安全构成严重威胁;2、改变飞艇受力状况,导致飞艇飞行高度波动,干扰飞艇执行任务。因此,准确获知飞艇平飞过程中的温度特性,对飞艇结构设计、材料选择、飞行试验规划、规避潜在危险等方面具有重要意义,而目前还没有一个系统性地计算带太阳能电池的平流层飞艇平飞过程平均温度的计算方法。
发明内容
(一)要解决的技术问题
本发明的目的在于,提供一种带太阳能电池的平流层飞艇平飞过程平均温度计算方法,可快速而准确地获得带太阳能电池的平流层飞艇平飞过程平均温度数据。
(二)技术方案
本发明提供一种带太阳能电池的平流层飞艇平飞过程平均温度计算方法,包括:
S1,根据飞艇飞行任务需求,计算飞艇飞行参数及飞艇设计参数;
S2,测量艇体材料特性参数、太阳能电池特性参数及电池隔热材料特性参数;
S3,计算飞艇大气环境参数及飞艇热环境参数;
S4,基于飞艇几何特征及传热模式,将飞艇划分为多个节点,建立各节点的能量微分方程;
S5,根据艇体材料特性参数和太阳能电池特性参数,联立求解飞艇多节点的能量微分方程组,计算飞艇平飞过程各节点平均温度数据。
(三)有益效果
本发明可以快速和准确地获知带太阳能电池的平流层飞艇平飞过程中的平均温度特性,在带太阳能电池的平流层飞艇设计、材料选择、飞行试验规划、规避潜在危险等方面具有指导意义,可以提高带太阳能电池的平流层飞艇设计一次成功率,缩短带太阳能电池的平流层飞艇设计周期,降低带太阳能电池的平流层飞艇设计成本。
附图说明
图1是本发明实施例提供的带太阳能电池的平流层飞艇结构示意图。
图2是本发明实施例提供的带太阳能电池的平流层飞艇平飞过程平均温度计算方法流程图。
具体实施方式
本发明提供一种带太阳能电池的平流层飞艇平飞过程平均温度计算方法,其根据飞艇飞行参数、飞艇设计参数、艇体材料特性参数、太阳能电池特性参数及电池隔热材料特性参数,计算大气环境参数及飞艇热环境参数,并基于飞艇几何特征及传热模式,将飞艇划分为多个节点,建立各节点的能量微分方程,通过求解飞艇多节点的能量微分方程组,计算飞艇平飞过程各节点平均温度数据。
根据本发明的一种实施方式,温度计算方法包括:
S1,根据飞艇飞行任务需求,计算飞艇飞行参数及飞艇设计参数;
S2,测量艇体材料特性参数、太阳能电池特性参数及电池隔热材料特性参数;
S3,计算飞艇大气环境参数及飞艇热环境参数;
S4,基于飞艇几何特征及传热模式,将飞艇划分为多个节点,建立各节点的能量微分方程;
S5,根据艇体材料特性参数和太阳能电池特性参数,联立求解飞艇多节点的能量微分方程组,计算飞艇平飞过程各节点平均温度数据。
根据本发明的一种实施方式,飞艇飞行参数包括飞艇飞行时间、飞艇飞行地点经度Lon、飞艇飞行地点纬度Lat、飞艇飞行海拔高度h和飞艇飞行空速v;
飞艇设计参数包括飞艇体积V、飞艇长度L、飞艇最大直径D、飞艇表面积A和太阳能电池面积AS。
根据本发明的一种实施方式,艇体材料特性参数包括艇体材料表面吸收率α、艇体材料表面发射率ε、艇体材料面密度ρ和艇体材料比热容c;
太阳能电池特性参数包括太阳能电池效率η、太阳能电池表面吸收率αS、太阳能电池表面发射率εS、太阳能电池面密度ρS和太阳能电池比热容cS;
电池隔热材料特性参数隔热材料特性参数包括隔热材料厚度δS_I和隔热材料导热系数λS_I。
根据本发明的一种实施方式,飞艇大气环境参数包括飞艇飞行海拔高度h处的大气温度TAtm、大气压力PAtm和大气密度ρAtm,
其中,大气温度TAtm的数学表达式为:
大气压力PAtm的数学表达式为:
大气密度ρAtm的数学表达式为:
飞艇热环境参数包括飞艇辐射热环境参数和对流换热环境参数,所述飞艇辐射热环境参数包括太阳直接辐射热流qD_S、大气散射太阳辐射热流qA_S、地面反射太阳辐射热流qG_S、大气长波辐射热流qA_IR和地面长波辐射热流qG_IR,
太阳直接辐射热流qD_S的数学表达式为:
qD_S=I0·τAtm,
其中,I0为大气层上界太阳辐射强度,τAtm为太阳直接辐射衰减系数;
所述大气散射太阳辐射热流qA_S的数学表达式为:
qA_S=k·qD_S,
其中,k为大气散射系数;
地面反射太阳辐射热流qG_S的数学表达式为:
qG_S=IGround·rGround·τIR_G,
其中,IGround为抵达地球表面太阳直接辐射强度,rGround为地球表面反射系数,τIR_G为地球表面辐射衰减系数;
所述大气长波辐射热流qA_IR的数学表达式为:
其中,σ为辐射常数,TAtm为大气温度;
地面长波辐射热流qG_IR的数学表达式为:
其中,TGround为地面温度,εGround为地面发射率;
对流换热环境参数包括飞艇与外部环境的对流换热系数hEx、飞艇与内部氦气的对流换热系数hIn,
飞艇与内部氦气的对流换热系数hIn的数学表达式为:
其中,NuEx为飞艇与外部空气的对流换热努赛尔数,λAir为空气导热系数;
飞艇与内部氦气的对流换热系数hIn的数学表达式为:
其中,NuIn为内部自然对流换热努赛尔数,λHe为氦气导热系数。
根据本发明的一种实施方式,多个节点包括太阳能电池、飞艇艇体上半部分被太阳能电池覆盖的部分、飞艇艇体上半部分未被太阳能电池覆盖的部分、飞艇艇体下半部分和飞艇内氦气,其中,
太阳能电池的能量微分方程为:
其中,TS为太阳能电池平均温度,t为时间,QS_D为吸收太阳直接辐射热量,QS_Atm为吸收大气散射辐射热量,QS_IR_Atm为吸收大气长波辐射热量,QS_IR为对外界环境长波辐射热量,QS_Conv为与外界环境对流换热热量,QS_Cond是通过隔热层与艇体上半部分的传导换热热量;
飞艇艇体上半部分被太阳能电池覆盖的部分的能量微分方程为:
其中,TEnup_S为艇体上半部分被太阳能电池遮盖的部分平均温度,AEnup_S为太阳能电池的面积,QEnup_S_IR为与艇体下半部分长波辐射换热热量,QEnup_S_ConvI为与飞艇内氦气对流换热热量,QEnup_S_Cond为通过隔热层与太阳能电池的传导换热热量;
飞艇艇体上半部分未被太阳能电池覆盖的部分的能量微分方程为:
其中,TEnup_R是艇体上半部分未被太阳能电池遮盖的部分平均温度,AEnup_R是艇体上半部分未被太阳能电池遮盖的部分面积,QER_D是吸收太阳直接辐射热量,QER_Atm是吸收大气散射辐射热量,QER_IR_Atm是吸收大气长波辐射热量,QER_IR_E是对外界环境长波辐射热量,QER_IR_I是与艇体下半部分长波辐射换热热量,QER_ConvE是与外界环境对流换热热量,QER_ConvI与飞艇内氦气对流换热热量;
飞艇艇体下半部分的能量微分方程为:
其中,TEndown是艇体下半部分平均温度,AEndown=A/2是艇体下半部分面积,QEnd_Atm是吸收大气散射辐射热量,QEnd_G是吸收地面反射辐射热量,QEnd_IR_Atm是吸收大气长波辐射热量,QEnd_IR_G是吸收地面长波辐射热量,QEnd_IR_E是对外界环境长波辐射热量,QEnd_IR_I是与艇体上半部分长波辐射换热热量,QEnd_ConvE是与外界环境对流换热热量,QEnd_ConvI是与飞艇内氦气对流换热热量;
飞艇内氦气的能量微分方程为:
其中,THe是飞艇内氦气平均温度,mHe是氦气质量,cp,He是氦气定压比热容,VHe是氦气体积,PHe是氦气压力。QHe_Enup_S是与艇体上半部分被太阳能电池遮盖的部分对流换热热量,QHe_Enup_R是与艇体上半部分未被太阳能电池遮盖的部分对流换热热量,QHe_Endown是与艇体下半部分对流换热热量。
根据本发明的一种实施方式,步骤S5中,利用四阶标准龙格-库塔法求解能量微分方程组。
综上所述,本发明可以快速和准确地获知带太阳能电池的平流层飞艇平飞过程中的平均温度特性,在带太阳能电池的平流层飞艇设计、材料选择、飞行试验规划、规避潜在危险等方面具有指导意义,可以提高带太阳能电池的平流层飞艇设计一次成功率,缩短带太阳能电池的平流层飞艇设计周期,降低带太阳能电池的平流层飞艇设计成本。
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
如图1所示,本发明实施例提供的带太阳能电池的平流层飞艇包括飞艇由艇体上半部分1、艇体下半部分2、太阳能电池3、太阳能电池隔热层4、尾翼5和推进装置6构成。
其中,飞艇主体由艇体上半部分1和艇体下半部分2构成,艇体上半部分顶部铺设有太阳能电池3,太阳能电池与艇体上半部分之间安装有隔热层4,尾翼5呈倒Y型安装于飞艇尾部,推进装置6左右对称安装于飞艇两侧。
如图2所示,带太阳能电池的平流层飞艇平飞过程平均温度计算方法,包括:
根据飞艇飞行任务需求,计算出本实施例中的飞艇主要飞行参数如表1所示,主要设计参数如表2所示。
表1飞艇主要飞行参数
表2飞艇主要设计参数
测量拟采用的飞艇艇体材料特性参数如表3所示;测量太阳能电池特性与太阳能电池隔热材料特性参数如表4所示。
表3艇体材料特性参数
表4太阳能电池与太阳能电池隔热材料特性参数
计算飞艇热环境:大气压力、温度、密度。其中,飞艇在海拔高度h处的大气温度TAtm(K)、大气压力PAtm(Pa)、大气密度ρAtm(kg/m3)可由公式计算:
大气温度随海拔高度h变化的数学表达式为:
大气压力随海拔高度h变化的数学表达式为:
大气密度随海拔高度h变化的数学表达式为:
计算太阳直接辐射热流qD_S,大气散射太阳辐射热流qA_S,地面反射太阳辐射热流qG_S,大气长波辐射热流qA_IR,地面长波辐射热流qG_IR;对流换热环境参数包括飞艇与外部环境的对流换热系数hEx,飞艇与内部氦气的对流换热系数hIn。
太阳直接辐射热流qD_S是大气层上界太阳辐射强度I0与太阳直接辐射衰减系数τAtm的乘积,计算式如下:
qD_S=I0·τAtm (4)
大气散射太阳辐射热流qA_S是太阳直接辐射热流qD_S与大气散射系数k的乘积,计算式如下:
qA_S=k·qD_S (5)
地面反射太阳辐射热流qG_S是抵达地球表面太阳直接辐射强度IGround、地球表面反射系数rGround与地球表面辐射衰减系数τIR_G的乘积,计算式如下:
qG_S=IGround·rGround·τIR_G (6)
大气长波辐射热流qA_IR计算式如下:
其中,σ是辐射常数,TAtm是大气温度。
地面长波辐射热流qG_IR计算式如下:
其中,TGround是地面温度。
对流换热环境参数包括飞艇与外部环境的对流换热系数,飞艇内部对流换热系数。
飞艇与外部环境的对流换热系数hEx的计算式:
其中,D是飞艇艇体与外界强迫对流换热特征长度,取飞艇最大直径。
其中,NuEx计算式为:
式中,NuEf是飞艇与外界强迫对流努赛尔数,NuEn是飞艇与外界自然对流努赛尔数。
飞艇与外界强迫对流换热努赛尔数NuEf的表达式为:
其中,Re是雷诺数,其计算式为:
飞艇与外界自然对流换热努赛尔数NuEn的表达式为:
式中,Ra是自然对流格拉晓夫数Gr与外界空气普朗特数PrAtm的乘积,其计算式为:
飞艇与内部氦气的对流换热系数hIn计算式为:
其中,NuIn为内部自然对流换热努赛尔数,其表达式为:
式中,Ra是自然对流格拉晓夫数Gr与内部氦气普朗特数PrHe的乘积,其计算式为:
建立飞艇各节点的瞬态能量方程,包括:太阳能电池、艇体上半部分被太阳能电池遮盖的部分、艇体上半部分未被太阳能电池遮盖的部分、艇体下半部分、内部氦气的瞬态能量方程。
太阳能电池瞬态能量方程表述如下:
其中,TS是太阳能电池平均温度,QS_D是吸收太阳直接辐射热量,QS_Atm是吸收大气散射辐射热量,QS_IR_Atm是吸收大气长波辐射热量,QS_IR是对外界环境长波辐射热量,QS_Conv是与外界环境对流换热热量,QS_Cond是通过隔热层与艇体上半部分的传导换热热量。
太阳能电池瞬态能量方程中各项热量计算式列述如下:
QS_D=αS·qD_S·AS·FS-S (19)
其中,FS-S是太阳能电池与太阳直接辐射的辐射角系数。
QS_Atm=αS·qIR_Atm·AS (20)
QS_IR_Atm=εS·qIR_Atm·AS (21)
QS_Conv=hEm·(TAtm-TS)·AS (23)
其中,hEm是太阳能电池与外界环境的对流换热系数,TAtm是外界环境温度。
通过隔热层与艇体上半部分的传导换热热量
其中,TEnup_S是艇体上半部分被太阳能电池遮盖的部分平均温度。
艇体上半部分被太阳能电池遮盖的部分瞬态能量方程表述如下:
其中,AEnup_S是艇体上半部分被太阳能电池遮盖的部分的面积,等于太阳能电池的面积。QEnup_S_IR是与艇体下半部分长波辐射换热热量,QEnup_S_ConvI是与飞艇内氦气对流换热热量,QEnup_S_Cond是通过隔热层与太阳能电池的传导换热热量。
艇体上半部分被太阳能电池遮盖的部分瞬态能量方程中各项热量计算式列述如下:
其中,σ是辐射常数,TEndown是艇体下半部分平均温度,FEnup_S-Endown是与艇体下半部分的辐射角系数。
QEnup_S_ConvI=hIn·(THe-TEnup_S)·AEnup_S (27)
其中,hIn是与内部氦气的对流换热系数。
艇体上半部分未被太阳能电池遮盖的部分瞬态能量方程表述如下:
其中,TEnup_R是艇体上半部分未被太阳能电池遮盖的部分平均温度,AEnup_R=A/2-AEnup_S是艇体上半部分未被太阳能电池遮盖的部分面积。QER_D是吸收太阳直接辐射热量,QER_Atm是吸收大气散射辐射热量,QER_IR_Atm是吸收大气长波辐射热量,QER_IR_E是对外界环境长波辐射热量,QER_IR_I是与艇体下半部分长波辐射换热热量,QER_ConvE是与外界环境对流换热热量,QER_ConvI与飞艇内氦气对流换热热量。
艇体上半部分未被太阳能电池遮盖的部分瞬态能量方程中各项热量计算式列述如下:
QER_D=α·qD_S·AEnup_R·FEnup_R-S (30)
其中,AEnup_R是艇体上半部分未被太阳能电池遮盖的部分面积,FEnup_R-S是艇体上半部分未被太阳能电池遮盖的部分与太阳直接辐射的辐射角系数。
QER_Atm=α·qA_S·AEnup_R (31)
QER_IR_Atm=ε·qA_IR·AEnup_R (32)
其中,ε是艇体材料发射率。
QER_ConvE=hEx·(TAtm-TEnup_R)·AEnup_R (35)
QER_ConvI=hIn·(THe-TEnup_R)·AEnup_R (36)
其中,THe是飞艇内氦气的温度。
艇体下半部分瞬态能量方程表述如下:
其中,TEndown是艇体下半部分平均温度,AEndown=A/2是艇体下半部分面积。QEnd_Atm是吸收大气散射辐射热量,QEnd_G是吸收地面反射辐射热量,QEnd_IR_Atm是吸收大气长波辐射热量,QEnd_IR_G是吸收地面长波辐射热量,QEnd_IR_E是对外界环境长波辐射热量,QEnd_IR_I是与艇体上半部分长波辐射换热热量,QEnd_ConvE是与外界环境对流换热热量,QEnd_ConvI是与飞艇内氦气对流换热热量。
艇体下半部分瞬态能量方程中各项热量计算式列述如下:
QEnd_Atm=α·qA_S·AEndown (38)
QEnd_G=α·qG_S·AEndown (39)
QEnd_IR_Atm=ε·qA_IR·AEndown (40)
QEnd_IR_G=ε·qG_IR·AEndown (41)
QEnd_ConvE=hEx·(TAtm-TEndown)·AEndown (44)
QEnd_ConvI=hIn·(THe-TEndown)·AEndown (45)
内部氦气瞬态能量方程表述如下:
其中,THe是飞艇内氦气平均温度,mHe是氦气质量,cp,He是氦气定压比热容,VHe是氦气体积,等于飞艇体积,PHe是氦气压力。QHe_Enup_S是与艇体上半部分被太阳能电池遮盖的部分对流换热热量,QHe_Enup_R是与艇体上半部分未被太阳能电池遮盖的部分对流换热热量,QHe_Endown是与艇体下半部分对流换热热量。
内部氦气瞬态能量方程中各项热量计算式列述如下:
QHe_Enup_S=hIn·(TEnup_S-THe)·AEnup_S (47)
QHe_Enup_R=hIn·(TEnup_R-THe)·AEnup_R (48)
QHe_Endown=hIn·(TEndown-THe)·AEndown (49)
氦气压力控制范围为:
0≤ΔPHe=PHe-PAtm≤300Pa (50)
其中,ΔPHe是氦气超压量,PHe是氦气绝对压力,PAtm是大气环境压力。
氦气质量控制方法为:当飞艇内部氦气超压超过300Pa时候,氦气阀门打开,排出部分氦气,至超压量等于300Pa时阀门关闭。
氦气质量流量计算式为:
其中,ρHe是氦气密度,Av_He是氦气阀门面积,kv_He是氦气阀门流量系数。
求解微分方程:
方程(18)、(25)、(28)、(36)、(25)、(45)和(50)含有一阶偏微分项,将其离散以便于编程计算。将方程(18)、(25)、(28)、(36)、(25)、(45)和(50)写成如下矢量形式:
y'=f(t,y) (52)
其中,
y=(TSTEnup_STEnup_RTEndownTHemHe)T (53)
f=(f1f2f3f4f5f6)T (54)
如果函数y具有连续的(n+1)阶导数,则函数(51)关于y的泰勒级数展开式为:
式中,Δt是时间步长,下标i表示当前时间节点,下标i+1表示下一个时间节点。将方程(54)转换成如下形式:
采用标准四阶龙格-库塔方法离散方程(55),得到
联立上述所有方程,开发仿真计算程序。
输入飞艇设计参数、飞行任务参数,计算得到飞艇各节点平均温度数据。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种带太阳能电池的平流层飞艇平飞过程平均温度计算方法,其特征在于,包括:
S1,根据飞艇飞行任务需求,计算飞艇飞行参数及飞艇设计参数;
S2,测量艇体材料特性参数、太阳能电池特性参数及电池隔热材料特性参数;
S3,计算飞艇大气环境参数及飞艇热环境参数;
S4,基于飞艇几何特征及传热模式,将飞艇划分为多个节点,建立各节点的能量微分方程,其中,所述多个节点包括太阳能电池、飞艇艇体上半部分被太阳能电池覆盖的部分、飞艇艇体上半部分未被太阳能电池覆盖的部分、飞艇艇体下半部分和飞艇内氦气;
S5,根据艇体材料特性参数和太阳能电池特性参数,联立求解飞艇多节点的能量微分方程组,计算飞艇平飞过程各节点平均温度数据。
2.根据权利要求1所述的温度计算方法,其特征在于,所述飞艇飞行参数包括飞艇飞行时间、飞艇飞行地点经度Lon、飞艇飞行地点纬度Lat、飞艇飞行海拔高度h和飞艇飞行空速v;
所述飞艇设计参数包括飞艇体积V、飞艇长度L、飞艇最大直径D、飞艇表面积A和太阳能电池面积AS。
3.根据权利要求2所述的温度计算方法,其特征在于,所述艇体材料特性参数包括艇体材料表面吸收率α、艇体材料表面发射率ε、艇体材料面密度ρ和艇体材料比热容c;
所述太阳能电池特性参数包括太阳能电池效率η、太阳能电池表面吸收率αS、太阳能电池表面发射率εS、太阳能电池面密度ρS和太阳能电池比热容cS;
所述电池隔热材料特性参数包括隔热材料厚度δS_I和隔热材料导热系数λS_I。
4.根据权利要求3所述的温度计算方法,其特征在于,所述飞艇大气环境参数包括飞艇飞行海拔高度h处的大气温度TAtm、大气压力PAtm和大气密度ρAtm,
其中,大气温度TAtm的数学表达式为:
大气压力PAtm的数学表达式为:
大气密度ρAtm的数学表达式为:
所述飞艇热环境参数包括飞艇辐射热环境参数和对流换热环境参数,所述飞艇辐射热环境参数包括太阳直接辐射热流qD_S、大气散射太阳辐射热流qA_S、地面反射太阳辐射热流qG_S、大气长波辐射热流qA_IR和地面长波辐射热流qG_IR,
所述太阳直接辐射热流qD_S的数学表达式为:
qD_S=I0·τAtm,
其中,I0为大气层上界太阳辐射强度,τAtm为太阳直接辐射衰减系数;
所述大气散射太阳辐射热流qA_S的数学表达式为:
qA_S=k·qD_S,
其中,k为大气散射系数;
所述地面反射太阳辐射热流qG_S的数学表达式为:
qG_S=IGround·rGround·τIR_G,
其中,IGround为抵达地球表面太阳直接辐射强度,rGround为地球表面反射系数,τIR_G为地球表面辐射衰减系数;
所述大气长波辐射热流qA_IR的数学表达式为:
其中,σ为辐射常数,TAtm为大气温度;
所述地面长波辐射热流qG_IR的数学表达式为:
其中,TGround为地面温度,εGround为地面发射率;
所述对流换热环境参数包括飞艇与外部环境的对流换热系数hEx、飞艇与内部氦气的对流换热系数hIn,
飞艇与外部环境的对流换热系数hEx的数学表达式为:
其中,NuEx为飞艇与外部空气的对流换热努赛尔数,λAir为空气导热系数;
飞艇与内部氦气的对流换热系数hIn的数学表达式为:
其中,NuIn为内部自然对流换热努赛尔数,λHe为氦气导热系数。
5.根据权利要求4所述的温度计算方法,其特征在于,所述太阳能电池的能量微分方程为:
其中,TS为太阳能电池平均温度,t为时间,QS_D为太阳能电池吸收太阳直接辐射热量,QS_Atm为太阳能电池吸收大气散射辐射热量,QS_IR_Atm为太阳能电池吸收大气长波辐射热量,QS_IR为太阳能电池对外界环境长波辐射热量,QS_Conv为太阳能电池与外界环境对流换热热量,QS_Cond为太阳能电池通过隔热层与艇体上半部分被太阳能电池覆盖的部分的传导换热热量;
飞艇艇体上半部分被太阳能电池覆盖的部分的能量微分方程为:
其中,TEnup_S为艇体上半部分被太阳能电池覆盖的部分平均温度,AEnup_S为太阳能电池的面积,QEnup_S_IR为与艇体上半部分被太阳能电池覆盖的部分与艇体下半部分长波辐射换热热量,QEnup_S_ConvI为艇体上半部分被太阳能电池覆盖的部分与飞艇内氦气对流换热热量,QEnup_S_Cond为艇体上半部分被太阳能电池覆盖的部分通过隔热层与太阳能电池的传导换热热量;
飞艇艇体上半部分未被太阳能电池覆盖的部分的能量微分方程为:
其中,TEnup_R是艇体上半部分未被太阳能电池覆盖的部分平均温度,AEnup_R是艇体上半部分未被太阳能电池覆盖的部分面积,QER_D是艇体上半部分未被太阳能电池覆盖的部分吸收太阳直接辐射热量,QER_Atm是艇体上半部分未被太阳能电池覆盖的部分吸收大气散射辐射热量,QER_IR_Atm是艇体上半部分未被太阳能电池覆盖的部分吸收大气长波辐射热量,QER_IR_E是艇体上半部分未被太阳能电池覆盖的部分对外界环境长波辐射热量,QER_IR_I是艇体上半部分未被太阳能电池覆盖的部分与艇体下半部分长波辐射换热热量,QER_ConvE是艇体上半部分未被太阳能电池覆盖的部分与外界环境对流换热热量,QER_ConvI与艇体上半部分未被太阳能电池覆盖的部分与飞艇内氦气对流换热热量;
飞艇艇体下半部分的能量微分方程为:
其中,TEndown是艇体下半部分平均温度,AEndown=A/2是艇体下半部分面积,QEnd_Atm是艇体下半部分吸收大气散射辐射热量,QEnd_G是艇体下半部分吸收地面反射辐射热量,QEnd_IR_Atm是艇体下半部分吸收大气长波辐射热量,QEnd_IR_G是艇体下半部分吸收地面长波辐射热量,QEnd_IR_E是艇体下半部分对外界环境长波辐射热量,QEnd_IR_I是艇体下半部分与艇体上半部分长波辐射换热热量,QEnd_ConvE是艇体下半部分与外界环境对流换热热量,QEnd_ConvI是艇体下半部分与飞艇内氦气对流换热热量;
飞艇内氦气的能量微分方程为:
其中,THe是飞艇内氦气平均温度,mHe是氦气质量,cp,He是氦气定压比热容,VHe是氦气体积,PHe是氦气压力;QHe_Enup_S是与艇体上半部分被太阳能电池遮盖的部分对流换热热量,QHe_Enup_R是与艇体上半部分未被太阳能电池遮盖的部分对流换热热量,QHe_Endown是与艇体下半部分对流换热热量。
6.根据权利要求5所述的温度计算方法,其特征在于,所述步骤S5中,利用四阶标准龙格-库塔法求解能量微分方程组。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510772490.1A CN105426606B (zh) | 2015-11-12 | 2015-11-12 | 带太阳能电池的平流层飞艇平飞过程平均温度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510772490.1A CN105426606B (zh) | 2015-11-12 | 2015-11-12 | 带太阳能电池的平流层飞艇平飞过程平均温度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105426606A CN105426606A (zh) | 2016-03-23 |
CN105426606B true CN105426606B (zh) | 2019-03-19 |
Family
ID=55504815
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510772490.1A Active CN105426606B (zh) | 2015-11-12 | 2015-11-12 | 带太阳能电池的平流层飞艇平飞过程平均温度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105426606B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106227060A (zh) * | 2016-07-19 | 2016-12-14 | 中国航天员科研训练中心 | 一种平流层飞艇载荷舱电子设备热特性模拟与仿真预测系统与方法 |
CN112487559B (zh) * | 2020-12-03 | 2023-04-28 | 中国人民解放军63660部队 | 参数自适应的单囊体平流层飞艇浮重平衡评估方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101976275A (zh) * | 2010-09-21 | 2011-02-16 | 北京航空航天大学 | 飞机红外辐射与大气透过率建模方法 |
-
2015
- 2015-11-12 CN CN201510772490.1A patent/CN105426606B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101976275A (zh) * | 2010-09-21 | 2011-02-16 | 北京航空航天大学 | 飞机红外辐射与大气透过率建模方法 |
Non-Patent Citations (7)
Title |
---|
《A comprehensive numerical model investigating the thermal-dynamic performance of scientific balloon》;Q. Liu等;《Advances in Space Research》;20131118;第53卷(第2期);第325-338页 * |
《A heat transient model for the thermal behavior prediction of stratospheric airships》;Wei Yao等;《Applied Thermal Engineering》;20140523;第70卷(第1期);第380-387页 * |
《Simplified analytical model for predicting the temperature of balloon on high-altitude》;Junjiang Xiong等;《International Journal of Thermal Sciences》;20140228;第76卷(第2期);第83-86页 * |
《临近空间飞艇外形优化设计与仿真》;杨燕初等;《中国空间科学技术》;20111025(第5期);第76-83页 * |
《平流层气球热动力学仿真》;刘强等;《北京航空航天大学学报》;20130911;第39卷(第12期);摘要,第1579-1583页第1-4节 * |
《平流层飞艇太阳能电池热特性研究》;李小建等;《第三届高分辨率对地观测学术年会优秀论文集》;20141201;第1-6页第1-3节 * |
《平流层飞艇热力学建模与仿真研究》;赵攀峰等;《合肥工业大学学报(自然科学版)》;20130428;第36卷(第4期);第501-505页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105426606A (zh) | 2016-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dai et al. | Performance simulation of high altitude scientific balloons | |
Zhang et al. | Thermal performance analysis of a high-altitude solar-powered hybrid airship | |
EP3845465B1 (en) | Aerostat icing characteristic numerical simulation and experimental verification system | |
CN109800488B (zh) | 关于液体火箭高空环境下底部热环境的数值计算方法 | |
Liu et al. | A comprehensive numerical model investigating the thermal-dynamic performance of scientific balloon | |
CN103616078B (zh) | 一种热红外载荷便携式野外定标及水面温度验证系统装置 | |
CN105426606B (zh) | 带太阳能电池的平流层飞艇平飞过程平均温度计算方法 | |
Zhu et al. | Thermal performance of high-altitude solar powered scientific balloon | |
CN108241769A (zh) | 高空气球航迹的预测方法及装置 | |
Pourbagian et al. | Surrogate-based optimization of electrothermal wing anti-icing systems | |
CN109145372A (zh) | 一种飞机油箱热计算建模方法及其计算模型 | |
CN105303053B (zh) | 带太阳能电池的高空气球平飞过程平均温度计算方法 | |
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 | |
Dai et al. | Numerical study of forced convective heat transfer around airships | |
CN117725794A (zh) | 一种柔性浮空器流热结构耦合计算方法 | |
Jiang et al. | An optimization approach for improving the solar array output power of stratospheric aerostat | |
CN105426604B (zh) | 带太阳能电池的平流层飞艇平飞过程分布温度计算方法 | |
Nosov et al. | Turbulence structure over heated surfaces: Numerical solutions | |
CN112560251B (zh) | 一种伪装网表面温度数值计算的简化方法 | |
Sun et al. | Membrane Material‐Based Rigid Solar Array Design and Thermal Simulation for Stratospheric Airships | |
CN105404774B (zh) | 带太阳能电池的高空气球平飞过程分布温度计算方法 | |
Mingalev et al. | A mechanism of formation of polar cyclones and possibility of their prediction using satellite observations | |
Yang et al. | Bionic Design and Experimental Validation of a Robotic Airship Inspired by the Physalia physalis | |
Liu et al. | Three dimensional simulation of stratospheric airship ice accretion in ascending process |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |