不同温度下含裂纹压气机叶片的模态特征和振动特性分析_周震霆.pdf

不同温度下含裂纹压气机叶片的模态特征和振动特性分析不同温度下含裂纹压气机叶片的模态特征和振动特性分析周震霆，贺星，刘永葆(海军工程大学动力工程学院，湖北武汉430033)摘 要:本文选取含裂纹悬臂梁模拟含裂纹压气机叶片，基于无质量扭转弹簧和呼吸式裂纹刚度模型，得到一种不同环境温度下含裂纹压气机叶片的模态特征和振动特性的分析方法。通过弹性模量引入温度模块，利用无质量扭转弹簧连续条件得到关于含裂纹梁的特征方程，分析环境温度和裂纹深度对含裂纹梁固有频率的影响；利用悬臂梁的强迫弯曲振动方程，引入呼吸式裂纹刚度模型，改变激振力频率，分析环境温度及激振频率对含裂纹梁振动位移响应的影响。结果表明，环境温度越高，含裂纹梁的固有频率越小，梁的振动位移响应越大。同时，激振频率的选取也具有一定影响。关键词：压气机叶片；呼吸式裂纹；温度；模态特征；振动特性中图分类号：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-bao(SchoolofPowerEngineering,NavalUniversityofEngineering,Wuhan430033,China)Abstract:Inthispaper,acrackedcantileverbeamisselectedtosimulateacrackedcompressorblade.Basedonamass-lesstorsionspringandabreathingcrackstiffnessmodel,amethodforanalyzingthemodalcharacteristicsandvibrationchar-acteristicsofacrackedcompressorbladeatdifferentambienttemperaturesisobtained.Theelasticmodulusisintroducedin-tothetemperaturemodule,andthecharacteristicequationaboutthecrackedbeamisobtainedbyusingthecontinuouscondi-tionofthemasslesstorsionspring,andtheinfluenceoftheambienttemperatureandcrackdepthonthenaturalfrequencyofthecrackedbeamisanalyzed;theforcedbendingvibrationequationofthecantileverbeamisusedtointroducethebreathingThecrackstiffnessmodelwasusedtochangetheexcitationforcefrequencytoanalyzetheinfluenceofambienttemperatureandexcitationfrequencyonthevibrationdisplacementresponseofthebeamwithcracks.Theresultsshowthatthehighertheambienttemperature,thesmallerthenaturalfrequencyofthebeamwithcracks,andthegreaterthevibrationdisplacementresponseofthebeam;atthesametime,theselectionoftheexcitationfrequencyalsohasacertaininfluence.Key words:compressorblade；breathingcrack；temperature；modalcharacteristics；vibrationcharacteristics0引言压气机叶片是燃气轮机运行过程中完成能量转换的基本单元，也是燃气轮机关键部件中出现故障的最典型部件之一。裂纹是压气机叶片出现故障的常见形式之一，裂纹的存在将会改变叶片的模态特征和振动特性。燃气轮机压气机叶片工作环境恶劣，承受着离心载荷和气动载荷等复杂载荷的影响，特别是高压压气机，随着压比的升高，高压压气机周围气流的温度也不断升高，同时燃烧室通过传热等方式也会导致高压压气机温度升高。在高温的影响下，压气机叶片的材料和力学性能会发生变化，也会对压气机叶片的模第 45卷第12期舰船科学技术Vol.45,No.122023年6月SHIPSCIENCEANDTECHNOLOGYJun.,2023收稿日期:20220507基金项目:国家科技重大专项作者简介:周震霆(1998)，男，硕士研究生，主要研究方向为动力机械及热力系统的设计、仿真与优化。态和振动特性产生影响。压气机叶片可被看作悬臂梁，通过裂纹悬臂梁模型模拟含裂纹的压气机叶片。对于含裂纹梁的模态特征和振动特性，20 世纪便有人开始进行了研究，并发展至今19，可见研究含裂纹梁的模态特征和振动特性具有实际的工程意义。分析裂纹深度、裂纹位置、结构阻尼、刚度等因素对梁的模态特征和振动特性的影响，却很少考虑温度的影响1012。然而在实际过程中，空气流过各级叶片时，由于受到压缩，压力和温度会逐级升高，位于后面级的压气机，温度产生的影响很难忽略，尤其是在迫切渴求实时监测和在线诊断的今天。若想在燃气轮机工作过程中得到压气机叶片的状态信息，其所处工况的环境温度应考虑在内。本文通过含裂纹悬臂梁模型模拟含裂纹的压气机叶片，采用无质量扭转弹簧等效代替裂纹。通过弹性模量引入温度模块，推导出不同温度下含裂纹悬臂梁的模态特征方程，并通过梁的受迫弯曲振动方程，得到温度影响下梁的振动特性曲线。1模型建立和模态分析 1.1 含裂纹悬臂梁模型建立如图 1(a)所示，建立一段长为 L，宽为 b，高为h 的等截面矩形悬臂梁模拟裂纹叶片，并忽略裂纹尺寸随温度的变化。在该悬臂梁表面有一条深度为 a 的裂纹，该裂纹距离固定端的距离为 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 为悬臂梁的长度，为裂纹的相对位置，为 泊 松 比，为 不 同 温 度 下 的 弹 性 模 量；为相对裂纹深度。通过文献 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)假定有分离变量解存在，令，代入式（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舰船科学技术第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(kR)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)cosh(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)(xlp)。(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*用呼吸式裂纹刚度 kbr替代，从而得到呼吸式裂纹梁的弯曲第45卷周震霆，等：不同温度下含裂纹压气机叶片的模态特征和振动特性95振动方程。对于呼吸式裂纹刚度 kbr，刘文光等14提出：kbr=ko+kko2(1+cos(t)。(15)其中：ko=kTkkT+k。(16)式中，ko为张开式裂纹梁模型。当 t=2n(n=1,2,)时，裂纹完全闭合，此时呼吸梁模型等效为闭合梁模型，即 kbr=k*；当 t=(2n-1)，(n=1,2,)时，裂纹完全张开，此时呼吸梁模型等效为张开梁模型，即kbr=ko。将式（10）的 k*用式（16）中的呼吸式裂纹刚度kbr替代，从而得到呼吸式裂纹梁的弯曲振动方程：m2T(t)t2+cT(t)t+kbrT(t)=F。(17)2.2 振动特性分析由式（13）可知，悬臂梁的广义刚度与弹性模量具有一定的关系，因此可知呼吸式裂纹刚度与弹性模量也具有关系。温度的变化通过对弹性模量的影响，导致呼吸式裂纹刚度发生变化，进而会对梁的振动位移响应造成影响。通过对梁的振动位移及其他振动特性与温度的变化之间的规律进行分析，得出相关结论。3结果与讨论以图 1 所示悬臂梁为例模拟含裂纹叶片，假设叶片的几何尺寸为：长度 L=0.09m，宽度 b=0.055m，高度 h=0.006m；材料参数为：密度=4480kg.m3，泊松比 v=0.33，阻尼比=0.05，弹性模量 ET随温度发生变化。3.1 环境温度对裂纹梁模态特征的影响112212312412512612712结合假设的含裂纹叶片的几何尺寸，取裂纹的相对位置为R/L=1/9，相对深度为a/h,环境温度的变化范围为 T20100200300400500600。通过数值计算得到环境温度和不同相对裂纹深度与含裂纹梁的第 1 阶固有频率之间的关系如图 3 所示。可以看出，在悬臂梁裂纹位置和深度一定的情况下，当环境温度不断增大时，悬臂梁的固有频率不断减少。通过比较，可知悬臂梁在环境温度为 600下的固有频率比在 20 下的降低 70Hz 左右，大约为悬臂梁固有频率的 12%左右。由此可见，环境温度对悬臂梁固有频率具有一定的影响。而当温度不变的情况下，梁的相对裂纹深度越大，固有频率越低。同时也可以看出，当梁的相对裂纹位置不变的情况下，环境温度对梁的一阶固有频率的影响明显比相对裂纹深度大。3.2 环境温度对裂纹梁振动特性的影响在假设的含裂纹叶片几何尺寸以及材料参数基础上，裂纹位置 R=0.01m，裂纹深度为 a=0.0005m，引入外部简谐激振力 F=pcos(t)，其中激振力幅值为100N，环境温度的变化为 T20300600，使梁分别在 400Hz、固有频率及 4000Hz 下振动，通过数值计算得到环境温度与含裂纹梁的振动特性之间的关系如图 4图 6 所示。由图 4 和图 5 中的位移响应图可以看出，在选取的低频和固有频率作用下，温度越高，梁的振动幅值越大。这是因为随着温度的升高，梁的刚度和阻尼会出现减小的情况，在相同激振力幅值的影响下，刚度和阻尼越小，梁的位移响应越大。从图 4 和图 5 的相图也可以看到相关现象，反映出选取的呼吸式裂纹刚度所作用的效果。同时，从图 5（b）可以看到在共振情况下阻尼产生明显的响应滞后效应。在选取的高频作用下，为温度越高，梁的振动幅值越大，与低频作用下的现象不同，如图 6 所示。随着温度的升高，梁的刚度和阻尼减小，梁的位移响应增大。在相同阶数弯曲频率下，随着温度的升高，梁的固有频率下降，且下降幅度较大。在每一阶频率附近都存在一个幅值峰，各阶频率构成了连续的幅值峰。当梁处于 300 和 600 时，4000Hz 可能位于该温度下某一高阶幅值峰的中段，而位于 20 下某一高阶与临近低阶间的谷底。4结语本文通过含裂纹悬臂梁模拟含裂纹压气机叶片，选取无质量扭转弹簧模拟裂纹，并引入呼吸式裂纹刚度模拟裂纹振动过程中的呼吸效应。通过弹性模量引图3固有频率对环境温度变化曲线Fig.3Thenaturalfrequencyvs.ambienttemperature96舰船科学技术第45卷入温度模块，研究环境温度对含裂纹压气机叶片的一阶固有频率和振动响应之间的关系，得出以下结论：1）环境温度对梁的一阶固有频率具有一定的影响，随着环境温度的升高，梁的一阶固有频率下降，而且下降的幅度较裂纹深度带来的影响更大。2）在低频及固有频率激振力作用下，随着环境温图4400Hz 频率下振动特性图Fig.4Thevibrationperformancewith400Hzfrequency图5固有频率下振动特性图Fig.5Thevibrationperformancewithnaturalfrequency第45卷周震霆，等：不同温度下含裂纹压气机叶片的模态特征和振动特性97度的升高，梁的位移响应越大，但在高频激振力作用下，出现了不同的现象。这是所选激振频率处于不同温度下各阶频率的距离远近导致的。由此可见，环境温度对叶片的模态特征和振动特性具有一定影响。若想在实际工作过程中实现实时监测和在线诊断，应选取合适的工作频率，需要后续进行研究得到。参考文献：OSTACHOWICZWM,KRAWCZUKM.Analysisoftheef-fectofcrackonthenaturalfrequenciesofacantileverbeamJ.JournalofSoundandVibration,1991,150(2):191201.1王术新,姜哲.裂缝悬臂梁的振动特性及其裂缝参数识别J.振动与冲击,2003,22(3):8385.2QIANGL,GUS.N,JIANGJ.S,ThedynamicbehaviorandcrackdetectionofabeamwithacrackJ.JournalofSoundandVibration,1990,138(2):233-243.3董广明,陈进,邹剑.裂纹转子的振动特性与裂纹参数识别 J.上海交通大学学报,2004,38(6):857861.4CHENLH,DUANJW,SUNY,etal.Thestudyofthevibra-tioncharacteristicsofthecantileverbeamwithasurfaceCJ.AppliedMechanicsandMaterials,2013(394):121127.5杨海燕,杨秉玉,刘启州.疲劳裂纹叶片振动的非线性特性研究 J.西北工业大学学报,1999,17(2):198204.6胡家顺,冯新,周晶.呼吸裂纹梁非线性动力特性研究 J.振动与冲击,2009,28(1):7680,87.7刘文光,陈国平.含裂纹悬臂梁的振动与疲劳耦合分析 J.振动与冲击,2011,30(5):140144.8XIEJ,ZIY,ZHANGM,etal.AnovelvibrationmodelingmethodforarotatingbladewithbreathingcracksJ.ScienceChina:TechnologySciences,2019,62(2):333348.9马一江,陈国平.不同温度下裂纹悬臂梁的模态和疲劳寿命分析 J.振动与冲击,2017,36(8):132137+163.10孙强,张忠平,柴桥,等.航空发动机压气机叶片振动频率与温度的关系 J.应用力学学报,2004,21(4):137139.11王振清,刘兵,韩玉来.高温下含裂纹铝合金梁自由振动频率分析 J.哈尔滨工程大学学报,2012,33(3):320324.12中国标准出版社.中国航空材料手册钛合金、铜合金 M.2001,2(3).13刘文光,陈国平.呼吸式裂纹梁的振动疲劳裂纹扩展耦合分析 J.中国机械工程,2010,21(23):27982802.14辛格雷苏.S.拉奥.机械振动 M.李欣业,杨理诚译.北京:清华大学出版社,2016.15图64000Hz 下振动特性图Fig.6Thevibrationperformancewith4000Hzfrequency98舰船科学技术第45卷