CN104215395B - 一种转子不平衡故障的检测方法及装置 - Google Patents
一种转子不平衡故障的检测方法及装置 Download PDFInfo
- Publication number
- CN104215395B CN104215395B CN201410455745.7A CN201410455745A CN104215395B CN 104215395 B CN104215395 B CN 104215395B CN 201410455745 A CN201410455745 A CN 201410455745A CN 104215395 B CN104215395 B CN 104215395B
- Authority
- CN
- China
- Prior art keywords
- rotor
- unbalance
- characteristic parameter
- function model
- fitting
- 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
Abstract
本发明提供了一种转子不平衡故障的检测方法及装置。所述包括:采集转子振动信号;根据转子振动信号,提取转子的特征参数并分析每个特征参数对于不平衡量的敏感度,提取出N个敏感的特征参数,作为N个不平衡特征参数;采用最佳拟合方法,对N个不平衡特征参数进行拟合,得到转子的N个不平衡特征参数的N个拟合函数模型;将N个转子不平衡特征参数的数值代入相应的拟合函数模型,计算对应的N个不平衡率;对N个不平衡率取平均值,生成最终的转子不平衡率;将转子不平衡率乘以转子质量,生成转子不平衡量。本发明使得转子在不停机的情况下即可获知转子的不平衡量,可有效减少停机检测带来的经济损失,同时减少因维护、校正不及时带来的安全隐患。
Description
技术领域
本发明涉及机械的检测领域,尤其涉及一种对转子的不平衡故障进行检测的方法及装置,具体的讲是一种在不停机的状态下对转子的不平衡故障进行检测的方法及装置。
背景技术
在石油和化工企业中,旋转机械的转子出现不平衡故障的几率非常高,约占旋转机械故障的30%以上。转子出现不平衡故障时,将引起整个机器的振动,并产生噪声,加速轴承的磨损,降低机器的寿命,甚至使机器控制失灵,发生严重事故。引起转子不平衡的原因是多方面的,例如转子的结构设计不合理、机械加工质量偏差、装配误差、材质不均匀、动平衡精度差;运行中联轴器相对位置的改变;转子部件缺损等等。
为了确知转子是否存在不平衡故障,通常解决方法是将设备停机,用专用检测仪器对转子系统进行不平衡测量。如果不平衡量超标,则需要对其进行维修或校正。但是对设备停机,会使生产中断,严重影响经济效益,因此通常在不平衡故障很严重时才被迫停下来检测,这时就往往错过了最佳的维修和校正时机,增加了维修费用,也给设备安全运行带来了严重的隐患。
发明内容
本发明的目的是提供一种转子不平衡故障的检测方法及装置,用以在不停机的情况下定量评估转子不平衡量。
为了达到上述目的,本发明实施例提供一种转子不平衡故障的检测方法,包括:采集转子振动信号;根据所述转子振动信号,提取转子的特征参数并分析每个所述特征参数对于不平衡量的敏感度,提取出N个敏感的特征参数,作为转子的N个不平衡特征参数;采用最佳拟合方法,对所述N个不平衡特征参数进行拟合,得到所述转子的N个不平衡特征参数的N个拟合函数模型;将所述N个转子不平衡特征参数的数值代入相应的拟合函数模型,计算对应的N个不平衡率;对所述N个不平衡率取平均值,生成最终的转子不平衡率;对所述转子不平衡率,乘以转子质量,即得到转子不平衡量。
为了达到上述目的,本发明实施例还提供一种转子不平衡故障的检测装置,包括:转子振动信号采集单元,用于采集转子振动信号;不平衡特征参数提取单元,用于根据所述转子振动信号,提取转子的特征参数并分析每个所述特征参数对于不平衡量的敏感度,提取出N个敏感的特征参数,作为转子的N个不平衡特征参数;拟合单元,用于采用最佳拟合方法,对所述N个不平衡特征参数进行拟合,得到所述转子的N个不平衡特征参数的N个拟合函数模型;不平衡率计算单元,用于将所述N个转子不平衡特征参数的数值代入相应的拟合函数模型,计算对应的N个不平衡率,并对所述N个不平衡率取平均值,生成最终的转子不平衡率;不平衡量计算单元,用于对所述转子不平衡率,乘以转子质量,生成转子不平衡量。
本发明实施例的转子不平衡故障的检测方法及装置的有益效果是:通过转子的不停机检测,确定转子不平衡的定量估算模型,从而使得转子在不停机的情况下即可获知转子的不平衡量。该方法可以及时地获得转子的平衡状态,为转子的维护、保养和校正提供依据。可以有效地减少停机检测带来的经济损失,同时减少因维护、校正不及时带来的安全隐患。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的转子不平衡故障的检测方法的流程图;
图2为本发明实施例的转子不平衡故障的检测装置的结构示意图;
图3为本发明具体实施例的1倍频随配重的变化趋势图;
图4为本发明具体实施例的0.5倍频随配重的变化趋势图;
图5为本发明具体实施例的2倍频与3倍频随配重的变化趋势图;
图6为本发明具体实施例的峰峰值随配重的变化趋势图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提供了一种在设备不停机情况下的对转子的不平衡故障进行检测的方法及装置,在设备正常运行中采集转子的振动信号,提取不平衡特征参数,然后用拟合方法挖掘特征参数与不平衡严重程度的函数关系和规律,从而定量估算转子的不平衡量。
图1为本发明实施例的转子不平衡故障的检测方法的流程图。如图1所示,所述检测方法包括:步骤S101,采集转子振动信号;步骤S102,根据所述转子振动信号,提取转子的特征参数并分析每个所述特征参数对于不平衡量的敏感度,提取出N个敏感的特征参数,作为转子的N个不平衡特征参数;步骤S103,采用最佳拟合方法,对所述N个不平衡特征参数进行拟合,得到所述转子的N个不平衡特征参数的N个拟合函数模型;步骤S104,将所述N个转子不平衡特征参数的数值代入相应的拟合函数模型,计算对应的N个不平衡率;步骤S105,对所述N个不平衡率取平均值,生成最终的转子不平衡率;步骤S106,计算不平衡量,不平衡量=不平衡率×转子质量。
在本实施例的步骤S101中,可以通过振动加速度传感器以及MDES(MachineDiagnosis Expert System,机器诊断专家系统)-5数据采集器采集转子振动信号。所述提取的转子的特征参数包括振动信号的时域特征参数和频域特征参数,其中:所述时域特征参数包括绝对值、峰峰值、方差值、有效值、偏斜度和方根幅值;所述频域特征参数包括1倍频幅值、2倍频幅值、3倍频幅值以及半倍频幅值。
在本实施例的步骤S102中,分析所述每个特征参数对于不平衡量的敏感度,包括:根据每一特征参数随不平衡量变化的趋势,分析所述每个特征参数对于不平衡量的敏感度。敏感度定义为:指标变化量/不平衡变化量,比值越大,敏感度越高。它是基于实验得出的结果,因为不同特征参数对故障敏感度是不一样的,选择对故障敏感度高的参数会提高故障诊断的准确率,因此通过大量实验模拟最终获得对故障敏感的故障特征。
在本实施例的步骤S103中,采用最佳拟合方法对所述N个不平衡特征参数进行拟合,得到所述转子的N个不平衡特征参数的N个拟合函数模型。所谓最佳拟合就是拟合结果能够在实现最大相关的同时还能反映因变量随自变量改变的变化趋势。为使估算结果更加符合设备的实际情况,提出了“不平衡率λ”这一无量纲概念。其拟合步骤如下:
步骤1,确定自变量和因变量:将提取出来的不平衡特征参数作为拟合函数模型的自变量,将不平衡量质量与转子质量之比Δ作为拟合函数的因变量
步骤2,选定拟合函数模型:一般四阶函数表达式足以表达上述因变量与自变量之间的关系,所以选定λ=a+b·x+c·x2+d·x3+e·x4为拟合结果的函数模型。式中因变量λ为不平衡率,它与Δ同为无量纲变量。
步骤3,拟合函数模型中系数的求解:用编程软件(1stOpt)编写曲线拟合程序代码;选定1stOpt软件中的共轭梯度法作为求取待定系数的方法;将Δ作为因变量、第1个特征参量为自变量代入程序代码中,获得系数a1、b1、c1、d1、e1,再代入拟合函数模型中,得到λ1=a1+b1·x+c1·x2+d1·x3+e1·x4;同理,将Δ作为因变量、第i个特征参量为自变量代入到程序代码中,得到λi=ai+bi·x+ci·x2+di·x3+ei·x4。最终获得全部的拟合公式。本步骤用均方差、相关系数和决定系数3个指标判断拟合效果,均方差越小,相关系数和决定系数越大,表示拟合结果越好。
在本实施例的步骤S104和S105中,在额定转速(本实施例为8610rpm)下进行振动信号采集,提取的特征参数分别代入相应的拟合函数模型,即N个λ=a+b·x+c·x2+d·x3+e·x4类型的多项式,其中N的个数取决于不平衡特征参数选取的个数,计算出对应的λ1,λ2,……,λN。并对所述N个不平衡率取平均值,生成最终的转子不平衡率λ=(λ1+λ2+……+λN)/N。
根据对转子不平衡故障的定量研究,确定了转子振动信号和不平衡量之间的一一对应关系。同时,将理论公式和实际设备检测的特征参数结合,可得到实际设备转子不平衡量与特征参数之间的关系,进而通过检测振动信号计算出转子的不平衡量。同时,考虑到现场设备功率、结构、转速等因素也会对不平衡率λ的造成影响,因此,需要进一步验证并增添修正系数k,最终得到完整的不平衡率估算公式:
λ=(λ1+λ2+……+λN)/N。
在本实施例的步骤S106中,计算不平衡量:不平衡量=λ×转子质量。
图2为本发明实施例的转子不平衡故障的检测装置的结构示意图。如图所示,本发明的检测装置包括:
转子振动信号采集单元101,用于采集转子振动信号;不平衡特征参数提取单元102,用于根据所述转子振动信号,提取转子的特征参数并分析每个所述特征参数对于不平衡量的敏感度,提取出N个敏感的特征参数,作为转子的N个不平衡特征参数;拟合单元103,用于采用最佳拟合方法,对所述N个不平衡特征参数进行拟合,得到所述转子的N个不平衡特征参数的N个拟合函数模型;不平衡率计算单元104,用于将所述N个转子不平衡特征参数的数值代入相应的拟合函数模型,计算对应的N个不平衡率,并对所述N个不平衡率取平均值,生成最终的转子不平衡率;不平衡量计算单元105,用于对所述转子不平衡率,乘以转子质量,生成转子不平衡量。
在本实施例中,所述转子振动信号采集单元101通过振动加速度传感器以及MDES-5数据采集器进行采集。
在本实施例中,所述不平衡特征参数提取单元102提取的转子的特征参数包括振动信号的时域特征参数和频域特征参数,其中:所述时域特征参数包括绝对值、峰峰值、方差值、有效值、偏斜度和方根幅值;所述频域特征参数包括1倍频幅值、2倍频幅值、3倍频幅值以及半倍频幅值。
在本实施例中,所述不平衡特征参数提取单元102分析所述每个特征参数对于不平衡量的敏感度,包括:根据每一特征参数随不平衡量变化的趋势,分析所述每个特征参数对于不平衡量的敏感度。
在本实施例中,所述拟合单元103采用最佳拟合方法,对所述N个不平衡特征参数进行拟合,得到所述转子的N个不平衡特征参数的N个拟合函数模型,包括:
步骤1,确定自变量和因变量:将所述不平衡特征参数作为所述拟合函数模型的自变量,将不平衡量质量与转子质量之比Δ作为拟合函数的因变量
步骤2,选定拟合函数模型:一般四阶函数表达式足以表达所述因变量与自变量之间的关系,所以选定λ=a+b·x+c·x2+d·x3+e·x4为拟合结果的函数模型,式中因变量λ为所述不平衡率,它与Δ同为无量纲变量;
步骤3,拟合函数模型中系数的求解:用编程软件1stOpt编写曲线拟合程序代码;选定1stOpt软件中的共轭梯度法作为求取待定系数的方法;将Δ作为因变量、将第1个特征参数作为自变量代入程序代码中,获得系数a1、b1、c1、d1、e1,再代入拟合函数模型中,得到λ1=a1+b1·x+c1·x2+d1·x3+e1·x4;同理,将Δ作为因变量、将第i个特征参数作为自变量代入到程序代码中,得到λi=ai+bi·x+ci·x2+di·x3+ei·x4;最终获得全部的拟合函数模型。
在本实施例中,所述不平衡率计算单元104将所述N个转子不平衡特征参数的数值代入相应的拟合函数模型,计算对应的N个不平衡率,包括:将提取的特征参数分别代入相应的拟合函数式,计算出对应的λ1,λ2,……,λN,并对所述N个不平衡率取平均值,生成最终的转子不平衡率λ=(λ1+λ2+……+λN)/N。
在本实施例中,所述不平衡量计算单元105对所述转子不平衡率,乘以转子质量,生成转子不平衡量:不平衡量=λ×转子质量。
具体实施例:
首先建立转子实验台,对转子不平衡程度进行定量的模拟和测算,实验采用本特利RK4转子实验台。转子实验台硬件系统包括加速度传感器、MDES-5数据采集器、便携式计算机及连接电缆,转子质量为500g。
(1)数据采集
采用MDES-5数据采集器在转子不停机情况下进行振动信号的采集和处理。测点统一为靠近联轴器端轴承座的垂直方向,在每种不平衡量情况下转速由591rpm到8610rpm,转速间隔约为500rpm,共计18种转速。在每一种转速下,不平衡量从0.1g到1.1g每0.1g采集一组数据,共计11组。因此总共采集了18×11=198组数据,如表1所示。
表1转子不平衡实验数据集
表1所示为所做实验保存的振动信号文件名。依照所选转速和不平衡量大小的不同进行采集。以i1g为例,表示转速为591rpm,不平衡量质量为0.1g情况下所采集的数据文件。
(2)特征提取
提取振动信号的特征值,包括0.5倍频幅值(简称0.5倍频),1倍频幅值(简称1倍频),2倍频幅值(简称2倍频),3倍频幅值(简称3倍频),峰峰值,有效值,方根幅值。首先分析每个特征参数对于不平衡量的敏感度,剔除不敏感的不平衡量。最终确定用于拟合的特征参数为:0.5倍频,1倍频,2倍频,3倍频,峰峰值。其相应幅值大小详细见表2,各参数随不平衡量变化趋势如图3至图6。
表2转速8610rpm下振动信号的5个拟合特征参数幅值
(3)拟合函数模型
用编程软件(1stOpt)编写曲线拟合程序代码;选定1stOpt软件中的共轭梯度法作为求取待定系数的方法;将表2中的Δ作为因变量、第1个特征参量(0.5倍频)作为自变量代入程序代码中,获得系数a1、b1、c1、d1、e1,再代入拟合函数模型中,得到λ1=-0.00088+0.00172·x-0.00046·x2+4.84684·x3-1.70883·x4。表3为1倍频幅值拟合结果。
表3 8610rpm转速下1倍频幅值拟合结果(单位:m/s2)
从表中可以看出,均方差接近于0,相关系数和决定系数都接近于1,表示拟合结果较好。
均方差、相关系数和决定系数,这三个指标是表示实验中计算结果和实际值的相关程度,是判定拟合程度好坏的标准,均方差越小,相关系数和决定系数越大,表示拟合结果越好。
同理,得到全部的拟合函数模型,如下:
1倍频:λ1=-0.00088+0.00172·x-0.00046·x2+4.84684·x3-1.70883·x4;
2倍频:λ2=0.00348-0.03462·x+0.14944·x2-0.2501·x3+0.14420·x4;
3倍频:λ3=0.00022+0.00638·x-0.00843·x2+0.00442·x3-0.00082·x4;
0.5倍频:λ4=0.00365-0.28217·x+9.13921·x2-115.73749·x3+522.04335·x4;
峰峰值:λ5=-0.00634+0.00089·x-3.5996·x2+5.858·x3-3.2734·x4;
(4)计算不平衡率
最终得到实验室转子不平衡率λ的估算表达式为:
λ=(λ1+λ2+λ3+λ4+λ5)/5。
为了验证上述公式的有效性,对不同不平衡量下转子实验台数据进行验证。如不平衡量为1.0g时,1倍频幅值为12.231m/s2,2倍频幅值为0.085m/s2,3倍频幅值为0.396m/s2,0.5倍频幅值为0.079m/s2,峰峰值为60.768m/s2,分别代入相应的拟合函数模型,得到λ1=0.0019,λ2=0.0015,λ3=0.0017,λ4=0.0017,λ5=0.0019,经计算得λ=0.00174,与实际不平衡量Δ相比误差为13%(Δ=1/500=0.002)。同理,对其他实验数据进行验证,结果如表4,整体实验误差为8.819%,误差都没有超过20%,与实际吻合较好。
表4转速8610rpm下振动信号的5个特征参数值的不平衡率
(5)计算不平衡量
不平衡量=λ×转子质量=0.00174×500=0.87g。
根据对转子不平衡故障的定量研究,确定了特征参数和不平衡量之间的一一对应关系。从表4中,可以看到随着不平衡率增大,设备振动也在逐渐增大,符合不平衡的故障机理,也证明了在不停机情况下利用定量评估的方法判断转子不平衡量的正确性。将理论公式和实际设备检测的特征参数结合,可得到实际设备转子不平衡量与特征参数之间的关系,进而通过检测振动信号估算出转子的不平衡量,为转子维护和校正提供依据。
本发明实施例的转子不平衡故障的检测方法及装置,通过转子的不停机检测,确定转子不平衡的定量估算模型,从而使得转子在不停机的情况下即可获知转子的不平衡量。该方法可以及时地获得转子的平衡状态,为转子的维护、保养和校正提供依据。可以有效地减少停机检测带来的经济损失,同时减少因维护、校正不及时带来的安全隐患。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (14)
1.一种转子不平衡故障的检测方法,其特征在于,所述检测方法包括:
采集转子振动信号;
根据所述转子振动信号,提取转子的特征参数并分析每个所述特征参数对于不平衡量的敏感度,提取出N个敏感的特征参数,作为转子的N个不平衡特征参数;
采用最佳拟合方法,对所述N个不平衡特征参数进行拟合,得到所述转子的N个不平衡特征参数的N个拟合函数模型;
将所述N个转子不平衡特征参数的数值代入相应的拟合函数模型,计算对应的N个不平衡率;
对所述N个不平衡率取平均值,生成最终的转子不平衡率;
对所述转子不平衡率,乘以转子质量,生成转子不平衡量。
2.根据权利要求1所述的转子不平衡故障的检测方法,其特征在于,所述采集转子振动信号通过振动加速度传感器以及MDES-5数据采集器进行采集。
3.根据权利要求1所述的转子不平衡故障的检测方法,其特征在于,所述提取的转子的特征参数包括振动信号的时域特征参数和频域特征参数,其中:
所述时域特征参数包括绝对值、峰峰值、方差值、有效值、偏斜度和方根幅值;
所述频域特征参数包括1倍频幅值、2倍频幅值、3倍频幅值以及半倍频幅值。
4.根据权利要求1所述的转子不平衡故障的检测方法,其特征在于,分析所述每个特征参数对于不平衡量的敏感度,包括:
根据每一特征参数随不平衡量变化的趋势,分析所述每个特征参数对于不平衡量的敏感度。
5.根据权利要求1所述的转子不平衡故障的检测方法,其特征在于,采用最佳拟合方法,对所述N个不平衡特征参数进行拟合,得到所述转子的N个不平衡特征参数的N个拟合函数模型,包括:
步骤1,确定自变量和因变量:将所述不平衡特征参数作为所述拟合函数模型的自变量,将不平衡量质量与转子质量之比Δ作为拟合函数的因变量
步骤2,选定拟合函数模型:一般四阶函数表达式足以表达所述因变量与自变量之间的关系,所以选定λ=a+b·x+c·x2+d·x3+e·x4为拟合结果的函数模型,式中因变量λ为所述不平衡率,它与Δ同为无量纲变量;
步骤3,拟合函数模型中系数的求解:用编程软件1stOpt编写曲线拟合程序代码;选定1stOpt软件中的共轭梯度法作为求取待定系数的方法;将Δ作为因变量、将第1个特征参数作为自变量代入程序代码中,获得系数a1、b1、c1、d1、e1,再代入拟合函数模型中,得到λ1=a1+b1·x+c1·x2+d1·x3+e1·x4;同理,将Δ作为因变量、将第i个特征参数作为自变量代入到程序代码中,得到λi=ai+bi·x+ci·x2+di·x3+ei·x4;最终获得全部的拟合函数模型。
6.根据权利要求5所述的转子不平衡故障的检测方法,其特征在于,将所述N个转子不平衡特征参数的数值代入相应的拟合函数模型,计算对应的N个不平衡率,包括:
将提取的特征参数分别代入相应的拟合函数式,计算出对应的λ1,λ2,……,λN,并对所述N个不平衡率取平均值,生成最终的转子不平衡率λ=(λ1+λ2+……+λN)/N。
7.根据权利要求6所述的转子不平衡故障的检测方法,其特征在于,对所述转子不平衡率,乘以转子质量,生成转子不平衡量:不平衡量=λ×转子质量。
8.一种转子不平衡故障的检测装置,其特征在于,所述检测装置包括:
转子振动信号采集单元,用于采集转子振动信号;
不平衡特征参数提取单元,用于根据所述转子振动信号,提取转子的特征参数并分析每个所述特征参数对于不平衡量的敏感度,提取出N个敏感的特征参数,作为转子的N个不平衡特征参数;
拟合单元,用于采用最佳拟合方法,对所述N个不平衡特征参数进行拟合,得到所述转子的N个不平衡特征参数的N个拟合函数模型;
不平衡率计算单元,用于将所述N个转子不平衡特征参数的数值代入相应的拟合函数模型,计算对应的N个不平衡率,并对所述N个不平衡率取平均值,生成最终的转子不平衡率;
不平衡量计算单元,用于对所述转子不平衡率,乘以转子质量,生成转子不平衡量。
9.根据权利要求8所述的转子不平衡故障的检测装置,其特征在于,所述转子振动信号采集单元通过振动加速度传感器以及MDES-5数据采集器进行采集。
10.根据权利要求8所述的转子不平衡故障的检测装置,其特征在于,所述不平衡特征参数提取单元提取的转子的特征参数包括振动信号的时域特征参数和频域特征参数,其中:
所述时域特征参数包括绝对值、峰峰值、方差值、有效值、偏斜度和方根幅值;
所述频域特征参数包括1倍频幅值、2倍频幅值、3倍频幅值以及半倍频幅值。
11.根据权利要求8所述的转子不平衡故障的检测装置,其特征在于,所述不平衡特征参数提取单元分析所述每个特征参数对于不平衡量的敏感度,包括:
根据每一特征参数随不平衡量变化的趋势,分析所述每个特征参数对于不平衡量的敏感度。
12.根据权利要求8所述的转子不平衡故障的检测装置,其特征在于,所述拟合单元采用最佳拟合方法,对所述N个不平衡特征参数进行拟合,得到所述转子的N个不平衡特征参数的N个拟合函数模型,包括:
步骤1,确定自变量和因变量:将所述不平衡特征参数作为所述拟合函数模型的自变量,将不平衡量质量与转子质量之比Δ作为拟合函数的因变量
步骤2,选定拟合函数模型:一般四阶函数表达式足以表达所述因变量与自变量之间的关系,所以选定λ=a+b·x+c·x2+d·x3+e·x4为拟合结果的函数模型,式中因变量λ为所述不平衡率,它与Δ同为无量纲变量;
步骤3,拟合函数模型中系数的求解:用编程软件1stOpt编写曲线拟合程序代码;选定1stOpt软件中的共轭梯度法作为求取待定系数的方法;将Δ作为因变量、将第1个特征参数作为自变量代入程序代码中,获得系数a1、b1、c1、d1、e1,再代入拟合函数模型中,得到λ1=a1+b1·x+c1·x2+d1·x3+e1·x4;同理,将Δ作为因变量、将第i个特征参数作为自变量代入到程序代码中,得到λi=ai+bi·x+ci·x2+di·x3+ei·x4;最终获得全部的拟合函数模型。
13.根据权利要求8所述的转子不平衡故障的检测装置,其特征在于,所述不平衡率计算单元将所述N个转子不平衡特征参数的数值代入相应的拟合函数模型,计算对应的N个不平衡率,包括:
将提取的特征参数分别代入相应的拟合函数式,计算出对应的λ1,λ2,……,λN,并对所述N个不平衡率取平均值,生成最终的转子不平衡率λ=(λ1+λ2+……+λN)/N。
14.根据权利要求12所述的转子不平衡故障的检测装置,其特征在于,所述不平衡量计算单元对所述转子不平衡率,乘以转子质量,生成转子不平衡量:不平衡量=λ×转子质量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410455745.7A CN104215395B (zh) | 2014-09-09 | 2014-09-09 | 一种转子不平衡故障的检测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410455745.7A CN104215395B (zh) | 2014-09-09 | 2014-09-09 | 一种转子不平衡故障的检测方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104215395A CN104215395A (zh) | 2014-12-17 |
CN104215395B true CN104215395B (zh) | 2016-11-23 |
Family
ID=52097118
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410455745.7A Active CN104215395B (zh) | 2014-09-09 | 2014-09-09 | 一种转子不平衡故障的检测方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104215395B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105021351B (zh) * | 2015-08-16 | 2017-09-26 | 北京航空航天大学 | 一种基于光纤光栅振动传感系统的直升机桨叶质量不平衡检测方法 |
CN105510010A (zh) * | 2015-11-26 | 2016-04-20 | 江南大学 | 一种针对旋转机械构造异常为不对中的特征参量模型 |
CN105756864B (zh) * | 2016-04-15 | 2018-03-13 | 浙江运达风电股份有限公司 | 基于双馈风力发电机组定子电流的桨叶不平衡故障诊断方法 |
CN106227909B (zh) * | 2016-06-17 | 2019-09-06 | 华电电力科学研究院 | 一种一次消除汽轮发电机组转子动态挠曲的方法 |
CN107025348B (zh) * | 2017-04-11 | 2020-04-03 | 西北工业大学 | 航空发动机转子早期积液故障识别方法 |
CN110276164B (zh) * | 2019-07-10 | 2023-04-07 | 湖南科技大学 | 环形零部件筋板焊接错位量引起转子结构不平衡量分析法 |
CN111579153A (zh) * | 2020-05-19 | 2020-08-25 | 深圳市唯真电机发展有限公司 | 一种转子不平衡量测试装置 |
CN112653306B (zh) * | 2020-12-18 | 2022-08-23 | 重庆忽米网络科技有限公司 | 基于工业ai的电动机动平衡工艺优化方法及系统 |
CN113295412B (zh) * | 2021-05-26 | 2022-10-11 | 华能澜沧江水电股份有限公司 | 一种检测立式水轮发电机组导轴承受力不平衡原因的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101451898A (zh) * | 2009-01-13 | 2009-06-10 | 华北电力大学 | 汽轮发电机组转子热弯曲不平衡故障实时诊断方法 |
CN101625277A (zh) * | 2008-07-07 | 2010-01-13 | 西门子公司 | 不平衡状态定量检测方法和装置及工件装夹状态检测方法 |
CN102353500A (zh) * | 2011-06-03 | 2012-02-15 | 上海师范大学 | 一种用于动平衡测量的不平衡信号提取方法 |
JP5212789B2 (ja) * | 2008-03-21 | 2013-06-19 | 株式会社Ihi | 回転体のアンバランス修正方法及び装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010151643A (ja) * | 2008-12-25 | 2010-07-08 | Yokohama Rubber Co Ltd:The | 回転体への挿入物のアンバランス測定方法及びアンバランス修正方法 |
-
2014
- 2014-09-09 CN CN201410455745.7A patent/CN104215395B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5212789B2 (ja) * | 2008-03-21 | 2013-06-19 | 株式会社Ihi | 回転体のアンバランス修正方法及び装置 |
CN101625277A (zh) * | 2008-07-07 | 2010-01-13 | 西门子公司 | 不平衡状态定量检测方法和装置及工件装夹状态检测方法 |
CN101451898A (zh) * | 2009-01-13 | 2009-06-10 | 华北电力大学 | 汽轮发电机组转子热弯曲不平衡故障实时诊断方法 |
CN102353500A (zh) * | 2011-06-03 | 2012-02-15 | 上海师范大学 | 一种用于动平衡测量的不平衡信号提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104215395A (zh) | 2014-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104215395B (zh) | 一种转子不平衡故障的检测方法及装置 | |
CN109655200B (zh) | 一种风力发电机组风轮不平衡的诊断方法及系统 | |
CN103608739A (zh) | 振动监测 | |
JP6849446B2 (ja) | 振動監視システム | |
CN103776480A (zh) | 基于多次移动平均的微小故障检测方法和装置 | |
Kusnick et al. | Wind turbine rotor imbalance detection using nacelle and blade measurements | |
CN102288283A (zh) | 汽轮发电机组振动故障的监测方法及装置 | |
CN110579312B (zh) | 一种无试重旋转机械多轮盘轴系动平衡故障检测方法 | |
CN106762343B (zh) | 基于在线数据的水轮发电机组推力轴承故障的诊断方法 | |
CN106052743B (zh) | 一种评估传感器质量对频响函数影响大小的方法 | |
CN110823576A (zh) | 基于生成对抗网络的机械异常检测方法 | |
CN114216640A (zh) | 用于检测工业设备故障状态的方法、设备和介质 | |
WO2011032606A1 (en) | Method for balancing a rotor mounted on a hub of a wind turbine | |
CN103900763B (zh) | 轮胎多级轮辋不平衡量标定方法 | |
CN104330214B (zh) | 一种简易现场动平衡方法 | |
CN106596110B (zh) | 基于在线数据的水轮发电机组水力不平衡故障的自动分析诊断方法 | |
CN100478650C (zh) | 改善发电设备故障检测的方法和装置 | |
CN107367353A (zh) | 在线动平衡处理转子弯曲故障校正可行性判定方法和系统 | |
CN110646138A (zh) | 一种旋转机械无键相无试重动平衡方法和分析装置 | |
CN109085375A (zh) | 一种旋转机械的键相脉冲信号监测方法及装置 | |
CN104180823A (zh) | 一种温度补偿方法及装置 | |
CN106596034B (zh) | 一种判断模型试验台旋翼天平试验数据有效性的方法 | |
JPH0337133B2 (zh) | ||
WO2023216451A1 (zh) | 一种超重力离心机动平衡长期监测方法 | |
CN112945535B (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 | ||
C53 | Correction of patent for invention or patent application | ||
CB03 | Change of inventor or designer information |
Inventor after: Duan Lixiang Inventor after: Zhang Laibin Inventor before: Duan Lixiang |
|
COR | Change of bibliographic data |
Free format text: CORRECT: INVENTOR; FROM: DUAN LIXIANG TO: DUAN LIXIANG ZHANG LAIBIN |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |