一种大型旋转机械的转子不平衡故障的实时诊断方法技术

本发明专利技术涉及一种大型旋转机械的转子不平衡故障的实时诊断方法，属于旋转机械振动监测技术领域。该方法首先通过监测的振动数据拟合转子系统的频响函数并辨识转子系统的相关物理参数；其次,根据该类转子系统规定的平衡允差及其物理参数，确定其产生不平衡故障的临界激励；最后,基于转子系统的频响函数和产生不平衡故障的临界激励，计算转子系统产生不平衡故障的临界振幅,该临界振幅可作为转子系统是否产生转子不平衡故障的依据。本发明专利技术反映了系统故障状态的振动信号与系统物理参数之间的联系，并以临界振幅作为产生不平衡故障的依据，可对大型旋转机械转子系统的原始不平衡、渐发性不平衡和突发性不平衡进行实时监测、诊断。断。断。

【技术实现步骤摘要】
一种大型旋转机械的转子不平衡故障的实时诊断方法


[0001]本专利技术属于旋转机械振动监测
，涉及一种大型旋转机械的转子不平衡故障的实时诊断方法。

技术介绍

[0002]旋转机械是一类被广泛使用的机械设备，在许多行业中，旋转机械造成的恶性事故会造成极大的损失，因此对旋转机械的运行实施有效地检测和诊断是非常必要的。常见的汽轮发电机组、风机等典型旋转机械常被作为特护机组进行监控和维护，其最常见的故障是转子不平衡，在现场发生的机组振动故障中，约80％是由于转子质量不平衡造成的。引起质量不平衡的原因是多种多样的，如旋转设备转子结构设计缺陷、安装制造精度低、旋转部件运行中，存在介质持续点蚀、磨损、介质结垢等，都可能成为转子质量不平衡的原因。
[0003]目前，大型旋转机械的振动分析与故障诊断工作主要还是依靠具有一定现场振动故障诊断经验的专家完成，诊断经济成本高，周期长，且对转子运行中逐步发展的渐发性不平衡和突发性不平衡，诊断难度大。因此，亟需一种大型旋转机械的转子不平衡故障实时的智能诊断方法来降低诊断成本，且保证诊断准确度和效率。

技术实现思路

