CN114036727B - 一种基于末制导的滑翔弹道快速规划方法 - Google Patents
一种基于末制导的滑翔弹道快速规划方法 Download PDFInfo
- Publication number
- CN114036727B CN114036727B CN202111235427.6A CN202111235427A CN114036727B CN 114036727 B CN114036727 B CN 114036727B CN 202111235427 A CN202111235427 A CN 202111235427A CN 114036727 B CN114036727 B CN 114036727B
- Authority
- CN
- China
- Prior art keywords
- glide
- attack angle
- trajectory
- angle alpha
- segment
- 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
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000013461 design Methods 0.000 claims abstract description 15
- 230000009189 diving Effects 0.000 claims abstract description 13
- 238000012407 engineering method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
本申请涉及一种基于末制导的滑翔弹道快速规划方法,属于导弹弹道设计技术领域,包括:根据中、末交班时刻弹目距离LJB划分滑翔段和俯冲下压段;根据目标点经纬度和滑翔高度要求,确定滑翔段预测目标点;滑翔段根据拟平衡滑翔攻角αHX和带落角约束的最优比例导引率求取的滑翔段攻角αJBQ,按照滑翔段剩余航程进行加权后形成新的滑翔攻角αHX';俯冲下压段采用带落角约束的最优比例导引律求取俯冲下压段攻角αJBH;判断终端落速vD是否满足期望终端落速vM,若是,则规划结束,若否,则对滑翔攻角αHX'进行迭代计算,求出满足期望终端落速vM的滑翔段设计攻角αsj,即可生成可行弹道。本申请基于末制导的滑翔弹道快速规划方法简单新颖,弹道生成速度快,非常适合工程应用,且适合在线规划。
Description
技术领域
本申请涉及导弹弹道设计技术领域,特别涉及一种基于末制导的滑翔弹道快速规划方法。
背景技术
滑翔导弹为了能够对特定目标实施精确打击,一般对终端打击精度、速度和角度都有一定的约束条件。国内外文献资料对于满足终端约束的滑翔弹道规划方法主要涉及很少,大多致力于末制导方法研究。国内有关满足终端约束的滑翔弹道规划工程方法大多是采用多变量迭代规划的设计方法,该方法不足之处是迭代变量多、迭代速度慢。因此,亟待找到一种满足终端约束的滑翔弹道快速规划方法,适合在线规划,能够较好地应用到工程设计中。
发明内容
本申请实施例提供一种基于末制导的滑翔弹道快速规划方法,以解决相关技术中滑翔弹道规划工程方法迭代变量多、迭代速度慢的问题。
本申请实施例提供了一种基于末制导的滑翔弹道快速规划方法,所述方法包括:
S1、根据中、末交班时刻弹目距离LJB划分滑翔段和俯冲下压段;
S2、根据目标点经纬度和滑翔高度要求,确定滑翔段预测目标点;
S3、滑翔段采用拟平衡滑翔攻角αHX与带落角约束的最优比例导引率求取的滑翔段攻角αJBQ按照滑翔段剩余航程进行加权后形成新的滑翔攻角αHX';
S4、俯冲下压段采用带落角约束的最优比例导引律求取俯冲下压段攻角αJBH;
S5、判断终端落速vD是否满足期望终端落速vM,若是,则进入S6,若否,则迭代计算出满足期望终端落速vM的滑翔段设计攻角αsj,并对滑翔攻角αHX'进行迭代更新,返回S3;
S6、规划结束,输出满足终端约束的滑翔弹道。
在一些实施例中:所述滑翔攻角αHX'的计算公式如下:
其中:
S:滑翔段剩余航程;
S0:滑翔初始点剩余航程。
在一些实施例中:所述拟平衡滑翔攻角αHX为滑翔段设计攻角αsj与滑翔段当地弹道倾角θd的差值,所述滑翔攻角αHX'的计算公式如下:
其中:
S:滑翔段剩余航程;
S0:滑翔初始点剩余航程。
在一些实施例中:所述滑翔段设计攻角αsj的计算公式如下:
αsj n+1=αsj n-(vD-vM)/H;
其中:
αsj n+1:第n+1次迭代的滑翔段设计攻角;
αsj n:第n次迭代的滑翔段设计攻角;
H:设计迭代步长。
在一些实施例中:步骤S3和步骤S4在标准条件下使用三自由度弹道模型计算弹道,标准条件包括:
a)地球模型为旋转椭球体;
b)大气条件为国家标准大气;
c)所有高度上无风;
d)弹道计算的初值为滑翔初始点的速度、位置参数;
e)符合瞬时平衡原理。
本申请提供的技术方案带来的有益效果包括:
本申请实施例提供了一种基于末制导的滑翔弹道快速规划方法,通过滑翔段与俯冲下压段的划分,以及最优比例导引律、攻角加权设计方法的应用,只需对滑翔设计攻角进行迭代计算,即可生成一条满足终端精度、角度、速度约束的滑翔弹道,将理论与工程设计进行了有效结合,设计新颖、简单,可靠性强,迭代速度快,适合在线规划,具有很好的工程应用价值。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例的方法流程图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请的一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本申请保护的范围。
本申请实施例提供了一种基于末制导的滑翔弹道快速规划方法,其能解决相关技术中滑翔弹道规划工程方法迭代变量多、迭代速度慢的问题。
参见图1所示,本申请实施例提供了一种基于末制导的滑翔弹道快速规划方法,所述方法包括:
S1、滑翔段和俯冲下压段划分:
根据中、末交班时刻弹目距离LJB划分滑翔段和俯冲下压段,中、末交班时刻弹目距离LJB=60km。
S2、滑翔段预测目标点的确定:
根据目标点经度λM、纬度和最低滑翔要求高程,最低滑翔要求高程hMIN=30km,确定滑翔段预测目标点经纬高为λHX、/>hHX;
其中:
λHX=λM
hHX=hMIN
S3、滑翔段弹道设计:
当滑翔段弹目距离>LJB时(例如LJB=60km),滑翔飞行器处于滑翔段,滑翔段垂直平面采用拟平衡滑翔攻角αHX与带落角约束的最优比例导引率求取的滑翔段攻角αJBQ按照滑翔段剩余航程进行加权,加权后形成新的滑翔攻角αHX'。
拟平衡滑翔攻角αHX为滑翔段设计攻角αsj与滑翔段当地弹道倾角θd的差值,滑翔攻角αHX'的计算公式如下:
其中:
S:滑翔段剩余航程;
S0:滑翔初始点剩余航程。
其中,αJBQ是根据带落角约束的最优比例导引律求取而来,求取过程这里不再赘述,θd为滑翔段当地弹道倾角,滑翔段预测目标点的期望落角设计为0°。
S4、俯冲下压段弹道设计:
当俯冲下压段弹目距离≤LJB时(例如LJB=60km),滑翔飞行器处于俯冲下压段,俯冲下压段采用带落角约束的最优比例导引律求取俯冲下压段攻角αJBH;其中,俯冲下压段的期望落角设计为80°。
S5、落速迭代:
根据上述步骤S3、步骤S4中的攻角,求解导弹运动方程,得到终端落速vD,判断终端落速vD是否满足期望终端落速vM=700m/s±10m/s,若是,进入下一步S6,若否,对滑翔攻角αHX'进行迭代,迭代出满足终端落速vM的滑翔段设计攻角αsj,并返回S3;
滑翔段设计攻角αsj的计算公式如下:
αsj n+1=αsj n-(vD-vM)/H;
其中:
αsj n+1:第n+1次迭代的滑翔段设计攻角;
αsj n:第n次迭代的滑翔段设计攻角;
H:设计迭代步长,H=200。
S6、规划结束,输出满足终端约束的滑翔弹道。
其中,步骤S3和步骤S4在标准条件下使用三自由度弹道模型计算弹道,标准条件包括:
a)地球模型为旋转椭球体;
b)大气条件为国家标准大气;
c)所有高度上无风;
d)弹道计算的初值为滑翔初始点的速度、位置参数;
e)符合瞬时平衡原理。
在本申请的描述中,需要说明的是,术语“上”、“下”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本申请的限制。除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本申请中的具体含义。
需要说明的是,在本申请中,诸如“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅是本申请的具体实施方式,使本领域技术人员能够理解或实现本申请。对这些实施例的多种修改对本领域的技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本申请的精神或范围的情况下,在其它实施例中实现。因此,本申请将不会被限制于本文所示的这些实施例,而是要符合与本文所申请的原理和新颖特点相一致的最宽的范围。
Claims (4)
1.一种基于末制导的滑翔弹道快速规划方法,其特征在于,所述方法包括:
S1、根据中、末交班时刻弹目距离LJB划分滑翔段和俯冲下压段;
S2、根据目标点经纬度和滑翔高度要求,确定滑翔段预测目标点;
S3、滑翔段采用拟平衡滑翔攻角αHX与带落角约束的最优比例导引率求取的滑翔段攻角αJBQ按照滑翔段剩余航程进行加权后形成新的滑翔攻角αHX';
S4、俯冲下压段采用带落角约束的最优比例导引律求取俯冲下压段攻角αJBH;
S5、判断终端落速vD是否满足期望终端落速vM,若是,则进入S6,若否,则迭代计算出满足期望终端落速vM的滑翔段设计攻角αsj,并对滑翔攻角αHX'进行迭代更新,返回S3;
所述滑翔段设计攻角αsj的计算公式如下:
其中:
第n+1次迭代的滑翔段设计攻角;
第n次迭代的滑翔段设计攻角;
H:设计迭代步长;
S6、规划结束,输出满足终端约束的滑翔弹道。
2.如权利要求1所述的一种基于末制导的滑翔弹道快速规划方法,其特征在于,所述滑翔攻角αHX'的计算公式如下:
其中:
S:滑翔段剩余航程;
S0:滑翔初始点剩余航程。
3.如权利要求1或2所述的一种基于末制导的滑翔弹道快速规划方法,其特征在于:
所述拟平衡滑翔攻角αHX为滑翔段设计攻角αsj与滑翔段当地弹道倾角θd的差值,所述滑翔攻角αHX'的计算公式如下:
其中:
S:滑翔段剩余航程;
S0:滑翔初始点剩余航程。
4.如权利要求1所述的一种基于末制导的滑翔弹道快速规划方法,其特征在于:
步骤S3和步骤S4在标准条件下使用三自由度弹道模型计算弹道,标准条件包括:
a)地球模型为旋转椭球体;
b)大气条件为国家标准大气;
c)所有高度上无风;
d)弹道计算的初值为滑翔初始点的速度、位置参数;
e)符合瞬时平衡原理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111235427.6A CN114036727B (zh) | 2021-10-22 | 2021-10-22 | 一种基于末制导的滑翔弹道快速规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111235427.6A CN114036727B (zh) | 2021-10-22 | 2021-10-22 | 一种基于末制导的滑翔弹道快速规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114036727A CN114036727A (zh) | 2022-02-11 |
CN114036727B true CN114036727B (zh) | 2024-06-18 |
Family
ID=80141750
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111235427.6A Active CN114036727B (zh) | 2021-10-22 | 2021-10-22 | 一种基于末制导的滑翔弹道快速规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114036727B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115629618B (zh) * | 2022-11-04 | 2024-10-01 | 航天科工火箭技术有限公司 | 一种基于落点选择和伪谱法的分离体最优弹道规划的方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105930550A (zh) * | 2016-04-01 | 2016-09-07 | 方洋旺 | 一种吸气式高超声速导弹助推-跳跃式弹道优化方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107145761B (zh) * | 2017-06-18 | 2019-12-03 | 湖北航天技术研究院总体设计所 | 带动力滑翔导弹多样式弹道规划方法 |
CN108398959A (zh) * | 2018-03-28 | 2018-08-14 | 北京临近空间飞行器系统工程研究所 | 一种高超声速滑翔飞行器快速下压制导控制方法 |
CN110471450B (zh) * | 2019-08-29 | 2020-07-14 | 大连理工大学 | 在高度速度剖面内直接规划再入轨迹的方法 |
-
2021
- 2021-10-22 CN CN202111235427.6A patent/CN114036727B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105930550A (zh) * | 2016-04-01 | 2016-09-07 | 方洋旺 | 一种吸气式高超声速导弹助推-跳跃式弹道优化方法 |
Non-Patent Citations (1)
Title |
---|
滑翔飞行器弹道在线规划及制导技术研究;罗志伟等;战术导弹技术;20150731(第4期);第66-71页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114036727A (zh) | 2022-02-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112613250B (zh) | 火星进入器表面流动转捩位置预测方法 | |
Richardson | Gulf Stream trajectories measured with free-drifting buoys | |
CN110595472B (zh) | 一种无人船双目标气象航线优化方法和系统 | |
Wang et al. | The increase of snowfall in Northeast China after the mid-1980s | |
CN114036727B (zh) | 一种基于末制导的滑翔弹道快速规划方法 | |
CN104298851B (zh) | 一种用于预报强降水天气的数据处理方法 | |
CN108894768B (zh) | 一种基于蝙蝠算法和井壁稳定的钻进轨迹设计方法与系统 | |
CN111678524B (zh) | 一种基于飞行安全的救援航空器路径规划方法及系统 | |
CN110750850B (zh) | 强约束复杂任务条件下的三维剖面优化设计方法、系统及介质 | |
CN112711899A (zh) | 一种蒸发波导高度的融合预测方法 | |
Wang | An inverse balance equation in sigma coordinates for model initialization | |
CN115860268B (zh) | 基于cnn-xgboost模型的短时强降水预测方法 | |
US20150355215A1 (en) | Approach to assess available wind resource distribution based on interpolation method | |
CN110147585A (zh) | 一种基于Kriging改进模型的防风栅气动优化方法 | |
CN113326624A (zh) | 一种沙漠地区大气边界层高度的预测方法及系统 | |
CN115795792A (zh) | 一种研究南大洋中尺度涡旋个数季节变化规律的方法 | |
CN114509074B (zh) | 一种在洋流补偿下进行航位推算的水下滑翔机定位方法 | |
CN107274012B (zh) | 基于云进化粒子群算法的短期风电功率预测方法 | |
Mashiko et al. | Super high-resolution simulation of the 6 May 2012 Tsukuba supercell tornado: Near-Surface structure and its evolution | |
Zhou et al. | Constrained least-squares method for computing equivalent static wind loads of large-span roofs | |
CN113032510B (zh) | 一种基于高斯权重的lbgm法初值扰动生成方法 | |
CN110231619A (zh) | 基于恩克法的雷达交接时刻预报方法及装置 | |
CN105160108A (zh) | 一种探空火箭气动外形优化目标函数的确定方法 | |
Cheong et al. | Tropical cyclone track and intensity prediction with a structure adjustable balanced vortex | |
Cao et al. | Preliminary assessment of the wind power resource around the thousand-meter scale megatall building |
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 |