CN102081350A - 静不稳定飞行器等价气动伺服弹性鲁棒稳定性的飞行试验确定方法 - Google Patents
静不稳定飞行器等价气动伺服弹性鲁棒稳定性的飞行试验确定方法 Download PDFInfo
- Publication number
- CN102081350A CN102081350A CN 201110030527 CN201110030527A CN102081350A CN 102081350 A CN102081350 A CN 102081350A CN 201110030527 CN201110030527 CN 201110030527 CN 201110030527 A CN201110030527 A CN 201110030527A CN 102081350 A CN102081350 A CN 102081350A
- Authority
- CN
- China
- Prior art keywords
- omega
- transfer function
- stability
- flight test
- loop transfer
- 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
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种静不稳定飞行器等价气动伺服弹性鲁棒稳定性的飞行试验确定方法,用于解决现有的飞行试验确定方法建模误差大的技术问题。技术方案是通过扫频飞行试验可以得到静不稳定飞行器多回路***的闭环传递函数频率特性矩阵,按照闭环传递函数与开环传递函数之间的对应关系,建立闭环传递函数与ASE稳定性之间的关系式,并将***闭环传递函数频率特性矩阵进行特征分解,根据***临界稳定的条件,得到多回路稳定裕度分析标量方程式,从而计算整体回路的ASE稳定性,大大简化了问题难度。
Description
技术领域
本发明涉及一种飞行试验确定方法,特别涉及静不稳定飞行器等价气动伺服弹性(ASE)鲁棒稳定性的飞行试验确定方法。
背景技术
气动伺服弹性稳定性分析的研究目的是依靠气动伺服弹性模型,并结合试飞试验数据,分析并预测气动伺服弹性***的稳定性,最终确定飞机稳定飞行的飞行包线。传统的稳定性分析方法是借助飞行试验,逐个分析***的每个回路,利用经典的Bode图或Nyquist图获得***的幅值裕度和相角裕度。然而,对于现代飞机而言,由于采用了多回路的电传操纵***,从而导致其气动伺服弹性***也存在较为复杂的耦合关系,上述的SISO***分析方法并不适用,因此,如何有效的设计、分析以及预测多回路ASE***的稳定性就成了当前亟待研究的问题。
多回路气动伺服弹性***稳定性分析是伴随我国自主研制新型飞机出现的新问题,也是当今国际上飞机设计及试飞中普遍关注的问题。研究多回路ASE***的分析方法可以满足我国新型电传操纵飞机的设计及试飞的要求,以便在飞机设计之初就能够对ASE稳定性进行有效评估,指导飞机总体设计;同时在飞机试飞定型阶段,运用该方法设计试飞方案,从试验角度验证实际ASE***的稳定性。在国际上同类研究中,美国最早采用鲁棒分析方法分析并预测ASE***的飞行包线,并在F18试验机的ASE试飞验证中得到应用,遗憾的是该方法对气动弹性模型的依赖度过高,需要成熟的建模方法和准确的标称模型;而且,鲁棒分析方法的保守性及多次迭代计算也是制约其应用的障碍;特别是由于现代战斗机放宽了静稳定性,在飞行试验时,必须使飞机带有增稳***,否则容易出现飞行事故;因此,通常都是通过闭环飞行试验并结合控制器模型的方法近似得到飞行器的开环特性,由此带来了较大的建模误差。
发明内容
为了克服现有的飞行试验确定方法建模误差大的不足,本发明提供一种静不稳定飞行器等价气动伺服弹性鲁棒稳定性的飞行试验确定方法,该方法通过扫频飞行试验可以得到静不稳定飞行器多回路***的闭环传递函数频率特性矩阵,按照闭环传递函数与开环传递函数之间的对应关系,建立闭环传递函数与ASE稳定性之间的关系式,并将***闭环传递函数频率特性矩阵进行特征分解,根据***临界稳定的条件,得到多回路稳定裕度分析标量方程式,从而计算整体回路的ASE稳定性。
本发明解决其技术问题采用的技术方案是,一种静不稳定飞行器等价气动伺服弹性鲁棒稳定性的飞行试验确定方法,其特点是包括以下步骤:
1、通过扫频飞行试验,获得静不稳定飞行器多回路***含有不确定量的闭环传递函数频率特性矩阵Φ(jω)=G(jω)[I+G(jω)]-1;为了确定***的相位裕度和幅值裕度,定义:Yp(jω)=Kpe-τjω,并串联接至等价开环***;
式中,j为虚数符号,ω表示频率,Φ(jω)为闭环传递函数频率特性矩阵,G(jω)为开环传递函数频率特性矩阵,I为单位矩阵,Kp为每条回路的附加的增益,τ为每条回路的附加时间滞后;
2、当***临界稳定时,复变量s的实部为零,令s=jω为纯虚数,满足以下行列式关系|(Kpe-τjω-1)Φ(jω)+I|=0,
对Φ(jω)进行特征分解,得:Φ(jω)=T(jω)D(jω)T-1(jω)
式中,T(jω)为线性变换矩阵,D(jω)为约当阵;
3、取τ=0,选定颤振频率段[ω0,ω1],求解满足以下等式:
时得到的最小ω值和回路幅值裕度Kp值,最小的Kp为整个***在选定颤振频率段[ω0,ω1]的幅值裕度,δ>0为不确定性影响的估计值,di(jω)为矩阵D(jω)的第i行第i列元素;
4、相位裕度转化为Kp=1时在选定颤振频率段[ω0,ω1]上计算
时得到的最小ω值和最小的τ,此时的ωτ为整体***的相位裕度,λ>0为不确定性影响的估计值。
本发明的有益效果是:由于通过扫频飞行试验直接获得了静不稳定飞行器多回路***的闭环传递函数频率特性矩阵,按照闭环传递函数与开环传递函数之间的对应关系,建立了闭环传递函数与ASE稳定性之间的关系式,根据***临界稳定的条件,得到了通过闭环传递函数频率特性矩阵分析ASE稳定性的两个等价行列式,并将***闭环传递函数频率特性矩阵进行特征分解,根据***临界稳定的条件,得到多回路稳定裕度分析标量方程式,从而计算整体回路的ASE稳定性,大大简化了问题难度。
下面结合附图和实施例对本发明作详细说明。
附图说明
附图是本发明静不稳定飞行器等价气动伺服弹性鲁棒稳定性的飞行试验确定方法的流程图。
具体实施方式
参照附图,详细说明本发明。
1通过扫频飞行试验,获得静不稳定飞行器多回路***的闭环传递函数频率特性
矩阵Φ(jω)=G(jω)[I+G(jω)]-1。
式中,j是虚数符号,ω表示频率,I是单位矩阵。为了便于说明,假定通过扫频飞行试验直接获得了静不稳定飞行器2输入2输出***的G(jω)假定为
为了确定***的相位裕度和幅值裕度,定义:Yp(jω)=Kpe-τjω,并串联接至等价开环***;
2、当***临界稳定时,复变量s的实部为零,令s=jω为纯虚数,满足以下行列式关系
|(Kpe-τjω-1)Φ(jω)+I|=0,对Φ(jω)进行特征分解,得:Φ(jω)=T(jω)D(jω)T-1(jω)
3、取τ=0,选定颤振频率段[5,50],求解满足以下等式:
时得到在ω=17.68000时,回路幅值裕度Kp=189.78为整个***在选定颤振频率段[5,50]的幅值裕度,δ>0.43为不确定性影响的估计值,
4、相位裕度转化为Kp=1时在选定颤振频率段[0.5,50]上计算
时得到的最小ω=1.225值和最小的τ=2.02,整体***的相位裕度为141.8°,λ>0.025为不确定性影响的估计值。
说明飞行器具有ASE稳定性。
Claims (1)
1.一种静不稳定飞行器等价气动伺服弹性鲁棒稳定性的飞行试验确定方法,其特征在于包括下述步骤:
(a)通过扫频飞行试验,获得静不稳定飞行器多回路***含有不确定量的闭环传递函数频率特性矩阵Φ(jω)=G(jω)[I+G(jω)]-1;为了确定***的相位裕度和幅值裕度,定义:Yp(jω)=Kpe-τjω,并串联接至等价开环***;
式中,j为虚数符号,ω表示频率,Φ(jω)为闭环传递函数频率特性矩阵,G(jω)为开环传递函数频率特性矩阵,I为单位矩阵,Kp为每条回路的附加的增益,τ为每条回路的附加时间滞后;
(b)当***临界稳定时,复变量s的实部为零,令s=jω为纯虚数,满足以下行列式关系|(Kpe-τjω-1)Φ(jω)+I|=0,
对Φ(jω)进行特征分解,得:Φ(jω)=T(jω)D(jω)T-1(jω)
式中,T(jω)为线性变换矩阵,D(jω)为约当阵;
(c)取τ=0,选定颤振频率段[ω0,ω1],求解满足以下等式:
时得到的最小ω值和回路幅值裕度Kp值,最小的Kp为整个***在选定颤振频率段[ω0,ω1]的幅值裕度,δ>0为不确定性影响的估计值,di(jω)为矩阵D(jω)的第i行第i列元素;
(d)相位裕度转化为Kp=1时在选定颤振频率段[ω0,ω1]上计算
时得到的最小ω值和最小的τ,此时的ωτ为整体***的相位裕度,λ>0为不确定性影响的估计值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110030527A CN102081350B (zh) | 2011-01-27 | 2011-01-27 | 静不稳定飞行器等价气动伺服弹性鲁棒稳定性的飞行试验确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110030527A CN102081350B (zh) | 2011-01-27 | 2011-01-27 | 静不稳定飞行器等价气动伺服弹性鲁棒稳定性的飞行试验确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102081350A true CN102081350A (zh) | 2011-06-01 |
CN102081350B CN102081350B (zh) | 2012-09-26 |
Family
ID=44087357
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110030527A Expired - Fee Related CN102081350B (zh) | 2011-01-27 | 2011-01-27 | 静不稳定飞行器等价气动伺服弹性鲁棒稳定性的飞行试验确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102081350B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102929128A (zh) * | 2012-10-10 | 2013-02-13 | 西北工业大学 | 飞行器不确定模型的控制器设计方法 |
CN102929129A (zh) * | 2012-10-10 | 2013-02-13 | 西北工业大学 | 飞行器多个时间滞后时变模型的控制器设计方法 |
CN102929140A (zh) * | 2012-10-10 | 2013-02-13 | 西北工业大学 | 飞行器时间滞后模型逼近及控制器设计方法 |
CN102929136A (zh) * | 2012-10-10 | 2013-02-13 | 西北工业大学 | 飞行器多个时间滞后模型逼近及控制器设计方法 |
CN102929145A (zh) * | 2012-10-10 | 2013-02-13 | 西北工业大学 | 飞行器多个时间滞后时变模型逼近及控制器设计方法 |
CN102929130A (zh) * | 2012-10-10 | 2013-02-13 | 西北工业大学 | 一种鲁棒飞行控制器设计方法 |
CN103777523A (zh) * | 2014-02-28 | 2014-05-07 | 西安费斯达自动化工程有限公司 | 飞行器多回路模型簇复合pid鲁棒控制器设计方法 |
CN103809448A (zh) * | 2014-02-28 | 2014-05-21 | 西安费斯达自动化工程有限公司 | 飞行器多回路模型簇颤振抑制复合根轨迹补偿鲁棒控制器设计方法 |
CN103823367A (zh) * | 2014-02-28 | 2014-05-28 | 西安费斯达自动化工程有限公司 | 纵向飞行模型簇颤振抑制复合频率鲁棒控制器设计方法 |
CN103823377A (zh) * | 2014-02-28 | 2014-05-28 | 西安费斯达自动化工程有限公司 | 纵向飞行模型簇颤振抑制复合根轨迹补偿鲁棒控制器设计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1477606A (zh) * | 2002-08-23 | 2004-02-25 | 中国科学技术大学 | 旋转式倒立摆 |
CN200976108Y (zh) * | 2006-11-24 | 2007-11-14 | 中国科学院沈阳自动化研究所 | 一种小型无人直升机自主飞行控制*** |
CN101192064A (zh) * | 2006-11-24 | 2008-06-04 | 中国科学院沈阳自动化研究所 | 小型无人直升机自主飞行控制*** |
-
2011
- 2011-01-27 CN CN201110030527A patent/CN102081350B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1477606A (zh) * | 2002-08-23 | 2004-02-25 | 中国科学技术大学 | 旋转式倒立摆 |
CN200976108Y (zh) * | 2006-11-24 | 2007-11-14 | 中国科学院沈阳自动化研究所 | 一种小型无人直升机自主飞行控制*** |
CN101192064A (zh) * | 2006-11-24 | 2008-06-04 | 中国科学院沈阳自动化研究所 | 小型无人直升机自主飞行控制*** |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102929136B (zh) * | 2012-10-10 | 2015-03-18 | 西北工业大学 | 飞行器多个时间滞后模型逼近及控制器设计方法 |
CN102929128B (zh) * | 2012-10-10 | 2015-04-08 | 西北工业大学 | 飞行器不确定模型的控制器设计方法 |
CN102929130B (zh) * | 2012-10-10 | 2015-06-17 | 西北工业大学 | 一种鲁棒飞行控制器设计方法 |
CN102929136A (zh) * | 2012-10-10 | 2013-02-13 | 西北工业大学 | 飞行器多个时间滞后模型逼近及控制器设计方法 |
CN102929145A (zh) * | 2012-10-10 | 2013-02-13 | 西北工业大学 | 飞行器多个时间滞后时变模型逼近及控制器设计方法 |
CN102929130A (zh) * | 2012-10-10 | 2013-02-13 | 西北工业大学 | 一种鲁棒飞行控制器设计方法 |
CN102929129B (zh) * | 2012-10-10 | 2015-04-29 | 西北工业大学 | 飞行器多个时间滞后时变模型的控制器设计方法 |
CN102929129A (zh) * | 2012-10-10 | 2013-02-13 | 西北工业大学 | 飞行器多个时间滞后时变模型的控制器设计方法 |
CN102929140A (zh) * | 2012-10-10 | 2013-02-13 | 西北工业大学 | 飞行器时间滞后模型逼近及控制器设计方法 |
CN102929128A (zh) * | 2012-10-10 | 2013-02-13 | 西北工业大学 | 飞行器不确定模型的控制器设计方法 |
CN102929145B (zh) * | 2012-10-10 | 2015-03-18 | 西北工业大学 | 飞行器多个时间滞后时变模型逼近及控制器设计方法 |
CN102929140B (zh) * | 2012-10-10 | 2015-03-11 | 西北工业大学 | 飞行器时间滞后模型逼近及控制器设计方法 |
CN103809448B (zh) * | 2014-02-28 | 2016-05-25 | 西安费斯达自动化工程有限公司 | 飞行器多回路模型簇颤振抑制复合根轨迹补偿鲁棒控制器设计方法 |
CN103809448A (zh) * | 2014-02-28 | 2014-05-21 | 西安费斯达自动化工程有限公司 | 飞行器多回路模型簇颤振抑制复合根轨迹补偿鲁棒控制器设计方法 |
CN103777523A (zh) * | 2014-02-28 | 2014-05-07 | 西安费斯达自动化工程有限公司 | 飞行器多回路模型簇复合pid鲁棒控制器设计方法 |
CN103823367A (zh) * | 2014-02-28 | 2014-05-28 | 西安费斯达自动化工程有限公司 | 纵向飞行模型簇颤振抑制复合频率鲁棒控制器设计方法 |
CN103777523B (zh) * | 2014-02-28 | 2016-01-06 | 西安费斯达自动化工程有限公司 | 飞行器多回路模型簇复合pid鲁棒控制器设计方法 |
CN103823367B (zh) * | 2014-02-28 | 2016-04-06 | 西安费斯达自动化工程有限公司 | 纵向飞行模型簇颤振抑制复合频率鲁棒控制器设计方法 |
CN103823377A (zh) * | 2014-02-28 | 2014-05-28 | 西安费斯达自动化工程有限公司 | 纵向飞行模型簇颤振抑制复合根轨迹补偿鲁棒控制器设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102081350B (zh) | 2012-09-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102081350B (zh) | 静不稳定飞行器等价气动伺服弹性鲁棒稳定性的飞行试验确定方法 | |
CN102081349B (zh) | 多输入-多输出等价气动伺服弹性鲁棒稳定性的飞行试验确定方法 | |
US20220309122A1 (en) | Iterative algorithm for aero-engine model based on hybrid adaptive differential evolution | |
CN103472820B (zh) | 一种基于偏最小二乘算法的推进***故障诊断方法 | |
Kobayashi et al. | Integration of on-line and off-line diagnostic algorithms for aircraft engine health management | |
CN101859146A (zh) | 一种基于预测滤波和经验模态分解的卫星故障预测方法 | |
He et al. | Preliminary exploration on digital twin for power systems: Challenges, framework, and applications | |
CN103345546B (zh) | 频率轨迹与粒子群算法相结合的调速器参数辨识方法 | |
CN103676941A (zh) | 基于运动学和动力学模型的卫星控制***故障诊断方法 | |
CN102176118A (zh) | 多输入-多输出飞行器等价稳定裕度鲁棒确定方法 | |
CN104571087B (zh) | 一种噪声影响下航天器控制***可诊断性确定方法 | |
CN104019817A (zh) | 一种用于卫星姿态估计的范数约束强跟踪容积卡尔曼滤波方法 | |
CN101813920B (zh) | 电站汽轮机组温度传感器虚拟冗余方法 | |
CN102081355B (zh) | 静不稳定飞行器等价稳定裕度的飞行试验鲁棒确定方法 | |
Grauer | Real-time data-compatibility analysis using output-error parameter estimation | |
Novara et al. | Control of MIMO nonlinear systems: A data-driven model inversion approach | |
Shimizu et al. | Output-based error estimation for chaotic flows using reduced-order modeling | |
Zhang et al. | Performance analysis of digital flight control systems with rollback error recovery subject to simulated neutron-induced upsets | |
Poudel et al. | Comparative analysis of state and parameter estimation techniques for power system frequency dynamics | |
Yan et al. | Simplified Markov model for reliability analysis of phased-mission system using states merging method | |
Liu et al. | Fault diagnosability qualitative analysis of spacecraft based on temporal fault signature matrix | |
Rahman et al. | Convergence of the fast state estimation for power systems | |
Jianwang et al. | Combing instrumental variable and variance matching for aircraft flutter model parameters identification | |
Suh et al. | Robust modal filtering and control of the X-56A model with simulated fiber optic sensor failures | |
CN106507990B (zh) | 静不稳定飞行器多回路ase稳定性飞行试验直接确定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120926 Termination date: 20200127 |
|
CF01 | Termination of patent right due to non-payment of annual fee |