曲线纤维变刚度壁板的亚音速气弹稳定性研究

段静波 ，徐步青

（石家庄铁道大学工程力学系，河北 石家庄 050043）

为追求列车设计的高速化与轻量化，新型复合材料越来越广泛地应用于高速轨道交通领域，复合材料曲线纤维层合板作为一种近年新出现的复合材料，与传统的直线纤维铺层相比，由于纤维方向随位置的变化而变化，曲线纤维铺层呈现出变刚度的特点.Gűrdal等[1]研究表明，刚度在空间上非均匀分布的变刚度壁板能大幅改善结构承载能力，适应更为复杂的力学环境.将曲线纤维变刚度复合材料层合板应用于高速列车蒙皮、裙板等壁板结构，有利于新一代复合材料高铁列车高速化与轻量化设计目标的顺利实现.

然而，伴随着高速列车行进速度的不断提高，车体蒙皮、裙板等壁板类结构所面临的力学环境越来越恶劣，设计不合理可能会导致壁板发生剧烈的气动弹性振动，严重影响结构的使用寿命甚至列车的运行安全.因此，亚音速曲线纤维变刚度壁板气动弹性稳定性问题将是车体新型复合材料应用中的一个关键问题.

曲线纤维复合材料最早由Hyper等[2]提出，Gűrdal等[1,3]则明确提出了变刚度复合材料的概念.随后，学者们对于曲线纤维变刚度壁板力学问题进行广泛研究：Hamed等[4]对曲线纤维变刚度壁板进行大变形和应力分析；马洪涛[5]开展了曲线纤维变刚度壁板的力学性能研究；聂国隽等[6]建立了曲线纤维增强复合材料层合板自由振动问题的基本方程，得到了层合板的自振频率及相应的振型；马成[7]通过对纤维曲线铺放层合板阻尼比的研究，探究纤维曲线铺放层合板的减振效果；Groh等[8]研究了横向剪切效应对变刚度层合板屈曲性能的影响；Hao等[9]研究了复合材料变刚度板的屈曲及优化问题，提出了复合材料变刚度结构的精确建模、分析与结构优化一体化设计框架；孙士平等[10]研究了复合载荷作用下变刚度复合材料回转壳的屈曲优化问题.

迄今为止，有关曲线纤维变刚度壁板气动弹性问题的研究较少，但已经引起了国内外学者的关注：Vahid等[11]研究了高速流场中的曲线纤维变刚度层合壁板颤振问题；欧阳小穗和刘毅[12]研究了高速流场中的曲线纤维变刚度层合壁板非线性颤振响应；Touraj等[13]将曲线纤维变刚度壁板应用于飞机机翼承力结构，研究了机翼-引擎系统的气动弹性特性.

基于此，本文以新一代新型复合材料高速列车壁板结构为背景，基于目前列车设计时速范围，研究曲线纤维变刚度壁板的亚音速气动弹性稳定性问题，考虑壁板可能出现的黏弹性效应，分析复合材料关键参数对复合变刚度壁板气动弹性稳定性的影响规律.

1 曲线纤维路径表征

图1为曲线纤维复合材料层合板示意，以层合板中第k层单层变刚度铺层为研究对象(图1（b）)，图中：zk和zk+1分别为该铺层上、下表面在坐标系xOz中的坐标值；hk为该铺层厚度；V∞为来流速度；L、H、h分别为矩形变刚度复合材料层合板的长度、宽度、厚度；在局部坐标系x′Oy′下，T0为纤维在板中心处与x′轴方向的夹角；T1为纤维在边界(x′= ±L/2)处与x′轴方向的夹角；ϑ为x′轴与x轴的夹角，表示该铺层的方向角.该铺层曲线纤维方向角用 ＜T0｜T1＞表示.假设纤维方向从板的中心位置开始呈线性变化，则该铺层在局部坐标系下任意位置处纤维方向角可表示为

图1 曲线纤维层合板Fig.1 Curvilinear fibers composite panels

2 复合变刚度壁板颤振方程

2.1 壁板变形几何关系

结构方面，考虑厚度方向上的剪切变形，忽略厚度方向上法向变形，采用Mindlin厚板理论描述曲线纤维变刚度层合壁板，其应变-位移关系为

