18燃气轮机涡轮转子叶片设计及动力学分析-潘永生(11)

826第二十八届（2012）全国直升机年会论文燃气轮机涡轮转子叶片设计及动力学分析潘永生1杨库2李成2（1.海军驻哈尔滨地区军事代表室；2.中航工业哈飞，哈尔滨150066）摘要：燃气轮机涡轮转子叶片是燃气轮机的主要零部件之一，它的工作状态、工作环境、结构形式十分复杂。它的可靠性直接影响燃气轮机能否正常工作。有关统计资料指出，叶片的损坏事故中，绝大部分是由于振动引起的。因而，长期以来，叶片的振动一直是燃气轮机的一个普遍存在的严重问题，必须予以充分重视。由于外载荷、几何尺寸、材料性能等方面存在着不确定因素使振动具有随机性，从而影响叶片的可靠度。针对叶片工作中诸多参数的不确定性，如何进行动力学分析是结构设计和实际应用中一个十分重要的研究课题。关键词：燃气轮机；叶片；振动；动力学1引言燃气轮机广泛应用于航空、航天、船舶和电力等对国家安全和国民经济起着重要作用的领域及行业，尤其是水面舰船和航空飞行器的重要动力设备，对作战载体的战斗力起着举足轻重的作用，因此作战环境中燃气轮机的可靠性是保证战斗力的必要条件，而转子叶片是燃气轮机的关键零件，涡轮转子叶片在高温条件下高速旋转，承受离心力、气动力等循环交变载荷与动载荷作用，它的工作状态、工作环境、结构行为十分复杂，容易产生振动，统计资料表明，叶片的损坏事故中，绝大部分是由于振动引起的。因而，长期以来，叶片的振动一直是燃气轮机的一个普遍而又严重的问题，必须重视。所以对燃气轮机结构进行动力学分析，是十分必要的。2燃气轮机涡轮转子叶片力学分析燃气轮机是一种连续回转的内燃、叶轮机械式的热机，主要由压气机、燃烧室和涡轮三大部分构成。涡轮是受力最为复杂的部分。涡轮由转子和静子组成，转子是涡轮转动部分的总称，由涡轮盘、涡轮轴、工作叶片（转子叶片）及联接零件等组成，参见图1。转子叶片是燃气轮机中受力、受热最严重的零件之一，是进行功能转换的重要零件。图1涡轮转子示意图燃气轮机涡轮转子叶片由叶身、矩形缘板和燕尾形榫头三部分组成，见图2。叶身几何形状由叶背与叶盆表面上数据点的坐标，以流面形式给出。坐标系为发动机总体坐标系，x为转轴方向（向827前)，z轴沿径向通过叶根截面重心，y轴沿周向，xyz组成右手系。图2叶片图发动机工作时，燃气轮机涡轮转子叶片主要承受自身质量产生的离心力、气流的横向气体力、热负荷和振动负荷。本章将计算转子叶片的静强度可靠性。在计算涡轮转子叶片上的应力时，做如下假设：把转子叶片看作根部完全固装的悬臂梁。忽略叶片承受各种负荷后产生的变形。转子叶片仅承受自身质量离心力和横向气体力，只计算由离心力产生的拉伸应力和弯曲应力，以及由气体力引起的弯曲应力，不计叶片上的扭转应力、热应力和振动应力。设转子叶片各截面的扭转中心（刚心）、气体力压力中心与重心三者重合，则离心力和气体力均作用于重心，叶片受力变形的情况将大为简化。2.1叶片自身质量产生的离心力由于发动机的转速很高，故叶片自身的质量离心力很大。例如一个质量0.3kg的转子叶片在最大转速（11150r/min）工作时产生的离心力达7000多牛顿，为自重的23000多倍。叶片在其自身离心力的作用下，将产生很大的拉伸应力和弯曲应力，还能产生扭转应力。2.1.1离心力作用下的拉伸应力图3叶片微元体所产生的离心力如图3所示，在叶片上取一个高度为dZ的微元段，并在该微元段上取一个微元体dXdYdZ，令微元体的面积dAdXdY，则微元体的离心力dP为：828dAdZZdP2式中：——叶片材料的密度；ω——转子旋转角速度；Z——微元体重心到旋转轴X的距离。注意到cosZZ，则离心力dP沿Z轴方向的分量为：因此，截面积为A(Z)的叶片微元段质量所产生的沿Z轴方向的离心力为：ZdZZAdAZdZdPZAZAZ)(2)(2)(这样，就可求出转子叶片某一截面iZZ以上的叶片质量的离心力沿Z轴方向的分量为：kiZZiZdZZAP)(2,离式中：kZ——叶尖处的Z坐标值；A(Z)——叶片横截面面积，是随Z坐标变化的函数。则叶片某一截面iZZ上的离心拉伸应力为：)()(2,iZZiZAZdZZAki离（1-1）对于等截面叶片，)()(iZAZA常数，故式（3-1）变为222,21ikiZZ离（1-2）由式（1-2）可知，截面所在的坐标iZ越小，则该截面上的离心拉伸应力越大，叶根截面上的离心拉伸应力最大，而在叶尖处为零。2.1.2离心力弯矩由于转子叶片各截面重心的连线不与Z轴重合，所以叶片旋转时产生的离心力将引起离心弯矩。通常采用数值积分法计算离心力弯矩。如图4所示，将叶片分成n段共n+1个截面，假定叶片各截面的面积A1和重心坐标iiiZYX、、为已知，则第i段叶片的离心力分量为：图4离心力弯矩计算简图82922220sincosXiYiimiimiZiimiimippVRVYpVRVZ（1-3）式中：——叶片材料的密度；——旋转角速度；iV——第i段叶片的体积；miR——第i段叶片重心（坐标为mimiZY、、miX）到旋转轴的距离。可近似地认为：)(21)(21)(21))((2111111iimiiimiiimiiiiiiZZZYYYXXXZZAAV（1-4）叶片第j截面上的离心力弯矩应等于第1段到第j段叶片的离心力分量Yip、Zip对截面x、y轴的力矩的总和，即：jijmiziyj,jijijmiyijmizixj-XXPM-ZZP-YYPM1离11,)()()(离（1-5）利用式（1-3）、（1-4）、（1-5）进行计算，就可得到叶片各截面上的离心力弯矩。对于接近于等截面的涡轮工作叶片，通常只需求叶根截面的弯矩，此时将整个叶身作为一段进行近似计算，而不必分段计算。2.2叶片弯曲应力2.2.1气体力弯矩当气流流过叶栅通道时，气流的轴向和周向速度会发生变化，也就是说气流的动量发生变化，说明气流受到了力的作用。图5流过平均半径处的气流情况830如图5所示，根据动量定量，这两个轴向力之和应等于每秒钟内气流的动量变化量。在叶片平均半径处取宽为一个栅距，高为1单位长度的窗口，则流过该窗口的气流每秒钟内的动量变化为：ammammammammxctcctcP111122221)1()1(式中：amammmcc2121、、、——进、出口截面处气流的密度和轴向速度。设进、出口截面处的栅距相等，即122mmmmtttZQ，则：)(21212221ammammmxccQZP（1-6）式中：mZ——叶片的平均半径，2尖跟rrZm；Q——叶片数；mt——平均半径处的栅距。由叶栅进、出口截面处气流压差引起的轴向力为：)(21)(21212mmmmmmxpPQZtppP（1-7）式中：mmpp21、——分别是叶栅进、出口截面平均半径处气流的静压。将（1-6）、（1-7）式，即得叶片给予流过窗口气流的轴向力为：mmammammmppccQZ121212222（1-8）根据作用力与反作用力定律，单位叶高上受到的气体力轴向分量为：mmammammmxmppccQZp212221212（1-9）同理可得到单位叶高上受到的气体力周向分量为：)(2222111umammumammmymccccQZp（1-10）式中：umumcc21、——叶栅进、出口截面平均半径处气流的周向速度。2.2.2叶片弯曲应力的计算作用在转子叶片j截面上的总弯矩，也称为合成弯矩，应等于作用在该截面上的气体力弯矩和离心力弯矩的代数和，即：，离，气，合，离，气，合yjyjjxjxjxjMMMMMMy（1-11）按式（1-11）求出对转子叶片某截面上X，Y轴的合成弯矩后，还不能立即计算弯曲应力。因为X，Y轴通常不是截面的主惯性轴。如图6所示，叶片截面与一个狭长的长方形截面相仿。这类截面的抗弯能力在各个方向是不一样的。通过截面重心，有一对惯性轴、，对于轴的惯性距最小，对轴的惯性矩最大。计算弯曲应力时，总是对惯性主轴、而言的。因此，必须把对X，Y轴的合成弯矩合xM、合yM转换为对，轴的合成弯矩。831图6叶片截面的主惯性轴图7合成弯矩的转换通过叶片截面重心的最小惯性轴近似地平行于叶弦，其正向指向叶片后缘，轴与y轴的夹角为。合成弯矩合xM、合yM在，轴上的投影之代数和，即为合,M、合,M，转换公式为（参见图7）：αMαMMαMαMMy,x,ξ,y,x,η,cossinsincos合合合合合合(1-12)根据叠加定理，这三个点在合,M、合,M作用下的弯曲应力计算式为：CηCξCBηBξBAηAAJMJMJMJJMJ合合弯合合弯合合，弯,,,,,,,,MM(1-13)式中：A、A、B、B、c、C——A,B,C三点的坐标，必须带有正负号；JJ、——叶片截面的主惯性矩。8323燃气轮机涡轮转子叶片动力学分析3.1叶片的自振频率影响因素叶片的自振频率（固有频率）是叶片本身的固有特性，叶片尺寸已定，其自振频率便已确定，但是随同叶片的旋转及所受温度环境和根部安装形式等不同影响，叶片的自振频率将有所改变。3.1.1离心力场的影响涡轮的工作叶片处在高速旋转状态下工作，转动使叶片产生很大的质量离心力。这个力通常以一个拉力的方式作用在叶身上，迫使振动的叶片有恢复到原平衡位置趋势。这相当增强了叶片的弹性恢复力，结果是叶片的各阶固有频率有所提高。因此，旋转着的叶片的固有频率随转速的增大而增高。通常转动下的叶片的自振频率称为动频，以Df表示。不转动时叶片的自振频率称为静频，以f表示。3.1.2温度的影响涡轮叶片处于高温状态下工作，由于温度升高，材料的弹性模量E要下降，叶片的各阶频率都会随E而改变。一般应用下式计算温度对频率的影响：00ttttEEff（1-15）式中：tf——叶片高温下的自振频率；0tf——片在常温下的自振频率；tE——高温下叶片材料的弹性模量；0tE——常温下叶片材料的弹性模量。涡轮叶片温度通常使频率降低，而转速又使频率增高，有时两者升降率相差不多，使叶片的动频没有多大变化。3.1.3一般研究叶片的振动特性都认为叶片根部与轮盘的联结是绝对刚性固装的。但大多数叶片榫头与盘槽的配合是留有少量间隙的，目的是为了消除工作时的热应力。并且叶片是弹性体，在两个零件的联结处也都有一定的弹性。严格地说，这种联结结构不能看作固装形式。因为随转速增加，叶身的巨大离心力作用在榫头上，结合面上产生巨大的压力，从而逐渐使叶片失去其在间隙内活动的可能，根部联结处的刚性随转速的增高而加大。最终导致成“完全固装”条件。图8中给出了具有纵树形榫头联结的涡轮叶片的振动频率随侧向力P的变化曲线。侧向力相当于叶片旋转时，叶片离心力负荷加给榫头上的紧固力。可见，在低转速范围，紧固力随转速变化使叶根刚性改变较大，振动频率增加较快。转速较高时，紧固力使叶根刚性改变较慢，振动频率变化也较缓慢。就是说叶片在低转速下，振动频率随叶根部刚性变化的影响较大，而高转速时则影响较小。833图8纵树形榫头频率随紧固力P的变化曲线实际上，当转速达到最大值的某一百分比时，就可以看成“完全固装”了。3.2叶片动力学分析及算例本节给出动力涡轮级数为4，研究的转子叶片为二级转子叶片，其前后的每级静叶片数都为68，因此尾流激振力的激振频率中，构造系数K=68或K=136。考虑四根支柱，K=4。本文利用ANSYS软件对叶片固有频率进行分析。得出固有频率均值。由于涡轮转子结构是由具有旋转周期对称性的N个重复扇区组成的。即绕其转轴转动N2(N为叶片数，本例中为71)角度后，结构的几何形状和转动前完全一样。根据叶片与轮盘结构的旋转对称