Analysis of blade vibration response based on contact stiffness distribution characteristics
-
摘要:
根据接触刚度在锯齿形叶冠结合面上非均匀分布的特性,提出一种基于定义和有限元计算相结合的接触刚度计算方法。在此基础上,将微-宏滑动摩擦模型作为叶冠结合面处的摩擦模型,推导微滑动和宏滑动状态下摩擦力表达式。利用谐波平衡法将非线性摩擦力转化为等效刚度和等效阻尼进行振动特性分析。针对叶冠结合面相对位移幅值动态变化的特点,提出了一种迭代求解的振动响应分析方法。与文献提供的带摩擦阻尼结构悬臂梁振动响应实验数据相比,共振峰附近的振幅误差为1.99 mm,相对误差为3.9%。振幅的最大误差为5.53 mm,出现在远离共振峰的位置,验证了响应分析方法的可行性。将该方法应用于带冠叶片上,结果表明当激振力频率为812.3 Hz时,振幅为0.56 mm。
Abstract:According to the non-uniform distribution characteristics of contact stiffness between interfaces of the zigzag shroud, a method for calculating contact stiffness was proposed based on the combination of definition and finite element calculation. On this basis, the micro-macro-slide friction model was applied to shroud contact interfaces, and the frictional force expressions were derived in micro-slide state and the macro-slide state. The harmonic balance method was used to convert the nonlinear friction force into the equivalent stiffness and equivalent damping. Considering the dynamic change of the relative displacement amplitude between interfaces, an iterative solution method for analyzing vibration response was proposed. Compared with experimental data of vibration response of cantilever beam with frictional damping structure provided in literature, the amplitude error near the resonance peak was 1.99 mm, and the relative error rate was 3.9%. The maximum error of the amplitude was 5.53 mm, which appeared far away from the resonance peak. It was believed that the accuracy of the response analysis method was relatively high. This method was applied to the shrouded blade, and the results showed that the amplitude was 0.56 mm when the excitation force frequency was 812.3 Hz.
-
表 1 叶冠静力学分析结果
Table 1. Statics analysis results of blade shroud
项目 结合面网格尺寸/mm 0.5 0.25 0.125 结合面最大应力/MPa 264.74 394.29 502.70 结合面最大位移/mm 0.34319 0.34317 0.34323 位移的相对误差/% −0.0058 0.0175 -
[1] 谢永慧,孟庆集. 汽轮机叶片疲劳寿命预测方法的研究[J]. 西安交通大学学报,2002,36(9): 912-915.XIE Yonghui,MENG Qingji. Numerical model for steam turbine blade fatigue life[J]. Journal of Xi’an Jiaotong University,2002,36(9): 912-915. (in Chinese) [2] 左焕. 整圈自带冠叶片系统接触碰撞特性研究[D]. 哈尔滨: 哈尔滨工程大学, 2014.ZUO Huan. Research on contact collision characteristics for the whole circle integrally shrouded blade[D]. Harbin: Harbin Engineering University, 2014. (in Chinese) [3] 洪杰,杨鑫,陈璐璐,等. 基于接触状态分析的叶冠阻尼评估与改进设计[J]. 北京航空航天大学学报,2012,38(5): 681-687.HONG Jie,YANG Xin,CHEN Lulu,et al. Damping assessment and improving design of blade shrouds based on contact states[J]. Journal of Beijing University of Aeronautics and Astronautics,2012,38(5): 681-687. (in Chinese) [4] OZOLI R,OZOLI I,RAMĀNA L. Turbine blade of gas turbine engine additional unloading by changing the layout of the gravity centre of the shroud shelf[J]. Transport and Aerospace Engineering,2019,7(2): 14-23. [5] 刘丽兰,刘宏昭,吴子英,等. 机械系统中摩擦模型的研究进展[J]. 力学进展,2008,38(2): 201-213.LIU Lilan,LIU Hongzhao,WU Ziying,et al. Research progress of friction model in mechanical system[J]. Advances in Mechanics,2008,38(2): 201-213. (in Chinese) [6] BURDEKIN M,COWLEY A,BACK N. An elastic mechanism for the micro-sliding characteristics between contacting machined surfaces[J]. Journal of Mechanical Engineering Science,2006,20(3): 121-127. [7] JOHNSON K L. 接触力学[M]. 徐秉业, 罗学富, 刘信声, 等, 译. 北京: 高等教育出版社, 1991. [8] MINDLIN R D. Compliance of elastic bodies in contact[J]. Journal of Applied Mechanics,1949,16(9): 259-268. [9] DERESIEWICZ H. Oblique contact of nonspherical elastic bodies[J]. Journal of Applied Mechanics,1957,24(4): 623-624. doi: 10.1115/1.4011612 [10] LIU Yalin,BO Shangguan,XU Zili. A friction contact stiffness model of fractal geometry in forced response analysis of a shrouded blade[J]. Nonlinear Dynamics,2012,70(3): 2247-2257. doi: 10.1007/s11071-012-0615-8 [11] ZHAO Jiazhe, WANG Yanrong, LUO Yanbin, et al. A vibration damping analysis method of turbine blade shroud dampers based on a given eigen-mode[C]//Proceedings of the ASME 2018 International Design Engineering. Quebec, Canada: ASME, 2018: 1-11. [12] 陈璐璐,张大义,文敏,等. 带凸肩风扇叶片振动特性及设计方法研究[J]. 推进技术,2015,36(9): 1389-1394.CHEN Lulu,ZHANG Dayi,WEN Min,et al. Dynaical effects of shrouds on fan blade vibration and its corresponding design method[J]. Journal of Propulsion Technology,2015,36(9): 1389-1394. (in Chinese) [13] BERTHILLIER M,DUPONT C,MONDAL R,et al. Blades forced response analysis with friction dampers[J]. Journal of Vibration and Acoustics,1998,120(2): 468-474. doi: 10.1115/1.2893853 [14] 黄平, 郭丹, 温诗铸. 界面力学[M]. 北京: 清华大学出版社, 2013. [15] 申秀丽,高鹏新,董少静,等. 带锯齿形叶冠叶片接触应力计算及静强度分析[J]. 航空动力学报,2020,35(4): 732-743. doi: 10.13224/j.cnki.jasp.2020.04.007SHEN Xiuli,GAO Pengxin,DONG Shaojing,et al. Contact stress calculation and static strength analysis of zig-zag shrouded blade[J]. Journal of Aerospace Power,2020,35(4): 732-743. (in Chinese) doi: 10.13224/j.cnki.jasp.2020.04.007 [16] ZHANG S,SHEN X,DONG S,et al. A method to analyze the contact stress of dovetail attachments in aeroengine[J]. Journal of Engineering for Gas Turbines and Power,2019,141(2): 022505.1-022505.10. [17] 靳顺. 摩擦阻尼建模及涡轮叶片结构减振分析[D]. 北京: 北京航空航天大学, 2015.JIN Shun. Friction damping modeling and vibration reduction research in gas-turbine structure[D]. Beijing: Beijing University of Aeronautics and Astronautics, 2015. (in Chinese) [18] 张建伟, 白海波, 李昕. ANSYS 14.0超级学习手册[M]. 北京: 人民邮电出版社, 2013: 242-243. [19] 康佳豪,徐超,李东武,等. 基于谐波平衡-时频转换法的摩擦振子稳态响应分析[J]. 振动与冲击,2020,39(12): 170-176, 205.KANG Jiahao,XU Chao,LI Dongwu,et al. Analysis of steady responses for a frictional oscillator based on the harmonic balance-alternating frequency/time domain method[J]. Journal of Vibration and Shock,2020,39(12): 170-176, 205. (in Chinese) [20] 朱永新. 带冠叶片摩擦减振特性研究[D]. 南京: 南京航空航天大学, 2007.ZHU Yongxin. Research on characteristic of the forced response of shrouded blades[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2007. (in Chinese)