本文考虑了大型汽轮机转子在高转速时,裂纹转子系统存在裂纹生长的情况,研究了陀螺效应对于大型汽轮机转子发展中裂纹的刚度的影响,并针对发展中的裂纹,建立了相应裂纹转子系统的三维陀螺力矩模型和三维刚度模型。进一步研究了质量偏心距、质量偏心角和陀螺力矩对裂纹转子的刚度的影响,得到了裂纹转子系统在三维陀螺力矩、质量偏心和偏心角耦合作用下的高转速转子在裂纹发展的情况下的刚度变化规律。 In this paper, the crack growth in the cracked rotor system of a large turbine rotor at high speed is considered, and this paper studies the influence of the gyro effect on the stiffness of a large steam turbine rotor with crack growth. For the developing cracks, the three-dimensional gyroscopic moment model and three-dimensional stiffness model of the corresponding cracked rotor system are formed. The coupling effects of mass eccentricity, unbalance angle and gyro moment on the stiffness of the heavy-duty rotor with damage evolution are studied respectively. The stiffness variation law of the high-speed rotor with three different types of cracks under the coupling of three- dimensional gyro moment, mass eccentricity and unbalance angle are obtained.
本文考虑了大型汽轮机转子在高转速时,裂纹转子系统存在裂纹生长的情况,研究了陀螺效应对于大型汽轮机转子发展中裂纹的刚度的影响,并针对发展中的裂纹,建立了相应裂纹转子系统的三维陀螺力矩模型和三维刚度模型。进一步研究了质量偏心距、质量偏心角和陀螺力矩对裂纹转子的刚度的影响,得到了裂纹转子系统在三维陀螺力矩、质量偏心和偏心角耦合作用下的高转速转子在裂纹发展的情况下的刚度变化规律。
转子,裂纹扩展,陀螺力矩,质量偏心角,刚度模型
Jiazhen Han, Huichun Peng
School of Mathematics and Physics, North China Electric Power University, Beijing
Received: Nov. 28th, 2020; accepted: Dec. 23rd, 2020; published: Dec. 30th, 2020
In this paper, the crack growth in the cracked rotor system of a large turbine rotor at high speed is considered, and this paper studies the influence of the gyro effect on the stiffness of a large steam turbine rotor with crack growth. For the developing cracks, the three-dimensional gyroscopic moment model and three-dimensional stiffness model of the corresponding cracked rotor system are formed. The coupling effects of mass eccentricity, unbalance angle and gyro moment on the stiffness of the heavy-duty rotor with damage evolution are studied respectively. The stiffness variation law of the high-speed rotor with three different types of cracks under the coupling of three- dimensional gyro moment, mass eccentricity and unbalance angle are obtained.
Keywords:Rotor, Crack Growth, Gyro Moment, Unbalance Angle, Stiffness Model
Copyright © 2020 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
涡轮机械转子通常在高温、高压及各种交变应力的作用下运行,导致了疲劳裂纹的产生,这可能会严重损害机器部件,甚至导致灾难性的故障。裂纹故障的研究已经有了很长的历史 [
刚度模型主要分为闭裂纹模型、开裂纹模型和呼吸裂纹模型 [
在研究转子横向运动时,旋转机械的陀螺效应就被考虑了进来。高转速是现代机械的一个发展趋势,对于大型重载(重力占优)的裂纹转子,在高转速的情况下,陀螺力矩对转子系统振动特性的影响也更加复杂 [
在本文的研究中,对大型重载(重力占优)的裂纹转子,考虑了高转速裂纹转子系统存在裂纹生长的情况,基于这三种类型的裂纹,建立了相应裂纹转子系统的陀螺力矩模型和刚度模型。考虑大型汽轮机转子的裂纹模型,并分别研究该模型对质量偏心角和陀螺力矩对裂纹转子的刚度的影响,最终得到了轴裂纹转子系统在三维陀螺力矩、质量偏心和偏心角耦合作用下的高转速转子在裂纹发展的情况下刚度变化规律。
当圆盘不装在两支承中间位置而偏向于一侧时,轴承变形后,圆盘轴线与两支点 A , B 的连线有一夹角 φ 。设圆盘的自转角速度为 Ω ,极转动惯量为 J p ,则圆盘对质心的动量矩为 H = J p Ω ,它与 A , B 的夹角也为 φ [
M g = J p Ω ω n sin φ (1)
这一惯性力矩称为陀螺力矩或回转力矩。由于 φ 角较小, sin φ ≈ φ ,上式可写作
M g = J p Ω ω n φ (2)
图1. 陀螺效应模型
设有如上图1的所示的转子,其转动角速度为 ω n ,极转动惯量为 J p ,直径转动惯量为 J d ,质量为m。建立如下模型:
{ Q ξ = m ε Ω 2 cos β + m g cos θ + J p Ω ω ξ Q η = m ε Ω 2 sin β + m g sin θ + J p Ω ω η Q ζ = J p Ω ω ζ (3)
{ ω ξ = θ ˙ ξ cos φ − θ ˙ y cos θ ξ sin φ ω η = − θ ˙ ξ sin φ + θ ˙ y cos θ ξ cos φ ω ζ = Ω − θ ˙ y sin θ ζ (4)
涡轮机械转子通常在高温、高压及各种交变应力的作用下运行,导致了疲劳裂纹的产生。本文研究的裂纹是受垂直于裂纹面的拉应力而产生的裂纹类型,因为其最危险容易引起低应力脆断。
本文只考虑转子产生的裂纹是张开型裂纹,由于 K ζ I = 0 ,故只考虑在 Q ξ , Q η 的作用,以上模型变为:
{ Q ξ = m ε Ω 2 cos β + m g cos θ + J p Ω ω ξ Q η = m ε Ω 2 sin β + m g sin θ + J p Ω ω η (5)
{ ω ξ = θ ˙ ξ cos φ − θ ˙ y cos θ ξ sin φ ω η = − θ ˙ η sin φ + θ ˙ y cos θ η cos φ (6)
由应力手册 [
{ σ ξ I ( ω ) = Q ξ L 4 α ′ 2 / I σ η I ( ω ) = Q η L 4 w / I (7)
其中, I = ( π D 4 / 64 ) , α ′ = D 2 − ( 2 w ) ,
总应力强度因子为
K I = K ξ I + K η I (8)
在这里,
{ K ξ I = σ ξ I π α F ( α / α ′ ) K η I = σ η I π α F ( α / α ′ ) (9)
其中
{ F ( α / α ′ ) = 2 α ′ π α tan ( π α 2 α ′ ) 0.923 + 0.199 [ 1 − sin ( π α / 2 α ′ ) ] 4 cos ( π α / 2 α ′ ) F ′ ( α / α ′ ) = 2 α ′ π α tan ( π α 2 α ′ ) 0.752 + 2.02 ( α / α ′ ) + 0.37 [ 1 − sin ( π α / 2 α ′ ) ] 3 cos ( π α / 2 α ′ ) F I I = 1.122 − 0.561 ( α / α ′ ) + 0.085 ( α / α ′ ) 2 + 0.81 ( α / α ′ ) 3 1 − α / α ′ F I I I = 2 α ′ π α tan ( π α 2 α ′ ) (10)
总的柔度计算公式
{ g 0 = L 3 48 E I g ξ ξ = g 0 + ∬ 128 L 2 α ′ 2 α E π D 8 F 2 ( α / α ′ ) d α d w g η η = g 0 + ∬ 512 L 2 w 2 α E π D 8 F 2 ( α / α ′ ) d α d w g ξ η = g η ξ = ∬ 256 L 2 α ′ 2 w E π D 8 α F ( α / α ′ ) F ′ ( α / α ′ ) d α d w (11)
L 3 / 48 E I 表示未开裂轴的灵活性。需要注意的是,上述表达式中裂纹深度的积分限为0到 α ,宽度的积分限没有指定。宽度积分的极限依赖于裂纹的开口部分,可以由裂纹边缘的应力强度因子得到。在找到裂纹边缘不同位置的总应力强度因子后,可以找出应力强度因子改变符号的位置。正应力强度因子表示张应力场和裂纹处于张开状态,负应力强度因子表示压缩应力场和裂纹处于闭合状态。一旦知道了裂纹张开部分的数量,就可以对裂纹张开部分进行积分。
利用上述柔度值,得到转动参照系刚度系数:
{ k 0 = 1 g 0 k ξ = g η g ξ g η − g ξ η 2 k η = g ξ g ξ g η − g ξ η 2 k ξ η = k η ξ = − g ξ η g ξ g η − g ξ η 2 (12)
上述刚度值随着转子的旋转而变化,这是由于受重力影响旋转时裂纹的开闭状态不同而导致的。由裂纹边缘的应力强度因子,找出裂纹的开口部分,进而得到新的柔度(式11),从而得到新的刚度值(式12)。
在式(12)中, k 0 为无裂纹转子的刚度,N/m; k ξ 、 k η 分别为 ξ 和 η 方向的刚度,N/m; k ξ η 、 k η ξ 为 ξ 和 η 两个方向相互作用的耦合柔度,N/m。通过矩阵变换,计算出定系下的刚度
[ k x k y x k x y k y ] = T − 1 [ k ξ k ξ η k η ξ k η ] T (13)
其中,矩阵T为
T = [ cos θ − sin θ sin θ cos θ ] (14)
由式(11)与式(12)即可求得转子在固定坐标系与旋转坐标系下的刚度。
本节考虑了裂纹的发展、轴的柔性以及盘的位置和质量偏心的不对称性,研究了陀螺力矩、质量偏心、不平衡角对转子系统时程刚度系数的耦合影响。对于时变刚度的计算,由于强度因子为零法有较高的计算精度且更加符合实际的试验情况,因而本文中采用此种方法基于Jeffcott转子进行计算,转子示意图如图2所示。
(1) 如图3,图4所示,相对刚度系数 k 1 = k ξ / k 0 , k 2 = k η / k 0 ,当值为1时裂纹的呼吸行为是完全封闭的,值小于1时裂纹的呼吸行为逐渐开放,裂纹扩展现象出现。各呼吸行为曲线的相对刚度系数到达最小值时裂纹几乎完全打开,此时 k 1 min = 0.941 , k 2 min = 0.953 。
图2. 转子示意图
(2) 取 ω ξ = 3 rad / s , ω η = 4 rad / s 时,改变偏心角的大小,所生成的相对刚度变化曲线如图3所示。在相同的低水平陀螺力矩下,质量偏心的方向对刚度的变化曲线有较大影响,且呈规律性变化。在一个周期内,当偏心角 β 由小增大时, ξ 轴方向和 η 轴方向上刚度的最值均在逐渐增大,范围则在逐渐缩小,说明在较低水平的陀螺力矩下,偏心角对于刚度的影响比较明显,此时陀螺效应对裂纹转子系统产生的影响较小。
图3. 陀螺力矩下不同偏心角的相对刚度
(3) 取 ω ξ = 9.7 rad / s , ω η = 11 rad / s 时,改变偏心角的大小,所生成的相对刚度变化曲线如图4所示。在相同的高水平陀螺力矩下,质量偏心角的大小对刚度的变化几乎没有影响。在一个周期内,当偏心角 β 由小增大时, ξ 轴方向和 η 轴方向上刚度的均无明显变化。因此,在高水平的陀螺力矩下,偏心角对于刚度的影响可忽略不计。
图4. 陀螺力矩下不同偏心角的相对刚度
(4) 综上,在陀螺力矩值偏小的情况下随着偏心角的增大,对裂纹转子系统的时变刚度影响逐渐增强。但当陀螺效应增大到一定程度时,偏心角的变化可忽略不计。此结论既符合工程实际,也在翟鹏程等 [
(1) 如4.1节的(1)所言,相对刚度系数 k 1 = k ξ / k 0 , k 2 = k η / k 0 ,当值为1时裂纹的呼吸行为是完全封闭的,值小于1时裂纹的呼吸行为逐渐开放,裂纹扩展现象出现。但在图5中能明显看出相对刚度系数值始终小于1,因此在偏心角为30。时,裂纹转子系统始终处于裂纹扩展的情况下,此时 k 1 , k 2 的最小值不变。
图5. 同一偏心角不同陀螺力矩下的相对刚度
(2) 取 β = 30 ∘ 时,改变陀螺力矩的大小,所生成的相对刚度变化曲线如图5所示。在相同的偏心角作用下,陀螺力矩对刚度的变化曲线有较大影响,且呈规律性变化。在一个周期内,当偏陀螺力矩由小增大时, ξ 轴方向和 η 轴方向上刚度的最值均在逐渐减小,范围则在逐渐增大,说明随着陀螺力矩的不断增大,刚度逐渐减小。
(3) 取 β = 150 ∘ 时,改变陀螺力矩的大小,所生成的相对刚度变化曲线如图6所示。对比图5,我们可以明显看出刚度曲线的变化,即当偏心角较大时,在相同的偏心角作用下,陀螺力矩对刚度的变化曲线有较大影响,且呈规律性变化。在一个周期内,当偏陀螺力矩由小增大时, ξ 轴方向和 η 轴方向上刚度的最值均在逐渐减小,范围逐渐增大,说明随着陀螺力矩的不断增大,刚度逐渐减小。
图6. 同一偏心角不同陀螺力矩下的相对刚度
(4) 综上,在偏心角不变的情况下,随着转速增加,陀螺效应逐渐增强,裂纹转子系统的时变刚度也逐渐减小,裂纹的呼吸行为也逐渐打开。早在蔡大文,胡腾等人 [
(1) 本文研究了陀螺效应对于大型汽轮机转子发展中轴裂纹的刚度的影响。对大型汽轮机转子,考虑了高转速裂纹转子系统存在裂纹生长的情况,并针对这种类型的裂纹,建立了相应裂纹转子系统的三维陀螺力矩模型和三维刚度模型。针对大型汽轮机转子的裂纹模型,研究了质量偏心距、质量偏心角和陀螺力矩对裂纹转子的刚度的影响,得到了裂纹转子系统在三维陀螺力矩、质量偏心和偏心角耦合作用下的高转速转子在裂纹发展的情况下的刚度变化规律。
(2) 当转速维持在一定范围内,陀螺力矩固定,偏心角对裂纹转子系统的刚度和稳定性都会产生一定的影响。当转速较低,陀螺力矩较小时,偏心角对刚度的影响较大。随着偏心角的增大,裂纹对转轴刚度的削弱程度逐渐降低,水平方向和竖直方向的影响效果并不相同。而当转速增加至一定值,陀螺力矩足够大时,偏心角的作用可以忽略不计。
(3) 同理,当偏心角保持不变,转速逐渐增加,陀螺效应逐渐增强。偏心角对裂纹转子系统的刚度和稳定性都会产生一定的影响。随着陀螺力矩的增大,刚度逐渐减小,裂纹扩展现象出现,直至达到一稳定值保持不变,裂纹的呼吸行为完全打开。
(4) 综上,我们在研究中得到了分别针对三种不同类型裂纹的转子系统在三维陀螺力矩、质量偏心和偏心角耦合作用下的高转速转子在裂纹发展的情况下刚度变化规律,这些研究结果能够为工程实践和故障转子系统失稳控制提供理论指导。
韩嘉桢,彭慧春. 大型汽轮机转子的发展中裂纹对刚度的影响Effect on the Stiffness of a Large Turbine Rotor with Developing Crack[J]. 力学研究, 2020, 09(04): 159-167. https://doi.org/10.12677/IJM.2020.94019