式中： ε 为线应变； ε0为壁板中面面内位移产生的应变向量，如式（3）； γ 为横向剪切应变向量，如式（3）；κ为弯曲时壁板的曲率向量，如式（3）；z为壁板厚度方向的坐标.

式中：u、 υ 分别为中面上的点沿x轴和y轴方向的位移；w为z轴方向的位移； θx、 θy分别为中面绕x轴和y轴方向的转角.

2.2 变刚度壁板弹性/黏弹性本构关系

对于正交各向异性铺层材料，采用kelvin-Voigt黏弹性本构模型，则其应力-应变关系[14]可写为

式中： η 为材料黏弹性阻尼系数，当 η =0 时，铺层材料就退化为弹性材料； σ 、τ、γ、˙ 分别为正应力、剪应力、切应变、线应变率；Qm、Qs分别为面内刚度和剪切刚度，其元素与弹性本构关系一致，在此不再赘述，如（5）.

由于曲线纤维角度影响，单层板内本构关系随壁板平面是变化的，当多个单层板形成层合板后，则曲线纤维角度和各层铺层角共同影响壁板本构关系，第k层本构关系为

式 中 ： σk、τk分 别 为 该 铺 层 的 正 应 力 和 切 应 力 ；εk、γk分别为该铺层的线应变和切应变；Qmk、Qsk分别为偏轴面内刚度矩阵和剪切刚度矩阵，=，其中，Tmk、Tsk如式（7）.

式中：Tk为曲线纤维方位角.

根据复合材料层合理论，对于n层的复合变刚度壁板本构关系可写为

式中：N、M、Fs分别为复合变刚度壁板的膜力、弯矩、横向剪力，矩阵A、B、D、As分别为

2.3 亚音速气动力模型

气动力方面，由于高速列车马赫数较小时，气体的可压缩性可以忽略，因此，采用Dowell基于线性无黏势流理论提出的亚音速气动力模型[15]，如式（9）.

式中：A(x) 为壁板内气动力相互影响效应，如式（10）；D1(x)、D2(x)为壁板边界条件对气动力的影响效应，如式（10）， ρ∞为空气密度.

2.4 气动弹性有限元方程

基于复合变刚度壁板黏弹性本构关系、几何关系以及气动力模型，可以写出其虚动能 δT、虚应变能 δU、外力虚功 δW为

采用板单元对复合变刚度壁板进行网格划分，根据Hamilton原理，通过变分可得到壁板的有限元颤振方程为

3 算法验证与结果分析

3.1 正确性验证

首先，对本文方法及程序进行正确性验证.采用文献[12]中曲线纤维变刚度复合材料层合壁板算例，通过退化气动力相关矩阵，选取铺层方式为[0/90＜0|45°＞/＜0|-45°＞]s的变刚度复合壁板进行了固有频率的网格收敛性分析，具体结果见表1.分析时壁板采用四边固支约束边界条件.为方便对比，固有频率按文献[12]进行无量纲处理.

表1 变刚度复合壁板固有频率的网格收敛性Tab.1 Grid convergence of natural frequencies of composite panels with variable stiffness

从表1中可以看到：随着壁板网格规模从5 ×5增加至30 × 30，变刚度复合壁板前六阶固有频率逐渐收敛，且收敛后的各阶固有频率与文献[12]的结果最大误差不超过3.56%.

图2给出了6种铺层方式下变刚度复合壁板的固有频率，6种铺层方式分别为铺层A1[0/90＜0|45°＞/＜0|-45°＞]s、铺层 A2[0/90＜15°|45°＞/＜-15°|-45°＞]s、铺层 A3[0/90＜30°|45°＞/＜-30°|-45°＞]s、 铺层 A4[0/90＜45°|40°＞/＜-45°|-60°＞]s、 铺层 A5 [0/90＜45°|75° ＞/＜-45°|-75°＞]s、 铺层A6[0/90＜90°|45°＞/＜-45°|-90°＞]s.壁板采用30 × 30的网格划分，四边固支约束边界条件.从图中对比可以看出：各种铺层方式下，本文前三阶固有频率结果与文献[12]的结果吻合得都比较好.因此，本文后续壁板稳定性分析均采用30 × 30的网格划分.

