CN107121015A - 一种快速弹上弹道在线规划方法 - Google Patents

一种快速弹上弹道在线规划方法 Download PDF

Info

Publication number
CN107121015A
CN107121015A CN201710461448.7A CN201710461448A CN107121015A CN 107121015 A CN107121015 A CN 107121015A CN 201710461448 A CN201710461448 A CN 201710461448A CN 107121015 A CN107121015 A CN 107121015A
Authority
CN
China
Prior art keywords
angle
attack
trajectory
gliding
number table
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
Application number
CN201710461448.7A
Other languages
English (en)
Other versions
CN107121015B (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.)
General Designing Institute of Hubei Space Technology Academy
Original Assignee
General Designing Institute of Hubei Space Technology Academy
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 General Designing Institute of Hubei Space Technology Academy filed Critical General Designing Institute of Hubei Space Technology Academy
Priority to CN201710461448.7A priority Critical patent/CN107121015B/zh
Publication of CN107121015A publication Critical patent/CN107121015A/zh
Application granted granted Critical
Publication of CN107121015B publication Critical patent/CN107121015B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F41WEAPONS
    • F41GWEAPON SIGHTS; AIMING
    • F41G3/00Aiming or laying means
    • F41G3/22Aiming or laying means for vehicle-borne armament, e.g. on aircraft

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)

Abstract

本发明公开了一种快速弹上弹道在线规划方法,对于再入式机动弹头,机动段采用射面内的长距离滑翔式弹道,包括如下步骤:确定机动段弹道的攻角、侧滑角变化模式;确定机动段弹道的对机动距离、落速均有重大影响的调节参数;调节参数插值数表制备;使用上一步骤制备的调节参数插值数表在线规划机动弹道。本发明的方法使弹上计算机性能较低的导弹也能具备弹上在线规划的能力。预先制备调节参数插值数表,通过对调节参数插值数表进行插值得到弹道在线规划调节参数,再进行弹道在线规划,计算量小,无需迭代计算,仅通过单次计算就能生成满足落速、落角要求以及落点精度要求的机动弹道。该方法适用于无侧向机动、能够长距离滑行的再入式机动弹头。

Description

