CN107341295B - 具有终端角度和速度约束的下压段弹道设计方法 - Google Patents
具有终端角度和速度约束的下压段弹道设计方法 Download PDFInfo
- Publication number
- CN107341295B CN107341295B CN201710456078.8A CN201710456078A CN107341295B CN 107341295 B CN107341295 B CN 107341295B CN 201710456078 A CN201710456078 A CN 201710456078A CN 107341295 B CN107341295 B CN 107341295B
- Authority
- CN
- China
- Prior art keywords
- angle
- trajectory
- push
- attack
- section
- 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
- 238000013461 design Methods 0.000 title claims abstract description 35
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims description 10
- XOFYZVNMUHMLCC-ZPOLXVRWSA-N prednisone Chemical compound O=C1C=C[C@]2(C)[C@H]3C(=O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 XOFYZVNMUHMLCC-ZPOLXVRWSA-N 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 230000000881 depressing effect Effects 0.000 claims description 2
- 230000000977 initiatory effect Effects 0.000 claims 2
- 230000007547 defect Effects 0.000 description 1
- 238000012407 engineering method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005381 potential energy Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001052 transient effect Effects 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/17—Mechanical parametric or variational design
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了具体涉及具有终端角度和速度约束的下压段弹道设计方法,包括确定下压起始点步骤、下压段攻角设计步骤、弹道计算步骤,首先设计一条基准滑翔弹道,从该弹道上确定出下压点起始时刻所对应的滑翔机动弹头所具有的机械能,定义此机械能为机械基能E0,每条弹道下压起始点的机械能都管控到所述机械基能E0附近,作为下压点起始点,并且采用满足终端角度约束的自适应设计方法来设计下压段攻角,设计出下压段弹道,最后计算弹道,结果均能满足精度要求。该方法简单新颖,弹道生成速度快,非常适合工程应用。
Description
技术领域
本发明属于航天飞行器总体设计技术领域,具体涉及具有终端角度和速度约束的下压段弹道设计方法。
背景技术
滑翔机动弹头为了能够对特定目标实施精确打击,一般对终端速度和角度都有一定的约束条件。目前国内外有关终端角度和速度约束的理论研究大多都是致力于具有终端角度和速度约束的制导方法。具有终端角度和速度约束的末端制导律方法研究这一新兴方向,目前仍处于探索阶段,而国内有关终端角度和速度约束的下压段弹道设计工程方法大多是采用迭代规划的设计方法,该方法不足之处是本身迭代速度不够快速,以及对硬件的要求也非常高,若想在同等时间内在飞行器上生成一条弹道,就要求硬件计算机计算速度更快。因此,亟待找到一种快速精准的弹道设计方法能够较好地应用到工程设计中,进而提高滑翔机动弹头对目标的打击精度。
发明内容
为解决工程设计中采用迭代规划的设计方法迭代速度不够快速以及对硬件的要求过高的技术问题,本发明提供了具体涉及具有终端角度和速度约束的下压段弹道设计方法,对于滑翔机动弹头,机动段采用射面内的滑翔式弹道,其特征在于,所述下压段弹道设计方法包括如下步骤:
确定下压起始点步骤,首先设计一条满足落速及落角要求的基准滑翔弹道,从该弹道上确定出下压点起始时刻所对应的滑翔机动弹头所具有的机械能E,定义此机械能E为机械基能E0;然后设计该滑翔机动弹头所有弹道时,每条弹道下压起始点的机械能都管控到所述机械基能E0附近;
下压段攻角设计步骤,下压段开始后,设计下压段攻角满足落角的约束;
弹道计算步骤,在标准条件下使用三自由度弹道模型计算弹道,标准条件包括:大气条件为国家标准大气、所有高度上无风、弹道计算的初值为滑翔初始点的速度和位置参数、符合瞬时平衡原理,弹道计算时程序角设计如下:
|α′xy|≤α′xy_max;
具体地,所述确定下压起始点步骤中所述每条弹道下压起始点的机械能都管控到所述机械基能E0附近,具体地:设定每条弹道的机械能E满足条件E0-ΔE≤E≤E0+ΔE时所对应的弹道上的点为下压起始点,ΔE为小量,不大于E0的0.1%。
优选地,所述下压段攻角设计步骤中所述设计下压段攻角满足终端角度约束,具体为下压段攻角采用满足终端角度约束的自适应设计方法,其表达式为α′xy﹦0‐(θd-θd_m),其中,θd为当地弹道倾角,θd_m为标准落角,通过将导弹当前时刻的当地弹道倾角θd与标准落角θd_m的差值负反馈到下压段攻角设计中,该下压攻角α′xy值通过自动调整来适应并满足落角的约束。
优选地,所述弹道计算步骤中所述的最大可用攻角α′xy_max根据下压段过载约束确定。
可选地,所述弹道计算步骤中所述的最大可用攻角α′xy_max根据下压段气动热流约束确定。
可选地,所述弹道计算步骤中所述的最大可用攻角α′xy_max根据下压段控制稳定约束确定。
本发明的具有终端角度和速度约束的下压段弹道设计方法,首先设计一条基准滑翔弹道,从该弹道上确定出下压点起始时刻所对应的滑翔机动弹头所具有的机械能,定义此机械能为机械基能E0,每条弹道下压起始点的机械能都管控到所述机械基能E0附近,作为下压点起始点,并且采用满足终端角度约束的自适应设计方法来设计下压段攻角,设计出下压段弹道,最后计算弹道,结果均能满足精度要求。该方法简单新颖,弹道生成速度快,非常适合工程应用。
附图说明
图1为本实施例的弹道设计方法设计的弹道的射程-速度曲线图;
图2为本实施例的弹道设计方法设计的弹道的射程-当地弹道倾角曲线。
具体实施方式
下面结合具体实施例对本发明作进一步说明,但不应该以此限制本发明的保护范围。
本发明提供的具体涉及具有终端角度和速度约束的下压段弹道设计方法,对于滑翔机动弹头,机动段采用射面内的滑翔式弹道,下压段弹道设计方法包括如下步骤:
第一步 确定下压起始点
首先设计一条满足落速及落角要求的基准滑翔弹道,本实施例中落速要求800m/s±10m/s,落角要求-80°±1°,从该弹道上确定出下压点开始时刻所对应的滑翔机动弹头具有的机械能E,定义此机械能E为机械基能E0,即动能Ev与势能Eh之和,本实施例中E0=m×1.1×160J,其中m为弹头质量,然后设计该滑翔机动弹头所有弹道时,每条弹道下压起始点的机械能E都管控到该机械基能E0附近,本实施例中设定每条弹道的机械能E满足条件E0-ΔE≤E≤E0+ΔE时所对应的弹道上的点为下压起始点,ΔE为小量,不大于E0的0.1%。
第二步 下压段攻角设计即弹道设计
下压段攻角设计时,设计下压段攻角满足终端角度约束,本实施例中具体为下压段攻角采用满足终端角度约束的自适应设计方法,其表达式为α′xy﹦0‐(θd-θd_m),其中,θd为当地弹道倾角,θd_m为标准落角,通过将导弹当前时刻的当地弹道倾角θd与标准落角θd_m的差值负反馈到下压段攻角设计中,该下压攻角α′xy值通过自动调整来适应并满足落角(即标准落角θd_m)的约束,本实施例中最大可用攻角α′xy_max根据下压段过载约束确定,本实施例中,下压段过载不超过12g,要求下压段最大可用攻角α′xy_max为18°。最大可用攻角α′xy_max也可以根据下压段气动热流约束确定,或者根据下压段控制稳定约束确定。
第三步 计算弹道
在标准条件下使用三自由度弹道模型计算弹道。标准条件包括:
a)大气条件为国家标准大气;
b)所有高度上无风;
c)弹道计算的初值为滑翔初始点的速度、位置参数;
d)符合瞬时平衡原理。
弹道计算时程序角设计如下:
|α′xy|≤α′xy_max;
具体给出一型滑翔机动弹头的弹道设计实例如下:
第一步:确定下压起始点。
首先设计一条射程为2500km满足落速及落角要求的基准滑翔弹道,本实施例中落速要求800m/s±10m/s,落角要求-80°±1°,从该弹道上确定出下压点开始时刻所对应的滑翔机动弹头具有的机械基能E0=590.7×106,ΔE=0.5907×106。
然后设计该滑翔机动弹头所有弹道(例如射程为2500km、2100km、1700km三条弹道)时,每条弹道下压起始点的机械能E都管控到该机械基能E0附近,即E0-ΔE≤E≤E0+ΔE。
第二步 下压段攻角设计即弹道设计。
下压段弹道设计为α′xy=0-(θd-θd_m)=0-(θd-(-80°))。
第三步 计算弹道。
结果如下:
弹道一:射程为2500km,落速802.5m/s,落角-80.4°;
弹道二:射程为2200km,落速799.5m/s,落角-80.4°;
弹道三:射程为1700km,落速798.1m/s,落角-80.4°。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.具有终端角度和速度约束的下压段弹道设计方法,对于滑翔机动弹头,机动段采用射面内的滑翔式弹道,其特征在于,所述下压段弹道设计方法包括如下步骤:
确定下压起始点步骤,首先设计一条满足落速及落角要求的基准滑翔弹道,从该弹道上确定出下压点起始时刻所对应的滑翔机动弹头所具有的机械能E,定义此机械能E为机械基能E0;然后设计该滑翔机动弹头所有弹道时,每条弹道下压起始点的机械能都管控到所述机械基能E0附近;
下压段攻角设计步骤,下压段开始后,设计下压段攻角满足落角的约束;
弹道计算步骤,在标准条件下使用三自由度弹道模型计算弹道,标准条件包括:大气条件为国家标准大气、所有高度上无风、弹道计算的初值为滑翔初始点的速度和位置参数、符合瞬时平衡原理,弹道计算时程序角设计如下:
|α′xy|≤α′xy_max;
2.根据权利要求1所述的下压段弹道设计方法,其特征在于,所述确定下压起始点步骤中所述每条弹道下压起始点的机械能都管控到所述机械基能E0附近,具体地:设定每条弹道的机械能E满足条件E0-ΔE≤E≤E0+ΔE时所对应的弹道上的点为下压起始点,ΔE为小量,不大于E0的0.1%。
3.根据权利要求1或2所述的下压段弹道设计方法,其特征在于,所述下压段攻角设计步骤中所述设计下压段攻角满足终端角度约束,具体为下压段攻角采用满足终端角度约束的自适应设计方法,其表达式为α′xy﹦0-(θd-θd_m),其中,θd为当地弹道倾角,θd_m为标准落角,通过将导弹当前时刻的当地弹道倾角θd与标准落角θd_m的差值负反馈到下压段攻角设计中,该下压攻角α′xy值通过自动调整来适应并满足落角的约束。
4.根据权利要求3所述的下压段弹道设计方法,其特征在于所述弹道计算步骤中所述的最大可用攻角α′xy_max根据下压段过载约束确定。
5.根据权利要求2所述的下压段弹道设计方法,其特征在于所述弹道计算步骤中所述的最大可用攻角α′xy_max根据下压段气动热流约束确定。
6.根据权利要求2所述的下压段弹道设计方法,其特征在于所述弹道计算步骤中所述的最大可用攻角α′xy_max根据下压段控制稳定约束确定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710456078.8A CN107341295B (zh) | 2017-06-16 | 2017-06-16 | 具有终端角度和速度约束的下压段弹道设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710456078.8A CN107341295B (zh) | 2017-06-16 | 2017-06-16 | 具有终端角度和速度约束的下压段弹道设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107341295A CN107341295A (zh) | 2017-11-10 |
CN107341295B true CN107341295B (zh) | 2020-09-01 |
Family
ID=60219980
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710456078.8A Active CN107341295B (zh) | 2017-06-16 | 2017-06-16 | 具有终端角度和速度约束的下压段弹道设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107341295B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109343563B (zh) * | 2018-10-15 | 2020-06-05 | 北京理工大学 | 考虑舵机失效及落角约束的飞行器制导系统及方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8130137B1 (en) * | 2005-07-26 | 2012-03-06 | Lockheed Martin Corporation | Template updated boost algorithm |
CN103838914A (zh) * | 2013-12-30 | 2014-06-04 | 北京航空航天大学 | 一种高超声速飞行器滑翔段弹道解析求解方法 |
CN104978489A (zh) * | 2015-07-02 | 2015-10-14 | 北京航天自动控制研究所 | 一种适用于滑翔飞行器的最小铰链力矩下压弹道计算方法 |
-
2017
- 2017-06-16 CN CN201710456078.8A patent/CN107341295B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8130137B1 (en) * | 2005-07-26 | 2012-03-06 | Lockheed Martin Corporation | Template updated boost algorithm |
CN103838914A (zh) * | 2013-12-30 | 2014-06-04 | 北京航空航天大学 | 一种高超声速飞行器滑翔段弹道解析求解方法 |
CN104978489A (zh) * | 2015-07-02 | 2015-10-14 | 北京航天自动控制研究所 | 一种适用于滑翔飞行器的最小铰链力矩下压弹道计算方法 |
Non-Patent Citations (3)
Title |
---|
Modeling a Low-Energy Ballistic Lunar Transfer Using Dynamical Systems Theory;Parker J S 等;《Journal of Spacecraft and Rockets》;20081231;第45卷(第6期);第1269-1281页 * |
一种滑翔飞行器弹道的参数优化方法研究;王方鹏 等;《航空制造技术》;20151106;第104-107页 * |
助推—滑翔式飞行器弹道设计与制导技术研究;刘欣;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20141015;第C032-7页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107341295A (zh) | 2017-11-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110471450B (zh) | 在高度速度剖面内直接规划再入轨迹的方法 | |
CN111473696B (zh) | 一种基于落点估计的制导火箭垂直攻击制导方法 | |
CN111306998A (zh) | 一种参数摄动自适应的制导火箭弹垂直攻击制导方法 | |
CN103983143B (zh) | 包含速度过程约束和多终端约束的空地导弹投放下滑段制导方法 | |
CN110717245B (zh) | 基于落角落速约束的拟滑翔弹道的设计方法 | |
CN104567545B (zh) | Rlv大气层内主动段的制导方法 | |
CN103307938A (zh) | 一种旋转弹气动参数获取方法 | |
CN107341295B (zh) | 具有终端角度和速度约束的下压段弹道设计方法 | |
CN114200828A (zh) | 一种超音速大机动靶标持续大过载防失速方法 | |
CN115657458A (zh) | 一种基于能量匹配的飞行器爬升轨迹控制方法 | |
CN111679687B (zh) | 一种带有落角约束的导引控制一体化方法 | |
CN112629339A (zh) | 一种基于直接法的火箭软着陆轨迹规划方法 | |
Ge et al. | Three-dimensional trajectory tracking guidance law based on linear quadratic regulator | |
CN116301028B (zh) | 基于吸气式高超声速平台的多约束在线飞行轨迹规划中段导引方法 | |
CN107145761A (zh) | 带动力滑翔导弹多样式弹道规划方法 | |
CN111737815B (zh) | 一种基于无量纲形式的飞行系统落点误差修正方法 | |
WO2011127338A2 (en) | Predictive roll capture | |
CN111026153A (zh) | 一种提高助推滑翔飞行器飞行距离的导引方法及导引装置 | |
CN106643298A (zh) | 一种基于预置命中点的大气层内反导拦截器中制导方法 | |
Zhang et al. | Trajectory optimization of a deflectable nose missile | |
CN118244804A (zh) | 一种采用附加侧滑角在线快速迭代速度控制方法 | |
CN111486755B (zh) | 一种可变气动外形制导装备的落速控制方法 | |
CN116361926B (zh) | 一种弹道导弹滑翔增程段初始机械能闭环调整方法 | |
CN114754628B (zh) | 基于落点预测和虚拟跟踪的飞行体弹道控制方法 | |
Xu et al. | Ballistic characteristics of rocket projectile with deflection nose |
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 |