1、不同温度下含裂纹压气机叶片的模态特征和振动特性分析不同温度下含裂纹压气机叶片的模态特征和振动特性分析周震霆，贺星，刘永葆(海军工程大学动力工程学院，湖北武汉430033)摘 要:本文选取含裂纹悬臂梁模拟含裂纹压气机叶片，基于无质量扭转弹簧和呼吸式裂纹刚度模型，得到一种不同环境温度下含裂纹压气机叶片的模态特征和振动特性的分析方法。通过弹性模量引入温度模块，利用无质量扭转弹簧连续条件得到关于含裂纹梁的特征方程，分析环境温度和裂纹深度对含裂纹梁固有频率的影响；利用悬臂梁的强迫弯曲振动方程，引入呼吸式裂纹刚度模型，改变激振力频率，分析环境温度及激振频率对含裂纹梁振动位移响应的影响。结果表明，环境温度越
2、高，含裂纹梁的固有频率越小，梁的振动位移响应越大。同时，激振频率的选取也具有一定影响。关键词：压气机叶片；呼吸式裂纹；温度；模态特征；振动特性中图分类号：TK221文献标识码：A文章编号：16727649(2023)12009306doi：10.3404/j.issn.16727619.2023.12.017Analysis of modal characteristics and vibration characteristics ofcracked compressor blades at different temperaturesZHOUZhen-ting,HEXing,LIUYong
3、-bao(SchoolofPowerEngineering,NavalUniversityofEngineering,Wuhan430033,China)Abstract:Inthispaper,acrackedcantileverbeamisselectedtosimulateacrackedcompressorblade.Basedonamass-lesstorsionspringandabreathingcrackstiffnessmodel,amethodforanalyzingthemodalcharacteristicsandvibrationchar-acteristicsofa
4、crackedcompressorbladeatdifferentambienttemperaturesisobtained.Theelasticmodulusisintroducedin-tothetemperaturemodule,andthecharacteristicequationaboutthecrackedbeamisobtainedbyusingthecontinuouscondi-tionofthemasslesstorsionspring,andtheinfluenceoftheambienttemperatureandcrackdepthonthenaturalfrequ
5、encyofthecrackedbeamisanalyzed;theforcedbendingvibrationequationofthecantileverbeamisusedtointroducethebreathingThecrackstiffnessmodelwasusedtochangetheexcitationforcefrequencytoanalyzetheinfluenceofambienttemperatureandexcitationfrequencyonthevibrationdisplacementresponseofthebeamwithcracks.Theresu
6、ltsshowthatthehighertheambienttemperature,thesmallerthenaturalfrequencyofthebeamwithcracks,andthegreaterthevibrationdisplacementresponseofthebeam;atthesametime,theselectionoftheexcitationfrequencyalsohasacertaininfluence.Key words:compressorblade；breathingcrack；temperature；modalcharacteristics；vibra
7、tioncharacteristics0引言压气机叶片是燃气轮机运行过程中完成能量转换的基本单元，也是燃气轮机关键部件中出现故障的最典型部件之一。裂纹是压气机叶片出现故障的常见形式之一，裂纹的存在将会改变叶片的模态特征和振动特性。燃气轮机压气机叶片工作环境恶劣，承受着离心载荷和气动载荷等复杂载荷的影响，特别是高压压气机，随着压比的升高，高压压气机周围气流的温度也不断升高，同时燃烧室通过传热等方式也会导致高压压气机温度升高。在高温的影响下，压气机叶片的材料和力学性能会发生变化，也会对压气机叶片的模第 45卷第12期舰船科学技术Vol.45,No.122023年6月SHIPSCIENCEANDTE
8、CHNOLOGYJun.,2023收稿日期:20220507基金项目:国家科技重大专项作者简介:周震霆(1998)，男，硕士研究生，主要研究方向为动力机械及热力系统的设计、仿真与优化。态和振动特性产生影响。压气机叶片可被看作悬臂梁，通过裂纹悬臂梁模型模拟含裂纹的压气机叶片。对于含裂纹梁的模态特征和振动特性，20 世纪便有人开始进行了研究，并发展至今19，可见研究含裂纹梁的模态特征和振动特性具有实际的工程意义。分析裂纹深度、裂纹位置、结构阻尼、刚度等因素对梁的模态特征和振动特性的影响，却很少考虑温度的影响1012。然而在实际过程中，空气流过各级叶片时，由于受到压缩，压力和温度会逐级升高，位于后面
9、级的压气机，温度产生的影响很难忽略，尤其是在迫切渴求实时监测和在线诊断的今天。若想在燃气轮机工作过程中得到压气机叶片的状态信息，其所处工况的环境温度应考虑在内。本文通过含裂纹悬臂梁模型模拟含裂纹的压气机叶片，采用无质量扭转弹簧等效代替裂纹。通过弹性模量引入温度模块，推导出不同温度下含裂纹悬臂梁的模态特征方程，并通过梁的受迫弯曲振动方程，得到温度影响下梁的振动特性曲线。1模型建立和模态分析 1.1 含裂纹悬臂梁模型建立如图 1(a)所示，建立一段长为 L，宽为 b，高为h 的等截面矩形悬臂梁模拟裂纹叶片，并忽略裂纹尺寸随温度的变化。在该悬臂梁表面有一条深度为 a 的裂纹，该裂纹距离固定端的距离为
10、 R。图1含裂纹叶片等效模型Fig.1AequivalentmodelofcrackedcompressorbladeKT如图 1（b）所示，该含裂纹悬臂梁可以等效成由1 条无质量扭转弹簧连接的 2 段无损弹性梁。对于该悬臂梁的弯曲振动，只考虑等效弹簧的刚度，该裂纹等效刚度的表达式为：1KT=72L(12)(12)ETbh2，(1)=0.629r21.047r3+4.602r49.975r5+20.285r632.993r7+47.041r840.693r9+19.6r10。(2)=l/LETr=a/h式中：L 为悬臂梁的长度，为裂纹的相对位置，为 泊 松 比，为 不 同 温 度 下 的 弹
11、性 模 量；为相对裂纹深度。通过文献 13 得到材料TC11 的弹性模量如表 1 所示。表 1 TC11 不同温度下的弹性模量Tab.1Theelasticmodulusatdifferenttemperature温度T/20100200300400500600弹性模量ET/GPa1231191141101049994 1.2 模态分析图 1 为各向同性的含裂纹等截面均质悬臂梁，则梁的自由振动微分方程为：ETI4y（x,t）x4+A2y(x,t)x2=0。(3)ETy(x,t)式中：为材料密度，为弹性模量，I 为截面惯性矩，为梁的挠度，A 为梁的截面面积。y(x,t)=Y(x)T(t)假定有分
12、离变量解存在，令，代入式（1）分离变量 t 可得：dydxk4y=0。(4)k4=Ag2ETI其中：，A 为截面面积，g为固有频率，则式（4）的通解可成如下形式：|Y1(x)=A1cos(kx)+B1sin(kx)+C1cosh(kx)+D1sinh(kx)，Y2(x)=A2cos(kx)+B2sin(kx)+C2cosh(kx)+D2sinh(kx)。(5)Y(x)表示位移，A1,B1,C1,D1表示裂纹左侧梁与边界条件有关的未知参数，A2,B2,C2,D2表示裂纹右侧梁与边界条件有关的未知参数。悬臂梁的边界条件为：固定端（x=0）的挠度和挠角Y1(0)=0，Y1(x)x=0。94舰船科学技
13、术第45卷自由端（x=l）的弯矩和剪力Y23(x)x2=0，Y32(x)x3=0。假设裂纹所在位置为 x=R，将裂纹等效成无质量扭转弹簧，因此，裂纹处的连续条件为：裂纹两侧挠度相等Y1(x)=Y2(x)。裂纹两侧挠角变化关系Y1(x)x+1KTY1(x)x2=Y2(x)x，裂纹两侧弯矩相等Y21(x)x2=Y22(x)x2，裂纹两侧剪力相等Y31(x)x3=Y32(x)x3。(6)式中，KT为扭转弹簧刚度。由梁的自由振动方程解和梁的边界条件及弹簧连续条件，可得：det(S)=0。(7)其中，矩阵 S 表达式为：S=|101010000000cos(kR)sin(kR)cosh(kR)sin(k
14、R)k/KTcos(kR)cos(kR)k/KTsin(kR)sinh(kR)+k/KTcosh(kR)cos(kR)sin(kR)cosh(kR)sin(kR)cos(kR)sinh(kR)00000100000cos(kl)sin(kl)cosh(kl)sinh(kl)0sin(kl)cos(kl)sinh(kl)cosh(kl)sinh(kR)cosh(kR)sin(kR)cosh(kR)sinh(kR)cosh(kR)+k/KTsinh(kR)sin(kR)cosh(kR)sinh(kR)cosh(kR)sinh(kR)cos(kR)sin(kR)cosh(kR)sinh(kR)cos
15、h(kR)sin(kR)cos(kR)sinh(kR)cosh(kR)。|。根据式（7）可得到 k 值，进而得到裂纹梁的固有频率 g，将得到的 k 值代入到位移方程中，可得到含裂纹梁的模态振型。2振动特性分析 2.1 建立呼吸式裂纹梁弯曲振动方程在图 1（a）的基础上，距悬臂梁固定端 lp处引入大小为 F 的余弦激振力，激振频率为 J，如图 2 所示。图2含裂纹悬臂梁模型Fig.2Thecrackedcantileverbeammodel忽略裂纹对悬臂梁振动的影响，横向力激励下梁的受迫弯曲振动方程为14：ETI4y(x,t)x4+cy(x,t)tdx+A2y(x,t)t2=pcos(Jt)(x
16、lp)。(8)式中：c 为阻尼系数，p 为横向激励力的幅值，lp为激励力离固支端的距离，J为外激励频率，为 Dirac 函数，x 为坐标，t 为时间。y(x,t)=Y(x)T(t)利 用 假 设 的 裂 纹 梁 横 向 固 有 振 动 形 式，根据悬臂梁的边界条件，得到悬臂梁的振型函数为15：Y(x)=(sinkxsinhkx)+n(coskxcoshkx)，(9)n=sinkL+sinhkLcoskL+coshkL式中，。通过文献 15 可知，无损悬臂梁的一阶模态为：kL=1.8751；二阶模态为kL=4.6941；三阶模态为 kL=7.8548。y(x,t)=Y(x)T(t)将代入式（8）中，左右两端同时乘 U(x)，并在 0L 上积分可得：m2T(t)t2+cT(t)t+kT(t)=F。(10)其中：m=AL0Y2(x)dx，(11)c=cL0Y2(x)dx，(12)k=ETIL04Y(x)x4Y(x)dx。(13)F=cos(t)L0p(xlp)Y(x)dx。(14)式中：m*为广义质量；c*为广义阻尼；k*为无裂纹梁的广义刚度，即无损梁刚度；F*为广义力。其中 k*用呼吸式裂
