CN109029689B - 一种基于转子两端运动轨迹的旋转机械振动分析方法 - Google Patents

一种基于转子两端运动轨迹的旋转机械振动分析方法 Download PDF

Info

Publication number
CN109029689B
CN109029689B CN201810893544.3A CN201810893544A CN109029689B CN 109029689 B CN109029689 B CN 109029689B CN 201810893544 A CN201810893544 A CN 201810893544A CN 109029689 B CN109029689 B CN 109029689B
Authority
CN
China
Prior art keywords
vibration
rotor
radius
amplitude
rotating
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
Application number
CN201810893544.3A
Other languages
English (en)
Other versions
CN109029689A (zh
Inventor
张文涛
房林铁
夏亚磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Datang Corp Science and Technology Research Institute Co Ltd Huazhong Branch
Original Assignee
China Datang Corp Science and Technology Research Institute Co Ltd Huazhong Branch
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Datang Corp Science and Technology Research Institute Co Ltd Huazhong Branch filed Critical China Datang Corp Science and Technology Research Institute Co Ltd Huazhong Branch
Priority to CN201810893544.3A priority Critical patent/CN109029689B/zh
Publication of CN109029689A publication Critical patent/CN109029689A/zh
Application granted granted Critical
Publication of CN109029689B publication Critical patent/CN109029689B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Testing Of Balance (AREA)

Abstract

本发明涉及一种基于转子两端运动轨迹的旋转机械振动分析方法，技术方案是，根据得到的同向和反向振动分量值，判断转子上不平衡力型式：若At大于旋转机械振动标准中规定的优秀值，转子上同向不平衡力大；若Af大于旋转机械振动标准中规定的优秀值，转子上反向不平衡力大；若At和Af都大于旋转机械振动标准中规定的优秀值，转子上混合型式的不平衡力大，本发明方法简单，振动分析和不平衡力型式分析时综合考虑了垂直和水平方向振动，而不是仅仅依据单一方向振动，所得到的分析结论全面、可靠，是旋转机械振动分析方法上的创新。

Description

一种基于转子两端运动轨迹的旋转机械振动分析方法
技术领域
本发明涉及一种旋转机械振动测试分析方法，可以帮助技术人员根据振动测试数据分析转子上不平衡力的分布型式，如力不平衡、力偶不平衡和混合不平衡，指导开展旋转机械动平衡工作，主要应用领域包括：各类旋转机械，如泵、风机、电动机、汽轮发电机组等。
背景技术
不平衡是旋转机械最常见的振动故障，转子动平衡是旋转机械振动故障治理的有效手段，也是经常要开展的工作。准确判定不平衡型式可以提高动平衡效率和精度、减少开机次数，对于提高旋转机械运行稳定性和可靠性具有重要意义。
不平衡包括静不平衡、动不平衡和混合式不平衡。静不平衡状态下，转子中部存在不平衡力，或者在转子两端存在一组大小相等，方向相同的力，又称为力不平衡。动不平衡状态下，在转子两端存在一组大小相等，方向相反的力，又称为力偶不平衡。当转子中同时存在力不平衡和力偶不平衡时，称为混合式不平衡。开展转子动平衡之前，必须首先了解转子上不平衡力的型式和分布。
机组实际运行状态下，无法直接测试到不平衡力，都是根据振动来诊断不平衡故障，并分析转子上是否有不平衡力。振动大，不平衡力大；振动小，不平衡力小。同样，也是根据转子两端的振动来判断不平衡力的型式。如转子两端振动信号中存在较大的同向分量，则认为转子上存在力不平衡；如转子两端振动信号中存在较大的反向分量，则认为转子上存在力偶不平衡；如转子两端振动信号中同时存在较大的同向分量和反向，则认为转子上存在混合型式的不平衡。
目前判断不平衡力型式都是依靠垂直或水平单一方向上的振动进行的。对于汽轮发电机组等大型旋转机械而言，由于机组垂直和水平方向上的支撑刚度、阻尼特性等存在较大差异，转子上不平衡力所引起的轴承振动在垂直和水平方向上也有较大差别。实际工程中，有时会出现根据垂直方向振动和根据水平方向振动所判断得到的不平衡力型式完全不同的情况。因此，其改进和创新势在必行。
发明内容
针对上述情况，为克服现有技术之缺陷，本发明之目的就是提供一种基于转子两端运动轨迹的旋转机械振动分析方法，可有效解决分析转子上的不平衡型式，提高动平衡的效率的问题。
本发明解决的技术方案是：
一种基于转子两端运动轨迹的旋转机械振动分析方法，包括以下步骤：
(1)在转轴上设置一个反光标记，将光电传感器对准反光标记，或者在转轴上开一键槽，将涡流传感器对准键槽，用光电传感器配合反光标记或涡流传感器配合键槽的方法，测取键相脉冲信号，转轴每旋转一周，产生一个这样的脉冲信号，记相邻两个脉冲之间的时间间隔为T，则实时转动频率f＝1/T；
(2)在转子两端待测部位的垂直和水平方向各布置一组振动传感器，以键相脉冲信号为触发基准，同步采集转子两端垂直和水平方向振动信号，记为x1(t),x2(t),y1(t),y2(t)；
(3)对采集到的振动信号进行快速傅里叶变换，得到与转动频率相同的工频谐波分量的幅值和相位，分别记为：Ax1,Ax2,Ay1,Ay2和
Figure GDA0002433556880000021
(4)根据一倍频谐波分量的幅值和相位，重构出与转子两端垂直和水平振动信号相对应的一倍频谐波分量信号：
Figure GDA0002433556880000022
Figure GDA0002433556880000023
Figure GDA0002433556880000024
Figure GDA0002433556880000025
(5)记与转动频率f相对应的周期为T，则有：T＝1/f；将[0,T]间划分为1000段，得到转动一周内不同时刻下各点振动值：x′1(ti),y′1(ti),x′2(ti),y′2(ti)；
(7)定义不同时刻下转子振动半径ri
Figure GDA0002433556880000026
Figure GDA0002433556880000027
(8)绘制出转子两端在不同时刻下由振动半径组成的振动曲线；
(9)找到t＝0时刻后转子两端振动半径曲线上第1个振动半径峰值点所对应的时间点及其幅值，分别记为：t′1,A1和t′2,A2；
(10)计算出与振动半径峰值点相对应的相位角ψ1,ψ2；
Figure GDA0002433556880000031
Figure GDA0002433556880000032
(11)以振动半径A1∠ψ1,A2∠ψ2作为转子两端振动，根据转子两端振动，将其分解为同向振动At以及反向振动Af；
Figure GDA0002433556880000033
Figure GDA0002433556880000034
(12)根据得到的同向和反向振动分量值，判断转子上不平衡力型式：
若At大于旋转机械振动标准中规定的优秀值，转子上同向不平衡力大；
若Af大于旋转机械振动标准中规定的优秀值，转子上反向不平衡力大；
若At和Af都大于旋转机械振动标准中规定的优秀值，转子上混合型式的不平衡力大。
本发明方法简单，认为转子两端振动是在垂直于轴线的平面内进行的，需要用垂直和水平两个方向上的振动来综合表示，据此，提出根据转子两端垂直和水平方向上的振动信号，求得转子两端在垂直轴线平面内的运动轨迹，定义不同时刻下的振动半径为轨迹上该时刻点到圆心的距离，求得转子两端振动半径变化曲线后，以t＝0时刻后第1个振动半径峰值所在点的相位角为该端相位角，以振动半径峰值作为该端振动。在此基础上，对转子两端振动信息进行分析，判断转子上不平衡型式。将转子振动看作是在垂直于轴线的平面内的平面轨迹运动，振动分析和不平衡力型式分析时综合考虑了垂直和水平方向振动，而不是仅仅依据单一方向振动，所得到的分析结论全面、可靠，是旋转机械振动分析方法上的创新。
附图说明
图1是转子***振动测试示意简图。
图中各个附图标记分别表示：1轴承、2转子、3振动传感器、4键相传感器、5振动数据采集仪、6轮盘
图2是转子测试得到的两端4个振动测点原始波形。
其中：x1和y1是一轴承侧水平和垂直方向振动，x2和y2是另一轴承侧水平和垂直方向振动；
图3是测试得到的键相脉冲信号图。
图4是测试得到的转子两端4个测点原始振动频谱图，图上标注了一倍频谐波分量的幅值和相位。
图5是利用本方法将测试信号重构后得到的与转子两端4个测点一倍频谐波分量相对应的波形图。
图6是利用本方法得到的转子两端振动半径波形图。
图7本发明方法的流程图
具体实施方式
以下结合附图对本发明的具体实施方式作进一步详细说明。
由图1-7给出，本发明一种基于转子两端运动轨迹的旋转机械振动分析方法，包括以下步骤：
(1)在转轴上设置一个反光标记，将光电传感器对准反光标记，或者在转轴上开一键槽，将涡流传感器对准键槽，用光电传感器配合反光标记或涡流传感器配合键槽的方法，测取键相脉冲信号，转轴每旋转一周，产生一个这样的脉冲信号，记相邻两个脉冲之间的时间间隔为T，则实时转动频率f＝1/T；
(2)在转子两端待测部位的垂直和水平方向各布置一组振动传感器，以键相脉冲信号为触发基准，同步采集转子两端垂直和水平方向振动信号，记为x1(t),x2(t),y1(t),y2(t)；
(3)对采集到的振动信号进行快速傅里叶变换，得到与转动频率相同的工频谐波分量的幅值和相位，分别记为：Ax1,Ax2,Ay1,Ay2和
Figure GDA0002433556880000041
(4)根据一倍频谐波分量的幅值和相位，重构出与转子两端垂直和水平振动信号相对应的一倍频谐波分量信号：
Figure GDA0002433556880000042
Figure GDA0002433556880000043
Figure GDA0002433556880000044
Figure GDA0002433556880000045
(5)记与转动频率f相对应的周期为T，则有：T＝1/f；将[0,T]间划分为1000段，得到转动一周内不同时刻下各点振动值：x′1(ti),y′1(ti),x′2(ti),y′2(ti)；
(7)定义不同时刻下转子振动半径ri
Figure GDA0002433556880000051
Figure GDA0002433556880000052
(8)绘制出转子两端在不同时刻下由振动半径组成的振动曲线；
(9)找到t＝0时刻后转子两端振动半径曲线上第1个振动半径峰值点所对应的时间点及其幅值，分别记为：t′1,A1和t′2,A2；
(10)计算出与振动半径峰值点相对应的相位角ψ1,ψ2；
Figure GDA0002433556880000053
Figure GDA0002433556880000054
(11)以振动半径A1∠ψ1,A2∠ψ2作为转子两端振动，根据转子两端振动，将其分解为同向振动At以及反向振动Af；
Figure GDA0002433556880000055
Figure GDA0002433556880000056
(12)根据得到的同向和反向振动分量值，判断转子上不平衡力型式：
若At大于旋转机械振动标准中规定的优秀值，转子上同向不平衡力大；
若Af大于旋转机械振动标准中规定的优秀值，转子上反向不平衡力大；
若At和Af都大于旋转机械振动标准中规定的优秀值，转子上混合型式的不平衡力大。
所述键相传感器、振动传感器和数据采集仪均为现有技术，如ROS-W型光电键相传感器、ZXP型振动传感器、EVM-8型振动数据采集仪。
由上述情况可以看出，转子振动和轴承振动作用在垂直于轴线的平面上，这是一个平面运动问题，仅仅依靠单一方向上的振动来进行分析是不完备的，一种发展趋势就是将垂直和水平方向上的振动综合起来分析。例如，传统的频谱分析技术是针对单一信号进行的，而全谱分析技术则是在传统频谱分析技术基础上发展起来的，综合利用了垂直和水平两个方向上的振动信号，可以在传统频谱基础上进一步得到正进动和反进动等深层次信息和特征。本发明认为转子两端振动是在垂直于轴线的平面内进行的，需要用垂直和水平两个方向上的振动来综合表示。据此，提出根据转子两端垂直和水平方向上的振动信号，求得转子两端在垂直轴线平面内的运动轨迹，定义不同时刻下的振动半径为轨迹上该时刻点到圆心的距离，求得转子两端振动半径变化曲线后，以t＝0时刻后第1个振动半径峰值所在点的相位角为该端相位角，以振动半径峰值作为该端振动。在此基础上，对转子两端振动信息进行分析，判断转子上不平衡型式。其方法简单，将转子振动看作是在垂直于轴线的平面内的平面轨迹运动，振动分析和不平衡力型式分析时综合考虑了垂直和水平方向振动，而不是仅仅依据单一方向振动，所得到的分析结论全面、可靠，是旋转机械振动分析方法上的创新经实际应用，取得了良好的技术效果：
如下表所示，是利用本方法得到的转子两端振动半径幅值和相位，以x向振动为依据，同向分量幅值小，反向分量幅值大，转子存在反向型式不平衡力。以Y向振动为依据，同向分量幅值大，反向分量幅值小，转子存在同向型式不平衡力。根据x和y向振动对不平衡型式的判断结论不同，存在矛盾。以本方法得到的振动半径为依据，同向分量幅值大于旋转机械振动标准中规定的优秀值，反向分量幅值小于旋转机械振动标准中规定的优秀值，转子存在同向型式不平衡力。本方法综合考虑了x和y两个方向上的振动，而不是仅仅依据单一方向振动，所得到的分析结论全面、可靠，是旋转机械振动分析方法上的创新经实际应用。
Figure GDA0002433556880000061

Claims (1)

1.一种基于转子两端运动轨迹的旋转机械振动分析方法，其特征在于，包括以下步骤：
(1)在转轴上设置一个反光标记，将光电传感器对准反光标记，或者在转轴上开一键槽，将涡流传感器对准键槽，用光电传感器配合反光标记或涡流传感器配合键槽的方法，测取键相脉冲信号，转轴每旋转一周，产生一个这样的脉冲信号，记相邻两个脉冲之间的时间间隔为T，则实时转动频率f＝1/T；
(2)在转子两端待测部位的垂直和水平方向各布置一组振动传感器，以键相脉冲信号为触发基准，同步采集转子两端垂直和水平方向振动信号，记为x1(t),x2(t),y1(t),y2(t)；
(3)对采集到的振动信号进行快速傅里叶变换，得到与转动频率相同的工频谐波分量的幅值和相位，分别记为：Ax1,Ax2,Ay1,Ay2和
Figure FDA0002433556870000011
(4)根据一倍频谐波分量的幅值和相位，重构出与转子两端垂直和水平振动信号相对应的一倍频谐波分量信号：
Figure FDA0002433556870000012
Figure FDA0002433556870000013
Figure FDA0002433556870000014
Figure FDA0002433556870000015
(5)记与转动频率f相对应的周期为T，则有：T＝1/f；将[0,T]间划分为1000段，得到转动一周内不同时刻下各点振动值：x'1(ti),y'1(ti),x'2(ti),y'2(ti)；
(7)定义不同时刻下转子振动半径ri
Figure FDA0002433556870000016
Figure FDA0002433556870000017
(8)绘制出转子两端在不同时刻下由振动半径组成的振动曲线；
(9)找到t＝0时刻后转子两端振动半径曲线上第1个振动半径峰值点所对应的时间点及其幅值，分别记为：t′1,A1和t′2,A2；
(10)计算出与振动半径峰值点相对应的相位角ψ1,ψ2；
Figure FDA0002433556870000021
Figure FDA0002433556870000022
(11)以振动半径A1∠ψ1,A2∠ψ2作为转子两端振动，根据转子两端振动，将其分解为同向振动At以及反向振动Af；
Figure FDA0002433556870000023
Figure FDA0002433556870000024
(12)根据得到的同向和反向振动分量值，判断转子上不平衡力型式：
若At大于旋转机械振动标准中规定的优秀值，转子上同向不平衡力大；
若Af大于旋转机械振动标准中规定的优秀值，转子上反向不平衡力大；
如果At和Af都大于旋转机械振动标准中规定的优秀值，转子上混合型式的不平衡力大。
CN201810893544.3A 2018-08-08 2018-08-08 一种基于转子两端运动轨迹的旋转机械振动分析方法 Active CN109029689B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810893544.3A CN109029689B (zh) 2018-08-08 2018-08-08 一种基于转子两端运动轨迹的旋转机械振动分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810893544.3A CN109029689B (zh) 2018-08-08 2018-08-08 一种基于转子两端运动轨迹的旋转机械振动分析方法

Publications (2)

Publication Number Publication Date
CN109029689A CN109029689A (zh) 2018-12-18
CN109029689B true CN109029689B (zh) 2020-05-22

Family

ID=64650071

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810893544.3A Active CN109029689B (zh) 2018-08-08 2018-08-08 一种基于转子两端运动轨迹的旋转机械振动分析方法

Country Status (1)

Country Link
CN (1) CN109029689B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110849414B (zh) * 2019-10-29 2021-05-04 润电能源科学技术有限公司 转子弯曲方向的识别方法、装置、设备及存储介质
CN110779438A (zh) * 2019-10-29 2020-02-11 润电能源科学技术有限公司 转子弯曲值测量方法、汽轮发电机及计算机可读存储介质
CN111473933B (zh) * 2020-04-27 2021-11-02 上海海事大学 多功能叶片与转子试验台
CN111550429B (zh) * 2020-05-14 2021-07-20 北京化工大学 一种离心压缩机转子***稳定性综合控制方法及装置
CN113191057B (zh) * 2021-05-12 2022-04-12 四川长虹空调有限公司 不平衡力方向的确定方法及装置
CN113358203B (zh) * 2021-06-02 2024-06-11 中国大唐集团科学技术研究院有限公司华东电力试验研究院 一种基于谐分量分解的转子固有频率识别方法和***

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1269504A (zh) * 1999-04-06 2000-10-11 黄震西 转子不平衡量两点测定法
CN101487756A (zh) * 2009-01-13 2009-07-22 东南大学 旋转机械振动分析中的谐分量转速平衡方法
CN107255549A (zh) * 2017-06-14 2017-10-17 西安交通大学 一种转子多点不平衡量识别方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2930489B1 (en) * 2014-04-09 2019-08-28 Balance Systems S.r.L. Balancing Process and Device for a Rotating Body

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1269504A (zh) * 1999-04-06 2000-10-11 黄震西 转子不平衡量两点测定法
CN101487756A (zh) * 2009-01-13 2009-07-22 东南大学 旋转机械振动分析中的谐分量转速平衡方法
CN107255549A (zh) * 2017-06-14 2017-10-17 西安交通大学 一种转子多点不平衡量识别方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
柔性转子转轴弯曲与不平衡耦合振动分析;夏亚磊;《动力工程学报》;20161130;第36卷(第11期);第877-882页 *

Also Published As

Publication number Publication date
CN109029689A (zh) 2018-12-18

Similar Documents

Publication Publication Date Title
CN109029689B (zh) 一种基于转子两端运动轨迹的旋转机械振动分析方法
Akin et al. A simple real-time fault signature monitoring tool for motor-drive-embedded fault diagnosis systems
CN104400560B (zh) 一种数控机床切削工况下主轴轴心轨迹在线测量方法
CN101556200B (zh) 基于全矢谱的柔性转子动平衡方法
Zhou et al. Study on intra-wave frequency modulation phenomenon in detection of rub-impact fault
CN111649921A (zh) 一种旋转机械故障诊断测试***及工作方法
CN110579312B (zh) 一种无试重旋转机械多轮盘轴系动平衡故障检测方法
Lu et al. A new method for characterizing axis of rotation radial error motion: Part 1. Two-dimensional radial error motion theory
Vafaei et al. Indicated repeatable runout with wavelet decomposition (IRR-WD) for effective determination of bearing-induced vibration
Thiele et al. Computationally efficient method for identifying manufacturing induced rotor and stator misalignment in permanent magnet brushless machines
Ahmed et al. Detection of eccentricity faults in machine usingfrequency spectrum technique
Chirindo et al. Analysis of non-intrusive rotor speed estimation techniques for inverter-fed induction motors
Rosero et al. Fault detection of eccentricity and bearing damage in a PMSM by means of wavelet transforms decomposition of the stator current
KR102224333B1 (ko) Cnc공작기계 스핀들주축의 런아웃 검사장치 및 검사방법
CN110646138A (zh) 一种旋转机械无键相无试重动平衡方法和分析装置
CN106441840A (zh) 一种基于模型的旋转机械伪次同步故障全息诊断方法
CN111307426A (zh) 一种基于FrFT-EWT原理的旋转机械故障特征提取方法
CN103712746B (zh) 一种转子动平衡试验中机械滞后角的求取方法
CN104165729B (zh) 一种高速转子的动平衡方法
CN115356108B (zh) 一种调制高阶水平挤压变换机械故障诊断方法与装置
Chai et al. Misalignment detection of rotor system based on adaptive input-output model identification of motor speed
CN112945535B (zh) 一种基于数值模拟的旋转机械故障检测方法及装置
Zhang et al. Order based modal analysis using Vold-Kalman filter
CN1103443C (zh) 非对称转子的全息动平衡方法
Dreher et al. Directional coordinates for the identification of backward and forward frequencies of rotating machines via OMA

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
GR01 Patent grant