首页> 中国专利> 一种基于攻角或倾侧角变频的高超平稳机动滑翔弹道设计方法

一种基于攻角或倾侧角变频的高超平稳机动滑翔弹道设计方法

摘要

一种基于攻角或倾侧角变频的高超平稳机动滑翔弹道设计方法,它有八大步骤:一、高超声速飞行器再入过程约束建模;二、平稳滑翔攻角走廊的建立;三、平稳滑翔倾侧角走廊的建立;四、平稳机动滑翔参考攻角曲线设计;五、平稳机动滑翔的参考倾侧角曲线设计;六、初始滑翔高度和弹道倾角求解;七、定阻尼微分反馈控制方案;八、平稳机动滑翔弹道的生成。本发明提出了平稳机动滑翔的概念,它是一种新的再入飞行模式,能够兼顾平衡滑翔便于跟踪制导和跳跃滑翔突防能力强的优点,降低了弹道规划的难度。

著录项

  • 公开/公告号CN105550402A

    专利类型发明专利

  • 公开/公告日2016-05-04

    原文格式PDF

  • 申请/专利权人 北京航空航天大学;

    申请/专利号CN201510888393.9

  • 发明设计人 陈万春;胡锦川;

    申请日2015-12-07

  • 分类号G06F17/50(20060101);

  • 代理机构11232 北京慧泉知识产权代理有限公司;

  • 代理人王顺荣;唐爱华

  • 地址 100191 北京市海淀区学院路37号

  • 入库时间 2023-12-18 15:54:55

法律信息

  • 法律状态公告日

    法律状态信息

    法律状态

  • 2018-08-31

    授权

    授权

  • 2016-06-01

    实质审查的生效 IPC(主分类):G06F17/50 申请日:20151207

    实质审查的生效

  • 2016-05-04

    公开

    公开

说明书

技术领域

本发明提供了一种基于攻角或倾侧角变频的高超平稳机动滑翔弹道设计方法,属 于航天技术、武器技术领域。

背景技术

高超声速飞行器一般是指飞行速度超过5倍音速的飞机、导弹、炮弹之类的有翼或 无翼飞行器,它具有飞行高度高、速度快、侧向机动性好的优点,已成为21世纪世界航空航 天事业发展的一个主要方向。滑翔弹道规划能有效的发掘高超声速飞行器的任务潜能,同 时保证飞行过程的安全,是其中的一项关键技术。

高超声速飞行器滑翔弹道通常被分为平衡滑翔和跳跃滑翔两种飞行模式,其中平 衡滑翔具有高度变化平缓、热流密度和动压峰值小、攻角和倾侧角曲线光滑等优点,便于再 入跟踪制导;而跳跃滑翔则具有飞行距离远、突防能力强的优点。平稳滑翔是在上述基础上 新提出的一种飞行模式,对于给定的连续可导的参考攻角和倾侧角曲线,均存在着相应的 平稳滑翔弹道。因此,可通过设计变化的参考攻角和倾侧角曲线,来获得兼顾平衡滑翔和跳 跃滑翔两者优点的平稳机动滑翔弹道。

发明内容

本发明的目的是提供一种基于攻角或倾侧角变频的高超平稳机动滑翔弹道设计 方法,它利用平衡滑翔条件将再入过程约束转化成为攻角走廊(倾侧角走廊),并给出了一 种走廊内参考攻角(倾侧角)的设计方案;在此基础上给出了该参考攻角和倾侧角对应平稳 机动滑翔弹道的初始高度和弹道倾角;并采用定阻尼微分反馈消除可能存在的状态偏差; 最终利用弹道积分即可获得相应的平稳机动滑翔弹道。

本发明为一种基于攻角或倾侧角变频的高超平稳机动滑翔弹道设计方法,包含以 下几个步骤:

步骤1:高超声速飞行器再入过程约束建模

高超声速飞行器再入过程中需要考虑最大热流密度、最大动压和最大过载等过程 约束,同时攻角和倾侧角也有边界约束,具体如下:

Q·=kρV3Q·maxq=12ρV2qmax

n=L2+D2nmaxαminααmaxσmin|σ|σmax

式中,为热流密度;为最大热流密度;ρ为大气密度;V为速度;k为飞行器的 热流密度常数。q为动压;n为过载;qmax和nmax分别为最大动压和最大过载。α为攻角;σ为倾侧 角;αmin、αmax分别为最小攻角、最大攻角;σmin和σmax分别为最小倾侧角和最大倾侧角。L和D分 别为归一化的升力和阻力。

步骤2:平稳滑翔攻角走廊的建立

