CN110287505A - 飞行器稳定性分析方法 - Google Patents
飞行器稳定性分析方法 Download PDFInfo
- Publication number
- CN110287505A CN110287505A CN201910212456.7A CN201910212456A CN110287505A CN 110287505 A CN110287505 A CN 110287505A CN 201910212456 A CN201910212456 A CN 201910212456A CN 110287505 A CN110287505 A CN 110287505A
- Authority
- CN
- China
- Prior art keywords
- aircraft
- stability
- mode
- aerodynamic
- rank
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提出一种飞行器稳定性分析方法,依据飞行器的外形建立空间网格,并对飞行器的选定状态进行定常流场计算,得到气动力系数和系数矩阵;获得飞行器的模态频率和振型;计算流场物面的振型;根据飞行器的结构特征设计模态训练信号;得到飞行器的弹性振动表达公式;得到各模态的广义气动力模型;得到离散空间内的气动力状态方程;提取出不同减缩频率的广义气动力模型;得到气动/结构/控制耦合的闭环系统运行方程;将闭环系统运行方程转换为传递函数,并进行稳定性分析,从而对飞行器稳定性进行分析。本发明的方法将飞行器的弹性变形进入到对控制系统的影响中,提高了控制系统稳定性分析的准确性。
Description
技术领域
本发明涉及飞行器气动/结构/控制耦合的稳定性分析,属于飞行器气动弹性 技术领域。
背景技术
高速飞行器特别是吸气式高速飞行器对于结构轻质化提出了很高要求,结构 刚度必然下降,同时长时高速飞行的热环境严酷,高温热环境使得结构刚度进一 步发生变化,飞行过程中机体在气动力、热作用下必然会发生弹性变形和振动, 装在飞行器上的传感器不仅接收到飞行器的刚体运动信号,同时接受到结构弹性 振动的信号,传感器中这些信号均会传递给控制系统,而控制系统通过伺服传动 给舵面的指令中就增加了这部分弹性振动频率的偏转运动,从而增加了一部分舵 面非定常气动力和惯性力,这是因控制系统的反馈而额外作用于飞行的控制力, 形成气动/结构/控制耦合的闭环回路系统。当闭环回路出现不稳定时,受到外界 非定常气动力的扰动,系统会迅速发散,造成飞行器破坏。
对于此问题的分析,以往采用简化的平板面元法的气动力模型与结构、控制 耦合进行频域范围内的稳定性分析,但平板面元法计算非定常气动力,具有局限 性,无法反应高速飞行器的复杂流动现象和粘性效应,导致对飞行器进行控制稳 定性的分析结果不准确。
发明内容
本发明的目的在于克服现有技术不足,提供一种控制耦合的高速飞行器稳定 性分析方法。
本发明的技术解决方案:
一种飞行器稳定性分析方法,步骤如下:
依据飞行器的外形建立空间网格,并对飞行器的选定状态进行定场流场计 算,得到气动力系数和系数矩阵;
依据飞行器结构和材料建立结构有限元模型,并进行给定飞行状态和温度场 下的热模态分析,获得飞行器的模态频率和振型;
将得到的模态振型插值到获得的空间网格的物面上,得到流场物面的振型;
根据飞行器的结构特征设计模态训练信号,所述的模态训练信号由白噪声进 行滤波后,根据一定的采样频率进行采样得到的。
将设计好的模态训练信号与得到的流场物面振型相乘并叠加,得到飞行器的 弹性振动表达公式;
利用得到的弹性振动公式模拟流场网格的运动,在获得的定常流场计算的基 础上进行非定常计算得到流场的压力分布,并根据压力分布得到各模态的广义气 动力模型;
根据得到的非定常广义气动力模型和多输入多输出的系统描述方程,得到连 续系统状态空间的模型为:式中:xa(t)表示状态向 量,xa(k)=[fa(k-1),…,fa(k-na),ξ(k-1),…,ξ(k-nb+1)]T,其中,fa(k-1)表示第k-1 个模态气动力系数,ξ(k-1)表示第k-1个模态位移,na表示fa(t)的延迟阶数,nb 表示ξ(t)的延迟阶数,fa(t)表示模态气动力系数,ξ(t)表示模态位移,Aa、Ba、Ca、 Da为系数矩阵,fa0为步骤一中定常流场计算得到的气动力系数。
将得到的连续系统状态空间模型转换为传递函数模型,获得气动结构耦合的 开环系统的幅相频特性,提取出不同减缩频率的广义气动力模型;
根据得到的各个广义气动力模型,得到气动/结构/控制耦合的闭环系统运行 方程为:
其中,q模态坐标向量,δ为控制面刚体偏转坐标向量,ρ为大气密度,V为 飞行速度,Mqq为各阶弹性模态对应的广义质量矩阵和Mqδ为舵偏1度产生的广 义质量矩阵,Cqq为各阶弹性模态对应的广义阻尼矩阵,Kqq为各阶弹性模态对应 的广义刚度矩阵,Qqq是各阶弹性模态对应的非定常气动力矩阵,Qqδ为舵偏1度 产生的非定常气动力矩阵。
将得到的方程转换为传递函数,并进行稳定性分析,得到闭环系统的稳定裕 度,依据稳定裕度对飞行器稳定性进行分析判断。
进一步的,所述的飞行器的模态频率包括飞行器垂直于翼展方向和平行于翼 展方向的前三阶弯曲频率、扭转频率和飞行器舵面的第一阶弯曲频率和第一阶旋 转频率。
进一步的,所述的白噪声的采样频率为其中b为飞行器的参考 长度,a为飞行状态下的当地声速,Ma为马赫数,dt为非定常气动力计算的时 间步长。
进一步的,所述的减缩频率为其中,V是飞行器飞行速度,ω是飞 行器振动的圆周频率。
进一步的,所述的每一个结构模态生成一组训练信号。
本发明与现有技术相比的有益效果:
(1)本发明通过模态训练信号的设计,采样频率根据飞行器的长度、飞行 器的速度和当地的声速等因素,使飞行器的模态训练信号更贴合当前飞行器的状 态,通过模态训练信号与流场物面振型共同计算得到飞行器的弹性振动表达公 式,使飞行器的弹性变形进入到对控制系统的影响中,提高了控制系统稳定性分 析的准确性;
(2)本发明通过将结构/气动力/控制耦合,将飞行器的结构、气动力和控 制系统等对飞行器的稳定性有影响的因素耦合到一起,对飞行器的稳定性进行分 析,这种分析方法不仅能适用于高超声速复杂飞行器的流场模拟、考虑了粘性效 应,同时可直接转换到状态空间方程进行稳定性分析,提高了对高速复杂外形的 适用性和控制系统稳定性分析结果的准确性。
附图说明
图1为本发明飞行器稳定性分析方法的流程图;
图2为本发明飞行器结构坐标系图;
图3为本发明某飞行器第一阶模态的滤波信号;
图4为本发明第一阶模态的广义气动力计算结果;
图5为本发明前4阶模态气动力在不同减缩频率下的幅值曲线;
图6为本发明前4阶模态气动力在不同减缩频率下的相位曲线;
图7为本发明飞行器稳定性的奈奎斯特图;
图8为本发明飞行器稳定性的伯德幅值图;
图9为本发明飞行器稳定性的伯德相角图。
具体实施方式
下面结合具体实例及附图对本发明进行详细说明。
如图1所示,一种飞行器稳定性分析方法,步骤如下:
步骤一,依据飞行器的外形划分流场CFD计算的空间网格,并对飞行器的选 定状态进行定场流场计算,此为本领域公知技术;
步骤二,依据飞行器结构和材料建立结构有限元模型,此为本领域公知技术, 优选的,如图2所示,结构有限元模型与流场计算的空间网格保持坐标系一致、 单位统一,流场物面与结构外表面的几何外形一致;对飞行器有限元模型进行给 定飞行状态和温度场下的热模态分析,获得飞行器的模态频率和振型,此为本领 域公知技术;
在一个实施例中,飞行器的模态频率包括飞行器垂直于翼展方向和平行于翼 展方向的前三阶弯曲频率、扭转频率和飞行器舵面的第一阶弯曲频率和第一阶旋 转频率,具体选择前几阶的频率,根据精度的需要进行选择。
步骤三,将步骤二得到的模态振型插值到步骤一获得的空间网格的物面上, 得到流场物面的振型,此方法为公知技术;
在一个实施例中,插值的方法采用RBF径向基函数插值的方法或其他现有的 插值计算方法。
步骤四,根据飞行器的结构特征设计非定常气动力降阶模型需要的训练信 号;
在一个实施例中,通过Matlab生成一段随机白噪声,对白噪声信号进行滤 波采样,采样频率为其中b为飞行器的参考长度,a为飞行状态下的 当地声速,dt为非定常气动力计算的时间步长;每一阶模态生成一组训练信号, 在一个实施例中,第一阶模态的训练信号如图3所示。
在一个实施例中,第一阶模态的频率为30Hz,参考长度b=1,飞行马赫数 Ma=0.8,飞行速度V=266m/s,则采样滤波的频率为0.70683。
步骤五,将步骤四设计好的模态训练信号与步骤三得到的流场物面振型相乘 并叠加,得到飞行器的弹性振动表达公式:其中qi为模态坐标向量, N代表模态的阶数,Φi代表各阶模态振型,此为本领域公知技术;
步骤六,计算步骤五得到的弹性振动公式模拟流场网格的运动,在步骤一获 得的定常流场计算的基础上进行非定常计算得到流场的压力分布,并根据压力分 布得到各模态的广义气动力模型其中p代表 流场的压力分布,ds代表流场物面单元面积,此为本领域公知技术,本发明第 一阶模态的广义气动力计算结果如图4所示;
步骤七,根据步骤六得到的非定常广义气动力模型和多输入多输出的系统描 述方程,通过最小二乘方法的差分方程的参数辨识过程,得到离散空间内的气动 力状态方程,此为本领域公知技术。
在一个实施例中,多输入多输出系统描述方程如下:
其中y(k)为系统输出量的第k次观测值,y(k-1)为系统输出量的第k-1次观 测值,以此类推;u(k)为系统第k个输入量,u(k-1)为系统第k-1个输入量,以 此类推;e(k)为零均值的随机噪声,Ai和Bi为待辨识的系数矩阵,na和nb分别 为输出和输入的延迟阶数。
采用基于最小二乘方法的差分方程进行参数辨识,对于非定常气动力模型, 令模态位移ξ=u,模态气动力系数fa=y,为了便于进行气动/结构/控制耦合的 稳定性分析,将上述多输入多输出系统描述方程转化为状态空间模型,定义状态 向量xa(k):
xa(k)=[fa(k-1),…,fa(k-na),ξ(k-1),…,ξ(k-nb+1)]T
则离散空间内的气动力状态方程和输出方程可以写为:
其中,为系数矩阵,fa0为第一步定常流场计算得到的气动 力系数。
式中,Ai和Bi为待辨识的系数矩阵,i=0、1、2、…、na,I为单 位矩阵。
步骤八,根据步骤七中得到的连续系统状态空间模型转换为传递函数模型, 获得气动结构耦合的开环系统的幅值和相位特性,提取出不同减缩频率的广义气 动力模型,减缩频率为其中b是参考长度,V是飞行速度,ω是振动的圆 周频率,本实施例前4阶模态气动力在不同减缩频率下的幅值曲线如图5所示, 前4阶模态气动力在不同减缩频率下的相位曲线如图6所示;
步骤九,根据步骤八中得到的各个广义气动力模型Q,得到气动/结构/控制 耦合的闭环系统运行方程为:
其中,q模态坐标向量,δ为控制面刚体偏转坐标向量,ρ为大气密度,V为 飞行速度,Mqq为各阶弹性模态对应的广义质量矩阵和Mqδ为舵偏1度产生的广 义质量矩阵,Cqq为各阶弹性模态对应的广义阻尼矩阵,Kqq为各阶弹性模态对应 的广义刚度矩阵,Qqq是各阶弹性模态对应的非定常气动力矩阵,Qqδ为舵偏1度 产生的非定常气动力矩阵。
步骤十,将步骤九得到的方程转换为频率表达:
其中,Qq,Δδ为由于弹性引起的附加舵偏Δδ产生的控制力矩阵,因此从飞行器的弹性振动到舵偏的传递函数可表示为:
进行稳定性分析,得到闭环系统的稳定裕度,从而对飞行器稳定性进行分析。
在一个实施例中,通过稳定性分析得到奈奎斯特图和伯德图,据此判断飞行 器的稳定性,如图7、8、9所示,由于奈奎斯特图中包围(-1,0)点,且伯德图 中幅值超过了0,因此飞行器气动/结构/控制耦合特性是不稳定的。
在其他的实施例中,步骤可以发生改变,但都在本发明的保护范围内。
本发明未详细说明部分为本领域技术人员公知技术。
Claims (6)
1.一种飞行器稳定性分析方法,其特征在于:步骤如下,
依据飞行器的外形建立空间网格,并对飞行器的选定状态进行定常流场计算,得到气动力系数和系数矩阵;
建立飞行器结构有限元模型,并进行给定飞行状态和温度场下的热模态分析,获得飞行器的模态频率和振型;
计算流场物面的振型;
根据飞行器的结构特征设计模态训练信号,所述的模态训练信号由白噪声进行滤波后,根据一定的采样频率进行采样得到;
将设计好的模态训练信号与流场物面振型相乘并叠加,得到飞行器的弹性振动表达公式;
根据得到的弹性振动公式模拟流场网格的运动,在获得的定常流场计算的基础上进行非定常计算得到流场的压力分布,并根据压力分布得到各模态的广义气动力模型;
根据得到的非定常广义气动力模型和多输入多输出的系统描述方程,得到离散空间内的气动力状态方程;
根据得到的连续系统状态空间模型转换为传递函数模型,获得气动结构耦合的开环系统的幅相频特性,提取出不同减缩频率的广义气动力模型;
根据各个广义气动力模型,得到气动/结构/控制耦合的闭环系统运行方程;
将得到的方程转换为传递函数,并进行稳定性分析,得到闭环系统的稳定裕度,依据稳定裕度对飞行器稳定性进行分析判断。
2.根据权利要求1所述的高超声速飞行器稳定性分析方法,其特征在于:所述的飞行器的模态频率包括飞行器垂直于翼展方向和平行于翼展方向的前三阶弯曲频率、扭转频率和飞行器舵面的第一阶弯曲频率和第一阶旋转频率。
3.根据权利要求1所述的高超声速飞行器稳定性分析方法,其特征在于:步骤三中所述的采样频率为其中b为飞行器的参考长度,a为飞行状态下的当地声速,dt为非定常气动力计算的时间步长。
4.根据权利要求1所述的高超声速飞行器稳定性分析方法,其特征在于:所述的每一个结构模态生成一组训练信号。
5.根据权利要求1所述的高超声速飞行器稳定性分析方法,其特征在于:所述的离散空间内的气动力状态方程为:式中:xa(k)表示状态向量,xa(k)=[fa(k-1),…,fa(k-na),ξ(k-1),…,ξ(k-nb+1)]T,其中,fa(k-1)表示第k-1个模态气动力系数,ξ(k-1)表示第k-1个模态位移,na表示fa(t)的延迟阶数,nb表示ξ(t)的延迟阶数,fa(k)表示模态气动力系数,ξ(k)表示模态位移,为系数矩阵,fa0为步骤一中定常流场计算得到的气动力系数。
6.根据权利要求1所述的高超声速飞行器稳定性分析方法,其特征在于:所述的闭环系统运行方程为:
其中,q模态坐标向量,δ为控制面刚体偏转坐标向量,ρ为大气密度,V为飞行速度,Mqq为各阶弹性模态对应的广义质量矩阵和Mqδ为舵偏1度产生的广义质量矩阵,Cqq为各阶弹性模态对应的广义阻尼矩阵,Kqq为各阶弹性模态对应的广义刚度矩阵,Qqq是各阶弹性模态对应的非定常气动力矩阵,Qqδ为舵偏1度产生的非定常气动力矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910212456.7A CN110287505B (zh) | 2019-03-20 | 2019-03-20 | 飞行器稳定性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910212456.7A CN110287505B (zh) | 2019-03-20 | 2019-03-20 | 飞行器稳定性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110287505A true CN110287505A (zh) | 2019-09-27 |
CN110287505B CN110287505B (zh) | 2020-12-25 |
Family
ID=68001248
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910212456.7A Active CN110287505B (zh) | 2019-03-20 | 2019-03-20 | 飞行器稳定性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110287505B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112528408A (zh) * | 2020-12-11 | 2021-03-19 | 中国直升机设计研究所 | 一种直升机旋翼与机身耦合稳定性建模方法 |
CN113111430A (zh) * | 2021-03-06 | 2021-07-13 | 北京航空航天大学 | 基于非线性气动力降阶的弹性飞机飞行动力学建模方法 |
CN114912202A (zh) * | 2022-05-24 | 2022-08-16 | 大连理工大学 | 宽速域吸气式动力飞行器机体推进一体化耦合控制方法 |
CN116956782A (zh) * | 2023-09-20 | 2023-10-27 | 西北工业大学 | 非线性颤振分析方法 |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7437014B2 (en) * | 2004-04-21 | 2008-10-14 | The Boeing Company | Methods and systems for model reduction and system identification of dynamic systems with multiple inputs |
CN102169046A (zh) * | 2010-12-22 | 2011-08-31 | 北京航空航天大学 | 一种磁悬浮机电设备弹性模态在线测试系统 |
CN103077259A (zh) * | 2011-10-26 | 2013-05-01 | 上海机电工程研究所 | 高超声速导弹多场耦合动力学一体化仿真分析方法 |
CN103400035A (zh) * | 2013-07-29 | 2013-11-20 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种高可信度快速预测飞行器滚转动导数的方法 |
CN103853157A (zh) * | 2014-03-19 | 2014-06-11 | 湖北蔚蓝国际航空学校有限公司 | 一种基于自适应滑模的飞行器姿态控制方法 |
CN104022742A (zh) * | 2014-06-19 | 2014-09-03 | 金陵科技学院 | 基于神经网络观测器的飞行器姿态鲁棒反演容错控制方法 |
CN104049533A (zh) * | 2013-03-11 | 2014-09-17 | 王岩 | 基于隐半马尔可夫模型的网络控制系统控制器设计方法 |
CN104063533A (zh) * | 2014-04-04 | 2014-09-24 | 西安电子科技大学 | 一种机翼蒙皮天线机电耦合分析方法 |
CN105629725A (zh) * | 2014-10-31 | 2016-06-01 | 北京临近空间飞行器系统工程研究所 | 一种后缘舵滑翔飞行器的弹性运动建模方法 |
CN105843073A (zh) * | 2016-03-23 | 2016-08-10 | 北京航空航天大学 | 一种基于气动力不确定降阶的机翼结构气动弹性稳定性分析方法 |
CN106444885A (zh) * | 2016-11-09 | 2017-02-22 | 南京航空航天大学 | 一种颤振主动抑制控制器构成及其模拟方法 |
CN106773782A (zh) * | 2016-12-15 | 2017-05-31 | 中国航空工业集团公司西安飞机设计研究所 | 一种气动伺服弹性混合建模方法 |
US9694918B1 (en) * | 2016-05-26 | 2017-07-04 | Beihang University | Method for disturbance compensation based on sliding mode disturbance observer for spacecraft with large flexible appendage |
CN107976908A (zh) * | 2017-12-06 | 2018-05-01 | 中国科学院力学研究所 | 一种飞行器耦合动稳定性特征分析方法 |
CN108256264A (zh) * | 2018-02-08 | 2018-07-06 | 北京航空航天大学 | 一种基于地面频响试验的气动伺服弹性稳定性预测方法 |
-
2019
- 2019-03-20 CN CN201910212456.7A patent/CN110287505B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7437014B2 (en) * | 2004-04-21 | 2008-10-14 | The Boeing Company | Methods and systems for model reduction and system identification of dynamic systems with multiple inputs |
CN102169046A (zh) * | 2010-12-22 | 2011-08-31 | 北京航空航天大学 | 一种磁悬浮机电设备弹性模态在线测试系统 |
CN103077259A (zh) * | 2011-10-26 | 2013-05-01 | 上海机电工程研究所 | 高超声速导弹多场耦合动力学一体化仿真分析方法 |
CN104049533A (zh) * | 2013-03-11 | 2014-09-17 | 王岩 | 基于隐半马尔可夫模型的网络控制系统控制器设计方法 |
CN103400035A (zh) * | 2013-07-29 | 2013-11-20 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种高可信度快速预测飞行器滚转动导数的方法 |
CN103853157A (zh) * | 2014-03-19 | 2014-06-11 | 湖北蔚蓝国际航空学校有限公司 | 一种基于自适应滑模的飞行器姿态控制方法 |
CN104063533A (zh) * | 2014-04-04 | 2014-09-24 | 西安电子科技大学 | 一种机翼蒙皮天线机电耦合分析方法 |
CN104022742A (zh) * | 2014-06-19 | 2014-09-03 | 金陵科技学院 | 基于神经网络观测器的飞行器姿态鲁棒反演容错控制方法 |
CN105629725A (zh) * | 2014-10-31 | 2016-06-01 | 北京临近空间飞行器系统工程研究所 | 一种后缘舵滑翔飞行器的弹性运动建模方法 |
CN105843073A (zh) * | 2016-03-23 | 2016-08-10 | 北京航空航天大学 | 一种基于气动力不确定降阶的机翼结构气动弹性稳定性分析方法 |
US9694918B1 (en) * | 2016-05-26 | 2017-07-04 | Beihang University | Method for disturbance compensation based on sliding mode disturbance observer for spacecraft with large flexible appendage |
CN106444885A (zh) * | 2016-11-09 | 2017-02-22 | 南京航空航天大学 | 一种颤振主动抑制控制器构成及其模拟方法 |
CN106773782A (zh) * | 2016-12-15 | 2017-05-31 | 中国航空工业集团公司西安飞机设计研究所 | 一种气动伺服弹性混合建模方法 |
CN107976908A (zh) * | 2017-12-06 | 2018-05-01 | 中国科学院力学研究所 | 一种飞行器耦合动稳定性特征分析方法 |
CN108256264A (zh) * | 2018-02-08 | 2018-07-06 | 北京航空航天大学 | 一种基于地面频响试验的气动伺服弹性稳定性预测方法 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112528408A (zh) * | 2020-12-11 | 2021-03-19 | 中国直升机设计研究所 | 一种直升机旋翼与机身耦合稳定性建模方法 |
CN113111430A (zh) * | 2021-03-06 | 2021-07-13 | 北京航空航天大学 | 基于非线性气动力降阶的弹性飞机飞行动力学建模方法 |
CN113111430B (zh) * | 2021-03-06 | 2022-09-30 | 北京航空航天大学 | 基于非线性气动力降阶的弹性飞机飞行动力学建模方法 |
CN114912202A (zh) * | 2022-05-24 | 2022-08-16 | 大连理工大学 | 宽速域吸气式动力飞行器机体推进一体化耦合控制方法 |
CN116956782A (zh) * | 2023-09-20 | 2023-10-27 | 西北工业大学 | 非线性颤振分析方法 |
CN116956782B (zh) * | 2023-09-20 | 2023-12-01 | 西北工业大学 | 非线性颤振分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110287505B (zh) | 2020-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110287505A (zh) | 飞行器稳定性分析方法 | |
CN104443427B (zh) | 飞行器颤振预测系统及方法 | |
Biedron et al. | Recent enhancements to the FUN3D flow solver for moving-mesh applications | |
Cunningham et al. | Modern wing flutter analysis by computational fluid dynamics methods | |
CN109141476B (zh) | 一种动态变形下传递对准过程中角速度解耦合方法 | |
CN107976908B (zh) | 一种飞行器耦合动稳定性特征分析方法 | |
CN113128085B (zh) | 一种基于状态空间形式涡格法的阵风减缓分析方法 | |
CN111736468B (zh) | 一种信息融合下的飞行器抗干扰控制方法 | |
CN109840349A (zh) | 一种固定翼飞机阵风响应建模分析方法 | |
Pfifer et al. | System Identification of a small flexible aircraft-invited | |
CN103400035A (zh) | 一种高可信度快速预测飞行器滚转动导数的方法 | |
CN103310060A (zh) | 一种跨音速极限环颤振分析方法 | |
CN106444885A (zh) | 一种颤振主动抑制控制器构成及其模拟方法 | |
CN113221237A (zh) | 一种基于降阶建模的大迎角颤振分析方法 | |
Chen et al. | Limit cycle oscillation control for transonic aeroelastic systems based on support vector machine reduced order model | |
Schirrer et al. | A comprehensive robust control design and optimization methodology for complex flexible-structure systems | |
Waite et al. | Aeroservoelastic Control Law Development for the Integrated Adaptive Wing Technology Maturation Wind-Tunnel Test | |
CN108388741B (zh) | 飞行器颤振分析网格模型沃尔什建模方法 | |
CN108446461A (zh) | 飞行器颤振分析网格模型建模方法 | |
Loubimov et al. | Measuring Atmospheric Boundary Layer Profiles Using UAV Control Data | |
Schütte et al. | Experimental and numerical aspects of simulating unsteady flows around the X-31 configuration | |
Gang et al. | Active control law design for Flutter/LCO suppression based on reduced order model method | |
Zeng et al. | Ground vibration test identified structure model for flutter envelope prediction | |
CN108399297B (zh) | 飞行器颤振分析网格模型勒让德建模方法 | |
Sundresan et al. | Review of aeroelasticity testing technology |
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 |