燃气轮机低压涡轮压气机转子动力学分析

本文研究的对象是该型燃气轮机低压涡轮压气机转子, 通过对转子-支承建模, 使用SAMCEF专业转子动力学分析软件, 采用有限元素法分析了其转子动力学特性, 包括转子的临界转速计算、稳态不平衡响应分析、转子稳定性分析等。验证了其在工程应用方面的可用性及可靠性, 同时得出了分析其动力学特性的基本方法及结论。

1　结构简介

该型燃气轮机低压涡轮压气机转子呈轴流轮毂式整体结构, 它由9级低压压气机和1级低压涡轮组成, 低压压气机与低压涡轮之间通过低压涡轮轴连接并以花键传递扭矩(见图1)。0 ～ 8级低压压气机轮盘、叶片材料为钛合金;低压涡轮盘、叶片材料为高温合金;低压涡轮轴材料为马氏体不锈钢。低压涡轮压气机转子呈3点支承结构, 前支承采用径向止推滚珠轴承, 支承点位于压气机0级轮盘前段, 由弹性支承、挤压油膜阻尼器、滚珠轴承组成;中间支承采用滚柱轴承, 由挤压油膜阻尼器、滚珠轴承组成, 支承点位于低压压气机后轴径后段;后支承采用滚柱轴承, 由弹性支承、挤压油膜阻尼器、滚珠轴承组成, 支承点位于低压涡轮轴后段。

2　计算模型

该型燃机低压涡轮压气机转子的3D模型较大,这将导致在网格划分以及计算过程中花费大量时间,因此我们对本机组的计算采用2D轴对称模型。在总体直角坐标系下建立二维轴对称单元, 其种类有3节点或高阶6 节点的三角形单元、4 节点或高阶8 节点的四边形单元。每个节点有9个自由度, 前6个自由度与梁单元一样, 分别为沿旋转轴线方向的拉伸和扭转, 以及由弯曲而引起的其他2个方向的线位移和角位移。另外, 3个自由度与旋转轴的横截面变形有关, 分别为拉伸引起的径向位移和弯曲引起的2个切向位移。同时, 使用这类单元可以很好地模拟转子的“涡动效应” 。

建立总体直角坐标系(X , Y , Z )和局部圆柱坐标系(er, z , eθ)。

设定ur, uz, uθ, r, z, θ为局部坐标系下的位移分量, 将它们沿θ向按傅立叶级数展开如下:

式中:n为波数;m为相数。每个谐函数(n, m)对应1个结构的基本变形。

低压涡轮压气机转子的压气机叶盘间的连接方式为焊接或螺栓紧固, 涡轮叶盘与涡轮轴之间也是以螺栓紧固的方式连接, 低压压气机与低压涡轮轴之间的连接为花键。因此可以认为压气机叶盘之间、涡轮盘与涡轮轴之间均为刚性连接;花键连接按铰链简化, 以此为依据进行整体建模。低压涡轮压气机转子叶片较长, 其质量及转动惯量都很大, 会产生较大的离心力和回转力矩, 这种回转效应会对转子的临界转速产生比较大的影响, 叶片等效简化为集中质量加到计算模型中。

3　计算内容

3.1　支承刚度和阻尼计算

现代航空发动机高速转子系统大都采用“柔轴”设计, 使发动机工作转速高于转子系统的临界转速,该型燃机低压涡轮压气机转子就采用了这种设计理念。由转子在临界转速时的挠度Y=meωc/c和传递率T≈ mωc/c可知, 此类转子工作在临界转速或者通过临界转速时, 其振幅Y正比于不平衡量(me)和转子的临界角速度(ωc), 反比于阻尼系数(c) 。显然, 对于“柔轴”类转子系统采取适当措施减小转子的不平衡量、降低临界转速和增大阻尼都是可以达到减小转子振动幅值的目的。而在发动机总体结构设计已经确定的前提下, 采用弹性支承或者阻尼器, 或者两者兼而有之, 就成为了最为有效和可行的方法。该型燃机低压涡轮压气机转子就采用了这种减振方式。机组所采用的弹性支承就是所谓低刚性弹性支承, 可以使转子在通过临界转速时的振幅和通过轴承座的外传载荷均很小, 这就可使转子工作在低阶临界转速以上。因此确定该型燃机低压压气机转子支承刚度是计算转子-支承系统临界转速所必需的前提条件。

3.1.1　前支承刚度和阻尼计算

低压压气机转子前支承位于0级轮盘前段, 由鼠笼式弹性支撑、挤压油膜阻尼器、AII76734BITI滚珠轴承组成。

1)AII76734BITI滚珠轴承的刚度

2)鼠笼刚度

以有限元法计算了鼠笼的静刚度。

3)挤压油膜阻尼器的刚度及阻尼计算

转子-支承系统采用挤压油膜阻尼器后, 便会出现油膜刚度和油膜阻尼, 在分析其动力特性时, 除转子-支承系统自身的刚度外, 还要计及油膜刚度。油膜阻尼远远大于转子-支承系统其他阻尼, 存在油膜阻尼时, 其他阻尼甚至可以忽略不计 。因此, 对于前支承来说, 油膜阻尼就是其支承系统的阻尼。依据其结构分析, 前支承的阻尼器应按短轴承半油膜计算。

阻尼器的8个参数计算

4)前支承总刚度、阻尼

3.1.2　中间支承刚度、阻尼计算

中间支承位于压气机后轴径后段, 由挤压油膜阻尼器、滚柱轴承组成。

1)滚柱轴承的刚度计算

2)中间支承总刚度及总阻尼计算

3.1.3　后支承刚度、阻尼计算

后支承位于低压涡轮轴后段, 由弹性支承、挤压油膜阻尼器、滚柱轴承组成(见图4)。

