CN107609307A - 一种考虑气弹和地球影响的远程火箭弹弹道分析方法 - Google Patents
一种考虑气弹和地球影响的远程火箭弹弹道分析方法 Download PDFInfo
- Publication number
- CN107609307A CN107609307A CN201710933732.XA CN201710933732A CN107609307A CN 107609307 A CN107609307 A CN 107609307A CN 201710933732 A CN201710933732 A CN 201710933732A CN 107609307 A CN107609307 A CN 107609307A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- rocket projectile
- 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
Landscapes
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
本发明公开一种考虑气弹和地球影响的远程火箭弹弹道分析方法,属于飞行器弹道领域。本发明实现方法为:对远程弹性体火箭弹进行模态分析,通过有限元前处理建模软件进行结构建模,用有限元后处理建模软件求解计算,得到模态参数;再建立坐标系,包括发射坐标系和弹体坐标系;建立火箭弹的气动弹性分析和动力学模型;考虑地球自转与燃料变化影响,将火箭弹的气动弹性分析和动力学模型进行外力分解;建立考虑气弹和地球影响的远程火箭弹的制导控制方案,提高制导精度。本发明提出完整的六自由度远程弹性体火箭弹全弹道分析方法,能够实现考虑气弹和地球影响的远程火箭弹弹道分析,解决远程火箭弹弹道领域相关实际工程问题,具有精度高的优点。
Description
技术领域
本发明涉及一种考虑气弹和地球影响的远程火箭弹弹道分析方法,属于飞行器弹道领域。
背景技术
远程火箭是陆军炮兵拥有的一类重要装备,主要包括火箭炮、火箭弹等战斗装备,侦察、指挥等信息装备,装填、运输等保障装备,可适应全地形、全天候作战,主要遂行面压制、小幅员压制及火力突击等作战任务。随着当今武器装备技术的不断发展和战场需求的不断提高,战场所需求的不仅仅是火箭弹的威力,也对火箭弹的打击精度提出了更高的要求,因此,如何建立高精度弹道仿真模型是未来重点研究方向之一。
随着远程火箭弹飞行速度等性能的不断提高,气动弹性对火箭弹动态品质的影响已成为火箭弹设计中必须考虑的问题。气动力、弹性力、惯性力以及控制系统之问的相互作用,产生了各种气动弹性问题。气动弹性不稳定性对于火箭弹结构来说是致命的,所以在火箭弹弹道仿真过程中,必须考虑气动弹性的影响。
目前对于远程弹性体火箭弹弹道分析问题,通常把地球当作平面大地,同时忽略地球曲率,地球自转以及燃料变化产生的附加力,若不考虑阵风等干扰且发射方向指向目标,就使得火箭弹空间六自由度运动变成铅垂面三自由度运动,从而造成仿真结果不准确。目前,尽管已有较多的远程火箭弹弹道分析研究,但现有工作都是从不同角度考虑了部分因素对弹道分析的影响,缺乏能够满足远程弹性体火箭弹高精度弹道仿真需求的研究。因此,有必要在弹性体火箭弹弹道仿真过程中,充分考虑地球曲率、地球自转以及燃料变化产生附加力的影响。
发明内容
为了解决考虑气弹和地球影响的远程火箭弹弹道分析问题,本发明公开一种考虑气弹和地球影响的远程火箭弹弹道分析方法,提出完整的六自由度远程弹性体火箭弹全弹道分析方法,能够实现考虑气弹和地球影响的远程火箭弹弹道分析,解决远程火箭弹弹道领域相关实际工程问题,具有精度高的优点。
本发明目的是通过下述技术方案实现。
本发明公开一种考虑气弹和地球影响的远程火箭弹弹道分析方法,首先对远程弹性体火箭弹进行模态分析,通过有限元前处理建模软件进行结构建模,用有限元后处理建模软件求解计算,得到模态参数,所述的模态参数包括模态振型和刚度矩阵;再建立坐标系,包括发射坐标系和弹体坐标系;之后建立火箭弹的气动弹性分析和动力学模型,动力学模型包括火箭弹质心运动、绕质心转动和弹性振动;然后考虑地球自转与燃料变化影响,将火箭弹的气动弹性分析和动力学模型外力分解为地球引力、推力、附加哥氏力、离心惯性力、哥氏惯性力和空气动力,其中空气动力包括定常气动力和弹性振动引起的非定常气动力;最后建立考虑气弹和地球影响的远程火箭弹的制导控制方案,提高制导精度。
本发明公开一种考虑气弹和地球影响的远程火箭弹弹道分析方法,包括如下步骤:
步骤1,远程弹性体火箭弹模态分析。
在火箭弹结构有限元前处理建模过程中,弹身采用梁模型构建,把弹翼、舵面和燃料当作非结构单元加在相应位置的有限元节点上,将燃料等规则体均匀分布在发动机位置的每个有限元节点上,将弹翼和舵面等不规则体按比例进行质量配置。远程火箭弹射程远,飞行时间长,且火箭弹主动段不断消耗燃料,导致火箭弹质量、质心位置、转动惯量、气动参数、模态振型和刚度矩阵等不断变化,即完成远程弹性体火箭弹模态分析。
作为优选的,为加快求解速度,步骤1中选取主动段的预设数量的特征点,用有限元分析软件求解计算,得到特征点时刻的模态振型和刚度矩阵,然后采用插值的方法,得到主动段每一时刻的模态振型和刚度矩阵,同理,火箭弹主动段每一时刻的质量、质心位置、转动惯量和气动参数也通过插值的方式求得。
步骤2,坐标系的建立。
发射坐标系O—xyz:坐标原点与发射点O固连,Ox轴在发射点水平面内指向发射瞄准方向,Oy轴垂直于发射点水平面指向上方,Oz轴与xOy面相垂直并构成右手坐标系。弹体坐标系O1—x1y1z1:坐标原点O1为火箭弹的质心,O1x1为火箭弹外壳对称轴,指向火箭弹的头部,O1y1在火箭弹的主对称面内,该平面在发射瞬时与发射坐标系xOy平面重合,Oy1轴垂直于Ox1轴,Oz1轴垂直于主对称面,顺着发射方向Oz1轴指向右方。
步骤3,建立远程火箭弹气动弹性分析模型。
由于弹性效应,火箭弹在主动段飞行过程中,喷口处发生弹性变形,导致推力Fp方向发生改变,火箭弹所受的推力Fp和推力矩Mp用式(1)表示
其中,Fp和Mp分别表示火箭弹所受的推力和推力矩;P为发动机推力;x(t)为推力作用点坐标;Фrz(xt)和Фry(xt)分别表示推力作用点沿y轴和z轴方向振型的转动分量。
此外,由于气动弹性和结构耦合的作用,火箭弹在飞行过程中会受到非定常气动力的影响。非定常气动力分别对弹身、弹翼和舵面三部分进行计算,其中弹身上的非定常气动力用细长体理论求解,将火箭弹沿纵向轴分成n份,得到由各分段上单位长度上的法向气动力组成向量为
其中,ΔP为火箭弹弹身单位长度法向气动力向量;ρ为空气密度;V为来流速度;s=diag(s1,s2,…,sn)为火箭弹各横截面积组成的对角矩阵。
运用细长体理论,弹身产生的非定常气动力Funst_body和力矩Munst_body的计算公式如式(3)所示
其中,E=[1 1 … 1]为n个1组成的矩阵;X=[x1 x2 … xn]为火箭弹弹身各段中点坐标矩阵;L=[l1 l2 … ln]为火箭弹弹身各段长度矩阵。
对于弹翼和舵面,采用气动导数法计算弹性变形产生的附加攻角αadd、侧滑角βadd和舵偏角δzadd、δyadd,如式(4)所示
其中,uy(xf)和uz(xf)分别表示火箭弹弹翼位置沿y轴和z轴方向的变形量;uy(xr)和uz(xr)分别表示火箭弹舵面位置沿y轴和z轴方向的变形量;Φty(xf)和Φtz(xf)分别表示火箭弹弹翼位置沿y轴和z轴方向振型的平动分量;Φty(xr)和Φtz(xr)分别表示火箭弹舵面位置沿y轴和z轴方向振型的平动分量;αadd和βadd分别表示附加攻角和附加侧滑角;δzadd和δyadd分别表示附加俯仰舵偏角和附加偏航舵偏角。
由弹翼和舵面产生的非定常气动力Funst_fin、Funst_rud和力矩Munst_fin、Munst_rud用式(5)表示
其中,和分别表示单独弹翼的法向力系数关于攻角的导数和侧向力系数关于侧滑角的导数;和分别表示单独舵面的法向力系数关于俯仰舵偏角的导数和侧向力系数关于偏航舵偏角的导数;x(f)和x(r)分别表示弹翼和舵面的坐标。
步骤4,建立远程火箭弹动力学模型。
弹性体火箭弹动力学方程写成矩阵形式,如下所示
其中,M为火箭弹质量;J为火箭弹转动惯量矩阵;v为火箭弹在发射坐标系下的质心运动速度,v=[vx vy vz]T;ω为火箭弹在弹体坐标系下的转动角速度,ω=[ωx ωyωz]T;q为弹性振动的模态坐标,q=[qz1 qz2 qz3 qy1 qy2 qy3]T,其中qz1、qz2和qz3分别表示火箭弹纵向前三阶弹性振动的模态坐标,qy1、qy2和qy3分别表示火箭弹侧向前三阶弹性振动的模态坐标;Φt为弹体模态振型的平动分量;Fall为火箭弹所受到的合外力;Mall为火箭弹所受到的合外力矩;Qall为沿火箭弹弹轴方向的分布力矩阵;Mq、Cq和Kq分别表示火箭弹结构质量矩阵、阻尼矩阵和刚度矩阵。
步骤5,对步骤4中建立的远程火箭弹动力学模型进行外力分解。
考虑地球自转与燃料变化,对作用在火箭弹上的外力进行分解,包括如下步骤:
步骤5.1,地球引力Mg。地球引力Mg用式(7)、式(8)和式(9)表示
其中,g为引力加速度矢量;g′r为引力加速度矢量在子午面内的分量;gωe为引力加速度矢量垂直于子午面的分量;r0为火箭弹与地心连线在子午面内投影的单位矢量;为火箭弹与地心连线垂直于子午面的单位矢量;f为万有引力常数;M1为地球质量;r为火箭弹距地心的距离;J为带谐系数;ae为地球赤道平均半径;φ为地心纬度。
步骤5.2,推力Fp。考虑到火箭弹喷口处的弹性变形,火箭弹所受的推力Fp和推力矩Mp用式(1)表示。
步骤5.3,附加哥氏力F′k。附加哥氏力F′k是火箭弹主动段阶段燃料的消耗,导致质心位置的变化,从而产生的一个相对力。火箭弹所受的附加哥氏力F′k用式(10)表示
其中,F′k为火箭弹所受的附加哥氏力;为燃料消耗的秒流量;ωT为火箭弹相对于惯性坐标系的转动角速度矢量;ρe为质心到喷管出口中心点的矢量。
步骤5.4,离心惯性力Fe。由于地球自转,火箭弹所受的离心惯性力Fe用式(11)表示
Fe=-Mωe×(ωe×r) (11)
其中,Fe为火箭弹所受的离心惯性力;m为火箭弹总质量;ωe为地球自转角速度矢量;r为火箭弹与地心之间的矢径。
步骤5.5,哥氏惯性力Fk。由于地球自转,火箭弹所受的哥氏惯性力Fk用式(12)表示
其中,Fk为火箭弹所受的哥氏惯性力;为火箭弹相对于发射坐标系的速度。
步骤5.6,空气动力。火箭弹受到的空气动力分为定常气动力Fst和非定常气动力Funst,其中,定常气动力Fst和力矩Mst用式(13)表示
其中,Q为动压;S为火箭弹特征面积;l为火箭弹特征长度;cA、cN和cZ分别表示火箭弹轴向力、法向力和侧向力系数;mx、my和mz分别表示火箭弹滚转力矩系数、偏航力矩系数和俯仰力矩系数。
弹身、弹翼和舵面产生的非定常气动力Funst_body、Funst_fin、Funst_rud和力矩Munst_body、Munst_fin、Munst_rud如式(3)和式(5)所示。
步骤6,根据步骤至步骤5建立考虑气弹和地球影响的远程火箭弹的制导控制方案,提高制导精度。
步骤6优选如下方法实现:
采用经典三回路驾驶仪,纵向先进行等弹道倾角爬升,中间过程保持无控飞行,末制导采用带落角约束的比例导引制导律。
俯仰方向舵偏角δz变化规律如式(14)所示,纵向等弹道倾角爬升和末制导指令信号θ*如式(15)所示
其中,kACT为舵机增益;kac为加速度计增益;KDC、KA、ωI以及kg为自驾仪设计参数;c为加速度计放置位置;t0为主动段结束时间;r0为进入末制导时的弹目距离;V为火箭弹的速度;α为攻角;为俯仰角;为俯仰角速度;θ*为纵向控制指令;θ0为火箭弹爬升段弹道倾角信号;NP、 分别表示位置项、落角项、目标机动项的导航系数;qz为弹目视线角垂直分量;为弹目视线角速度垂直分量;qF为落角约束;θ为火箭弹弹道倾角;tgo为火箭弹此刻到击中目标的时间。
火箭弹侧向采用比例导引制导律,侧向偏航方向舵偏角δy变化规律如式(16)所示,侧向末制导指令信号如式(17)所示
其中,β为侧滑角;为偏航角;为偏航角速度;为侧向控制指令;N为侧向比例导引系数;为弹目视线角速度水平分量;为弹道偏角。
式(14)、(15)、(16)、(17)即为根据步骤一至步骤五建立考虑气弹和地球影响的远程火箭弹的制导控制方案,提高制导精度。
有益效果:
1、本发明公开一种考虑气弹和地球影响的远程火箭弹弹道分析方法,建立火箭弹的气动弹性分析和动力学模型充分考虑气动弹性、地球曲率,地球自转以及燃料变化产生的附加力影响,在原有的三自由度基础上增加地球自转产生侧向运动的三个自由度,即建立完整的六自由度远程弹性体火箭弹全弹道分析模型,能够提高分析精度。
2、本发明公开一种考虑气弹和地球影响的远程火箭弹弹道分析方法,采用经典三回路驾驶仪,相比于已有技术中两回路驾驶仪能够进一步提高制导控制方案的精度。
附图说明
图1为本发明公开一种考虑气弹和地球影响的远程火箭弹弹道分析方法流程图;
图2为火箭弹结构与一维梁等效结果模型示意图;
图3为火箭弹主动段一阶振型的平动分量;
图4为火箭弹主动段二阶振型的平动分量;
图5为火箭弹主动段三阶振型的平动分量;
图6为火箭弹被动段一阶振型的平动分量;
图7为火箭弹被动段二阶振型的平动分量;
图8为火箭弹被动段三阶振型的平动分量;
图9为发射坐标系示意图;
图10为弹体坐标系示意图;
图11为刚体和弹性体火箭弹射高曲线;
图12为刚体和弹性体火箭弹侧偏曲线。
具体实施方式
为了更好地说明本发明的目的,下面通过仿真试验,结合表格和附图对本发明做进一步说明。
火箭弹基本参数定义如表1所示。
表1火箭弹基本参数定义
本实施例公开一种考虑气弹和地球影响的远程火箭弹弹道分析方法,具体的实施步骤如下:
步骤1,远程弹性体火箭弹模态分析。在火箭弹结构有限元建模过程中,弹身采用梁模型构建,把弹翼、舵面和燃料当作非结构单元加在相应位置的有限元节点上,将燃料等规则体均匀分布在发动机位置的每个有限元节点上,将弹翼和舵面等不规则体按比例进行质量配置,如图2所示。通过NASTRAN求解计算,可得到火箭弹主动段前三阶频率分别为15.7Hz,43.3Hz和70.8Hz,前三阶振型的平动分量分别如图3、图4和图5所示;被动段前三阶频率分别为22.9Hz,55.1Hz和104.8Hz,前三阶振型的平动分量分别如图6、图7和图8所示。
步骤2,坐标系的建立。发射坐标系O—xyz:坐标原点与发射点O固连,Ox轴在发射点水平面内指向发射瞄准方向,Oy轴垂直于发射点水平面指向上方,Oz轴与xOy面相垂直并构成右手坐标系,如图9所示。弹体坐标系O1—x1y1z1:坐标原点O1为火箭弹的质心,O1x1为火箭弹外壳对称轴,指向火箭弹的头部,O1y1在火箭弹的主对称面内,该平面在发射瞬时与发射坐标系xOy平面重合,Oy1轴垂直于Ox1轴,Oz1轴垂直于主对称面,顺着发射方向Oz1轴指向右方,如图10所示。
步骤3,建立远程火箭弹气动弹性分析模型。
由于弹性效应,火箭弹在主动段飞行过程中,喷口处发生弹性变形,导致推力Fp方向发生改变,火箭弹所受的推力Fp和推力矩Mp用式(18)表示
其中,Fp和Mp分别表示火箭弹所受的推力和推力矩;P为发动机推力,主动段推力为136670N,被动段推力为64462N;x(t)为推力作用点坐标,x(t)=8m;Фrz(xt)和Фry(xt)分别表示推力作用点沿y轴和z轴方向振型的转动分量。
此外,由于气动弹性和结构耦合的作用,火箭弹在飞行过程中会受到非定常气动力的影响。非定常气动力分别对弹身、弹翼和舵面三部分进行计算,其中弹身上的非定常气动力可用细长体理论求解,将火箭弹沿纵向轴分成95份,可得到由各分段上单位长度上的法向气动力组成向量为
其中,ΔP为火箭弹弹身单位长度法向气动力向量;ρ为空气密度;V为来流速度;s=diag(s1,s2,…,sn)为火箭弹各横截面积组成的对角矩阵。
运用细长体理论,弹身产生的非定常气动力Funst_body和力矩Munst_body的计算公式如式(20)所示
其中,E=[1 1 … 1]为n个1组成的矩阵;X=[x1 x2 … xn]为火箭弹弹身各段中点坐标矩阵;L=[l1 l2 … ln]为火箭弹弹身各段长度矩阵。
对于弹翼和舵面,采用气动导数法计算弹性变形产生的附加攻角αadd、侧滑角βadd和舵偏角δzadd、δyadd,如式(21)所示
其中,uy(xf)和uz(xf)分别表示火箭弹弹翼位置沿y轴和z轴方向的变形量;uy(xr)和uz(xr)分别表示火箭弹舵面位置沿y轴和z轴方向的变形量;Φty(xf)和Φtz(xf)分别表示火箭弹弹翼位置沿y轴和z轴方向振型的平动分量;Φty(xr)和Φtz(xr)分别表示火箭弹舵面位置沿y轴和z轴方向振型的平动分量;αadd和βadd分别表示附加攻角和附加侧滑角;δzadd和δyadd分别表示附加俯仰舵偏角和附加偏航舵偏角。
由弹翼和舵面产生的非定常气动力Funst_fin、Funst_rud和力矩Munst_fin、Munst_rud用式(22)表示
其中,和分别表示单独弹翼的法向力系数关于攻角的导数和侧向力系数关于侧滑角的导数;和分别表示单独舵面的法向力系数关于俯仰舵偏角的导数和侧向力系数关于偏航舵偏角的导数;x(f)和x(r)分别表示弹翼和舵面的坐标,x(f)=7.7176m x(r)=0.7880m。
步骤4,建立远程火箭弹的动力学模型。
弹性体火箭弹动力学方程写成矩阵形式,如下所示
其中,M为火箭弹质量,主动段质量为1280kg,被动段质量为560kg;J为火箭弹转动惯量矩阵,主动段沿x、y、z三个方向的惯性张量矩阵为[25.3 4526.46 4526.46]Tkg·m2,被动段沿x、y、z三个方向的惯性张量矩阵为[13.03 3084.2 3084.2]Tkg·m2;v为火箭弹在发射坐标系下的质心运动速度,v=[vx vy vz]T;ω为火箭弹在弹体坐标系下的转动角速度,ω=[ωx ωy ωz]T;q为弹性振动的模态坐标,q=[qz1 qz2 qz3 qy1 qy2 qy3]T,其中qz1、qz2和qz3分别表示火箭弹纵向前三阶弹性振动的模态坐标,qy1、qy2和qy3分别表示火箭弹侧向前三阶弹性振动的模态坐标;Φt为弹体模态振型的平动分量;Fall为火箭弹所受到的合外力;Mall为火箭弹所受到的合外力矩;Qall为沿火箭弹弹轴方向的分布力矩阵;Mq、Cq和Kq分别表示火箭弹结构质量矩阵、阻尼矩阵和刚度矩阵。
步骤5,对步骤4中建立的远程火箭弹动力学模型进行外力分解。
考虑地球自转与燃料变化,对作用在火箭弹上的外力进行分解,包括如下步骤:
步骤51,地球引力Mg。地球引力Mg可用式(24)、式(25)和式(26)表示
其中,g为引力加速度矢量;g′r为引力加速度矢量在子午面内的分量;gωe为引力加速度矢量垂直于子午面的分量;r0为火箭弹与地心连线在子午面内投影的单位矢量;为火箭弹与地心连线垂直于子午面的单位矢量;f为万有引力常数;M1为地球质量,fM1=3.986005×1014;r为火箭弹距地心的距离;J为带谐系数,J=1.5×1.08263×10-3;ae为地球赤道平均半径,ae=6378140m;φ为地心纬度。
步骤52,推力Fp。考虑到火箭弹喷口处的弹性变形,火箭弹所受的推力Fp和推力矩Mp用式(18)表示。
步骤53,附加哥氏力F′k。附加哥氏力F′k是火箭弹主动段阶段燃料的消耗,导致质心位置的变化,从而产生的一个相对力。火箭弹所受的附加哥氏力F′k用式(27)表示
其中,F′k为火箭弹所受的附加哥氏力;为燃料消耗的秒流量,=26.7kg/s;ωT为火箭弹相对于惯性坐标系的转动角速度矢量;ρe为质心到喷管出口中心点的矢量。
步骤54,离心惯性力Fe。由于地球自转,火箭弹所受的离心惯性力Fe用式(28)表示
Fe=-Mωe×(ωe×r) (28)
其中,Fe为火箭弹所受的离心惯性力;m为火箭弹总质量;ωe为地球自转角速度矢量,ωe=7.292115×10-5rad/s;r为火箭弹与地心之间的矢径。
步骤55,哥氏惯性力Fk。由于地球自转,火箭弹所受的哥氏惯性力Fk用式(29)表示
其中,Fk为火箭弹所受的哥氏惯性力;为火箭弹相对于发射坐标系的速度。
步骤56,空气动力。火箭弹受到的空气动力分为定常气动力Fst和非定常气动力Funst,其中,定常气动力Fst和力矩Mst用式(30)表示
其中,Q为动压;S为火箭弹特征面积,S=0.1075m2;l为火箭弹特征长度,l=8m;cA、cN和cZ分别表示火箭弹轴向力、法向力和侧向力系数;mx、my和mz分别表示火箭弹滚转力矩系数、偏航力矩系数和俯仰力矩系数。
弹身、弹翼和舵面产生的非定常气动力Funst_body、Funst_fin、Funst_rud和力矩Munst_body、Munst_fin、Munst_rud如式(20)和式(22)所示。
步骤6,根据步骤一至步骤5建立考虑气弹和地球影响的远程火箭弹的制导控制方案,提高制导精度。
步骤6选如下方法实现:
采用经典三回路驾驶仪,纵向先进行等弹道倾角爬升,中间过程保持无控飞行,末制导采用带落角约束的比例导引制导律。
火箭弹初始发射位置为北纬37°,东经0°,初始发射速度为0,发射方向指向北,目标在发射坐标系中的坐标为(280000,-6163,0)m。采用经典三回路驾驶仪,纵向先进行等弹道倾角爬升,中间过程保持无控飞行,末制导采用带落角约束的比例导引制导律。
俯仰方向舵偏角δz变化规律如式(31)所示,纵向等弹道倾角爬升和末制导指令信号θ*如式(32)所示
其中,kACT为舵机增益,kACT=1;kac为加速度计增益,kac=1;KDC、KA、ωI以及kg为自驾仪设计参数;c为加速度计放置位置,c=0;t0为主动段结束时间,t0=21.54s;r0为进入末制导时的弹目距离,r0=40000m;V为火箭弹的速度;α为攻角;为俯仰角;为俯仰角速度;θ*为纵向控制指令;θ0为火箭弹爬升段弹道倾角信号,θ0=30.7°;NP、分别表示位置项、落角项、目标机动项的导航系数;qz为弹目视线角垂直分量;为弹目视线角速度垂直分量;qF为落角约束;θ为火箭弹弹道倾角;tgo为火箭弹此刻到击中目标的时间。
火箭弹侧向采用比例导引制导律,侧向偏航方向舵偏角δy变化规律如式(33)所示,侧向末制导指令信号如式(34)所示
其中,β为侧滑角;为偏航角;为偏航角速度;为侧向控制指令;N为侧向比例导引系数,N=3;为弹目视线角速度水平分量;为弹道偏角。
经过上述步骤的操作,能够得到刚体火箭弹和弹性体火箭弹脱靶量,如表2所示,同时可得到火箭弹的射高和侧偏曲线分别如图11和图12所示。从图11和图12可以看出,在弹道仿真过程中,充分考虑了气弹、地球曲率、地球自转以及燃料变化产生的附加力影响,刚体和弹性体火箭弹射高和侧偏曲线有明显偏差,且刚体火箭弹脱靶量为11.6m,弹性体火箭弹脱靶量为8.5m,说明刚体和弹性体火箭弹脱靶量相近,在初步设计中,火箭弹的结构刚度是足够的。本发明所提出的一种考虑气弹和地球影响的远程火箭弹弹道分析方法,模型准确,能够为火箭弹在真实场景中能否精确命中目标提供依据。
表2刚体和弹性体火箭弹脱靶量
模型 | 脱靶量 |
刚体火箭弹 | 11.6m |
弹性体火箭弹 | 8.5m |
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例,用于解释本发明,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种考虑气弹和地球影响的远程火箭弹弹道分析方法,其特征在于:包括如下步骤:
步骤1,远程弹性体火箭弹模态分析;
在火箭弹结构有限元前处理建模过程中,弹身采用梁模型构建,把弹翼、舵面和燃料当作非结构单元加在相应位置的有限元节点上,将燃料等规则体均匀分布在发动机位置的每个有限元节点上,将弹翼和舵面等不规则体按比例进行质量配置;远程火箭弹射程远,飞行时间长,且火箭弹主动段不断消耗燃料,导致火箭弹质量、质心位置、转动惯量、气动参数、模态振型和刚度矩阵不断变化,即完成远程弹性体火箭弹模态分析;
步骤2,坐标系的建立;
发射坐标系O—xyz:坐标原点与发射点O固连,Ox轴在发射点水平面内指向发射瞄准方向,Oy轴垂直于发射点水平面指向上方,Oz轴与xOy面相垂直并构成右手坐标系;弹体坐标系O1—x1y1z1:坐标原点O1为火箭弹的质心,O1x1为火箭弹外壳对称轴,指向火箭弹的头部,O1y1在火箭弹的主对称面内,该平面在发射瞬时与发射坐标系xOy平面重合,Oy1轴垂直于Ox1轴,Oz1轴垂直于主对称面,顺着发射方向Oz1轴指向右方;
步骤3,建立远程火箭弹气动弹性分析模型;
由于弹性效应,火箭弹在主动段飞行过程中,喷口处发生弹性变形,导致推力Fp方向发生改变,火箭弹所受的推力Fp和推力矩Mp用式(1)表示
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>P</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P&Phi;</mi>
<mrow>
<mi>r</mi>
<mi>z</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>q</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>P&Phi;</mi>
<mrow>
<mi>r</mi>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>q</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>M</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P&Phi;</mi>
<mrow>
<mi>r</mi>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>qx</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P&Phi;</mi>
<mrow>
<mi>r</mi>
<mi>z</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>qx</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Fp和Mp分别表示火箭弹所受的推力和推力矩;P为发动机推力;x(t)为推力作用点坐标;Фrz(xt)和Фry(xt)分别表示推力作用点沿y轴和z轴方向振型的转动分量;
此外,由于气动弹性和结构耦合的作用,火箭弹在飞行过程中会受到非定常气动力的影响;非定常气动力分别对弹身、弹翼和舵面三部分进行计算,其中弹身上的非定常气动力用细长体理论求解,将火箭弹沿纵向轴分成n份,得到由各分段上单位长度上的法向气动力组成向量为
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>P</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;p</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&Delta;p</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&Delta;p</mi>
<mi>n</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mi>&rho;</mi>
<mo>(</mo>
<mi>s</mi>
<mi>&Phi;</mi>
<mover>
<mi>q</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mo>+</mo>
<mn>2</mn>
<msup>
<mi>Vs&Phi;</mi>
<mo>&prime;</mo>
</msup>
<mover>
<mi>q</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>+</mo>
<msup>
<mi>Vs</mi>
<mo>&prime;</mo>
</msup>
<mi>&Phi;</mi>
<mover>
<mi>q</mi>
<mo>&CenterDot;</mo>
</mover>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msup>
<mi>V</mi>
<mn>2</mn>
</msup>
<msup>
<mi>s</mi>
<mo>&prime;</mo>
</msup>
<msup>
<mi>&Phi;</mi>
<mo>&prime;</mo>
</msup>
<mi>q</mi>
<mo>+</mo>
<msup>
<mi>V</mi>
<mn>2</mn>
</msup>
<msup>
<mi>s&Phi;</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
<mi>q</mi>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,ΔP为火箭弹弹身单位长度法向气动力向量;ρ为空气密度;V为来流速度;s=diag(s1,s2,…,sn)为火箭弹各横截面积组成的对角矩阵;
运用细长体理论,弹身产生的非定常气动力Funst_body和力矩Munst_body的计算公式如式(3)所示
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>s</mi>
<mi>t</mi>
<mo>_</mo>
<mi>b</mi>
<mi>o</mi>
<mi>d</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mi>&Delta;</mi>
<mi>P</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>s</mi>
<mi>t</mi>
<mo>_</mo>
<mi>b</mi>
<mi>o</mi>
<mi>d</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mi>&Delta;</mi>
<mi>P</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,E=[1 1 ... 1]为n个1组成的矩阵;X=[x1 x2 ... xn]为火箭弹弹身各段中点坐标矩阵;L=[l1 l2 ... ln]为火箭弹弹身各段长度矩阵;
对于弹翼和舵面,采用气动导数法计算弹性变形产生的附加攻角αadd、侧滑角βadd和舵偏角δzadd、δyadd,如式(4)所示
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>y</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>&Phi;</mi>
<mrow>
<mi>t</mi>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>q</mi>
<mo>,</mo>
<msub>
<mi>u</mi>
<mi>z</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>&Phi;</mi>
<mrow>
<mi>t</mi>
<mi>z</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>q</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>y</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>&Phi;</mi>
<mrow>
<mi>t</mi>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>q</mi>
<mo>,</mo>
<msub>
<mi>u</mi>
<mi>z</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>&Phi;</mi>
<mrow>
<mi>t</mi>
<mi>z</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>q</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>a</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mi>y</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>y</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mi>V</mi>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>a</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mi>z</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>z</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mi>V</mi>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>z</mi>
<mi>a</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mi>y</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>y</mi>
<mi>a</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mi>z</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,uy(xf)和uz(xf)分别表示火箭弹弹翼位置沿y轴和z轴方向的变形量;uy(xr)和uz(xr)分别表示火箭弹舵面位置沿y轴和z轴方向的变形量;Φty(xf)和Φtz(xf)分别表示火箭弹弹翼位置沿y轴和z轴方向振型的平动分量;Φty(xr)和Φtz(xr)分别表示火箭弹舵面位置沿y轴和z轴方向振型的平动分量;αadd和βadd分别表示附加攻角和附加侧滑角;δzadd和δyadd分别表示附加俯仰舵偏角和附加偏航舵偏角;
由弹翼和舵面产生的非定常气动力Funst_fin、Funst_rud和力矩Munst_fin、Munst_rud用式(5)表示
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>s</mi>
<mi>t</mi>
<mo>_</mo>
<mi>f</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>c</mi>
<mi>N</mi>
<mi>&alpha;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>a</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>c</mi>
<mi>Z</mi>
<mi>&beta;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>a</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>Q</mi>
<mi>S</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>s</mi>
<mi>t</mi>
<mo>_</mo>
<mi>r</mi>
<mi>u</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>c</mi>
<mi>N</mi>
<msub>
<mi>&delta;</mi>
<mi>z</mi>
</msub>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>z</mi>
<mi>a</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>c</mi>
<mi>Z</mi>
<msub>
<mi>&delta;</mi>
<mi>y</mi>
</msub>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>y</mi>
<mi>a</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>Q</mi>
<mi>S</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>s</mi>
<mi>t</mi>
<mo>_</mo>
<mi>f</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>c</mi>
<mi>Z</mi>
<mi>&beta;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>a</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>c</mi>
<mi>N</mi>
<mi>&alpha;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>a</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<msup>
<mi>QSx</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>s</mi>
<mi>t</mi>
<mo>_</mo>
<mi>r</mi>
<mi>u</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>c</mi>
<mi>Z</mi>
<msub>
<mi>&delta;</mi>
<mi>y</mi>
</msub>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>y</mi>
<mi>a</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>c</mi>
<mi>N</mi>
<msub>
<mi>&delta;</mi>
<mi>z</mi>
</msub>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>z</mi>
<mi>a</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<msup>
<mi>QSx</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,和分别表示单独弹翼的法向力系数关于攻角的导数和侧向力系数关于侧滑角的导数;和分别表示单独舵面的法向力系数关于俯仰舵偏角的导数和侧向力系数关于偏航舵偏角的导数;x(f)和x(r)分别表示弹翼和舵面的坐标;
步骤4,建立远程火箭弹动力学模型;
弹性体火箭弹动力学方程写成矩阵形式,如下所示
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>M</mi>
<mover>
<mi>v</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>a</mi>
<mi>l</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>J</mi>
<mover>
<mi>&omega;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>M</mi>
<mrow>
<mi>a</mi>
<mi>l</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&omega;</mi>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mi>J</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>M</mi>
<mi>q</mi>
</msub>
<mover>
<mi>q</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mo>+</mo>
<msub>
<mi>C</mi>
<mi>q</mi>
</msub>
<mover>
<mi>q</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>+</mo>
<msub>
<mi>K</mi>
<mi>q</mi>
</msub>
<mi>q</mi>
<mo>=</mo>
<msup>
<msub>
<mi>&Phi;</mi>
<mi>t</mi>
</msub>
<mi>T</mi>
</msup>
<msub>
<mi>Q</mi>
<mrow>
<mi>a</mi>
<mi>l</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,M为火箭弹质量;J为火箭弹转动惯量矩阵;v为火箭弹在发射坐标系下的质心运动速度,v=[vx vy vz]T;ω为火箭弹在弹体坐标系下的转动角速度,ω=[ωx ωy ωz]T;q为弹性振动的模态坐标,q=[qz1 qz2 qz3 qy1 qy2 qy3]T,其中qz1、qz2和qz3分别表示火箭弹纵向前三阶弹性振动的模态坐标,qy1、qy2和qy3分别表示火箭弹侧向前三阶弹性振动的模态坐标;Φt为弹体模态振型的平动分量;Fall为火箭弹所受到的合外力;Mall为火箭弹所受到的合外力矩;Qall为沿火箭弹弹轴方向的分布力矩阵;Mq、Cq和Kq分别表示火箭弹结构质量矩阵、阻尼矩阵和刚度矩阵;
步骤5,对步骤4中建立的远程火箭弹动力学模型进行外力分解;
步骤6,根据步骤一至步骤5建立考虑气弹和地球影响的远程火箭弹的制导控制方案,提高制导精度。
2.如权利要求1所述的一种考虑气弹和地球影响的远程火箭弹弹道分析方法,其特征在于:步骤5具体实现方法为,
考虑地球自转与燃料变化,对作用在火箭弹上的外力进行分解,包括如下步骤,
步骤5.1,地球引力Mg;地球引力Mg用式(7)、式(8)和式(9)表示
<mrow>
<mi>M</mi>
<mi>g</mi>
<mo>=</mo>
<msubsup>
<mi>Mg</mi>
<mi>r</mi>
<mo>&prime;</mo>
</msubsup>
<msup>
<mi>r</mi>
<mn>0</mn>
</msup>
<mo>+</mo>
<msub>
<mi>Mg</mi>
<mrow>
<mi>&omega;</mi>
<mi>e</mi>
</mrow>
</msub>
<msubsup>
<mi>&omega;</mi>
<mi>e</mi>
<mn>0</mn>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>g</mi>
<mi>r</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>fM</mi>
<mn>1</mn>
</msub>
</mrow>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<mi>J</mi>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>a</mi>
<mi>e</mi>
</msub>
<mi>r</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mn>5</mn>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>g</mi>
<mrow>
<mi>&omega;</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<mrow>
<msub>
<mi>fM</mi>
<mn>1</mn>
</msub>
</mrow>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
</mfrac>
<mi>J</mi>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>a</mi>
<mi>e</mi>
</msub>
<mi>r</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&phi;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,g为引力加速度矢量;g′r为引力加速度矢量在子午面内的分量;gωe为引力加速度矢量垂直于子午面的分量;r0为火箭弹与地心连线在子午面内投影的单位矢量;为火箭弹与地心连线垂直于子午面的单位矢量;f为万有引力常数;M1为地球质量;r为火箭弹距地心的距离;J为带谐系数;ae为地球赤道平均半径;φ为地心纬度;
步骤5.2,推力Fp;考虑到火箭弹喷口处的弹性变形,火箭弹所受的推力Fp和推力矩Mp用式(1)表示;
步骤5.3,附加哥氏力Fk′;附加哥氏力Fk′是火箭弹主动段阶段燃料的消耗,导致质心位置的变化,从而产生的一个相对力;火箭弹所受的附加哥氏力Fk′用式(10)表示
<mrow>
<msubsup>
<mi>F</mi>
<mi>k</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mo>-</mo>
<mn>2</mn>
<mover>
<mi>M</mi>
<mo>&CenterDot;</mo>
</mover>
<msub>
<mi>&omega;</mi>
<mi>T</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>&rho;</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Fk′为火箭弹所受的附加哥氏力;为燃料消耗的秒流量;ωT为火箭弹相对于惯性坐标系的转动角速度矢量;ρe为质心到喷管出口中心点的矢量;
步骤5.4,离心惯性力Fe;由于地球自转,火箭弹所受的离心惯性力Fe用式(11)表示
Fe=-Mωe×(ωe×r) (11)
其中,Fe为火箭弹所受的离心惯性力;m为火箭弹总质量;ωe为地球自转角速度矢量;r为火箭弹与地心之间的矢径;
步骤5.5,哥氏惯性力Fk;由于地球自转,火箭弹所受的哥氏惯性力Fk用式(12)表示
<mrow>
<msub>
<mi>F</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>M&omega;</mi>
<mi>e</mi>
</msub>
<mo>&times;</mo>
<mfrac>
<mrow>
<mi>&delta;</mi>
<mi>r</mi>
</mrow>
<mrow>
<mi>&delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Fk为火箭弹所受的哥氏惯性力;为火箭弹相对于发射坐标系的速度;
步骤5.6,空气动力;火箭弹受到的空气动力分为定常气动力Fst和非定常气动力Funst,其中,定常气动力Fst和力矩Mst用式(13)表示
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>l</mi>
<mi>b</mi>
</mrow>
</msub>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>c</mi>
<mi>A</mi>
</msub>
</mrow>
</mtd>
<mtd>
<msub>
<mi>c</mi>
<mi>N</mi>
</msub>
</mtd>
<mtd>
<msub>
<mi>c</mi>
<mi>Z</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mi>Q</mi>
<mi>S</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
</mrow>
</mtd>
<mtd>
<msub>
<mi>m</mi>
<mi>y</mi>
</msub>
</mtd>
<mtd>
<msub>
<mi>m</mi>
<mi>z</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mi>Q</mi>
<mi>S</mi>
<mi>l</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Q为动压;S为火箭弹特征面积;l为火箭弹特征长度;cA、cN和cZ分别表示火箭弹轴向力、法向力和侧向力系数;mx、my和mz分别表示火箭弹滚转力矩系数、偏航力矩系数和俯仰力矩系数;
弹身、弹翼和舵面产生的非定常气动力Funst_body、Funst_fin、Funst_rud和力矩Munst_body、Munst_fin、Munst_rud如式(3)和式(5)所示。
3.如权利要求1或2所述的一种考虑气弹和地球影响的远程火箭弹弹道分析方法,其特征在于:步骤6具体实现方法为,
采用经典三回路驾驶仪,纵向先进行等弹道倾角爬升,中间过程保持无控飞行,末制导采用带落角约束的比例导引制导律;
俯仰方向舵偏角δz变化规律如式(14)所示,纵向等弹道倾角爬升和末制导指令信号θ*如式(15)所示
<mrow>
<msup>
<mi>&theta;</mi>
<mo>*</mo>
</msup>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo><</mo>
<mi>t</mi>
<mo><</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mi>P</mi>
</msub>
<mi>V</mi>
<msub>
<mover>
<mi>q</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>z</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>N</mi>
<mi>&theta;</mi>
</msub>
<mi>V</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mi>z</mi>
</msub>
<mo>-</mo>
<msub>
<mi>q</mi>
<mi>F</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>t</mi>
<mrow>
<mi>g</mi>
<mi>o</mi>
</mrow>
</msub>
</mfrac>
<mo>+</mo>
<msub>
<mi>N</mi>
<msub>
<mi>a</mi>
<mi>T</mi>
</msub>
</msub>
<mi>g</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>)</mo>
<mo>/</mo>
<mi>V</mi>
<mo>+</mo>
<mi>&theta;</mi>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>r</mi>
<mo><</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,kACT为舵机增益;kac为加速度计增益;KDC、KA、ωI以及kg为自驾仪设计参数;c为加速度计放置位置;t0为主动段结束时间;r0为进入末制导时的弹目距离;V为火箭弹的速度;α为攻角;为俯仰角;为俯仰角速度;θ*为纵向控制指令;θ0为火箭弹爬升段弹道倾角信号;NP、分别表示位置项、落角项、目标机动项的导航系数;qz为弹目视线角垂直分量;为弹目视线角速度垂直分量;qF为落角约束;θ为火箭弹弹道倾角;tgo为火箭弹此刻到击中目标的时间;
火箭弹侧向采用比例导引制导律,侧向偏航方向舵偏角δy变化规律如式(16)所示,侧向末制导指令信号如式(17)所示
其中,β为侧滑角;为偏航角;为偏航角速度;为侧向控制指令;N为侧向比例导引系数;为弹目视线角速度水平分量;为弹道偏角;
式(14)、(15)、(16)、(17)即为根据步骤一至步骤五建立考虑气弹和地球影响的远程火箭弹的制导控制方案,提高制导精度。
4.如权利要求3所述的一种考虑气弹和地球影响的远程火箭弹弹道分析方法,其特征在于:为加快求解速度,步骤1中选取主动段的预设数量的特征点,用有限元分析软件求解计算,得到特征点时刻的模态振型和刚度矩阵,然后采用插值的方法,得到主动段每一时刻的模态振型和刚度矩阵,同理,火箭弹主动段每一时刻的质量、质心位置、转动惯量和气动参数也通过插值的方式求得。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710933732.XA CN107609307B (zh) | 2017-10-10 | 2017-10-10 | 一种考虑气弹和地球影响的远程火箭弹弹道分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710933732.XA CN107609307B (zh) | 2017-10-10 | 2017-10-10 | 一种考虑气弹和地球影响的远程火箭弹弹道分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107609307A true CN107609307A (zh) | 2018-01-19 |
CN107609307B CN107609307B (zh) | 2019-09-10 |
Family
ID=61068148
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710933732.XA Expired - Fee Related CN107609307B (zh) | 2017-10-10 | 2017-10-10 | 一种考虑气弹和地球影响的远程火箭弹弹道分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107609307B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109165411A (zh) * | 2018-07-26 | 2019-01-08 | 西安航天动力技术研究所 | 采用偏置且斜切结构喷管的固体发动机内弹道计算方法 |
CN109343563A (zh) * | 2018-10-15 | 2019-02-15 | 北京理工大学 | 考虑舵机失效及落角约束的飞行器制导系统及方法 |
CN110162818A (zh) * | 2018-11-26 | 2019-08-23 | 北京机电工程研究所 | 伞弹系统弹道计算方法 |
CN110516410A (zh) * | 2019-10-21 | 2019-11-29 | 湖南德雅坤创科技有限公司 | 一种飞行器总体设计参数的确定方法及装置 |
CN110717239A (zh) * | 2019-07-30 | 2020-01-21 | 蓝箭航天空间科技股份有限公司 | 一种箭(弹)体分布气动特性计算方法 |
CN111578793A (zh) * | 2020-05-07 | 2020-08-25 | 北京星途探索科技有限公司 | 一种有风时火箭整流罩分离的侧滑角控制方法 |
CN111985169A (zh) * | 2020-08-26 | 2020-11-24 | 北京中科宇航技术有限公司 | 一种近轴对称飞行器气动特性建模方法 |
CN114415509A (zh) * | 2022-01-11 | 2022-04-29 | 西北工业大学 | 一种针对远程火箭弹的伺服弹性主动抑制方法 |
CN115469546A (zh) * | 2022-09-20 | 2022-12-13 | 北京理工大学 | 大长细比火箭弹视线角速率提取方法 |
CN115615261A (zh) * | 2022-09-20 | 2023-01-17 | 北京理工大学 | 大长细比火箭弹弹性辨识与视线角速率提取融合方法 |
CN116362163A (zh) * | 2023-06-01 | 2023-06-30 | 西安现代控制技术研究所 | 一种非奇异多约束弹道快速优化方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120211593A1 (en) * | 2008-11-12 | 2012-08-23 | General Dynamics Ordnance And Tactical Systems, Inc. | Trajectory modification of a spinning projectile |
CN103604316A (zh) * | 2013-11-22 | 2014-02-26 | 北京机械设备研究所 | 一种用于多弹发射的弹道校正方法 |
CN104750948A (zh) * | 2015-04-22 | 2015-07-01 | 北京理工大学 | 一种处理飞行器设计中多极值多约束问题的优化方法 |
CN105631125A (zh) * | 2015-12-28 | 2016-06-01 | 北京理工大学 | 基于降阶模型的气动-热-结构耦合分析方法 |
CN105629725A (zh) * | 2014-10-31 | 2016-06-01 | 北京临近空间飞行器系统工程研究所 | 一种后缘舵滑翔飞行器的弹性运动建模方法 |
CN106021786A (zh) * | 2016-05-31 | 2016-10-12 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机六自由度非线性方程组配平方法 |
CN106294280A (zh) * | 2016-08-22 | 2017-01-04 | 麻毅威 | 一种弹道规划方法 |
-
2017
- 2017-10-10 CN CN201710933732.XA patent/CN107609307B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120211593A1 (en) * | 2008-11-12 | 2012-08-23 | General Dynamics Ordnance And Tactical Systems, Inc. | Trajectory modification of a spinning projectile |
CN103604316A (zh) * | 2013-11-22 | 2014-02-26 | 北京机械设备研究所 | 一种用于多弹发射的弹道校正方法 |
CN105629725A (zh) * | 2014-10-31 | 2016-06-01 | 北京临近空间飞行器系统工程研究所 | 一种后缘舵滑翔飞行器的弹性运动建模方法 |
CN104750948A (zh) * | 2015-04-22 | 2015-07-01 | 北京理工大学 | 一种处理飞行器设计中多极值多约束问题的优化方法 |
CN105631125A (zh) * | 2015-12-28 | 2016-06-01 | 北京理工大学 | 基于降阶模型的气动-热-结构耦合分析方法 |
CN106021786A (zh) * | 2016-05-31 | 2016-10-12 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机六自由度非线性方程组配平方法 |
CN106294280A (zh) * | 2016-08-22 | 2017-01-04 | 麻毅威 | 一种弹道规划方法 |
Non-Patent Citations (1)
Title |
---|
刘莉等: "基于物理规划的弹道多目标优化", 《北京理工大学学报》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109165411A (zh) * | 2018-07-26 | 2019-01-08 | 西安航天动力技术研究所 | 采用偏置且斜切结构喷管的固体发动机内弹道计算方法 |
CN109343563A (zh) * | 2018-10-15 | 2019-02-15 | 北京理工大学 | 考虑舵机失效及落角约束的飞行器制导系统及方法 |
CN110162818A (zh) * | 2018-11-26 | 2019-08-23 | 北京机电工程研究所 | 伞弹系统弹道计算方法 |
CN110717239A (zh) * | 2019-07-30 | 2020-01-21 | 蓝箭航天空间科技股份有限公司 | 一种箭(弹)体分布气动特性计算方法 |
CN110717239B (zh) * | 2019-07-30 | 2020-07-03 | 蓝箭航天空间科技股份有限公司 | 一种箭体或弹体分布气动特性计算方法 |
CN110516410A (zh) * | 2019-10-21 | 2019-11-29 | 湖南德雅坤创科技有限公司 | 一种飞行器总体设计参数的确定方法及装置 |
CN110516410B (zh) * | 2019-10-21 | 2020-01-31 | 湖南德雅坤创科技有限公司 | 一种飞行器总体设计参数的确定方法及装置 |
CN111578793B (zh) * | 2020-05-07 | 2022-08-23 | 北京星途探索科技有限公司 | 一种有风时火箭整流罩分离的侧滑角控制方法 |
CN111578793A (zh) * | 2020-05-07 | 2020-08-25 | 北京星途探索科技有限公司 | 一种有风时火箭整流罩分离的侧滑角控制方法 |
CN111985169A (zh) * | 2020-08-26 | 2020-11-24 | 北京中科宇航技术有限公司 | 一种近轴对称飞行器气动特性建模方法 |
CN111985169B (zh) * | 2020-08-26 | 2024-04-02 | 北京中科宇航技术有限公司 | 一种近轴对称飞行器气动特性建模方法 |
CN114415509A (zh) * | 2022-01-11 | 2022-04-29 | 西北工业大学 | 一种针对远程火箭弹的伺服弹性主动抑制方法 |
CN114415509B (zh) * | 2022-01-11 | 2023-04-28 | 西北工业大学 | 一种针对远程火箭弹的伺服弹性主动抑制方法 |
CN115469546A (zh) * | 2022-09-20 | 2022-12-13 | 北京理工大学 | 大长细比火箭弹视线角速率提取方法 |
CN115615261A (zh) * | 2022-09-20 | 2023-01-17 | 北京理工大学 | 大长细比火箭弹弹性辨识与视线角速率提取融合方法 |
CN116362163A (zh) * | 2023-06-01 | 2023-06-30 | 西安现代控制技术研究所 | 一种非奇异多约束弹道快速优化方法 |
CN116362163B (zh) * | 2023-06-01 | 2023-09-01 | 西安现代控制技术研究所 | 一种多约束弹道快速优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107609307B (zh) | 2019-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107609307B (zh) | 一种考虑气弹和地球影响的远程火箭弹弹道分析方法 | |
Petrescu et al. | History of aviation-a short review | |
CN110471456B (zh) | 高超声速飞行器俯冲段制导、姿控、变形一体化控制方法 | |
CN101625571B (zh) | 一种模拟自旋式飞行器六自由度运动的方法 | |
CN111721291A (zh) | 一种发射系下捷联惯组导航的工程算法 | |
CN103631153B (zh) | 可视化双模制导航空时敏炸弹半实物仿真方法 | |
CN108803639A (zh) | 一种基于反步法的四旋翼飞行器飞行控制方法 | |
CN106444836B (zh) | 一种无控探空火箭抗干扰设计方法 | |
CN109703768B (zh) | 一种基于姿态/轨迹复合控制的软式空中加油对接方法 | |
CN110334439B (zh) | 一种利用气动力辅助的超轨道速度飞行技术方法 | |
CN114756959A (zh) | 一种飞行器近距空战机动智能决策机模型设计方法 | |
CN114509946A (zh) | 一种基于预设时间滑模的飞行器制导控制一体化设计方法 | |
CN111176315B (zh) | 一种基于l1自适应控制的变体巡飞弹控制方法 | |
CN111008488B (zh) | 一种螺旋桨无人机发射过程反扭矩建模方法 | |
CN116432310A (zh) | 一种六自由度非完全可观空战机动智能决策模型设计方法 | |
CN106125571A (zh) | 一种巡航导弹的速度控制建模方法 | |
CN109634110B (zh) | 一种基于非线性干扰观测器的尾座式无人机过渡段鲁棒控制方法 | |
CN109445283B (zh) | 一种用于欠驱动浮空器在平面上定点跟踪的控制方法 | |
CN115390590B (zh) | 一种轴对称飞行器大机动控制方法及相关设备 | |
Yao et al. | Trajectory tracking controller based on PID-NLADRC | |
Baba et al. | Air combat guidance law for an UCAV | |
Fu et al. | Partial integrated guidance and control method for the interception of nearspace hypersonic target | |
Goncharenko et al. | A method to maneuver aeroballistic missiles under counteractions | |
Wang et al. | Research and Analysis on modeling method of Tandem wing UVA | |
CN116296218A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190910 Termination date: 20211010 |