CN113847981B - 一种基于机械特性的水电机组保护性振动监测方法 - Google Patents
一种基于机械特性的水电机组保护性振动监测方法 Download PDFInfo
- Publication number
- CN113847981B CN113847981B CN202111089115.9A CN202111089115A CN113847981B CN 113847981 B CN113847981 B CN 113847981B CN 202111089115 A CN202111089115 A CN 202111089115A CN 113847981 B CN113847981 B CN 113847981B
- Authority
- CN
- China
- Prior art keywords
- bearing
- load
- generating set
- hydroelectric generating
- axis
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000012544 monitoring process Methods 0.000 title claims abstract description 24
- 230000001681 protective effect Effects 0.000 title claims abstract description 10
- 238000013461 design Methods 0.000 claims abstract description 15
- 238000006073 displacement reaction Methods 0.000 claims description 46
- 238000013016 damping Methods 0.000 claims description 23
- 230000003068 static effect Effects 0.000 claims description 16
- 230000033001 locomotion Effects 0.000 claims description 13
- 238000001228 spectrum Methods 0.000 claims description 10
- 230000003993 interaction Effects 0.000 claims description 4
- 230000001360 synchronised effect Effects 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 4
- 238000005303 weighing Methods 0.000 claims description 2
- 230000000704 physical effect Effects 0.000 abstract description 3
- 230000008569 process Effects 0.000 abstract description 3
- 238000005259 measurement Methods 0.000 description 10
- 230000001133 acceleration Effects 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 6
- 238000006880 cross-coupling reaction Methods 0.000 description 5
- 230000000670 limiting effect Effects 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 239000013307 optical fiber Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 229910000831 Steel Inorganic materials 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000000314 lubricant Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 238000004513 sizing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/04—Bearings
- G01M13/045—Acoustic or vibration analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M99/00—Subject matter not provided for in other groups of this subclass
- G01M99/005—Testing of complete machines, e.g. washing-machines or mobile phones
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M99/00—Subject matter not provided for in other groups of this subclass
- G01M99/007—Subject matter not provided for in other groups of this subclass by applying a load, e.g. for resistance or wear testing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Hydraulic Turbines (AREA)
Abstract
一种基于机械特性的水电机组保护性振动监测方法，包括以下步骤：步骤1，构建水电机组径向力和振动模型；步骤2，获取水电机组设计标准中载荷约束；步骤3，使用步骤1中的模型和步骤2中获得的约束，设置警报和跳闸水平。与现有技术相比，根据水电机组的物理性质，即根据水电机组的设计标准和运行过程中测量的径向负荷，提出振动监测方法。有关状态监测的方法基于关键部件的机械特性和作用在这些部件上的测量径向载荷。使用这种方法，可以根据与正常运行期间预期载荷水平相关的径向载荷水平以及关键部件的疲劳极限来确定监控系统的警报和跳闸水平。
Description
技术领域
本发明属于水电机组保护领域，更具体地，涉及一种基于机械特性的水电机组保护性振动监测方法。
背景技术
现有技术文件1(CN104330149A)公开了一种水电站机组振动测量监测方法，压电式加速度传感器，置于水电站机组机架或顶盖上；前置处理模块，置于压电式加速度传感器后端，将压电式加速度传感器的数据信号经过两次硬件积分为振动位移信号，最终输出到振摆保护数据采集装置完成A/D数据采集，同时负责抑制传感器共振信号；压电式加速度传感器和前置处理模块之间采用低噪声、低衰减电缆，以保证电信号在传输过程中不带来衰减和更大的噪声。
现有技术文件2(CN105628173B)公开了一种基于光纤传感的水电站机组振动测量监测装置，包括光纤光栅加速度传感器组、1×N光纤耦合器、环形器、M个前置波分复用器、后置波分复用器、宽带光源、光电转换滤波模块、非平衡迈克尔逊干涉仪解调模块、现场可编程门阵列模块、净化电源箱、微处理器工控主板模块。其中采用光纤光栅加速度传感器安装在水电机组的顶盖或机架上，通过设置固定差别长度的铠装光缆连接进行光信号传输，然后对光信号干涉法解调。
包括ISO标准在内，例如，ISO-7919-5标准中关于水电机组中轴和轴承座径向运动的现有建议，ISO 10816-5：机械振动——通过测量非机械振动来评估机器振动，振动监测为基于统计数据的警报级别，不考虑机器的机械特性。根据这些标准，水电机组的同步速度决定了推荐的最大轴位移和外壳加速度。然而，为了保护机器免受有害的径向载荷，监测水电机组中的径向载荷非常重要。
发明内容
为解决现有技术中存在的不足，本发明的目的在于，提供一种基于机械特性的水电机组保护性振动监测方法。
本发明采用如下的技术方案。一种基于机械特性的水电机组保护性振动监测方法，其特征在于，包括以下步骤：
步骤1，构建水电机组径向力和振动模型；
步骤2，获取水电机组设计标准中载荷约束；
步骤3，使用步骤1中的模型和步骤2中获得的约束，设置警报和跳闸水平。
优选地，步骤1进一步包括：
步骤1.1，构建轴承与支架之间的相互作用模型；
步骤1.2，测量或计算获得水电机组径向力；
步骤1.3，使用轴承特性和轴运动确定水电机组中导向轴承中的径向载荷；
步骤1.4，实测轴承载荷和根据轴位移和轴承特性计算的轴承载荷。
优选地，步骤1.1中，
轴承油膜特性的力平衡模型以如下公式表示，
式中：
f表示由转子引起作用在轴承上的径向载荷，fx表示在x轴的分量，fy表示在y轴的分量；
k表示轴承的刚度，kxx表示其x轴的分量，kyy表示其y轴的分量，kxy，kyx表示其交叉耦合刚度，
c表示轴承的阻尼，cxx表示其x轴的分量，cyy表示其y轴的分量，cxy，cyx表示其交叉耦合阻尼，
uSx表示轴的绝对位移矢量us的x轴分量，uSy表示y轴分量；
uHx表示轴承箱的绝对位移矢量uH的x轴分量，uHy表示y轴分量。
优选地，位移与组合刚度和阻尼特性之间的关系以如下公式表示，
D＝[I+ZH-1]-1 (4)
式中：
D表示轴承和支架的组合特性，D＝Kc+iΩCc；
Kc是组合刚度，Cc是组合阻尼。
优选地，步骤1.3中，通过计算轴在x和y方向的偏心距，轴的总偏心距εt和相位θ以如下公式表示，
式中：
cbx表示径向轴承游隙x轴的分量；
cby表示径向轴承游隙y轴的分量。
优选地，步骤2进一步包括：
步骤2.1，根据水电机组组件的载荷谱获得水电机组载荷约束；
步骤2.2，根据平衡等级获得载荷约束；
步骤2.3，根据发电机偏心率和转轮不均匀流动特性获得载荷约束。
优选地，步骤2.3中，弗朗西斯水轮机的模型以如下公式表示，
式中：
fr表示转轮的径向载荷；
P表示额定功率；
Ω表示水轮机组同步转速；
D是转轮直径。
优选地，步骤2.1中，循环数和事件数为载荷谱提供输入数据，用作单元疲劳设计的支持。
优选地，步骤3进一步包括：
步骤3.1，识别水电机组中机械部件的特性；
步骤3.2，振动与载荷水平之间的关系；
步骤3.3，确定正常运行期间的最大载荷水平；
步骤3.4，基于运行模式、平衡等级和疲劳极限设置报警和跳闸阈值。
优选地，步骤3.2进一步包括：
(a)在轴承内部或后面安装称重传感器，测量静态和动态载荷；
(b)根据计算出的轴承特性和测得的轴位移计算轴承载荷；
(c)根据轴承支架的刚度特性和测得的轴承座运动计算轴承载荷
本发明的有益效果在于，与现有技术相比，根据水电机组的物理性质，即根据水电机组的设计标准和运行过程中测量的径向负荷，提出振动监测方法。有关状态监测的方法基于关键部件的机械特性和作用在这些部件上的测量径向载荷。使用这种方法，可以根据与正常运行期间预期载荷水平相关的径向载荷水平以及关键部件的疲劳极限来确定监控系统的警报和跳闸水平。本发明还提供了一种使用计算出的轴承参数和轴位移测量值来测量水电机组轴承载荷的替代方法。这种方法预计不会像使用称重传感器的载荷测量方法具有相同的高分辨率，但该方法所需的安装快速、简单，并且不需要对水电机组中的组件进行任何修改。
附图说明
图1为水电机组轴承和支架示意图；
图2为轴承和支架中的总阻尼和运动分布图；
图3为水电机组可倾瓦轴承的刚度与偏心率的关系图；
图4为水电机组可倾瓦轴承的阻尼特性与偏心率的关系图；
图5为作为偏心函数的轴承载荷曲线图；
图6为水力发电机组可倾瓦轴承间隙和偏心距的函数计算轴承载荷示意图。
具体实施方式
下面结合附图对本申请作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本申请的保护范围。
本发明提供了一种基于机械特性的水电机组保护性振动监测方法，包括以下步骤：
步骤1，构建水电机组径向力和振动模型，具体包括：
在立式水电机组中，影响径向轴承和支架的径向力理想情况下为零。作用在轴承段上的压力在轴承箱内部平衡，进而不会影响周围的结构。但在工程实践当中，理想情况并不存在。所有水电机组或多或少都会受到几何偏差、转子平衡、定子和转子偏心度、水轮机中的动态流动特性等的影响。
步骤1.1，构建轴承与支架之间的相互作用模型。
通常，使用安装在与刚性支架相连的轴承箱上的加速度计进行振动测量，并且忽略轴承和支架的刚度特性，这种情况下轴的振动水平将被严重低估。水电机组的轴承具备刚性，最大的位移将在支架中，如果使用位移传感器测量轴和轴承箱之间的距离将严重低估轴的移动。因此，本发明的一个创新点在于，组合刚度和阻尼特性，而不仅仅是轴承或支架特性。
如图1所示，建立轴承与支架之间的相互作用模型，在垂直于轴的平面，建立坐标系xOy，坐标原点O位于静态状态下的轴承几何中心，以kij表示轴承的刚度参数，以cij表示轴承的阻尼参数，以hij表示支架的刚度参数，i，j＝x，y；以u表示位移矢量，是指轴承中的轴位移，即轴相对于轴承中心的位移，位移矢量u指向轴承的几何中心以频率Ω旋转，以如下公式表示轴承油膜特性的力平衡模型表示为，
式中：
u表示位移矢量，ux表示其x轴的分量，uy表示其y轴的分量，u＝uS-uH，uS表示轴的绝对位移矢量，uH表示轴承箱的绝对位移矢量；
f表示由转子引起作用在轴承上的径向载荷，fx表示在x轴的分量，fy表示在y轴的分量；
k表示轴承的刚度，kxx表示其x轴的分量，kyy表示其y轴的分量，kxy，kyx表示其交叉耦合刚度，
c表示轴承的阻尼，cxx表示其x轴的分量，cyy表示其y轴的分量，cxy，cyx表示其交叉耦合阻尼，
m表示轴承的质量，mxx表示其x轴的分量，myy表示其y轴的分量，mxy，myx表示其交叉耦合质量。
在大多数传统轴承模型中，流体惯性力可以被忽略。对于旋转频率为Ω的纯谐波运动，轴承油膜特性的力平衡模型进一步可以表示为如下公式，
式中：
uSx表示轴的绝对位移矢量uS的x轴分量，uSy表示y轴分量；
uHx表示轴承箱的绝对位移矢量uH的x轴分量，uHy表示y轴分量；
作用在轴承上的径向载荷f也作用在轴承座上，轴承油膜特性的力平衡模型进一步可以表示为如下公式，
式中：
uH表示轴承箱的绝对位移矢量，
H表示支架的刚度，
hij表示支架的刚度参数，i，j＝x，y；
轴承支架中的交叉耦合阻尼和刚度被忽略，支架由钢梁制成，位移与组合刚度和阻尼特性之间的关系可以简写为，
f＝HuH
uH＝H-1f
f＝Z(uS-uH)＝ZuS-ZuH＝ZuS-ZH-1f
f＝[I+ZH-1]-1ZuS＝DuS
D＝[I+ZH-1]-1 (4)
式中：
D表示轴承和支架的组合特性，D＝Kc+iΩCc；
Kc是组合刚度，Cc是组合阻尼。
轴承和支架的刚度均为线性，则轴的位移将分布在轴承和支架之间，如图2所示，组合阻尼相对于轴承阻尼的变化取决于轴承和支架之间的刚度关系，其中，横坐标表示轴承与支架之间的刚度关系，水电机组中，轴承支架的刚度在0.2GN/m到4GN/m之间变化，具体取决于支架内外部分的连接方式以及钢的尺寸结构体。
水电机组在正常运行时的轴承特性为kxx，kyy～0.4-2GN/m，交叉耦合刚度项kxy和kyx明显小于可倾瓦轴承的kxx和kyy。在图2所示的实例中，轴承中的阻尼cxx和cyy设置为0.1GNs/m。这些轴承特性的假定值很大程度上取决于轴承游隙和运行条件，即轴承载荷。图2中的函数显示了轴在轴承中的运动和轴承支架运动之间轴的总位移分布，即分别为usj/uj和uhj/uj。图3还显示了总阻尼与轴承阻尼之间的关系，即ccjj/cjj。上述模型非常简单，并没有考虑到轴颈轴承是非线性的，并且刚度变化在很大程度上取决于计算刚度的偏心距这一事实。这表明在不知道轴承和支架的特性的情况下评估从位移传感器或加速度计收集的振动数据是很困难的。
步骤1.2，测量或计算获得水电机组径向力。
水电机组中的转子会以静态和动态的力影响径向轴承。由于静力不会在结构中产生振动，因此无法使用加速度计检测轴承的静载荷。然而，例如由大定子偏心引起的静力会在轴承、支撑结构和定子上产生大的力。在装置启动和停机时，力会发生变化，从而导致水电机组中的部件疲劳。结构的动态载荷很容易用振动设备检测，因为变化的力会产生位移，可以使用加速度计或位移传感器检测。
可以使用不同的方法来确定立式水电机组的径向力。对于直接力测定，一个优选但非限制性的实施方式为，将称重传感器安装在枢轴销后面，将应变计安装在轴承枢轴销内，将应变计安装在轴承支架上。
力也可以通过轴承座运动(用加速度计测量)乘以将轴承座连接到混凝土结构的支架的刚度来确定。根据加速度计测量确定的径向力仅包括动态载荷。使用加速度计无法测量静态径向载荷。
步骤1.3，使用轴承特性和轴运动确定水电机组中导向轴承中的径向载荷。
大型水电机组中的径向导向轴承是流体动力径向轴承，通常是可倾瓦类型。轴承中的力是由轴和轴承座之间的相对运动引起的。轴承的动态特性，即K和C，取决于轴承的几何形状、润滑剂的特性、转速、偏心率等。
径向轴承中的轴偏心距ε是径向轴位移与轴承径向游隙之间的关系。在小偏心距(ε<0.6)下，水电机组中使用的轴颈轴承的轴承特性可以被认为是线性的，为了确定可倾瓦轴承的动态特性，需要进行数值计算。一个优选但非限制性的实施方式为，轴承动力学软件通常可以在规定的轴承载荷或轴颈偏心率下计算轴承参数。
图3和图4显示了直接刚度和阻尼特性(kxx和cxx)的计算结果，作为水电机组中可倾瓦轴承在固定轴承游隙和转速下的偏心距的函数。可以理解的是，可以使用任意可用的分析软件执行计算，一个优选但非限制性的实施方式为，使用商业转子动力学分析软件RAPID-RDATM。
如果计算所有轴承特性，即还计算kxy，kyx，kyy，cxy...cyy，mxx...myy，载荷施加在x方向，则根据方程1计算出的轴承载荷作为偏心距的函数将如图5所示。
要根据测量的轴位移和计算的轴承参数确定轴承载荷，必须知道当前轴承间隙的轴承参数以及轴相对于轴承中心的位移。轴承间隙随轴承、周围结构和轴的温度而变化。一个优选但非限制性的实施方式为，在每个轴承上使用四个位移传感器，以90度间隔安装，可以补偿热变化并确定轴承的中心。
轴承中心和间隙是通过使用液压千斤顶在+x、-x、+y和-y方向推动轴并使用传感器记录中心位置和轴承间隙来确定的。当在水电机组上测量轴承游隙并假设由于温度和外力引起的对称形状变化时，则可以描述轴距中心的位移和轴承游隙的变化。cmx表示在x方向上测得的径向轴承游隙。当轴承内的温度发生变化时，由于轴和轴承直径的变化，轴承间隙也会发生变化。当前径向轴承游隙cbx可以通过将轴承Δdb和轴Δds的几何变化与测量的轴承游隙相加或减去来确定，即，cbx＝cmx+Δdb-Δds。
当轴在x方向上位移x距离时轴承表面和轴表面之间的距离是传感器1处的cbx-x和传感器3处的cbx+x。传感器3和传感器1总和的一半代表径向轴承间隙，传感器3和传感器1之间的一半差异提供了轴距轴承中心的位移/>使用相应的推理，也可以得到y轴的位移。
由于热变化引起的轴承游隙的变化通常小于轴承游隙的15％，但这种变化仍会影响轴承参数，需要加以考虑。
通过了解在这些偏心下轴承中的径向轴位移、当前轴承游隙以及轴承特性，可以确定载荷，前提是偏心的大小对轴承特性起决定性作用，而非偏心的动态部分与静态部分之间的关系。
通过计算轴在x和y方向的偏心距，轴的总偏心距εt和相位θ以如下公式5表示，
由此可以绘制轴承载荷与偏心距和轴承间隙的函数关系图，如图6所示，借助该图，可以完成载荷的确定。
步骤1.4，实测轴承载荷和根据轴位移和轴承特性计算的轴承载荷。
实测轴承载荷是使用安装在枢轴销内的应变计测量，轴位移是使用位移传感器测量，测量时还需要同步对轴承游隙进行查表，即依照测量时间获得对应的轴承游隙。为了根据测得的轴位移确定轴承载荷，依照载荷-偏心距-间隙关系乘以轴的偏心距。
步骤2，获取水电机组设计标准中载荷约束。
在制造新水电机组时，部件根据特定的设计标准设计。水力发电厂中的一些机械部件为根据功率、温度、启动和停止、液压载荷、不平衡、自身重量和可能发生的故障等因素而变化的载荷构建。
步骤2.1，根据水电机组组件的载荷谱获得水电机组载荷约束。
水电机组的设计标准规定了机组的使用方式以及在其技术寿命期间可以启动和停机的次数。循环数和事件数为载荷谱提供输入数据，用作单元疲劳设计的支持。包含在谱中的载荷分量示例如表2所示。
使用水电厂确定的载荷谱，设计水电机组时会考虑其结构强度和疲劳。为了避免损坏机械部件，了解部件的设计载荷并监控关键部件中出现的载荷非常重要。对于旋转结构产生的高载荷造成的损坏，损坏多发生在轴承、轴承支架的焊接、联轴器和相互连接处。在设计新水电机组时，制造商会根据水电厂的载荷谱执行有限元分析(FEA)和疲劳分析。为了保护水电机组免受有害载荷的影响，水电厂保留有关这些载荷的信息非常重要。对于无法进行疲劳计算的旧水电机组，应使用FEA确定关键部件，并使用疲劳分析确定关键载荷。水电机组中应该出现的正常载荷水平，包括静态和动态，可以通过平衡标准、发电机形状偏差的限制以及水电机组中测量载荷的经验来估计。
表2水电机组组件的载荷谱示例
步骤2.2，根据平衡等级获得载荷约束。
根据ISO-1940，最大允许径向载荷可以根据所选的平衡质量等级、转子质量和转速进行计算。平衡质量等级决定了允许的残余特定不平衡eper与转速Ω的乘积的最大允许幅度。根据ISO-1940，G6.3是水电机组推荐的平衡质量等级，这意味着eper·Ω＝6.3mm/s。根据方程(6)，平衡质量等级G6.3和G16在不同转子质量和转子速度下的最大允许轴承载荷见表3
式中：
fu表示径向力；
m是转子质量；
Ω是转子速度；
eper是允许的残余特定不平衡。
表3对称定位转子的轴承载荷示例
| Rotormass | Rotorspeed | BearingloadG16 | BearingloadG6.3 |
| (kg) | (RPM) | (kN) | (kN) |
| 200,000 | 167 | 28 | 11 |
| 500,000 | 83 | 35 | 13 |
| 100,000 | 500 | 42 | 16 |
步骤2.3，根据发电机偏心率和转轮不均匀流动特性获得载荷约束。
除了疲劳标准和平衡等级之外，制造商/水电厂还对发电机中允许的最大形状偏差提出了建议。由于不平衡磁拉力(UMP)，相对于转子的偏心定子会在发电机轴承上产生静载荷。转子还将承受频率为1×Ω的循环载荷，这是由偏心定子引起的UMP引起的。相对于轴心偏心的发电机转子中心将导致发电机轴承上的动态载荷和转子上的静态载荷。与蒸汽和燃气发电机组中的发电机相比，水电机组中的UMP力较高；1mm偏心率下的正常UMP力为200–400kN。由水电厂决定发电机允许的最大偏心率；常用值为3％的最大定子偏心和1.5％的转子偏心。
对于转轮中的径向载荷，更难以确定正常运行时的最大允许载荷。由于质量不平衡，转轮的平衡等级设置了最大允许的动态载荷。水轮机中的静载荷是由水轮机中不均匀的流动特性引起的。弗朗西斯水轮机的模型以如下公式表示，
式中：
fr表示转轮的径向载荷；
P表示额定功率；
Ω表示水轮机组同步转速；
D是转轮直径。
使用疲劳计算、平衡等级、发电机形状偏差的建议以及确定的转轮特性和径向载荷之间的关系，可以确定水电机组的径向载荷以及正常运行情况下可以预期的径向载荷。
步骤3，使用步骤1中的模型和步骤2中获得的约束，设置警报和跳闸水平。
监测水电机组振动的方法和允许振动水平的设置应基于机器的机械特性、机器尺寸确定的载荷以及振动引起的载荷。根据测得的载荷水平和机器的设计特性来确定警报和跳闸水平。
步骤3.1，识别水电机组中机械部件的特性。需要对机器的尺寸标注先决条件进行分析，并确定关键部件。必须使用制造商设计数据或执行新的FEA和疲劳计算来确定关键部件的疲劳极限。
步骤3.2，振动与载荷水平之间的关系。(a)在轴承内部或后面安装称重传感器。测量静态和动态载荷。(b)根据计算出的轴承特性和测得的轴位移计算轴承载荷。(c)根据轴承支架的刚度特性和测得的轴承座运动计算轴承载荷(只能使用此方法确定动态载荷)。
步骤3.3，确定正常运行期间的最大载荷水平。如果使用平衡标准(例如ISO1940)，则可以计算满足标准的最大载荷水平。由不平衡引起的载荷是动态载荷。
对于发电机，定子和转子的最大偏心度决定了发电机的最大载荷。最大允许不平衡拉力是通过将气隙(正常气隙：15-25毫米)、偏心率和发电机的磁刚度(通常为200-400106N/m)相乘计算得出的。如果有上、下发电机导向轴承，载荷分布在上、下发电机轴承之间；否则，单个发电机轴承将承担发电机的所有载荷。涡轮导向轴承将承受由转轮引起的所有载荷。
步骤3.4，基于运行模式、平衡等级和疲劳极限设置报警和跳闸阈值。警报和跳闸级别应基于操作模式、平衡等级和疲劳极限。在正常工作条件下，径向载荷不应超过允许动载荷加上最大允许静载荷；如果超过这些载荷水平(通常在75到125kN之间的载荷水平)，机器应发出警报。机器的跳闸等级应为疲劳极限设置较大的裕度。由于转轮的动态行为，在启动期间应允许更高的载荷。如果仅使用加速度计，则无法确定静态载荷(来自定子偏心等)。不建议只使用加速度计。
本发明的有益效果在于，与现有技术相比，根据水电机组的物理性质，即根据水电机组的设计标准和运行过程中测量的径向负荷，提出振动监测方法。有关状态监测的方法基于关键部件的机械特性和作用在这些部件上的测量径向载荷。使用这种方法，可以根据与正常运行期间预期载荷水平相关的径向载荷水平以及关键部件的疲劳极限来确定监控系统的警报和跳闸水平。本发明还提供了一种使用计算出的轴承参数和轴位移测量值来测量水电机组轴承载荷的替代方法。这种方法预计不会像使用称重传感器的载荷测量方法具有相同的高分辨率，但该方法所需的安装快速、简单，并且不需要对水电机组中的组件进行任何修改。
本发明申请人结合说明书附图对本发明的实施示例做了详细的说明与描述，但是本领域技术人员应该理解，以上实施示例仅为本发明的优选实施方案，详尽的说明只是为了帮助读者更好地理解本发明精神，而并非对本发明保护范围的限制，相反，任何基于本发明的发明精神所作的任何改进或修饰都应当落在本发明的保护范围之内。
Claims (3)
1.一种基于机械特性的水电机组保护性振动监测方法，其特征在于，包括以下步骤：
步骤1，构建水电机组径向力和振动模型；包括：
步骤1.1，构建轴承与支架之间的相互作用模型；轴承油膜特性的力平衡模型以如下公式表示，
式中：
f表示由转子引起作用在轴承上的径向载荷，fx表示在x轴的分量，fy表示在y轴的分量；k表示轴承的刚度，kxx表示其x轴的分量，kyy表示其y轴的分量，kxy，kyx表示其交叉耦合刚度，c表示轴承的阻尼，cxx表示其x轴的分量，cyy表示其y轴的分量，cxy，cyx表示其交叉耦合阻尼，uSx表示轴的绝对位移矢量uS的x轴分量，uSy表示y轴分量；uHx表示轴承箱的绝对位移矢量uH的x轴分量，uHy表示y轴分量；
位移与组合刚度和阻尼特性之间的关系以如下公式表示，
D＝[I+ZH-1]-1 (4)
式中：
D表示轴承和支架的组合特性，D＝Kc+iΩCc；Kc是组合刚度，Cc是组合阻尼；
步骤1.2，测量或计算获得水电机组径向力；
步骤1.3，使用轴承特性和轴运动确定水电机组中导向轴承中的径向载荷；通过计算轴在x和y方向的偏心距，轴的总偏心距εt和相位θ以如下公式表示，
式中：
cbx表示径向轴承游隙x轴的分量；cby表示径向轴承游隙y轴的分量；
步骤1.4，实测轴承载荷和根据轴位移和轴承特性计算的轴承载荷；
步骤2，获取水电机组设计标准中载荷约束；包括：
步骤2.1，根据水电机组组件的载荷谱获得水电机组载荷约束；
步骤2.2，根据平衡等级获得载荷约束；
步骤2.3，根据发电机偏心率和转轮不均匀流动特性获得载荷约束；弗朗西斯水轮机的模型以如下公式表示，
式中：
fr表示转轮的径向载荷P表示额定功率；Ω表示水轮机组同步转速D是转轮直径；
步骤3，使用步骤1中的模型和步骤2中获得的约束，设置警报和跳闸水平。
2.根据权利要求1所述的一种基于机械特性的水电机组保护性振动监测方法，其特征在于：
步骤2.1中，循环数和事件数为载荷谱提供输入数据，用作单元疲劳设计的支持。
3.根据权利要求1或2所述的一种基于机械特性的水电机组保护性振动监测方法，其特征在于：
步骤3进一步包括：
步骤3.1，识别水电机组中机械部件的特性；
步骤3.2，振动与载荷水平之间的关系；包括：
(a)在轴承内部或后面安装称重传感器，测量静态和动态载荷；
(b)根据计算出的轴承特性和测得的轴位移计算轴承载荷；
(c)根据轴承支架的刚度特性和测得的轴承座运动计算轴承载荷；
步骤3.3，确定正常运行期间的最大载荷水平；
步骤3.4，基于运行模式、平衡等级和疲劳极限设置报警和跳闸阈值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111089115.9A CN113847981B (zh) | 2021-09-16 | 2021-09-16 | 一种基于机械特性的水电机组保护性振动监测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111089115.9A CN113847981B (zh) | 2021-09-16 | 2021-09-16 | 一种基于机械特性的水电机组保护性振动监测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113847981A CN113847981A (zh) | 2021-12-28 |
| CN113847981B true CN113847981B (zh) | 2024-05-24 |
Family
ID=78974391
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111089115.9A Active CN113847981B (zh) | 2021-09-16 | 2021-09-16 | 一种基于机械特性的水电机组保护性振动监测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113847981B (zh) |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202869806U (zh) * | 2012-08-24 | 2013-04-10 | 国电联合动力技术有限公司 | 一种测量风力发电机组振动与载荷对应关系的装置 |
| CN105865783A (zh) * | 2016-03-23 | 2016-08-17 | 湖南大学 | 一种测量滑动轴承油膜特性的方法 |
| CN106194598A (zh) * | 2016-07-15 | 2016-12-07 | 山东中车风电有限公司 | 一种风电机组状态的监控方法 |
| CN107843429A (zh) * | 2016-09-19 | 2018-03-27 | 舍弗勒技术股份两合公司 | 轴承状态监测控制方法及控制装置、监测设备、监测方法 |
| CN107893682A (zh) * | 2016-10-03 | 2018-04-10 | 通用电气公司 | 用于检测润滑轴承状态的系统的方法 |
| CN108225771A (zh) * | 2017-12-26 | 2018-06-29 | 新昌县江北轴承有限公司 | 一种滚子轴承振动测量装置 |
| CN108426714A (zh) * | 2018-03-15 | 2018-08-21 | 国家电网公司 | 一种可倾式滑动轴承工作状态的监测系统的监测方法 |
| CN109635438A (zh) * | 2018-12-12 | 2019-04-16 | 哈动国家水力发电设备工程技术研究中心有限公司 | 一种水轮机顶盖振动数值的确定方法 |
| CN111504587A (zh) * | 2020-05-30 | 2020-08-07 | 西安热工研究院有限公司 | 一种立式水轮发电机导轴承支架振动特性评估系统及方法 |
| CN111597651A (zh) * | 2020-04-30 | 2020-08-28 | 上海工程技术大学 | 一种基于hwpso-svdd模型的滚动轴承性能退化评估方法 |
| CN111608842A (zh) * | 2020-05-06 | 2020-09-01 | 哈尔滨电机厂有限责任公司 | 一种水轮机顶盖径向刚度检验的方法 |
| CN112284575A (zh) * | 2020-09-18 | 2021-01-29 | 西安交通大学 | 一种基于电涡流位移传感器的轴承载荷在线监测方法 |
| CN112464529A (zh) * | 2020-11-16 | 2021-03-09 | 西安热工研究院有限公司 | 柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法 |
| CN112834222A (zh) * | 2021-02-02 | 2021-05-25 | 厦门物之联智能科技有限公司 | 动态监测列车轴承使用寿命的方法及电子设备 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| IT201900024030A1 (it) * | 2019-12-16 | 2021-06-16 | Skf Ab | A system and method for determining bearing preload by vibration measurement |
-
2021
- 2021-09-16 CN CN202111089115.9A patent/CN113847981B/zh active Active
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202869806U (zh) * | 2012-08-24 | 2013-04-10 | 国电联合动力技术有限公司 | 一种测量风力发电机组振动与载荷对应关系的装置 |
| CN105865783A (zh) * | 2016-03-23 | 2016-08-17 | 湖南大学 | 一种测量滑动轴承油膜特性的方法 |
| CN106194598A (zh) * | 2016-07-15 | 2016-12-07 | 山东中车风电有限公司 | 一种风电机组状态的监控方法 |
| CN107843429A (zh) * | 2016-09-19 | 2018-03-27 | 舍弗勒技术股份两合公司 | 轴承状态监测控制方法及控制装置、监测设备、监测方法 |
| CN107893682A (zh) * | 2016-10-03 | 2018-04-10 | 通用电气公司 | 用于检测润滑轴承状态的系统的方法 |
| CN108225771A (zh) * | 2017-12-26 | 2018-06-29 | 新昌县江北轴承有限公司 | 一种滚子轴承振动测量装置 |
| CN108426714A (zh) * | 2018-03-15 | 2018-08-21 | 国家电网公司 | 一种可倾式滑动轴承工作状态的监测系统的监测方法 |
| CN109635438A (zh) * | 2018-12-12 | 2019-04-16 | 哈动国家水力发电设备工程技术研究中心有限公司 | 一种水轮机顶盖振动数值的确定方法 |
| CN111597651A (zh) * | 2020-04-30 | 2020-08-28 | 上海工程技术大学 | 一种基于hwpso-svdd模型的滚动轴承性能退化评估方法 |
| CN111608842A (zh) * | 2020-05-06 | 2020-09-01 | 哈尔滨电机厂有限责任公司 | 一种水轮机顶盖径向刚度检验的方法 |
| CN111504587A (zh) * | 2020-05-30 | 2020-08-07 | 西安热工研究院有限公司 | 一种立式水轮发电机导轴承支架振动特性评估系统及方法 |
| CN112284575A (zh) * | 2020-09-18 | 2021-01-29 | 西安交通大学 | 一种基于电涡流位移传感器的轴承载荷在线监测方法 |
| CN112464529A (zh) * | 2020-11-16 | 2021-03-09 | 西安热工研究院有限公司 | 柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法 |
| CN112834222A (zh) * | 2021-02-02 | 2021-05-25 | 厦门物之联智能科技有限公司 | 动态监测列车轴承使用寿命的方法及电子设备 |
Non-Patent Citations (2)
| Title |
|---|
| 振动测试齿轮箱轴承处动态载荷逆向识别研究;姜雨霆等;设备管理与维修(第03期);第31-35页 * |
| 深沟球轴承三维非线性时变振动特性研究;邵毅敏等;振动工程学报(第06期);第29-36页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113847981A (zh) | 2021-12-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| De Castro et al. | Whirl and whip instabilities in rotor-bearing system considering a nonlinear force model | |
| Muszynska | Vibrational diagnostics of rotating machinery malfunctions | |
| Muszynska | Whirl and whip—rotor/bearing stability problems | |
| Zhang et al. | Failure investigation of a Kaplan turbine blade | |
| Hou et al. | Measurement of load distribution in a cylindrical roller bearing with an instrumented housing: Finite element validation and experimental study | |
| Kumar et al. | Dynamic analysis and identification of unbalance and misalignment in a rigid rotor with two offset discs levitated by active magnetic bearings: a novel trial misalignment approach | |
| Egusquiza et al. | Failure investigation of a Pelton turbine runner | |
| Bin et al. | Development of whole-machine high speed balance approach for turbomachinery shaft system with N+ 1 supports | |
| JP2012255688A (ja) | 軸受試験装置 | |
| Nässelqvist et al. | A methodology for protective vibration monitoring of hydropower units based on the mechanical properties | |
| Tiwari et al. | An innovative virtual trial misalignment approach for identification of unbalance, sensor and active magnetic bearing misalignment along with its stiffness parameters in a magnetically levitated flexible rotor system | |
| Safizadeh et al. | Prediction of oil whirl initiation in journal bearings using multi-sensors data fusion | |
| Nässelqvist et al. | Bearing load measurement in a hydropower unit using strain gauges installed inside pivot pin | |
| Shamsah et al. | Estimating rotor unbalance from a single run-up and using reduced sensors | |
| Blaut et al. | Application of the Teager-Kaiser energy operator in diagnostics of a hydrodynamic bearing | |
| CN113847981B (zh) | 一种基于机械特性的水电机组保护性振动监测方法 | |
| Shen et al. | Experimental and numerical analysis of nonlinear dynamics of rotor–bearing–seal system | |
| CN106989879A (zh) | 一种汽轮发电机机座模态测试方法 | |
| KR101977761B1 (ko) | 회전기기의 불평형 감시방법 및 그 장치 | |
| CN216381720U (zh) | 变形评估设备和风力发电机组 | |
| Sinha et al. | Estimating the static load on the fluid bearings of a flexible machine from run-down data | |
| Brito Junior et al. | Using simplified models to assist fault detection and diagnosis in large hydrogenerators | |
| Brito Jr et al. | Experimental estimation of journal bearing stiffness for damage detection in large hydrogenerators | |
| Jamil et al. | Study the Dynamic Behavior of Rotor Supported on a Worn Journal Bearings | |
| Shen et al. | Effect of the seal force on nonlinear dynamics and stability of the rotor-bearing-seal system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant |