CN106323449B - 一种变频调速驱动下的旋转机械轴系扭振故障监测方法 - Google Patents

一种变频调速驱动下的旋转机械轴系扭振故障监测方法 Download PDF

Info

Publication number
CN106323449B
CN106323449B CN201610664399.2A CN201610664399A CN106323449B CN 106323449 B CN106323449 B CN 106323449B CN 201610664399 A CN201610664399 A CN 201610664399A CN 106323449 B CN106323449 B CN 106323449B
Authority
CN
China
Prior art keywords
shafting
coordinate system
frequency
torsional oscillation
formula
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
CN201610664399.2A
Other languages
English (en)
Other versions
CN106323449A (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.)
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN201610664399.2A priority Critical patent/CN106323449B/zh
Publication of CN106323449A publication Critical patent/CN106323449A/zh
Application granted granted Critical
Publication of CN106323449B publication Critical patent/CN106323449B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H11/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties
    • G01H11/06Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties by electric means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts

Abstract

本发明公开了一种变频调速驱动下的旋转机械轴系扭振故障监测方法,包括如下步骤:采集异步电机三相定子电流信号ia、ib、ic;对三相定子电流进行Park变换和Clarke变换,得到d‑q坐标系下的d轴和q轴定子电流id、iq;对q轴定子电流分量iq做快速傅里叶变换得到其频谱;设定轴系关键参数计算轴系扭振固有频率;计算q轴定子电流分量iq中频率在轴系扭振固有频率附近且绝对值最大的分量和频率在0Hz附近的绝对值最大的分量的比值;对该比值进行判定,超过一定阈值时视为发生扭振故障。本发明的有益效果为:从电气侧着手,对三相定子电流进行分析,监测指标简单,便于操作;不需要借助现有或者额外的机械设备,不影响转轴正常运转,适合进行长期监测。

Description

