[发明专利]一种跨音速极限环颤振分析方法

权利要求书

1.一种跨音速极限环颤振分析方法，其特征在于：采用以下步骤：

步骤1：应用拉格朗日方程建立二元机翼颤振运动微分方程：

mh··+Sαα··+Khh=-12ρV2(2b)c1Sαh··+Iαα··+Kαα=12ρV2(2b)2cm]]>

其中b为半弦长；m为质量；ρ为空气密度；V为速度；c_l为升力系数；c_m为力矩系数；S_α为机翼对刚心的质量静矩，S_α＝mx_αb；I_α为机翼对刚心的质量惯性矩，K_h为线弹簧刚度，ω_h为沉浮自由度的“部分”频率；K_α为扭转弹簧刚度，ω_α为俯仰自由度的“部分”频率；二元机翼刚心位于弦线中点后ab处，a为系数；

步骤2：引入无量纲质量将步骤1中的二元机翼颤振运动微分方程转化矩阵形式：

[M]{ξ··}+ωα2[K]{ξ}=V2πμb2{f}]]>

其中{ξ}=hbα]]>为广义位移，为沉浮位移，α为俯仰位移；{f}=-cl2cm]]>为广义气动力；[M]=1xαxαrα2]]>为质量阵；[K]=(ωhωα)200rα2]]>为刚度阵；

步骤3：计算线性颤振边界：

取沉浮位移和俯仰位移为

h/b=0α=α0sin(ωt)]]>

其中α₀不大于0.5°；通过最小二乘法拟合得到由俯仰位移引起的线性广义气动力的形式为

为待定系数；则由俯仰位移引起的跨音速动态线性频域气动力系数为

取沉浮位移和俯仰位移为

h/b=(h/b)0sin(ωt)α=0]]>

其中(h/b)₀不大于0.05；通过最小二乘法拟合得到由沉浮位移引起的线性广义气动力的形式为

为待定系数；则由沉浮位移引起的跨音速动态线性频域气动力系数为

机翼的跨音速线性广义气动力由俯仰位移和沉浮位移所引起的线性广义气动力叠加得到

(-cl)(2cm)L=(-cl)L,h(2cm)L,h+(-cl)L,α(2cm)L,α=QlL,hQmL,h(h/b)0+QlL,αQmL,αα0=QlL,hQlL,αQmL,hQmL,α(h/b)0α0=[Q0](h/b)0α0]]>

其中频域气动力系数矩阵[Q₀]为马赫数及减缩频率的函数[Q₀(Ma,k)]；将机翼的跨音速线性广义气动力代入步骤2中二元机翼颤振运动微分方程的矩阵形式内，并将矩阵形式转化到频域内，得到频域颤振方程

-ω2[M]{ξ0}+ωα2[K]{ξ0}=V2πμb2[Q0(Ma,k)]{ξ0}]]>

对频域颤振方程通过频域颤振分析方法直接求解得到无量纲化的线性颤振速度无量纲化的线性颤振频率及线性颤振模态{uL}=(hbα)fL1T;]]>

步骤4：计算跨音速极限环特性：

步骤4.1：根据步骤3得到的线性颤振模态和给定的俯仰极限环幅值α_0,NL，得到相应的沉浮极限环幅值

步骤4.2：取沉浮位移和俯仰位移为

h/b=0α=α0,NLsin(ωt)]]>

通过最小二乘法拟合得到由俯仰位移引起的非线性广义气动力的形式为

为待定系数；相应的由俯仰位移引起的跨音速气动力描述函数为

取沉浮位移和俯仰位移为

h/b=(h/b)0,NLsin(ωt)α=0]]>

通过最小二乘法拟合得到由沉浮位移引起的非线性广义气动力的形式为

为待定系数；相应的由沉浮位移引起的跨音速气动力描述函数为

机翼的跨音速非线性广义气动力由俯仰位移和沉浮位移所引起的非线性广义气动力叠加得到

(-cl)(2cm)=(-cl)NL,h(2cm)NL,h+(-cl)NL,α(2cm)NL,α=QlNL,hQmNL,h(h/b)0,NL+QlNL,αQmNL,αα0,NL=QlNL,hQlNL,αQmNL,hQmNL,α(h/b)0,NLα0,NL=[QD](h/b)0,NLα0,NL]]>

其中频域气动力系数矩阵[Q_D]为马赫数、减缩频率及给定俯仰位移幅值的函数[Q_D(Ma,k,α_0,NL)]；将机翼的跨音速非线性广义气动力代入步骤2中二元机翼颤振运动微分方程的矩阵形式内，并将矩阵形式转化到频域内，得到跨音速气动力非线性颤振方程

-ω2[M]{ξ0,NL}+ωα2[K]{ξ0,NL}=V2πμb2[QD(Ma,k,α0,NL)]{ξ0,NL}]]>

对跨音速气动力非线性颤振方程通过频域颤振分析方法直接求解得到给定极限环幅值α0,NL(h/b)0,NL]]>下的无量纲化的颤振速度无量纲化的颤振频率

步骤5：重复步骤4，得到不同给定极限环幅值下的无量纲化颤振速度和无量纲化颤振频率，得到考虑气动力非线性的跨音速极限环颤振特性。