CN116576737A - 减小高度散布的垂直发射远程导弹转弯程序角生成方法 - Google Patents
减小高度散布的垂直发射远程导弹转弯程序角生成方法 Download PDFInfo
- Publication number
- CN116576737A CN116576737A CN202310600762.4A CN202310600762A CN116576737A CN 116576737 A CN116576737 A CN 116576737A CN 202310600762 A CN202310600762 A CN 202310600762A CN 116576737 A CN116576737 A CN 116576737A
- Authority
- CN
- China
- Prior art keywords
- missile
- ballistic
- interpolation
- angle
- reference axis
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 241000764238 Isis Species 0.000 claims 1
- 239000006185 dispersion Substances 0.000 abstract description 14
- 230000009194 climbing Effects 0.000 abstract description 5
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000005457 optimization Methods 0.000 description 3
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 230000003121 nonmonotonic effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F42—AMMUNITION; BLASTING
- F42B—EXPLOSIVE CHARGES, e.g. FOR BLASTING, FIREWORKS, AMMUNITION
- F42B15/00—Self-propelled projectiles or missiles, e.g. rockets; Guided missiles
- F42B15/01—Arrangements thereon for guidance or control
-
- 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)
- General Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Aviation & Aerospace Engineering (AREA)
- Combustion & Propulsion (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种减小高度散布的垂直发射远程导弹转弯程序角生成方法，以垂直发射远程导弹为背景，对初始爬升段的弹道倾角指令插值方法进行研究，定义新的弹道倾角指令参考轴，并用射前优化得到的倾角指令与参考轴生成插值表，利用导弹位置在参考轴上的投影作为插值变量进行一维插值得到弹道倾角指令，计算结果表明本发明能够非常有效得解决时间插值弹道倾角指令带来的弹道高度散布问题。
Description
技术领域
本发明属于导弹技术领域，具体涉及一种垂直发射远程导弹转弯程序角生成方法。
背景技术
垂直发射远程导弹发射后，飞控装置需要在爬升段对导弹的弹道倾角进行控制，保证导弹按照程序弹道倾角指令进行转弯，同时保证弹道最高点散布较小。获得弹道倾角指令传统方法是利用弹上时间插射前优化的方案角指令得到。通过仿真分析，在不同拉偏条件下，弹道最高点散布较大，这对后续的高度控制、射程和落速等均非常不利。
发明内容
为了克服现有技术的不足，本发明提供了一种减小高度散布的垂直发射远程导弹转弯程序角生成方法，以垂直发射远程导弹为背景，对初始爬升段的弹道倾角指令插值方法进行研究，定义新的弹道倾角指令参考轴，并用射前优化得到的倾角指令与参考轴生成插值表，利用导弹位置在参考轴上的投影作为插值变量进行一维插值得到弹道倾角指令，计算结果表明本发明能够非常有效得解决时间插值弹道倾角指令带来的弹道高度散布问题。
本发明解决其技术问题所采用的技术方案包括如下步骤：
步骤1：定义弹道倾角参考轴与发射系XOY在同一平面且与X轴夹角为，其中X为前向位置，Y为纵向位置；
步骤2：利用导弹位置在弹道倾角参考轴上的投影插值得到弹道倾角；
步骤3：在发射前根据倾角指令与参考轴角度生成插值表，在发射中利用导弹位置在参考轴上的投影作为插值变量进行一维插值得到弹道倾角指令，然后控制弹道倾角跟随弹道倾角指令。
优选地，所述°，则点火时刻和弹道最高点时导弹速度与参考轴夹角均为45°。
优选地，所述为随时间变化的时变值/>。
本发明的有益效果如下：
本发明方法由于采用了位置在参考轴上的投影作为插值变量，可以有效屏蔽掉由于温度不同带来的发动机推力在时间轴上散布的影响，较好地解决了弹道高度散布较大问题。
附图说明
图1 是本发明实施例时间插值弹道倾角指令变化曲线。
图2 是本发明实施例时间插值对应不同温度弹道高度变化曲线。
图3 是本发明实施例导弹位置x插值弹道倾角及局部放大曲线，（a）弹道倾角放大曲线，（b）局部放大曲线。
图4 是本发明实施例导弹位置y插值弹道倾角及局部放大曲线，（a）弹道倾角放大曲线，（b）局部放大曲线。
图5 是本发明实施例弹道倾角指令随xy45°插值曲线。
图6 是本发明实施例不同温度发动机推力曲线。
图7 是本发明实施例不同温度弹道倾角指令曲线。
图8 是本发明实施例弹道高度随时间变化曲线。
图9 是本发明实施例弹道高度随导弹位置x变化曲线。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
某垂直发射远程导弹初始转弯段采用弹上时间插值弹道倾角控制转弯，由于发动机工作温度不同推力曲线差异较大等原因会带来爬升段高度散布较大，采用时间插值得到弹道倾角指令表时最高点散布约15km，这对后续的高度控制、射程和落速等均有较大的影响。指令曲线如图1所示，高度散布如图2所示。本发明通过优化弹道倾角插值方法，解决弹道高度散布问题。
为解决上述问题，本发明提供一种弹道倾角指令插值方法，保证插值变量单调的同时，能够解决传统时间插值方法带来的高度散布问题。
为解决高度散布问题，首先联想到使用位置插值，保证在位置到规定值的时候，倾角指令为0，则可大幅度减小高度散布。尝试弹道倾角随导弹位置x或者位置y插值得到，绘制随x插值得到的弹道倾角曲线和局部放大如图3所示，从图中可以看出垂直发射初始段导弹位置可能会不单调且垂直发射初始段弹道倾角指令随x变化斜率过大，这都不利于得到平滑的弹道倾角指令。绘制随y插值得到的弹道倾角曲线和局部放大如图4所示。从图中可以看出弹道最高点也是导弹倾角指令随y变化斜率由正变负的点，也就是为保证插值表的单调只能取到最高点，这样微小的高度正散布也会带来外插问题。
一种减小高度散布的垂直发射远程导弹转弯程序角生成方法，主要包括以下内容：
如下图所示，定义弹道倾角参考轴，与发射系XOY（X为前向位置，Y为纵向位置）在同一平面且与X轴夹角为固定值或者时变值/>，下面以固定值为例进行说明。
利用导弹位置在参考轴上的投影插值得到弹道倾角。
考虑到垂直发射情况取°，这样可以保证点火时刻和弹道最高点导弹速度与参考轴夹角均为45°，绘制弹道倾角指令变化曲线如图5所示，从图中可以看出，弹道倾角过零点对应弹道最高点，重新定义参考轴之后，很好解决了初始段弹道倾角随导弹位置x或y变化斜率过大甚至不单调的问题。
由于采用了位置在某参考轴上的投影作为插值变量，可以有效屏蔽掉由于温度不同带来的发动机推力在时间轴上散布的影响。
具体实施例：
1. 系统
垂直发射远程导弹发射后弹载计算机根据装定的爬升方案角插值表进行插值得到弹道倾角指令，通过弹道倾角指令与弹道倾角偏差计算得到舵偏角指令。电动舵机执行弹载计算机发出的舵偏角指令，产生控制力矩，控制导弹转向目标方向。为增强弹道倾角指令的适应性，避免因为发动机推力散布等原因带来的弹道高度散布过大问题，需要对弹道倾角指令生成方法进行设计。
2. 方法
本发明改变传统的时间插值得到弹道倾角指令方法，定义新的弹道倾角指令参考轴，并用射前优化得到的倾角指令与参考轴生成插值表，利用导弹位置在参考轴上的投影作为插值变量进行一维插值得到弹道倾角指令，然后控制弹道倾角跟随弹道倾角指令。较好地解决了弹道高度散布较大问题。
采用相同的发射前优化结果，只改变插值方法，对比不同弹道倾角插值方法对弹道高度散布的影响。
在不同温度（-10℃、20℃、45℃）条件下发动机推力曲线如图6所示，可以看出发动机工作时间和推力大小差异较大。采用时间插值正如上面所述，带来最高点高度散布约15km。
利用本发明方法插值得到弹道倾角，考虑到垂直发射情况取°，这样可以保证点火时刻和弹道最高点导弹速度与参考轴夹角均为45°。
绘制不同温度下弹道倾角指令变化曲线如图7所示，从图中可以看出弹道倾角指令会有一定的差异以适应温度变化。绘制弹道高度随时间和导弹位置x变化曲线如图8和图9所示，可以看出不同温度下使用新方法虽然达到最高点的时间有所不同，但弹道高度散布极小，达到了预期的目的。
Claims (3)
1.一种减小高度散布的垂直发射远程导弹转弯程序角生成方法，其特征在于，包括如下步骤：
步骤1：定义弹道倾角参考轴与发射系XOY在同一平面且与X轴夹角为，其中X为前向位置，Y为纵向位置；
步骤2：利用导弹位置在弹道倾角参考轴上的投影插值得到弹道倾角；
步骤3：在发射前根据倾角指令与参考轴角度生成插值表，在发射中利用导弹位置在参考轴上的投影作为插值变量进行一维插值得到弹道倾角指令，然后控制弹道倾角跟随弹道倾角指令。
2.根据权利要求1所述的一种减小高度散布的垂直发射远程导弹转弯程序角生成方法，其特征在于，所述°，则点火时刻和弹道最高点时导弹速度与参考轴夹角均为45°。
3.根据权利要求1所述的一种减小高度散布的垂直发射远程导弹转弯程序角生成方法，其特征在于，所述为随时间变化的时变值/>。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310600762.4A CN116576737A (zh) | 2023-05-25 | 2023-05-25 | 减小高度散布的垂直发射远程导弹转弯程序角生成方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310600762.4A CN116576737A (zh) | 2023-05-25 | 2023-05-25 | 减小高度散布的垂直发射远程导弹转弯程序角生成方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116576737A true CN116576737A (zh) | 2023-08-11 |
Family
ID=87535603
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310600762.4A Pending CN116576737A (zh) | 2023-05-25 | 2023-05-25 | 减小高度散布的垂直发射远程导弹转弯程序角生成方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116576737A (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107121015A (zh) * | 2017-06-16 | 2017-09-01 | 湖北航天技术研究院总体设计所 | 一种快速弹上弹道在线规划方法 |
| CN111306998A (zh) * | 2020-03-24 | 2020-06-19 | 北京理工大学 | 一种参数摄动自适应的制导火箭弹垂直攻击制导方法 |
| KR102488431B1 (ko) * | 2022-10-26 | 2023-01-13 | 한화시스템(주) | 사격 제원 초기 파라미터 추정장치 및 그 방법 |
| CN115932778A (zh) * | 2023-02-23 | 2023-04-07 | 南京大学 | 一种低数据率下利用hrrp序列进行弹道中段目标识别的方法 |
| CN116045744A (zh) * | 2023-01-09 | 2023-05-02 | 航天科工火箭技术有限公司 | 一种固体运载火箭分离体残骸落区的控制方法和装置 |
-
2023
- 2023-05-25 CN CN202310600762.4A patent/CN116576737A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107121015A (zh) * | 2017-06-16 | 2017-09-01 | 湖北航天技术研究院总体设计所 | 一种快速弹上弹道在线规划方法 |
| CN111306998A (zh) * | 2020-03-24 | 2020-06-19 | 北京理工大学 | 一种参数摄动自适应的制导火箭弹垂直攻击制导方法 |
| KR102488431B1 (ko) * | 2022-10-26 | 2023-01-13 | 한화시스템(주) | 사격 제원 초기 파라미터 추정장치 및 그 방법 |
| CN116045744A (zh) * | 2023-01-09 | 2023-05-02 | 航天科工火箭技术有限公司 | 一种固体运载火箭分离体残骸落区的控制方法和装置 |
| CN115932778A (zh) * | 2023-02-23 | 2023-04-07 | 南京大学 | 一种低数据率下利用hrrp序列进行弹道中段目标识别的方法 |
Non-Patent Citations (1)
| Title |
|---|
| 聂川义;郭建国;周军;呼卫军;: "空射弹道导弹助推段弹道设计与优化", 飞行力学, no. 04, 15 August 2012 (2012-08-15), pages 357 - 361 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108168381B (zh) | 一种多枚导弹协同作战的控制方法 | |
| CN108362171A (zh) | 一种具有攻击时间和攻击角度约束的制导策略 | |
| CN111442697A (zh) | 一种基于伪谱法修正的过重补制导方法和弹道整形制导方法 | |
| CN111594322B (zh) | 一种基于Q-Learning的变循环航空发动机推力控制方法 | |
| CN114200828B (zh) | 一种超音速大机动靶标持续大过载防失速方法 | |
| CN115079565B (zh) | 变系数的带落角约束制导方法、装置和飞行器 | |
| CN114035616A (zh) | 一种飞行器对移动目标打击控制方法及系统 | |
| CN104656659B (zh) | 舰载机滑跃起飞自动飞行控制方法 | |
| CN105159311A (zh) | 一种用于捷联导引头的自动驾驶仪的设计方法 | |
| CN116576737A (zh) | 减小高度散布的垂直发射远程导弹转弯程序角生成方法 | |
| CN116592708A (zh) | 适用于多种时空约束耦合下的飞行器弹道成型制导方法 | |
| CN115657458A (zh) | 一种基于能量匹配的飞行器爬升轨迹控制方法 | |
| CN111680426B (zh) | 一种变系数比例导引参数设计方法 | |
| CN116301028B (zh) | 基于吸气式高超声速平台的多约束在线飞行轨迹规划中段导引方法 | |
| CN107145761B (zh) | 带动力滑翔导弹多样式弹道规划方法 | |
| CN107341295B (zh) | 具有终端角度和速度约束的下压段弹道设计方法 | |
| CN111737815B (zh) | 一种基于无量纲形式的飞行系统落点误差修正方法 | |
| CN113901645A (zh) | 一种应用于近程超声速巡航弹的定高飞行技术 | |
| CN108388272B (zh) | 一种基于时间最优的无人机前飞速度计算方法 | |
| CN117537670B (zh) | 一种提高远程导弹舵资源利用效率的舵分配方法 | |
| CN115342683B (zh) | 一种具有自动避障功能的调炮方法 | |
| CN117932793B (zh) | 采用耦合迭代优化的制导火箭弹内外弹道联合设计方法 | |
| KR102136267B1 (ko) | 지대지 유도탄의 궤적성형을 위한 가표적 선정 방법 및 이를 수행하는 장치 | |
| CN116540772A (zh) | 一种攻击角度和收敛时间控制的滑模制导律设计方法 | |
| CN116047906A (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 |