一种变频调速驱动下的旋转机械轴系扭振故障监测方法
技术领域
本发明涉及旋转机械故障诊断方法领域,尤其是一种变频调速驱动下的旋转机械轴系扭振故障监测方法。
背景技术
旋转机械在电力、能源、化工、国防等领域内发挥着无可替代的重要作用。其中,泵与风机等旋转机械多由异步电机带动。近年来,从调速节能的角度出发,变频器驱动电机已经成为一种应用日益广泛的电机驱动方式。
使用变频器驱动电机,可以通过调整电机转速适应负载变化的情况,有效减小节流损失,带来可观的经济效益。但由于异步电机由带有频率丰富的谐波的交流电进行驱动,其中频率与轴系扭振固有频率重合的谐波可能导致扭振故障。扭振故障将使得轴系承受大幅交变应力,严重影响设备安全性。在由变频器驱动异步电机带动其他旋转机械工作时,应对轴系扭振故障进行监测,确保设备安全运行。
目前监测扭振的技术主要可以分为接触式和非接触式方法。接触式扭振监测方法在转轴上架设应变片等形式的传感器,使得应变片测量轴上的剪应变,通过剪应变的变化情况来反应轴的扭振情况。接触式扭振监测方法需要再轴上额外架设设备,因此不一定适用于所有场合。非接触式扭振监测方法一般采用测齿法,利用轴上的测速齿轮盘或其他现有设备,通过测量角度的不均匀性来实现扭振监测。该方法无需在轴上另外架设设备,不影响轴的正常运行,因此适合长期监测。但该方法存在对等分机构要求很高、信号的提取和分析复杂等问题。
发明内容
本发明所要解决的技术问题在于,提供一种变频调速驱动下的旋转机械轴系扭振故障监测方法,不需要在转轴上额外架设设备,有效对轴系扭振故障进行监测,适合进行长期监测。
为解决上述技术问题,本发明提供一种变频调速驱动下的旋转机械轴系扭振故障监测方法,包括如下步骤:
(1)采集异步电机三相定子电流信号ia、ib、ic
(2)对三相定子电流进行Park变换和Clarke变换,得到d-q坐标系下的d轴和q轴定子电流id、iq
Park变换将三相a-b-c坐标系中的电流、电压、磁链等电磁量转换到静止的两相d-q坐标系中,ds和qs为静止d-q坐标系坐标轴,电磁量变换如式(1)所示:
Clarke变换将静止的两相d-q坐标系中电磁量转换到旋转的两相d-q坐标系中,坐标系旋转速度为电机的同步转速,因此该旋转坐标系也称为同步坐标系;ds和qs为静止d-q坐标系坐标轴,de和qe为旋转d-q坐标系坐标轴,θ=wt,w为同步转速,电磁量变换如式(2)所示:
将Park变换和Clarke变换相结合,将三相a-b-c坐标系中的电磁量转换到以同步转速旋转的两相d-q坐标系中,电磁量变换如式(3)所示:
(3)对q轴定子电流分量iq做快速傅里叶变换得到其频谱;
有限长离散信号x(n),n=0,1…N-1的DFT定义为:
DFT需要计算大约N2次乘法和N2次加法,当N较大时计算量很大;FFT将x(n)分成一个偶数序列和一个奇数序列的和:
x(n)=x1(n)+x2(n) (5)
这两个序列的长度都是N/2,x1(n)是偶数序列,x2(n)是奇数序列,则:
由于则有:
其中X1(k)和X2(k)分别为x1(n)和x2(n)的N/2点DFT,经过多次分解,最终将N点DFT分解为N/2个两点DFT,最终使计算量减少为次乘法和次加法;
对q轴定子电流分量iq做快速傅里叶变换FFT得到其频谱,即将iq信号分解为振幅和频率不同的分量。
(4)设定轴系关键参数计算轴系扭振固有频率;
a、将轴系模化为有附加转动惯量的等截面阶梯轴,将轴系进行详细分段,取转轴面积变化和存在附加转动惯量变化的位置为分段点,附加转动惯量指汽轮机叶片等零件;在该模型中一个等截面轴段单元的主要参数有:内径Di、外径D0、长度L、密度ρ、扭转弹性模量G以及附加转动惯量AI;汽轮机叶片是一个形状复杂的几何体,计算转动惯量时采用线性处理的有限积分法将叶片分为许多小的简单几何体,然后求出这些有限小的几何体对同一转轴的转动惯量之和;计算时,将叶片不同半径ri处的截面积fi用一个宽度为bi、半径为ri、表面积与叶片在该处的截面积及转动惯量相同的圆柱面来替代,即
式中Zd为叶片排的动叶数目;bi曲线所包络平面的转动惯量就等于叶片排的转动惯量
式中,Zd为叶片排的动叶数目;
b、轴系模化为阶梯轴后,需求出每个等截面轴段的转动惯量I′和抗扭刚度k′,其计算公式为
式中,ρ为材料密度,G为材料扭转弹性模量,D0为轴段外径,Di为轴段内径,l为轴段长度,AI为轴段上的附加转动惯量;
c、将轴系进一步模化成多段集中质量扭振模型,对n段等截面轴进行模化时,将每个等截面轴段沿着长度等分,以每个等截面轴段的抗扭刚度作为轻质弹簧的抗扭刚度,将两个相邻轴段等分线间轴段的转动惯量作为刚性轮盘的转动惯量;将n段等截面轴段组成的轴系模化为由n个理想弹簧和n+1个刚性轮盘组成的弹簧-质量块系统;
d、应用Holzer法计算扭振固有频率,首先假定一个试算频率ω,然后令θ1=1,然后按照式(12)递推计算扭角:
式中:Ii为第i段转动惯量(kg·m2),kj为扭转刚度(N·m);
最后将所有的θi代入下式计算剩余扭矩Mi
以一定的步长改变试算频率ω,代入式(13),即可得到剩余扭矩Mi与ω2的关系曲线,当Mi=0时,对应的ω值就是一个扭振固有频率。
(5)计算q轴定子电流分量iq中频率在轴系扭振固有频率附近且绝对值最大的分量和频率在0Hz附近的绝对值最大的分量的比值;
a、取q轴定子电流分量iq中频率在轴系扭振固有频率周边一定范围内,且绝对值最大的分量;优选地,该范围取为±1Hz,视实际需求可进行调整;取q轴定子电流分量iq中频率在轴系扭振固有频率周边±1Hz范围内绝对值最大的分量,记其绝对值为a;
b、取q轴定子电流分量iq中频率在0Hz周边一定范围内,且绝对值最大的分量;优选地,该范围取为±1Hz,视实际需求可进行调整;取q轴定子电流分量i.中频率在轴系扭振固有频率周边±1Hz范围内绝对值最大的分量,记其绝对值为b;
c、计算比值a/b;
(6)对该比值进行判定,超过一定阈值时视为发生扭振故障;该阈值取为1,视实际需求可进行调整。
优选的,步骤(4)中,轴系参数包括转轴铸造材料的密度、弹性模量、泊松比、各轴段长度、内外径、质量、刚度、转动惯量和阻尼。
优选的,步骤(6)中,初始阈值取为1。
本发明的有益效果为:从电气侧着手,对三相定子电流进行分析,监测指标简单,便于操作;不需要借助现有或者额外的机械设备,不影响转轴正常运转,适合进行长期监测。
附图说明
图1为本发明的方法流程示意图。
图2为本发明进行Park变换的坐标系示意图。
图3为本发明进行Clarke变换的坐标系示意图。
图4为本发明结合Park变换和Clarke变换两种变换后的坐标系示意图。
图5为本发明将带附加转动惯量的等截面阶梯轴示意图。
图6为本发明的模化的多段集中质量模型的示意图。
图7为本发明的轴系未发生扭振故障时,q轴定子电流分量iq的频谱分析图。
图8为本发明的轴系发生扭振故障时,q轴定子电流分量iq的频谱分析图。
具体实施方式
如图1所示,一种变频调速驱动下的旋转机械轴系扭振故障监测方法,包括如下步骤:
(1)采集异步电机三相定子电流信号ia、ib、ic
(2)对三相定子电流进行Park变换和Clarke变换,得到d-q坐标系下的d轴和q轴定子电流id、iq
如图2、3和4所示,Park变换将三相a-b-c坐标系中的电流、电压、磁链等电磁量转换到静止的两相d-q坐标系中,ds和qs为静止d-q坐标系坐标轴,电磁量变换如式(1)所示:
Clarke变换将静止的两相d-q坐标系中电磁量转换到旋转的两相d-q坐标系中,坐标系旋转速度为电机的同步转速,因此该旋转坐标系也称为同步坐标系;ds和qs为静止d-q坐标系坐标轴,de和qe为旋转d-q坐标系坐标轴,θ=wt,w为同步转速,电磁量变换如式(2)所示:
将Park变换和Clarke变换相结合,将三相a-b-c坐标系中的电磁量转换到以同步转速旋转的两相d-q坐标系中,电磁量变换如式(3)所示:
(3)对q轴定子电流分量iq做快速傅里叶变换得到其频谱;
快速傅里叶变换FFT算法是一种改进的傅里叶变换DFT算法,FFT在利用离散信号x(n)来计算信号X(k)的频谱时所需计算步骤远少于DFT;
有限长离散信号x(n),n=0,1…N-1的DFT定义为:
DFT需要计算大约N2次乘法和N2次加法,当N较大时计算量很大;FFT将x(n)分成一个偶数序列和一个奇数序列的和:
x(n)=x1(n)+x2(n) (5)
这两个序列的长度都是N/2,x1(n)是偶数序列,x2(n)是奇数序列,则:
由于则有:
其中X1(k)和X2(k)分别为x1(n)和x2(n)的N/2点DFT,经过多次分解,最终将N点DFT分解为N/2个两点DFT,最终使计算量减少为次乘法和次加法;
对q轴定子电流分量iq做快速傅里叶变换FFT得到其频谱,即将iq信号分解为振幅和频率不同的分量;
(4)设定轴系关键参数计算轴系扭振固有频率;轴系参数包括转轴铸造材料的密度、弹性模量、泊松比、各轴段长度、内外径、质量、刚度、转动惯量和阻尼,如图5和图6所示;
a、将轴系模化为有附加转动惯量的等截面阶梯轴,将轴系进行详细分段,取转轴面积变化和存在附加转动惯量变化的位置为分段点,附加转动惯量指汽轮机叶片等零件;在该模型中一个等截面轴段单元的主要参数有:内径Di、外径D0、长度L、密度ρ、扭转弹性模量G以及附加转动惯量AI;汽轮机叶片是一个形状复杂的几何体,计算转动惯量时采用线性处理的有限积分法将叶片分为许多小的简单几何体,然后求出这些有限小的几何体对同一转轴的转动惯量之和;计算时,将叶片不同半径ri处的截面积fi用一个宽度为bi、半径为ri、表面积与叶片在该处的截面积及转动惯量相同的圆柱面来替代,即
式中Zd为叶片排的动叶数目;bi曲线所包络平面的转动惯量就等于叶片排的转动惯量
式中,Zd为叶片排的动叶数目;
b、轴系模化为阶梯轴后,需求出每个等截面轴段的转动惯量I′和抗扭刚度k′,其计算公式为
式中,ρ为材料密度,G为材料扭转弹性模量,D0为轴段外径,Di为轴段内径,l为轴段长度,AI为轴段上的附加转动惯量;
c、将轴系进一步模化成多段集中质量扭振模型,对n段等截面轴进行模化时,将每个等截面轴段沿着长度等分,以每个等截面轴段的抗扭刚度作为轻质弹簧的抗扭刚度,将两个相邻轴段等分线间轴段的转动惯量作为刚性轮盘的转动惯量;将n段等截面轴段组成的轴系模化为由n个理想弹簧和n+1个刚性轮盘组成的弹簧-质量块系统;
d、应用Holzer法计算扭振固有频率,首先假定一个试算频率ω,然后令θ1=1,然后按照式(12)递推计算扭角:
式中:Ii为第i段转动惯量(kg·m2),kj为扭转刚度(N·m);
最后将所有的θi代入下式计算剩余扭矩Mi
以一定的步长改变试算频率ω,代入式(13),即可得到剩余扭矩Mi与ω2的关系曲线,当Mi=0时,对应的ω值就是一个扭振固有频率。
(5)计算q轴定子电流分量iq中频率在轴系扭振固有频率附近且绝对值最大的分量和频率在0Hz附近的绝对值最大的分量的比值;
a、取q轴定子电流分量iq中频率在轴系扭振固有频率周边一定范围内,且绝对值最大的分量;优选地,该范围取为±1Hz,视实际需求可进行调整;取q轴定子电流分量iq中频率在轴系扭振固有频率周边±1Hz范围内绝对值最大的分量,记其绝对值为a;
b、取q轴定子电流分量iq中频率在0Hz周边一定范围内,且绝对值最大的分量;优选地,该范围取为±1Hz,视实际需求可进行调整;取q轴定子电流分量i.中频率在轴系扭振固有频率周边±1Hz范围内绝对值最大的分量,记其绝对值为b;
c、计算比值a/b;
(6)对该比值进行判定,超过一定阈值时视为发生扭振故障;该阈值取为1,视实际需求可进行调整。
如图7和8所示,为轴系未发生扭振故障和发生扭振故障时的q轴定子电流分量iq的频谱,轴系固有扭振频率已经计算得到17.2Hz。
轴系未发生扭振故障时,取q轴定子电流分量iq中频率在17.2Hz周边±1Hz内,且绝对值最大的分量,为17.2Hz的分量,绝对值近似0。取q轴定子电流分量iq中频率在0Hz周边±1Hz内,且绝对值最大的分量。取得0.2Hz的分量,绝对值为535。计算出17.2Hz的分量的绝对值与0.2Hz的分量的绝对值的比值,接近于0,比值小于1,判定为未发生扭振故障。
轴系发生扭振故障时,取q轴定子电流分量iq中频率在17.2Hz周边±1Hz内,且绝对值最大的分量,为16.8Hz的分量,绝对值为784。取q轴定子电流分量iq中频率在0Hz周边±1Hz内,且绝对值最大的分量。取得0.35Hz的分量,绝对值为542。计算出16.8Hz的分量的绝对值与0.35Hz的分量的绝对值的比值为1.45,比值小于1,判定为发生扭振故障。
尽管本发明就优选实施方式进行了示意和描述,但本领域的技术人员应当理解,只要不超出本发明的权利要求所限定的范围,可以对本发明进行各种变化和修改。

