CN110781579B - 一种空投物体下落过程仿真方法 - Google Patents
一种空投物体下落过程仿真方法 Download PDFInfo
- Publication number
- CN110781579B CN110781579B CN201910900309.9A CN201910900309A CN110781579B CN 110781579 B CN110781579 B CN 110781579B CN 201910900309 A CN201910900309 A CN 201910900309A CN 110781579 B CN110781579 B CN 110781579B
- Authority
- CN
- China
- Prior art keywords
- air
- aerial delivery
- flow field
- simulating
- speed
- 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
- 238000000034 method Methods 0.000 title claims abstract description 63
- 230000008569 process Effects 0.000 title claims abstract description 40
- 238000004088 simulation Methods 0.000 claims abstract description 25
- 230000003068 static effect Effects 0.000 claims abstract description 3
- 238000006073 displacement reaction Methods 0.000 claims description 12
- 230000005484 gravity Effects 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 230000001052 transient effect Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 description 5
- 238000012827 research and development Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 1
- 229910052753 mercury Inorganic materials 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种空投物体下落过程仿真方法,采用静态网格技术通过考虑流场来流和物体的相对关系模拟空投物体的下落过程,包括以下步骤:1)确定空投物体下落时的初始状态和初始流场边界条件;2)根据空投物体的受力情况确定下一时刻物体的运动状态;3)根据空投物体和来流的相对运动状态确定新时刻的流场边界条件;4)更新流场边界条件;5)判断空投物体是否落地,若空投物体已落地,则终止仿真,如未落地,则重复步骤2)‑4)。与现有技术相比,本发明具有计算量小、仿真时间短、结果准确等优点。
Description
技术领域
本发明涉及计算机仿真领域,尤其是涉及一种空投物体下落过程仿真方法。
背景技术
空投物体,如炸弹在研发过程中需要按照国家军用标准GJB 349.22-1989通过实际投弹测量空投物体的“标准落下时间”进行评估其下落性能,试验方法为在标准大气条件下(地面气压为760mm水银柱高,地面气温为15℃,地面大气密度为1.225kg/m3,气温递增率为0.0065℃/m)从2000m高度、航速为40m/s作水平飞行的飞机上投下的炸弹的落下时间。但上述过程执行周期长,成本高,由于炸弹在设计阶段中需要不断优化外形以提高其下落性能,因此上述试验方法难以满足实际研发的急切需要,为此可通过仿真方法研究炸弹的下落过程,缩短炸弹的研发设计时间。
利用仿真手段研究物体下落动态过程,一般可以通过动网格技术模拟物体在流场中的运动轨迹。但采用此种方法,需要创建一个较大的流场,相比于流场,物体本身的尺寸要小得多。如对于空投物体下落过程的仿真,需要创建一个高度达到2000m,长度数百米的流场,而空投物体的相对尺寸则小得多,需要大量的网格以保证仿真精度,计算量和仿真时间难以预测,普通计算机很难满足。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种空投物体下落过程仿真方法。
本发明的目的可以通过以下技术方案来实现:
一种空投物体下落过程仿真方法,采用静态网格技术通过考虑流场来流和物体的相对关系模拟空投物体的下落过程,包括以下步骤:
1)确定空投物体下落时的初始状态和初始流场边界条件;
2)根据空投物体的受力情况确定下一时刻物体的运动状态;
3)根据空投物体和来流的相对运动状态确定新时刻的流场边界条件;
4)更新流场边界条件;
5)判断空投物体是否落地,若空投物体已落地,则终止仿真,如未落地,则重复步骤2)-4)。
空投物体在下落过程中的自由度包括相对地面的竖直下落、水平方向的前进以及绕着空投物体重心所在轴线的俯仰。
所述的流场边界条件包括流场的速度矢量、大气密度、温度和压强,其由空投物体的运动状态确定,运动状态包括水平速度、垂直速度、垂直下落距离以及空投物体轴线与水平方向的夹角。
所述的步骤1)中,初始流场边界条件根据炸弹下落时的初始状态确定,则有:
其中,ux0为初始位置时空投物体的对地水平方向速度,uy0为初始位置时空投物体的对地竖直方向速度,vx0为初始位置时流场对空投物体的水平方向风速,vy0为初始位置时流场对空投物体竖直方向速度风速。
在仿真空投物体下落的过程中,以空投物体为参考系,则流场相对于空投物体的速度为:
其中,与分别为t时刻流场相对于空投物体的水平速度和垂直速度,α为空投物体速度矢量与空投物体轴线的夹角,即空投物体的迎风角,为t时刻空投物体在绝对坐标系中的水平方向速度,为t时刻空投物体在绝对坐标系中的垂直方向速度。
在仿真空投物体下落的过程中,空投物体从t时刻到t+Δt时刻的速度变化通过牛顿第二定律获取,则有:
其中,m为空投物体的质量,J为空投物体绕z轴的转动惯量,g为重力加速度,分别为t时刻空投物体沿运动方向的受到的法向力、切向力以及绕z轴的转矩,分别为t时刻和t+Δt时刻空投物体在绝对坐标系中的水平速度、垂直速度、绕z轴的转动速度,θ0为t时刻空投物体轴线与水平方向间的夹角。
空投物体在Δt时间内的水平位移Δx、竖直位移Δy和转角Δθ分别为:
空投物体从初始位置到新位置时的总位移通过位移求和得到。
该仿真为瞬态仿真,物体和流场的相对位置不发生改变。
空投物体下落的动态过程通过改变流场边界条件实现。
与现有技术相比,本发明具有以下优点:
本仿真方法误差小,采用本发明后,流场的计算域缩小,对计算机的计算性能要求降低,减少了计算量和仿真时间,缩短了研发周期,降低了研发成本。
附图说明
图1为本发明中空投物体下落过程动力学模型。
图2为本发明中流场相对空投物体速度示意图
图3为本发明中空投物体落过程的仿真求解路线图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
本实施例为通过相对运动和坐标变换的概念描述本发明的原理,并用此仿真方法对炸弹按照国家军用标准的规定进行“标准落下时间”的计算。
如图1所示,考虑炸弹在三个方向的自由度,分别是相对地面的竖直下落,水平方向的前进以及绕着炸弹重心所在轴线的俯仰,定义相对于地面静止的绝对坐标系oxy,初始位置时炸弹的对地水平方向速度为ux0,对地竖直方向速度为uy0;此时流场对炸弹的水平方向风速为vx0,对炸弹竖直方向速度风速为vy0,由相对位置关系,可知在初始位置
如图1所示,下落过程中,炸弹从t时刻到t+Δt时刻的速度变化通过牛顿第二定律求得
其中,m为炸弹的质量,J为炸弹绕z轴转动惯量,g为重力加速度, 分别为t时刻炸弹沿运动方向的受到的法向力、切向力以及受到的绕z轴转矩,以及分别为t时刻和t+Δt时刻炸弹在绝对坐标系中的水平速度、垂直速度及绕z轴的转动速度,θ0为t时刻炸弹轴线和水平方向的夹角。
如图2所示,以炸弹为参考系,t时刻流场相对于炸弹的速度为:
则炸弹在Δt时间内的水平位移、竖直位移和转角分别为
炸弹水平和竖直方向的总位移s和h即可通过上述位移求和得到。
炸弹下落过程中,周边的气压和气温及密度都会随着海拔降低而升高,视空气为理想气体。以海平面为基准,气温为15℃,大气压强为1.01325×105Pa,海拔高度每上升1m,大气压强减小11.1Pa,气温减小0.0065℃。则大气压强p(h)和气温T(h)可表示为炸弹竖直下落距离的函数,如式(5)所示,密度可通过理想气体状态方程求出:
如图2所示,根据上述动力学模型,可以采用fluent流场仿真对炸弹下落过程进行仿真。其中,炸弹在任意时刻的受力可通过对弹体表面压力求和得到,虚线框内的过程利用fluent的UDF功能编程实现。
基于上述仿真方法选择为层流模型,介质为理想气体考虑空气的可压缩性,边界条件设置为压力远场,其中流场的初始边界条件为:水平方向马赫数为0.1203(取当地音速为332.5m/s),垂直方向分量为0。大气压强为79102Pa,大气温度为2℃,仿真时间步长取0.01s,即每隔0.01s调整一次流场的边界条件。经过仿真计算得到该炸弹的标落下时间为20.72s,与实际投弹测量得到的标准落下时间误差较小,证实了仿真可以有效准确地模拟出炸弹投弹时下落过程。
本发明阐明了采用fluent流场仿真软件对某型炸弹下落过程进行仿真的原理和过程,并和实验对比分析了其计算标准落下时间的精确性。但本领域技术人员应很容易理解的是,本发明对于基于本仿真方法的而采用其他仿真软件如CFX等进行空投物体下落过程的仿真同样应纳入本发明的权利要求范围。
Claims (9)
1.一种空投物体下落过程仿真方法,采用静态网格技术通过考虑流场来流和物体的相对关系模拟空投物体的下落过程,其特征在于,包括以下步骤:
1)确定空投物体下落时的初始状态和初始流场边界条件;
2)根据空投物体的受力情况确定下一时刻物体的运动状态,在仿真空投物体下落的过程中,空投物体从t时刻到t+Δt时刻的速度变化通过牛顿第二定律获取,则有:
其中,m为空投物体的质量,J为空投物体绕z轴的转动惯量,g为重力加速度,分别为t时刻空投物体沿运动方向的受到的法向力、切向力以及绕z轴的转矩,分别为t时刻和t+Δt时刻空投物体在绝对坐标系中的水平速度、垂直速度、绕z轴的转动速度,θ0为t时刻空投物体轴线与水平方向间的夹角;
3)根据空投物体和来流的相对运动状态确定新时刻的流场边界条件;
4)更新流场边界条件;
5)判断空投物体是否落地,若空投物体已落地,则终止仿真,如未落地,则重复步骤2)-4)。
2.根据权利要求1所述的一种空投物体下落过程仿真方法,其特征在于,空投物体在下落过程中的自由度包括相对地面的竖直下落、水平方向的前进以及绕着空投物体重心所在轴线的俯仰。
3.根据权利要求1所述的一种空投物体下落过程仿真方法,其特征在于,所述的流场边界条件包括流场的速度矢量、大气密度、温度和压强,其由空投物体的运动状态确定,运动状态包括水平速度、垂直速度、垂直下落距离以及空投物体轴线与水平方向的夹角。
7.根据权利要求6所述的一种空投物体下落过程仿真方法,其特征在于,空投物体从初始位置到新位置时的总位移通过位移求和得到。
8.根据权利要求1所述的一种空投物体下落过程仿真方法,其特征在于,该仿真为瞬态仿真,物体和流场的相对位置不发生改变。
9.根据权利要求1所述的一种空投物体下落过程仿真方法,其特征在于,物体下落的动态过程通过改变流场边界条件实现。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910900309.9A CN110781579B (zh) | 2019-09-23 | 2019-09-23 | 一种空投物体下落过程仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910900309.9A CN110781579B (zh) | 2019-09-23 | 2019-09-23 | 一种空投物体下落过程仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110781579A CN110781579A (zh) | 2020-02-11 |
CN110781579B true CN110781579B (zh) | 2021-07-20 |
Family
ID=69383781
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910900309.9A Active CN110781579B (zh) | 2019-09-23 | 2019-09-23 | 一种空投物体下落过程仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110781579B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111428434B (zh) * | 2020-04-10 | 2022-07-01 | 山东大学 | 一种物体在洋流中下落的数值仿真方法及*** |
CN114608384A (zh) * | 2020-12-08 | 2022-06-10 | 江苏数字鹰科技股份有限公司 | 一种空投物落点在瞄具中位置的确定方法以及相应的瞄具 |
CN113155393B (zh) * | 2021-03-03 | 2023-03-28 | 中国人民解放军95795部队 | 一种空投缓冲气囊试验装置 |
CN113688543A (zh) * | 2021-07-21 | 2021-11-23 | 昆山丘钛微电子科技股份有限公司 | 一种跌落仿真方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103500144A (zh) * | 2013-09-29 | 2014-01-08 | 北京航空航天大学 | 一种基于动态边界条件的助推器再入运动分析方法 |
CN103559335A (zh) * | 2013-10-16 | 2014-02-05 | 南京航空航天大学 | 一种基于流固耦合的伞降仿真方法 |
CN104504174A (zh) * | 2014-11-28 | 2015-04-08 | 南京航空航天大学 | 一种适于伞降空投过程流固耦合数值模拟的自适应网格控制方法 |
WO2017004419A1 (en) * | 2015-06-30 | 2017-01-05 | Unmanned Innovations, Inc. | Systems and methods for multi-mode unmanned vehicle mission planning and control |
CN109165423A (zh) * | 2018-08-03 | 2019-01-08 | 北京航空航天大学 | 一种基于流函数的绕物流场建模方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ATE534538T1 (de) * | 2008-02-14 | 2011-12-15 | Sumitomo Rubber Ind | Verfahren zur parameterbestimmung, das beim luftdruckabnahme-erkennungsverfahren verwendet wird |
CN103559338B (zh) * | 2013-10-21 | 2016-06-15 | 南京航空航天大学 | 一种降落伞工作特性仿真方法 |
CN104166771B (zh) * | 2014-08-26 | 2017-07-18 | 中国人民解放军装甲兵工程学院 | 复杂环境下空投着陆仿真模拟与评价方法 |
CN105468808B (zh) * | 2015-09-10 | 2019-02-19 | 中国人民解放军空军空降兵学院 | 空降空投并行蒙特卡洛仿真的实现方法 |
-
2019
- 2019-09-23 CN CN201910900309.9A patent/CN110781579B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103500144A (zh) * | 2013-09-29 | 2014-01-08 | 北京航空航天大学 | 一种基于动态边界条件的助推器再入运动分析方法 |
CN103559335A (zh) * | 2013-10-16 | 2014-02-05 | 南京航空航天大学 | 一种基于流固耦合的伞降仿真方法 |
CN104504174A (zh) * | 2014-11-28 | 2015-04-08 | 南京航空航天大学 | 一种适于伞降空投过程流固耦合数值模拟的自适应网格控制方法 |
WO2017004419A1 (en) * | 2015-06-30 | 2017-01-05 | Unmanned Innovations, Inc. | Systems and methods for multi-mode unmanned vehicle mission planning and control |
CN109165423A (zh) * | 2018-08-03 | 2019-01-08 | 北京航空航天大学 | 一种基于流函数的绕物流场建模方法 |
Non-Patent Citations (3)
Title |
---|
Flow fields of graupel falling in air;Pao K. Wang,et al;《Atmospheric Research》;20130428;正文第158-169页 * |
炸弹投放模型的建立与数值模拟;陈少君 等;《电光与控制》;20120731;正文第37-40、56页 * |
降落伞充满稳定阶段流场数值模拟;刘巍 等;《中国空间科学技术》;20070228;正文第61-64页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110781579A (zh) | 2020-02-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110781579B (zh) | 一种空投物体下落过程仿真方法 | |
CN107679655B (zh) | 一种航天发射火箭落点预测*** | |
CN106918439B (zh) | 一种基于风洞的***模拟器、其运行方法及其所得***模型 | |
CN102393200B (zh) | 基于飞行仿真的通用惯导测试方法 | |
CN109710961B (zh) | 一种基于gps数据的高空无人机升限数据处理方法 | |
CN105137804B (zh) | 一种针对飞行姿态扰动的实验室模拟方法 | |
CN106933104A (zh) | 一种基于dic‑pid的四旋翼飞行器姿态与位置的混合控制方法 | |
CN207456730U (zh) | 一种基于风洞的***模拟器 | |
CN103868648B (zh) | 三轴气浮仿真实验平台的质心测量方法 | |
CN106096091B (zh) | 一种飞机运动模拟方法 | |
CN113848963B (zh) | 一种飞行控制***的控制律参数设计方法 | |
CN103234729A (zh) | 风洞常规测力试验中气动刚度与气动阻尼的视频测量方法 | |
CN114186508A (zh) | 一种基于cfd软件的水下航行器水动力系数测算方法 | |
CN104197955B (zh) | 全自动三轴气浮台微干扰力矩测量方法 | |
CN106570242B (zh) | 低空大动压整体式整流罩高速分离流固耦合仿真方法 | |
CN115791075A (zh) | 多旋翼风阻系数标定方法、设备及计算机可读存储介质 | |
CN106290969A (zh) | 一种考虑减速伞气动力影响的风速风向探测方法 | |
CN107315419A (zh) | 一种实现飞艇远距离定点悬停的三段式控制策略 | |
CN105973237B (zh) | 基于实际飞行数据插值的仿真动态轨迹解析生成方法 | |
CN102566446A (zh) | 基于线性模型组的无人直升机全包线数学模型构建方法 | |
CN112668107B (zh) | 基于立式转台的航天运载器试验方法、装置及存储介质 | |
CN116992700B (zh) | 一种物流无人机导航精度确定的方法及设备 | |
CN114139797A (zh) | 光伏板表面积灰预测***及方法、存储介质及电子设备 | |
CN108073180A (zh) | 无人机的控制方法、装置和*** | |
CN116702439B (zh) | 基于再入落点预测模型的仿真用飞行器再入落点预测方法 |
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 |