一种快速弹上弹道在线规划方法
技术领域
本发明属于航天飞行器总体设计技术领域,具体涉及一种快速弹上弹道在线规划方法。
背景技术
随着导弹武器系统对射程提出了更高的要求,以及耗尽关机闭路制导技术的逐渐成熟,发动机耗尽关机方式开始在越来越多的地地弹道导弹上采用。与精确关机方式相比,发动机耗尽关机需要面对的一个重要问题是发动机工作过程中由于装药生产工艺、环境温度等因素产生的总冲、燃速偏差。这些偏差体现在弹道上,主要表现为发动机关机后的速度偏差和再入点位置偏差。近年来得益于耗尽关机闭路制导技术的发展,耗尽关机产生的偏差在很大程度上得到控制。尤其是带有变轨发动机的导弹,可以在变轨段通过交变姿态的方式耗散多余的能量,但当变轨段的待增视速度方向与最佳速度修正方向(对射程影响最大的方向)比较接近时,仍会产生较大的再入点偏差,可达10km量级,对于制导与控制非常不利,即便制导系统强行追随标准弹道,亦会造成落速、落角与设计值相差较大。
为解决再入点位置偏差过大的问题,一种有效的解决方式是以新的再入点为起点重新生成一条机动弹道,即在线规划方式。现有的许多在线规划技术需要进行弹上迭代计算,这对弹上计算机的性能要求较高,对于那些弹上计算机性能较低的导弹无法满足条件。
发明内容
为解决再入点位置偏差过大而弹上计算机性能较低的导弹以新的再入点为起点重新生成一条机动弹道的技术问题,本发明提供一种快速弹上弹道在线规划方法,对于再入式机动弹头,机动段采用射面内的长距离滑翔式弹道,其特征在于,所述弹道在线规划方法包括如下步骤:
步骤1、确定机动段弹道的攻角、侧滑角变化模式,由于机动段在固定的射面内飞行、不作横向机动,侧滑角取值为0,机动段飞行分为四段,攻角按下面的方式变化:
第一段,从再入时刻起,经一定时间,攻角按固定的线性速率增加到最大正攻角;
第二段,从到最大正攻角时刻起到滑翔结束时刻为滑翔段,两者时间差值为滑翔时间,滑翔段的攻角数值为滑翔攻角减去当前的当地弹道倾角,其中滑翔攻角设计为常值,当计算得到的攻角数值大于最大正攻角时则限幅为最大正攻角,记录滑翔结束时刻的攻角数值,以进行下一阶段的计算;
第三段,从滑翔结束时刻到最小负攻角时刻,攻角以固定的线性速率从滑翔结束时刻的攻角数值减小到最小负攻角;
第四段,从最小负攻角时刻到落地为最后一段,该阶段弹头将以负攻角下压弹道,攻角数值为用于限定落角的设计角度值减去当前的当地弹道倾角,其中用于限定落角的设计角度值为常值,当计算得到的攻角数值小于最小负攻角时限幅为最小负攻角,当计算得到的攻角数值大于0时则限幅为0;
步骤2、确定机动段弹道的对机动距离、落速均有重大影响的调节参数用于弹上在线规划生成机动弹道;
步骤3、调节参数插值数表制备,先制备步骤2确定的调节参数的调节参数插值数表;
步骤4、使用步骤3制备的调节参数插值数表在线规划机动弹道。
具体地,当弹道的机动距离与再入速度可建立一一对应关系时,所述步骤2中所述的机动段弹道的调节参数为滑翔攻角和滑翔时间。
进一步地,所述步骤3中所述的调节参数插值数表制备具体为滑翔攻角插值数表制备和滑翔时间插值数表制备,其方法包括如下步骤:
步骤a1、首先确定再入点射面内纵向位置差;
步骤b1、再确定再入速度、目标点高程;
步骤c1、最后确定的调节参数插值数表为两个三维数表,滑翔攻角插值数表,表内三维数据分别为再入速度、目标点高程和纵向位置差;滑翔时间插值数表,表内三维数据分别为再入速度、目标点高程和纵向位置差;
步骤d1、计算每组再入速度、目标点高程和纵向位置差组合下的滑翔攻角和滑翔时间,使用滑翔攻角和滑翔时间数据计算出的每条弹道的落速、落角满足设计值,即生成了一组用于弹道在线规划的调节参数插值数表。
具体地,当弹道的机动距离与再入速度不再相关时,所述步骤2中所述的机动段弹道的调节参数为滑翔攻角和机动射程随滑翔时间的变化率。
进一步地,所述步骤3中所述的调节参数插值数表制备具体为滑翔攻角插值数表制备和机动射程随滑翔时间的变化率插值数表制备,其方法包括如下步骤:
步骤a2、确定再入速度、目标点高程、滑翔时间;
步骤b2、再确定机动射程随滑翔时间的变化率;
步骤c2、最后确定的调节参数插值数表为两个三维数表,滑翔攻角插值数表,表内三维数据分别为再入速度、目标点高程和滑翔时间;机动射程随滑翔时间的变化率插值数表表内三维数据分别为再入速度、目标点高程和滑翔时间;
步骤d2、计算每组再入速度、目标点高程和滑翔时间组合下的滑翔攻角和机动射程随滑翔时间的变化率,使用滑翔攻角和机动射程随滑翔时间的变化率数据计算出的每条弹道的落速、落角满足设计值,即生成了一组用于弹道在线规划的调节参数插值数表。
优选地,所述步骤1确定滑翔时间,当再入点射面内纵向位置出现偏差时,可由机动射程变化量除以机动射程随滑翔时间的变化率首先求得滑翔时间增量,再将原先的滑翔时间用滑翔时间加上滑翔时间增量代替,利用新的滑翔时间通过滑翔攻角数表插值确定滑翔攻角。
本发明的一种快速弹上弹道在线规划方法,使弹上计算机性能较低的导弹也能具备弹上在线规划的能力。其主要特点是预先制备调节参数插值数表,通过对调节参数插值数表进行插值得到弹道在线规划调节参数,再进行弹道在线规划,计算量小,无需迭代计算,仅通过单次计算就能生成满足落速、落角要求以及落点精度要求的机动弹道。该方法适用于无侧向机动、能够长距离滑行的再入式机动弹头。
附图说明
图1为机动段弹道攻角分四段变化示意图;
图2为不同偏差条件下的在线规划弹道与标准弹道的弹道形状对比;
图3为不同偏差条件下的在线规划弹道与标准弹道的攻角对比;
图4为不同偏差条件下的在线规划弹道与标准弹道的速度曲线对比。
具体实施方式
下面结合具体实施例和附图对本发明作进一步说明,但不应该以此限制本发明的保护范围。
实施例一:
本发明以的一种快速弹上弹道在线规划方法,对于再入式机动弹头,机动段采用射面内的长距离滑翔式弹道,机动弹道主要由攻角、侧滑角的设计决定,所述弹道在线规划方法包括如下步骤:
步骤1、确定机动段弹道的攻角ALP、侧滑角BET变化模式,由于机动段在固定的射面内飞行、不作横向机动,侧滑角BET取值为0,机动段飞行分为四段,攻角ALP按下面的分段情况变化,如图1所示:
第一段,从再入时刻T1起,经一定时间到T2时刻,攻角ALP按固定的线性速率(例如6°/s)增加到最大正攻角ZALP;
第二段,从到最大正攻角时刻T2起到滑翔结束时刻T3为滑翔段,两者时间差值(T3-T2)为滑翔时间TT2,滑翔段的攻角ALP数值为滑翔攻角ALP0减去当前的当地弹道倾角THED(即ALP=ALP0-THED),其中滑翔攻角ALP0设计为常值,当计算得到的攻角ALP数值大于最大正攻角ZALP时则限幅为最大正攻角ZALP数值,记录滑翔结束时刻攻角ALP32数值,以进行下一阶段的计算;
第三段,从滑翔结束时刻T3到最小负攻角时刻T4,ALP攻角以固定的线性速率(例如-7°/s)从滑翔结束时刻攻角ALP32数值减小到最小负攻角FALP;
第四段,从最小负攻角时刻T4到落地为最后一段,该阶段弹头将以负攻角下压弹道,攻角ALP数值为用于限定落角的设计角度ALP1值减去当前的当地弹道倾角THED(即ALP=ALP1-THED),例如设计角度是-80°则令ALP1=-80°,其中用于限定落角的设计角度值为常值,当计算得到的攻角ALP数值小于最小负攻角FALP时限幅为最小负攻角FALP数值,当计算得到的攻角ALP数值大于0时则限幅为0;
步骤2、确定机动段弹道的对机动距离、落速均有重大影响的调节参数,当弹道的机动距离与再入速度可建立一一对应关系时,确定机动段弹道的对机动距离、落速均有重大影响的调节参数,由步骤1可知,最大正攻角ZALP、最小负攻角FALP、滑翔攻角ALP0、滑翔时间TT2都是可变参量,由于用于限定落角的设计角度ALP1值是固定值可以预先确定;为了简化设计,应当以尽量减少可调参量为原则,所以攻角ALP的上下限ZALP、FALP也可以设计为固定值。剩下还有滑翔攻角ALP0和滑翔时间TT2两个参量,计算表明这两个参量对机动距离、落速都有影响,故这两个量将共同决定机动弹道的形式,用于弹上弹道在线规划生成机动弹道时首先要计算的也是这两个参量;
步骤3、调节参数插值数表制备,由于该在线规划方法不进行迭代计算,故需要一次性确定机动弹道的两个参数滑翔攻角ALP0和滑翔时间TT2,预先制备步骤2确定的滑翔攻角插值数表制备和滑翔时间插值数表制备,由于要应对再入点位置偏差,滑翔攻角ALP0、滑翔时间TT2首先应当和再入点位置偏差相关,考虑到中段闭路制导后再入点的位置偏差主要在纵向偏差,所以纵向位置差ΔL是确定滑翔攻角ALP0、滑翔时间TT2的一个参量,另外机动段弹道还与再入速度zrV、目标点高程mbdH相关。最后确定的调节参数插值数表为两个3维数表:ALP0(zrV,mbdH,ΔL)、TT2(zrV,mbdH,ΔL),其具体制备方法顺次包括如下步骤:
步骤a1、首先确定再入点射面内纵向位置差ΔL;
步骤b1、再确定再入速度zrV、目标点高程mbdH;
步骤c1、最后确定的调节参数插值数表为两个三维数表,滑翔攻角插值数表,表内三维数据分别为再入速度、目标点高程和纵向位置差,即ALP0(zrV,mbdH,ΔL);滑翔时间插值数表,表内三维数据分别为再入速度、目标点高程和纵向位置差,即TT2(zrV,mbdH,ΔL);
步骤d1、计算每组再入速度zrV、目标点高程mbdH和纵向位置差ΔL组合下的滑翔攻角ALP0和滑翔时间TT2,使用滑翔攻角ALP0和滑翔时间TT2数据计算出的每条弹道的落速、落角满足设计值,即生成了一组用于弹道在线规划的滑翔攻角插值数表和滑翔时间插值数表。
步骤4、使用步骤3制备的调节参数插值数表即滑翔攻角插值数表和滑翔时间插值数表在线规划机动弹道。具体在线规划机动弹道实现方法顺次包括如下步骤:
第1步:读入射表装订诸元中的标准机动弹道参数:规划算法标志位GHI(整数)、最大正攻角ZALP(单位为°)、最小负攻角FALP(单位为°)、滑翔攻角ALP0(单位为°)、滑翔时间TT2(单位为s),以及标准再入点在发射系内的坐标(X0,Y0,Z0)、发射系内的再入速度(VX0,VY0,VZ0)、发射系再入弹道倾角THE0、再入时刻T0。GHI=0代表不进行在线规划(射程较小时,闭路制导能力足够修正偏差,且在线规划规律性不强故不引入在线规划);GHI=1表示在线规划采用算法1;GHI=2代表在线规划采用算法2。出现的坐标均为发射系下的坐标,下同。
第2步:变轨发动机工作结束后,根据当时的速度位置算出实际的发射系内再入点坐标(X45,Y45,Z45)。弹上规划只针对纵向位置偏差,所以Z45实际并没有用处。
第3步:计算再入点纵向偏差ΔL,计算方法如下:
定义ΔX=X45-X0,ΔY=Y45-Y0,ΔZ=ΔX*VZ0/VX0,dL=sqrt(ΔX^2+ΔY^2+ΔZ^2),sqrt为开平方根。
则ΔL=0.993*dL*SIGN(ΔX),(常数0.993是为了将45km高度处的位置偏差转换到0km高程处)。
注意dL为大于等于0的数,ΔL则有正负号之分,SIGN(ΔX)表示取ΔX的符号;
第4步:确定在线规划计算的初始条件(T,VX,VY,VZ,X,Y,Z),即再入时刻、发射系内的速度分量、发射系内的位置分量。各项按下式确定:
定义dTHE=ΔL/6371000.79,(6371000.79为地球平均半径),dTHE为再入点位置产生偏差后,弹道倾角随之产生的变化量。
T=T0
VX=VX0-dTHE*sin(THE0)*sqrt(VX0^2+VY0^2)
VY=VY0+dTHE*cos(THE0)*sqrt(VX0^2+VY0^2)
VZ=VZ0
X=X0+ΔX=X45
Y=Y0+ΔY=Y45
Z=Z0+ΔZ;
第5步:3维数表ALP0(mbdH,Vzr,ΔL)和TT2(mbdH,Vzr,ΔL)线性插值得到一组新的ALP0和TT2。例如可以固定两个数表形状为6*19*11,则三个插值参数的数组如下(为和实际的高程、再入速度、再入点偏差加以区别插值数组的名称前加了“s”;插值参数此处为一个具体例子,这些参数的个数、具体数值可依据实际情况设计):
smbdH/0,500,1000,1500,2000,2500/(m)
sVzr/1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400/(m/s)
sΔL/-10000,-8000,-6000,-4000,-2000,0,2000,4000,6000,8000,10000/(m)
第6步:利用第4步与第5步结果计算机动弹道。
实施例二:
本发明以的一种快速弹上弹道在线规划方法,对于再入式机动弹头,机动段采用射面内的长距离滑翔式弹道,机动弹道主要由攻角、侧滑角的设计决定,所述弹道在线规划方法包括如下步骤:
步骤1、同实施例一;
步骤2、确定机动段弹道的对机动距离、落速均有重大影响的调节参数,当弹道的机动距离与再入速度不再相关时,在大射程条件下,不同射程可以对应同一个再入速度zrV,实施例一中的三维数组(zrV,mbdH,ΔL)不能完全确定机动弹道参数,如果将ΔL换成机动段距离L,可以部分解决该问题,但这时的机动段参数还和射向相关,数表会扩大为四维,为使计算简单,考虑该情况下弹道射程主要靠机动段的滑翔时间TT2确定,所以可以考虑由数组(zrV,mbdH,TT2)确定ALP0和机动射程L随滑翔时间TT2的变化率D(D=ΔL/ΔTT2),即确定两个数表ALP0(zrV,mbdH,TT2)和D(zrV,mbdH,TT2)。当再入点射面内纵向位置出现偏差时,可由ΔL/D首先求得TT2增量ΔTT2,再将原先的TT2用TT2+ΔTT2代替,利用新的TT2通过数表插值确定ALP0和D。使用这种高阶求偏导数的方式可将射向的影响大幅缩小,从而将数表维持在3维。本实施例中再入点射面内纵向位置为出现偏差的情况;因此,调节参数确定为滑翔攻角ALP0和机动射程L随滑翔时间TT2的变化率D;
步骤3、调节参数插值数表制备,两个3维数表:滑翔攻角插值数表ALP0(zrV,mbdH,TT2)和机动射程随滑翔时间的变化率插值数表D(zrV,mbdH,TT2),其具体制备方法顺次包括如下步骤:
步骤a2、确定再入速度zrV、目标点高程mbdH、滑翔时间TT2;
步骤b2、再确定机动射程随滑翔时间的变化率D;
步骤c2、最后确定的调节参数插值数表为两个三维数表,滑翔攻角插值数表,表内三维数据分别为再入速度、目标点高程和滑翔时间,即ALP0(zrV,mbdH,TT2);机动射程随滑翔时间的变化率插值数表,表内三维数据分别为再入速度、目标点高程和滑翔时间,即D(zrV,mbdH,TT2);
步骤d2、计算每组再入速度zrV、目标点高程mbdH和滑翔时间TT2组合下的滑翔攻角ALP0和机动射程随滑翔时间的变化率D,使用滑翔攻角ALP0和机动射程随滑翔时间的变化率D数据计算出的每条弹道的落速、落角满足设计值,即生成了一组用于弹道在线规划的调节参数插值数表。
步骤4、使用步骤3制备的调节参数插值数表即滑翔攻角插值数表ALP0(zrV,mbdH,TT2)和机动射程随滑翔时间的变化率插值数表D(zrV,mbdH,TT2)在线规划机动弹道。具体在线规划机动弹道实现方法顺次包括如下步骤:
第1步:读入射表装订诸元中的标准机动弹道参数:规划算法标志位GHI(整数)、最大正攻角ZALP(单位为°)、最小负攻角FALP(单位为°)、滑翔攻角ALP0(单位为°)、滑翔时间TT2(单位为s),以及标准再入点在发射系内的坐标(X0,Y0,Z0)、发射系内的再入速度(VX0,VY0,VZ0)、发射系再入弹道倾角THE0、再入时刻T0。GHI=0代表不进行在线规划(射程较小时,闭路制导能力足够修正偏差,且在线规划规律性不强故不引入在线规划);GHI=1表示在线规划采用算法1;GHI=2代表在线规划采用算法2。出现的坐标均为发射系下的坐标,下同。
第2步:变轨发动机工作结束后,根据当时的速度位置算出实际的发射系内再入点坐标(X45,Y45,Z45)。弹上规划只针对纵向位置偏差,所以Z45实际并没有用处。
第3步:计算再入点纵向偏差ΔL,计算方法如下:
定义ΔX=X45-X0,ΔY=Y45-Y0,ΔZ=ΔX*VZ0/VX0,dL=sqrt(ΔX^2+ΔY^2+ΔZ^2),sqrt为开平方根。
则ΔL=0.993*dL*SIGN(ΔX),(常数0.993是为了将45km高度处的位置偏差转换到0km高程处)。
注意dL为大于等于0的数,ΔL则有正负号之分,SIGN(ΔX)表示取ΔX的符号。
第4步:确定在线规划计算的初始条件(T,VX,VY,VZ,X,Y,Z)。即再入时刻、发射系内的速度分量、发射系内的位置分量。各项按下式确定:
定义dTHE=ΔL/6371000.79,(6371000.79为地球平均半径),dTHE为再入点位置产生偏差后,弹道倾角随之产生的变化量。
T=T0
VX=VX0-dTHE*sin(THE0)*sqrt(VX0^2+VY0^2)
VY=VY0+dTHE*cos(THE0)*sqrt(VX0^2+VY0^2)
VZ=VZ0
X=X0+ΔX=X45
Y=Y0+ΔY=Y45
Z=Z0+ΔZ;
第5步:3维数表D(mbdH,Vzr,TT2)插值得到机动段射程对TT2的偏导数D。得到TT2修正值dTT2=ΔL/D,则新的TT2=TT2+dTT2。
用新的TT2插值数表ALP0(mbdH,Vzr,TT2)得到新的ALP0。
以下作为例子将两个数表D(mbdH,Vzr,TT2)和ALP0(mbdH,Vzr,TT2)形状设计为6*19*15,其三个插值参数的数组如下(为和实际的高程、再入速度、滑翔时间加以区别,插值参数的数组名称前加了“s”;插值参数此处为一个具体例子,这些参数的个数、具体数值可依据实际情况设计):
smbdH/0,500,1000,1500,2000,2500/(m)
sVzr/1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400/(m/s)
sTT2/50,60,70,80,90,100,110,120,130,140,150,160,170,180,190/(s)
第6步:利用第4步与第5步结果计算机动弹道。
经测试,以上在线规划方法可以保证重新规划的机动弹道落点和实际目标点的位置偏差在大约2km的范围内。
图2~4为某一射程下的标准机动弹道以及再入点纵向偏差±10.km时的在线规划弹道参数的对比。图2为高度-航程曲线,横轴代表航程、纵轴代表高度,表现了弹道在射面内的几何形状,其中Hbz表示标准状态下的弹道、Hgh+表示再入点偏差10km时的规划弹道、Hgh-表示再入点偏差-10km时的规划弹道;图3为攻角-时间曲线,横轴代表时间、纵轴代表攻角,ALPbz表示标准弹道的攻角、ALPgh+表示再入点偏差10km时的规划弹道、ALPgh-表示再入点偏差-10km时的规划弹道;图4为速度-时间曲线,横轴代表时间、纵轴代表速度,Vbz表示标准弹道、Vhg+表示再入点偏差10km时的规划弹道、Vhg-表示再入点偏差-10km时的规划弹道。从图中可以看到偏差条件下的规划弹道在落点以及落速上均保持和标准弹道的一致性,说明了规划算法的有效性。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (6)

