2. 沈阳鼓风机集团股份有限公司, 辽宁 沈阳 110869
2. Shenyang Blower Works Group Corporation, Shenyang 110869, China
研究叶片与柔性机匣的碰摩特性, 合适的碰摩模型是非常重要的, 很多学者对此进行了深入研究.最简单的为线性弹簧模型, Padova和Choy等[1-2]推导出叶片沿径向的变形和碰摩力之间关系的碰摩模型.考虑到叶片旋转产生的离心刚化影响, Jiang等[3]对Padova模型进行了修正.刘昕等[4]采用Hertz弹性接触模型研究碰摩力与侵入深度的关系.Chen[5]考虑叶片数和转静间隙变化对碰摩力的影响, 以弹性碰摩模型为基础, 模拟单点、多点、局部,以及全周等多种碰摩形式.Ma等[6]研究了叶片形状和弹性机匣整体变形对碰摩力的影响, 提出了新的叶片机匣碰摩力模型.
以上所有碰摩模型中, 对机匣的考虑过于简单, 文献[1-3]完全没有考虑机匣变形对碰摩力的影响, 文献[4]直接引用Hertz接触模型计算法向接触力, 对于机匣这种薄壁结构, 无法反映出其厚度的影响.文献[5-6]将机匣假设为线弹性支承结构, 文献[7-8]研究了机匣的柔性响应, 但没有将其参数纳入碰摩力的计算中, 依旧采用线性碰摩刚度的算法.
本文将机匣看作薄壁圆环结构[7-9], 其半径、厚度、碰摩宽度、支承刚度都是影响碰摩力大小的重要参数.叶片简化为Timoshenko梁模型[10], 并考虑离心刚化、旋转软化和碰摩软化效应, 研究其在碰摩过程中的弹性变形.
1 旋转叶片-柔性机匣碰摩力模型叶片和柔性机匣的碰摩力取决于两者共同的变形:叶片考虑弯曲变形, 机匣考虑柔性变形及与支承刚度相关的整体位移.在推导法向碰摩力的过程中, 叶片-机匣碰摩看作为准静态接触过程, 假定轮盘为刚性, 不考虑碰摩过程中摩擦热效应.
1.1 旋转叶片柔性变形将叶片简化为Timoshenko梁模型, 如图 1所示.图中给出了叶片受力情况, 包括离心力、气动力、法向接触力, 以及摩擦力.离心力会使叶片产生离心刚化和旋转软化效应, 而法向接触力会使叶片产生软化现象.微元体离心力计算公式为
(1) |
Timoshenko梁存在剪切变形, 如果剪切变形为0, 则中心线的切线将与截面的法线重合, y为中心线的挠度, φ为弯曲引起的斜率, dy/dx-φ为斜率损失, 等于剪切角, Q为截面的剪力, 梁的弹性方程为
(2) |
(3) |
法向接触力沿着半径指向圆心, 摩擦力与其垂直并与相对运动方向相反, 离心载荷沿着半径方向背离圆心, 将这三个力沿水平和竖直方向分解, 由图中关系, 可得
(4) |
式中:L为叶片长度,y(L)为叶尖的弯曲位移.
离心力水平方向分量的力矩与叶片运动方向相反, 产生离心刚化效应.而竖直方向的分量与叶片运动方向相同, 产生软化效应.文献[11]在微元体受力分析中描述了旋转软化力, 与式(4)表达式完全一致.
在梁的任意截面x0处, 分析其力矩平衡关系, 各分力及其力臂如图 2所示, 可得叶片弯曲位移方程:
(5) |
式(5)是二阶非齐次微分方程, 直接求解比较困难, 采用幂级数方法展开求解, 并略去高阶小量, 设定叶片机匣摩擦系数为μ, 叶片安装角为β, 则有Ft=μFn, 可得
(6) |
式中
如果忽略离心刚化、旋转软化和碰摩软化效应, 直接引用材料力学公式, 可得叶片弯曲位移为
(7) |
文献[6]在推导叶片-机匣碰摩力模型中考虑了叶片的离心势能, 基于功能互等定理推导出叶片的弯曲位移.为了验证本文叶片挠度计算式(6)的准确性, 将文献[6]的计算模型和式(7)的计算结果与有限元结果进行对比, 有限元仿真在ANSYS中进行, 采用Beam188单元, 叶盘系统参数见文献[6].
综合考虑离心刚化、旋转软化和碰摩软化影响, 将转速设定为5 000 r/min, 图 3a中碰摩力为2 kN, 此时软化效应大于离心刚化效应, 实际弯曲位移大于3个影响因素都不考虑的结果.图 3b碰摩力设定为200 N, 此时软化效应小于离心刚化效应, 实际弯曲位移也小于材料力学模型计算结果.图 3a中文献[6]结果与本文相差较大, 而图 3b中则极为接近.通过对比可以看出, 文献[6]充分考虑了离心刚化效应, 而未过多考虑碰摩软化效应.
航空发动机机匣的简化模型如图 4所示, 叶片的宽度为b, 安装角为β, 则机匣沿轴向与叶片的碰摩宽度为Lc=bcosβ.取位于碰摩宽度内的这段机匣作为研究对象, 其等效支承刚度可以通过梁函数法[5]计算或实验测得, 柔性变形则采用圆环来模拟, 其任意点的最终位移为机匣总体位移和柔性位移的向量和, 共包括水平位移uc, 垂直位移vc, 径向柔性位移us(向心为正), 切向柔性位移ws(沿角度增大方向为正).
圆环的弹性势能为
(8) |
式中:Ec为弹性模量; Rc为机匣半径; Ic为圆环截面惯性矩, 与圆环厚度和碰摩宽度有关.依据功能互等定理, 当径向碰摩力为Fn、作用角度为φ时, Fn所做的功为
(9) |
圆环的柔性变形通常采用节径振动形式来表达, 以n≥2的所有节径振动作为模态振型, 其切向位移和径向位移可以展开为
(10) |
式中n代表节径数, 也称为周向波数.将式(9)按照节径位移展开, 可得
(11) |
并得机匣结构刚度为
(12) |
为了验证任意点的总体位移为平动位移和柔性位移的向量和, 本文将机匣展开至4阶.某机匣半径为224 mm, 厚度为3 mm, 碰摩长度50 mm, 施加恒定径向激励F=1 N, 与水平方向夹角为60°, 在kcx=kcy=∞和kcx=kcy=20 kN/m时求取最大位移, 并与有限元仿真进行对比, 如表 1所示.两种工况下绘制机匣形状, 所有位移均放大1 000倍.
叶片挠度y(x)是碰摩力的函数, 由参考文献[6], 叶片径向位移可表达为
(13) |
根据弹性协调相容条件, 叶片-机匣在碰摩过程中侵入深度δ可以分解为3部分:叶片径向位移uL、沿碰摩法向机匣的平动位移ucn和弹性位移usn, 如图 5所示.可得
(14) |
将式(12)和式(13)代入式(14), 略去高阶小量, 可得法向碰摩力的表达式为
(15) |
其中
为了验证上述模型的准确性, 搭建了叶片转子碰摩实验台.实验台由动力系统、转子系统、进给系统, 以及测试系统组成.实验台详细参数及支承刚度见文献[6].由于条件所限, 本实验不能实现机匣的柔性变形, 只能考虑叶片以及机匣支承刚度的影响, 故将式(14)中Ks设为无穷大.
实验在1 000, 1 500和2 000 r/min三种转速下进行.机匣采用钢制和铝制两种材料, 铝制机匣的支承刚度为2×107 N/m, 钢制机匣的支承刚度为3.5×107 N/m.叶片分为厚度为3 mm的薄叶片和5 mm的厚叶片两种, 本文与文献[6]模型的对比如图 6所示.各工况下本文的精度略高于文献[6]模型.对比图 6a和6b可知, 增加叶片厚度将提高抗弯刚度; 侵入量相同时, 碰摩力会大幅增加.由图 6c和6d以及图 6e和6f的对比可知, 侵入量相同时, 由于支承刚度大, 钢制机匣的法向碰摩力大于铝制机匣, 对比图 6c和6e以及图 6a, 6d和6f可以看出, 随着转速的提高, 离心刚化作用增加叶片的抗弯刚度, 保持侵入量不变, 碰摩力会少许增加.
文献[12]对航空发动机整机进行建模, 并通过力传感器测得碰摩力.为进一步验证机匣柔性变形与相关参数的影响, 本文与其进行对比.固定转速和支承刚度不变, 随着碰摩程度逐渐加重, 碰摩力逐渐增加, 实验结果如图 7a所示.本文在机匣厚度hc分别为3,4和5 mm三种工况下计算碰摩力与侵入深度的关系, 如图 7b所示, 可以看出, 由于机匣较薄, 其结构刚度远小于支承刚度, 法向碰摩力与侵入深度呈现出线性变化的关系.当机匣厚度hc为4 mm时, 本文模型中法向碰摩力的计算结果与实验结果最为接近.
1) 叶片的弯曲变形同时受到离心刚化、旋转软化和碰摩软化的作用, 不同转速及碰摩力条件下, 叶片的弯曲特性不同:提高转速会增加叶片的抗弯刚度, 碰摩力的增加会降低叶片的抗弯刚度;叶片实际工作中的状态由二者共同决定.
2) 叶片-机匣碰摩中, 碰摩部位机匣的位移由两部分组成:整体机匣位移和局部柔性位移.
3) 叶片-机匣的碰摩力受叶片刚度、机匣支承刚度和结构刚度的影响:侵入量相同时, 机匣支承刚度和结构刚度的提高会增加碰摩力的值, 其中机匣的结构刚度, 主要由机匣的厚度、半径和碰摩宽度决定.
[1] |
Padova C, Barton J, Dunn M, et al.
Experimental results from controlled blade tip/shroud rubs at engine speed[J]. Journal of Turbomachinery, 2007, 129(4): 1165–1177.
|
[2] |
Padovan J, Choy F K.
Nonlinear dynamics of rotor/blade/casing rub interactions[J]. Journal of Turbomachinery, 1987, 109(4): 527–534.
DOI:10.1115/1.3262143 |
[3] |
Jiang J, Ahrens J, Ulbrich H, et al.A contact model of a rotating rubbing blade[C]// Proceedings of the 5th International Conference on Rotor Dynamics.Darmstadt: Technische Universität Darmstadt, 1998: 478-489.
|
[4] |
刘昕, 张华彪, 孙小磊, 等.
叶轮转子碰摩的非线性动力学响应[J]. 振动与冲击, 2016, 35(20): 172–177.
( Liu Xin, Zhang Hua-biao, Sun Xiao-lei, et al. Nonlinear dynamics on rubbing of a blade-rotor system[J]. Journal of Vibration and Shock, 2016, 35(20): 172–177. ) |
[5] |
Chen G.
Simulation of casing vibration resulting from blade-casing rubbing and its verifications[J]. Journal of Sound and Vibration, 2016, 361: 190–209.
DOI:10.1016/j.jsv.2015.09.041 |
[6] |
Ma H, Tai X Y, Han Q K, et al.
A revised model for rubbing between rotating blade and elastic casing[J]. Journal of Sound and Vibration, 2015, 337: 244–262.
DOI:10.1016/j.jsv.2014.10.010 |
[7] |
Lesaffre N, Sinou J J, Thouverez F.
Model and stability analysis of a flexible bladed rotor[J]. International Journal of Rotating Machinery, 2006, 2006: 1–16.
|
[8] |
Parent M O, Thouverez F.
Phenomenological model for stability analysis of bladed rotor-to-stator contacts[J]. International Symposium on Transport Phenomena and Dynamics of Rotating Machinery, 2016, 10/11/12/13/14/15: 1–13.
|
[9] |
Lesaffre N, Sinou J J, Thouverez F.
Stability analysis of rotating beams rubbing on an elastic circular structure[J]. Journal of Sound and Vibration, 2007, 299(4/5): 1005–1032.
|
[10] |
Chen G.
Vibration modeling and verifications for whole aero-engine[J]. Journal of Sound and Vibration, 2015, 349: 163–176.
DOI:10.1016/j.jsv.2015.03.029 |
[11] |
Chu S M, Cao D Q, Sun S P, et al.
Impact vibration characteristics of a shrouded blade with asymmetric gaps under wake flow excitations[J]. Nonlinear Dynamics, 2013, 72(3): 539–554.
DOI:10.1007/s11071-012-0732-4 |
[12] |
李勇, 姜广义, 王德友, 等.
转静件碰摩状态下的叶片振动载荷和振动特性测试分析[J]. 航空动力学报, 2008, 23(11): 1988–1992.
( Li Yong, Jiang Guang-yi, Wang De-you, et al. Test analysis of blades vibration load and vibration characteristics on the condition of rotor and stator rubbing[J]. Journal of Aerospace Power, 2008, 23(11): 1988–1992. ) |