CN113656896A - 翻转平板在陨落过程中的气动力模型建立方法 - Google Patents
翻转平板在陨落过程中的气动力模型建立方法 Download PDFInfo
- 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
Links
Images
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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling 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,且满足的矩形,如航天器的太阳能帆板。片状物上的流动可以简单的视为平板边界层流动,这类片状模型也比较容易推广到多个面的部件上去,有助于复杂流动的分析。这里假定碎片是绕过其最小惯量主对称轴的翻滚运动,因此其运动可以分解为质心的平动和自旋的翻转运动,如图2所示。这时其质心的平动速度V沿绝对坐标系的X轴,翻滚的角速度指向Z轴。
在平板的截面内看其陨落过程,如图3所示,那么其短轴b方向上的微元运动可以分解为质心的平动叠加上绕质心的滚转运动,由于各个微元的周向速度不一样,因此各微元相对于气流的当地攻角也不同。对于如图3中的微元dr,首先具有与质心相同的平移速度,然后有由于翻转产生的相对速度,两者速度的合成速度才是微元的真实速度,同时该速度相对于平板的夹角才是流动的真实攻角,具体可以计算为:
二、全流域下压力系数C P的计算
飞行器在陨落过程中,当进入大气层时经历的气体从稀薄进入稠密,其流动也从稀薄流逐步演变为连续流,这过程用Kn数来衡量,它表示为分子自由程与参考长度之比。当Kn>10时,表明气体的分子自由程大于特征长度,为典型的稀薄流条件。即从物体表面反射的空气分子要运动到距离物体很远的地方才会和其它气体分子碰撞,此时分子间的碰撞过程可以忽略,只考虑气体分子与物体壁面的碰撞作用即可。当Kn<0.001时,表明气体的分子自由程远小于特征长度,满足连续流假设,其流动控制方程为Navier-Stokes方程,但是考虑到该方程难以求解,且陨落过程一般属于典型的高超声速流动问题,所以这里采用牛顿流模型来计算平板上的压力分布。
在稀薄流区(Kn≥10):假定Maxwell平衡气体分布,通过求解无碰撞的Boltzmann方程得到物体壁面的压力,可得相应的稀薄流区压力系数为:
在连续流区(Kn≤0.001):采用牛顿流模型,其基本思想是将流体质点看作是直线运动的粒子,碰撞到壁面后法向动量完全损失而切线动量保留下来。但是这种原始的牛顿流模型没有考虑有限马赫数的影响,这里采用修正牛顿公式,则连续流区的表面压力系数为:
过渡流区(0.001<Kn<10),位于连续流区和稀薄流区之间,由于流场比较稀薄使得传统的连续流理论不适用,而分子间的自由碰撞既不能忽略也不能采用简单的稀薄流的理论公式。通常采用工程上常用的桥函数法进行加权叠加来连接两个区间,得到过渡区的气动参数。这里采用sine-squared桥函数作为加权函数,其中权函数的自变量就是Kn数,将连续流区和稀薄流区进行加权叠加得到过渡区的压力系数C P,trans为:
三、翻滚运动的气动力模型
由陨落运动形态的假定和速度的分解,可知当片状物翻滚再入时,其角速度矢量w总是垂直速度方向。平板局部微元dr的真实速度V为相对的翻滚转动速度叠加上来流速度,并以此来计算当地攻角α,如式(1),可知它是名义攻角α m、平板短轴方向坐标r、翻转角速度w、来流速度V∞,可记为。
根据当前陨落的Kn数,确定当前是属于稀薄流、连续流还是过渡流,然后由(2~4)式可以计算出该局部微元dr上的压力分布C p ,如(5)式;
然后沿整个平板的迎风面对压力系数CP(αm,r)进行积分,同时沿X方向(阻力方向)和Y方向(升力方向)进行分解,可以得到当前名义攻角下瞬时的升力系数和阻力系数,如(6)式:
由此可见,整个平板的瞬时气动升力和阻力与名义攻角α m密切相关。平板陨落翻滚时,其名义攻角不断变化,并认为翻滚的角速度为常数,即。因此其气动力是随时间做周期性变化,将其在一个时间周期内积分后,可得相应的平均气动系数:
这样,就建立了平板在翻滚陨落时的气动力模型:
Claims (8)
1.一种翻转平板在陨落过程中的气动力模型建立方法,其特征在于:包括如下步骤:
步骤一、对平板的翻滚运动进行分解;
步骤二、计算当地攻角;
步骤三、计算全流域下的压力系数;
步骤四、计算当前名义攻角下瞬时的升力系数和阻力系数;
步骤五、建立平板在翻滚陨落时的气动力模型。
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116306185A (zh) * | 2023-02-14 | 2023-06-23 | 北京理工大学 | 基于不规则离散元的陨星进入地球大气层后运动仿真方法及系统 |
Citations (2)
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 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种航天器再入陨落的分析预报方法 |
-
2021
- 2021-09-10 CN CN202111062927.4A patent/CN113656896B/zh active Active
Patent Citations (2)
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)
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 |