设滑翔段的参考倾侧角绝对值如下所示:

ref|=fσ(V)

上式中,σref为参考倾侧角,fσ(V)为速度相关函数。利用平衡滑翔条件可得最大热 流密度、最大动压和最大过载对应的攻角下界分别为αQ、αq和αn,则攻角走廊的下界为,

αlow=max(αQqnmin)

上式中,αlow为攻角走廊的下界。而攻角走廊的上界则为αup=αmax

步骤3:平稳滑翔倾侧角走廊的建立

设滑翔段的攻角曲线如下所示,

αref=fα(V)

上式中,αref为参考攻角,fα(V)为速度相关函数。利用平衡滑翔条件可得最大热流 密度、最大动压和最大过载对应的倾侧角绝对值的上界分别为|σQ|、|σq|和|σn|,则倾侧角 走廊的上界为,

σup=min(|σQ|,|σq|,|σn|,σmax)

上式中,σup为倾侧角走廊的上界。而倾侧角走廊的下界则为σlow=σmin

步骤4:平稳机动滑翔参考攻角曲线设计

本发明采用多项式函数和三角函数组合设计参考攻角曲线,如下所示:

αref=fα1(V)+Ca1sin(Ca2π(V-Vf)Ca3+Ca4)

式中,αref为参考攻角曲线;fα1(V)为以速度为自变量的多项式函数,它决定了参考 攻角的均值;Ca1为控制攻角振幅的常数,通常由攻角走廊决定其大小;Ca2为控制攻角振荡 频率的常数,是调整平稳机动滑翔飞行模态的主要参数;Ca3为速度常数,Ca3=V0-Vf,其中V0和Vf分别为滑翔弹道的初始速度和终端速度;Ca4为相位调整常数。

步骤5:平稳机动滑翔的参考倾侧角曲线设计

与步骤4类似,采用多项式函数和三角函数组合设计参考倾侧角曲线,如下所示:

σref=fσ1(V)+Cs1sin(Cs2π(V-Vf)Cs3+Cs4)

式中,σref为参考倾侧角曲线;fσ1(V)为以速度为自变量的多项式函数,它决定了参 考倾侧角的均值,进而影响了规划弹道的射程;Cs1为控制倾侧角振幅的常数,通常由倾侧角 走廊决定其大小;Cs2为控制倾侧角振荡频率的常数,是调整平稳机动滑翔飞行模态的主要 参数;Cs3为速度常数,Cs3=V0-Vf,其中V0和Vf分别为滑翔弹道的初始速度和终端速度;Cs4为 相位调整常数。

步骤6:初始滑翔高度和弹道倾角求解

根据步骤2和步骤4(或者步骤3和步骤5)确定的参考攻角和倾侧角,求解平稳机动 滑翔弹道的初始高度和弹道倾角,具体如下:

hsg=-1βln(g-V2/r-asgkcV2CL1*)

γsg=-(fa+fV)(g+asg-V2/r)K1*βV+K1*g(fa+fV)---(1)

上式中,hsg、γsg和asg分别为平稳滑翔高度、弹道倾角和纵向加速度;kc为飞行器 的组合常数,kc=ρ0S/(2m);和分别为参考攻角和倾侧角对应的升力系数纵向分量和 升阻比纵向分量;fa和fV均为速度相关函数;asg为平稳滑翔纵向加速度。

步骤7:定阻尼微分反馈控制方案

定阻尼微分反馈控制能够使得滑翔弹道快速收敛到平稳滑翔状态,具体如下:

k2=2ζc(g-V2/r)β-(g-V2/r)/(VK1*)

CL1(real)=CL1*-2mk2CL1*ρVScosσ*(γ-γsg)

上式中,ζc为给定的阻尼,通常取ζc=0.707;ρ、g、V、r、γ分别为飞行器当前的大气 密度、重力加速度、速度、地心距和弹道倾角;为实际的升力系数纵向分量;γsg为参考 弹道倾角,由式(1)获得;σ*为参考倾侧角;为由参考升力系数纵向分量;k2为微分反馈 系数。

步骤8:平稳机动滑翔弹道的生成

以步骤6获得的初始高度和弹道倾角为积分初值,以步骤7获得的实际升力系数纵 向分量为控制变量,进行弹道积分;即可获得平稳机动滑翔弹道。

本发明的优点在于:

提出了平稳机动滑翔的概念,它是一种新的再入飞行模式,能够兼顾平衡滑翔便 于跟踪制导和跳跃滑翔突防能力强的优点。

在进行平稳机动滑翔弹道设计时,能够采用攻角走廊或者倾侧角走廊来处理过程 约束,降低了弹道规划的难度。

