[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN116361926B - 一种弹道导弹滑翔增程段初始机械能闭环调整方法 - Google Patents

一种弹道导弹滑翔增程段初始机械能闭环调整方法 Download PDF

Info

Publication number
CN116361926B
CN116361926B CN202310638532.7A CN202310638532A CN116361926B CN 116361926 B CN116361926 B CN 116361926B CN 202310638532 A CN202310638532 A CN 202310638532A CN 116361926 B CN116361926 B CN 116361926B
Authority
CN
China
Prior art keywords
mechanical energy
range
missile
steps
normalized
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
Application number
CN202310638532.7A
Other languages
English (en)
Other versions
CN116361926A (zh
Inventor
苗昊春
高登巍
栗金平
邓海鹏
李琪
常江
潘瑞
范中平
裴培
张一杰
张哲明
李雅君
张梦典
张建松
李昊远
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Institute of Modern Control Technology
Original Assignee
Xian Institute of Modern Control Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xian Institute of Modern Control Technology filed Critical Xian Institute of Modern Control Technology
Priority to CN202310638532.7A priority Critical patent/CN116361926B/zh
Publication of CN116361926A publication Critical patent/CN116361926A/zh
Application granted granted Critical
Publication of CN116361926B publication Critical patent/CN116361926B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling 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)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)

Abstract

本发明公开了一种弹道导弹滑翔增程段初始机械能闭环调整方法,首先建立剩余射程和机械能的函数关系,从而达到预测任意剩余射程下机械能的目的;其次利用对爬升方案角的闭环修正,使得导弹低空爬升时间增加或者缩短来达到附加阻力的目的进行闭环能量修正。该策略计算效率高,矫正正能力强,工程易实现。

Description