图2 曲线纤维壁板固有频率结果对比Fig.2 Comparison of natural frequencies of curved fiber composite laminates

其次，采用文献[16]中的两边固支铝制壁板，将本文程序退化为各向同性材料壁板进行气动弹性稳定性计算，并分析验证，结果如表2.从表2中可以看到，本文方法计算结果与文献[16]解吻合很好.

表2 两端固支壁板发散稳定性结果对比Tab.2 Divergence stability of the plate with bilateral fixation

3.2 复合变刚度壁板发散稳定性分析

1） 失稳速度计算

采用复模态理论求解复合变刚度弹性、黏弹性壁板的气动弹性发散失稳速度.图3给出了铺层方式为[0/90/＜0|45＞/＜0|-45＞]s的弹性/黏弹性变刚度壁板在四边固支、四边简支条件下，前两阶频率随来流速度变化情况.黏弹性材料的黏性阻尼系数η=0.002，下文中不特别注明，均取此值.

从图3中可以看到：无论变刚度壁板是弹性材料还是黏弹性材料，其一阶固有频率随着来流速度增大均是减小的，二阶固有频率随着来流速度先减小后增大；变刚度壁板材料的黏弹性效应会导致壁板二阶固有频率明显降低；当一阶固有频率减小至0时，壁板就出现发散失稳，据此可判断，四边固支的变刚度弹性、黏弹性壁板发散失稳速度分别为108 m/s和91 m/s，四边简支弹性、黏弹性变刚度壁板发散失稳速度分别为72 m/s和66 m/s.由此可知，与弹性变刚度壁板相比，壁板材料的黏弹性效应会导致其失稳速度降低.从图中失稳前壁板固有频率值变化剧烈程度可以看到，黏弹性壁板发散失稳出现得更突兀.

图3 黏弹性和弹性变刚度壁板频率随气流速度变化Fig.3 Change of the two frequencies of viscoelastic and elastic variable-stiffness panels with air velocity

2） 曲线铺设和直线铺设的比较

图4为弹性、黏弹性壁板分别采用直线和曲线纤维铺设情况下，四边固支纤维复合弹性壁板的发散稳定特性.直线铺设与曲线铺设的铺层角考虑了7种组合方式，分别为铺层 B1[0/15/45/60]、铺层B2[0/30/45/75]、铺层 B3[0/15/30/60]、铺层 B4[30/45/60/75]、铺层 B5[45/15/30/60]、铺层 B6[30/15/45/60]、铺层B7[45/30/60/75].其中，曲线纤维铺设时，纤维方向角T0=0°，T1=45°.从图中可以看到：无论是弹性壁板还是黏弹性壁板，在纤维曲线铺设下的失稳速度相比直线纤维铺设都有提高.采用铺层B7时，曲线纤维壁板失稳速度相比直线纤维壁板可提高将近50%.

图4 直线铺设与曲线铺设对变刚度壁板发散临界速度的影响Fig.4 Comparison of variable-stiffness panel instability with straight and curved fiber laying

3） 铺层主方向角的影响

图5为铺层主方向对曲线纤维变刚度壁板发散临界速度的影响.7种不同的铺层方式为铺层C1[0/0]s、 铺 层 C2[0/15]s、 铺 层 C3[0/30]s、 铺 层C4[0/45]s、 铺 层 C5[0/60]s、 铺 层 C6[0/75]s、 铺 层C7[0/90]s等，壁板采用四边固支约束.从图中可以看出，随着铺层主方向角从 0°～90° 增大，无论是变刚度弹性壁板还是变刚度黏弹性壁板，壁板发散失稳速度呈增大的趋势.由此可见，通过增大铺层主方向角可以提高复合材料变刚度壁板发散速度.

图5 铺层主方向对曲线纤维变刚度壁板发散临界速度影响Fig.5 Influence of ply direction on divergence critical speed of variable-stiffness panels with curved fibers

4） 曲线纤维方向角的影响

进一步研究曲线纤维方向角对弹性、黏弹性变刚度复合壁板发散稳定性的影响.壁板为四边固支约束，其铺层形式为[0/90/＜T0|T1＞/＜-T0|-T1＞]s.考虑T0和T1单独变化的两种情况.图6为T1不变，T0=-90°～90°（间隔 15°）变刚度壁板发散失稳特性随T0的变化规律.图7 为T0不变，T1= -90°～90°（间隔15°）变刚度壁板发散失稳特性随T1的变化规律.

