CN116305590A - 基于灵敏度因子分析的超远程制导火箭射程能力评估方法 - Google Patents
基于灵敏度因子分析的超远程制导火箭射程能力评估方法 Download PDFInfo
- Publication number
- CN116305590A CN116305590A CN202310572536.XA CN202310572536A CN116305590A CN 116305590 A CN116305590 A CN 116305590A CN 202310572536 A CN202310572536 A CN 202310572536A CN 116305590 A CN116305590 A CN 116305590A
- Authority
- CN
- China
- Prior art keywords
- range
- rocket
- guided rocket
- glide
- mechanical energy
- 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
- 230000035945 sensitivity Effects 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000000556 factor analysis Methods 0.000 title claims abstract description 16
- 238000004458 analytical method Methods 0.000 claims abstract description 24
- 238000011156 evaluation Methods 0.000 claims description 22
- 239000003380 propellant Substances 0.000 claims description 19
- 238000002347 injection Methods 0.000 claims description 10
- 239000007924 injection Substances 0.000 claims description 10
- 239000002245 particle Substances 0.000 claims description 8
- 238000012512 characterization method Methods 0.000 claims description 5
- 230000010354 integration Effects 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 abstract description 8
- 230000002596 correlated effect Effects 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 235000015842 Hesperis Nutrition 0.000 description 1
- 235000012633 Iberis amara Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
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
-
- 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)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Testing Of Engines (AREA)
Abstract
本发明公开了一种基于灵敏度因子分析的超远程制导火箭射程能力评估方法,首先建立制导火箭滑翔段射程与机械能之间的微分表达式,然后引入准平衡滑翔条件和发动机内弹道特性,建立关于发动机内弹道特性的制导火箭滑翔段射程解析评估公式,紧接着确定制导火箭射程影响因素并对其灵敏度因子进行计算,最后基于灵敏度因子分析各因素对制导火箭射程的影响关系。该基于灵敏度因子分析的超远程制导火箭射程能力评估方法具有形式简单,通用性强等特性,未来可广泛应用于各类超远程制导弹药的方案设计及优化过程中。
Description
技术领域
本发明属于火箭技术领域,具体涉及一种超远程制导火箭射程能力评估方法。
背景技术
超远程制导火箭采用全程大气层内滑翔弹道方案,能够利用气动力大幅改变弹道飞行轨迹,具有机动能力强、隐蔽性好等特点,得到了世界各军事强国的广泛关注。超远程制导火箭总体设计过程是一个关于多学科、多变量、多约束的复杂优化问题的求解过程。在采用传统的数值优化方法时,射程能力作为其中重要的战技指标之一,通常被选为目标函数,整个优化过程涉及质量特性、气动特性、力-热载荷特性、推力特性等多种复杂因素的联合寻优,这导致优化速度慢、效率低,难以快速分析并筛选出对制导火箭射程影响较大的因素。
发明内容
为了克服现有技术的不足,本发明提供了一种基于灵敏度因子分析的超远程制导火箭射程能力评估方法,在纵向平面质点动力学模型基础上,结合准平衡滑翔条件和内弹道特性,建立了关于发动机内弹道特性的制导火箭滑翔段射程解析评估公式,确定影响制导火箭射程的主要因素,进而求取每个因素的射程灵敏度因子并进行分析。
本发明针对超远程制导火箭射程能力评估过程中分析变量多、学科耦合复杂、离散数值分析效率低的问题,首先建立制导火箭滑翔段射程与机械能之间的微分表达式,然后引入准平衡滑翔条件和发动机内弹道特性,建立关于发动机内弹道特性的制导火箭滑翔段射程解析评估公式,紧接着确定制导火箭射程影响因素并对其灵敏度因子进行计算,最后基于灵敏度因子分析各因素对制导火箭射程的影响关系。该基于灵敏度因子分析的超远程制导火箭射程能力评估方法具有形式简单,通用性强等特性,未来可广泛应用于各类超远程制导弹药的方案设计及优化过程中。
本发明解决其技术问题所采用的技术方案包括如下步骤:
步骤1:建立制导火箭滑翔段射程与机械能之间的微分关系;
首先基于纵向平面制导火箭质点动力学模型,建立制导火箭被动段机械能表达公式,并求取机械能的时间导数,进而结合超远程制导火箭滑翔段射程关于时间的导数,得到制导火箭滑翔段射程关于机械能的微分表达式;
步骤2:基于机械能表征的制导火箭滑翔段射程的解析评估公式;
基于准平衡滑翔条件,对制导火箭射程关于机械能的微分表达式进行化简,然后对其进行积分,得到基于机械能表征的制导火箭滑翔段射程的解析评估公式;
步骤3:建立关于发动机内弹道特性的制导火箭滑翔段射程解析评估公式;
根据齐奥尔科夫斯基公式,获得制导火箭推进剂喷射后的理想机械能,在此基础上,通过对理想机械能进行工程修正,从而近似获得主动段结束后的实际机械能,并将其带入到步骤2中得到的制导火箭滑翔段射程评估公式中,得到关于发动机内弹道特性的制导火箭滑翔段射程解析评估公式;
步骤4:影响因素确定及射程灵敏度因子计算;
根据步骤3中建立的制导火箭滑翔段射程解析评估公式,分析得到关于超远程制导火箭滑翔段射程能力的影响因素,进而求取各个因素的射程灵敏度因子,即步骤3中的射程解析评估公式关于各个因素的偏导数;
步骤5:各影响因素与制导火箭滑翔段射程的相关性及优先级分析;
对步骤4中得到各因素射程灵敏度因子进行分析,得到各影响因素与超远程制导火箭滑翔段射程的相关关系;同时,针对具体项目,将初步方案数据带入灵敏度因子公式,得到各因素对射程的具体灵敏度数值,并进行优先级排序,筛选出对制导火箭射程影响敏感的因素。
进一步地,所述步骤1具体为:
纵向平面质点动力学模型如式(1)所示:
机械能定义如下:
根据公式(1)和(3)得到制导火箭滑翔段射程关于机械能的微分表达式如下
进一步地,所述步骤2具体为:
将公式(5)带入公式(4)可得
对公式(7)进行积分得到:
进一步地,所述步骤3具体为:
进一步地,所述步骤4具体为:
进一步地,所述步骤5具体为:
由公式(12)-(16)可知,升阻比、发动机比冲、满载质量与射程呈正相关关系,空载质量和落速与射程呈负相关关系;
为提升滑翔段射程从以下角度进行优化:
a) 增加升阻比;
将数据带入公式(12)-(16)中,从而计算出每个因素的射程灵敏度因子数值,并进行优先级排序,筛选出对制导火箭射程影响敏感的因素。
本发明的有益效果如下:
本发明提出的基于灵敏度因子分析的制导火箭射程能力评估方法能够快速分析出各影响因素与射程之间的相关性,并筛选出对射程影响最为敏感的因素,可加快制导火箭总体方案优化效率。该方法具有广阔的应用前景。
附图说明
图1为本发明方法流程图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
超远程制导火箭普遍采用全大气层内机动滑翔弹道,具有弹道特性复杂、约束限制多样、影响因素众多等特征。传统的数值积分方法在评估射程能力时存在时间长、工况多、迭代慢等问题,导致难以快速分析得到众多影响因素对制导火箭射程的影响程度。
为解决上述问题,本发明提出一种基于灵敏度因子分析的超远程制导火箭射程能力评估方法。考虑到超远程制导火箭主动段持续时间较短,同时主动段普遍采用程序角控制思路,因此主动段飞行射程较为确定,且占总射程比例较低,所以超远程制导火箭射程能力主要由被动段滑翔距离决定。针对上述分析,本发明根据质点动力学方程,首先建立制导火箭滑翔段射程与机械能之间的微分关系式;其次,在准平衡滑翔条件下,建立制导火箭滑翔段射程关于机械能的解析评估公式;然后,引入发动机内弹道特性,建立关于发动机内弹道特性的滑翔段射程解析评估公式;接着,求取各个影响因素的射程灵敏度因子;最后,对各个因素的灵敏度因子进行分析。
如图1所示,一种基于灵敏度因子分析的超远程制导火箭射程能力评估方法,具体步骤如下:
步骤一,建立制导火箭滑翔段射程与机械能之间的微分关系;
首先基于纵向平面制导火箭质点动力学模型,建立制导火箭被动段机械能表达公式,并求取机械能的时间导数,进而结合超远程制导火箭滑翔段射程关于时间的导数,可得到制导火箭滑翔段射程关于机械能的微分表达式。
针对纵向平面质点动力学模型,如下所示
机械能定义如下:
根据公式(1)和(3)得到制导火箭滑翔段射程关于机械能的微分表达式如下
步骤二,基于机械能表征的制导火箭滑翔段射程的解析评估公式;
将公式(5)带入公式(4)可得
对公式(7)进行积分得到:
步骤三,建立关于发动机内弹道特性的制导火箭滑翔段射程解析评估公式;
根据齐奥尔科夫斯基公式,可以获得制导火箭推进剂喷射后的理想机械能,再此基础上,考虑到制导火箭主动段通常按给定方案角飞行,规律较为一致,因此通过对理想机械能进行工程修正,从而近似获得主动段结束后的实际机械能,并将其带入到步骤二中得到的制导火箭滑翔段射程评估公式中,可得关于发动机内弹道特性的制导火箭滑翔段射程解析评估公式。
在方案论证阶段,在分析制导火箭的最大射程时,通常假设发射点和目标点高度为0m,因此由公式(9)可得,推进剂喷射完时制导火箭的理想机械能为。上述分析中未考虑主动段大气阻力消耗,因此比主动段结束后的实际机械能偏大。由于主动段通常按给定方案角飞行,规律较为一致,因此根据气动特性和发动机工作特性,可选取一个工程系数/>近似气动阻力耗散,从而近似得到推进剂喷射完的实际机械能,即
步骤四,影响因素确定及射程灵敏度因子计算;
根据步骤三中建立的制导火箭滑翔段射程解析评估公式,可分析得到关于超远程制导火箭滑翔段射程能力的主要影响因素,进而求取各个因素的射程灵敏度因子,即步骤三中的射程解析评估公式关于各个因素的偏导数。
步骤五,各影响因素与制导火箭滑翔段射程的相关性及优先级分析;
对步骤四中得到各因素射程灵敏度因子进行分析,可得到各影响因素与超远程制导火箭滑翔段射程的相关关系;同时,针对具体项目,可将初步方案数据带入灵敏度因子公式,得到各因素对射程的具体灵敏度数值,并进行优先级排序,快速筛选出对制导火箭射程影响较为敏感的因素。
由公式(12)-(16)可知,升阻比、发动机比冲、满载质量与射程呈正相关关系,空载质量和落速与射程呈负相关关系。通常,落速作为战技指标的关键项,必须要处于指定范围内,因此为提升滑翔段射程可以从以下角度进行优化:
a) 增加升阻比,需要对制导火箭的气动外形进行优化;
若要分析上述各个因素对射程的灵敏度大小排序,可将初步方案数据带入公式(12)-(16)中,从而计算出每个因素的灵敏度因子数值,并进行优先级排序,快速筛选出对制导火箭射程影响较为敏感的因素。
Claims (7)
1.一种基于灵敏度因子分析的超远程制导火箭射程能力评估方法,其特征在于,包括以下步骤:
步骤1:建立制导火箭滑翔段射程与机械能之间的微分关系;
首先基于纵向平面制导火箭质点动力学模型,建立制导火箭被动段机械能表达公式,并求取机械能的时间导数,进而结合超远程制导火箭滑翔段射程关于时间的导数,得到制导火箭滑翔段射程关于机械能的微分表达式;
步骤2:基于机械能表征的制导火箭滑翔段射程的解析评估公式;
基于准平衡滑翔条件,对制导火箭射程关于机械能的微分表达式进行化简,然后对其进行积分,得到基于机械能表征的制导火箭滑翔段射程的解析评估公式;
步骤3:建立关于发动机内弹道特性的制导火箭滑翔段射程解析评估公式;
根据齐奥尔科夫斯基公式,获得制导火箭推进剂喷射后的理想机械能,在此基础上,通过对理想机械能进行工程修正,从而近似获得主动段结束后的实际机械能,并将其带入到步骤2中得到的制导火箭滑翔段射程评估公式中,得到关于发动机内弹道特性的制导火箭滑翔段射程解析评估公式;
步骤4:影响因素确定及射程灵敏度因子计算;
根据步骤3中建立的制导火箭滑翔段射程解析评估公式,分析得到关于超远程制导火箭滑翔段射程能力的影响因素,进而求取各个因素的射程灵敏度因子,即步骤3中的射程解析评估公式关于各个因素的偏导数;
步骤5:各影响因素与制导火箭滑翔段射程的相关性及优先级分析;
对步骤4中得到各因素射程灵敏度因子进行分析,得到各影响因素与超远程制导火箭滑翔段射程的相关关系;同时,针对具体项目,将初步方案数据带入灵敏度因子公式,得到各因素对射程的具体灵敏度数值,并进行优先级排序,筛选出对制导火箭射程影响敏感的因素。
4.根据权利要求3所述的一种基于灵敏度因子分析的超远程制导火箭射程能力评估方法,其特征在于,所述步骤3具体为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310572536.XA CN116305590B (zh) | 2023-05-22 | 2023-05-22 | 基于灵敏度因子分析的超远程制导火箭射程能力评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310572536.XA CN116305590B (zh) | 2023-05-22 | 2023-05-22 | 基于灵敏度因子分析的超远程制导火箭射程能力评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116305590A true CN116305590A (zh) | 2023-06-23 |
CN116305590B CN116305590B (zh) | 2023-09-01 |
Family
ID=86827283
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310572536.XA Active CN116305590B (zh) | 2023-05-22 | 2023-05-22 | 基于灵敏度因子分析的超远程制导火箭射程能力评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116305590B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104656450A (zh) * | 2015-01-20 | 2015-05-27 | 北京航空航天大学 | 一种高超声速飞行器平稳滑翔再入弹道设计方法 |
CN106029185A (zh) * | 2014-05-06 | 2016-10-12 | 深圳市大疆创新科技有限公司 | 检测表面上的投射物击中的设备、系统和方法 |
CN109190248A (zh) * | 2018-09-03 | 2019-01-11 | 中国运载火箭技术研究院 | 一种用于滑翔飞行器的滑翔射程解析方法及解析系统 |
CN112329137A (zh) * | 2020-10-27 | 2021-02-05 | 深圳三零三防务科技有限公司 | 一种基于平衡飞行理论的运载火箭在线入轨能力评估方法 |
EP4068042A1 (en) * | 2021-03-24 | 2022-10-05 | Beijing Interstellar Glory Space Technology Co., Ltd. | Recovery control method and apparatus for rocket, electronic device and storage medium |
-
2023
- 2023-05-22 CN CN202310572536.XA patent/CN116305590B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106029185A (zh) * | 2014-05-06 | 2016-10-12 | 深圳市大疆创新科技有限公司 | 检测表面上的投射物击中的设备、系统和方法 |
CN104656450A (zh) * | 2015-01-20 | 2015-05-27 | 北京航空航天大学 | 一种高超声速飞行器平稳滑翔再入弹道设计方法 |
CN109190248A (zh) * | 2018-09-03 | 2019-01-11 | 中国运载火箭技术研究院 | 一种用于滑翔飞行器的滑翔射程解析方法及解析系统 |
CN112329137A (zh) * | 2020-10-27 | 2021-02-05 | 深圳三零三防务科技有限公司 | 一种基于平衡飞行理论的运载火箭在线入轨能力评估方法 |
EP4068042A1 (en) * | 2021-03-24 | 2022-10-05 | Beijing Interstellar Glory Space Technology Co., Ltd. | Recovery control method and apparatus for rocket, electronic device and storage medium |
Non-Patent Citations (3)
Title |
---|
CHUANG ZHOU 等: "Analysis of dynamic characteristics and sensitivity of hydrogen-oxygen expansion cycle rocket engine system", ACTA ASTRONAUTICA * |
张可菊 等: "弹箭弹道参数对气动参数灵敏度分析", 沈阳理工大学学报 * |
杨云刚 等: "基于最小二乘法的弹体参数辨识工程算法", 弹箭与制导学报 * |
Also Published As
Publication number | Publication date |
---|---|
CN116305590B (zh) | 2023-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110083971B (zh) | 一种基于作战推演的自爆式无人机集群作战兵力分配方法 | |
CN105956281A (zh) | 固体火箭发动机装药设计方法 | |
CN105608251B (zh) | 直升机火控系统精度敏感性分析的BNSobol法 | |
CN108830023B (zh) | 一种炮射导弹增速发动机点火过程的预测方法 | |
CN111506113A (zh) | 飞行器制导指令计算方法、侧滑角计算方法及制导方法 | |
CN111176315B (zh) | 一种基于l1自适应控制的变体巡飞弹控制方法 | |
CN117892559A (zh) | 一种超远程制导火箭总体协调多学科分级优化方法 | |
CN117892558B (zh) | 一种超远程制导火箭多学科动态优化模型构建方法 | |
CN118094779B (zh) | 基于蚁群算法的制导火箭大空域滑翔增程弹道优化方法 | |
CN115828412A (zh) | 考虑外弹道约束的单室双推固体火箭发动机配比参数设计方法 | |
CN116305590B (zh) | 基于灵敏度因子分析的超远程制导火箭射程能力评估方法 | |
CN109325288B (zh) | 一种基于不确定性优化的固体运载器总体参数确定方法及系统 | |
CN106643298A (zh) | 一种基于预置命中点的大气层内反导拦截器中制导方法 | |
CN112286059B (zh) | 一种具有攻击角和视场角约束的二阶前置角重塑制导方法 | |
Xu et al. | Artillery structural dynamic responses uncertain optimization based on robust Nash game method | |
CN114819055A (zh) | 一种基于lstm网络的导弹拦截点预测方法 | |
SAVASTRE et al. | Theoretical evaluation of drag coefficient for different geometric configurations of ballistic caps for an experimental 30× 165 mm AP-T projectile | |
CN118133585B (zh) | 一种考虑多维力热约束的远程火箭弹道优化方法 | |
Cai et al. | Research on the Internal Ballistic Performance Prediction of Solid Rocket Motor Based on Monte Carlo | |
Zhou et al. | An integrate VWS-TOPSIS-PSO approach and application in overall design parameter optimization of Ballistic Missile | |
CN118133634B (zh) | 一种固体火箭发动机装药结构完整性评估方法 | |
CN109204828A (zh) | 一种用于步枪后坐缓冲的小型无人机低后坐前冲机构 | |
Soham et al. | Trajectory simulations by the numerical solution of the point-mass equations of motion for 7.62 mm/. 308” rifle bullets | |
Aboelkhair et al. | Effect of Projectile Shape on the Power of Fire in Personal Defense Hand-Held Weapons | |
CN117932793B (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 |