一种弹道导弹滑翔增程段初始机械能闭环调整方法
技术领域
本发明属于导弹技术领域,具体涉及一种弹道导弹滑翔增程段初始机械能闭环调整方法。
背景技术
现有的弹道导弹能量管理技术通常在滑翔増程段附加侧滑角(STT)或附加倾侧角(BTT)进行能量修正。但由于弹道导弹的能量散布主要来自于固体发动机的推力散布,且弹道导弹的初始能量全来自于发动机推力,因此推力对初始能量的影响非常大,较大的初始推力散布限制了附加侧滑角和附加倾侧角减速的能力。此外,在滑翔増程段通过附加大侧滑角和倾侧角会导致导弹飞行稳定性的下降,限制了减速能力有限。故如何修正滑翔增程段初始能量使得其初始能量处于后续能量管理能力范围内有重要意义。
发明内容
为了克服现有技术的不足,本发明提供了一种弹道导弹滑翔增程段初始机械能闭环调整方法,首先建立剩余射程和机械能的函数关系,从而达到预测任意剩余射程下机械能的目的;其次利用对爬升方案角的闭环修正,使得导弹低空爬升时间增加或者缩短来达到附加阻力的目的进行闭环能量修正。该策略计算效率高,矫正能力强,工程易实现。
本发明解决其技术问题所采用的技术方案包括如下步骤:
步骤1:计算归一化机械能和剩余射程;
归一化机械能:
(1)
其中,/>分别为导弹归一化速度和归一化地心距,/>为导弹速度;/>分别为速度和距离缩放比例;/>,/>分别为水平面地心引力加速度和地球平均半径;/>为发射系下的地心矢量;
归一化剩余射程:
(2)
其中为剩余弹目距离,/>分别为发射系下的导弹位置和目标位置三分量;
步骤2:建立归一化剩余射程和归一化机械能函数关系;
当不考虑地球旋转和推力对弹道的影响时得到如下关于时间的导数关系:
(3)
其中为当地弹道倾角,/>为气动阻力,/>为气动升力;
将式(3)带入机械能和剩余射程表达式(1)和(2)得:
(4)
时主动段机械能时剩余射程的单调递减函数;如果考虑推力则:
(5)
其中为推力在速度方向上的投影;
步骤3:爬升方案角闭环修正;
最高点处满足,此时机械能最大限度的以重力势能的形式存在并包含部分动能,机械能和高度的函数关系为:
(6)
设标称情况下主动段结束时到最高点处的机械能为,当到达最大高度的机械能为/>,标称的爬升方案角为/>,主动段结束后到最高点处的爬升方案角为:
(7)
优选地,所述,/> m。
本发明的有益效果如下:
本发明方法在线计算效率高且闭环修正效果好,可以满足导弹在线能量闭环调整的要求,同时具有很强的通用性。设计的滑翔增程段初始能量散布小,保证了后续精确能量管理的适应性,具有广阔的应用前景。
附图说明
图1是本发明实施例滑翔增程段初始能量调整示意图。
图2是本发明实施例蒙特卡洛打靶弹道剩余射程和机械能曲线。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
为扩展已有弹道导弹能量管理能力实现主动段散布大的情况下能量管理的适应范围,需研究一种针对滑翔増程段初始能量闭环修正的方法。本发明提供了一种基于剩余射程机械能预测和爬升方案角闭环校正的算法进行弹道导弹初始能量校准,适用于固体火箭推力散布较大导致的导弹进入滑翔段时的初始能量散布大的情况,主要用于缩小弹道导弹滑翔増程段的初始能量散布。该方法包含给定剩余射程下的机械能预测以及剩余爬升方案角矫正系数计算两个部分组成。该算法在给定剩余射程处判定机械能,通过闭环修正主动段结束到最高点处的爬升方案角,通过变更爬升方案角改变气动阻力实现在最高点处的速度趋于一致。
步骤一,计算归一化机械能和剩余射程
归一化机械能:
其中为发射系下的地心矢量,可通过相对发射点的方位直接计算获得。
归一化剩余射程:
其中为剩余弹目距离,/>为发射系下的导弹位置和目标位置三分量。
步骤二,建立归一化剩余射程和归一化机械能函数关系
当不考虑地球旋转和推力对弹道的影响时可得到如下关于时间的导数关系:
其中为当地弹道倾角,/>为气动阻力,/>为气动升力。将上式带入机械能和剩余射程表达式可得:
从表达式可以看出,当时主动段机械能时剩余射程的单调递减函数;而如果考虑推力则:
其中为推力在速度方向上的投影,由于/>,显然主动段机械能相对剩余射程保持单调递增,故主动段结束时即为机械能的最值点。因此通过主动段结束时的机械能峰值可以判断在主动段推力散布影响下的机械能散布。
步骤三,爬升方案角闭环修正
最高点处满足,此时机械能最大限度的以重力势能的形式存在并包含部分动能。机械能和高度的函数关系为:
从上式可知通过调整爬升方案角修正机械能具有很大的调节效率。设标称情况下主动段结束时到最高点处的机械能为,当到达最大高度的机械能为/>,标称的爬升方案角为/>,设计主动段结束后到最高点处的爬升方案角为:
具体实施例:
结合某滑翔式弹道导弹算例对本发明做进一步描述。导弹发射点海拔0m,对射程x=1200km处,侧偏z=0m处h=0m海拔的导弹进行打击。
一,计算发射点地心矢径
首先转换角度单位:
其中卯酉圆半径为;偏心率为/>;/>—参考椭球的长、短半轴,分别为6378137m、6356752.3141m;/>为发射点和目标点所在位置的椭球高。
二,根据标称弹道计算弹道地心距
根据发射点地心矢径和标称弹道数据,计算弹道地心矢径/>
计算归一化地心矢径:
三,计算标称机械能
根据归一化缩放标称速度:,然后计算机械能:
四,计算闭环爬升方案角
根据归一化机械能和地心矢径建立的插值关系,结合标称爬升方案角/>,计算闭环爬升方案角:
为了验证本发明的可行性,分别进行了本发明提出的方案角对应机械能修正后的仿真结果,通过数学仿真对比得图1和图2曲线。由仿真结果可知:本发明设计的爬升方案角闭环修正方法可有效消除主动段结束后的能量误差,为滑翔增程段初始能量的归一化建立合理初值,配合增程段速度控制使得绝大部分蒙特卡洛打靶仿真弹道剩余能量最终趋近于一致。

Claims (2)