1.一种快速弹上弹道在线规划方法,对于再入式机动弹头,机动段采用射面内的长距离滑翔式弹道,其特征在于,所述弹道在线规划方法包括如下步骤:
步骤1、确定机动段弹道的攻角、侧滑角变化模式,由于机动段在固定的射面内飞行、不作横向机动,侧滑角取值为0,机动段飞行分为四段,攻角按下面的方式变化:
第一段,从再入时刻起,经一定时间,攻角按固定的线性速率增加到最大正攻角;
第二段,从到最大正攻角时刻起到滑翔结束时刻为滑翔段,两者时间差值为滑翔时间,滑翔段的攻角数值为滑翔攻角减去当前的当地弹道倾角,其中滑翔攻角设计为常值,当计算得到的攻角数值大于最大正攻角时则限幅为最大正攻角,记录滑翔结束时刻的攻角数值,以进行下一阶段的计算;
第三段,从滑翔结束时刻到最小负攻角时刻,攻角以固定的线性速率从滑翔结束时刻的攻角数值减小到最小负攻角;
第四段,从最小负攻角时刻到落地为最后一段,该阶段弹头将以负攻角下压弹道,攻角数值为用于限定落角的设计角度值减去当前的当地弹道倾角,其中用于限定落角的设计角度值为常值,当计算得到的攻角数值小于最小负攻角时限幅为最小负攻角,当计算得到的攻角数值大于0时则限幅为0;
步骤2、确定机动段弹道的对机动距离、落速均有重大影响的调节参数用于弹上在线规划生成机动弹道;
步骤3、调节参数插值数表制备,先制备步骤2确定的调节参数的调节参数插值数表;
步骤4、使用步骤3制备的调节参数插值数表在线规划机动弹道。
2.根据权利要求1所述的弹道在线规划方法,其特征在于,当弹道的机动距离与再入速度可建立一一对应关系时,所述步骤2中所述的机动段弹道的调节参数为滑翔攻角和滑翔时间。
3.根据权利要求2所述的弹道在线规划方法,其特征在于所述步骤3中所述的调节参数插值数表制备具体为滑翔攻角插值数表制备和滑翔时间插值数表制备,其方法包括如下步骤:
步骤a1、首先确定再入点射面内纵向位置差;
步骤b1、再确定再入速度、目标点高程;
步骤c1、最后确定的调节参数插值数表为两个三维数表,滑翔攻角插值数表,表内三维数据分别为再入速度、目标点高程和纵向位置差;滑翔时间插值数表,表内三维数据分别为再入速度、目标点高程和纵向位置差;
步骤d1、计算每组再入速度、目标点高程和纵向位置差组合下的滑翔攻角和滑翔时间,使用滑翔攻角和滑翔时间数据计算出的每条弹道的落速、落角满足设计值,即生成了一组用于弹道在线规划的调节参数插值数表。
4.根据权利要求1所述的弹道在线规划方法,其特征在于,当弹道的机动距离与再入速度不再相关时,所述步骤2中所述的机动段弹道的调节参数为滑翔攻角和机动射程随滑翔时间的变化率。
5.根据权利要求4所述的弹道在线规划方法,其特征在于所述步骤3中所述的调节参数插值数表制备具体为滑翔攻角插值数表制备和机动射程随滑翔时间的变化率插值数表制备,其方法包括如下步骤:
步骤a2、确定再入速度、目标点高程、滑翔时间;
步骤b2、再确定机动射程随滑翔时间的变化率;
步骤c2、最后确定的调节参数插值数表为两个三维数表,滑翔攻角插值数表,表内三维数据分别为再入速度、目标点高程和滑翔时间;机动射程随滑翔时间的变化率插值数表表内三维数据分别为再入速度、目标点高程和滑翔时间;
步骤d2、计算每组再入速度、目标点高程和滑翔时间组合下的滑翔攻角和机动射程随滑翔时间的变化率,使用滑翔攻角和机动射程随滑翔时间的变化率数据计算出的每条弹道的落速、落角满足设计值,即生成了一组用于弹道在线规划的调节参数插值数表。
6.根据权利要求5所述的弹道在线规划方法,其特征在于所述步骤1确定滑翔时间,当再入点射面内纵向位置出现偏差时,可由机动射程变化量除以机动射程随滑翔时间的变化率首先求得滑翔时间增量,再将原先的滑翔时间用滑翔时间加上滑翔时间增量代替,利用新的滑翔时间通过滑翔攻角数表插值确定滑翔攻角。
CN201710461448.7A 2017-06-16 2017-06-16 一种快速弹上弹道在线规划方法 Active CN107121015B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710461448.7A CN107121015B (zh) 2017-06-16 2017-06-16 一种快速弹上弹道在线规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710461448.7A CN107121015B (zh) 2017-06-16 2017-06-16 一种快速弹上弹道在线规划方法

