摘要
针对航空发动机整体叶盘设计在进一步降低发动机质量、减少连接件复杂性的同时,也带来阻尼不足引起的振动问题,首先,提出了一种基于叶盘切缝并辅以填充物的方法,使系统在振动过程中产生非弹性碰撞能量损耗以及填充材料的阻尼损耗,从而为整体叶盘振动提供额外的阻尼;其次,建立了切缝整体叶盘与切缝填充后的时变参数动力学简化模型;最后,分别从仿真与实验的角度对切缝整体叶盘与切缝填充后的自由振动响应进行了分析,并利用希尔伯特变换法,识别了整体叶盘结构的瞬时频率以及瞬时阻尼比。结果表明,整体叶盘切缝的碰撞与填充物可以增加系统阻尼,能够起到抑制振动的作用。
随着材料科学的发展以及相关加工工艺的进步,航空发动机整体叶盘逐渐取代榫接叶盘。与榫接叶盘相比,整体叶盘具有轮盘薄、质量轻等优点,且没有榫头与榫槽的摩擦与碰撞,因此降低了材料的磨损。但是,整体叶盘的轮盘及叶片质量与刚度相差较小,容易发生耦合振
针对整体叶盘振动抑制,最常用的减振方法是从增加摩擦阻尼或碰撞增加额外阻尼入
许多学者对基于碰撞增加额外阻尼振动抑制的方法也进行了研究。Kaneko
笔者提出了一种新的整体叶盘阻尼减振方法,通过在整体叶盘的设计过程中进行切缝处理,并补以高阻尼材料,使得叶片受外激励下振动时,叶盘间隙发生碰撞和摩擦,从而引起振动能量的耗散,为整体叶盘提供了额外阻尼,实现整体叶盘振动抑制的目的。首先,根据分段非线性模型,建立了叶盘扇区的动力学模型,并进行数值仿真模拟分析;其次,开展扇区、整体叶盘模态测试实验,并基于希尔伯特变换法进行参数辨识;然后,分析了不同碰撞条件下的结构阻尼变化规律,并依据模态测试结果修正仿真模型参数;最后,总结了所提出的整体叶盘减振方法的可行性,为后续整体叶盘结构减振设计提供思路。
叶盘切缝如

图1 叶盘切缝
Fig.1 Model of a blisk with slitting

图2 带间隙非线性振动模型
Fig.2 Dynamics model of a blisk with slitting
叶盘模型的运动微分方程如下
(1) |
(2) |
(3) |
其中:m为质量;k1,k2为系统刚度;c1,c2为系统阻尼;k3,k4分别为系统左右两边的碰撞刚度;δ0为细缝间隙;F为作用于系统的外力;a0为外力的幅值;f00为外力的激励频率; k (x)为结构恢复力;x为沿轮盘切线方向位移,即叶盘周向位移。
通过调整参数大小,可以模拟分析叶片在不同参数下的振动响应规律。
对上述运动微分方程做如下归一化处理
(4) |
(5) |
(6) |
(7) |
(8) |
(9) |
其中:f0为固有频率;k0=2k1=2k2;c0=2c1=2c2;β,γ为无量纲参数。
(10) |
(11) |
(12) |
其中:为分段线性每一段的固有圆频率;Fi为分段线性每一段的等效外激励N。
基于希尔伯特变换的非线性参数识别方
具有等价黏性阻尼的线性时变单自由度系统的自由振动微分方程为
(13) |
若已知结构在自由振动条件下的位移响应、速度响应、加速度响应以及对应的希尔伯特变换,就能识别出结构的时变模态参数,其计算式为
(14) |
(15) |
其中:f0(t) 为无碰撞条件下的线性系统的固有频率;(t)为阻尼比;m(t)为质量;和为加速度响应及其希尔伯特变换;和为速度响应及其希尔伯特变换;x(t)和为位移响应及其希尔伯特变换。
为简化表达,后文中速度响应表示为v(t),加速度响应表示为a(t)。
对上述单自由度系统,分别计算线性与非线性系统的自由衰减响应,为后续分析做参考。令外激励Fi=0 N,固有频率f0=3.5 Hz,阻尼比=1.55%,u0=0.021,γ=1。设系统的初始位移x(0)=5u0,初始速度v(0)=0。先取α=0,即碰撞刚度为0,此时系统为单自由度线性系统,即无碰撞引起的刚度非线性。线性与非线性状态下瞬时响应与参数识别结果对比如

图3 线性与非线性状态下瞬时响应与参数识别结果对比
Fig.3 Identified parameters from linear and nonlinear systems
根据
由于参数识别假定了非时变参数的线性结构,所以图中所识别的系统参数(如瞬时频率以及瞬时阻尼)均为时变函数。另外,由于位移与时间的关系已知,所以所识别的时变参数也可以转化成位移的函数,如
为了验证碰撞减振应用在整体叶盘振动抑制中的可行性,设计了模态测试实验来获取切缝叶盘的模态参数,通过不同结构、材料的叶盘模态参数比较,判断碰撞减振的实际效果。
实验用力捶对切缝整体叶盘结构施加脉冲激振力,采用NI数据采集仪对叶盘加速度数据进行采集,实验测试实况及部分仪器如

