超声速物面振动气动力的非线性效应研究

刘 成, 叶正寅

(西北工业大学 航空学院, 陕西 西安 710072)

0 引 言

本文对比了活塞理论、考虑激波膨胀波的活塞理论、Euler方程对气动力的计算结果，研究了在颤振计算过程中气动力非线性的主要来源，并对影响非线性效应的主要因素进行了分析。

1 考虑激波膨胀波的活塞理论

活塞理论(Piston Theory, PT)是由动量定理和等熵关系式推导得到的气动力表达式。在超声速、高超声速非定常气动力计算中，活塞理论具有良好的精度。其表达式如下：

(1)

(2)

式中，p为当地压强，p∞、a∞、v∞分别为来流压强、声速、速度，vn为非定常下洗速度，w为位移。

在超声速流场中，物面的振动变形会产生一系列的激波膨胀波，从而导致气动力计算出现非线性特征。从活塞理论的推导过程来看，并没有显式地考虑激波和膨胀波的存在，为此，本文建立了考虑激波膨胀波的活塞理论(Shock Expansion-Piston Theory, SE-PT)。

激波角β与气流偏转角θ的关系为：

(3)

式中Ma1为波前马赫数。求得β后再根据斜激波关系式可求得波后气动参数：

(4)

(5)

(6)

式中，Man,1、ρ1、p1为波前马赫数的垂直分量、波前密度、波前压强，Man,2、ρ2、p2为波后马赫数的垂直分量、波后密度、波后压强。

对于膨胀波来说，结合普朗特-迈耶函数v(Ma)和气流偏转角θ即可求得波后马赫数，公式可写为：

θ=v(Ma2)-v(Ma1)

(7)

式中，Ma1为波前马赫数，Ma2为波后马赫数。再由膨胀波关系式求出波后气动参数：

(8)

(9)

式中，ρ1、p1为波前密度、波前压强，ρ2、p2为波后密度、波后压强。

这样，在壁板的每一步振动中，都调用激波膨胀波关系式，得到当地的压强，并结合活塞理论计算气动力，即为考虑激波膨胀波的活塞理论(SE-PT)。

在颤振分析中，壁板模态m一般选取为如下形式：

(10)

式中，A表示广义位移，i表征模态阶数，2l为壁板长度。可以看到，模态振型为高阶光滑的连续函数，满足正交性，符合几何边界条件的要求。这种模态分解方法，在壁板颤振分析中取得极大成功。但对本文的研究内容来说，正弦形式的模态却不是最合适的，原因有两点：

1) 本文首先建立了考虑激波膨胀波的活塞理论，而激波膨胀波关系式不适应于求解光滑曲面的气动载荷；

2) 本文探究的是气动力的非线性效应，与壁板结构的真实模态形状无关。

因此，需要重新设计模态振型。新模态主要需考虑两点因素：便于激波膨胀波关系式的使用，满足几何边界条件。基于这样的考虑,我们将模态振型构造为折线形状。设壁板长度为2l，模态表达式如下：

(11)

采用式(11)表征的模态，可近似应用激波膨胀波关系式计算变形后的定常气动载荷，且该模态保留了壁板结构的主要几何特征，得到的气动力分析结果是具有普遍性意义的。

取计算状态为Ma=3，输入如图1所示的扫频位移信号，比较SE-PT理论与Euler方程的气动力系数F求解结果，如图2所示，可见二者符合良好，证明所采用的考虑激波膨胀波的活塞理论在气动力预测方面具有良好的精度。

图1 扫频位移信号Fig.1 Sweep signal

图2 扫频信号下的气动力响应Fig.2 Aerodynamic response of sweep signal

2 一阶活塞理论在壁板振动中的应用

计算状态取为Ma=3，研究对象为壁板，一阶模态为前述的折线形状，输入频率50 Hz的正弦位移信号，如图3所示。分别采用一阶活塞理论(PT)、考虑激波膨胀波的一阶活塞理论(SE-PT)、Euler方程计算气动力系数，结果如图4所示。可以清晰地看出，SE-PT理论与Euler方程的求解结果符合良好，气动力响应的非线性特征较为明显，而一阶活塞理论只能预测气动力的线性部分。将气动力响应做频谱分析，得到的幅频特性曲线如图5所示，从中可以看出，气动力响应的一阶谐波分量与二阶谐波分量的幅值比约为1:0.5，而一阶分量与三阶谐波分量的比值约为1:0.02，基本可以忽略。也就是说，气动力响应的非线性特征主要是出现了二阶谐波分量。