Publications (2)

Publication Number Publication Date
CN107121015A true CN107121015A (zh) 2017-09-01
CN107121015B CN107121015B (zh) 2018-10-16

Family

ID=59718459

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710461448.7A Active CN107121015B (zh) 2017-06-16 2017-06-16 一种快速弹上弹道在线规划方法

Country Status (1)

Country Link
CN (1) CN107121015B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110717245A (zh) * 2019-09-03 2020-01-21 湖北航天技术研究院总体设计所 基于落角落速约束的拟滑翔弹道的设计方法
CN111859526A (zh) * 2020-06-04 2020-10-30 中国人民解放军国防科技大学 助推滑翔导弹总体参数的快速确定方法
CN114167887A (zh) * 2021-11-16 2022-03-11 湖北航天技术研究院总体设计所 飞行器落速控制方法、装置、设备及可读存储介质
CN114237287A (zh) * 2021-12-10 2022-03-25 中国空间技术研究院 基于滑翔弹能力约束的天基中制导信息应用方法
CN116576737A (zh) * 2023-05-25 2023-08-11 西安现代控制技术研究所 减小高度散布的垂直发射远程导弹转弯程序角生成方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5537909A (en) * 1995-04-17 1996-07-23 Hughes Missile System Company All-aspect bomb damage assessment system
CN103488814A (zh) * 2013-08-16 2014-01-01 北京航空航天大学 一种适用于再入飞行器姿态控制的闭环仿真系统
CN104035335A (zh) * 2014-05-27 2014-09-10 北京航空航天大学 基于高精度纵、横程解析预测方法的平稳滑翔再入制导律
CN104648695A (zh) * 2014-12-18 2015-05-27 北京控制工程研究所 一种基于倾侧角可用性的再入走廊最优规划方法
CN105223962A (zh) * 2015-10-27 2016-01-06 湖北航天技术研究院总体设计所 一种减小外形跳变时弹体过载的动态攻角限幅方法
CN106294280A (zh) * 2016-08-22 2017-01-04 麻毅威 一种弹道规划方法
WO2017014085A1 (ja) * 2015-07-17 2017-01-26 三菱重工業株式会社 航空機管理装置、航空機、及び航空機の軌道算出方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5537909A (en) * 1995-04-17 1996-07-23 Hughes Missile System Company All-aspect bomb damage assessment system
CN103488814A (zh) * 2013-08-16 2014-01-01 北京航空航天大学 一种适用于再入飞行器姿态控制的闭环仿真系统
CN104035335A (zh) * 2014-05-27 2014-09-10 北京航空航天大学 基于高精度纵、横程解析预测方法的平稳滑翔再入制导律
CN104648695A (zh) * 2014-12-18 2015-05-27 北京控制工程研究所 一种基于倾侧角可用性的再入走廊最优规划方法
WO2017014085A1 (ja) * 2015-07-17 2017-01-26 三菱重工業株式会社 航空機管理装置、航空機、及び航空機の軌道算出方法
CN105223962A (zh) * 2015-10-27 2016-01-06 湖北航天技术研究院总体设计所 一种减小外形跳变时弹体过载的动态攻角限幅方法
CN106294280A (zh) * 2016-08-22 2017-01-04 麻毅威 一种弹道规划方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110717245A (zh) * 2019-09-03 2020-01-21 湖北航天技术研究院总体设计所 基于落角落速约束的拟滑翔弹道的设计方法
CN110717245B (zh) * 2019-09-03 2022-04-29 湖北航天技术研究院总体设计所 基于落角落速约束的拟滑翔弹道的设计方法
CN111859526A (zh) * 2020-06-04 2020-10-30 中国人民解放军国防科技大学 助推滑翔导弹总体参数的快速确定方法
CN114167887A (zh) * 2021-11-16 2022-03-11 湖北航天技术研究院总体设计所 飞行器落速控制方法、装置、设备及可读存储介质
CN114167887B (zh) * 2021-11-16 2023-08-15 湖北航天技术研究院总体设计所 飞行器落速控制方法、装置、设备及可读存储介质
CN114237287A (zh) * 2021-12-10 2022-03-25 中国空间技术研究院 基于滑翔弹能力约束的天基中制导信息应用方法
CN116576737A (zh) * 2023-05-25 2023-08-11 西安现代控制技术研究所 减小高度散布的垂直发射远程导弹转弯程序角生成方法