图4 实验测试实况及部分仪器
Fig.4 Experimental setup and instrumentations
利用带通滤波以及希尔伯特变换,可以对分析频率范围内任意阶模态实施非线性参数识别。本研究仅对叶盘的第1阶模态开展分析与讨论。
由于现有工艺的限制,叶盘加工缝隙能实现的最小尺寸为0.2 mm,所以仅依靠锤击激励产生的振动很难发生碰撞,因此在锤击实验时叶盘缝隙间填充了不同材料以增加叶盘的等效阻尼。
对铝合金填充的扇区叶盘进行实验,得到原始加速度响应信号及其傅里叶频谱如

图5 原始加速度响应信号及其傅里叶频谱
Fig.5 Raw acceleration response and its Fourier spec trum

图6 滤波后加速度、速度及位移响应信号
Fig.6 Acceleration, velocity and displacement response signal after filtering fil tering
扇区叶盘非线性参数识别结果如

图7 扇区叶盘非线性参数识别结果
Fig.7 Sector blisk nonlinear parameter identification results
因前述仿真模型并未考虑叶盘缝隙中填充材料的影响,为了尽可能使仿真模型更加贴近实验条件下的真实叶盘结构,本节基于模态实验结果对仿真模型进行修正。
缝隙填充后叶盘仍然可以看作是一个非线性系统,其非线性主要来自于以下方面:①填充材料本身的非线性;②填充材料与缝隙接触面在相对运动过程中的拉压非对称性。依据文献[

图8 修正后叶盘模型
Fig.8 The revised blisk model
修正叶盘模型的运动微分方程如下
(16) |
(17) |
(18) |
(19) |
(20) |
其中:c (x)为结构时变阻尼;k (x)为修正前结构恢复力;x为沿轮盘切线方向位移;为阻尼比;kc(x)为修正后结构刚度。
将参数识别结果代入修正后的模型中并与实验结果进行对比,如

图9 修正模型后仿真与实验结果对比
Fig.9 Comparison of the parameter identifications between the modified simulation and the test
为验证不同的填充材料对等效阻尼的影响,本次实验采用了3种不同的填充材料:铝膜、纸以及橡胶。时域响应如

图10 时域响应
Fig.10 Time domain response

图11 叶盘非线性参数识别结果
Fig.11 Nonlinear parameter identification results of blisk
状态 | 频率 均值/Hz | 频率均值变化率/% | 阻尼 均值/% | 阻尼均值变化率/% |
---|---|---|---|---|
未切缝 未填充 |
850.00 821.00 |
0 ‒3.41 |
0.48 0.52 |
0 8.33 |
铝合金 | 829.50 | ‒2.41 | 0.75 | 56.25 |
纸制品 | 826.50 | ‒2.76 | 1.13 | 135.42 |
橡胶 | 823.00 | ‒3.18 | 1.41 | 193.75 |
对比
1) 所设计的叶盘切缝对所分析的1阶频率影响不大(小于4%)。
2) 可以通过切缝的碰撞增加振动的能量耗散,从而提高叶盘的阻尼比,但所增加的阻尼跟切缝的大小有关。
3) 切缝填充能进一步提高叶盘的阻尼,且阻尼大小与填充物的材料性质有关。本研究测试的切缝填充材料中,阻尼性能最好的橡胶将叶盘的等价阻尼提高到了未作处理时整体叶盘的约3倍。
4) 由于目前实验室条件的限制,实验仅在叶盘静态时进行,未考虑旋转状态下离心力、转速等指标,后续研究可以考虑模拟更真实的叶盘工作状态进行测试验证。
5) 整体叶盘切缝后,在缝隙底部会造成较为严重的应力集中,需要优化设计切缝以避免此现象。
6) 进一步优化切缝与切缝填充物,以优化叶盘的多模态阻尼。
参考文献
陈玉刚.整体叶盘振动特性分析与硬涂层阻尼减振设计研究[D].大连:大连理工大学,2017. [百度学术]
张亮,刘铁箭,李欣,等.失谐叶盘结构减振问题研究综述[J].河北科技大学学报,2016,37(2):109-117. [百度学术]
ZHANG Liang, LIU Tiejian, LI Xin, et al. Research overview on vibration damping of mistuned bladed disk assemblies[J].Journal of Hebei University of Science and Technology,2016, 37(2):109-117.(in Chinese) [百度学术]
许愕俊,田忠贤. 论叶片摩擦阻尼减振[J].航空发动机, 2000(1):49-54. [百度学术]
XU Ejun, TIAN Zhongxian. On friction damping of blades [J]. Aeroengine, 2000(1):49-54. (in Chinese) [百度学术]
CIGEROGLU E, OZGUVEN H N. Nonlinear vibration analysis of bladed disks with dry friction dampers [J]. Journal of Sound and Vibration,2006,295(3/4/5):1028-1043. [百度学术]
LAXALDE D, THOUVEREZ F, LOMBARD J P. Vibration control for integrally bladed disks using friction ring dampers[C]∥ASME Turbo Expo 2007: Power for Land, Sea and Air. Montreal, Canada:ASME,2007. [百度学术]
LAXALDE D, THOUVEREZ F, SINOU J J, et al. Qualitative analysis of forced response of blisks with friction ring dampers[J]. European Journal of Mechanics-A/Solids, 2007, 26(4):676-687. [百度学术]
PATSIAS S, SAXTON C, SHIPTON M. Hard damping coatings: an experimental procedure for extraction of damping characteristics and modulus of elasticity[J]. Materials Science and Engineering: A, 2004, 370(1/2):412-416. [百度学术]
ALLARA M. A model for the characterization of friction contacts in turbine blades[J]. Journal of Sound and Vibration,2009, 320(3):527-544. [百度学术]
WEI S T, PIERRE C. Effects of dry friction damping on the occurrence of localized forced vibrations in nearly cyclic structures[J]. Journal of Sound and Vibration,1989,129(3):397-416. [百度学术]
陈香,朱靖,李光辉,等. 涡轮带冠叶片干摩擦阻尼减振实验[J].航空动力学报,2012,27(4):817-823. [百度学术]
CHEN Xiang,ZHU Jing,LI Guanghui,et al. Experiment on dry friction damping of shrouded turbine blades[J]. Journal of Aerospace Power, 2012,27(4):817-823. (in Chinese) [百度学术]
陈璐璐,张大义,文敏,等.带凸肩风扇叶片振动特性及设计方法研究[J].推进技术,2015,36(9):1389-1394. [百度学术]
CHEN Lulu, ZHANG Dayi, WEN Min, et al. Dynamical effects of shrouds on fan blade vibration and its corresponding design method[J]. Tournal of Propulsion Technology, 2015,36(9):1389-1394. (in Chinese) [百度学术]
KANEKO Y, MORI K, OHYAMA H. Development and verification of 3000 rpm 48 inch integral shroud blade for steam turbine[J]. JSME International Journal, 2006,49(2):205-211. [百度学术]
LU Z, LU X L, MASRI S F. Studies of the performance of particle dampers under dynamic loads [J]. Journal of Sound and Vibration,2010,26(20):5415-5433. [百度学术]
LU Z, MASRI S F, LU X L. Studies of the performance of particle dampers attached to a two-degrees-of-freedom system under random excitation [J]. Journal of Vibration and Control,2011,17(10):1454-1471. [百度学术]
卢绪祥,黄树红,刘正强,等.汽轮机自带冠叶片碰撞减振的研究现状与发展[J].振动与冲击,2010,29(2) :11-16. [百度学术]
LU Xuxiang, HUANG Shuhong, LIU Zhengqiang, et al. Impact-damping studies for integral shroud blades of a turbine machinery[J]. Journal of Vibration and Shock, 2010,29(2):11-16. (in Chinese) [百度学术]
卢绪祥,刘正强,黄树红,等.含间隙碰撞振动系统的非线性振动特性[J].动力工程学报,2012,32(5) :388-393. [百度学术]
LU Xuxiang, LIU Zhengqiang, HUANG Shuhong, et al. Nonlinear vibration characteristics of a vibro-impact system with clearance[J]. Journal of Chinese Society of Power Engineering, 2012,32(5):388-393. (in Chinese) [百度学术]
赵亚英,陈长秀,赵熹.汽轮机带冠叶片阻尼振动特性的有限元分析[J].安徽理工大学学(自然科学版), 2014,34(1):61-67. [百度学术]
ZHAO Yaying, CHEN Changxiu, ZHAO Xi. The finite element analysis on damping vibration dynamic characteristics for shrouded blades of steam turbine[J]. Journal of Anhui University of Science and Technology( Natural Science),2014,34 (1):61-67. (in Chinese) [百度学术]
FELDMAN M. Non-linear system vibration analysis Hilbert transform-I: free vibration analysis method 'freevib'[J]. Mechanical Systems and Signal Process⁃ing, 1994,8(2):119-127. [百度学术]
FELDMAN M. Non-linear system vibration analysis Hilbert transform-II:forced vibration analysis method 'forcevib'[J]. Mechanical Systems and Signal Process⁃ing, 1994, 8(3):309-318. [百度学术]
FELDMAN M, BRAUN S. Identification of non-linear system parameters via the instantaneous frequency: application of the Hilbert transform and Wigner-Ville technique[C]∥Proceedings of the 13th IMAC. Nashville,TN,USA:SEA,1995:637-642. [百度学术]
FELDMAN M. Identification of weakly nonlinearities in multiple coupled oscillators[J]. Journal of Sound and Vibration, 2007,303:357-370. [百度学术]
FELDMAN M. Time-varying vibration decomposition and analysis based on the Hilbert transform [J]. Journal of Sound and Vibration, 2006, 295:518-530. [百度学术]