CN114139428B - 一种转子现场动平衡多目标优化方法及系统 - Google Patents
一种转子现场动平衡多目标优化方法及系统 Download PDFInfo
- Publication number
- CN114139428B CN114139428B CN202111223314.4A CN202111223314A CN114139428B CN 114139428 B CN114139428 B CN 114139428B CN 202111223314 A CN202111223314 A CN 202111223314A CN 114139428 B CN114139428 B CN 114139428B
- Authority
- CN
- China
- Prior art keywords
- rotor
- signal
- vibration
- balance
- optimization
- 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
- 238000005457 optimization Methods 0.000 title claims abstract description 76
- 238000000034 method Methods 0.000 title claims abstract description 56
- 239000002245 particle Substances 0.000 claims abstract description 40
- 238000006073 displacement reaction Methods 0.000 claims abstract description 36
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 35
- 238000004364 calculation method Methods 0.000 claims abstract description 30
- 238000012360 testing method Methods 0.000 claims abstract description 18
- 230000006870 function Effects 0.000 claims description 30
- 230000002159 abnormal effect Effects 0.000 claims description 22
- 238000012545 processing Methods 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 16
- 238000001914 filtration Methods 0.000 claims description 14
- 238000012937 correction Methods 0.000 claims description 7
- 230000008602 contraction Effects 0.000 claims description 6
- 238000011156 evaluation Methods 0.000 claims description 6
- 230000000630 rising effect Effects 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 238000005520 cutting process Methods 0.000 claims description 2
- 230000010363 phase shift Effects 0.000 claims description 2
- 239000000758 substrate Substances 0.000 claims 3
- 230000007774 longterm Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 7
- 238000011065 in-situ storage Methods 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000005299 abrasion Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005272 metallurgy Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000013450 outlier detection Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Computational Linguistics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Hardware Design (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Geometry (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Testing Of Balance (AREA)
Abstract
本发明公开一种转子现场动平衡多目标优化方法及系统,方法包括获取转子原始振动信号、键相信号、加试重后振动信号和键相信号,利用互相关算法计算得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位;构建以各平衡转速下各测点残余振动位移幅值平方和最小、最大残余振动位移幅值最小为优化目标的优化函数;通过基于速度限制的多目标粒子群算法对所述目标函数进行转子动平衡优化计算;以各平衡转速下各测点的平衡效率和最大为决策目标,得到最优的配重方案。本发明可以得到最优的动平衡配重方案,保证转子长期平稳运行,并可用于各类旋转机械的刚性转子和柔性转子的现场动平衡。
Description
技术领域
本发明涉及旋转机械转子动平衡领域,尤其涉及一种转子现场动平衡多目标优化方法及系统。
背景技术
旋转机械如离心压缩机、汽轮发电机、风机、燃气轮机等,是石油、化工、冶金、电力等工业领域的关键设备,支撑着工业技术的发展。转子是旋转机械最重要和核心的部件,随着工业的发展,要求转子更加高速化、精密化,转子的转速也要求越来越高,有的转子甚至达到了十几万转每分钟。转子的工作转速大于一阶临界转速时,在运转过程中转子自挠曲变形不可忽略,挠曲变形会产生质量偏心,导致运行过程中产生与转频同频的离心力,进而产生转子自身的不平衡振动;转子在设计时的误差,制造过程中的材质不均匀、偏心,装配过程中的不对称、不同心以及运行过程中零件的飞出、磨损、腐蚀、积灰等都会产生不平衡量。不平衡是引起旋转机械振动的主要原因,会引起转子挠曲和内应力,加速轴承和轴封等零件的磨损,降低机器的工作效率和使用寿命,严重时会造成重大的经济损失和人员伤亡。
不平衡不仅是旋转机械主要的激振源,也是许多自激振动的诱发因素。统计资料表明:约有80%的振动可直接归因于转子不平衡。不平衡质量越大、机械运转速度越高,机械振动越剧烈。一旦平衡状况得以改善,其他一些故障现象也将随之消失,特别是在如今转子系统朝着超细长、超高速的发展趋势下,转子不平衡问题必须快速、高效的解决,才能保证旋转机械安全稳定运行。因此转子动平衡技术是现代工业中的一项关键技术。
转子动平衡方法可分为平衡机法和现场动平衡法,平衡机法需要将转子拆卸到动平衡机上进行动平衡,不仅操作繁琐,而且由于动平衡机和现场工况的不同会导致平衡效果难以满足实际工程应用要求。基于最小二乘影响系数法的现场动平衡法是转子动平衡领域应用最广泛的方法,但该方法会因影响系数矩阵的病态问题导致最终平衡的结果误差较大。基于单目标优化的转子动平衡方法虽然能够避免最小二乘影响系数法影响系数矩阵的病态问题,但优化目标为各平衡转速下各测点残余振动位移幅值平方和最小,难以保证残余振动的均匀性。
发明内容
针对现有转子现场动平衡方法中影响系数矩阵病态的缺陷以及残余振动均匀性难以保证的问题,本发明提出一种转子现场动平衡多目标优化方法,以实现现场动平衡时转子最优配重的快速,精确计算,提高动平衡效率。
为了实现上述目的,本发明提供了如下的技术方案。
一种转子现场动平衡多目标优化方法,包括以下步骤:
获取转子原始振动信号、键相信号、加试重后振动信号和键相信号,利用互相关算法计算得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位;
构建以各平衡转速下各测点残余振动位移幅值平方和最小、最大残余振动位移幅值最小为优化目标的优化函数;
通过基于速度限制的多目标粒子群算法对所述目标函数进行转子动平衡优化计算;
以各平衡转速下各测点的平衡效率均方根值最大为决策目标,选择最优的配重方案。
作为本发明的进一步改进,所述获取转子原始振动信号、键相信号、加试重后振动信号和键相信号,具体包括:
获取各平衡转速下的键相信号、各测点的转子原始振动位移信号以及在校正面上加重后的键相信号、各测点的转子振动位移信号;
对采集的转子键相信号进行中值滤波处理,去除噪声干扰,得到处理后的转子键相信号;
对采集的转子振动位移信号进行振动异常值检测和替换处理;
对振动异常值处理后的转子振动信号进行滤波处理,去除噪声干扰;
利用预处理后的转子键相信号和转子振动信号进行分析计算,得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位。
作为本发明的进一步改进,所述对采集的转子振动位移信号进行振动异常值检测和替换处理具体是采用基于移动窗口的中位数绝对偏差法进行转子振动信号的振动异常值检测,检测到振动异常值后利用分段三次样条插值对振动异常值进行替换。
作为本发明的进一步改进,所述对采集的转子振动位移信号进行振动异常值检测和替换处理具体包括:
确定移动窗口的大小,并以移动窗口的大小将采集的振动信号划分为若干子集;
计算每个子集内各个数据点之间的距离并进行标准化,得到各个子集内各个数据点之间的标准化距离
计算各个子集的稳健距离RKi,RKi为子集内各个数据点之间的标准化距离减去子集中位数的绝对值;
计算临界值C,C为子集的稳健距离RDi的中位数乘以正态分布的比例因子常量;
判断RDi是否大于3倍的C,若大于3倍的C说明第i个点为振动异常值;
利用分段三次样条插值对振动异常值进行插值替换。
作为本发明的进一步改进,所述对振动异常值处理后的转子振动信号进行滤波处理,去除噪声干扰,具体包括:
采用零相移滤波器对振动信号进行带通滤波,滤除转子转频以外的频率分量,得到纯净的转子振动信号。
作为本发明的进一步改进,所述利用互相关算法计算得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位,具体包括:
以中值滤波后转子键相信号的第一个上升沿和最后一个上升沿为基准将转子振动信号整周期截断,利用互相关算法计算得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位。
作为本发明的进一步改进,所述通过基于速度限制的多目标粒子群算法对所述目标函数进行转子动平衡优化计算,其中,基于速度限制的多目标粒子群算法引入速度收缩和多项式变异,构建外部精英档案来存储非支配解,并通过计算拥挤距离来确定粒子是否被放入外部精英档案;转子动平衡优化计算,具体包括:
a)初始化粒子种群和外部精英档案;
b)利用速度收缩计算每个粒子的速度;
c)根据粒子速度计算并更新粒子的位置;
d)对部分粒子进行多项式变异;
e)对粒子进行拥挤距离计算并评估,根据评估结果更新粒子种群和外部精英档案;
f)判断算法迭代数是否小于最大迭代数,若小于最大迭代数则继续进行步骤a)至步骤f),否则退出循环,得到的外部精英档案为转子动平衡优化计算的结果。
一种转子现场动平衡多目标优化系统,包括:
幅值及相位计算模块,用于获取转子原始振动信号、键相信号、加试重后振动信号和键相信号,利用互相关算法计算得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位;
优化函数构建模块,用于构建以各平衡转速下各测点残余振动位移幅值平方和最小、最大残余振动位移幅值最小为优化目标的优化函数;
平衡优化计算模块,用于通过基于速度限制的多目标粒子群算法对所述目标函数进行转子动平衡优化计算;
配重方案选择模块,用于以各平衡转速下各测点的平衡均方根最大为决策目标,选择最优的配重方案。
一种电子设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现所述转子现场动平衡多目标优化方法的步骤。
一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现所述转子现场动平衡多目标优化方法的步骤。
与现有技术相比,本发明具有以下有益效果:
本发明公开了一种转子现场动平衡多目标优化方法,通过对转子原始振动信号和键相信号以及加试重后振动信号和键相信号的采集、处理、分析、计算,得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位。构建以各平衡转速下各测点残余振动位移幅值平方和最小、最大残余振动位移幅值最小为优化目标的优化函数,并基于速度限制的多目标粒子群算法进行转子动平衡优化计算,以各平衡转速下各测点的平衡效率和最大为决策目标,能够快速、精确地得到最优的配重方案。根据最优配重方案对转子各校正面进行加重校正,达到转子稳定运行要求。本发明避免了传统的最小二乘影响系数法因影响系数矩阵病态导致动平衡精度变差以及单目标优化算法难以保证残余振动均匀性的问题,可以得到最优的动平衡配重方案,保证转子长期平稳运行,并可用于各类旋转机械的刚性转子和柔性转子的现场动平衡。
附图说明
在此描述的附图仅用于解释目的,而不意图以任何方式来限制本发明公开的范围。另外,图中的各部件的形状和比例尺寸等仅为示意性的,用于帮助对本发明的理解,并不是具体限定本发明各部件的形状和比例尺寸。在附图中:
图1是本发明涉及的一种转子现场动平衡多目标优化方法总体流程图;
图2是本发明涉及的一种转子现场动平衡多目标优化方法的振动异常值检测、替换流程图;
图3是本发明涉及的一种转子现场动平衡多目标优化方法的不平衡信号相位、幅值互相关提取流程图;
图4是本发明涉及的一种转子现场动平衡多目标优化方法的多目标优化计算流程图。
图5为本发明优选实施例转子现场动平衡多目标优化系统结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明中的技术方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
需要说明的是,当元件被称为“设置于”另一个元件,它可以直接在另一个元件上或者也可以存在居中的元件。当一个元件被认为是“连接”另一个元件,它可以是直接连接到另一个元件或者可能同时存在居中元件。本文所使用的术语“垂直的”、“水平的”、“左”、“右”以及类似的表述只是为了说明的目的,并不表示是唯一的实施例。
除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的,不是旨在于限制本发明。本文所使用的术语“和/或”包括一个或多个相关的所列项目的任意的和所有的组合。
本发明第一个目的是提供一种转子现场动平衡多目标优化方法,包括以下步骤:
获取转子原始振动信号、键相信号、加试重后振动信号和键相信号,利用互相关算法计算得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位;
构建以各平衡转速下各测点残余振动位移幅值平方和最小、最大残余振动位移幅值最小为优化目标的优化函数;
通过基于速度限制的多目标粒子群算法对所述目标函数进行转子动平衡优化计算;
以各平衡转速下各测点的平衡效率均方根值最大为决策目标,选择最优的配重方案。
通过采用上述方法,无需将转子拆卸至动平衡机进行动平衡,能够快速准确的实现旋转机械转子现场动平衡的优化计算,得到最优的配重方案,保证在该配重方案下各测振面的振动均得到很好的消除,提高动平衡效率,保证转子长期平稳运行。
以下结合具体附图和实施例对本发明的方法进行详细说明。
实施例
一种转子现场动平衡多目标优化方法,如图1所示,步骤如下:
S1:通过数据采集设备得到各平衡转速下的键相信号、各测点的转子原始振动位移信号以及在校正面上加重后的键相信号、各测点的转子振动位移信号;
S2:对采集的转子键相信号进行中值滤波处理,去除噪声干扰,得到处理后的转子键相信号;
S3:对采集的转子振动位移信号进行振动异常值检测和替换处理;
S4:对振动异常值处理后的转子振动信号进行滤波处理,去除噪声干扰;
S5:利用预信号处理后的转子键相信号和转子振动信号进行分析计算,得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位;
S6:构建优化目标函数;
S7:使用速度限制的多目标粒子群算法对步骤S6中得到的目标函数进行优化计算;
S8:构建决策目标,从步骤S7计算得到的解集中选择最优的配重方案;
具体地,如图2所示,所述步骤S3中,采用基于移动窗口的中位数绝对偏差法进行转子振动信号的振动异常值检测,检测到离群点后利用分段三次样条插值对离群点进行替换。步骤S3具体为:
S3.1确定移动窗口的大小,并以移动窗口的大小将采集的振动信号划分为若干子集;
S3.2计算每个子集内各个数据点之间的距离并进行标准化,得到各个子集内各个数据点之间的标准化距离
S3.3计算各个子集的稳健距离RDi,RDi为子集内各个数据点之间的标准化距离减去子集中位数的绝对值;
S3.4计算临界值C,C为子集的稳健距离RKi的中位数乘以正态分布的比例因子常量1.4826;
S3.5判断RKi是否大于3倍的C,若大于3倍的C说明第i个点为离群点;
S3.6利用分段三次样条插值对离群点进行插值替换;
进一步地,步骤S4中,采用零相位滤波器对振动信号进行带通滤波,滤除转子转频以外的频率分量,得到纯净的转子振动信号。
具体地,如图3所示,所述步骤S5中,以中值滤波后转子键相信号的第一个上升沿和最后一个上升沿为基准将转子振动信号整周期截断,利用互相关算法计算得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位。
进一步地,所述步骤S6中,以各平衡转速下各测点残余振动位移幅值平方和最小、最大残余振动位移幅值最小为目标构建优化函数。
具体地,如图4所示,所述步骤S7中,使用速度限制的多目标粒子群算法对步骤S6中的目标函数进行优化计算。该算法引入速度收缩和多项式变异保证解的快速收敛和多样性,构建一个外部精英档案来存储非支配解,并通过计算拥挤距离来确定粒子是否被放入外部精英档案。步骤S7具体为:
S7.1初始化粒子种群(位置、速度、最佳个体)和外部精英档案;
S7.2利用速度收缩计算每个粒子的速度;
S7.3根据粒子速度计算并更新粒子的位置;
S7.4对部分粒子进行多项式变异;
S7.5对粒子进行拥挤距离计算并评估,根据评估结果更新粒子种群和外部精英档案;
S7.6判断算法迭代数是否小于最大迭代数,若小于最大迭代数则继续进行步骤7.2至步骤7.5,否则退出循环,得到的外部精英档案为转子动平衡优化计算的结果;
进一步地,步骤S8中,以各平衡转速下各测点的平衡效率和最大为决策目标,从步骤S7中的外部精英档案中选择最优的配重方案。
根据最优配重方案对转子各校正面进行加重校正,达到转子稳定运行要求。
如图4所示,本发明的另一目的在于提出一种转子现场动平衡多目标优化系统,包括:
幅值及相位计算模块,用于获取转子原始振动信号、键相信号、加试重后振动信号和键相信号,利用互相关算法计算得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位;
优化函数构建模块,用于构建以各平衡转速下各测点残余振动位移幅值平方和最小、最大残余振动位移幅值最小为优化目标的优化函数;
平衡优化计算模块,用于通过基于速度限制的多目标粒子群算法对所述目标函数进行转子动平衡优化计算;
配重方案选择模块,用于以各平衡转速下各测点的平衡均方根最大为决策目标,选择最优的配重方案。
本发明第三个目的是提供一种电子设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现所述转子现场动平衡多目标优化方法的步骤。
所述转子现场动平衡多目标优化方法包括以下步骤:
获取转子原始振动信号、键相信号、加试重后振动信号和键相信号,利用互相关算法计算得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位;
构建以各平衡转速下各测点残余振动位移幅值平方和最小、最大残余振动位移幅值最小为优化目标的优化函数;
通过基于速度限制的多目标粒子群算法对所述目标函数进行转子动平衡优化计算;
以各平衡转速下各测点的平衡效率均方根值最大为决策目标,选择最优的配重方案。
本发明第四个目的是提供一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现所述转子现场动平衡多目标优化方法的步骤。
所述转子现场动平衡多目标优化方法包括以下步骤:
获取转子原始振动信号、键相信号、加试重后振动信号和键相信号,利用互相关算法计算得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位;
构建以各平衡转速下各测点残余振动位移幅值平方和最小、最大残余振动位移幅值最小为优化目标的优化函数;
通过基于速度限制的多目标粒子群算法对所述目标函数进行转子动平衡优化计算;
以各平衡转速下各测点的平衡效率均方根值最大为决策目标,选择最优的配重方案。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求保护范围之内。
Claims (7)
1.一种转子现场动平衡多目标优化方法,其特征在于,包括以下步骤:
获取转子原始振动信号、键相信号、加试重后振动信号和键相信号,利用互相关算法计算得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位;
构建以各平衡转速下各测点残余振动位移幅值平方和最小、最大残余振动位移幅值最小为优化目标的优化函数;
通过基于速度限制的多目标粒子群算法对所述目标函数进行转子动平衡优化计算;
以各平衡转速下各测点的平衡效率均方根值最大为决策目标,选择最优的配重方案;
所述获取转子原始振动信号、键相信号、加试重后振动信号和键相信号,具体包括:
获取各平衡转速下的键相信号、各测点的转子原始振动位移信号以及在校正面上加重后的键相信号、各测点的转子振动位移信号;
对采集的转子键相信号进行中值滤波处理,去除噪声干扰,得到处理后的转子键相信号;
对采集的转子振动位移信号进行振动异常值检测和替换处理;
对振动异常值处理后的转子振动信号进行滤波处理,去除噪声干扰;
利用预处理后的转子键相信号和转子振动信号进行分析计算,得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位;
所述利用互相关算法计算得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位,具体包括:
以中值滤波后转子键相信号的第一个上升沿和最后一个上升沿为基准将转子振动信号整周期截断,利用互相关算法计算得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位;
所述通过基于速度限制的多目标粒子群算法对所述目标函数进行转子动平衡优化计算,其中,基于速度限制的多目标粒子群算法引入速度收缩和多项式变异,构建外部精英档案来存储非支配解,并通过计算拥挤距离来确定粒子是否被放入外部精英档案;转子动平衡优化计算,具体包括:
a)初始化粒子种群和外部精英档案;
b)利用速度收缩计算每个粒子的速度;
c)根据粒子速度计算并更新粒子的位置;
d)对部分粒子进行多项式变异;
e)对粒子进行拥挤距离计算并评估,根据评估结果更新粒子种群和外部精英档案;
f)判断算法迭代数是否小于最大迭代数,若小于最大迭代数则继续进行步骤a)至步骤f),否则退出循环,得到的外部精英档案为转子动平衡优化计算的结果。
2.根据权利要求1所述的方法,其特征在于,
所述对采集的转子振动位移信号进行振动异常值检测和替换处理具体是采用基于移动窗口的中位数绝对偏差法进行转子振动信号的振动异常值检测,检测到振动异常值后利用分段三次样条插值对振动异常值进行替换。
3.根据权利要求1所述的方法,其特征在于,
所述对采集的转子振动位移信号进行振动异常值检测和替换处理具体包括:
确定移动窗口的大小,并以移动窗口的大小将采集的振动信号划分为若干子集;
计算每个子集内各个数据点之间的距离并进行标准化,得到各个子集内各个数据点之间的标准化距离
计算各个子集的稳健距离RDi,RDi为子集内各个数据点之间的标准化距离减去子集中位数的绝对值;
计算临界值C,C为子集的稳健距离RDi的中位数乘以正态分布的比例因子常量;
判断RDi是否大于3倍的C,若大于3倍的C说明第i个点为振动异常值;
利用分段三次样条插值对振动异常值进行插值替换。
4.根据权利要求1所述的方法,其特征在于,
所述对振动异常值处理后的转子振动信号进行滤波处理,去除噪声干扰,具体包括:
采用零相移滤波器对振动信号进行带通滤波,滤除转子转频以外的频率分量,得到纯净的转子振动信号。
5.一种转子现场动平衡多目标优化系统,基于权利要求1至4任一项所述的方法其特征在于,包括:
幅值及相位计算模块,用于获取转子原始振动信号、键相信号、加试重后振动信号和键相信号,利用互相关算法计算得到转子原始不平衡信号和加试重后不平衡信号的幅值及相位;
优化函数构建模块,用于构建以各平衡转速下各测点残余振动位移幅值平方和最小、最大残余振动位移幅值最小为优化目标的优化函数;
平衡优化计算模块,用于通过基于速度限制的多目标粒子群算法对所述目标函数进行转子动平衡优化计算;
配重方案选择模块,用于以各平衡转速下各测点的平衡均方根最大为决策目标,选择最优的配重方案。
6.一种电子设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现权利要求1-4任一项所述转子现场动平衡多目标优化方法的步骤。
7.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现权利要求1-4任一项所述转子现场动平衡多目标优化方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111223314.4A CN114139428B (zh) | 2021-10-20 | 2021-10-20 | 一种转子现场动平衡多目标优化方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111223314.4A CN114139428B (zh) | 2021-10-20 | 2021-10-20 | 一种转子现场动平衡多目标优化方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114139428A CN114139428A (zh) | 2022-03-04 |
CN114139428B true CN114139428B (zh) | 2024-05-10 |
Family
ID=80395198
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111223314.4A Active CN114139428B (zh) | 2021-10-20 | 2021-10-20 | 一种转子现场动平衡多目标优化方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114139428B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116317781B (zh) * | 2023-04-17 | 2024-01-30 | 北京化工大学 | 基于微速差的同轴双转子系统不平衡振动控制方法和装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5818135A (ja) * | 1981-07-27 | 1983-02-02 | Hitachi Ltd | 回転機のバランス方法及び該バランス方法に使用する装置 |
CN101556200A (zh) * | 2008-06-24 | 2009-10-14 | 郑州大学 | 基于全矢谱的柔性转子动平衡方法 |
CN109632188A (zh) * | 2019-01-29 | 2019-04-16 | 北京化工大学 | 虑及多转速状态的柔性转子动平衡双重优化方法及系统 |
-
2021
- 2021-10-20 CN CN202111223314.4A patent/CN114139428B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5818135A (ja) * | 1981-07-27 | 1983-02-02 | Hitachi Ltd | 回転機のバランス方法及び該バランス方法に使用する装置 |
CN101556200A (zh) * | 2008-06-24 | 2009-10-14 | 郑州大学 | 基于全矢谱的柔性转子动平衡方法 |
CN109632188A (zh) * | 2019-01-29 | 2019-04-16 | 北京化工大学 | 虑及多转速状态的柔性转子动平衡双重优化方法及系统 |
Non-Patent Citations (2)
Title |
---|
张西宁 ; 赵明 ; 温广瑞 ; .一种轴系全息动平衡配重的多目标优化方法.中国机械工程.2010,(16),全文. * |
运侠伦 ; 梅雪松 ; 姜歌东 ; 胡振邦 ; 张尊浩 ; .辨识主轴转子系统能量传递特性的矢量匹配平衡方法研究.振动与冲击.2019,(10),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN114139428A (zh) | 2022-03-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6297555B2 (ja) | 回転ブレードの健康状態を監視するための方法及びシステム | |
EP2746541B1 (en) | A method of detecting and localizing partial rotor-stator rubbing during the operation of a turbine | |
CA2054598C (en) | Apparatus and method for removing common mode vibration data from digital turbine blade vibration data | |
CN103608739A (zh) | 振动监测 | |
CN114139428B (zh) | 一种转子现场动平衡多目标优化方法及系统 | |
CN101622457B (zh) | 用于确定涡轮气泵的泵转子的疲劳的方法 | |
JP2009506261A (ja) | タービンの回転軸の動的挙動を監視するための方法と装置 | |
CN104215395B (zh) | 一种转子不平衡故障的检测方法及装置 | |
EP3241987B1 (en) | Rotor balancing | |
JP4772594B2 (ja) | 回転機器の低速バランス法及び低速バランス実施装置 | |
Pavlenko et al. | Using Computer Modeling and Artificial Neural Networks for Ensuring the Vibration Reliability of Rotors. | |
CN110646138A (zh) | 一种旋转机械无键相无试重动平衡方法和分析装置 | |
CN112329357A (zh) | 一种清水离心泵振动故障的简易诊断方法及诊断系统 | |
Zhang et al. | Vibration control of tie rod rotors with optimization of unbalanced force and unbalanced moment | |
CN109085375B (zh) | 一种旋转机械的键相脉冲信号监测方法及装置 | |
CN112491228B (zh) | 基于振动波形检测引起定子低频振动的关键磁极的方法 | |
US20150066389A1 (en) | Systems and methods to monitor a rotating component | |
CN111898238B (zh) | 一种有约束的高速动平衡力学解算方法 | |
RU2372595C1 (ru) | Способ балансировки сборного ротора | |
RU2449180C1 (ru) | Способ балансировки ротора | |
RU2426014C1 (ru) | Расчетно-имитационный способ балансировки вала | |
CN116209827A (zh) | 风力涡轮机的风力转子的不平衡估计 | |
Mashiachidi et al. | Prediction of fatigue life of mistuned steam turbine blades subjected to variations in blade geometry | |
CN113312730A (zh) | 一种双驱汽轮机转子应力监测方法 | |
CN111695206A (zh) | 一种旋转态叶盘的失谐识别和模型更新方法 |
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 |