Claims (6)

1.一种变频调速驱动下的旋转机械轴系扭振故障监测方法,其特征在于,包括如下步骤:
(1)采集异步电机三相定子电流信号ia、ib、ic
(2)对三相定子电流进行Park变换和Clarke变换,得到d-q坐标系下的d轴和q轴定子电流id、iq
(3)对q轴定子电流分量iq做快速傅里叶变换得到其频谱;
(4)设定轴系关键参数计算轴系扭振固有频率;
(5)计算q轴定子电流分量iq中频率在轴系扭振固有频率附近且绝对值最大的分量和频率在0Hz附近的绝对值最大的分量的比值;
(6)对该比值进行判定,超过一定阈值时视为发生扭振故障。
2.如权利要求1所述的变频调速驱动下的旋转机械轴系扭振故障监测方法,其特征在于,步骤(2)的具体过程为:Park变换将三相a-b-c坐标系中的电流、电压、磁链等电磁量转换到静止的两相d-q坐标系中,ds和qs为静止d-q坐标系坐标轴,电磁量变换如式(1)所示:
Clarke变换将静止的两相d-q坐标系中电磁量转换到旋转的两相d-q坐标系中,坐标系旋转速度为电机的同步转速,因此该旋转坐标系也称为同步坐标系;ds和qs为静止d-q坐标系坐标轴,de和qe为旋转d-q坐标系坐标轴,θ=wt,w为同步转速,电磁量变换如式(2)所示:
将Park变换和Clarke变换相结合,将三相a-b-c坐标系中的电磁量转换到以同步转速旋转的两相d-q坐标系中,电磁量变换如式(3)所示:
3.如权利要求1所述的变频调速驱动下的旋转机械轴系扭振故障监测方法,其特征在于,步骤(3)的具体过程为:有限长离散信号x(n),n=0,1…N-1的DFT定义为:
DFT需要计算大约N2次乘法和N2次加法,当N较大时计算量很大;FFT将x(n)分成一个偶数序列和一个奇数序列的和:
x(n)=x1(n)+x2(n) (5)
这两个序列的长度都是N/2,x1(n)是偶数序列,x2(n)是奇数序列,则:
由于则有:
其中X1(k)和X2(k)分别为x1(n)和x2(n)的N/2点DFT,经过多次分解,最终将N点DFT分解为N/2个两点DFT,最终使计算量减少为次乘法和次加法;
对q轴定子电流分量iq做快速傅里叶变换FFT得到其频谱,即将iq信号分解为振幅和频率不同的分量。
4.如权利要求1所述的变频调速驱动下的旋转机械轴系扭振故障监测方法,其特征在于,步骤(4)的具体过程为:
a、将轴系模化为有附加转动惯量的等截面阶梯轴,将轴系进行详细分段,取转轴面积变化和存在附加转动惯量变化的位置为分段点,附加转动惯量指汽轮机叶片;一个等截面轴段单元的主要参数有:内径Di、外径D0、长度L、密度ρ、扭转弹性模量G以及附加转动惯量AI;汽轮机叶片是一个形状复杂的几何体,计算转动惯量时采用线性处理的有限积分法将叶片分为许多小的简单几何体,然后求出这些有限小的几何体对同一转轴的转动惯量之和;计算时,将叶片不同半径ri处的截面积fi用一个宽度为bi、半径为ri、表面积与叶片在该处的截面积及转动惯量相同的圆柱面来替代,即
式中Zd为叶片排的动叶数目;bi曲线所包络平面的转动惯量就等于叶片排的转动惯量
式中,Zd为叶片排的动叶数目;
b、轴系模化为阶梯轴后,需求出每个等截面轴段的转动惯量I′和抗扭刚度k′,其计算公式为
式中,ρ为材料密度,G为材料扭转弹性模量,D0为轴段外径,Di为轴段内径,l为轴段长度,AI为轴段上的附加转动惯量;
c、将轴系进一步模化成多段集中质量扭振模型,对n段等截面轴进行模化时,将每个等截面轴段沿着长度等分,以每个等截面轴段的抗扭刚度作为轻质弹簧的抗扭刚度,将两个相邻轴段等分线间轴段的转动惯量作为刚性轮盘的转动惯量;将n段等截面轴段组成的轴系模化为由n个理想弹簧和n+1个刚性轮盘组成的弹簧-质量块系统;
d、应用Holzer法计算扭振固有频率,首先假定一个试算频率ω,然后令θ1=1,然后按照式(12)递推计算扭角:
式中:Ii为第i段转动惯量(kg·m2),kj为扭转刚度(N·m);
最后将所有的θi代入下式计算剩余扭矩Mi
以一定的步长改变试算频率ω,代入式(13),即可得到剩余扭矩Mi与ω2的关系曲线,当Mi=0时,对应的ω值就是一个扭振固有频率。
5.如权利要求1所述的变频调速驱动下的旋转机械轴系扭振故障监测方法,其特征在于,步骤(4)中,轴系参数包括转轴铸造材料的密度、弹性模量、泊松比、各轴段长度、内外径、质量、刚度、转动惯量和阻尼。
6.如权利要求1所述的变频调速驱动下的旋转机械轴系扭振故障监测方法,其特征在于,步骤(6)中,初始阈值取为1。
CN201610664399.2A 2016-08-12 2016-08-12 一种变频调速驱动下的旋转机械轴系扭振故障监测方法 Active CN106323449B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610664399.2A CN106323449B (zh) 2016-08-12 2016-08-12 一种变频调速驱动下的旋转机械轴系扭振故障监测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610664399.2A CN106323449B (zh) 2016-08-12 2016-08-12 一种变频调速驱动下的旋转机械轴系扭振故障监测方法