提出了一种满足过程约束的要求的平稳机动滑翔参考攻角和倾侧角设计方案。在 此基础上,利用平稳滑翔弹道动态特性,通过弹道积分即可获得所需的平稳机动滑翔弹道。

附图说明

图1是平稳机动滑翔弹道规划流程示意图。

图2是参考倾侧角为0deg时的平稳滑翔攻角走廊示意图。

图3是不同参考倾侧角下的平稳滑翔攻角走廊示意图。

图4是参考攻角为15deg时的平稳滑翔倾侧角走廊示意图。

图5是不同参考攻角下的平稳滑翔倾侧角走廊示意图。

图6是平稳机动滑翔参考攻角曲线示意图。

图7是平稳机动滑翔参考倾侧角曲线示意图。

图8是变化攻角曲线获得的平稳机动滑翔纵向弹道示意图。

图9是变化攻角曲线获得的速度曲线示意图。

图10是变化攻角曲线获得的横向弹道示意图。

图11是变化攻角曲线获得的法向加速度示意图。

图12是变化倾侧角曲线获得的平稳机动滑翔纵向弹道示意图。

图13是变化倾侧角曲线获得的速度曲线示意图。

图14是变化倾侧角曲线获得的横向弹道示意图。

图15是变化倾侧角曲线获得的法向加速度示意图。

图中符号说明如下:

σref:为参考倾侧角;αref为参考攻角。

具体实施方案

见图1—图15,下面将结合附图和实施案例对本发明作进一步的详细说明。

本发明提出了一种平稳机动滑翔的再入弹道飞行模式;并给出了相应的弹道规划 方法。首先利用平衡滑翔条件将再入过程约束转化成为相应的攻角走廊(或者倾侧角走 廊);然后在走廊内设计参考攻角(或倾侧角)曲线;最后利用平稳滑翔弹道动态特性积分获 得所需的平稳机动滑翔弹道。

本发明一种基于攻角或倾侧角变频的高超平稳机动滑翔弹道设计方法,它包括以 下步骤:步骤1:高超声速飞行器再入过程约束建模

高超声速飞行器再入过程中需要考虑最大热流密度、最大动压和最大过载等过程 约束,同时攻角和倾侧角也有边界约束,具体如下:

Q·=kρV3Q·max---(2)

q=12ρV2qmax---(3)

n=L2+D2nmax---(4)

αmin≤α≤αmax(5)

σmin≤|σ|≤σmax(6)

式中,为热流密度;为最大热流密度;ρ为大气密度;V为速度;k为飞行器的 热流密度常数。q为动压;n为过载;qmax和nmax分别为最大动压和最大过载。α为攻角;σ为倾侧 角;αmin、αmax分别为最小攻角、最大攻角;σmin和σmax分别为最小倾侧角和最大倾侧角。L和D分 别为归一化的升力和阻力,具体如下:

L=ρV2SCL2mg0D=ρV2SCD2mg0

式中,S为飞行器的气动参考面积,m为飞行器的质量;g0为海平面重力加速度, CL和CD分别为升力系数和阻力系数。

高超声速滑翔段飞行过程中若符合平稳滑翔条件,则纵向受力如下所示,

ρV2SCLcosσ2m+(V2r-g)cosγ-aϵ=0

上式中,g为重力加速度,g=μ/r2,μ为地球引力常数,r为飞行器的地心距;aε为平 稳滑翔坡度率加速度,与参考攻角有关,为一小量;γ为弹道倾角。在进行再入走廊分析时, 通常可忽略aε的影响,并且假定cosγ≈1,则可获得平衡滑翔条件如下:

ρV2SCLcosσ2m+(V2r-g)=0---(7)

在进行滑翔弹道规划时,若指定了参考倾侧角曲线,则可利用平衡滑翔条件将式 (2)至式(6)给出的过程约束转化为攻角走廊约束;反之,若指定了参考攻角曲线,则可利用 平衡滑翔条件将式(2)至式(6)给出的过程约束转化为倾侧角走廊约束。

步骤2:平稳滑翔攻角走廊的建立

设滑翔段的参考倾侧角绝对值如下所示:

ref|=fσ(V)

上式中,σref为参考倾侧角,fσ(V)为速度相关函数。则利用式(7)可求得最大热流 密度、最大动压和最大过载对应的升力系数下界分别为,

CLQ=2m(gQ-V2/rQ)S>cosσref(kV2Q·max)2

CLq=m(gq-V2/rq)qmaxS>cosσref

CLnCLn2+CDn2=(gn-V2/rn)nmaxg0cosσref

