CN106125571B - 一种巡航导弹的速度控制建模方法 - Google Patents
一种巡航导弹的速度控制建模方法 Download PDFInfo
- Publication number
- CN106125571B CN106125571B CN201610157141.3A CN201610157141A CN106125571B CN 106125571 B CN106125571 B CN 106125571B CN 201610157141 A CN201610157141 A CN 201610157141A CN 106125571 B CN106125571 B CN 106125571B
- Authority
- CN
- China
- Prior art keywords
- model
- missile
- speed
- speed control
- angle
- 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 27
- 238000013461 design Methods 0.000 claims abstract description 15
- 230000033001 locomotion Effects 0.000 claims description 34
- 230000003068 static effect Effects 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 8
- 238000002485 combustion reaction Methods 0.000 claims description 6
- 238000012795 verification Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 2
- 239000007789 gas Substances 0.000 claims 5
- 238000005259 measurement Methods 0.000 claims 2
- 239000002737 fuel gas Substances 0.000 claims 1
- 238000011156 evaluation Methods 0.000 abstract description 3
- 238000004080 punching Methods 0.000 abstract 3
- 238000010586 diagram Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000009194 climbing Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B17/00—Systems involving the use of models or simulators of said systems
- G05B17/02—Systems involving the use of models or simulators of said systems electric
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提出的巡航导弹速度控制建模方法,结合冲压发动机的动态特性的描述,给出了速度控制回路的组成及建模流程,具体包括冲压发动机的建模、面向仿真的速度回路全量模型、面向控制的小扰动模型,为巡航导弹的性能评估和控制器设计奠定了基础。本发明提出的面向仿真的速度控制回路建模方法,有效反映了弹体、冲压发动机、大气环境、测量装置和控制器之间的接口关系,为发动机性能和总体性能评估、技术指标实现提供基础。本发明提出的面向控制的小扰动模型建模方法,为速度控制器设计提供合理的设计模型,可有效提高控制品质。
Description
技术领域
本发明属于控制技术领域,具体涉及一种吸气式巡航导弹的速度控制建模方法。
背景技术
巡航导弹采用冲压发动机为动力,其优点是比冲高且推力连续可调,这使得速度精确控制成为可能,从而满足各种飞行任务要求。
传统导弹的速度控制是通过舵偏调节攻角来实现的,而攻角剖面又受弹道的约束,因此无法实现速度的精确控制。巡航导弹相比传统导弹,其多了一维控制变量,通过燃气流量(或余气系数)来调节推力,同时通过舵偏调节攻角,两者匹配实现速度的精确控制,完成加速爬升、巡航和减速下压攻击等任务。
巡航导弹的推力性能与飞行状态、弹体姿态、飞行环境等密切相关,故其速度控制的约束强,难度大。在总体性能一定的基础上,如何驾驭好冲压发动机以充分发挥其性能,如何设计高效合理的控制方法提高速度控制精度,建立准确的速度控制模型是解决上述问题的关键技术之一,速度控制模型是评估发动机性能、设计控制器、完成地面仿真、实现技术指标的基础与核心。
本发明给出了冲压发动机动态特性的描述方法,根据使用目的的不同提出了面向仿真的巡航导弹的速度控制全量模型建模方法、面向控制的小扰动模型建模方法,为巡航导弹的性能评估和控制器设计奠定坚实基础。
发明内容
本发明提出的巡航导弹速度控制建模方法,给出了速度控制回路的组成及建模流程,具体包括冲压发动机的建模方法、面向仿真的速度回路全量模型、速度回路与高度回路、姿态回路的解耦方法、面向控制的小扰动模型。
一种巡航导弹的速度控制建模方法,包括如下步骤:
(1)建立冲压发动机模型
所述冲压发动机模型包括静态模型和动态模型,静态模型包括推力模型和静压裕度模型,以插值表的形式表示;动态模型包括燃气流量调节特性和燃烧室特性,可通过理论建模和试验验证的方法获得;
(2)建立面向仿真的速度控制模型
所述巡航导弹的速度控制模型由冲压发动机、大气环境、测量装置、控制器和弹体五部分组成,大气环境模型根据飞行状态输出来流参数,从而获得空气流量计算推力和喘振裕度,测量装置模型测量敏感弹体运动给出导弹的运动信息,控制器模型则根据测量装置模型给出的运动信息按照预定速度进行控制,给出燃气流量指令,并输送给固冲发动机,固冲发动机响应燃气流量指令,并与当前的空气流量混合燃烧产生推力,从而实现导弹的速度控制;
(3)建立速度控制的小扰动模型
所述巡航导弹小扰动状态方程为
其中,e1、e2、e4、e5、e6分别表征速度、燃气流量、攻角、舵偏角和弹道倾角对导弹速度的影响系数;ci,i=1,2,3分别表征攻角、弹道倾角、舵偏角改变导弹速度方向的能力;bi,i=1,2,3分别表征弹体角速度、攻角和舵偏角产生的气动力矩改变绕质心运动的能力,ΔX表示某参数X的偏差量,X的含义如下:V为导弹速度,θ为当地弹道倾角,为俯仰角,ωz为弹体角速度,为俯仰舵偏,为燃气流量。
当e4,e5,e6远小于e2时,速度控制器可采用单入单出的控制理论进行单独设计;当e4,e5,e6不可忽略时,则采用多入多出的控制理论进行设计。
所述大气环境模型用以高度为自变量的插值表的形式表示,测量装置用测量精度和二阶传递函数表示,弹体用飞行力学中的运动学和动力学方程表示。
本发明的有益效果如下:本发明给出了冲压发动机动态特性的描述方法,有效考虑了冲压发动机动态特性对速度控制的影响,提高了建模精度。本发明提出的面向仿真的速度控制回路建模方法,有效反映了弹体、冲压发动机、大气环境、测量装置和控制器之间的接口关系,为发动机性能和总体性能评估、、技术指标实现提供基础。本发明提出的面向控制的小扰动模型建模方法,为速度控制器设计提供合理的设计模型,可有效提高控制品质。
附图说明
图1是冲压发动机设计模型的原理图;
图2是本发明解耦后的速度控制回路框图;
图3是传统导弹的速度控制原理图;
图4本发明提出的巡航导弹的速度控制原理图。
具体实施方式
下面结合附图和具体实施方式对本发明的技术方案做进一步详细说明。显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明要求保护的范围。
首先进行冲压发动机的建模,考虑到冲压发动机是一个闭环控制的小系统,其动态特性不能忽略。冲压发动机模型分别用静态模型和动态模型两部分来描述,如图1所示。静态模型是指在不同来流状态和燃气流量(或余气系数)下的稳态推力性能和静压裕度情况,以插值表的形式表示。动态特性描述了从燃气流量(或余气系数)指令到冲压发动机产生推力过程中的延迟、响应时间、超调量等,可通过理论建模和试验验证的方法获得。冲压发动机的动态特性由两部分组成,即燃气流量调节特性和燃烧室特性GP。
根据速度控制原理可将巡航导弹的速度控制回路由冲压发动机、大气环境、测量装置、控制器和弹体五部分组成,其接口关系如图4所示,大气环境模型主要用于建立飞行状态与来流参数之间的关系,从而获得空气流量计算推力和喘振裕度。测量装置敏感弹体运动给出导弹的运动信息,如高度、马赫数、攻角和侧滑角等;控制器则根据测量装置测量到的运动信息按照预定速度进行控制,给出合理的燃气流量指令,并输送给固冲发动机,固冲发动机响应燃气流量指令,并与当前的空气流量混合燃烧产生推力,从而实现导弹的速度控制。该速度控制回路有效的反映了各组成部分的接口关系,根据图4可以快速的建立速度控制模型,便于数学仿真分析。
设导弹相对于纵向平面对称,且以接近于垂直平面弹道飞行(侧向参数很小),则导弹的速度运动与高度、俯仰姿态运动共同构成纵向运动,全量运动方程如下:
其中,为考虑了冲压发动机动态特性的推力;为燃气流量指令;H为高度;V为导弹速度;Fqx1,Fqy1分别为沿弹道x轴和y轴的气动力;m为弹体质量;g为重力加速度;为俯仰角;θ为当地弹道倾角;α为攻角;Mz1为俯仰力矩;Jz1为转动惯量;ω为弹体转动角速度;为俯仰舵偏。
方程(1)与传统导弹的不同在于两方面:一是推力考虑了冲压发动机的动态特性;二是输入变量由一维(舵偏)变为两维(舵偏和燃气流量)。
下面进行面向控制的小扰动模型建模:
采用“小扰动”的假设,在平衡点处对方程(1)进行泰勒展开,去掉高次项得到对应的小扰动方程如下:
其中,d1、d2分别表征速度和弹道倾角改变高度的能力;ei,i=1,2…,6分别表征速度、燃气流量、高度、攻角、舵偏角和弹道倾角对导弹速度的影响系数;ci,i=1,2…,4分别表征攻角、弹道倾角、舵偏角和速度改变导弹速度方向的能力;bi,i=1,2…,4分别表征弹体角速度、攻角、舵偏角、速度产生的气动力矩改变绕质心运动的能力,Δ表示某参数的偏差量,如ΔV为表示导弹速度V的偏差量,为燃气流量的偏差量。
由上式知,高度运动、速度运动与姿态运动耦合,故需分析三个运动的交连情况,分析解耦设计的可能性。由于高度变化和速度变化对应于弹体的长周期,对导弹姿态运动的影响较小,故可忽略高度运动和速度运动对姿态运动的影响,姿态运动回路可独立设计。同理,可忽略高度运动对速度运动的影响。对于采用冲压发动机为动力的导弹速度运动,其阻力和推力均与弹体姿态紧密相关,由式(2)中的e4、e5可知,攻角和舵偏角对速度运动的影响较大,且姿态运动快于速度运动,故速度运动需考虑姿态运动的影响。当姿态运动剧烈时,对速度回路的响应特性要求高,从而提出对冲压发动机的快响应特性,以满足速度精度控制的需求。
式(2)可写为状态方程如式(3)。
其中,e1、e2、e4、e5、e6分别表征速度、燃气流量、攻角、舵偏角和弹道倾角对导弹速度的影响系数;ci,i=1,2,3分别表征攻角、弹道倾角、舵偏角改变导弹速度方向的能力;bi,i=1,2,3分别表征弹体角速度、攻角和舵偏角产生的气动力矩改变绕质心运动的能力,ΔX表示某参数X的偏差量,X的含义如下:V为导弹速度,θ为当地弹道倾角,为俯仰角,ωz为弹体角速度,为俯仰舵偏,为燃气流量。
当e4,e5,e6远小于e2时,则可将姿态运动对速度运动的耦合作为干扰处理,从而速度回路可采用单入单出的经典控制理论进行单独设计,速度回路对应的控制框图如图2所示。
当e4,e5,e6不可忽略时,则需进行速度回路与姿态回路的联合设计,采用多入多出的控制理论进行设计。
依据本发明提供的巡航导弹的速度控制建模方法,可以有效开展速度控制仿真模型和小扰动模型的建立,进行速度控制器设计与仿真,评价总体性能并完成飞行任务。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专利技术人员来说是显而易见的,本文中所定义的一般原理可以在不脱离本发明范围的情况下,在其他实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽范围。
Claims (2)
1.一种巡航导弹的速度控制建模方法,其特征在于,包括如下步骤:
(1)建立冲压发动机模型
所述冲压发动机模型包括静态模型和动态模型,静态模型包括推力模型和静压裕度模型,以插值表的形式表示;动态模型包括燃气流量调节特性和燃烧室特性,可通过理论建模和试验验证的方法获得;
(2)建立面向仿真的速度控制模型
所述巡航导弹的速度控制模型由冲压发动机、大气环境、测量装置、控制器和弹体五部分组成,大气环境模型根据飞行状态输出来流参数,从而获得空气流量计算推力和喘振裕度,测量装置模型测量敏感弹体运动给出导弹的运动信息,控制器模型则根据测量装置模型给出的运动信息按照预定速度进行控制,给出燃气流量指令,并输送给固冲发动机,固冲发动机响应燃气流量指令,并与当前的空气流量混合燃烧产生推力,从而实现导弹的速度控制;
(3)建立速度控制的小扰动模型
所述巡航导弹小扰动状态方程为
其中,e1、e2、e4、e5、e6分别表征速度、燃气流量、攻角、舵偏角和弹道倾角对导弹速度的影响系数;ci,i=1,2,3分别表征攻角、弹道倾角、舵偏角改变导弹速度方向的能力;bi,i=1,2,3分别表征弹体角速度、攻角和舵偏角产生的气动力矩改变绕质心运动的能力,ΔX表示某参数X的偏差量,X的含义如下:V为导弹速度,θ为当地弹道倾角,为俯仰角,ωz为弹体角速度,为俯仰舵偏,为燃气流量,
当e4,e5,e6远小于e2时,速度控制器可采用单入单出的控制理论进行单独设计;当e4,e5,e6不可忽略时,则采用多入多出的控制理论进行设计。
2.根据权利要求1所述的一种巡航导弹的速度控制建模方法,其特征在于,所述大气环境模型用以高度为自变量的插值表的形式表示,测量装置用测量精度和二阶传递函数表示,弹体用飞行力学中的运动学和动力学方程表示。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610157141.3A CN106125571B (zh) | 2016-03-18 | 2016-03-18 | 一种巡航导弹的速度控制建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610157141.3A CN106125571B (zh) | 2016-03-18 | 2016-03-18 | 一种巡航导弹的速度控制建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106125571A CN106125571A (zh) | 2016-11-16 |
CN106125571B true CN106125571B (zh) | 2019-05-07 |
Family
ID=57270349
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610157141.3A Active CN106125571B (zh) | 2016-03-18 | 2016-03-18 | 一种巡航导弹的速度控制建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106125571B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106950982B (zh) * | 2017-02-16 | 2020-02-14 | 北京临近空间飞行器系统工程研究所 | 再入飞行器姿控动力系统高空力矩特性辨识方法 |
CN109085847B (zh) * | 2018-09-03 | 2021-12-07 | 中国运载火箭技术研究院 | 一种大静不稳定度飞行器弹性稳定控制系统及控制方法 |
CN110765404B (zh) * | 2019-06-14 | 2023-08-25 | 上海航天控制技术研究所 | 一种机弹干扰辨识算法 |
CN111046497B (zh) * | 2019-12-24 | 2023-04-07 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种高空高速飞机突防生存力快速评估装置 |
CN112149292B (zh) * | 2020-09-13 | 2022-05-31 | 中国运载火箭技术研究院 | 发动机干扰力矩确定方法、设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
USH159H (en) * | 1984-12-17 | 1986-11-04 | The United States Of America As Represented By The Secretary Of The Navy | Shaped trajectory cruise missile launch mode |
CN103587718A (zh) * | 2013-11-07 | 2014-02-19 | 北京临近空间飞行器系统工程研究所 | 一种高超声速飞行器姿态稳定的频率管理方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7800032B1 (en) * | 2006-11-30 | 2010-09-21 | Raytheon Company | Detachable aerodynamic missile stabilizing system |
-
2016
- 2016-03-18 CN CN201610157141.3A patent/CN106125571B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
USH159H (en) * | 1984-12-17 | 1986-11-04 | The United States Of America As Represented By The Secretary Of The Navy | Shaped trajectory cruise missile launch mode |
CN103587718A (zh) * | 2013-11-07 | 2014-02-19 | 北京临近空间飞行器系统工程研究所 | 一种高超声速飞行器姿态稳定的频率管理方法 |
Non-Patent Citations (3)
Title |
---|
冲压发动机控制路径研究;陈新民等;《航天控制》;20130630;第31卷(第3期);第11-14,96页 |
基于标准遗传算法的巡航导弹造型_气动_隐身一体化设计;王旭刚等;《中国航空学会2007年学术年会》;20071231;第1-8页 |
高超声速飞行器动力系统研究进展;王振国等;《力学进展》;20091125;第39卷(第6期);第716-735页 |
Also Published As
Publication number | Publication date |
---|---|
CN106125571A (zh) | 2016-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106125571B (zh) | 一种巡航导弹的速度控制建模方法 | |
Garza et al. | A collection of nonlinear aircraft simulations in matlab | |
CN107942651B (zh) | 一种近空间飞行器控制系统 | |
CN109426146A (zh) | 高超声速飞行器的高阶非奇异Terminal滑模控制方法 | |
CN102073755A (zh) | 近空间高超声速飞行器运动控制仿真方法 | |
CN109709978B (zh) | 一种高超声速飞行器制导控制一体化设计方法 | |
Ohanian | Ducted fan aerodynamics and modeling, with applications of steady and synthetic jet flow control | |
CN110568765A (zh) | 面向攻角跟踪的高超声速飞行器非对称输出受限控制方法 | |
CN113868758B (zh) | 一种基于计算力学的弹性高超声速飞行器建模方法 | |
CN114564045A (zh) | 一种考虑降雨和阵风条件的无人机飞行控制律设计方法 | |
CN106896722A (zh) | 采用状态反馈与神经网络的高超飞行器复合控制方法 | |
CN111506113A (zh) | 飞行器制导指令计算方法、侧滑角计算方法及制导方法 | |
CN107870563A (zh) | 一种旋转弹全阶反馈控制器的插值增益调度方法 | |
CN109857146A (zh) | 一种基于前馈和权值分配的分层无人机跟踪控制方法 | |
Li et al. | Flight-propulsion integration dynamic analysis and adaptive control of the hypersonic vehicle at wide-range Mach numbers | |
Raju et al. | Empirical virtual sliding target guidance law design: An aerodynamic approach | |
Dou et al. | Analysis of the impact of translating cowl on hypersonic vehicle performance | |
Sridharan et al. | Impact of control specifications on vehicle design for scramjet-powered hypersonic vehicles | |
CN107489531B (zh) | 基于半积分与增益自适应的高超声速发动机供油规律设计方法 | |
Dong et al. | Maneuvering Strategy and Own Aircraft Movement Prediction in Trial Input Method----Low Angle of Attack | |
CN105674804B (zh) | 一种包含法向加速度导数的空射巡航弹下滑段多约束制导方法 | |
Amoignon et al. | Impact of aerodynamic modeling on control system design for a military aircraft concept | |
Wang et al. | Virtual flight simulation of the basic finner projectile based on fuzz control | |
Chen et al. | Modeling and Analysis of Fluid-Thermal-Structure Coupling Problems for Hypersonic Vehicles | |
Shabani et al. | Cost Function Determination for a WIG in Predefined Path and Height Using Conjugate Gradient Method |
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 |