Also Published As

Publication number Publication date
CN107121015B (zh) 2018-10-16

Similar Documents

Publication Publication Date Title
CN107121015B (zh) 一种快速弹上弹道在线规划方法
CN106681348B (zh) 考虑全捷联导引头视场约束的制导控制一体化设计方法
CN103090728B (zh) 一种基于滑模控制的带末角约束制导方法
CN107941087B (zh) 一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法
CN106997208B (zh) 一种面向不确定条件下的高超声速飞行器的控制方法
CN106020231B (zh) 基于再入点参数的高超声速飞行器再入轨迹优化方法
CN104077469B (zh) 基于速度预测的分段迭代剩余时间估计方法
CN106586033A (zh) 自适应分段的多段线性伪谱广义标控脱靶量再入制导方法
CN106444838A (zh) 一种自主水下航行器的精确航路跟踪控制方法
CN105159308A (zh) 一种可重复使用飞行器着陆段制导与控制律一体化耦合设计方法
CN111591470A (zh) 一种适应推力可调模式的飞行器精确软着陆闭环制导方法
CN110764523B (zh) 基于抗饱和平滑变换的比例积分前置导引攻击目标的方法
CN105550402B (zh) 一种基于攻角或倾侧角变频的高超平稳机动滑翔弹道设计方法
CN104317305A (zh) 一种面向复杂战场威胁的射前航迹确定方法
CN102654772A (zh) 一种基于控制力受限情况下飞行器航迹倾角反演控制方法
Bao et al. Integrated guidance and control for hypersonic morphing missile based on variable span auxiliary control
CN103744057A (zh) 基于输出相关自适应卡尔曼滤波的弹道轨迹形成方法
CN110750850A (zh) 强约束复杂任务条件下的三维剖面优化设计方法、系统及介质
CN111897223B (zh) 一种考虑自动驾驶仪动态特性的速度追踪制导方法
CN109827478B (zh) 一种带落角约束与过载约束的制导方法
CN110065649A (zh) 采用虚拟瞄准点的临近空间高超声速飞行器弹道设计方法
CN115585814A (zh) 一种基于沉降机制的飞行器变尺度地形跟随方法
CN113342044B (zh) 一种可重复使用运载器末端能量管理段地面轨迹设计方法
Raju et al. Empirical virtual sliding target guidance law design: An aerodynamic approach
CN111474948B (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