一种浮空器辐射热环境特性估算方法技术

本发明专利技术属于浮空器热环境分析领域，尤其涉及一种浮空器辐射热环境特性估算方法。包括步骤计算大气层外边界辐射强度Isun；计算地面处太阳直射辐射强度Id，0；计算地面处大气透过率τatm，0，τatm，0＝Id，0/Isun；根据地面处大气透过率τatm，0计算高度H处大气透过率τatm；根据H处大气透过率τatm、大气层外边界辐射强度Isun确定高度H处太阳直射辐射强度Id，Id＝Isun·τatm；计算高度H处地面长波辐射大气透过率τIR,atm；确定地表温度Tground，并根据地面长波辐射透过率和发射率及玻尔兹曼定律确定地面长波辐射强度。本发明专利技术根据较少的测量数据估算目标场地放飞的浮空器在不同高度H处所受辐射热环境特性参数，减小了高度和地理环境给浮空器热特性计算带来的误差，在浮空器热设计、蒙皮选材、飞行试验安全等方面具有指导意义。

A Method for Estimating Radiant Thermal Environment Characteristics of Aerostat

The invention belongs to the field of thermal environment analysis of a floater, in particular to a method for estimating the radiation thermal environment characteristics of a floater. Including steps to calculate the atmospheric external boundary radiation intensity Isun; calculate the solar direct radiation intensity Id, 0; calculate the atmospheric transmittance atm, 0, atm, 0=Id, 0/Isun; calculate the atmospheric transmittance ATM at ground level, 0 at altitude H; determine the solar direct radiation at altitude H according to the atmospheric transmittance ATM at H and the radiation intensity Isun at atmospheric external boundary. Radiation intensity Id, Id = Isun atm; Atmospheric transmittance_IR, ATM of ground long-wave radiation at altitude H; Surface temperature Tground is determined, and ground long-wave radiation intensity is determined according to ground long-wave radiation transmittance and emissivity and Boltzmann law. The invention estimates the radiation thermal environment characteristic parameters of the aerostat released from the target site at different altitudes H based on less measured data, reduces the error of the calculation of the thermal characteristics of the aerostat caused by altitude and geographical environment, and has guiding significance in thermal design of the aerostat, material selection of skin, flight test safety, etc.