上式中,CLQ和CLq分别为最大热流密度和最大动压对应的升力系数;CLn和CDn分别为 最大过载对应的阻力系数;k为与再入飞行器热流密度相关的常数。上述四者与当前速度和 攻角相关,采用secant法可快速求解出最大热流密度对应的攻角下界αQ、最大动压对应的 攻角下界αq和最大过载对应的攻角下界αn。gQ、gq和gn分别为最大热流密度、最大动压和最大 过载对应的重力加速度;rQ、rq和rn分别为最大热流密度、最大动压和最大过载对应的地心 距,如下所示:

rQ=R0-1βln[1ρ0(Q·maxkV3)2]

rq=R0-1βln(2qmaxp0V2)

rn=R0-1βln(2nmaxmg0ρ0V2SCLn2+CDn2)---(8)

上式中,R0为地球半径;ρ0为海平面大气密度;β为指数大气模型常数,通常取β=1/ 7200m-1。在求解rn时,可先假定rn≈R0,然后求解式(12)估算的取值,最后利用 式(8)获得rn的解。最后,可得攻角走廊的下界为,

αlow=max(αQqnmin)

上式中,αlow为攻角走廊的下界。而攻角走廊的上界则为αup=αmax

步骤3:平稳滑翔倾侧角走廊的建立

设滑翔段的攻角曲线如下所示,

αref=fα(V)(9)

上式中,αref为参考攻角,fα(V)为速度相关函数。利用式(7)可求得最大热流密度、 最大动压和最大过载对应的倾侧角上界分别为,

|σQ|=arccos[2m(gQ-V2/rQ)SCL*(kV2Q·max)2]

|σq|=arccos[m(gq-V2/rq)qmaxSCL*]

|σn|=arccos[(gn-V2/rn)1+(CD*/CL*)2nmaxg0]

上式中,和分别为参考攻角对应的升力系数和阻力系数;|σQ|、|σq|和|σn| 分别为最大热流密度、最大动压和最大过载对应的倾侧角绝对值的上边界。最终可得倾侧 角走廊的上界为,

σup=min(|σQ|,|σq|,|σn|,σmax)

上式中,σup为倾侧角走廊的上界。而倾侧角走廊的下界则为σlow=σmin

步骤4:平稳机动滑翔参考攻角曲线设计

在利用步骤2获得平稳滑翔攻角再入走廊之后,还需在再入走廊内设计机动滑翔 所需的连续可导的参考攻角曲线,才能最终确定平稳机动滑翔弹道。本发明采用多项式函 数和三角函数组合设计参考攻角曲线,如下所示:

αref=fα1(V)+Ca1sin(Ca2π(V-Vf)Ca3+Ca4)

式中,αref为参考攻角曲线;fα1(V)为以速度为自变量的多项式函数,它决定了参考 攻角的均值;Ca1为控制攻角振幅的常数,通常由攻角走廊决定其大小;Ca2为控制攻角振荡 频率的常数,是调整平稳机动滑翔飞行模态的主要参数;Ca3为速度常数,Ca3=V0-Vf,其中V0和Vf分别为滑翔弹道的初始速度和终端速度;Ca4为相位调整常数。

步骤5:平稳机动滑翔的参考倾侧角曲线设计

与步骤4类似,在步骤3确定了平稳滑翔倾侧角再入走廊的基础上,采用多项式函 数和三角函数组合设计参考倾侧角曲线,如下所示:

σref=fσ1(V)+Cs1sin(Cs2π(V-Vf)Cs3+Cs4)

式中,σref为参考倾侧角曲线;fσ1(V)为以速度为自变量的多项式函数,它决定了参 考倾侧角的均值,进而影响了规划弹道的射程;Cs1为控制倾侧角振幅的常数,通常由倾侧角 走廊决定其大小;Cs2为控制倾侧角振荡频率的常数,是调整平稳机动滑翔飞行模态的主要 参数;Cs3为速度常数,Cs3=V0-Vf,其中V0和Vf分别为滑翔弹道的初始速度和终端速度;Cs4为 相位调整常数。

步骤6:初始滑翔高度和弹道倾角求解

根据步骤2和步骤4(或者步骤3和步骤5)确定的参考攻角和倾侧角,求解平稳机动 滑翔弹道的初始高度和弹道倾角,具体如下:

hsg=-1βln(g-V2/r-asgkcV2CL1*)---(10)

γsg=-(fa+fV)(g+asg-V2/r)K1*βV+K1*g(fa+fV)---(11)

