CN104615863A - 一种带控制面的三自由度机翼的颤振边界预测方法 - Google Patents
一种带控制面的三自由度机翼的颤振边界预测方法 Download PDFInfo
- Publication number
- CN104615863A CN104615863A CN201510019577.1A CN201510019577A CN104615863A CN 104615863 A CN104615863 A CN 104615863A CN 201510019577 A CN201510019577 A CN 201510019577A CN 104615863 A CN104615863 A CN 104615863A
- Authority
- CN
- China
- Prior art keywords
- flutter
- wind speed
- value
- louth
- list
- 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
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
发明公开了一种带控制面的三自由度机翼的颤振边界预测方法,首先通过模态参数识别方法识别机翼的沉浮、俯仰与控制面偏转三个主要模态的频率与阻尼比,然后采用这些模态参数构造特征方程式并求解其特征系数,由特征系数计算一个稳定性判据F3;在每个不同风速下根据响应信号计算相应的F3,采用直线拟合外推颤振边界。该发明适用于带控制面的三自由度机翼的颤振边界预测,能够在较低风速下预测颤振边界,有助于提高试验的准确性与安全性。
Description
技术领域
本发明涉颤振边界预测技术领域,具体涉及一种带控制面的三自由度机翼的颤振边界预测方法。
背景技术
颤振试验是飞行器设计工作中必需的一个重要环节,通过试验得到的颤振边界能够确定飞行器的飞行包线,进而评判飞行器性能与保障飞行安全。
颤振试验中常用的颤振边界预测方法是速度-阻尼比外推法,在不同风速下采集振动响应信号识别不同模态的频率与阻尼比,拟合危险模态的阻尼比相对于风速的变化曲线,外推曲线得到阻尼比降为零风速即为颤振速度。
然而阻尼比与风速并非线性的关系,通常情况下,危险模态的阻尼比只在靠近颤振边界的风速时才有下降趋势,这将不利于颤振边界的提前预测,特别地,对于突发型颤振,速度-阻尼比外推法可能会失去预测颤振边界的能力。
此外,因为受到噪声的影响,阻尼比的识别不可避免地将具有一定的误差,而速度-阻尼法仅以阻尼比作为外推依据,这种误差将在外推过程中放大,很容易造成外推失真。
Zimmerman and Weissenburger提出了一种Z-W方法,可以进行典型机翼弯扭耦合的颤振边界预测,然而该方法只适用于二元机翼,对带控制面的三自由度机翼并不适用。
发明内容
发明目的:为了解决现有技术中,速度-阻尼法无法较早预测三自由度颤振边界的问题,本发明通过识别模态参数,并通过模态参数获得稳定判据,对稳定判据进行直线拟合外推,获得误差较小的颤振预测边界。
技术方案:一种带控制面的三自由度机翼的颤振边界预测方法,其特征在于,该方法包括以下步骤:
步骤1)在固定风速值下,通过模态参数识别方法识别机翼的沉浮、俯仰与控制面偏转三个模态的频率ωi参数与阻尼比gi参数;
步骤2)采用步骤1)中获得的三个模态的参数构造特征方程式s6+A5s5+A4s4+A3s3+A2s2+A1s+A0=0,特征方程的特征系数A的计算公式如下:
其中:
A0、A1、A2、A3、A4、A5分别是特征方程的系数;
i、j、k分别是叠加运算的系数,具体数值为叠加运算符上的数值;
步骤3)根据劳斯稳定判据,由以上特征系数A计算劳斯表前两列系数,其计算公式如下:
其中,
P11为劳斯列表第一行第一列;
P21为劳斯列表第二行第一列;P22为劳斯列表第二行第二列;
P31为劳斯列表第三行第一列;P32为劳斯列表第三行第二列;
P41为劳斯列表第四行第一列;P42为劳斯列表第四行第二列;
P51为劳斯列表第五行第一列;
步骤4)由劳斯表系数计算出稳定性判据F3:
步骤5)给定判定值N,0<N<1;
选取两个不同数值的风速,数值分别为X、Y,选取测量风速值为Z;
步骤6)利用步骤1)到步骤4)分别计算风速值X、Y、Z的稳定性判据;根据步骤5)中三个风速值和其分别对应的稳定性判据值,绘制稳定性判据值相对于风速的变化曲线,进行直线拟合并外推出F3降为0时的风速S,设定该风速S作为此刻测量风速下的预测颤振边界速度,即预测颤振边界速度为S;
步骤7)设定替换风速T,T值大于X、Y、Z中的最大值;用替换风速T替换更新X、Y、Z中的最小值,重复步骤5)-7),直到预测颤振速度S与替换风速T的差值小于判定值N;此时认定预测颤振边界速度S为临界点。
进一步的,判定值N=0.5,选择可以根据判定值控制误差,判定值越小,误差越小。判定值可以根据需要进行调节,这里的N=0.5为优选判定值。
有益效果:本发明基于带控制面的三自由度机翼颤振边界预测方法能够在风速较低的时候给出较接近的颤振边界,改善传统的颤振判断方法,有助于提高试验的准确性与安全性。
附图说明
图1是风速为50m/s的颤振边界预测结果;
图2是风速为80m/s的颤振边界预测结果;
图3是风速为110m/s的颤振边界预测结果;
图4是风速为130m/s的颤振边界预测结果;
图5是临近边界的逐步直线拟合结果;
图6是算例的V-f图与V-g图;
具体实施方式
下面结合附图对本发明做更进一步的解释。一种带控制面的三自由度机翼的颤振边界预测方法,其特征在于,该方法包括以下步骤:
步骤1)在固定风速值下,通过模态参数识别方法识别机翼的沉浮、俯仰与控制面偏转三个模态的频率ωi参数与阻尼比gi参数;
步骤2)采用步骤1)中获得的三个模态的参数构造特征方程式s6+A5s5+A4s4+A3s3+A2s2+A1s+A0=0,特征方程的特征系数A的计算公式如下:
其中:
A0、A1、A2、A3、A4、A5分别是特征方程的系数;
i、j、k分别是叠加运算的系数,具体数值为叠加运算符上的数值;
步骤3)根据劳斯稳定判据,由以上特征系数A计算劳斯表前两列系数,其计算公式如下:
其中,
P11为劳斯列表第一行第一列;
P21为劳斯列表第二行第一列;P22为劳斯列表第二行第二列;
P31为劳斯列表第三行第一列;P32为劳斯列表第三行第二列;
P41为劳斯列表第四行第一列;P42为劳斯列表第四行第二列;
P51为劳斯列表第五行第一列;
步骤4)由劳斯表系数计算出稳定性判据F3:
步骤5)给定判定值N,0<N<1;
选取两个不同数值的风速,数值分别为X、Y,选取测量风速值为Z;
步骤6)利用步骤1)到步骤4)分别计算风速值X、Y、Z的稳定性判据;根据步骤5)中三个风速值和其分别对应的稳定性判据值,绘制稳定性判据值相对于风速的变化曲线,进行直线拟合并外推出F3降为0时的风速S,设定该风速S作为此刻测量风速下的预测颤振边界速度,即预测颤振边界速度为S;
步骤7)设定替换风速T,T值大于X、Y、Z中的最大值;用替换风速T替换更新X、Y、Z中的最小值,重复步骤5)-7),直到预测颤振速度S与替换风速T的差值小于判定值N;此时认定预测颤振边界速度S为临界点。
进一步的,判定值N=0.5,选择可以根据判定值控制误差,判定值越小,误差越小。判定值可以根据需要进行调节,这里的N=0.5为优选判定值。
实施例:
(1)测得风速30m/s的三个模态频率与阻尼比分别为[46.75Hz,9.14%]、[104.04Hz,4.82%]、[434.10Hz,0.61%],根据前面所提出的方法计算得到F3=9.20E30;继续测得风速40m/s与50m/s的模态参数,计算得到判据F3分别为8.65E30与8.14E30,对F3与风速进行三点直线拟合,如图1所示,外推F3降为零时的风速为185.9m/s;
(2)根据上一步外推的颤振风速,选取适当的风速加载,替换上一步中最小的风速30m/s,本次选择风速为80m/s,计算得到判据F3=5.71E30,更新三点直线拟合所取数据为最新的三点,即40m/s、50m/s与80m/s的计算结果,对F3与风速进行三点直线拟合,如图2所示,外推F3降为零时的风速为157.1m/s;
(3)根据上一步外推的颤振风速,选取适当的风速加载,本次选择风速为110m/s,替换上一步中最小的风速40m/s,计算得到判据F3=3.12E30,更新三点直线拟合所取数据为最新的三点,即50m/s、80m/s与110m/s的计算结果,对F3与风速进行三点直线拟合,如图3所示,外推F3降为零时的风速为148.89m/s;
(4)根据上一步外推的颤振风速,选取适当的风速加载,本次选择风速为130m/s,替换上一步中最小的风速50m/s,计算得到判据F3=1.44E30,更新三点直线拟合所取数据为最新的三点,即80m/s、110m/s与130m/s的计算结果,对F3与风速进行三点直线拟合,如图4所示,外推F3降为零时的风速为146.70m/s;
(5)本例中判定值取为0.5m/s,逐步更新风速并外推新的预测风速,如图5所示,当风速达到171m/s时,预测的颤振边界为171.4m/s,二者之差小于判定值,认定颤振边界为171.4m/s,而实际颤振速度为172.1m/s,预测结果十分接近。改变判定值的大小可以调整预测结果的精度。
本例的V-g(风速-阻尼比)、V-f(风速-频率)图见图6,如果采用速度-阻尼方法进行颤振边界预测,危险模态的阻尼比只在140m/s之后才有下降趋势,因此在该风速之前无法进行预测,不利于试验安全。
本发明基于带控制面的三自由度机翼颤振边界预测方法能够在风速较低的时候给出较接近的颤振边界,采用逐步直线拟合能够准确的外推颤振边界,改善传统的颤振判断方法,有助于提高试验的准确性与安全性。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (2)
1.一种带控制面的三自由度机翼的颤振边界预测方法,其特征在于,该方法包括以下步骤:
1)在固定风速值下,通过模态参数识别方法识别机翼的沉浮、俯仰与控制面偏转三个模态的频率ωi参数与阻尼比gi参数;
2)采用步骤1)中获得的三个模态的参数构造特征方程式s6+A5s5+A4s4+A3s3+A2s2+A1s+A0=0,特征方程的特征系数A的计算公式如下:
其中:
A0、A1、A2、A3、A4、A5分别是特征方程的系数;
i、j、k分别是叠加运算的系数,具体数值为叠加运算符上的数值;
3)根据劳斯稳定判据,由以上特征系数A计算劳斯表前两列系数,其计算公式如下:
其中,
P11为劳斯列表第一行第一列;
P21为劳斯列表第二行第一列;P22为劳斯列表第二行第二列;
P31为劳斯列表第三行第一列;P32为劳斯列表第三行第二列;
P41为劳斯列表第四行第一列;P42为劳斯列表第四行第二列;
P51为劳斯列表第五行第一列;
4)由劳斯表系数计算出稳定性判据F3:
5)给定判定值N,0<N<1;
选取两个不同数值的风速,数值分别为X、Y,选取测量风速值为Z;
6)利用步骤1)到步骤4)分别计算风速值X、Y、Z的稳定性判据;根据步骤5)中三个风速值和其分别对应的稳定性判据值,绘制稳定性判据值相对于风速的变化曲线,进行直线拟合并外推出F3降为0时的风速S,设定该风速S作为此刻测量风速下的预测颤振边界速度,即预测颤振边界速度为S;
7)设定替换风速T,T值大于X、Y、Z中的最大值;用替换风速T替换更新X、Y、Z中的最小值,重复步骤5)-7),直到预测颤振速度S与替换风速T的差值小于判定值N;此时认定预测颤振边界速度S为临界点。
2.如权利要求1所述的一种带控制面的三自由度机翼的颤振边界预测方法,其特征在于,判定值N=0.5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510019577.1A CN104615863B (zh) | 2015-01-14 | 2015-01-14 | 一种带控制面的三自由度机翼的颤振边界预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510019577.1A CN104615863B (zh) | 2015-01-14 | 2015-01-14 | 一种带控制面的三自由度机翼的颤振边界预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104615863A true CN104615863A (zh) | 2015-05-13 |
CN104615863B CN104615863B (zh) | 2017-12-26 |
Family
ID=53150304
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510019577.1A Active CN104615863B (zh) | 2015-01-14 | 2015-01-14 | 一种带控制面的三自由度机翼的颤振边界预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104615863B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106599491A (zh) * | 2016-12-19 | 2017-04-26 | 厦门大学 | 一种基于qmu的颤振裕度评估方法 |
CN106777696A (zh) * | 2016-12-19 | 2017-05-31 | 厦门大学 | 基于qmu的颤振设计方法 |
CN112329274A (zh) * | 2020-12-29 | 2021-02-05 | 星河动力(北京)空间科技有限公司 | 火箭模态参数的确定方法、装置、设备及存储介质 |
CN116611175A (zh) * | 2023-07-18 | 2023-08-18 | 北京航空航天大学 | 一种大展弦比飞机体自由度颤振的预测方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103617323A (zh) * | 2013-11-27 | 2014-03-05 | 中国飞行试验研究院 | 一种基于飞行试验的颤振数据处理架构 |
-
2015
- 2015-01-14 CN CN201510019577.1A patent/CN104615863B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103617323A (zh) * | 2013-11-27 | 2014-03-05 | 中国飞行试验研究院 | 一种基于飞行试验的颤振数据处理架构 |
Non-Patent Citations (2)
Title |
---|
HAIXIN CHEN ET AL.: "Navier-Stokes simulations for transport aircraft wing body high-lift configurations", 《JOURNAL OF AIRCRAFT》 * |
S. J. PRICE ET AL.: "Evaluation and Extension of the Flutter-Margin Method for Flight Flutter Prediction", 《JOURNAL OF AIRCRAFT》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106599491A (zh) * | 2016-12-19 | 2017-04-26 | 厦门大学 | 一种基于qmu的颤振裕度评估方法 |
CN106777696A (zh) * | 2016-12-19 | 2017-05-31 | 厦门大学 | 基于qmu的颤振设计方法 |
CN106777696B (zh) * | 2016-12-19 | 2019-05-17 | 厦门大学 | 基于qmu的颤振设计方法 |
CN106599491B (zh) * | 2016-12-19 | 2020-05-15 | 厦门大学 | 一种基于qmu的颤振裕度评估方法 |
CN112329274A (zh) * | 2020-12-29 | 2021-02-05 | 星河动力(北京)空间科技有限公司 | 火箭模态参数的确定方法、装置、设备及存储介质 |
CN112329274B (zh) * | 2020-12-29 | 2021-03-23 | 星河动力(北京)空间科技有限公司 | 火箭模态参数的确定方法、装置、设备及存储介质 |
CN116611175A (zh) * | 2023-07-18 | 2023-08-18 | 北京航空航天大学 | 一种大展弦比飞机体自由度颤振的预测方法 |
CN116611175B (zh) * | 2023-07-18 | 2023-09-12 | 北京航空航天大学 | 一种大展弦比飞机体自由度颤振的预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104615863B (zh) | 2017-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104615863A (zh) | 一种带控制面的三自由度机翼的颤振边界预测方法 | |
US10690117B2 (en) | Method of controlling a wind turbine and related system | |
CN109190283A (zh) | 一种考虑高湍流自由来流效应的风力机翼型气动优化方法 | |
EP2637046B1 (en) | Method for predicting wind conditions in a wind farm | |
CN111079228B (zh) | 一种基于流场预测的气动外形优化方法 | |
US9821903B2 (en) | Closed loop control of aircraft control surfaces | |
CN102561118B (zh) | 一种基于轨迹偏差的高铁无砟轨道快速精调方法 | |
CN112362291B (zh) | 一种飞翼布局飞机纵向气动力系数的雷诺数效应修正方法 | |
CN105069476A (zh) | 基于两阶段集成学习的风电异常数据识别方法 | |
CN104881585A (zh) | 一种两自由度机翼的颤振边界预测方法 | |
CN101859146A (zh) | 一种基于预测滤波和经验模态分解的卫星故障预测方法 | |
CN103940608A (zh) | 一种提高风电机组齿轮箱故障等级判断精度的方法 | |
CN114912301B (zh) | 一种低速风洞全机模型测力试验数据处理与修正系统 | |
CN112052522B (zh) | 一种基于疲劳强度的船体结构优化简化计算方法 | |
CN103530511A (zh) | 紊流激励条件下的风洞颤振试验颤振边界预测方法 | |
CN112763180B (zh) | 声学风洞内高速列车受电弓模型气动噪声快速预测方法 | |
Wu et al. | Effects of lateral wind gusts on vertical axis wind turbines | |
CN113642068A (zh) | 一种基于决策树和循环神经网络的桥梁涡激振动幅值预测方法 | |
CN109540459B (zh) | 一种气动特性数值计算结果修正方法 | |
CN105134482A (zh) | 大型智能风机叶片系统灰色组合建模与优化振动控制的方法 | |
CN104713730B (zh) | 一种根据振动信号确定飞机发动机退化率的方法 | |
CN111693726B (zh) | 基于邻域粗糙集的通风系统故障诊断风速传感器布置方法 | |
CN105093932B (zh) | 一种确定lpv变增益控制器的鲁棒性的方法 | |
CN105508147A (zh) | 风电叶片单点疲劳加载试验弯矩匹配方法 | |
CN104627388A (zh) | 一种再入飞行器的自适应弹道预测方法 |
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 |