【技术实现步骤摘要】
一种浮空器辐射热环境特性估算方法
本专利技术属于浮空器热环境分析领域，尤其涉及一种浮空器辐射热环境特性估算方法。
技术介绍
浮空器是一种依靠浮升气体升空、体积巨大的低动态飞行器，在通信中继、对地观测、空间探测等领域具有广阔应用前景。浮空器升降及驻空过程中，外界热环境变化导致内部浮升气体温度变化，会影响到其内部压力、飞行控制及安全。浮空器的环境热效应主要受对流和热辐射综合影响，其中热辐射主要受太阳直射辐射和地面长波辐射的影响。不同地区热环境迥异，尤其是辐射热环境存在较大差异，以往研究中多采用通用经验模型进行估算，难以反映地区热辐射环境差异给浮空器热特性带来的影响。戴秋敏(浮空器热环境与热特性研究，南京航空航天大学博士学位论文，2014年)根据气溶胶、水汽、臭氧、二氧化碳等参数的垂直分布情况拟合了新的太阳直射辐射强度和地面长波辐射强度模型，计算精度得到了提高，但该方法测量数据种类较多，对于一些测量条件有限的区域难以实现。
技术实现思路
本专利技术目的在于，提供一种浮空器辐射热环境特性估算方法，可根据较少的测量数据估算目标场地放飞的浮空器在不同高度处所受辐射热环境参数。为达到上述专利技术目的，本专利技术所述浮空器辐射热环境特性估算方法，包括如下步骤：S1，计算大气层外边界辐射强度Isun；S2，计算地面处太阳直射辐射强度Id，0；S3，计算地面处大气透过率τatm，0，τatm，0＝Id，0/Isun；S4，根据地面处大气透过率τatm，0计算高度H处大气透过率τatm；S41，高度H处大气透过率τatm与大气质量m呈指数关系，τatm＝0.5[exp(-am)+exp(-bm)]；其中，a，b均为待拟合系数；S42，将高度0处大气透过率，即τatm地面处大气透过率τatm，0带入公式τatm＝0.5[exp(-am)+exp(-bm)]，利用最小二乘法对地面处大气透过率τatm，0与大气质量m进行拟合，确定系数a，b的值；S34，将确定的系数a，b带入公式τatm＝0.5[exp(-am)+exp(-bm)]，得到高度H处大气透过率τatm＝0.5[exp(-am)+exp(-bm)]；S5，根据H处大气透过率τatm、大气层外边界辐射强度Isun确定高度H处太阳直射辐射强度Id，Id＝Isun·τatm；S6，计算高度H处地面长波辐射大气透过率τIR,atm；式中，P0为地面处大气压强，Patm为高度H处大气压强；S7，确定地表温度Tground，并根据地面长波辐射透过率和发射率及玻尔兹曼定律确定地面长波辐射强度。S71，地表温度根据气温Tatm确定，Tground＝cTatm+d；其中，c，d均为拟合系数；S72，根据地表温度Tground、地面长波发射率εground和高度H处地面长波辐射大气透过率τIR,atm计算高度H处地面长波辐射强度IIR,ground，式中，σ为斯蒂芬-玻尔兹曼常数，其值为σ＝5.67×10-8。进一步，步骤S41中根据高度H处大气压强Patm以及所对应的太阳高度角h计算大气质量m，大气质量m计算公式为：进一步，所述不同高度H的范围为地面至海拔32km处。本专利技术根据较少的测量数据估算目标场地放飞的浮空器在不同高度H处所受辐射热环境特性参数，减小了高度和地理环境给浮空器热特性计算带来的误差，在浮空器热设计、蒙皮选材、飞行试验安全等方面具有指导意义。附图说明图1为本专利技术实施例中一种浮空器辐射热环境特性估算方法的流程示意图。具体实施方式下面结合附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术的一部分实施例，而不是全部的实施例、基于本专利技术中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。实施例中所采用的太阳直射辐射表等设备均为市售。实施例1一种浮空器辐射热环境特性估算方法，如图1所示，包括如下步骤：S1，计算大气层外边界辐射强度Isun；根据不同日期的天序号n确定大气层外边界辐射强度Isun，大气层外边界辐射强度Isun计算公式为其中I0为太阳常数，其值为1367W/m2。S2，利用太阳直射辐射表计算地面处太阳直射辐射强度Id，0；S3，根据大气层外边界辐射强度Isun及地面处太阳直射辐射强度Id，0计算地面处大气透过率τatm，0，τatm，0＝Id，0/Isun；S4，根据地面处大气透过率τatm，0计算高度H处大气透过率τatm；S41，不同高度处大气压不同，从而大气质量不同，先根据测量时刻的高度H处大气压强Patm以及所对应的太阳高度角h计算大气质量m，大气质量m计算公式为：高度H处大气透过率τatm与大气质量m呈指数关系，τatm＝0.5[exp(-am)+exp(-bm)]；其中，a，b均为待拟合系数；S42，将高度0处大气透过率，即τatm地面处大气透过率τatm，0带入公式τatm＝0.5[exp(-am)+exp(-bm)]，利用最小二乘法对地面处大气透过率τatm，0与大气质量m进行拟合，确定系数a，b的值；S34，将确定的系数a，b带入公式τatm＝0.5[exp(-am)+exp(-bm)]，得到高度H处大气透过率τatm＝0.5[exp(-am)+exp(-bm)]；所述不同高度H的范围为地面至海拔32km处。S5，根据H处大气透过率τatm、大气层外边界辐射强度Isun确定高度H处太阳直射辐射强度Id，Id＝Isun·τatm；S6，计算高度H处地面长波辐射大气透过率τIR,atm；地面长波辐射透过率在地面处为1，透过率随海拔上升呈指数递减，与太阳直射辐射的大气透过率形式类似，故高度H处地面长波辐射大气透过率τIR,atm采用下式计算：式中，系数a，b的值与步骤S4中确定的值相同，P0为地面处大气压强，Patm为高度H处大气压强；优选地，所在高度H处的大气压强Patm根据下式确定：S7，确定地表温度Tground，并根据H处地面长波辐射大气透过率τIR,atm和地面长波发射率εground及玻尔兹曼定律确定地面长波辐射强度；S71，地表温度根据气温Tatm确定，Tground＝cTatm+d；其中，c，d均为拟合系数；S72，根据地表温度Tground、地面长波发射率εground和高度H处地面长波辐射大气透过率τIR,atm计算高度H处地面长波辐射强度IIR,ground，式中，σ为斯蒂芬-玻尔兹曼常数，其值为σ＝5.67×10-8。优选地，地面长波发射率εground根据不同地表类型确定：草地0.98，土壤0.93，城市0.92，冰雪0.99，水体0.95，沙漠0.88。综上可知，本专利技术中的浮空器辐射热环境特性估算方法实现了利用较少的测量参数对目标地区放飞的浮空器在不同高度处受到的太阳直射辐射强度和地面长波辐射强度进行估算的效果。本文档来自技高网...