上式中,hsg、γsg和asg分别为平稳滑翔高度、弹道倾角和纵向加速度;kc为飞行器 的组合常数,kc=ρ0S/(2m);和分别为参考攻角和倾侧角对应的升力系数纵向分量和 升阻比纵向分量,表达式如下,

CL1*=CL*cosσrefK1*=CL*cosσref/CD*

fa和fV均为速度相关函数,表达式分别如下,

fa=(CL1*/V)/CL1*

fV=2(g-asg)/(gV-asgV-V3/r)

asg可采用下式进行估算:

asg=(gr-V2)(K1*)2βr[-gfaV+4g+V3fa/rV3+(g-V2r)faV-(gfa-V2fa/r+2g/VK1*)K1*V]

式中,和分别为fa和对速度的一阶偏导数。将式(10)和(11)中 的速度V取为平稳机动滑翔弹道的初始速度V0,即可获得相应的初始高度和弹道倾角。

步骤7:定阻尼微分反馈控制方案

定阻尼微分反馈控制能够使得滑翔弹道快速收敛到平稳滑翔状态,具体如下:

k2=2ζc(g-V2/r)β-(g-V2/r)/(VK1*)---(12)

CL1(real)=CL1*-2mk2CL1*ρVScosσ*(γ-γsg)---(13)

上式中,ζc为给定的阻尼,通常取ζc=0.707;ρ、g、V、r分别为飞行器当前的大气密 度、重力加速度、速度和地心距;为实际的升力系数纵向分量;γsg为参考弹道倾角, 由式(11)获得;σ*为参考倾侧角;为参考升力系数纵向分量;k2为微分反馈系数。

步骤8:平稳机动滑翔弹道的生成

以步骤6获得的初始高度和弹道倾角为积分初值,以步骤7获得的实际升力系数纵 向分量为控制变量,进行弹道积分;即可获得平稳机动滑翔弹道。

实施案例

为了检验解析求解算法的精度,选用CAV为仿真计算的飞行器模型,仿真中参数设 置如表格1。若取参考倾侧角|σref|=0deg,则攻角走廊如图2所示。可以看出,攻角走廊的上 界由最大攻角决定,而攻角走廊的下界则主要由最大热流密度约束和最大动压约束决定。 图3则给出了不同参考倾侧角对应的攻角走廊,可以看出,参考倾侧角越大,则攻角走廊的 下界越高。若取参考攻角αref=15deg,则倾侧角走廊如图4所示。可以看出,倾侧角走廊的下 界由最小倾侧角决定,而倾侧角的上界则由最大热流密度约束、最大倾侧角和最大动压决 定。图5则给出了不同参考攻角对应的倾侧角走廊,可以看出参考攻角越大则倾侧角走廊的 上界越高。表格2给出了平稳机动滑翔所需的参考攻角和参考倾侧角,它们在再入走廊中的 位置如图6和图7所示。

表格1CAV过程约束限制取值

表格2参考攻角曲线和参考倾侧角

编号 参考攻角(deg) 参考倾侧角(deg) 算例1 17-0.0004V+4sin[0.0004π(V-2000)] 45 算例2 17-0.0004V+4sin[0.0006π(V-2000)] 45 算例3 17-0.0004V+4sin[0.0008π(V-2000)] 45 算例4 17-0.0004V+4sin[0.001π(V-2000)] 45

算例5 15 45+15sin[0.0004π(V-2000)] 算例6 15 45+15sin[0.0006π(V-2000)] 算例7 15 45+15sin[0.0008π(V-2000)] 算例8 15 45+15sin[0.001π(V-2000)]

图8至图11给出了不同攻角变化频率对应的平稳机动滑翔弹道;而图12至图15给 出了不同倾侧角变化频率对应的平稳机动滑翔弹道。由图8和图12比较可知,两者均能活得 高度起伏变化的平稳滑翔纵向弹道,但由变化攻角曲线获得的结果更为平滑;由图9和图13 对比可知,变化攻角曲线对速度的影响更小;由图10和图14对比可知,变化攻角曲线对横向 弹道的影响也更小;由图11和图15对比可知,变化攻角曲线所需的法相加速度也更小。综上 可得,在进行平稳机动滑翔弹道规划时,通过在攻角走廊中规划连续可导的变化攻角是设 计平稳机动滑翔的最佳方式。

去获取专利,查看全文>

相似文献

  • 专利
  • 中文文献
  • 外文文献
获取专利

客服邮箱:kefu@zhangqiaokeyan.com

京公网安备:11010802029741号 ICP备案号:京ICP备15016152号-6 六维联合信息科技 (北京) 有限公司©版权所有
  • 客服微信

  • 服务号