后支承总刚度及总阻尼计算如下：

3.2　低压涡轮压气机转子临界转速计算

低压涡轮压气机转子是由低压压气机转子、低压涡轮转子组成的轴系, 掌握轴系中单独转子的动力学特性对分析整体轴系的动力学特性是十分必要的。由于低压涡轮前端以花键与低压压气机后轴径相连接, 并依靠低压后轴径上的中间轴承作为其前支承点, 即低压涡轮并没有独立的前支承, 也就没有单独计算低压涡轮转子临界转速的必要。因此在计算整体轴系的临界转速之前, 首先对低压压气机转子进行了临界转速计算分析。

3.2.1　材料属性

低压涡轮气机转子是由几种不同的主体材料结构组成, 对2种不同材料的属性分别进行定义。低压压气机轮盘材料为钛合金:1.19e11 N/m, 0.33, 4 480kg/m3;低压压气机后轴径为马氏体不锈钢:2.01e11N/m, 0.3, 7 840 kg/m3 ;低压涡轮轴为马氏体不锈钢:2.01e11 N/m, 0.3, 7 800 kg/m3;低压涡轮为高温合金:2.01e11 N/m, 0.3, 8 320 kg/m3;分别计算低压压

气机叶片与低压涡轮叶片的等效质量和等效转动惯量后, 以集中质量的方式对其加载(见图5)。

3.2.2　计算结果及分析

1)计算结果

低压涡轮压气机转子临界转速计算结果如表1所示。1 ～ 3阶临界转速的振型如图6 ～图8所示。

　2)结果分析

由于低压涡轮压气机的慢车转速N1 =2 400r/min, 工作转速为N0 =7 700 r/min, 根据本文的计算结果, 可知

临界转速的安全系数:

超速时的临界转速的安全系数:

3.3　稳态谐波响应分析

稳态不平衡响应的计算是转子动力学分析中与临界转速计算同等重要的基本任务。稳态不平衡响应分析也可以用来确定系统的临界转速, 但它更重要的任务是不是用来求解在转子系统中可能存在的不平衡量作用下, 转子-支承系统的稳态不平衡响应,分析研究如何采取措施, 限制最大不平衡响应及减小不平衡响应。

在任一转速下的不平衡响应, 可能是若干个主振型的线性叠加, 随着转速的改变可由某一阶振型向另一阶振型过渡。不平衡量的分布规律的改变可引起稳态不平衡响应的变化, 不仅仅是量的变化, 也会引起相应“振型”的变化。因此, 稳态不平衡响应计算比临界转速计算能提供更多的有用信息。转子的激振力就是转子的转速, 振动幅值与转速的平方成比例。分析该转子的不平衡响应对于掌握其在不平衡量作用下的动态响应特性是非常重要的。

3.3.1　计算模型

按G1级精度确定低压涡轮转子的单位质量允许的剩余不平衡量。加重位置位于第0级压气机轮盘及涡轮盘上, 第0级压气机轮盘上的不平衡量为2.54 g;涡轮盘上的不平衡量为1.37 g。

3.3.2　计算结果

在不平衡响应计算结果中, 可以得到转子在计算频域范围内的最大响应值, 同时也可以得到工作转速范围内的最大响应值。最大响应值包括最大相对变形以及最大相对旋转角度。

3.4　瞬态响应分析

转子系统的瞬态响应分析主要是转子系统不平衡突然变化, 作用在转子系统上的外载荷突然变化或转子系统在变转速下工作等情况下, 转子系统的响应分析, 包括转子系统的位移、变形以及支承结构的传递载荷分布。对于带挤压油膜阻尼器的转子系统, 在机动载荷作用下的响应需进行瞬态分析。飞机机动飞行时引起的机动载荷(陀螺力矩和惯性力), 通常按静载处理, 即只计及机动载荷作用下转子结构的静位移。对于带挤压油膜阻尼器的转子系统, 需要用瞬态响应分析方法分析在该冲击载荷作用下阻尼器的承载能力和抑制失稳的能力。此外, 对于支承各向异性或带非同心型挤压油膜阻尼器的转子系统, 需要采用瞬态响应分析方法才可能求得稳态响应特性。

对于Jeffcott转子, 其计算模型为1只计刚性不计质量的柔性轴, 轴中央置1具有质量偏心的单盘,轴两端对称地安装在刚性支承上。盘的质心坐标为xc和yc, 它与固定坐标x和y的关系为:

　转子除受不平衡力作用外, 还受由加速或减速引起的惯性力作用, 在固定坐标系统内转子系统的运动方程为:

式中:m, c, k分别为盘的质量、系统外阻尼、轴刚度;θ=ω· t, ω为自转角速度, α为角加速度, e为质量偏心距。

4　结　语

1)临界转速计算结论

本文计算的该型燃机低压涡轮压气机转子的慢车转速的临界转速安全系数为1.53;全工况工作转速的临界转速安全系数为1.60;超速时的临界转速安全系数为1.39。可见, 此结构方案能满足低压涡轮压气机转子稳定运行的要求, 且有较高裕度。

2)不平衡响应计算结论

在工作频率范围内, 转子的动态响应较小, 最大变形为0.055 mm(62.5 Hz)因此, 低压涡轮压气机转子选取的平衡精度是合理的, 同时具有较小的振动幅值。

3)瞬态响应计算结论

分析了转子由0转速升至工作转速过程中, 第0级叶片丢失状态下转子的瞬态响应。由转子的轴心轨迹判断, 转子系统是稳定的, 已形成了稳态的圆轨迹;转子的振动幅值响应是在系统阻尼间隙范围内。

经以上分析证明, 该转子系统的设计是合理的,此转子结构方案能保证低压涡轮压气机的稳定运行。