CN105468808B - 空降空投并行蒙特卡洛仿真的实现方法 - Google Patents
空降空投并行蒙特卡洛仿真的实现方法 Download PDFInfo
- Publication number
- CN105468808B CN105468808B CN201510573559.8A CN201510573559A CN105468808B CN 105468808 B CN105468808 B CN 105468808B CN 201510573559 A CN201510573559 A CN 201510573559A CN 105468808 B CN105468808 B CN 105468808B
- Authority
- CN
- China
- Prior art keywords
- drop
- air
- airborne
- wind
- monte carlo
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000000342 Monte Carlo simulation Methods 0.000 title claims abstract description 16
- 238000004458 analytical method Methods 0.000 claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 239000006185 dispersion Substances 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 6
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000002948 stochastic simulation Methods 0.000 description 1
- 238000013076 uncertainty analysis Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供了一种空降空投并行蒙特卡洛仿真的实现方法,先给出初始条件的期望值,再给出概率分布规律,蒙特卡洛分析模块将根据以上变量对应的概率分布情况随机生成大量初始条件,并导出到主程序的输入文件中供主程序进行计算;进行大量的循环计算后,分析空中风向和风速,空投初始时刻的重量、高度、速度对计算结果的影响。本发明充分考虑了每次空降空投初始时刻的重量、高度、速度和空降空投过程中的风向、风速大小都的不确定性,使空降空投仿真计算更接近真实情况。
Description
技术领域
本发明涉及空降空投过程仿真实验领域,具体是一种空降空投并行蒙特卡洛仿真的实现方法。
背景技术
空降空投过程仿真,其本质为常微分初值问题,初始条件和空降空投过程中的风干扰情况唯一地决定了空降空投过程中空投件的轨迹、姿态、落点的情况。但是在实际空降空投过程中,初始条件不是恒定不变的,而是在一定范围内变化的随机变量。随机因素的存在使得每次空降空投初始时刻的重量、高度、速度和空降空投过程中的风向、风速大小都是不确定的,要对这些不同情况的影响进行模拟和统计,需要采用蒙特卡洛(Monte Carlo)方法。
蒙特卡洛法也称随机模拟,它主要依据概率分布对随机变量进行抽样,然后将样本带入数学模型进行计算得到应变量。虽然蒙特卡洛模拟技术只给出的是统计估计而非精确的结果,且应用其研究问题需要花费大量的计算时间,但它对问题的维数不敏感,对求解对象是线性问题与否也没有原则性要求,因此在复杂系统的不确定分析中,蒙特卡洛方法成为不可或缺的手段。而且对于那些无法得到解析结果的复杂问题来说,这种手段可能是唯一有效的结果。
蒙特卡洛法的实质是通过大量随机试验,利用概率论解决问题的一种数值方法,基本思想是基于概率和体积间的相似性。国内外对空降空投过程进行仿真的过程中,由于每次空降空投初始时刻的重量、高度、速度和空降空投过程中的风向、风速大小都具有的不确定性,仿真结果偏差较大。如何实现利用蒙特卡洛法实现对空降空投过程的仿真,是一个亟需解决的问题。
发明内容
本发明为了对空降空投过程进行仿真,提供了一种空降空投并行蒙特卡洛仿真的实现方法,使空降空投仿真计算更接近真实情况,并且大大缩短了计算时间
本发明提供了一种空降空投并行蒙特卡洛仿真的实现方法,包括以下步骤:
1)给出初始条件的期望值,再给出概率分布规律,将空降空投中的风向和风速大小在一定范围内随机均匀分布,而空降空投初始时刻的重量、高度、速度服从正态分布,其中各随机变量概率分布密度如下:
风速大小:
式中:Vmin和Vmax分别为某一高度上最小和最大风速;
风速与初始航向夹角:
式中:当θ风=π时为逆风,θ风=-π时为顺风;
空降空投重量:
式中:σm为空降空投重量分布方差,mμ为期望空降空投重量;对于某型空投件,默认取mμ=150kg,σm=5kg。
空降空投高度:
式中:σh为空降空投高度分布方差,hμ为期望空降空投高度;对于某型空投件,默认取hμ=400m,σh=17m。
空降空投速度:
式中:σMa为空降空投初始马赫数分布方差,Maμ为期望空降空投马赫数;对于某型空投件,默认取Maμ=0.35,σMa=0.02。
2)蒙特卡洛分析模块将根据以上变量对应的概率分布情况随机生成大量初始条件,并导出到主程序的输入文件中供主程序进行计算;
3)进行大量的循环计算,分析空投落点散布时取500~5000次,分析空中风向和风速,空投初始时刻的重量、高度、速度对计算结果的影响。
在进行蒙特卡洛仿真时,由于每一次计算的情况与之前设置固定参数进行计算时无异,因此蒙特卡洛分析模块独立于主程序作为外循环存在。为了提高计算效率,在算法实现上,蒙特卡洛模块利用了MATLAB中的Parallel Computing Toolbox(并行计算工具箱),采用支持并行计算的parfor语句进行编写,能够在多核的计算机或计算机集群上进行并行计算,以缩短计算时间。
本发明有益效果在于:
1.充分考虑了每次空降空投初始时刻的重量、高度、速度和空降空投过程中的风向、风速大小都具有的不确定性,使空降空投仿真计算更接近真实情况。
2.利用了MATLAB中的Parallel Computing Toolbox(并行计算工具箱),采用支持并行计算的parfor语句进行编写,能够在多核的计算机或计算机集群上进行并行计算,大大缩短了需要进行大量的循环计算所带来的计算时间。
附图说明
图1为蒙特卡洛分析模块实现示意图。
具体实施方式
下面结合附图对本发明作进一步说明。
本发明提供了一种空降空投并行蒙特卡洛仿真的实现方法,包括以下步骤:
1)给出初始条件的期望值,再给出概率分布规律,将空降空投中的风向和风速大小在一定范围内随机均匀分布,而空降空投初始时刻的重量、高度、速度服从正态分布,其中各随机变量概率分布密度如下:
风速大小:
式中:Vmin和Vmax分别为某一高度上最小和最大风速;
风速与初始航向夹角:
式中:当θ风=π时为逆风,θ风=-π时为顺风;
空降空投重量:
式中:σm为空降空投重量分布方差,mμ为期望空降空投重量;对于某型空投件,默认取mμ=150kg,σm=5kg。
空降空投高度:
式中:σh为空降空投高度分布方差,hμ为期望空降空投高度;对于某型空投件,默认取hμ=400m,σh=17m。
空降空投速度:
式中:σMa为空降空投初始马赫数分布方差,Maμ为期望空降空投马赫数;对于某型空投件,默认取Maμ=0.35,σMa=0.02。
2)蒙特卡洛分析模块将根据以上变量对应的概率分布情况随机生成大量初始条件,并导出到主程序的输入文件中供主程序进行计算;的
3)进行大量的循环计算,分析空投落点散布时取500~5000次,分析空中风向和风速,空投初始时刻的重量、高度、速度对计算结果的影响。
在进行蒙特卡洛仿真时,由于每一次计算的情况与之前设置固定参数进行计算时无异,因此蒙特卡洛分析模块独立于主程序作为外循环存在,这从图1中可以看出。同时,出于数理统计的精度要求考虑,蒙特卡洛分析需要进行大量的循环计算(一般的,分析空投落点散布时可取500~5000次),为了提高计算效率,在算法实现上,蒙特卡洛模块利用了MATLAB中的Parallel Computing Toolbox(并行计算工具箱),采用支持并行计算的parfor语句进行编写,能够在多核的计算机或计算机集群上进行并行计算,以缩短计算时间。该模块的实现原理如图1所示,该模块的实现方法为成熟的现有技术,由计算机现有程序实现全部计算过程,对于具体计算过程本发明中不再累述。
本发明具体应用途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进,这些改进也应视为本发明的保护范围。
Claims (2)
1.一种空降空投并行蒙特卡洛仿真的实现方法,其特征在于包括以下步骤:
1)给出初始条件的期望值,再给出概率分布规律,将空降空投中的风向和风速大小在一定范围内随机均匀分布,而空降空投初始时刻的重量、高度、速度服从正态分布,其中各随机变量概率分布密度如下:
风速大小:
式中:Vmin和Vmax分别为某一高度上最小和最大风速;
风速与初始航向夹角:
式中:当θ风=π时为逆风,θ风=-π时为顺风;
空降空投重量:
式中:σm为空降空投重量分布方差,mμ为期望空降空投重量;
空降空投高度:
式中:σh为空降空投高度分布方差,hμ为期望空降空投高度;
空降空投速度:
式中:σMa为空降空投初始马赫数分布方差,Maμ为期望空降空投马赫数;
2)蒙特卡洛分析模块将根据以上变量对应的概率分布情况随机生成初始条件,并导出到主程序的输入文件中供主程序进行计算;
3)进行循环计算,分析空中风向和风速,空投初始时刻的重量、高度、速度对计算结果的影响。
2.根据权利要求1所述的空降空投并行蒙特卡洛仿真的实现方法,其特征在于:步骤3)中所述的循环计算中,分析空投落点散布时取500~5000次。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510573559.8A CN105468808B (zh) | 2015-09-10 | 2015-09-10 | 空降空投并行蒙特卡洛仿真的实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510573559.8A CN105468808B (zh) | 2015-09-10 | 2015-09-10 | 空降空投并行蒙特卡洛仿真的实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105468808A CN105468808A (zh) | 2016-04-06 |
CN105468808B true CN105468808B (zh) | 2019-02-19 |
Family
ID=55606505
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510573559.8A Expired - Fee Related CN105468808B (zh) | 2015-09-10 | 2015-09-10 | 空降空投并行蒙特卡洛仿真的实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105468808B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110781579B (zh) * | 2019-09-23 | 2021-07-20 | 同济大学 | 一种空投物体下落过程仿真方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104166771A (zh) * | 2014-08-26 | 2014-11-26 | 中国人民解放军装甲兵工程学院 | 复杂环境下空投着陆仿真模拟与评价方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110024569A1 (en) * | 2009-05-27 | 2011-02-03 | Conax Florida Corporation | Attenuated landing system |
-
2015
- 2015-09-10 CN CN201510573559.8A patent/CN105468808B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104166771A (zh) * | 2014-08-26 | 2014-11-26 | 中国人民解放军装甲兵工程学院 | 复杂环境下空投着陆仿真模拟与评价方法 |
Non-Patent Citations (2)
Title |
---|
复杂条件下军用装备定点投放的建模与仿真;李薇等;《兵工自动化》;20110930;第30卷(第9期);全文 |
空降着陆散布概率模型及仿真;张武生等;《火力与指挥控制》;20130331;第38卷(第3期);全文 |
Also Published As
Publication number | Publication date |
---|---|
CN105468808A (zh) | 2016-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105843073B (zh) | 一种基于气动力不确定降阶的机翼结构气动弹性稳定性分析方法 | |
IN2015DN04035A (zh) | ||
Abd-Elaal et al. | An analytical model for simulating steady state flows of downburst | |
CN107886158A (zh) | 一种基于迭代局部搜索和随机惯性权重的蝙蝠优化算法 | |
Yang et al. | LES study of turbulent flow fields over hilly terrains—Comparisons of inflow turbulence generation methods and SGS models | |
Guichard | Assessment of an improved Random Flow Generation method to predict unsteady wind pressures on an isolated building using Large-Eddy Simulation | |
Zahir et al. | Variable-fidelity optimization with design space reduction | |
CN105468808B (zh) | 空降空投并行蒙特卡洛仿真的实现方法 | |
CN104182511A (zh) | 一种簇特征加权的模糊紧致散布聚类方法 | |
CN109033487A (zh) | 一种基于蒙特卡洛仿真的飞行器总体概率设计方法 | |
CN105490836B (zh) | 一种复杂网络可靠度的蒙特卡罗评估方法 | |
Lei et al. | Variable step euler method for real-time simulation | |
CN102819237B (zh) | 一种实体法铣削仿真过程中仿真目标毛坯的产生方法 | |
Turkiyyah et al. | Fast vortex methods for predicting wind-induced pressures on buildings | |
CN105955029B (zh) | 一种保鲁棒性的pid控制参数优化方法 | |
US20140343905A1 (en) | Vorticity-refinement based numerical method for simulating aircraft wing-tip vortex flows | |
CN106815075A (zh) | 建筑火灾数值仿真的区域分解优化方法 | |
CN102779288A (zh) | 一种基于场理论的本体分析方法 | |
US8170813B2 (en) | Determining effects of turbine blades on fluid motion | |
CN106934729A (zh) | 建筑检测鉴定方法及装置 | |
Wu et al. | The Monte Carlo calculation method of multiple integration | |
CN104036098A (zh) | 一种并行同步扰动随机近似的气动优化设计方法 | |
CN104657595A (zh) | 一种单颗粒曳力模型系数标定方法 | |
RU2611892C1 (ru) | Способ трехмерного моделирования заданного гидрогеологического объекта, реализуемый в вычислительной системе | |
CN109583007A (zh) | 一种火星进入飞行状态不确定性量化方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190219 Termination date: 20200910 |
|
CF01 | Termination of patent right due to non-payment of annual fee |