CN113656896A - 翻转平板在陨落过程中的气动力模型建立方法 - Google Patents

翻转平板在陨落过程中的气动力模型建立方法 Download PDF

Info

Publication number
CN113656896A
CN113656896A CN202111062927.4A CN202111062927A CN113656896A CN 113656896 A CN113656896 A CN 113656896A CN 202111062927 A CN202111062927 A CN 202111062927A CN 113656896 A CN113656896 A CN 113656896A
Authority
CN
China
Prior art keywords
flat plate
aerodynamic model
aerodynamic
model
coefficient
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.)
Granted
Application number
CN202111062927.4A
Other languages
English (en)
Other versions
CN113656896B (zh
Inventor
张兆
李志辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Original Assignee
High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center filed Critical High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Priority to CN202111062927.4A priority Critical patent/CN113656896B/zh
Publication of CN113656896A publication Critical patent/CN113656896A/zh
Application granted granted Critical
Publication of CN113656896B publication Critical patent/CN113656896B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Mathematical Optimization (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明公开了一种翻转平板在陨落过程中的气动力模型建立方法,包括如下步骤:步骤一、对平板的翻滚运动进行分解;步骤二、计算当地攻角;步骤三、计算全流域下的压力系数;步骤四、计算当前名义攻角下瞬时的升力系数和阻力系数;步骤五、建立平板在翻滚陨落时的气动力模型。与现有技术相比,本发明的积极效果是:本发明方法针对典型的平板模型,通过速度分解充分考虑了平板模型绕长轴滚转时不同姿态角下滚转速度分量对当地攻角的影响,导出了该条件下的瞬时压力分布和气动力,然后积分得出了相应的平均气动力和气动力矩,从而建立了平板翻滚运动下的气动力模型,这对平板陨落轨迹的预测和落点分析、以及降低地面风险具有积极的意义。

Description

翻转平板在陨落过程中的气动力模型建立方法
技术领域
本发明涉及空气动力学与航天技术领域,具体涉及一种翻转平板在陨落过程中的气动力模型建立方法。
背景技术
大型航天器在轨服役多年后,为了避免形成太空垃圾,影响现役的航天器安全,需要对其进行离轨再入返回地面,整个陨落过程涉及弹道、气动力/热,烧蚀等多个学科,系统复杂度极高。这些大型航天器在再入过程中与大气层发生摩擦,在强烈的气动加热和气动力作用下解体产生大量碎片,由于高温作用部分碎片会在到达地面前完全烧蚀熔化,但最终仍有一些碎片落到地面。这些陨落的碎片一方面会对地面的人群、建筑以及生态系统造成威胁,另一方面一些有价值的碎片还需要进行回收研究。因此准确预报空间碎片陨落时间、地点对于降低航天器碎片的危险度和及时回收残骸部件是十分重要的。
对于大型航天器而言,其外部壳层及内部的舱室等空间结构含有大量的薄壁结构,当航天器服役期满后再入返回大气层时,这些结构将发生解体、破碎,必然产生大量的片状碎片。从形状上讲,这是真正意义上的太空碎“片”。因此在再入过程中,采用与实际碎片形状更为接近的片状模型,以取得准确性更高的再入弹道数据和预测落点位置似乎是一种必然。但在实际的应用中,由于片状模型的气动力、气动热计算复杂,计算周期长,且再入碎片一般存在的自旋现象又进一步增加了再入弹道预测的难度,故在再入弹道计算程序中实现较为困难。
目前常见的再入预测方法可以分为面向物体法和面向飞行器法两大类,如图1所示:(a)面向物体法将飞行器及其部件简化为若干基本外形物体,利用基本外形的气动力/热和烧蚀模型进行再入仿真,具有建模简单、计算速度快的优点。(b)面向飞行器法则直接对飞行器真实外形进行建模并采用工程算法计算飞行器的气动力和气动热,能够更准确地模拟再入烧蚀和解体等过程,但是这类方法计算量大,不能快速获得,影响了预测的实效性。
对于再入运动,飞行器的质量和惯性张量、空间姿态和运动状态作为与气动力模型密切相关的重要输入参数,决定了弹道方程中气动力的计算可靠性,影响着再入轨迹预测的准确度。本发明提供了一种典型的平板碎片在陨落过程中处于翻滚状态时的气动力模型。
发明内容
为了克服现有技术的上述缺点,本发明提出了一种翻转平板在陨落过程中的气动力模型建立方法。
本发明解决其技术问题所采用的技术方案是:一种翻转平板在陨落过程中的气动力模型建立方法,包括如下步骤:
步骤一、对平板的翻滚运动进行分解;
步骤二、计算当地攻角;
步骤三、计算全流域下的压力系数;
步骤四、计算当前名义攻角下瞬时的升力系数和阻力系数;
步骤五、建立平板在翻滚陨落时的气动力模型。
与现有技术相比,本发明的积极效果是:
典型的空间平板碎片进入大气层后的陨落过程,平板碎片在陨落过程中会产生翻滚运动,这种情况下的气动力如果采用计算流体力学的方法直接求解流场来获得平板所受的气动力和气动力矩,由于物理过程十分复杂且计算量巨大,难以快速获得相应的参数。工程中的做法就是采用等效攻角的经验模型来直接估算物体所受的气动力,这类方法并没有充分考虑物体空间姿态和运动形态(如进动、章动、翻滚等旋转运动)对气动力的影响,导致等效攻角模型预测的气动力/力矩的精度不够。
对于大型航天器而言,其太阳能帆板、外部壳层及内部的舱室等空间结构含有大量的薄壁结构,因此片状平板模型是航天器中最常见的部件。另外根据刚体动力学的一般性理论可知,对于薄的长宽比较大的长方体刚体而言,其绕长轴的翻滚运动是常见的且稳定的运动形态。基于此本发明方法针对典型的平板模型,通过速度分解充分考虑了平板模型绕长轴滚转时不同姿态角下滚转速度分量对当地攻角的影响,导出了该条件下的瞬时压力分布和气动力,然后积分得出了相应的平均气动力和气动力矩,从而建立了平板翻滚运动下的气动力模型,这对平板陨落轨迹的预测和落点分析、以及降低地面风险具有积极的意义。
附图说明
本发明将通过例子并参照附图的方式说明,其中:
图1为快速工程算法原理图,其中:(a)为面向物体法,(b)为面向航天器法;
图2为片状物的运动分解示意图;
图3为翻滚运动的当地攻角计算。
具体实施方式
空间飞行器中含有大量的板壳结构,在其陨落后都变成了薄板进入大气层。若要预测这类片状陨落物的落点,则需要对其陨落过程中的气动力进行建模,然后带入弹道动力学方程中分析其陨落轨迹。基于此,本发明以典型的平板为对象,提出了建立其在陨落过程中的典型翻滚运动状态下的一种气动力模型,用于其陨落轨迹的预测。
具体包含以下内容:
一、平板翻滚运动的分解与当地攻角的计算
这里为简单起见假定碎片为长宽厚分别为a,b,w,且满足
Figure 817288DEST_PATH_IMAGE001
的矩形,如航天器的太阳能帆板。片状物上的流动可以简单的视为平板边界层流动,这类片状模型也比较容易推广到多个面的部件上去,有助于复杂流动的分析。这里假定碎片是绕过其最小惯量主对称轴的翻滚运动,因此其运动可以分解为质心的平动和自旋的翻转运动,如图2所示。这时其质心的平动速度V沿绝对坐标系的X轴,翻滚的角速度
Figure 547477DEST_PATH_IMAGE002
指向Z轴。
在平板的截面内看其陨落过程,如图3所示,那么其短轴b方向上的微元运动可以分解为质心的平动叠加上绕质心的滚转运动,由于各个微元的周向速度不一样,因此各微元相对于气流的当地攻角也不同。对于如图3中的微元dr,首先具有与质心相同的平移速度,然后有由于翻转产生的相对速度
Figure 327215DEST_PATH_IMAGE003
,两者速度的合成速度才是微元的真实速度,同时该速度相对于平板的夹角才是流动的真实攻角,具体可以计算为:
Figure 784741DEST_PATH_IMAGE004
Figure 599113DEST_PATH_IMAGE006
二、全流域下压力系数C P的计算
飞行器在陨落过程中,当进入大气层时经历的气体从稀薄进入稠密,其流动也从稀薄流逐步演变为连续流,这过程用Kn数来衡量,它表示为分子自由程与参考长度之比。当Kn>10时,表明气体的分子自由程大于特征长度,为典型的稀薄流条件。即从物体表面反射的空气分子要运动到距离物体很远的地方才会和其它气体分子碰撞,此时分子间的碰撞过程可以忽略,只考虑气体分子与物体壁面的碰撞作用即可。当Kn<0.001时,表明气体的分子自由程远小于特征长度,满足连续流假设,其流动控制方程为Navier-Stokes方程,但是考虑到该方程难以求解,且陨落过程一般属于典型的高超声速流动问题,所以这里采用牛顿流模型来计算平板上的压力分布。
在稀薄流区(Kn≥10):假定Maxwell平衡气体分布,通过求解无碰撞的Boltzmann方程得到物体壁面的压力,可得相应的稀薄流区压力系数为:
Figure 234625DEST_PATH_IMAGE007
Figure 501658DEST_PATH_IMAGE009
在连续流区(Kn≤0.001):采用牛顿流模型,其基本思想是将流体质点看作是直线运动的粒子,碰撞到壁面后法向动量完全损失而切线动量保留下来。但是这种原始的牛顿流模型没有考虑有限马赫数的影响,这里采用修正牛顿公式,则连续流区的表面压力系数为:
Figure 700558DEST_PATH_IMAGE010
过渡流区(0.001<Kn<10),位于连续流区和稀薄流区之间,由于流场比较稀薄使得传统的连续流理论不适用,而分子间的自由碰撞既不能忽略也不能采用简单的稀薄流的理论公式。通常采用工程上常用的桥函数法进行加权叠加来连接两个区间,得到过渡区的气动参数。这里采用sine-squared桥函数作为加权函数,其中权函数的自变量就是Kn数,将连续流区和稀薄流区进行加权叠加得到过渡区的压力系数C P,trans为:
Figure 431754DEST_PATH_IMAGE011
三、翻滚运动的气动力模型
由陨落运动形态的假定和速度的分解,可知当片状物翻滚再入时,其角速度矢量w总是垂直速度方向。平板局部微元dr的真实速度V为相对的翻滚转动速度叠加上来流速度,并以此来计算当地攻角α,如式(1),可知它是名义攻角α m、平板短轴方向坐标r、翻转角速度w、来流速度V,可记为
Figure 425118DEST_PATH_IMAGE012
根据当前陨落的Kn数,确定当前是属于稀薄流、连续流还是过渡流,然后由(2~4)式可以计算出该局部微元dr上的压力分布C p ,如(5)式;
Figure 989567DEST_PATH_IMAGE013
Figure 726578DEST_PATH_IMAGE014
然后沿整个平板的迎风面对压力系数CPm,r)进行积分,同时沿X方向(阻力方向)和Y方向(升力方向)进行分解,可以得到当前名义攻角下瞬时的升力系数和阻力系数,如(6)式:
Figure 312281DEST_PATH_IMAGE015
由此可见,整个平板的瞬时气动升力和阻力与名义攻角α m密切相关。平板陨落翻滚时,其名义攻角不断变化,并认为翻滚的角速度为常数,即
Figure 742125DEST_PATH_IMAGE016
。因此其气动力是随时间做周期性变化,将其在一个时间周期内积分后,可得相应的平均气动系数:
Figure 983750DEST_PATH_IMAGE017
这样,就建立了平板在翻滚陨落时的气动力模型:
Figure 71923DEST_PATH_IMAGE018
该模型与滚翻角速度
Figure 715394DEST_PATH_IMAGE019
、流域的Kn数,陨落M数相关。

Claims (8)

1.一种翻转平板在陨落过程中的气动力模型建立方法,其特征在于:包括如下步骤:
步骤一、对平板的翻滚运动进行分解;
步骤二、计算当地攻角;
步骤三、计算全流域下的压力系数;
步骤四、计算当前名义攻角下瞬时的升力系数和阻力系数;
步骤五、建立平板在翻滚陨落时的气动力模型。
2.根据权利要求1所述的翻转平板在陨落过程中的气动力模型建立方法,其特征在于:按如下公式计算当地攻角:
Figure 47799DEST_PATH_IMAGE001
式中
Figure 839037DEST_PATH_IMAGE002
是质心的平动速度,
Figure 527639DEST_PATH_IMAGE003
是平板相对于来流的名义攻角,r是平板短轴b的坐标。
3.根据权利要求1所述的翻转平板在陨落过程中的气动力模型建立方法,其特征在于:按如下公式计算全流域下的压力系数:
Figure 291195DEST_PATH_IMAGE005
式中:C p,con为连续流区的表面压力系数;C p,trans为过渡区的压力系数;C p,fm为稀薄流区压力系数。
4.根据权利要求3所述的翻转平板在陨落过程中的气动力模型建立方法,其特征在于:按如下公式计算连续流区的表面压力系数:
Figure 62842DEST_PATH_IMAGE007
式中,M是来流马赫数,
Figure 826530DEST_PATH_IMAGE008
是气体分子的比热比,
Figure 709035DEST_PATH_IMAGE009
是当地攻角。
5.根据权利要求4所述的翻转平板在陨落过程中的气动力模型建立方法,其特征在于:按如下公式计算稀薄流区压力系数:
Figure 61519DEST_PATH_IMAGE011
式中,
Figure 613855DEST_PATH_IMAGE012
为来流和壁面温度,
Figure 51789DEST_PATH_IMAGE013
是壁面法向动量适应系数,
Figure 737985DEST_PATH_IMAGE014
是误差函数,S按下式计算得到:
Figure 210555DEST_PATH_IMAGE015
式中,
Figure 386321DEST_PATH_IMAGE016
为Boltzmann常数,m为气体分子质量。
6.根据权利要求5所述的翻转平板在陨落过程中的气动力模型建立方法,其特征在于:按如下公式计算过渡区的压力系数C P,trans
Figure 311552DEST_PATH_IMAGE017
式中,Kn为表征流动稀薄程度的努森数。
7.根据权利要求3所述的翻转平板在陨落过程中的气动力模型建立方法,其特征在于:按如下公式计算当前名义攻角下瞬时的升力系数和阻力系数:
Figure 345980DEST_PATH_IMAGE018
8.根据权利要求7所述的翻转平板在陨落过程中的气动力模型建立方法,其特征在于:步骤五所述建立平板在翻滚陨落时的气动力模型的方法为:
(1)按如下公式计算翻滚运动时的平均升力系数和平均阻力系数:
Figure 469793DEST_PATH_IMAGE019
(2)代入全流域下的压力系数得到与滚翻角速度
Figure 363931DEST_PATH_IMAGE020
、流域的Kn数,陨落M数相关的气动力模型:
Figure 776458DEST_PATH_IMAGE021
CN202111062927.4A 2021-09-10 2021-09-10 翻转平板在陨落过程中的气动力模型建立方法 Active CN113656896B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111062927.4A CN113656896B (zh) 2021-09-10 2021-09-10 翻转平板在陨落过程中的气动力模型建立方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111062927.4A CN113656896B (zh) 2021-09-10 2021-09-10 翻转平板在陨落过程中的气动力模型建立方法

Publications (2)

Publication Number Publication Date
CN113656896A true CN113656896A (zh) 2021-11-16
CN113656896B CN113656896B (zh) 2023-05-05

Family

ID=78483659

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111062927.4A Active CN113656896B (zh) 2021-09-10 2021-09-10 翻转平板在陨落过程中的气动力模型建立方法

Country Status (1)

Country Link
CN (1) CN113656896B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116306185A (zh) * 2023-02-14 2023-06-23 北京理工大学 基于不规则离散元的陨星进入地球大气层后运动仿真方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10118696B1 (en) * 2016-03-31 2018-11-06 Steven M. Hoffberg Steerable rotating projectile
CN111241634A (zh) * 2019-11-19 2020-06-05 中国空气动力研究与发展中心超高速空气动力研究所 一种航天器再入陨落的分析预报方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10118696B1 (en) * 2016-03-31 2018-11-06 Steven M. Hoffberg Steerable rotating projectile
CN111241634A (zh) * 2019-11-19 2020-06-05 中国空气动力研究与发展中心超高速空气动力研究所 一种航天器再入陨落的分析预报方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116306185A (zh) * 2023-02-14 2023-06-23 北京理工大学 基于不规则离散元的陨星进入地球大气层后运动仿真方法及系统
CN116306185B (zh) * 2023-02-14 2023-11-03 北京理工大学 基于不规则离散元的陨星进入地球大气层后运动仿真方法及系统

Also Published As

Publication number Publication date
CN113656896B (zh) 2023-05-05

Similar Documents

Publication Publication Date Title
Kinney Aero-thermodynamics for conceptual design
Shinan et al. Numerical investigation on impingement dynamics and freezing performance of micrometer-sized water droplet on dry flat surface in supercooled environment
CN109558660B (zh) 一种航天器碎片陨落落区预报方法
CN114580224B (zh) 一种分布式气动融合轨道耦合姿态摄动分析方法
Tong et al. Eulerian simulations of icing collection efficiency using a singularity diffusion model
CN116610137B (zh) 一种基于干扰预示的高超声速飞行器强抗扰控制方法
Xie et al. Robust and efficient prediction of the collection efficiency in icing accretion simulation for 3D complex geometries using the Lagrangian approach I: an adaptive interpolation method based on the restricted radial basis functions
CN113656896A (zh) 翻转平板在陨落过程中的气动力模型建立方法
Li et al. Hydrodynamic analysis of the energy dissipation of droplets on vibrating superhydrophobic surfaces
Scott et al. Navier-stokes solution to the flowfield over ice accretion shapes
Tewari Entry trajectory model with thermomechanical breakup
Xiao et al. A multiple model particle filter for maneuvering target tracking based on composite sampling
Shinan et al. Numerical simulation of ice shedding motion characteristic on airfoil surface
Harish et al. Store separation dynamics using grid-free Euler solver
Lofthouse et al. Static and dynamic simulations of a generic UCAV geometry using the kestrel flow solver
Lim et al. Prediction of HART II rotor BVI loading and wake system using CFD/CSD loose coupling
Dash et al. Computational modelling of freezing of water droplet impacting on an ice surface
Zhang et al. Aerodynamic thermal analysis of a spinning winged projectile
Bragg Effect of geometry on airfoil icing characteristics
Wada et al. Evaluation of aerodynamic performance of airfoil using the e-mps method after icing
Sathyanarayana et al. Trajectory prediction using coupled CFD-RBD with dynamic meshing
Iuliano et al. An eulerian approach to three-dimensional droplet impingement simulation in icing environment
Boudreau Artificially induced boundary-layer transition on blunt-slender cones at hypersonic speeds
Bragg Predicting rime ice accretion on airfoils
Bragg et al. An analytical evaluation of the icing properties of several low and medium speed airfoils

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
GR01 Patent grant
GR01 Patent grant