从图6中可以看出：当T1不变、T0变化时，变刚度弹性、黏弹性壁板发散速度随着T0的增大均呈先减小后增大的趋势，呈“凹”形分布；对变刚度弹性壁板，在T0从-90°～-60° 时，壁板在亚音速已不发生失稳，T0从 15°～45° 时，壁板发散失稳速度出现极小值范围，基本关于T0= 30° 对称分布；对于变刚度黏弹性壁板，其发散失稳速度分布规律与弹性壁板类似，只是变刚度黏弹性壁板发散失稳速度极小值相比变刚度弹性壁板的小.

图6 曲线纤维变刚度壁板失稳速度随T0的变化（T1不变）Fig.6 Changes of divergence characteristics with T0 for variable-stiffness plate with curved fibers (T1 unchanged)

从图7中可以看出：当T0不变、T1变化时，变刚度弹性、黏弹性壁板壁板发散速度随着T1从-90°～90° 也呈“凹”形分布；壁板发散失稳速度关于T0=0° 对称分布，发散失稳速度变化相对平缓.

图7 曲线纤维变刚度壁板失稳速度随T1的变化（T0不变）Fig.7 Changes of flutter characteristics with T1 for variablestiffness plate with curved fibers (T0 unchanged)

由此可见，复合变刚度壁板曲线纤维参数T0、T1对变刚度弹性、黏弹性壁板发散稳定性均具有明显影响，而且T1影响更为明显一些.此外，曲线纤维复合变刚度壁板在发散稳定性方面存在较大的优化空间，通过调整曲线纤维的路径可以提高复合材料变刚度壁板的发散速度.

5） 黏弹性材料阻尼系数的影响

图8给出了黏弹性材料阻尼系数对曲线纤维复合壁板发散稳定性的影响.从图可以看到：随着黏弹性材料阻尼系数增大，曲线纤维复合壁板的失稳速度减小.

图8 黏弹性阻尼对曲线纤维变刚度壁板发散临界速度影响Fig.8 Influence of viscoelastic damp on divergence critical speed of variable-stiffness panels with curved fibers

6） 边界条件影响

图9给出了不同约束情况下，变刚度弹性、黏弹性壁板的发散稳定性情况.壁板铺层形式为[0/90/＜T0|T1＞/＜-T0|-T1＞]s.从图中可以看出：约束条件对变刚度壁板失稳速度影响很明显，边界约束越强，壁板的稳定性越好，发散失稳临界速度越大.

图9 约束对曲线纤维变刚度壁板发散失稳临界速度影响Fig.9 Influence of boundary conditions on divergence critical speed of variable-stiffness panels with curved fibers

4 结 论

本文建立了亚音速流场中曲线纤维变刚度弹性、黏弹性壁板气动弹性稳定性分析模型，与文献[12]、文献[16]结果进行对比，验证了本文模型及算法的正确性，进而采用复模态理论求解复合变刚度壁板发散临界速度，在此基础上进行了关键参数影响分析，得出的主要结论如下：

1） 变刚度壁板材料的黏弹性效应会导致其失稳速度降低，与弹性壁板相比，黏弹性壁板失稳前固有频率值变化更剧烈，发散失稳出现得更突兀；

2） 通过调整铺层主方向角可以实现复合材料变刚度弹性、黏弹性壁板的发散速度的提高，通过调整曲线纤维的路径也可以提高复合材料变刚度壁板的发散速度，因此，曲线纤维复合变刚度壁板在发散稳定性方面存在较大的优化设计空间；

3） 弹性、黏弹性壁板发散速度随着曲线纤维路径参数T0、T1从 -90° ～ 90° 变化呈现“凹”形分布，且参数T1对壁板发散稳定性的影响相比T0更明显；

4） 边界条件对曲线纤维复合材料弹性、黏弹性壁板的发散临界速度有较大影响，边界约束越强，壁板的稳定性越好，发散失稳临界速度越大.

致谢：民航航空器适航审定技术重点实验室开放基金（SH2020112705）.