CN103412145A - 一种转子系统轴心轨迹进动方向自动识别方法 - Google Patents

一种转子系统轴心轨迹进动方向自动识别方法 Download PDF

Info

Publication number
CN103412145A
CN103412145A CN2013103615253A CN201310361525A CN103412145A CN 103412145 A CN103412145 A CN 103412145A CN 2013103615253 A CN2013103615253 A CN 2013103615253A CN 201310361525 A CN201310361525 A CN 201310361525A CN 103412145 A CN103412145 A CN 103412145A
Authority
CN
China
Prior art keywords
axis
rotor
arc
trajectory
alpha
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.)
Granted
Application number
CN2013103615253A
Other languages
English (en)
Other versions
CN103412145B (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.)
North China Electric Power University
Original Assignee
North China Electric Power University
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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN201310361525.3A priority Critical patent/CN103412145B/zh
Publication of CN103412145A publication Critical patent/CN103412145A/zh
Application granted granted Critical
Publication of CN103412145B publication Critical patent/CN103412145B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Operation Control Of Excavators (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明公开了轨迹检测技术领域的一种转子系统轴心轨迹进动方向自动识别方法。首先,设定采样频率,利用振动传感器采集转子的径向振动信号,并通过快速傅立叶变换FFT得到幅值谱;其次,采用相位差校正方法对转子振动信号中的主要频率成分进行校正后重构信号;最后,利用重构信号合成转子提纯轴心轨迹,并由提纯轴心轨迹通过局部弧段的进动方向识别得到整个轴心轨迹的进动方向。本发明方法从局部角度出发,对轴心轨迹各个局部弧段的进动方向进行识别,不受轴心轨迹复杂程度的影响,通用性更强。

Description

一种转子系统轴心轨迹进动方向自动识别方法
技术领域
本发明属于轨迹检测技术领域,尤其涉及一种转子系统轴心轨迹进动方向自动识别方法。
背景技术
轴心轨迹是判断旋转机械运行状态的一个重要依据。轴心轨迹的形状以及轴心轨迹的进动方向是轴心轨迹的两个重要特征。当转子处于不同运行状态时,轴心轨迹的形状将发生变化,进动方向也可能改变,如不平衡故障转子的轴心轨迹呈椭圆形,正进动;不对中故障轴心轨迹呈香蕉形或8字形,正进动与反进动并存。因此,通过识别轴心轨迹的形状以及进动方向可判断转子的运行状态。在轴心轨迹进动方向识别方面,不少学者已进行了研究。
有的学者将椭圆轴心轨迹视为两个正圆轨迹的合成,其中一个正进动,一个反进动。当正进动圆比反进动圆大时,为正进动;反之,为反进动。该方法仅能识别椭圆形轴心轨迹进动方向,不能识别复杂的轴心轨迹进动方向。有的学者利用轴心轨迹上连续的几个点进行判断,即选取连续的几个点(xi,yi)(i=1,2,3…N),计算点(xi,yi)到原点的连线与x轴的夹角αi=arctan(yi/xi),若αii+1i+2,则轴心轨迹的进动方向为正进动;反之,为反进动。该方法可识别包含坐标原点的简单轴心轨迹进动方向,但是不能识别复杂的轴心轨迹进动方向。还有一些学者通过截取轴心轨迹中某一旋转周期内的采样点构成复杂多边形,利用平面图形的平移变换判断多边形各顶点的凹凸性和多边形在各顶点处的旋转方向,进而得到轴心轨迹的进动方向。该方法计算较为复杂,并且对于复杂的轴心轨迹,只能识别轴心轨迹部分区域的进动方向,其他区域的进动方向需要人为给出该区域的凹凸特性进行判断,效率较低。
发明内容
针对背景技术中提到的现有的轴心轨迹识别方法在识别复杂的轴心轨迹时,只能识别轴心轨迹部分区域的进动方向,识别算法复杂且效率低下的问题,本发明提出了一种转子系统轴心轨迹进动方向自动识别方法。
一种转子系统轴心轨迹进动方向自动识别方法,其特征在于,所述方法具体包括以下步骤:
步骤1:设定采样频率,利用振动传感器采集转子X轴和Y轴方向的径向振动信号,将采集的电压信号保存到计算机中;
步骤2:对转子X轴和Y轴方向的振动信号通过快速傅立叶变换FFT得到幅值谱;
步骤3:采用相位差频谱校正法提取X轴和Y轴方向振动信号的主要频率分量;
步骤4:将步骤3得到的X轴和Y轴方向各主要频率分量的幅值、频率和相位合成波形,并叠加得到X轴和Y轴方向重构信号,再利用重构信号合成转子提纯轴心轨迹;
步骤5:利用步骤4得到的转子提纯轴心轨迹确定转子轴心轨迹的进动方向。
步骤3中,X轴和Y轴方向振动信号的主要频率分量包括信号的基频、分频及倍频。
步骤5中,利用步骤4得到的转子提纯轴心轨迹确定转子轴心的进动方向的具体过程为:
步骤501:设P1、P2、P3为转子提纯轴心轨迹上连续三点,其坐标分别为(x1,y1)、(x2,y2)、(x3,y3);通过这三个连续的点确定一段圆弧,圆弧圆心为O1,圆弧圆心坐标为(a,b);
步骤502:将原坐标系XOY进行平移,平移后坐标系为XO1Y;分别计算
Figure BDA0000368624140000031
Figure BDA0000368624140000032
Figure BDA0000368624140000034
轴的夹角角度α1、α2、α3
步骤503:将XO1Y坐标系绕点O1逆时针旋转α1,使
Figure BDA0000368624140000035
轴和
Figure BDA0000368624140000036
重合,即令α1=0,进而得新坐标系X1O1Y1;此时,
Figure BDA0000368624140000038
Figure BDA0000368624140000039
轴的夹角分别为 α 1 ′ = 0 α 2 ′ = α 2 - α 1 α 3 ′ = α 3 - α 1 ;
步骤504:在X1O1Y1坐标系中,以P1为起始点按P1→P2→P3的顺序沿圆弧旋转;比较α′2与α′3的大小来判断圆弧的旋转方向;当α′3>α′2时,圆弧的旋转方向为逆时针;当α′3<α′2时,圆弧的旋转方向为顺时针;
步骤505:以转子提纯轴心轨迹的起始点开始,将前一段圆弧的终点作为后一段圆弧的起始点,重复执行步骤501-504;最终确定转子轴心轨迹的进动方向。
本发明的有益效果是,所提出的技术方案从局部角度出发,对轴心轨迹各个局部弧段的进动方向进行识别,不受轴心轨迹复杂程度的影响,通用性更强。
附图说明
图1是本发明提供的圆弧示意图;
图2是本发明提供的原始振动信号的波形及转子轴心轨迹图;其中,(a)是X轴方向振动信号波形图;(b)是Y轴方向振动信号波形图;(c)是原始转子轴心轨迹图;
图3是本发明提供的方法对转子轴心轨迹进动方向识别结果。
具体实施方式
下面结合附图,对优选的实施例作详细说明。应该强调的是,下述说明仅仅是示例性的,而不是为了限制本发明的范围及其应用。
图1是本发明提供的圆弧示意图。图1中,P1、P2、P3为转子提纯轴心轨迹上连续三点,其坐标分别为(x1,y1)、(x2,y2)、(x3,y3);O1为通过这三个连续的点确定的一段圆弧的圆心;XO1Y为原坐标系XOY进行平移后坐标系为;X1O1Y1为XO1Y坐标系绕点O1逆时针旋转得到的新的坐标系。
根据本发明提出的方法,在Bently RK4转子实验台上进行转子碰磨实验,碰磨发生时转子转速为1500r/min。
步骤1:设定采样频率2560Hz,转子X轴和Y轴方向的径向振动信号,将采集的电压信号保存到计算机中;数据波形及原始轴心轨迹如图2;
步骤2:对转子X轴和Y轴方向的振动信号通过快速傅立叶变换FFT得到幅值谱;
步骤3:采用相位差频谱校正法提取X轴和Y轴方向振动信号的主要频率分量;本实施例中X轴方向振动信号的主要频率成分为1倍频、
Figure BDA00003686241400000511
倍频、2倍频、3倍频和4倍频;Y轴方向振动信号的主要频率成分为1倍频、2倍频、3倍频和4倍频。
步骤4:将步骤3得到的X轴和Y轴方向各主要频率分量的幅值、频率和相位合成波形,并叠加得到X轴和Y轴方向重构信号,再利用重构信号合成转子提纯轴心轨迹;
步骤5:利用步骤4得到的转子提纯轴心轨迹确定转子轴心的进动方向的具体过程为:
步骤501:设P1、P2、P3为转子提纯轴心轨迹上连续三点,其坐标分别为(x1,y1)、(x2,y2)、(x3,y3);通过这三个连续的点确定一段圆弧,圆弧圆心为O1,圆弧圆心坐标为(a,b);
步骤502:将原坐标系XOY进行平移,平移后坐标系为XO1Y;分别计算
Figure BDA0000368624140000052
Figure BDA0000368624140000053
Figure BDA0000368624140000054
Figure BDA0000368624140000055
轴的夹角角度α1、α2、α3
步骤503:将XO1Y坐标系绕点O1逆时针旋转α1,使轴和
Figure BDA0000368624140000057
重合,即令α1=0,进而得新坐标系X1O1Y1;此时,
Figure BDA00003686241400000510
Figure BDA0000368624140000061
轴的夹角分别为 α 1 ′ = 0 α 2 ′ = α 2 - α 1 α 3 ′ = α 3 - α 1 ;
步骤504:在X1O1Y1坐标系中,以P1为起始点按P1→P2→P3的顺序沿圆弧旋转;比较α′2与α′3的大小来判断圆弧的旋转方向;当α′3>α′2时,圆弧的旋转方向为逆时针;当α′3<α′2时,圆弧的旋转方向为顺时针;
步骤505:以转子提纯轴心轨迹的起始点开始,将前一段圆弧的终点作为后一段圆弧的起始点,重复执行步骤501-504;从而确定转子轴心轨迹的进动方向;轴心轨迹进动方向识别结果如图3,图3中“○”表示正进动,“*”表示反进动。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应该涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。

Claims (3)

1.一种转子系统轴心轨迹进动方向自动识别方法,其特征在于,所述方法具体包括以下步骤:
步骤1:设定采样频率,利用振动传感器采集转子X轴和Y轴方向的径向振动信号,将采集的电压信号保存到计算机中;
步骤2:对转子X轴和Y轴方向的振动信号通过快速傅立叶变换FFT得到幅值谱;
步骤3:采用相位差频谱校正法提取X轴和Y轴方向振动信号的主要频率分量;
步骤4:将步骤3得到的X轴和Y轴方向各主要频率分量的幅值、频率和相位合成波形,并叠加得到X轴和Y轴方向重构信号,再利用重构信号合成转子提纯轴心轨迹;
步骤5:利用步骤4得到的转子提纯轴心轨迹确定转子轴心轨迹的进动方向。
2.根据权利要求1所述的方法,其特征在于,所述步骤3中,X轴和Y轴方向振动信号的主要频率分量包括信号的基频、分频及倍频。
3.根据权利要求1所述的方法,其特征在于,步骤5中,利用步骤4得到的转子提纯轴心轨迹确定转子轴心的进动方向的具体过程为:
步骤501:设P1、P2、P3为转子提纯轴心轨迹上连续三点,其坐标分别为(x1,y1)、(x2,y2)、(x3,y3);通过这三个连续的点确定一段圆弧,圆弧圆心为O1,圆弧圆心坐标为(a,b);
步骤502:将原坐标系XOY进行平移,平移后坐标系为XO1Y;分别计算
Figure FDA0000368624130000022
Figure FDA0000368624130000023
轴的夹角角度α1、α2、α3
步骤503:将XO1Y坐标系绕点O1逆时针旋转α1,使轴和
Figure FDA0000368624130000026
重合,即令α1=0,进而得新坐标系X1O1Y1;此时,
Figure FDA0000368624130000027
Figure FDA0000368624130000028
Figure FDA0000368624130000029
Figure FDA00003686241300000210
轴的夹角分别为 α 1 ′ = 0 α 2 ′ = α 2 - α 1 α 3 ′ = α 3 - α 1 ;
步骤504:在X1O1Y1坐标系中,以P1为起始点按P1→P2→P3的顺序沿圆弧旋转;比较α′2与α′3的大小来判断圆弧的旋转方向;当α′3>α′2时,圆弧的旋转方向为逆时针;当α′3<α′2时,圆弧的旋转方向为顺时针;
步骤505:以转子提纯轴心轨迹的起始点开始,将前一段圆弧的终点作为后一段圆弧的起始点,重复执行步骤501-504;从而确定转子轴心轨迹的进动方向。
CN201310361525.3A 2013-08-19 2013-08-19 一种转子系统轴心轨迹进动方向自动识别方法 Expired - Fee Related CN103412145B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310361525.3A CN103412145B (zh) 2013-08-19 2013-08-19 一种转子系统轴心轨迹进动方向自动识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310361525.3A CN103412145B (zh) 2013-08-19 2013-08-19 一种转子系统轴心轨迹进动方向自动识别方法

Publications (2)

Publication Number Publication Date
CN103412145A true CN103412145A (zh) 2013-11-27
CN103412145B CN103412145B (zh) 2014-12-17

Family

ID=49605171

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310361525.3A Expired - Fee Related CN103412145B (zh) 2013-08-19 2013-08-19 一种转子系统轴心轨迹进动方向自动识别方法

Country Status (1)

Country Link
CN (1) CN103412145B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104165686A (zh) * 2014-06-17 2014-11-26 中州大学 一种基于二元经验模态分解的转子轴心轨迹提纯方法
CN104483118A (zh) * 2014-12-08 2015-04-01 西安交通大学 基于瞬时频率轴心轨迹的转子动静碰磨故障诊断方法
CN104834936A (zh) * 2015-04-28 2015-08-12 杨露 一种基于特征有效性分析的轴心轨迹智能识别方法
CN106017956A (zh) * 2016-05-18 2016-10-12 重庆大学 基于进动能量差密度谱分析的旋转机械故障诊断新方法
CN110119764A (zh) * 2019-04-16 2019-08-13 北京天泽智云科技有限公司 一种变转速工况下轴心轨迹的提纯方法
CN113340640A (zh) * 2021-07-22 2021-09-03 格创东智(深圳)科技有限公司 旋转机械轴心轨迹提纯方法、装置、设备和存储介质
CN113358307A (zh) * 2021-06-02 2021-09-07 西安西热节能技术有限公司 一种根据轴振信号确定转子涡动方向的判断方法
CN119147783A (zh) * 2024-11-13 2024-12-17 一汽中车电驱动系统有限公司 一种转子系统轴心轨迹进动方向自动识别方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101368870A (zh) * 2008-10-13 2009-02-18 南京大学 一种用于机械转子单截面轴振分析的幅值频谱
CN101451882A (zh) * 2008-12-24 2009-06-10 南京大学 一种用于机械转子单截面轴振分析的短时幅值频谱阵列
CN102322942A (zh) * 2011-09-13 2012-01-18 大连圣力来监测技术有限公司 旋转机械振动信号矢椭谱分析方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101368870A (zh) * 2008-10-13 2009-02-18 南京大学 一种用于机械转子单截面轴振分析的幅值频谱
CN101451882A (zh) * 2008-12-24 2009-06-10 南京大学 一种用于机械转子单截面轴振分析的短时幅值频谱阵列
CN102322942A (zh) * 2011-09-13 2012-01-18 大连圣力来监测技术有限公司 旋转机械振动信号矢椭谱分析方法

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104165686A (zh) * 2014-06-17 2014-11-26 中州大学 一种基于二元经验模态分解的转子轴心轨迹提纯方法
CN104165686B (zh) * 2014-06-17 2017-10-03 中州大学 一种基于二元经验模态分解的转子轴心轨迹提纯方法
CN104483118A (zh) * 2014-12-08 2015-04-01 西安交通大学 基于瞬时频率轴心轨迹的转子动静碰磨故障诊断方法
CN104483118B (zh) * 2014-12-08 2017-04-19 西安交通大学 基于瞬时频率轴心轨迹的转子动静碰磨故障诊断方法
CN104834936B (zh) * 2015-04-28 2017-11-10 杨露 一种基于特征有效性分析的轴心轨迹智能识别方法
CN104834936A (zh) * 2015-04-28 2015-08-12 杨露 一种基于特征有效性分析的轴心轨迹智能识别方法
CN106017956A (zh) * 2016-05-18 2016-10-12 重庆大学 基于进动能量差密度谱分析的旋转机械故障诊断新方法
CN110119764A (zh) * 2019-04-16 2019-08-13 北京天泽智云科技有限公司 一种变转速工况下轴心轨迹的提纯方法
CN110119764B (zh) * 2019-04-16 2021-03-02 北京天泽智云科技有限公司 一种变转速工况下轴心轨迹的提纯方法
CN113358307A (zh) * 2021-06-02 2021-09-07 西安西热节能技术有限公司 一种根据轴振信号确定转子涡动方向的判断方法
CN113340640A (zh) * 2021-07-22 2021-09-03 格创东智(深圳)科技有限公司 旋转机械轴心轨迹提纯方法、装置、设备和存储介质
CN119147783A (zh) * 2024-11-13 2024-12-17 一汽中车电驱动系统有限公司 一种转子系统轴心轨迹进动方向自动识别方法
CN119147783B (zh) * 2024-11-13 2025-01-28 一汽中车电驱动系统有限公司 一种转子系统轴心轨迹进动方向自动识别方法

Also Published As

Publication number Publication date
CN103412145B (zh) 2014-12-17

Similar Documents

Publication Publication Date Title
CN103412145B (zh) 一种转子系统轴心轨迹进动方向自动识别方法
CN103884502B (zh) 一种变转速下风力发电机行星齿轮系统故障诊断方法
CN106769033A (zh) 基于阶次包络时频能量谱的变转速滚动轴承故障识别方法
CN102998110B (zh) 一种基于阶比-全息谱原理的旋转机械故障特征提取方法
CN106644467A (zh) 一种齿轮箱非平稳信号故障特征提取方法
CN108871742B (zh) 一种改进的无键相故障特征阶次提取方法
CN106370403A (zh) 一种基于边缘检测的瞬时频率估计方法
CN110243605A (zh) 多源时频脊线提取方法
CN110926594B (zh) 一种旋转机械信号时变频率特征提取方法
CN106863320B (zh) 一种用于智能机器人的语音交互数据获取方法及装置
BR112014017641B1 (pt) Sistema e processo de aquisição de um sinal vibratório, e, sistema de monitoração de um motor rotativo
CN104034412A (zh) 一种基于分数阶全息原理的旋转机械故障特征提取方法
CN105758404A (zh) 智能设备的实时定位方法及系统
CN105388012A (zh) 基于非线性调频小波变换的阶次跟踪方法
CN107395086A (zh) 一种内置式永磁同步电机转子初始位置检测方法
CN107359903A (zh) 扩频通信码捕获方法及装置
CN107665337A (zh) 一种轴承滚动体故障增强诊断方法
Wang et al. Tacholess order-tracking approach for wind turbine gearbox fault detection
CN105447467A (zh) 一种用户行为模式的识别系统和识别方法
CN116380451A (zh) 基于多谐波时频脊线增强提取的变转速齿轮箱故障诊断方法
CN102269803A (zh) 基于时间延迟的离散频谱低频成分的校正方法
CN104677619A (zh) 基于复信号双边谱的旋转机械故障特征提取方法
CN106053070A (zh) 基于分离信号包络谱特征的轴承滚动体故障增强诊断方法
CN102004165B (zh) 同步发电机瞬时转速测量方法
CN110346594A (zh) 一种便携式全封闭压缩机高精度转速测量方法及装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20141217

Termination date: 20160819