[0004]有鉴于此，本专利技术的目的在于提供一种大型旋转机械的转子不平衡故障的实时诊断方法，对转子不平衡故障进行实时自动监测、分析、诊断，以提高故障诊断效率和准确度。
[0005]为达到上述目的，本专利技术提供如下技术方案：
[0006]一种大型旋转机械的转子不平衡故障的实时诊断方法，具体包括以下步骤：
[0007]S1：通过监测的振动数据拟合转子系统的频响函数，并辨识转子系统的物理参数，包括质量、阻尼和刚度特性等；
[0008]S2：根据相关(国家标准GB/T 9239.1
‑
2006《恒态(刚性)转子平衡品质要求》)规定的平衡允差及其物理参数，确定转子系统产生不平衡故障的临界激励；
[0009]S3：根据上述确定的转子系统的频响函数和产生不平衡故障的临界激励，计算转子系统产生不平衡故障的临界振幅，根据临界振幅判定转子系统是否产生不平衡故障。
[0010]进一步，步骤S1具体包括以下步骤：
[0011]S11：基于实测的振动数据拟合频响函数幂基多项式，并采用正交多项式拟合法计算频响函数幂基多项式的系数；
[0012]S12：根据动力学方程构建转子系统的加速度频响函数，由于该加速度频响函数与S11拟合的频响函数幂基多项式相等，则可计算转子系统的物理参数。
[0013]进一步，步骤S11具体包括：对于一个N自由度的转子系统，其加速度频响函数矩阵的第l行第p列元素表示为：
[0014][0015]其中，n＝2N，m＝2N
‑
2，N(s)、D(s)分别为传递函数的分子、分母多项式，a0,a1,
…
,a
m
、b0,b1,
…
,b
n
表示分子、分母多项式系数，且令s＝jω，ω表示频率；p0(s)＝1，p1(s)＝s，
…
，p
m
(s)＝s
m
；q0(s)＝1，q1(s)＝s，
…
，q
n
(s)＝s
n
。
[0016]采用正交多项式拟合法计算加速度频响函数幂基多项式系数；
[0017]首先，构造如下矩阵：
[0018][0019][0020]W＝[q
n
(jω
‑
L
),
…
,q
n
(jω
‑1),q
n
(jω1),
…
,q
n
(jω
L
)]T
，A＝[a0,a1,
…
,a
m
]T
，B＝[b0,b1,
…
,b
n
]T
；其中，A和B分别为分子和分母幂基多项式的系数向量，E为加速度频响的附加项矩阵，D为实测加速度频响序列的扩展矩阵，P、Q和W分别为分子、分母的扩展幂基矩阵和向量，L为被拟合频响数据的测量频率点数。本专利技术引入负频率的概念，使ω＝ω
‑
L
,
…
,ω
‑1,ω1,
…
,ω
L
共2L个点，且令ω
‑
i
＝
‑
ω
i
，频响函数的数学模型值H(jω
‑
i
)＝H
*
(jω
i
)。
[0021]其次，计算频响函数的幂基多项式系数，基于频响函数的实测值与理论值之间的总方差最小的理论，采用下式求得频响函数的幂基多项式系数，即矩阵A和B；
[0022][0023]其中，P
a
＝EP，Q
a
＝DQ，W
a
＝DW。
[0024]进一步，根据实测加速度振动数据对矩阵D进行构造，具体包括以下步骤：
[0025](1)计算转子径向动载荷：采用斯捷潘若夫公式计算各测量转速下转子的径向动载荷，该载荷即为激励力F；
[0026][0027]其中，P1为风机的全压，Pa；D2为叶轮出口直径，cm；B2为包括前后盖板的叶轮宽度，cm；q
v
为风机的工况流量，m3/s；q
vd
为风机设计工况流量，m3/s；
[0028](2)加速度频响序列构造：在相应转速下，通过计算实测振动加速度与激励力的比值，构造加速度频响序列，即矩阵D。
[0029]进一步，步骤S12中，根据动力学方程构建转子系统的加速度频响函数，具体包括：
[0030]下面以二自由度转子系统为例对本案进行说明，其方法可扩展到多自由度，二自由度转子系统的等效模型如附图2所示。
[0031]根据动力学方程，二自由度转子系统的加速度频响函数矩阵为：
[0032]H(ω)＝(jω)2[K+jωC+(jω)2M]‑1[0033]其中，表示质量矩阵，表示阻尼矩阵，表示刚度矩阵；m1、c1、k1为转子的等效质量、阻尼和刚度，m2、c2、k2为非旋转部件的质量、阻尼和刚度。
[0034]考虑实际测量的方便，常在轴承座的径向方向安装振动传感器测量其振动，即需要通过计算频响函数矩阵的第1行第2列元素H
12
(s)辨识转子系统的相关物理参数。根据动力学方程，该元素表示为：
[0035][0036]根据步骤S11，二自由度转子系统的加速度频响函数矩阵中第1行第2列元素为：
[0037][0038]由于实测振动数据拟合的频响函数幂基多项式与动力学方程计算的频响函数相等，则根据上述两式即可计算转子系统的相关物理参数。
[0039]进一步，步骤S2具体包括以下步骤：
[0040]S21：根据国家标准确定转子本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种大型旋转机械的转子不平衡故障的实时诊断方法，其特征在于，该方法具体包括以下步骤：S1：通过监测的振动数据拟合转子系统的频响函数，并辨识转子系统的物理参数；S2：根据转子系统规定的平衡允差及其物理参数，确定转子系统产生不平衡故障的临界激励；S3：根据转子系统的频响函数和产生不平衡故障的临界激励，计算转子系统产生不平衡故障的临界振幅，从而依据临界振幅判定转子系统是否产生不平衡故障。2.根据权利要求1所述的大型旋转机械的转子不平衡故障的实时诊断方法，其特征在于，步骤S1具体包括以下步骤：S11：基于实测的振动数据拟合频响函数幂基多项式，并采用正交多项式拟合法计算频响函数的幂基多项式系数；S12：根据动力学方程构建转子系统的加速度频响函数，由于该加速度频响函数与S11拟合的频响函数幂基多项式相等，即能计算转子系统的物理参数。3.根据权利要求2所述的大型旋转机械的转子不平衡故障的实时诊断方法，其特征在于，步骤S11具体包括：对于一个N自由度的转子系统，其加速度频响函数矩阵的第l行第p列元素表示为：其中，n＝2N，m＝2N
‑
2，N(s)、D(s)分别为传递函数的分子、分母多项式，a0,a1,
…
,a
m
、b0,b1,
…
,b
n
分别表示分子、分母多项式系数，且令s＝jω，ω表示频率；p0(s)＝1，p1(s)＝s，
…
，p
m
(s)＝s
m
；q0(s)＝1，q1(s)＝s，
…
，q
n
(s)＝s
n
；采用正交多项式拟合法计算加速度频响函数幂基多项式系数；首先，构造如下矩阵：首先，构造如下矩阵：W＝[q
n
(jω
L
),
…
,q
n
(jω1),q
n
(jω1),
…
,q
n
(jω
L
)]
T
，A＝[a0,a1,
…
,a
m
]
T
，B＝[b0,
b1,
…
,b
n
]
T
；其中，A和B分别为分子、分母幂基多项式的系数向量，E为加速度频响的附加项矩阵，D为实测加速度频响序列的扩展矩阵，P、Q和W分别为分子、分母的扩展幂基矩阵和向量，L为被拟合频响数据的测量频率点数；引入负频率的概...

【专利技术属性】
技术研发人员：张霞，陈惠，张鹏，石明全，王晓东，
申请(专利权)人：中国科学院重庆绿色智能技术研究院，
类型：发明
国别省市：