【技术保护点】
1.一种浮空器辐射热环境特性估算方法，其特征在于，包括如下步骤：S1，计算大气层外边界辐射强度Isun；S2，计算地面处太阳直射辐射强度Id，0；S3，计算地面处大气透过率τatm，0，τatm，0＝Id，0/Isun；S4，根据地面处大气透过率τatm，0计算高度H处大气透过率τatm；S41，高度H处大气透过率τatm与大气质量m呈指数关系，τatm＝0.5[exp(‑am)+exp(‑bm)]；其中，a，b均为待拟合系数；S42，将高度0处大气透过率，即τatm地面处大气透过率τatm，0带入公式τatm＝0.5[exp(‑am)+ex(p‑bm)]，利用最小二乘法对地面处大气透过率τatm，0与大气质量m进行拟合，确定系数a，b的值；S34，将确定的系数a，b带入公式τatm＝0.5[exp(‑am)+exp(‑bm)]，得到高度H处大气透过率τatm＝0.5[exp(‑am)+exp(‑bm)]；S5，根据H处大气透过率τatm、大气层外边界辐射强度Isun确定高度H处太阳直射辐射强度Id，Id＝Isun·τatm；S6，计算高度H处地面长波辐射大气透过率τIR,atm；

【技术特征摘要】
1.一种浮空器辐射热环境特性估算方法，其特征在于，包括如下步骤：S1，计算大气层外边界辐射强度Isun；S2，计算地面处太阳直射辐射强度Id，0；S3，计算地面处大气透过率τatm，0，τatm，0＝Id，0/Isun；S4，根据地面处大气透过率τatm，0计算高度H处大气透过率τatm；S41，高度H处大气透过率τatm与大气质量m呈指数关系，τatm＝0.5[exp(-am)+exp(-bm)]；其中，a，b均为待拟合系数；S42，将高度0处大气透过率，即τatm地面处大气透过率τatm，0带入公式τatm＝0.5[exp(-am)+ex(p-bm)]，利用最小二乘法对地面处大气透过率τatm，0与大气质量m进行拟合，确定系数a，b的值；S34，将确定的系数a，b带入公式τatm＝0.5[exp(-am)+exp(-bm)]，得到高度H处大气透过率τatm＝0.5[exp(-am)+exp(-bm)]；S5，根据H处大气透过率τatm、大气层外边界辐射强度Isun确定高度H处太阳...

【专利技术属性】
技术研发人员：李敏，宁辉，孟小君，叶虎，燕道华，张宇，李勇，
申请(专利权)人：中国人民解放军六三六五三部队，
类型：发明
国别省市：新疆,65