图3 50 Hz位移信号Fig.3 50 Hz displacement signal

图4 气动力响应对比图Fig.4 Comparison between aerodynamic responses

图5 气动力响应幅频特性曲线Fig.5 Amplitude-frequency curve of aerodynamic response

将一阶活塞理论中加入激波膨胀波关系式，得到的非定常气动力计算结果与Euler方程相同，因此得出，气动力的非线性特征是激波膨胀波带来的。

上述设计的壁板模态形状前后具有对称性，因此将壁板从中间位置分为两部分——前半部分和后半部分，同时也将壁板表面气动力分为前后两部分，分别作出两部分的气动力响应曲线以及幅频特性曲线，如图6、图7所示。可见，壁板前后两部分的气动力一阶、二阶谐波分量均相同，二阶分量仅占一阶分量的9%左右。而根据图5，整个壁板的气动力二阶谐波分量约占一阶分量的50%，非线性分量占比显著增加。

图6 前后表面气动力响应对比图Fig.6 Comparison between aerodynamic responses for former and latter parts

图7 前后表面气动力幅频特性曲线Fig.7 Amplitude-frequency curves of aerodynamic responses for former and latter parts

表1展示了气动力一阶、二阶谐波分量的相位大小。对一阶谐波分量来说，前后两部分的相位差较大，近似于反向，而前后表面的二阶分量几乎同相位。正是因为壁板前后部分的一阶分量反相、二阶分量同相，所以两个表面9%的非线性分量在求和后占比达到了50%，非线性效应显著增强。进一步探究一阶谐波分量反相的原因。活塞理论的一阶项可写为如下形式：

(12)

该式的第一项为变形分量，第二项为振动分量，由此将物面振动产生的气动力分解为变形导致的气动力分量和振动导致的气动力分量两部分，前者正比于壁板相对于来流的迎角，后者正比于壁板的横向振动速度。对于壁板的前后两部分来说，气流的迎角相反而横向振动速度相同，所以总的气动力中变形分量前后抵消，而振动分量则加强，同时在低速振动中，变形气动力分量占主导地位。因此前后部分的一阶气动力几乎是反相位的。

表1 前后表面气动力分量的相位Table 1 Phrases of aerodynamic harmonic components for former and latter part

3 一阶活塞理论在具有双侧来流的物面振动中的应用

取长度为l的无厚度板，绕前端点做俯仰运动，俯仰位移如图8所示，得到的气动力系数响应如图9所示。一阶活塞理论的求解结果与SE-PT理论、Euler方程的结果符合良好，说明此时的气动力响应中不存在高阶谐波分量。输入单一频率(50 Hz)的广义位移，得到的气动力响应如图10所示，频谱分析的幅频特性曲线如图11。可见，气动力响应只有一阶谐波分量，一阶活塞理论的计算结果是准确的。

图8 扫频位移信号Fig.8 Sweep signal

图9 扫频信号下气动力响应Fig.9 Aerodynamic response of sweep signal

图10 正弦位移信号Fig.10 Sinusoidal signal

图11 正弦信号下气动力响应Fig.11 Aerodynamic response of sinusoidal signal

对壁板这种只有单侧来流的物面振动来说，气动力的非线性特征明显。而在双侧来流的物面振动中，却没有非线性分量。下面对这一现象做进一步的研究。把双侧来流的物面分为上表面和下表面，分别计算它们的气动力响应。图12展示了气动力的时域响应曲线，图13展示了气动力的幅频特性曲线。

图12 上下表面气动力响应对比图Fig.12 Comparison between aerodynamic responses for upper and lower surfaces

图13 上下表面气动力幅频曲线Fig.13 Amplitude-frequency curves of aerodynamic responses for upper and lower surfaces

对上、下单个表面来说，气动力出现了明显的二阶谐波分量，其与一阶谐波分量的幅值比约为9%。表2展示了上、下表面各阶谐波分量的相位关系，可以看到，一阶分量相位相同，二阶分量相位相反，这就使得在计算总的气动力时，二阶谐波分量相互抵消，呈现出良好的线性特征。激波膨胀波带来了气动力的非线性效应，在物面振动过程中上表面出现激波时，下表面往往出现膨胀波，反之亦然。激波与膨胀波对气动力的影响相反，这是导致总的非线性效应相互抵消的物理原因。

表2 上、下表面气动力分量的相位Table 2 Phrases of aerodynamic harmonic components for upper and lower surfaces

4 高阶活塞理论