Publications (2)

Publication Number Publication Date
CN106323449A CN106323449A (zh) 2017-01-11
CN106323449B true CN106323449B (zh) 2019-02-01

Family

ID=57740364

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610664399.2A Active CN106323449B (zh) 2016-08-12 2016-08-12 一种变频调速驱动下的旋转机械轴系扭振故障监测方法

Country Status (1)

Country Link
CN (1) CN106323449B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109917170B (zh) * 2019-04-04 2020-03-10 西南交通大学 一种单相脉冲整流器无锁相环的dq电流检测方法
CN114114068A (zh) * 2020-08-28 2022-03-01 台达电子工业股份有限公司 马达连接故障检测方法
CN112857723B (zh) * 2021-01-21 2022-01-11 华中科技大学 一种转轴扭振的实验测试方法及装置
CN114545228B (zh) * 2022-01-26 2022-11-22 图湃(北京)医疗科技有限公司 基于电机电流反馈的oct系统滚珠丝杠磨损状态监测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4148222A (en) * 1977-12-22 1979-04-10 General Electric Company Apparatus and method for measuring torsional vibration
CN1695048A (zh) * 2002-09-10 2005-11-09 阿尔斯托姆科技有限公司 检测在电机上轴系振动的方法和装置
CN101603855A (zh) * 2009-07-16 2009-12-16 华北电力大学 汽轮发电机组轴系扭振的分析方法
CN104236704A (zh) * 2014-09-19 2014-12-24 华北电力大学 用于汽轮发电机组轴系次同步振荡扭振监测的方法及系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4148222A (en) * 1977-12-22 1979-04-10 General Electric Company Apparatus and method for measuring torsional vibration
CN1695048A (zh) * 2002-09-10 2005-11-09 阿尔斯托姆科技有限公司 检测在电机上轴系振动的方法和装置
CN101603855A (zh) * 2009-07-16 2009-12-16 华北电力大学 汽轮发电机组轴系扭振的分析方法
CN104236704A (zh) * 2014-09-19 2014-12-24 华北电力大学 用于汽轮发电机组轴系次同步振荡扭振监测的方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
变频电机驱动轴系扭矩脉动特性试验研究;李年仔 等;《动力工程学报》;20140930;第34卷(第9期);第731-735页