1.一种弹道导弹滑翔增程段初始机械能闭环调整方法,其特征在于,包括如下步骤:
步骤1:计算归一化机械能和剩余射程;
归一化机械能:
(1)
其中,/>分别为导弹归一化速度和归一化地心距,为导弹速度;/>分别为速度和距离缩放比例;/>,/>分别为水平面地心引力加速度和地球平均半径;/>为发射系下的地心矢量;
归一化剩余射程:
(2)
其中为剩余弹目距离,/>分别为发射系下的导弹位置和目标位置三分量;
步骤2:建立归一化剩余射程和归一化机械能函数关系;
当不考虑地球旋转和推力对弹道的影响时得到如下关于时间的导数关系:
(3)
其中为当地弹道倾角,/>为气动阻力,/>为气动升力;
将式(3)带入机械能和剩余射程表达式(1)和(2)得:
(4)
时主动段机械能时剩余射程的单调递减函数;如果考虑推力则:
(5)
其中为推力在速度方向上的投影;
步骤3:爬升方案角闭环修正;
最高点处满足,此时机械能最大限度的以重力势能的形式存在并包含部分动能,机械能和高度的函数关系为:
(6)
设标称情况下主动段结束时到最高点处的机械能为,当到达最大高度的机械能为,标称的爬升方案角为/>,主动段结束后到最高点处的爬升方案角为:
(7)。
2.根据权利要求1所述的一种弹道导弹滑翔增程段初始机械能闭环调整方法,其特征在于,所述,/> m。
CN202310638532.7A 2023-06-01 2023-06-01 一种弹道导弹滑翔增程段初始机械能闭环调整方法 Active CN116361926B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310638532.7A CN116361926B (zh) 2023-06-01 2023-06-01 一种弹道导弹滑翔增程段初始机械能闭环调整方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310638532.7A CN116361926B (zh) 2023-06-01 2023-06-01 一种弹道导弹滑翔增程段初始机械能闭环调整方法

Publications (2)

Publication Number Publication Date
CN116361926A CN116361926A (zh) 2023-06-30
CN116361926B true CN116361926B (zh) 2023-08-18

Family

ID=86923828

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310638532.7A Active CN116361926B (zh) 2023-06-01 2023-06-01 一种弹道导弹滑翔增程段初始机械能闭环调整方法

Country Status (1)

Country Link
CN (1) CN116361926B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104035335A (zh) * 2014-05-27 2014-09-10 北京航空航天大学 基于高精度纵、横程解析预测方法的平稳滑翔再入制导律
WO2020222250A1 (en) * 2019-04-30 2020-11-05 Sivaprakasam P Modified re-entry vehicle design with dynamic trajectory glide control system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104035335A (zh) * 2014-05-27 2014-09-10 北京航空航天大学 基于高精度纵、横程解析预测方法的平稳滑翔再入制导律
WO2020222250A1 (en) * 2019-04-30 2020-11-05 Sivaprakasam P Modified re-entry vehicle design with dynamic trajectory glide control system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
弹道导弹在被动段轴向加速变轨弹道设计与优化;张斌;任军学;王一白;刘宇;;北京航空航天大学学报(第06期);全文 *

Also Published As

Publication number Publication date
CN116361926A (zh) 2023-06-30

Similar Documents

Publication Publication Date Title
CN111306998B (zh) 一种参数摄动自适应的制导火箭弹垂直攻击制导方法
CN108168381B (zh) 一种多枚导弹协同作战的控制方法
CN108362171B (zh) 一种具有攻击时间和攻击角度约束的制导策略
CN111336871B (zh) 一种基于迂回式飞行的垂直攻击制导方法
CN111692919B (zh) 超近射程的飞行器精确制导控制方法
CN109737812B (zh) 空对地制导武器侧向攻击方法和装置
CN111580547B (zh) 一种高超声速飞行器编队控制方法
CN111473696A (zh) 一种基于落点估计的制导火箭垂直攻击制导方法
CN110425943B (zh) 面向变质心飞行器的工程化再入制导方法
CN104567545B (zh) Rlv大气层内主动段的制导方法
US20090212108A1 (en) Method of Determining a Fire Guidance Solution
CN114035616A (zh) 一种飞行器对移动目标打击控制方法及系统
CN116361926B (zh) 一种弹道导弹滑翔增程段初始机械能闭环调整方法
CN115657458A (zh) 一种基于能量匹配的飞行器爬升轨迹控制方法
CN115079565A (zh) 变系数的带落角约束制导方法、装置和飞行器
CN114153143A (zh) 一种导弹非奇异固定时间滑模制导律的设计方法
CN111680426B (zh) 一种变系数比例导引参数设计方法
CN113834385B (zh) 一种针对初始航迹角自由的控制时间的二维协同制导方法
CN117419609A (zh) 一种用于有翼火箭的电磁发射方法
CN116576735A (zh) 一种超远程制导火箭气动热主动缓解控制方法
CN113608547B (zh) 一种滑翔制导炸弹基于风补偿的机载火控方法及系统
CN115993772A (zh) 一种基于贝塞尔曲线的四阶段二维制导方法
CN107341295B (zh) 具有终端角度和速度约束的下压段弹道设计方法
Nobahari et al. Integrated optimization of guidance and control parameters in a dual spin flying vehicle
CN115342683B (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