图14对比了50 Hz广义位移输入信号下不同阶活塞理论的气动力计算结果，可见，二阶活塞理论基本能预测出气动力的非线性特征。Cheng和Mei[15]研究了不同阶活塞理论计算壁板振动极限环幅值的区别，从计算结果来看，一阶活塞理论不够精确，二阶与三阶活塞理论的差别在0.5%以内。Fridmann[2]对比了微幅振动下不同阶活塞理论得到的壁板表面气动力系数，发现二阶活塞理论、三阶活塞理论、Euler方程的计算结果一致。

图15对比了不同阶考虑激波膨胀波的活塞理论气动力计算结果，由于激波膨胀波关系式的引入，一阶理论就已经能够很好地呈现出气动力的非线性特征。

图14 不同阶活塞理论的气动力响应对比Fig.14 Aerodynamic response comparison between different-order piston theories

图15 考虑激波膨胀波活塞理论的气动力响应对比Fig.15 Aerodynamic response comparison between piston theories considering shock and expansion waves

对比二阶活塞理论、考虑激波膨胀波的一阶活塞理论与Euler方程在不同位移信号下的气动力求解结果，如图16所示，三者符合良好，证明壁板颤振中气动力的非线性因素主要是由激波膨胀波带来的，而活塞理论的二阶项基本可以反映出这种非线性。

实际工程中的壁板颤振分析可能面临着很宽的结构频率范围，而活塞理论在不同的频率范围下呈现出不同的计算效果。为了澄清这方面的机制，分析了频率对气动力非线性的影响。

取三个不同频率(50 Hz、200 Hz、500 Hz)的位移信号，分别采用一阶、二阶活塞理论计算气动力系数。图17是气动力的时域响应，图18是对应的幅频特性曲线。三个频率下气动力出现了不同程度的非线性效应，二阶谐波分量与一阶谐波分量幅值比分别为0.503、0.042、0.172，可以看到，气动力的非线性程度与频率有很大关系。进一步做出二阶谐波与一阶谐波分量幅值比随频率变化的曲线图，如图19(a)，图中另外画出了幅值比为0.05的虚线。可以看出，当振动频率较高或者较低时，气动力的非线性效应变强，而在中间某段频率(200～300 Hz)内，非线性分量占比很低，此时气动力可以近似看做是线性的。

探究在中频段内气动力的非线性分量占比很低的原因。观察气动力的一阶、二阶谐波分量幅值随频率的变化(图19b)。由前面的分析可知，一阶谐波可分为变形导致的气动力分量和振动导致的气动力分量，前者与振动频率无关，后者与振动频率成线性增大关系，二者的综合作用，使得气动力的一阶谐波分量幅值随频率线性增大。图19(b)中的虚线展示了这一规律。实线代表的是二阶谐波分量幅值随频率的变化关系，可以看到是先减小后增大的。因此在中间频段二阶分量的占比降低，气动力的线性效应增强。

图16不同位移信号下气动力响应对比图
Fig.16Serodynamicresponsecomparisonfordifferentdisplacementsignals

(a) 50 Hz (b) 200 Hz (c) 500 Hz

(a) Amplitude ratio vs. frequency (b) Amplitude vs. frequency

进一步探究二阶谐波的幅值随频率先减小后增大的原因。下式为活塞理论的二阶分量：

图20 气动力二阶谐波变形分量和振动分量幅值随频率变化Fig.20 Deformation andvibration components of 2nd aerodynamic harmonics vs. frequency

5 结 论

本文选取超声速来流的计算状态，忽略壁面黏性效应，对比了活塞理论、考虑激波膨胀波的活塞理论、Euler方程的气动力计算结果，探究了壁板颤振中气动力非线性的主要来源，同时研究了具有双侧气流的物面振动气动力的线性性质，解释了振动频率影响气动力非线性效应的原因，得到的主要结论如下：

(1) 超声速物面振动过程中产生的激波膨胀波是气动力非线性的主要来源；

(2) 一阶活塞理论适应于高超声速翼面气动弹性分析，其原因是上下表面的非线性分量相互抵消，气动力的线性效应增强；而在只有单侧来流的壁板颤振分析中，模态变形量的对称性使得气动力的线性分量降低，非线性分量升高，此时需采用高阶活塞理论来模拟这种非线性效应；

(3) 活塞理论中的高阶项表征的是激波膨胀波带来的非线性效应，它随振动频率的增加先减弱后增强，其原因是非线性项中的变形气动力分量和振动气动力分量相位相反，前者幅值不随频率变化，后者随频率增加而增大，在中间某段频域内二者的幅值近似相等，非线性效用相互抵消，只有此时气动力才可以看作是线性的。