Also Published As

Publication number Publication date
CN106323449A (zh) 2017-01-11

Similar Documents

Publication Publication Date Title
CN106323449B (zh) 一种变频调速驱动下的旋转机械轴系扭振故障监测方法
Gong et al. Current-based mechanical fault detection for direct-drive wind turbines via synchronous sampling and impulse detection
CN105262394A (zh) 一种内置式永磁同步电机的mtpa控制方法及其控制系统
CN104865400A (zh) 一种风电机组转速的检测识别方法及系统
CN110046418B (zh) 一种永磁电机周期定子的振动特性分析方法
CN104935222B (zh) 一种无刷双馈感应发电机转速估计系统
Shang et al. The effects of stator and rotor eccentricities on measurement accuracy of axial flux variable-reluctance resolver with sinusoidal rotor
Chen et al. Proportional selection scheme: A frequency band division tool for rolling element bearing diagnostics
CN102564765A (zh) 一种汽轮发电机组动静碰摩轴向位置定位方法
Sheng et al. Blade aerodynamic asymmetry fault analysis and diagnosis of wind turbines with doubly fed induction generator
CN104935229A (zh) 储能涡卷弹簧实时转动惯量的获取方法
Ye et al. Bearing fault diagnosis under time-varying speed and load conditions via speed sensorless algorithm and angular resample
Huang et al. Instantaneous amplitude-frequency feature extraction for rotor fault based on BEMD and Hilbert transform
CN108258949A (zh) 一种无感永磁同步电机的转子位置自适应估计方法
CN112327957B (zh) 转子低频振动多阶线谱的控制方法及装置
CN101510693B (zh) 一种获得次同步振荡转速信号的方法和装置
CN114112396B (zh) 一种采用电流信号分析的转速波动准平稳工况下轴承故障诊断方法
Xiang et al. Torsional vibration of a shafting system under electrical disturbances
CN105978430A (zh) 一种电机电阻自整定方法及系统
CN113037162B (zh) 无轴承永磁同步电机神经网络带通滤波器振动补偿控制器
Jiang et al. Electromagnetic Force and Mechanical Response of Turbo-Generator End Winding under Electromechanical Faults
Chen et al. A study on automatic machine condition monitoring and fault diagnosis for bearing and unbalanced rotor faults
CN102420434B (zh) 一种交流系统次同步谐振分析方法
Zhou et al. Application of intermediate bearing failure diagnosis based on computational order analysis
CN104931129A (zh) 一种发电机组转子扭振测量系统及测量方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: 210009 No. 87 Dingjiaqiao, Gulou District, Nanjing City, Jiangsu Province

Applicant after: Southeast University

Address before: No. 2, four archway in Xuanwu District, Nanjing, Jiangsu

Applicant before: Southeast University

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant