整体叶盘在气体激励下承受宽频激励, 而且轮盘与叶片的刚度非常相近, 这会促使其振动局部超标, 严重威胁发动机的安全性与可靠性^[1].阻尼减振是提高叶片寿命的较为可行的方案.由于整体叶盘的一体化结构特性, 对其实施干摩擦阻尼减振^[2]非常困难.为此, 必须设计一种全新的阻尼器实现整体叶盘的减振.

硬涂层主要是指由金属基、陶瓷基或两者混合制成的涂层材料, 具有抗高温、抗摩擦和抗腐蚀的能力.近年来发现, 硬涂层还能够降低构件的振动应力^[3-5].所以, 本文考虑在叶片涂敷硬涂层来实现整体叶盘的减振.事实上, 这种减振理念已经被NASA的研究人员初步应用并取得很好的效果^[6].为了有效实施整体叶盘的硬涂层阻尼减振, 需要对涂层叶盘进行系统化研究, 主要包括分析模型的建立与动力学分析.

集中和连续参数模型^[7-8]实际结构具有较大差距, 所以大量研究集中于有限元模型^[9].航空发动机整体叶盘是大型复杂的机械结构件, 如果直接对其全尺寸有限元模型进行分析通常需要耗费大量计算资源.对于具有失谐特征的整体叶盘, 其有限元模型的减缩已成为结构振动特性分析的首要问题^[10].子结构模态综合法以瑞利-里兹法为理论基础, 其基本思想是依据先分解后综合的原则^[11].

本文选取叶片双面涂敷NiCoCrAlY+YSZ硬涂层的整体叶盘为研究实例, 推导了固定界面动刚度模态综合法的数学表达式, 并对涂层叶盘的有限元模型进行了减缩; 利用有限元减缩模型获取涂层叶盘的模态参数, 并从计算精度与效率两方面作出评估; 对涂敷硬涂层前后的整体叶盘固有频率、模态损耗因子以及频响作出比较, 探究了硬涂层对整体叶盘的动力学特性的影响.

1 涂层叶盘的减缩建模

在不考虑系统阻尼影响的情况下, 则涂层叶盘的动力学方程可以表示为

(1)

其中:K_cs与M_cs分别表示涂层叶盘的刚度矩阵与质量矩阵; X_cs表示位移向量; F_cs表示激励力.

根据模态综合法的基本思想, 将涂层叶盘分离为轮盘子结构和涂层叶片子结构, 则子结构的刚度矩阵K^r、质量矩阵M^r与位移向量X^r分别表示为

(2)

(3)

(4)

式中, 右上标r=b或d, 且b和d分别表示叶片-硬涂层子结构和轮盘子结构; 右下标i和j分别代表子结构的内部节点与界面节点.

为便于后续描述, 考虑将子结构的K^r，M^r与X^r分别用K，M与X表示.当任一无阻尼子结构处于频率为ω的自由振动时, 其仅在子结构界面力作用f下的运动方程为

(5)

将式(5)展开并进行傅里叶变换, 可得

(6)

子结构在固定界面的情况下, 外部位移向量X_j与内部位移向量X_i分别满足以下条件:

(7)

式中, p和x^p分别表示模态阶次和模态坐标.

由式(7)可得固有频率Λ_ii与模态振型Φ_ii, 分别表示为

(8)

式中, l与h表示低阶模态与高阶模态.

固有频率与模态振型满足正交归一化条件, 即

(9)

式中, I表示单位矩阵.则由式(9)可求得

(10)

对式(10)第三行两端分别求逆, 可得内部自由度的动柔度表达式为

(11)

由式(6)得子结构界面力与位移间的关系, 即

(12)

式中, D_jj^*(λ)表示子结构界面动刚度, 且有

(13)

将式(8)和式(10)代入式(11), 整理可得

(14)

将式(14)代入式(13), 并利用模态正交关系可以得到

(15)

在实际应用中仅仅需要考虑部分阶次的模态特性, 当需要考虑的频率为ω时, 可将Δ_ih^-1(λ)展开为关于特征值λ的泰勒级数, 即

(16)

考虑到子结构固有频率与模态振型满足关系

(17)

将式(16)与式(17)代入式(15), 整理可得

(18)

式中,

将式(18)代入式(12), 整理可得

(19)

(20)

为了表述方便, 考虑两个子结构的模态综合并分别用下标α和β表示.由式(19)可以得到

(21)

为了弥补单一综合协调条件的不足, 本文同时采用界面力与界面位移的双重综合协调条件, 即

(22)

利用式(22), 由式(19)和(20)可得

(23)

(24)

将式(23)和式(24)联立可得

(25)

式中,

将式(25)的第一行进行展开, 并在展开式的两端同时乘以λⁿ(1≤n≤z, n∈Z), 整理可得

