一种带控制面的三自由度机翼的颤振边界预测方法技术

发明专利技术公开了一种带控制面的三自由度机翼的颤振边界预测方法，首先通过模态参数识别方法识别机翼的沉浮、俯仰与控制面偏转三个主要模态的频率与阻尼比，然后采用这些模态参数构造特征方程式并求解其特征系数，由特征系数计算一个稳定性判据F3；在每个不同风速下根据响应信号计算相应的F3，采用直线拟合外推颤振边界。该发明专利技术适用于带控制面的三自由度机翼的颤振边界预测，能够在较低风速下预测颤振边界，有助于提高试验的准确性与安全性。

【技术实现步骤摘要】
一种带控制面的三自由度机翼的颤振边界预测方法
本专利技术涉颤振边界预测
，具体涉及一种带控制面的三自由度机翼的颤振边界预测方法。
技术介绍
颤振试验是飞行器设计工作中必需的一个重要环节，通过试验得到的颤振边界能够确定飞行器的飞行包线，进而评判飞行器性能与保障飞行安全。颤振试验中常用的颤振边界预测方法是速度-阻尼比外推法，在不同风速下采集振动响应信号识别不同模态的频率与阻尼比，拟合危险模态的阻尼比相对于风速的变化曲线，外推曲线得到阻尼比降为零风速即为颤振速度。然而阻尼比与风速并非线性的关系，通常情况下，危险模态的阻尼比只在靠近颤振边界的风速时才有下降趋势，这将不利于颤振边界的提前预测，特别地，对于突发型颤振，速度-阻尼比外推法可能会失去预测颤振边界的能力。此外，因为受到噪声的影响，阻尼比的识别不可避免地将具有一定的误差，而速度-阻尼法仅以阻尼比作为外推依据，这种误差将在外推过程中放大，很容易造成外推失真。ZimmermanandWeissenburger提出了一种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，更新三本文档来自技高网...

【技术保护点】
一种带控制面的三自由度机翼的颤振边界预测方法，其特征在于，该方法包括以下步骤：1)在固定风速值下，通过模态参数识别方法识别机翼的沉浮、俯仰与控制面偏转三个模态的频率ωi参数与阻尼比gi参数；2)采用步骤1)中获得的三个模态的参数构造特征方程式s6+A5s5+A4s4+A3s3+A2s2+A1s+A0＝0，特征方程的特征系数A的计算公式如下：A5=-2Σi=13giA4=Σi=13(gi2+ωi2+4Πj=1:3j≠igj)A3=-2Σi=13gi(Σj=1:3j≠i(gj2+ωj2))+8Πk=13gkA2=Σi=13(4(gi2+ωi2)Πj=1:3j≠igj+Πj=1:3j≠i(gj2+ωj2))A1=-2Σi=13(giΠj=1:3j≠i(gj2+ωj2))A0=Πi=13(gi2+ωi2)]]>其中：A0、A1、A2、A3、A4、A5分别是特征方程的系数；i、j、k分别是叠加运算的系数，具体数值为叠加运算符上的数值；3)根据劳斯稳定判据，由以上特征系数A计算劳斯表前两列系数，其计算公式如下：P11=A5P21=A4-A3/A5,P22=A2-A1/A5P31=A3-A5P22/P21,P32=A1-A0A5/P21P41=P22-P21P32/P31,P42=A0P51=P32-A0P31/P41]]>其中，P11为劳斯列表第一行第一列；P21为劳斯列表第二行第一列；P22为劳斯列表第二行第二列；P31为劳斯列表第三行第一列；P32为劳斯列表第三行第二列；P41为劳斯列表第四行第一列；P42为劳斯列表第四行第二列；P51为劳斯列表第五行第一列；4)由劳斯表系数计算出稳定性判据F3：F3=P51P41P31P212/P114]]>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为临界点。...

【技术特征摘要】
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)根据劳斯稳定判据...

【专利技术属性】
技术研发人员：周丽，李扬，杨秉才，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏;32