CN105718660A - 临近空间大范围机动弹道三维包络计算方法 - Google Patents
临近空间大范围机动弹道三维包络计算方法 Download PDFInfo
- Publication number
- CN105718660A CN105718660A CN201610039940.0A CN201610039940A CN105718660A CN 105718660 A CN105718660 A CN 105718660A CN 201610039940 A CN201610039940 A CN 201610039940A CN 105718660 A CN105718660 A CN 105718660A
- Authority
- CN
- China
- Prior art keywords
- cos
- phi
- sin
- lambda
- trajectory
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 38
- 230000008859 change Effects 0.000 claims abstract description 9
- 238000000034 method Methods 0.000 claims description 17
- 241001442234 Cosa Species 0.000 claims description 9
- 244000089409 Erythrina poeppigiana Species 0.000 claims description 9
- 235000009776 Rathbunia alamosensis Nutrition 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 9
- 238000013178 mathematical model Methods 0.000 claims description 9
- 230000014509 gene expression Effects 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 20
- 238000004088 simulation Methods 0.000 description 6
- 230000009466 transformation Effects 0.000 description 3
- 239000003814 drug Substances 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 206010048669 Terminal state Diseases 0.000 description 1
- 239000008186 active pharmaceutical agent Substances 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000002904 solvent Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种临近空间大范围机动弹道三维包络计算方法,包括以下步骤:建立弹道包络侧向宽度与纵程的数学模型;建立换极坐标系;建立换极坐标系中飞行器动力学模型;计算换极坐标系弹道参数;计算换极坐标系弹道侧向包络;计算一般坐标系弹道三维包络。本发明所述临近空间大范围机动弹道三维包络计算方法立足于飞行器机动能力分析,兼顾射向、纵程及建模误差等对弹道包络的影响,可广泛应用于具有相同动力学描述的飞行器,计算得到的包络可覆盖具有不同机动特性的弹道,能适应考虑禁飞区弹道及侧向最大范围机动弹道等极端情况,为飞行器机动能力分析和射前空间信息获取提供了重要的方法基础。
Description
技术领域
本发明涉及飞行器弹道规划领域,特别涉及一种临近空间大范围机动弹道三维包络计算方法。
背景技术
本发明涉及的临近空间的大范围机动弹道特指高超声速滑翔飞行器包含初始下降段弹道和滑翔段弹道在内的弹道,其起点为再入点,终点为滑翔段弹道终点。该段弹道是高超声速滑翔飞行器全程弹道的主要组成部分,具有长时间处于临近空间并进行大范围侧向机动的特性,弹道特征显著区别于弹道式弹道。在再入点和目标点确定的条件下,临近空间大范围机动弹道的侧向包络可达几千公里,弹道形态因飞行任务不同而迥异。
进行弹道包络计算的意义包括但不限于如下两个方面。首先,计算弹道包络是分析飞行器机动能力的重要手段,为弹道规划和飞行任务制定提供先验知识。其次,计算弹道覆盖区域为射前高精度引力模型构建提供必要的空间位置信息,是修正引力高阶项影响的必要条件,对于快速发射和提升飞行器导航、制导精度具有重要意义。目前,尚无公开资料提出有关临近空间大范围机动弹道三维包络的计算方法或准则。
发明内容
本发明针对临近空间大范围机动弹道覆盖区域确定问题,首次提出一种临近空间大范围机动弹道三维包络计算方法。方法基于飞行器能力分析结果,能够根据再入点和目标点快速确定弹道三维包络。计算得到的包络可覆盖具有不同机动特性的弹道,为飞行器机动能力分析和射前空间信息获取提供了重要的方法基础。
本发明通过以下技术方案来实现上述目的:
一种临近空间大范围机动弹道三维包络计算方法,所述临近空间大范围机动弹道包括初始下降段弹道和滑翔段弹道,弹道起点为再入点,弹道终点为滑翔段弹道终点,所述计算方法包括以下步骤:
(1)建立弹道包络侧向宽度与纵程的数学模型;
(2)建立换极坐标系;
(3)建立换极坐标系中飞行器动力学模型;
(4)计算换极坐标系弹道参数;
(5)计算换极坐标系弹道侧向包络;
(6)计算一般坐标系弹道三维包络。
具体地,
(1)建立弹道包络侧向宽度与纵程的数学模型
令弹道再入点为I,其经度为λI、地心纬度为φI;弹道落点为T,其经度为λT、地心纬度为φT;以由再入点和落点确定的再入大圆弧平面为对称面,称沿射向方向对称面以左弹道为左侧机动弹道,沿射向方向对称面以右弹道为右侧机动弹道;记最大左侧机动弹道距对称面的最大距离为左边界Bl,记最大右侧机动弹道距对称面的最大距离为右边界Br;
根据给定的飞行器再入点飞行状态参数、弹道终端飞行状态参数、飞行过程约束条件及终端约束条件,针对再入点为(λI,φI)=(0°,0°),落点为(λki,0°)(i=1,2,3,4,…)的低空大范围机动弹道进行弹道计算,记各落点对应的弹道纵程Li,由式(1)计算Li,
Li=Rearccos(sinφIsinφI+cosφIcosφIcos(λT-λT))(1)
式中:Re为地球半径;
通过改变再入点初始速度方位角获得最大左侧机动弹道和最大右侧向机动弹道,记各落点对应的弹道侧向包络宽度Wio,由式(2)计算Wio,
Wio=Bl-Br(2)
考虑建模误差,取n倍Wio为实际侧向包络宽度Wi,对Li(i=1,2,3,4,…)和Wi(i=1,2,3,4,…)进行拟合,建立如式(3)所示的弹道包络侧向宽度与纵程的数学模型,确定模型系数ai,
(2)建立换极坐标系
A.定义一个再入大圆弧平面作为换极赤道平面:
a)对目标点确定的情况,将再入点和目标点地心矢径构成的再入大圆弧平面作为换极赤道平面;
b)对于目标点未确定的情况,根据再入点位置及速度方位角确定的再入大圆弧平面作为换极赤道面;
B.基于换极赤道平面定义换极坐标系OE为地心,轴沿再入点地心矢径方向,轴在换极赤道面内垂直于轴指向目标点方向,轴与轴、轴构成右手系。
(3)建立换极坐标系中飞行器动力学模型
根据步骤(2)建立的换极坐标系建立换极坐标系中飞行器的动力模型,且用表示换极坐标系中各物理量,用x表示一般坐标系中各物理量:
在换极坐标系中建立以时间为自变量的滑翔飞行器动力学方程,其飞行状态量为换极后的经度地心纬度航迹偏航角速度速度倾角和地心距
其中,Cσ、Cθ为哥氏加速度项,和为牵连加速度项,
其中,
其中,ωe为地球旋转加速度矢量,λp和φp为换极后极点P的经度和地心纬度,Ap为P的方位角;
根据换极坐标系定义,一般坐标系与换极坐标系中地心距、当地速度倾角及速度的定义一致,
定义
其中,ψf为点F的方位角。
由一般坐标系中λ和φ确定换极坐标系中和的表达式为,
由换极坐标系中和确定一般坐标系中λ和φ的表达式为,
由一般坐标系中σ确定换极坐标系中的表达式为
其中,
(4)计算换极坐标系弹道参数
根据弹道再入点经度λI、再入点地心纬度φI、落点经度λT、落点地心纬度φT计算得到换极弹道参数,由计算得记换极弹道落点经度为
(5)计算换极坐标系弹道侧向包络
由式(15)计算换极弹道纵程
将代入式(16),可计算得到换极弹道侧向包络宽度
将换极弹道侧向包络描述为长度为宽度为的矩形,其边界由式(17)确定,
其中,为换极系侧向包络东向下边界,为换极系侧向包络东向上边界;为换极系侧向包络北向下边界,为换极系侧向包络北向上边界。
(6)计算一般坐标系弹道三维包络
按照东向dλ、北向dφ的间隔将由式(17)确定的换极坐标系弹道侧向包络剖分为互不重叠的若干网格,记网格节点坐标为
由步骤(3)所述坐标转换关系,计算一般坐标系网格节点坐标N(λn,φn),由N(λn,φn)确定的区域即为一般坐标系弹道侧向包络区域。
优选地,考虑机动和偏差进行蒙特卡罗打靶,以弹道高度范围的m倍作为弹道包络的垂向范围ΔH。
本发明临近空间大范围机动弹道三维包络计算方法计算得到的弹道包络可覆盖具有各种机动特性的弹道,对考虑禁飞区的中程机动弹道和侧向最大范围机动弹道等极端情况均有较好的适应性,且具有以下优点:
1)首次提出了一种针对临近空间大范围机动弹道的三维包络计算方法。方法能够根据再入点和目标点快速确定弹道三维包络,计算得到的包络可覆盖具有不同机动特性的弹道,为飞行器机动能力分析和射前空间信息获取提供了重要的方法基础。
2)基于极点变换的思想,首次将一种换极坐标系引入弹道三维包络计算方法的建立中。通过将弹道再入大圆弧平面转移到换极赤道面附近,使得换极坐标系中的弹道包络关于换极赤道面对称,从而简化了确定包络边界的计算流程。
3)引入分片定位的思想,对换极坐标系弹道包络进行空域剖分,确定包络特征点位置。根据换极系特征点位置,通过坐标转换求解一般系特征点位置,进而确定一般系弹道包络。根据计算时间和计算量要求,可灵活调节用于特征点确定的空域剖分参数。
4)方法计算流程简洁、计算速度快,确定模型参数后即可完成弹道三维包络的快速计算。模型参数确定后,无需通过多次弹道规划和弹道计算来确定飞行器可达域,因此能够满足射前任务快速规划和空间信息快速获取的要求。
5)方法可广泛适用于具有相同动力学描述的飞行器,不受飞行器特征参数改变及飞行任务改变的影响。方法对考虑禁飞区的中程机动弹道以及最大侧向机动弹道等极端情况均有较好的适应性,且可根据模型偏差调整包络宽度系数,具有较强的鲁棒性。
附图说明
图1为实施例中不同射向弹道侧向包络示意图;
图2为实施例中东向纵程5560km弹道侧向包络示意图;
图3为实施例中东向纵程6671km弹道侧向包络示意图;
图4为实施例中东向纵程7784km弹道侧向包络示意图;
图5为实施例中东向纵程8896km弹道侧向包络示意图;
图6为实施例中东向不同纵程弹道侧向包络示意图;
图7为实施例中换极坐标系及飞行状态量定义示意图;
图8为实施例中弹道再入点I和换极后极点P的位置关系示意图;
图9为实施例中换极坐标系航迹偏航角的位置关系示意图;
图10为实施例中换极系弹道侧向包络节点示意图,其中n=1.5;
图11为实施例中一般系弹道侧向包络节点示意图,其中n=1.5;
图12为实施例中实施例中考虑禁飞区弹道与侧向包络位置关系示意图,其中n=1.5;
图13为实施例中最大左侧机动弹道与侧向包络位置关系示意图,其中n=1.5;
图14为实施例中最大右侧机动弹道与侧向包络位置关系示意图,其中n=1.5;
图15为实施例中考虑禁飞区弹道与侧向包络位置关系示意图,其中n=1.1;
图16为实施例中最大左侧机动弹道与侧向包络位置关系示意图其中n=1.1;
图17为实施例中最大右侧机动弹道与侧向包络位置关系示意图其中n=1.1。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明:
基于CAV-H飞行器模型进行仿真,基本仿真条件设置为:①质量m=908kg,参考面积Sm=0.48375m2;②再入点状态参数:速度VI=6500m/s,速度倾角θI=0°,高度HI=80km;③滑翔段终端状态参数:速度Vk=2500m/s,速度倾角θk=0°,高度Hk=30km,终端航迹偏航角偏差不大于±5°;④飞行约束条件:最大热流密度最大动压qmax=100kPa,最大过载nmax=3g;⑤终端结束时距目标点的待飞航程:Stogo=100km。
仿真计算机配置为:Intel(R)Core(TM)i5-3470CPU3.20GHz,内存为3.46GB。软件环境为WindowXP操作系统,计算程序基于VC++6.0开发。其具体步骤如下:
第一步,建立弹道包络侧向宽度与纵程的数学模型
令弹道再入点为I,其经度为λI、地心纬度为φI;弹道落点为T,其经度为λT、地心纬度为φT;以由再入点和落点确定的再入大圆弧平面为对称面,称沿射向方向对称面以左弹道为左侧机动弹道,沿射向方向对称面以右弹道为右侧机动弹道;记最大左侧机动弹道距对称面的最大距离为左边界Bl,记最大右侧机动弹道距对称面的最大距离为右边界Br;
针对不同射向的临近空间大范围机动弹道进行仿真,其再入点经度为λI=0°、再入点地心纬度为φI=0°,落点经度λT和落点地心纬度φT见表1。
表1不同射向弹道落点坐标
射向 | 正北 | 正东 | 正南 | 正西 |
经度(°) | 0 | 60 | 0 | -60 |
地心纬度(°) | 60 | 0 | -60 | 0 |
计算得到不同射向弹道侧向包络示意图见图1,不同射向弹道侧向包络边界见表2。由图1和表2可见,由于地球自转和非球形因素的影响,东向弹道包络最大、南北向弹道包络次之,西向弹道包络最小;南北向弹道包络向东有小范围偏移。
表2不同射向弹道侧向包络边界
射向 | 正东 | 正南 | 正西 | 正北 |
左边界Bl(km) | 1597 | 1072 | 1072 | 1072 |
右边界Br(km) | 1597 | 1513 | 1072 | 1528 |
针对射向为正东、纵程不同的临近空间大范围机动弹道进行仿真。弹道再入点经度为λI=0°,再入点地心纬度为φI=0°;落点地心纬度为0°,落点经度分别为50°、60°、70°和80°。记各落点对应的弹道纵程为Li,由式(1)计算Li,
Li=Rearccos(sinφIsinφI+cosφIcosφIcos(λT-λT))(1)
其中,Re为地球半径。
计算得到不同落点对应的纵程Li见表3。
表3不同落点对应的纵程
落点经度(°) | 50 | 60 | 70 | 80 |
纵程Li(km) | 5560 | 6671 | 7784 | 8896 |
通过改变再入点初始速度方位角获得由飞行器能力决定的最大左侧机动弹道和最大右侧向机动弹道。不同纵程弹道侧向包络示意图见图2、图3、图4、图5和图6。
记各落点对应的弹道侧向包络宽度Wio,由式(2)计算Wio,
Wio=Bl-Br(2)
考虑建模误差,取n倍Wio为实际侧向包络宽度Wi,计算得到不同纵程弹道侧向包络宽度见表4。
表4不同纵程弹道侧向包络宽度
对Li(i=1,2,3,4,…)和Wi(i=1,2,3,4,…)进行拟合,建立如式(3)所示的弹道包络侧向宽度与纵程的数学模型,
具体为W=-1.064813L+0.000879L2-1.225047e-7L3+4.607571e-12L4
第二步,建立换极坐标系
A.定义一个再入大圆弧平面作为换极赤道平面:
对目标点确定的情况,将再入点和目标点地心矢径构成的再入大圆弧平面作为换极赤道平面;
对于目标点未确定的情况,根据再入点位置及速度方位角确定的再入大圆弧平面作为换极赤道面;
B.基于换极赤道平面定义换极坐标系OE为地心,轴沿再入点地心矢径方向,轴在换极赤道面内垂直于轴指向目标点方向,轴与轴、轴构成右手系。见图7。
第三步,建立换极坐标系中飞行器动力学模型
在换极坐标系中建立以时间为自变量的滑翔飞行器动力学方程,其飞行状态量为换极后的经度地心纬度航迹偏航角速度速度倾角和地心距
其中,Cσ、Cθ为哥氏加速度项,和为牵连加速度项,
其中,
其中,ωe为地球旋转加速度矢量,λp和φp为换极后极点P的经度和地心纬度,AP为P的方位角,见图8。
根据换极坐标系定义,一般坐标系与换极坐标系中地心距、当地速度倾角及速度的定义一致,
定义
其中,ψf为点F的方位角。
由一般坐标系中λ和φ确定换极坐标系中和的表达式为,
由换极坐标系中和确定一般坐标系中λ和φ的表达式为,
由一般坐标系中σ确定换极坐标系中的表达式为
其中,
第四步,计算换极坐标系弹道参数
根据弹道再入点经度λI=100°、再入点地心纬度φI=20°、落点经度λT=165°、落点地心纬度φT=50°计算得到换极弹道参数,由计算得
第五步,计算换极坐标系弹道侧向包络
由式(15)计算换极弹道纵程
计算带入式(16),可计算得到换极弹道侧向包络宽度
其中,模型系数ai由第一步确定。
由计算的将换极弹道侧向包络描述为长度为宽度为的矩形,其边界由式(17)确定,
其中,为换极系侧向包络东向下边界,为换极系侧向包络东向上边界;为换极系侧向包络北向下边界,为换极系侧向包络北向上边界。
由计算得,
第六步,计算一般坐标系弹道三维包络
按照东向dλ=1°、北向dφ=1°的间隔将由式(17)确定的换极坐标系弹道侧向包络剖分为互不重叠的若干网格,记网格节点坐标为换极系弹道侧向包络节点见图10。
由第三步所述坐标转换关系,计算一般坐标系网格节点坐标N(λn,φn),由N(λn,φn)确定的区域即为一般坐标系弹道侧向包络区域,一般系弹道侧向包络节点见图11。
考虑机动和偏差进行蒙特卡罗打靶,取弹道包络的垂向范围为沿参考弹道上下浮动10km,即ΔH=20km。经过以上计算,最终可得给定再入点和落点弹道的三维弹道包络。
针对再入点与落点与第四步所述相同的临近空间大范围机动弹道进行仿真。图12、图13和图14分别给出了考虑禁飞区弹道、最大左侧机动弹道和最大右侧机动弹道与计算所得侧向包络的位置关系示意图。其中,禁飞区中心经度λnf=138°、中心地心纬度φnf=40°、高度Hnf=100km、半径Rnf=1000km。由图12、图13和图14可见,极端情况下的弹道均位于计算所得的包络内。
取n=1.1,即以1.1倍Wio为实际侧向包络宽度Wi建立三维弹道包络数学模型。图15、图16和图17分别给出了考虑禁飞区弹道、最大左侧机动弹道和最大右侧机动弹道与计算所得包络的位置关系示意图,弹道仿真条件与图12、图13和图14中弹道相同。由图15、图16和图17可见,极端情况下的弹道均位于计算所得的包络内。
综合上述仿真结果可获得以下结论:
1)由本发明确定的方法进行临近空间大范围机动弹道三维包络计算,可使计算得到的包络覆盖具有不同机动特性的弹道,对考虑禁飞区的中程机动弹道以及最大侧向机动弹道等极端情况均有较好的适应性;
2)由本发明确定的方法进行弹道三维包络计算,模型参数一经固定,则无需再进行多次弹道计算来确定弹道可达域,避免了多次弹道积分带来的巨大运算量,因此具有计算速度快和轻量化的特征,能够满足射前任务快速规划和空间信息快速获取的要求;
3)建立模型所涉及的包络侧向宽度系数,确定包络特征点所涉及的空域剖分参数均可灵活调节,因此能够适应不同建模误差、不同计算时间要求和不同计算量要求下的弹道包络计算;
4)方法可广泛适用于具有相同动力学模型描述的飞行器,不受飞行器特征参数改变及飞行任务改变的影响,计算不出现奇点,具有适应范围广的特征。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种临近空间大范围机动弹道三维包络计算方法,所述临近空间大范围机动弹道包括初始下降段弹道和滑翔段弹道,弹道起点为再入点,弹道终点为滑翔段弹道终点,其特征在于:所述计算方法包括以下步骤:
(1)建立弹道包络侧向宽度与纵程的数学模型;
(2)建立换极坐标系;
(3)建立换极坐标系中飞行器动力学模型;
(4)计算换极坐标系弹道参数;
(5)计算换极坐标系弹道侧向包络;
(6)计算一般坐标系弹道三维包络。
2.根据权利要求1所述的临近空间大范围机动弹道三维包络计算方法,其特征在于:
(1)建立弹道包络侧向宽度与纵程的数学模型
令弹道再入点为I,其经度为λI、地心纬度为φI;弹道落点为T,其经度为λT、地心纬度为φT;以由再入点和落点确定的再入大圆弧平面为对称面,称沿射向方向对称面以左弹道为左侧机动弹道,沿射向方向对称面以右弹道为右侧机动弹道;记最大左侧机动弹道距对称面的最大距离为左边界Bl,记最大右侧机动弹道距对称面的最大距离为右边界Br;
根据给定的飞行器再入点飞行状态参数、弹道终端飞行状态参数、飞行过程约束条件及终端约束条件,针对再入点为(λI,φI)=(0°,0°),落点为(λki,0°)(i=1,2,3,4,…)的低空大范围机动弹道进行弹道计算,记各落点对应的弹道纵程Li,由式(1)计算Li,
Li=Rearccos(sinφIsinφI+cosφIcosφIcos(λT-λT))(1)
式中:Re为地球半径;
通过改变再入点初始速度方位角获得最大左侧机动弹道和最大右侧向机动弹道,记各落点对应的弹道侧向包络宽度Wio,由式(2)计算Wio,
Wio=Bl-Br(2)
取n倍Wio为实际侧向包络宽度Wi,对Li(i=1,2,3,4,…)和Wi(i=1,2,3,4,…)进行拟合,建立如式(3)所示的弹道包络侧向宽度与纵程的数学模型,确定模型系数ai,
(2)建立换极坐标系
A.定义一个再入大圆弧平面作为换极赤道平面:
a)对目标点确定的情况,将再入点和目标点地心矢径构成的再入大圆弧平面作为换极赤道平面;
b)对于目标点未确定的情况,根据再入点位置及速度方位角确定的再入大圆弧平面作为换极赤道面;
B.基于换极赤道平面定义换极坐标系OE为地心,轴沿再入点地心矢径方向,轴在换极赤道面内垂直于轴指向目标点方向,轴与轴、轴构成右手系。
(3)建立换极坐标系中飞行器动力学模型
根据步骤(2)建立的换极坐标系建立换极坐标系中飞行器的动力模型,且用表示换极坐标系中各物理量,用x表示一般坐标系中各物理量:
在换极坐标系中建立以时间为自变量的滑翔飞行器动力学方程,其飞行状态量为换极后的经度地心纬度航迹偏航角速度速度倾角和地心距
其中,Cσ、Cθ为哥氏加速度项,和为牵连加速度项,
其中,
其中,ωe为地球旋转加速度矢量,λp和φp为换极后极点P的经度和地心纬度,Ap为P的方位角;
根据换极坐标系定义,一般坐标系与换极坐标系中地心距、当地速度倾角及速度的定义一致,
定义
其中,ψf为点F的方位角。
由一般坐标系中λ和φ确定换极坐标系中和的表达式为,
由换极坐标系中和确定一般坐标系中λ和φ的表达式为,
由一般坐标系中σ确定换极坐标系中的表达式为
其中,
(4)计算换极坐标系弹道参数
根据步骤(3)的坐标转换关系计算换极坐标系弹道参数,根据弹道再入点经度λI、再入点地心纬度φI、落点经度λT、落点地心纬度φT计算得到换极弹道参数,由计算得记换极弹道落点经度为
(5)计算换极坐标系弹道侧向包络
由式(15)计算换极弹道纵程
其中,Re为地球半径。
将代入式(16),可计算得到换极弹道侧向包络宽度
其中,模型系数ai由第一步确定。
将换极弹道侧向包络描述为长度为宽度为的矩形,其边界由式(17)确定,
其中,为换极系侧向包络东向下边界,为换极系侧向包络东向上边界;为换极系侧向包络北向下边界,为换极系侧向包络北向上边界。
(6)计算一般坐标系弹道三维包络
按照东向dλ、北向dφ的间隔将由式(17)确定的换极坐标系弹道侧向包络剖分为互不重叠的若干网格,记网格节点坐标为
由步骤(3)所述坐标转换关系,计算一般坐标系网格节点坐标N(λn,φn),由N(λn,φn)确定的区域即为一般坐标系弹道侧向包络区域。
3.根据权利要求2所述的临近空间大范围机动弹道三维包络计算方法,其特征在于:考虑机动和偏差进行蒙特卡罗打靶,以弹道高度范围的m倍作为弹道包络的垂向范围ΔH。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610039940.0A CN105718660B (zh) | 2016-01-21 | 2016-01-21 | 临近空间大范围机动弹道三维包络计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610039940.0A CN105718660B (zh) | 2016-01-21 | 2016-01-21 | 临近空间大范围机动弹道三维包络计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105718660A true CN105718660A (zh) | 2016-06-29 |
CN105718660B CN105718660B (zh) | 2019-03-01 |
Family
ID=56153644
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610039940.0A Active CN105718660B (zh) | 2016-01-21 | 2016-01-21 | 临近空间大范围机动弹道三维包络计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105718660B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106547991A (zh) * | 2016-11-25 | 2017-03-29 | 中国工程物理研究院总体工程研究所 | 沿滑翔弹道的扰动引力重构模型优化方法 |
CN107832528A (zh) * | 2017-11-13 | 2018-03-23 | 中国运载火箭技术研究院 | 一种空天飞行器三维空间可达域分析方法 |
CN109507877A (zh) * | 2018-11-22 | 2019-03-22 | 北京航空航天大学 | 飞行器弹道规划安全性评估方法 |
CN109800506A (zh) * | 2019-01-22 | 2019-05-24 | 南京航空航天大学 | 一种飞行器的性能评估方法及系统 |
CN113377128A (zh) * | 2021-06-10 | 2021-09-10 | 北京空天技术研究所 | 一种飞行器可达区域估算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103838914A (zh) * | 2013-12-30 | 2014-06-04 | 北京航空航天大学 | 一种高超声速飞行器滑翔段弹道解析求解方法 |
CN104035335A (zh) * | 2014-05-27 | 2014-09-10 | 北京航空航天大学 | 基于高精度纵、横程解析预测方法的平稳滑翔再入制导律 |
CN104392047A (zh) * | 2014-11-25 | 2015-03-04 | 北京航空航天大学 | 一种基于平稳滑翔弹道解析解的快速弹道规划方法 |
CN104634182A (zh) * | 2014-12-16 | 2015-05-20 | 北京控制工程研究所 | 一种跳跃式再入标准弹道在线修正的跟踪制导方法 |
CN105115512A (zh) * | 2015-09-23 | 2015-12-02 | 北京理工大学 | 一种火星大气进入段侧向预测校正制导方法 |
-
2016
- 2016-01-21 CN CN201610039940.0A patent/CN105718660B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103838914A (zh) * | 2013-12-30 | 2014-06-04 | 北京航空航天大学 | 一种高超声速飞行器滑翔段弹道解析求解方法 |
CN104035335A (zh) * | 2014-05-27 | 2014-09-10 | 北京航空航天大学 | 基于高精度纵、横程解析预测方法的平稳滑翔再入制导律 |
CN104392047A (zh) * | 2014-11-25 | 2015-03-04 | 北京航空航天大学 | 一种基于平稳滑翔弹道解析解的快速弹道规划方法 |
CN104634182A (zh) * | 2014-12-16 | 2015-05-20 | 北京控制工程研究所 | 一种跳跃式再入标准弹道在线修正的跟踪制导方法 |
CN105115512A (zh) * | 2015-09-23 | 2015-12-02 | 北京理工大学 | 一种火星大气进入段侧向预测校正制导方法 |
Non-Patent Citations (4)
Title |
---|
HUAN,XUEYING AN, ETC.: "Fast local representation of gravity anomaly along flight trajectory", 《PROCEEDINGS OF THE INSTITUTION OF MECHANICAL ENGINEERS PART G JOURNAL OF AEROSPACE ENGINEERING》 * |
YU XIE, LUHUA LIU,ETC.: "Rapid generation of entry trajectories with waypoint and no-fly zone constraints", 《ACTA ASTRONAUTICA》 * |
郑书娥: "滚转发射机弹分离轨迹的仿真研究", 《现代防御技术》 * |
钱山,郑伟,张士峰,蔡洪: "一种弹道导弹再入弹道解析方法", 《飞行力学》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106547991A (zh) * | 2016-11-25 | 2017-03-29 | 中国工程物理研究院总体工程研究所 | 沿滑翔弹道的扰动引力重构模型优化方法 |
CN106547991B (zh) * | 2016-11-25 | 2020-11-03 | 中国工程物理研究院总体工程研究所 | 沿滑翔弹道的扰动引力重构模型优化方法 |
CN107832528A (zh) * | 2017-11-13 | 2018-03-23 | 中国运载火箭技术研究院 | 一种空天飞行器三维空间可达域分析方法 |
CN109507877A (zh) * | 2018-11-22 | 2019-03-22 | 北京航空航天大学 | 飞行器弹道规划安全性评估方法 |
CN109800506A (zh) * | 2019-01-22 | 2019-05-24 | 南京航空航天大学 | 一种飞行器的性能评估方法及系统 |
CN113377128A (zh) * | 2021-06-10 | 2021-09-10 | 北京空天技术研究所 | 一种飞行器可达区域估算方法 |
CN113377128B (zh) * | 2021-06-10 | 2022-12-09 | 北京空天技术研究所 | 一种飞行器可达区域估算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105718660B (zh) | 2019-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111306989B (zh) | 一种基于平稳滑翔弹道解析解的高超声速再入制导方法 | |
CN105718660B (zh) | 临近空间大范围机动弹道三维包络计算方法 | |
CN107966156B (zh) | 一种适用于运载火箭垂直回收段的制导律设计方法 | |
CN104035335B (zh) | 基于高精度纵、横程解析预测方法的平稳滑翔再入制导方法 | |
CN109740198B (zh) | 一种基于解析预测的滑翔飞行器三维再入制导方法 | |
CN106547991B (zh) | 沿滑翔弹道的扰动引力重构模型优化方法 | |
WO2019109407A1 (zh) | 一种基于飞行路径角规划的再入轨迹设计方法 | |
EP2817689B1 (en) | Safe emergency landing of a uav | |
CN104392047B (zh) | 一种基于平稳滑翔弹道解析解的快速弹道规划方法 | |
CN107679655A (zh) | 一种航天发射火箭落点预测系统 | |
CN104615144A (zh) | 基于目标规划的高超声速飞行器再入轨迹在线优化方法 | |
CN106054604B (zh) | 基于模型预测控制理论的再入飞行器鲁棒最优制导方法 | |
CN108549785B (zh) | 一种基于三维飞行剖面的高超声速飞行器精准弹道快速预测方法 | |
Liang et al. | Lateral entry guidance with no-fly zone constraint | |
CN105740506B (zh) | 沿临近空间大范围机动弹道空间包络的扰动引力逼近方法 | |
CN113900448B (zh) | 一种基于滑模干扰观测器的飞行器预测校正复合制导方法 | |
Coombes et al. | Landing site reachability in a forced landing of unmanned aircraft in wind | |
CN105138808A (zh) | 基于摄动理论的滑翔弹道误差传播分析方法 | |
CN108766035B (zh) | 一种点密度引导下的无人机地形匹配飞行控制系统 | |
Dalle et al. | Flight envelope calculation of a hypersonic vehicle using a first principles-derived model | |
CN102508492B (zh) | 一种飞行器在等高航路点间的定高度大圆飞行实现方法 | |
CN105760573B (zh) | 沿临近空间大范围机动弹道的扰动引力延拓逼近方法 | |
CN109977543B (zh) | 基于侧向优先的三维剖面覆盖区域计算方法、系统及介质 | |
CN105354380B (zh) | 面向摄动因素影响补偿的滑翔弹道快速修正方法 | |
CN113093789B (zh) | 一种基于路径点优选的飞行器禁飞区规避轨迹规划方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |