CN109085847A - 一种大静不稳定度飞行器弹性稳定控制系统及控制方法 - Google Patents
一种大静不稳定度飞行器弹性稳定控制系统及控制方法 Download PDFInfo
- Publication number
- CN109085847A CN109085847A CN201811021187.8A CN201811021187A CN109085847A CN 109085847 A CN109085847 A CN 109085847A CN 201811021187 A CN201811021187 A CN 201811021187A CN 109085847 A CN109085847 A CN 109085847A
- Authority
- CN
- China
- Prior art keywords
- aircraft
- aerodynamic moment
- trajectory
- coefficient
- elastic
- 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
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000013461 design Methods 0.000 claims abstract description 71
- 230000006641 stabilisation Effects 0.000 claims abstract description 39
- 238000011105 stabilization Methods 0.000 claims abstract description 39
- 230000000087 stabilizing effect Effects 0.000 claims abstract description 36
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 claims abstract description 18
- 230000001133 acceleration Effects 0.000 claims description 10
- 230000001186 cumulative effect Effects 0.000 claims description 10
- 230000003068 static effect Effects 0.000 claims description 8
- 238000009434 installation Methods 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 abstract description 10
- 238000007796 conventional method Methods 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000003786 synthesis reaction Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000033001 locomotion Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/08—Control of attitude, i.e. control of roll, pitch, or yaw
- G05D1/0808—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明涉及一种大静不稳定度飞行器弹性稳定控制系统及控制方法,通过建立飞行器控制系统设计模型,选取气动力矩系数优化变量;通过弹道与姿态联合优化,降低气动力矩系数;并根据气动力矩系数设定值进行幅值相位稳定网络切换系数的选择,选取气动力矩系数较小区域,进行弹性幅值稳定网络设计,气动力矩系数较大区域,进行弹性相位稳定网络设计,实现飞行器幅值稳定和相位稳定的控制切换,保证飞行器的全频段频域稳定,该方法解决了大静不稳定度细长飞行器的弹性稳定控制难题,可广泛应用于飞行器姿态控制系统设计中。
Description
技术领域
本发明涉及一种大静不稳定度细长飞行器弹性稳定控制系统及控制方法,属于飞行器姿态控制系统设计领域,可广泛应用于细长飞行器的弹性稳定设计中。
背景技术
大静不稳定度细长飞行器,由于静不稳定度较大,要求刚体控制截频较高,以保证刚体稳定控制,而细长体飞行器的弹体弹性频率较低,从而导致飞行器的刚体控制截频与弹体一阶弹性频率接近,刚体稳定与弹性稳定矛盾突出。如采用传统飞行器频域设计方法,进行飞行器弹性幅值稳定,则由于弹性滤波对刚体相位影响过大,将导致刚体相位裕度过低甚至刚体不稳定。因此,需要开展大静不稳定度细长飞行器弹性稳定控制方法的研究,寻求能实现大静不稳定度细长飞行器稳定控制的方法。
发明内容
本发明的目的在于克服现有技术的上述缺陷,提供一种大静不稳定度飞行器弹性稳定控制方法,该方法以解决大静不稳定度细长飞行器弹性稳定控制难题为目的,开展弹道与姿控联合优化设计,并进行姿态控制方法优化设计,解决了大静不稳定度细长飞行器的弹性稳定控制难题,该方法可广泛应用于飞行器姿态控制系统设计中。
本发明的另外一个目的在于提供一种大静不稳定度飞行器弹性稳定控制系统。
本发明的上述目的主要是通过如下技术方案予以实现的:
一种大静不稳定度飞行器弹性稳定控制方法,包括:
建立飞行器控制系统设计模型,得到气动力矩系数模型,将气动力矩系数作为弹道优化参数;
调整弹道程序攻角,根据所述气动力矩系数模型计算气动力矩系数,使气动力矩系数最大值达到第一气动力矩系数设定值,并取气动力矩系数最大值达到第一气动力矩系数设定值的弹道作为标准弹道;
根据所述标准弹道建立包含幅值相位稳定网络切换系数的控制方程,并根据第二气动力矩系数设定值进行所述幅值相位稳定网络切换系数的选择,实现飞行器幅值稳定和相位稳定的控制切换。
在上述大静不稳定度飞行器弹性稳定控制方法中,所述建立的飞行器控制系统设计模型如下:
其中:Δθ为弹道倾角偏差,为弹道倾角偏差变化率,Δα为攻角偏差,为俯仰等效舵偏角,αw为风攻角,为合成干扰力系数,为俯仰角偏差,为合成干扰力矩系数,为弹性振动第i次振型广义坐标,为弹性振动第i次振型广义速度;为弹性振动第i次振型广义加速度;ωfi为弹性振动第i次振型固有频率,ξi为弹性振动第i次振型阻尼系数,为惯组敏感的俯仰角偏差,为俯仰偏差角速度;为俯仰偏差角加速度;Wi'(Xg)为惯组安装处第i次振型斜率;c1f为升力系数;c2f为引力系数;c3f为控制力系数;c′1f为风干扰对应的升力系数。
在上述大静不稳定度飞行器弹性稳定控制方法中,根据建立的飞行器控制系统设计模型得到气动力矩系数模型如下:
其中:b2f为气动力矩系数;Cmz为俯仰力矩系数,lk为参考长度,α为弹道程序攻角,q为动压,Sm为参考面积,Jz1为转动惯量。
在上述大静不稳定度飞行器弹性稳定控制方法中,所述第一气动力矩系数设定值-b2fmax1取值为25~35。
在上述大静不稳定度飞行器弹性稳定控制方法中,根据所述标准弹道建立含幅值相位稳定网络切换系数的控制方程如下:
其中:为俯仰舵偏角,为俯仰通道静态增益,Kfx为幅值相位稳定网络切换系数,为俯仰通道幅值稳定网络,为俯仰通道相位稳定网络,为俯仰角偏差。
在上述大静不稳定度飞行器弹性稳定控制方法中,所述含幅值相位稳定网络切换系数Kfx取值为1或0。
在上述大静不稳定度飞行器弹性稳定控制方法中,所述第二气动力矩系数设定值-b2fmax2取值为10~20。
在上述大静不稳定度飞行器弹性稳定控制方法中,根据第二气动力矩系数设定值进行所述幅值相位稳定网络切换系数的选择的具体方法为:
若-b2f<-b2fmax2,则Kfx=1,
若-b2f≥-b2fmax2,则Kfx=0。
在上述大静不稳定度飞行器弹性稳定控制方法中,
当Kfx=1时,实现飞行器幅值稳定控制,所述飞行器幅值稳定的控制方程为:
当Kfx=0时,实现飞行器相位稳定控制,所述飞行器相位稳定的控制方程为:
在上述大静不稳定度飞行器弹性稳定控制方法中,所述俯仰通道幅值稳定网络与俯仰通道相位稳定网络的具体表示形式如下:
其中:ω1,ω2,ω3,ω4,ω5,ω6,ω7,ω8,ω9,ω10,ω11,ω12为转折频率;
ξ3,ξ4,ξ5,ξ6,ξ7,ξ8,ξ9,ξ10,ξ11,ξ12为阻尼系数;s为拉式算子。
一种大静不稳定度飞行器弹性稳定控制系统,包括:设计模型建立模块、弹道优化参数选择模块、弹道优化模块和控制参数设计模块,其中:
设计模型建立模块,建立飞行器控制系统设计模型,根据所述飞行器控制系统设计模型得到气动力矩系数模型,并将所述气动力矩系数模型发送给弹道优化参数选择模块;
弹道优化参数选择模块,接收设计模型建立模块发送的气动力矩系数模型,选择气动力矩系数作为弹道优化参数,将所述气动力矩系数模型发送给弹道优化模块;
弹道优化模块,接收弹道优化参数选择模块发送的气动力矩系数模型,调整弹道程序攻角,根据所述气动力矩系数模型计算气动力矩系数,使气动力矩系数最大值达到第一气动力矩系数设定值,并取气动力矩系数最大值达到第一气动力矩系数设定值的弹道作为标准弹道,并将所述标准弹道发送给控制参数设计模块;
控制参数设计模块,接收弹道优化模块发送的标准弹道,根据所述标准弹道建立含幅值相位稳定网络切换系数的控制方程,并根据第二气动力矩系数设定值进行所述幅值相位稳定网络切换系数的选择,实现飞行器幅值稳定和相位稳定的控制切换。
在上述大静不稳定度飞行器弹性稳定控制系统中,所述设计模型建立模块建立的飞行器控制系统设计模型如下:
其中:Δθ为弹道倾角偏差,为弹道倾角偏差变化率,Δα为攻角偏差,为俯仰等效舵偏角,αw为风攻角,为合成干扰力系数,为俯仰角偏差,为合成干扰力矩系数,为弹性振动第i次振型广义坐标,为弹性振动第i次振型广义速度;为弹性振动第i次振型广义加速度;ωfi为弹性振动第i次振型固有频率,ξi为弹性振动第i次振型阻尼系数,为惯组敏感的俯仰角偏差,为俯仰偏差角速度;为俯仰偏差角加速度;Wi'(Xg)为惯组安装处第i次振型斜率;c1f为升力系数;c2f为引力系数;c3f为控制力系数;c′1f为风干扰对应的升力系数。
在上述大静不稳定度飞行器弹性稳定控制系统中,所述设计模型建立模块根据飞行器控制系统设计模型得到的气动力矩系数模型如下:
其中:b2f为气动力矩系数;Cmz为俯仰力矩系数,lk为参考长度,α为弹道程序攻角,q为动压,Sm为参考面积,Jz1为转动惯量。
在上述大静不稳定度飞行器弹性稳定控制系统中,所述控制参数设计模块根据标准弹道建立含幅值相位稳定网络切换系数的控制方程如下:
其中:为俯仰舵偏角,为俯仰通道静态增益,Kfx为幅值相位稳定网络切换系数,为俯仰通道幅值稳定网络,为俯仰通道相位稳定网络,为俯仰角偏差。
在上述大静不稳定度飞行器弹性稳定控制系统中,所述第一气动力矩系数设定值-b2fmax1取值为25~35;所述第二气动力矩系数设定值-b2fmax2取值为10~20。
在上述大静不稳定度飞行器弹性稳定控制系统中,所述含幅值相位稳定网络切换系数Kfx取值为1或0。
在上述大静不稳定度飞行器弹性稳定控制系统中,所述控制参数设计模块根据第二气动力矩系数设定值进行幅值相位稳定网络切换系数的选择的具体方法为:
若b2f<-b2fmax2,则Kfx=1,
若b2f≥-b2fmax2,则Kfx=0。
本发明与现有技术相比的有益效果如下:
(1)、本发明通过采用基于稳定裕度需求的弹道与姿态联合优化方法,通过理论分析及仿真试验,选取气动力矩系数作为弹道优化参数,选取的气动力矩系数对姿态控制系统设计稳定裕度影响最大,通过弹道与姿控的迭代优化设计,有效减小气动力矩系数,降低姿态控制系统设计难度,该设计方法可推广应用于各类飞行器的姿态控制;
(2)、本发明针对大静不稳定度细长飞行器,刚体控制截频与弹体弹性频率接近,姿态控制系统频域设计困难,通过采用大静不稳定度细长飞行器幅相稳定结合控制方法,合理设计包含幅值相位稳定网络切换系数的控制方程,控制方程中包含了幅值稳定和相位稳定控制网络,能够实现大静不稳定度细长飞行器的弹性稳定控制;
(3)、本发明根据标准弹道建立含幅值相位稳定网络切换系数的控制方程,并根据气动力矩系数设定值进行幅值相位稳定网络切换系数的选择,可以实现飞行器幅值稳定和相位稳定的控制切换,本发明采用幅相稳定相结合的方式,在静不稳定度较大区域采用相位稳定方案,解决了刚体与弹性稳定矛盾,在静不稳定度较小区域,采用幅值稳定,显著提高弹性稳定的可靠性;
(4)、本发明方法可以推广应用于各类大静不稳定度细长飞行器的弹性稳定控制,提高姿态控制系统适应性;
(5)、本发明通过理论分析结合仿真试验,合理确定气动力矩系数最大值的优化设定值,以及幅相位稳定切换系数对应的气动力矩系数优化设定值,使得幅值稳定与相位稳定均具有足够的稳定裕度。
(6)、本发明在细长飞行器的研制中进行了成功应用,取得了较好的效果,证明大静不稳定度细长飞行器弹性稳定控制方法的有效性,可以解决大静不稳定度细长飞行器姿态控制系统设计难题。
附图说明
图1为本发明弹道程序攻角与传统方法对比的曲线图;
图2为本发明弹道动压与传统方法对比的曲线图;
图3为本发明飞行器弹性稳定控制方法流程图;
图4为本发明飞行器弹性稳定控制控制系统结构示意图。
图5为本发明实施例中飞行器相位稳定频域综合图;
图6为本发明实施例中飞行器幅值稳定频域综合图。
具体实时方式
下面结合附图和具体实施例对本发明作进一步详细的描述:
如图3所示为本发明飞行器弹性稳定控制方法流程图,本发明大静不稳定度飞行器弹性稳定控制方法具体包括如下步骤:
一、建立飞行器控制系统设计模型,得到气动力矩系数模型,将气动力矩系数作为弹道优化参数。
建立飞行器俯仰通道的刚体运动及弹性运动数学模型如下:
其中:Δθ为弹道倾角偏差,为弹道倾角偏差变化率,Δα为攻角偏差,为俯仰等效舵偏角,αw为风攻角,为合成干扰力系数,为俯仰角偏差,为合成干扰力矩系数,为弹性振动第i次振型广义坐标,为弹性振动第i次振型广义速度;为弹性振动第i次振型广义加速度;ωfi为弹性振动第i次振型固有频率,ξi为弹性振动第i次振型阻尼系数,为惯组敏感的俯仰角偏差,为俯仰偏差角速度;为俯仰偏差角加速度;Wi'(Xg)为惯组安装处第i次振型斜率;c1f为升力系数;c2f为引力系数;c3f为控制力系数;c′1f为风干扰对应的升力系数。
大静不稳定度细长飞行器弹性稳定控制难点主要在于刚体控制截频与弹体一阶弹性频率较近,在进行弹性滤波网络设计时,将对刚体相位裕度产生较大影响,甚至造成刚体运动的不稳定。而直接影响刚体控制截频的为气动力矩系数b2f,气动力矩系数模型如下:
其中:Cmz为俯仰力矩系数,lk为参考长度,α为攻角,q为动压,Sm为参考面积,Jz1为转动惯量。
可见,影响-b2f大小的主要为全弹的静不稳定度、飞行动压及转动惯量,由于飞行器规模确定后,转动惯量很难改变。因此,可以通过弹道与姿控联合优化,设计弹道程序攻角,改变飞行器的静不稳定度和飞行动压,降低-b2f系数,从而降低姿态控制系统设计难度。
二、采用基于稳定裕度需求的弹道与姿态联合优化方法,调整弹道程序攻角,根据气动力矩系数模型计算气动力矩系数,使气动力矩系数最大值达到第一气动力矩系数设定值,并取气动力矩系数最大值达到第一气动力矩系数设定值的弹道作为标准弹道。
开展弹道与姿控的联合优化,针对飞行攻角进行优化设计,采用全程下压攻角,降低最大动压,并进行姿控与弹道的迭代优化设计,即调整弹道设计攻角α,姿态控制系统进行六自由度数学仿真,统计最大-b2f系数,并根据频域设计要求和控制力需求提出-b2f系数改进意见,即通过降低或提高程序攻角α,从而减小-b2f系数,使得最大-b2f系数达到第一气动力矩系数设定值-b2fmax,本发明中第一气动力矩系数设定值-b2fmax1取值为25~35,具体地,本发明实施例选取-b2fmax=30。
以某飞行器为例,弹道优化设计结果如图1和图2所示,图1为本发明弹道程序攻角与传统方法对比的曲线图;图2为本发明弹道动压与传统方法对比的曲线图;其中,实线为采用传统飞行器攻角一次下压方式,虚线为采用本发明方法设计的攻角。
由图1可知,本发明采用了常值攻角下压弹道,传统方法采用一次攻角下压弹道。由图2可知,采用本发明设计的弹道动压相比传统方法设计的弹道动压有较大幅度降低。
三、采用大静不稳定度细长飞行器幅相稳定结合控制方法,根据标准弹道建立包含幅值相位稳定网络切换系数的控制方程,并根据第二气动力矩系数设定值进行幅值相位稳定网络切换系数的选择,实现飞行器幅值稳定和相位稳定的控制切换。
针对弹性稳定控制,工程上有两种方法,幅值稳定及相位稳定。幅值稳定的实质为执行机构对弹性振动的激励小于其在固有阻尼下产生的衰减,而相位稳定的实质是把弹性振动信号作为控制信号的一部分,通过控制装置得到合适的相位,进行弹性运动的主动抑制。
针对大静不稳定度细长飞行器,刚体控制截频提高,与飞行器一阶弹性频率接近,进行弹性滤波网络设计时对刚体相位裕度影响较大,甚至造成刚体不稳定,刚体稳定与弹性稳定无法兼顾的难题,本发明进行姿控网络的分区设计和在线切换,在-b2f较大区域,对进行一阶弹性的相位稳定,减小对刚体的相位滞后,解决刚体与弹性频域稳定的矛盾;在-b2f较小区域,切换为幅值稳定,提高弹性稳定可靠性。
本发明根据标准弹道建立含幅值相位稳定网络切换系数的控制方程如下:
其中:为俯仰舵偏角,为俯仰通道静态增益,Kfx为幅值相位稳定网络切换系数,为俯仰通道幅值稳定网络,为俯仰通道相位稳定网络,为俯仰角偏差。
绘制弹体频域特性曲线,进行姿态控制系统的初步设计,综合评估设计难度及稳定裕度后选取Kfx系数,Kfx以马赫数进行分区,当飞行器飞行实时马赫数小于Maqh,Kfx=1,当飞行器飞行实时马赫数大于Maqh,Kfx=0。本发明中含幅值相位稳定网络切换系数Kfx取值为1或0。第二气动力矩系数设定值-b2fmax2取值为10~20。
本发明根据第二气动力矩系数设定值进行幅值相位稳定网络切换系数的选择的具体方法为:
若-b2f<-b2fmax2,则Kfx=1,
若-b2f≥-b2fmax2,则Kfx=0。
即当气动力矩系数-b2f小于第二气动力矩系数设定值-b2fmax2时,含幅值相位稳定网络切换系数Kfx取值为1,当气动力矩系数-b2f大于或等于第二气动力矩系数设定值-b2fmax2时,含幅值相位稳定网络切换系数Kfx取值为0。
具体地,本发明实施例中选取了-b2f大于或等于15(对应马赫数1.5)进行相位稳定设计,-b2f小于15进行幅值稳定设计,即当马赫数小于1.5时,Kfx=1,马赫数大于或等于1.5时,Kfx=0。
当Kfx=1时,实现飞行器幅值稳定控制,飞行器幅值稳定的控制方程为;
当Kfx=0时,实现飞行器相位稳定控制,飞行器相位稳定的控制方程为:
姿控网络及的连续网络形式,可以采用:
其中:ω1,ω2,ω3,ω4,ω5,ω6,ω7,ω8,ω9,ω10,ω11,ω12为转折频率;
ξ3,ξ4,ξ5,ξ6,ξ7,ξ8,ξ9,ξ10,ξ11,ξ12为阻尼系数;s为拉式算子。
俯仰通道静态增益姿控网络及参数通过如下方式确定:根据飞行器传递函数,采用MATLAB绘制伯德图,设计俯仰通道静态增益姿控网络及各参数,满足稳定裕度要求,保证系统频域稳定。
可见,针对大静不稳定区域进行相位稳定,小静不稳定区域采用幅值稳定的姿控网络分区设计方法,能够实现全飞行段姿态稳定控制,并具有一定的稳定裕度,从而解决大静不稳定度细长飞行器的弹性稳定控制难题。
如图5所示为本发明实施例中飞行器相位稳定频域综合图,图6为本发明实施例中飞行器幅值稳定频域综合图。由图5、图6可知,采用本发明方法设计的全弹道控制参数能够满足刚体幅值裕度大于6dB,相位裕度大于30°,弹性幅值裕度大于6dB,弹性相位裕度大于30°的稳定裕度要求,保证系统频域稳定。
如图4所示为本发明飞行器弹性稳定控制控制系统结构示意图,由图可知本发明大静不稳定度飞行器弹性稳定控制系统,包括设计模型建立模块、弹道优化参数选择模块、弹道优化模块和控制参数设计模块。
设计模型建立模块,建立飞行器控制系统设计模型,根据所述飞行器控制系统设计模型得到气动力矩系数模型,并将所述气动力矩系数模型发送给弹道优化参数选择模块。
弹道优化参数选择模块,接收设计模型建立模块发送的气动力矩系数模型,选择气动力矩系数作为弹道优化参数,将所述气动力矩系数模型发送给弹道优化模块。
弹道优化模块,接收弹道优化参数选择模块发送的气动力矩系数模型,调整弹道程序攻角,根据所述气动力矩系数模型计算气动力矩系数,使气动力矩系数最大值达到第一气动力矩系数设定值,并取气动力矩系数最大值达到第一气动力矩系数设定值的弹道作为标准弹道,并将所述标准弹道发送给控制参数设计模块。
控制参数设计模块,接收弹道优化模块发送的标准弹道,根据所述标准弹道建立含幅值相位稳定网络切换系数的控制方程,并根据第二气动力矩系数设定值进行所述幅值相位稳定网络切换系数的选择,实现飞行器幅值稳定和相位稳定的控制切换。
具体地,本发明实施例中设计模型建立模块根据飞行器控制系统设计模型得到的气动力矩系数模型如下:
其中:b2f为气动力矩系数;Cmz为俯仰力矩系数,lk为参考长度,α为弹道程序攻角,q为动压,Sm为参考面积,Jz1为转动惯量。
具体地,本发明实施例中,控制参数设计模块根据标准弹道建立含幅值相位稳定网络切换系数的控制方程如下:
其中:为俯仰舵偏角,为俯仰通道静态增益,Kfx为幅值相位稳定网络切换系数,为俯仰通道幅值稳定网络,为俯仰通道相位稳定网络,为俯仰角偏差。
控制参数设计模块根据第二气动力矩系数设定值进行幅值相位稳定网络切换系数的选择的具体方法为:
若b2f<-b2fmax2,则Kfx=1,
若b2f≥-b2fmax2,则Kfx=0;
当Kfx=1时,实现飞行器幅值稳定控制,所述飞行器幅值稳定的控制方程为;
当Kfx=0时,实现飞行器相位稳定控制,所述飞行器相位稳定的控制方程为:
具体地,本发明实施例中,第一气动力矩系数设定值-b2fmax1取值为25~35;所述第二气动力矩系数设定值-b2fmax2取值为10~20。
以上所述,仅为本发明一个具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
Claims (17)
1.一种大静不稳定度飞行器弹性稳定控制方法,其特征在于:包括:
建立飞行器控制系统设计模型,得到气动力矩系数模型,将气动力矩系数作为弹道优化参数;
调整弹道程序攻角,根据所述气动力矩系数模型计算气动力矩系数,使气动力矩系数最大值达到第一气动力矩系数设定值,并取气动力矩系数最大值达到第一气动力矩系数设定值的弹道作为标准弹道;
根据所述标准弹道建立包含幅值相位稳定网络切换系数的控制方程,并根据第二气动力矩系数设定值进行所述幅值相位稳定网络切换系数的选择,实现飞行器幅值稳定和相位稳定的控制切换。
2.根据权利要求1所述的大静不稳定度飞行器弹性稳定控制方法,其特征在于:所述建立的飞行器控制系统设计模型如下:
其中:Δθ为弹道倾角偏差,为弹道倾角偏差变化率,Δα为攻角偏差,为俯仰等效舵偏角,αw为风攻角,为合成干扰力系数,为俯仰角偏差,为合成干扰力矩系数,为弹性振动第i次振型广义坐标,为弹性振动第i次振型广义速度;为弹性振动第i次振型广义加速度;ωfi为弹性振动第i次振型固有频率,ξi为弹性振动第i次振型阻尼系数,为惯组敏感的俯仰角偏差,为俯仰偏差角速度;为俯仰偏差角加速度;Wi'(Xg)为惯组安装处第i次振型斜率;c1f为升力系数;c2f为引力系数;c3f为控制力系数;c′1f为风干扰对应的升力系数。
3.根据权利要求2所述的大静不稳定度飞行器弹性稳定控制方法,其特征在于:根据建立的飞行器控制系统设计模型得到气动力矩系数模型如下:
其中:b2f为气动力矩系数;Cmz为俯仰力矩系数,lk为参考长度,α为弹道程序攻角,q为动压,Sm为参考面积,Jz1为转动惯量。
4.根据权利要求1所述的大静不稳定度飞行器弹性稳定控制方法,其特征在于:所述第一气动力矩系数设定值-b2fmax1取值为25~35。
5.根据权利要求1所述的大静不稳定度飞行器弹性稳定控制方法,其特征在于:根据所述标准弹道建立含幅值相位稳定网络切换系数的控制方程如下:
其中:为俯仰舵偏角,为俯仰通道静态增益,Kfx为幅值相位稳定网络切换系数,为俯仰通道幅值稳定网络,为俯仰通道相位稳定网络,为俯仰角偏差。
6.根据权利要求1所述的大静不稳定度飞行器弹性稳定控制方法,其特征在于:所述含幅值相位稳定网络切换系数Kfx取值为1或0。
7.根据权利要求6所述的大静不稳定度飞行器弹性稳定控制方法,其特征在于:所述第二气动力矩系数设定值-b2fmax2取值为10~20。
8.根据权利要求1~7之一所述的大静不稳定度飞行器弹性稳定控制方法,其特征在于:根据第二气动力矩系数设定值进行所述幅值相位稳定网络切换系数的选择的具体方法为:
若-b2f<-b2fmax2,则Kfx=1,
若-b2f≥-b2fmax2,则Kfx=0。
9.根据权利要求8所述的大静不稳定度飞行器弹性稳定控制方法,其特征在于:
当Kfx=1时,实现飞行器幅值稳定控制,所述飞行器幅值稳定的控制方程为:
当Kfx=0时,实现飞行器相位稳定控制,所述飞行器相位稳定的控制方程为:
10.根据权利要求5~7之一所述的大静不稳定度飞行器弹性稳定控制方法,其特征在于:所述俯仰通道幅值稳定网络与俯仰通道相位稳定网络的具体表示形式如下:
其中:ω1,ω2,ω3,ω4,ω5,ω6,ω7,ω8,ω9,ω10,ω11,ω12为转折频率;
ξ3,ξ4,ξ5,ξ6,ξ7,ξ8,ξ9,ξ10,ξ11,ξ12为阻尼系数;s为拉式算子。
11.一种大静不稳定度飞行器弹性稳定控制系统,其特征在于:包括:设计模型建立模块、弹道优化参数选择模块、弹道优化模块和控制参数设计模块,其中:
设计模型建立模块,建立飞行器控制系统设计模型,根据所述飞行器控制系统设计模型得到气动力矩系数模型,并将所述气动力矩系数模型发送给弹道优化参数选择模块;
弹道优化参数选择模块,接收设计模型建立模块发送的气动力矩系数模型,选择气动力矩系数作为弹道优化参数,将所述气动力矩系数模型发送给弹道优化模块;
弹道优化模块,接收弹道优化参数选择模块发送的气动力矩系数模型,调整弹道程序攻角,根据所述气动力矩系数模型计算气动力矩系数,使气动力矩系数最大值达到第一气动力矩系数设定值,并取气动力矩系数最大值达到第一气动力矩系数设定值的弹道作为标准弹道,并将所述标准弹道发送给控制参数设计模块;
控制参数设计模块,接收弹道优化模块发送的标准弹道,根据所述标准弹道建立含幅值相位稳定网络切换系数的控制方程,并根据第二气动力矩系数设定值进行所述幅值相位稳定网络切换系数的选择,实现飞行器幅值稳定和相位稳定的控制切换。
12.根据权利要求11所述的大静不稳定度飞行器弹性稳定控制系统,其特征在于:所述设计模型建立模块建立的飞行器控制系统设计模型如下:
其中:Δθ为弹道倾角偏差,为弹道倾角偏差变化率,Δα为攻角偏差,为俯仰等效舵偏角,αw为风攻角,为合成干扰力系数,为俯仰角偏差,为合成干扰力矩系数,为弹性振动第i次振型广义坐标,为弹性振动第i次振型广义速度;为弹性振动第i次振型广义加速度;ωfi为弹性振动第i次振型固有频率,ξi为弹性振动第i次振型阻尼系数,为惯组敏感的俯仰角偏差,为俯仰偏差角速度;为俯仰偏差角加速度;Wi'(Xg)为惯组安装处第i次振型斜率;c1f为升力系数;c2f为引力系数;c3f为控制力系数;c′1f为风干扰对应的升力系数。
13.根据权利要求11所述的大静不稳定度飞行器弹性稳定控制系统,其特征在于:所述设计模型建立模块根据飞行器控制系统设计模型得到的气动力矩系数模型如下:
其中:b2f为气动力矩系数;Cmz为俯仰力矩系数,lk为参考长度,α为弹道程序攻角,q为动压,Sm为参考面积,Jz1为转动惯量。
14.根据权利要求11所述的大静不稳定度飞行器弹性稳定控制系统,其特征在于:所述控制参数设计模块根据标准弹道建立含幅值相位稳定网络切换系数的控制方程如下:
其中:为俯仰舵偏角,为俯仰通道静态增益,Kfx为幅值相位稳定网络切换系数,为俯仰通道幅值稳定网络,为俯仰通道相位稳定网络,为俯仰角偏差。
15.根据权利要求11所述的大静不稳定度飞行器弹性稳定控制系统,其特征在于:所述第一气动力矩系数设定值-b2fmax1取值为25~35;所述第二气动力矩系数设定值-b2fmax2取值为10~20。
16.根据权利要求11所述的大静不稳定度飞行器弹性稳定控制系统,其特征在于:所述含幅值相位稳定网络切换系数Kfx取值为1或0。
17.根据权利要求11~16之一所述的大静不稳定度飞行器弹性稳定控制系统,其特征在于:所述控制参数设计模块根据第二气动力矩系数设定值进行幅值相位稳定网络切换系数的选择的具体方法为:
若b2f<-b2fmax2,则Kfx=1,
若b2f≥-b2fmax2,则Kfx=0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811021187.8A CN109085847B (zh) | 2018-09-03 | 2018-09-03 | 一种大静不稳定度飞行器弹性稳定控制系统及控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811021187.8A CN109085847B (zh) | 2018-09-03 | 2018-09-03 | 一种大静不稳定度飞行器弹性稳定控制系统及控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109085847A true CN109085847A (zh) | 2018-12-25 |
CN109085847B CN109085847B (zh) | 2021-12-07 |
Family
ID=64840636
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811021187.8A Active CN109085847B (zh) | 2018-09-03 | 2018-09-03 | 一种大静不稳定度飞行器弹性稳定控制系统及控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109085847B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109407690A (zh) * | 2018-12-27 | 2019-03-01 | 湖北航天飞行器研究所 | 一种飞行器稳定控制方法 |
CN110750836A (zh) * | 2019-10-15 | 2020-02-04 | 北京电子工程总体研究所 | 一种基于稳态锥动的飞行器俯仰阻尼力矩计算方法 |
CN111273056A (zh) * | 2020-04-02 | 2020-06-12 | 山东创惠电子科技有限责任公司 | 一种不采用高度测量的高速飞行器的攻角观测方法 |
CN112068581A (zh) * | 2020-09-11 | 2020-12-11 | 中国运载火箭技术研究院 | 飞行器复合控制方法、控制装置及存储介质 |
CN112666959A (zh) * | 2020-11-30 | 2021-04-16 | 航天科工火箭技术有限公司 | 一种运载火箭姿态失稳状态下的姿态稳定控制方法 |
CN115202389A (zh) * | 2022-06-27 | 2022-10-18 | 中国航天空气动力技术研究院 | 一种减弱大流量轨控喷流干扰力矩的控制方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5170969A (en) * | 1988-11-23 | 1992-12-15 | The Boeing Company | Aircraft rudder command system |
JP2972731B1 (ja) * | 1998-07-09 | 1999-11-08 | 三菱電機株式会社 | 誘導飛しょう体 |
CN102081355A (zh) * | 2011-01-27 | 2011-06-01 | 西北工业大学 | 静不稳定飞行器等价稳定裕度的飞行试验鲁棒确定方法 |
CN104267732A (zh) * | 2014-09-29 | 2015-01-07 | 哈尔滨工业大学 | 基于频域分析的挠性卫星高稳定度姿态控制方法 |
CN105629725A (zh) * | 2014-10-31 | 2016-06-01 | 北京临近空间飞行器系统工程研究所 | 一种后缘舵滑翔飞行器的弹性运动建模方法 |
CN106125571A (zh) * | 2016-03-18 | 2016-11-16 | 中国运载火箭技术研究院 | 一种巡航导弹的速度控制建模方法 |
CN108445895A (zh) * | 2018-02-05 | 2018-08-24 | 天津大学 | 用于倾转式三旋翼无人机位置控制的鲁棒控制方法 |
-
2018
- 2018-09-03 CN CN201811021187.8A patent/CN109085847B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5170969A (en) * | 1988-11-23 | 1992-12-15 | The Boeing Company | Aircraft rudder command system |
JP2972731B1 (ja) * | 1998-07-09 | 1999-11-08 | 三菱電機株式会社 | 誘導飛しょう体 |
CN102081355A (zh) * | 2011-01-27 | 2011-06-01 | 西北工业大学 | 静不稳定飞行器等价稳定裕度的飞行试验鲁棒确定方法 |
CN104267732A (zh) * | 2014-09-29 | 2015-01-07 | 哈尔滨工业大学 | 基于频域分析的挠性卫星高稳定度姿态控制方法 |
CN105629725A (zh) * | 2014-10-31 | 2016-06-01 | 北京临近空间飞行器系统工程研究所 | 一种后缘舵滑翔飞行器的弹性运动建模方法 |
CN106125571A (zh) * | 2016-03-18 | 2016-11-16 | 中国运载火箭技术研究院 | 一种巡航导弹的速度控制建模方法 |
CN108445895A (zh) * | 2018-02-05 | 2018-08-24 | 天津大学 | 用于倾转式三旋翼无人机位置控制的鲁棒控制方法 |
Non-Patent Citations (1)
Title |
---|
王青 等: "降低频带祸合的强静不稳定弹性飞行器控制", 《弹箭与制导学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109407690A (zh) * | 2018-12-27 | 2019-03-01 | 湖北航天飞行器研究所 | 一种飞行器稳定控制方法 |
CN110750836A (zh) * | 2019-10-15 | 2020-02-04 | 北京电子工程总体研究所 | 一种基于稳态锥动的飞行器俯仰阻尼力矩计算方法 |
CN111273056A (zh) * | 2020-04-02 | 2020-06-12 | 山东创惠电子科技有限责任公司 | 一种不采用高度测量的高速飞行器的攻角观测方法 |
CN112068581A (zh) * | 2020-09-11 | 2020-12-11 | 中国运载火箭技术研究院 | 飞行器复合控制方法、控制装置及存储介质 |
CN112068581B (zh) * | 2020-09-11 | 2023-11-17 | 中国运载火箭技术研究院 | 飞行器复合控制方法、控制装置及存储介质 |
CN112666959A (zh) * | 2020-11-30 | 2021-04-16 | 航天科工火箭技术有限公司 | 一种运载火箭姿态失稳状态下的姿态稳定控制方法 |
CN115202389A (zh) * | 2022-06-27 | 2022-10-18 | 中国航天空气动力技术研究院 | 一种减弱大流量轨控喷流干扰力矩的控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109085847B (zh) | 2021-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109085847A (zh) | 一种大静不稳定度飞行器弹性稳定控制系统及控制方法 | |
CN109543271B (zh) | 一种应用于共轴带推力桨高速直升机的配平操纵设计方法 | |
Salazar-Cruz et al. | Trajectory tracking for a four rotor mini-aircraft | |
CN103307938A (zh) | 一种旋转弹气动参数获取方法 | |
CN102915037B (zh) | 基于混合控制的四旋翼无人机增稳控制方法 | |
CN108763772B (zh) | 一种固定翼无人机纵向姿态控制器设计方法 | |
CN108037662A (zh) | 一种基于积分滑模障碍李雅普诺夫函数的四旋翼飞行器输出受限反步控制方法 | |
CN107515612A (zh) | 基于侧喷流控制的弹性振动抑制方法 | |
CN107783429A (zh) | 一种通用飞行模拟器自动飞行控制仿真系统 | |
CN104656661B (zh) | 一种公务机下降着陆控制方法 | |
CN110673623B (zh) | 一种基于双环pd控制算法控制的四旋翼无人机着陆方法 | |
CN107066640A (zh) | 基于保护映射的变构型飞行器模型迭代设计方法 | |
CN108750120A (zh) | 一种植保无人飞机喷头主动平衡装置及调节方法 | |
CN113885543B (zh) | 一种复杂海况干扰下的航行体自抗扰横滚控制器的设计方法 | |
CN108958038A (zh) | 一种适应飞行器推力离散特性的控制参数调整方法 | |
CN110109473A (zh) | 可应用于旋转飞行器的非线性滚转稳定控制系统及方法 | |
EP3931086A1 (en) | A coaxial rotorcraft system and a method for controlling the same | |
CN113777932A (zh) | 一种基于Delta算子的四旋翼自适应滑模容错控制方法 | |
CN107992071A (zh) | 尾座式无人机纵向姿态双模糊控制系统及方法 | |
CN107831653B (zh) | 一种抑制参数摄动的高超声速飞行器指令跟踪控制方法 | |
Thai et al. | Robust design for the roll-channel autopilot of a canard-guided dual-spin projectile | |
JPH0933197A (ja) | 混合ミサイル自動操縦器 | |
CN109992003A (zh) | 鲁棒性滚转角速率控制方法及系统 | |
CN109484676A (zh) | 一种垂直起降火箭在线轨迹规划的等效姿态控制处理方法 | |
Roy et al. | Robust control for longitudinal and lateral dynamics of small scale helicopter |
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 |