(26)

略去式(26)中的λ^m(m>z)项, 整理可得

(27)

(28)

式中,

将式(24)与式(28)联立并整理可得

(29)

式中,

将式(27)代入式(29)并通过整理可得

(30)

式中,

2 涂层叶盘的振动分析

为了精确描述涂层叶盘的动力学特性, 必须考虑系统阻尼对它的影响.经典模态应变能法是含阻尼结构有限元模态分析的能量分析方法, 以耗散能量与总的变形能的比值来确定涂层叶盘的损耗因子^[12].

硬涂层材料损耗因子为η_c, 涂层叶盘有限元模型共有N个单元, 其中硬涂层有M个单元; 叶盘材料损耗因子为η_b, 涂层叶盘耗散能量表示为

(31)

式中, U_{k, p}与U_{g, p}分别表示硬涂层与叶盘的每个单元的模态应变能.由于叶盘的材料损耗因子远远小于硬涂层的材料损耗因子, 所以忽略叶盘材料的阻尼, 则涂层叶盘的耗散能量可以表示为

(32)

涂层叶盘的总变形能V_{z, p}可以表示为

(33)

根据模态应变能理论的定义, 可得到涂层叶盘的模态损耗因子η_{cs, p}, 即

(34)

式中, MSE_{c, p}与MSE_{cs, p}分别表示第p阶的硬涂层的模态应变能与涂层叶盘的模态应变能.

本文利用比例阻尼来近似表示涂层叶盘的系统阻尼c, 并通过如下公式求解, 即

(35)

式中:δ与ξ分别表示与质量矩阵和刚度矩阵相关的比例系数; λ_{cs, p}与ω_{cs, p}分别表示涂层叶盘在第p阶的特征值以及对应的固有频率.

考虑系统阻尼c, 则由式(1)的傅里叶变换可得涂层叶盘的频域动力学方程, 即

(36)

由式(36)可得到涂层叶盘的频响函数, 即

(37)

式中:χ_p表示与固有频率ω_{cs, p}对应的特征向量; ω表示激励频率; i表示虚数单位.

3 涂层叶盘的基本参数

图 1所示为涂层叶盘的实验试件.图 2所示为涂层叶盘的有限元模型(全尺寸和单扇区模型).几何参数、材料参数见表 1和表 2.

3.1 涂层叶盘的模态分析

图 3所示为涂层叶盘的前15阶固有频率.可以看到, 固有频率FC与FB及FA的对应性比较好.其中, 固有频率FA与FC之间的偏差略大, 这是实验干扰造成的, 但在允许范围内.这说明减缩模型具有较高的计算精度.此外发现, 当发生失谐时, 整体叶盘在低频段具有丰富的密集模态.

图 4所示为基于减缩模型与全尺寸模型的运算时间.可以看出, 减缩模型的运算时间较全尺寸模型大大减少, 最大可节省33.3%, 说明减缩模型计算效率比较高.可预见的是, 随着模态阶次的增大, 减缩模型的这种优势变得更加明显.

3.2 硬涂层对整体叶盘的影响

表 3所列为涂敷硬涂层前后整体叶盘前15阶固有频率.可以看到, 涂敷硬涂层后整体叶盘的固有频率较涂敷前略微增大, 变化率分别在0.51%~2.25%与0.71%~2.67%.这说明硬涂层不会对整体叶盘的固有频率造成显著影响.

图 5所示为涂敷硬涂层前后叶盘前15阶模态损耗因子.可以看到, 涂敷硬涂层后的整体叶盘模态损耗因子均明显增大.由于实验受到干扰, 实验测得的模态损耗因子普遍大于理论分析的模态损耗因子.即, A>B, C>D.此外, 同一结构的模态损耗因子有所差异, 但是它们都有相似的变化规律.即, A与C相似, B与D相似.

图 6所示为涂敷硬涂层前后的整体叶盘频响函数曲线.可以直观看到, 在相同阶次下, 整体叶盘在涂层后的频响幅值较涂层前明显减小, 说明叶盘在共振区的共振响应受到了显著抑制.

4 结论

1) 提出了一种叶片涂敷硬涂层的减振方法.这种方法可以有效提高整体叶盘在恶劣环境的工作能力, 而且具有刚度可调节的优点.

2) 建立了涂层叶盘的有限元模型, 推导了固定界面动刚度模态综合法公式对其进行减缩.涂层叶盘的有限元减缩模型具有较高的计算精度与效率, 能够快速获得涂层叶盘的动力学特性.

3) NiCoCrAlY+YSZ硬涂层表现出较强的阻尼能力, 在保证整体叶盘的固有频率不发生较大变化的前提下, 能够显著抑制其在共振区的共振响应, 有效提高失谐整体叶盘的工作可靠性.