CN110793615B - 风机振动状态检测方法、装置及配合端等效激励确定方法 - Google Patents
风机振动状态检测方法、装置及配合端等效激励确定方法 Download PDFInfo
- Publication number
- CN110793615B CN110793615B CN201911067941.6A CN201911067941A CN110793615B CN 110793615 B CN110793615 B CN 110793615B CN 201911067941 A CN201911067941 A CN 201911067941A CN 110793615 B CN110793615 B CN 110793615B
- Authority
- CN
- China
- Prior art keywords
- vibration
- signal
- fan
- amplitude
- monitoring part
- 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
- 230000005284 excitation Effects 0.000 title claims abstract description 131
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000001514 detection method Methods 0.000 title claims description 33
- 238000012544 monitoring process Methods 0.000 claims abstract description 168
- 238000005070 sampling Methods 0.000 claims abstract description 40
- 230000013011 mating Effects 0.000 claims description 72
- 230000000737 periodic effect Effects 0.000 claims description 25
- 238000004088 simulation Methods 0.000 claims description 24
- 230000036541 health Effects 0.000 claims description 17
- 238000012545 processing Methods 0.000 claims description 9
- 238000012952 Resampling Methods 0.000 claims description 5
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 12
- 238000013016 damping Methods 0.000 description 15
- 238000004458 analytical method Methods 0.000 description 14
- 230000008569 process Effects 0.000 description 10
- 230000009286 beneficial effect Effects 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000010248 power generation Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004806 packaging method and process Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 101100394003 Butyrivibrio fibrisolvens end1 gene Proteins 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
- G01H1/12—Measuring characteristics of vibrations in solids by using direct conduction to the detector of longitudinal or not specified vibrations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
- G01H1/12—Measuring characteristics of vibrations in solids by using direct conduction to the detector of longitudinal or not specified vibrations
- G01H1/16—Amplitude
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
- G01M7/02—Vibration-testing by means of a shake table
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种风机振动状态检测方法、装置及配合端等效激励确定方法,包括:按照设定采样周期对n个振动传感器在风机振动监测部位采集到的振动信号进行采样;根据采样得到的风机振动监测部位的振动信号,计算振动成分信号;根据振动成分信号的幅值以及振动传感器采集到振动信号的时间,计算风机振动监测部位配合端的等效激励信号;根据风机振动监测部位配合端的等效激励信号的周期性,确定风机振动监测部位的故障区域。本发明通过计算配合端等效激励信号来确定风机监测部件的故障情况,计算方法简单易行,能够准确地确定故障区域。
Description
技术领域
本发明涉及一种风机振动状态检测方法、装置及配合端等效激励确定方法,属于风力发电技术领域。
背景技术
作为一种清洁、可再生的新型能源,风电发展十分迅速,特别是海上风电开发已经成为全球的热点与前沿。风力发电机(简称风机)作为风力发电的关键设备,由于运行环境复杂,不可避免地会出现故障,如果因此导致风机的非正常停机,会严重影响发电效率。所以,监测风机的关键运行信号,分析诊断风机故障并确定故障区域是非常必要和迫切的。其中,风机的故障主要包括电气故障和机械故障,目前关于风机电气信号的监测和分析诊断系统已较为成熟,而对于风机机械信号的监测分析及故障诊断技术还相对匮乏。风机的机械故障中以振动引起的故障最为常见,因此监测并分析风机关键部件的振动影响的机械健康状态是十分必要和迫切的。
在现有技术中,对于风机部件振动影响的机械健康状态的监测和分析是通过在振动监测部位设置振动传感器,并将采集的振动数据通过电缆或光纤传输给主控PLC以及与PLC连接的服务器。现有技术中对风机部位进行振动监测和分析的对象是振动传感器采集的振动信号,而对于振动对监测部件的整体影响未涉及。实际上,振动对于监测部件的影响因为不同部位的承受能力不同而存在差异,因此,仅根据采集到的部分振动信号,会导致风机监测部分最终确定的故障区域不准确。由于对于风机部件健康状态评估是基于传感器采集的振动信号,因此会造成评估的滞后性,通常是在故障发生时或发生后,严重影响风机发电的效率。另外,振动传感器采集的振动信号会附有干扰信号,导致基于采集信号的分析存在较大误差,也会影响到故障区域的判断。此外,现有技术中对于风机监测部位健康状态的展示多属于数据和信号型的,不直观。并且现有对于振动造成影响的分析大多采用振动定位的方法,而在风机这种部件复杂配合的装置中明显不适用,并且振动定位方法通常计算过程复杂。
发明内容
本发明的目的是提供一种风机振动状态检测方法、装置及配合端等效激励确定方法,用于解决现有技术确定风机检测部件故障区域不准确的问题。
为解决上述技术问题,本发明提供了一种风机振动状态检测方法,步骤如下:
按照设定采样周期对n个振动传感器在风机振动监测部位采集到的振动信号进行再采样;
根据得到的风机振动监测部位的振动信号,计算振动成分信号;
根据振动成分信号的幅值以及振动传感器采集到振动信号的时间,计算风机振动监测部位配合端的等效激励信号;
根据风机振动监测部位配合端的等效激励信号的周期性,确定风机振动监测部位的故障区域。
本发明的有益效果是:通过将风机振动监测部位的振动信号均视为是通过在监测部位的配合端施加激励引起的,这样根据采集到的振动信号计算振动监测部位配合端的等效激励信号,根据配合端的等效激励信号的周期性,就可以获知振动监测部位发生故障的位置,以确定故障范围。本发明通过计算配合端等效激励信号来确定风机监测部件的故障情况,计算方法简单易行,能够准确地确定故障区域,在多部件复杂配合的风机装置振动影响分析中尤为有效。
进一步的,为了可靠获取等效激励信号以准确得到风机振动监测部位的健康状态,计算风机振动监测部位配合端的等效激励信号的步骤为:
若n个振动传感器在每个设定采样周期内均采集到1个振动信号X1(t)、X2(t)……Xn-1(t)、Xn(t),并且满足V1>V2>……>Vn-1>Vn且t1<t2<……<tn-1<tn或满足V1<V2<……<Vn-1<Vn且t1>t2>……>tn-1>tn,其中,Vi为第i个振动传感器采集到的振动信号计算出的振动成分信号的幅值,ti为第i个振动传感器采集到振动信号的时间,则根据其中两个振动成分信号的幅值以及对应两个振动传感器之间的距离,计算振动信号单位距离的衰减率;
根据选择的其中一个振动成分信号的幅值、振动信号单位距离的衰减率以及该选择的其中一个振动成分信号的幅值对应的振动传感器与目标配合端之间的距离,计算该目标配合端的等效激励信号,目标配合端为距离最大的一个振动成分信号的幅值对应的振动传感器最近的配合端。
进一步的,为了可靠获取等效激励信号以准确得到风机振动监测部位的健康状态,计算风机振动监测部位配合端的等效激励信号的步骤为:
若n个振动传感器在每个设定采样周期内均采集到2个振动信号X1,1(t)、X2,1(t)……Xn-1,1(t)、Xn,1(n)和X1,2(t)、X2,2(t)……Xn-1,2(t)、Xn,2(t),并且满足V1,1>V2,1>……>Vn-1,1>Vn,1,t1,1<t2,1<……<tn-1,1<tn,1且V1,2<V2,2<……<Vn-1,2<Vn,2,t1,2>t2,2>……>tn-1,2>tn,2或满足V1,1<V2,1<……<Vn-1,1<Vn,1,t1,1>t2,1>……>tn-1,1>tn,1且V1,2>V2,2>……>Vn-1,2>Vn,2,t1,2<t2,2<……<tn-1,2<tn,2,其中,Vi,1为第i个振动传感器采集到的第一振动信号计算出的振动成分信号的幅值,ti,1为第i个振动传感器采集到第一振动信号的时间,Vi,2为第i个振动传感器采集到的第二振动信号计算出的振动成分信号的幅值,ti,2为第i个振动传感器采集到第二振动信号的时间,则根据各振动成分信号的幅值V1,1、V2,1……Vn-1,1、Vn,1中的两个振动成分信号的幅值以及对应两个振动传感器之间的距离,计算第一振动信号单位距离的衰减率,根据各振动成分信号的幅值V1,2、V2,2……Vn-1,2、Vn,2中的两个振动成分信号的幅值以及对应两个振动传感器之间的距离,计算第二振动信号单位距离的衰减率;
根据各振动成分信号的幅值V1,1、V2,1……Vn-1,1、Vn,1中选择的一个振动成分信号的幅值、第一振动信号单位距离的衰减率以及该各振动成分信号的幅值V1,1、V2,1……Vn-1,1、Vn,1中选择的一个振动成分信号的幅值对应的振动传感器与第一目标配合端之间的距离,计算该第一目标配合端的等效激励信号,第一目标配合端为距离各振动成分信号的幅值V1,1、V2,1……Vn-1,1、Vn,1中最大的一个振动成分信号的幅值对应的振动传感器最近的配合端;根据各振动成分信号的幅值V1,2、V2,2……Vn-1,2、Vn,2中选择的一个振动成分信号的幅值、第二振动信号单位距离的衰减率以及该各振动成分信号的幅值V1,2、V2,2……Vn-1,2、Vn,2中选择的一个振动成分信号的幅值对应的振动传感器与第二目标配合端之间的距离,计算该第二目标配合端的等效激励信号,第二目标配合端为距离各振动成分信号的幅值V1,2、V2,2……Vn-1,2、Vn,2中最大的一个振动成分信号的幅值对应的振动传感器最近的配合端。
进一步的,为了可靠判断出影响振动监测部位的健康状态的故障区域,若振动监测部位的两个配合端的等效激励信号的幅值均大于设定阈值,且两个配合端的等效激励信号均为周期性信号,则判定该振动监测部位所属部件存在故障。
进一步的,为了可靠判断出影响振动监测部位的健康状态的故障区域,若振动监测部位的至少一个配合端的等效激励信号的幅值大于设定阈值,且在大于设定阈值的等效激励信号中只有第一配合端的等效激励信号为周期性信号,则判断是否满足与该第一配合端直接相连的部件的两个等效激励信号是周期性信号;若满足,则与该第一配合端直接相连的部件存在故障,否则,该第一配合端和与其直接相连的部件的配合区域存在故障。
进一步的,为了便于直观地看出振动监测部位的振动影响情况,还包括:将振动监测部位的配合端的等效激励信号施加于有限元仿真模型的对应位置,仿真得到振动监测部位的应变图像;若振动监测部位的配合端的等效激励信号的幅值大于设定阈值,且大于设定阈值的等效激励信号均不为周期性信号,则根据仿真得到振动监测部位各有限元点仿真应变与振动监测部位各点可承受额定应变的关系确定振动监测部位的健康状态。
进一步的,为了保证同一个振动信号在采样周期内只能在振动监测部位的传感器间传播一次,振动信号在振动监测部位相距最远的两个振动传感器间的传播时间<所述设定采样周期<振动信号在振动监测部位相距最远的两个振动传感器间的传播时间+振动信号在风机振动监测部位的不同配合端和距离该配合端最近的振动传感器之间的传播时间中的最小值。
为解决上述技术问题,本发明还提供了一种风机振动状态检测装置,包括处理器和存储器,所述处理器用于处理存储在存储器中的指令,以实现上述的风机振动状态检测方法。
本发明的有益效果是:通过将风机振动监测部位的振动信号均视为是通过在监测部位的配合端施加激励引起的,这样根据采集到的振动信号计算振动监测部位配合端的等效激励信号,根据配合端的等效激励信号的周期性,就可以获知振动监测部位发生故障的位置,以确定故障范围。本发明通过计算配合端等效激励信号来确定风机监测部件的故障情况,计算方法简单易行,能够准确地确定故障区域,在多部件复杂配合的风机装置振动影响分析中尤为有效。
为解决上述技术问题,本发明还提供了一种风机振动状态配合端等效激励确定方法,步骤如下:
按照设定采样周期对n个振动传感器在风机振动监测部位采集到的振动信号进行再采样;
根据得到的风机振动监测部位的振动信号,计算振动成分信号;
根据振动成分信号的幅值以及振动传感器采集到振动信号的时间,计算风机振动监测部位配合端的等效激励信号。
本发明的有益效果是:通过将风机振动监测部位的振动信号均视为是通过在振动监测部位的配合端施加激励引起的,根据采集到的振动信号,可以计算出振动监测部位配合端的等效激励信号,利用该等效激励信号可以实现很多功能,例如利用等效激励信号的周期性,就能够准确地确定风机监测部件的故障区域。
为解决上述技术问题,本发明还提供了一种风机振动状态检测系统,包括处理器、存储器以及设置在风机振动监测部位的振动传感器,所述处理器连接存储器和振动传感器,用于:按照设定采样周期接收n个振动传感器发送过来的在风机振动监测部位采集到的振动信号,并将其存储在存储器中;根据存储器中存储的振动信号,计算振动成分信号;根据振动成分信号的幅值以及振动传感器采集到振动信号的时间,计算风机振动监测部位配合端的等效激励信号;根据风机振动监测部位配合端的等效激励信号的周期性,确定风机振动监测部位的故障区域。
本发明的有益效果是:通过将风机振动监测部位的振动信号均视为是通过在监测部位的配合端施加激励引起的,这样根据采集到的振动信号计算振动监测部位配合端的等效激励信号,根据配合端的等效激励信号的周期性,就可以获知振动监测部位发生故障的位置,以确定故障范围。本发明通过计算配合端等效激励信号来确定风机监测部件的故障情况,计算方法简单易行,能够准确地确定故障区域,在多部件复杂配合的风机装置振动影响分析中尤为有效。
附图说明
图1为本发明风机振动状态监测分析系统的整体结构示意图;
图2为本发明风机振动状态检测方法的流程图;
图3为本发明振动监测部位振动传感器布置示意图;
图4为本发明配合端等效激励确定方法的流程示意图;
图5为本发明振动监测部位有限元仿真示意图;
图6为本发明振动状态监测分析的流程示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及具体实施例对本发明进行进一步详细说明。
风机振动状态检测方法实施例:
如图1所示,本实施例提供了一种风机振动状态监测分析系统,该系统包括:依次通讯连接的振动传感单元U1、数据处理单元U2、有限元仿真单元U3以及监测分析单元U4。其中,振动传感单元U1为在风机振动监测部位非均匀设置的多个振动传感器,以3个为例,其布置情况如图3所示,其中振动传感器1和2之间的距离L2>振动传感器2和3之间的距离L3。该基于有限元仿真分析的风机振动状态监测分析系统的工作原理如下:
振动传感单元U1将采集到的振动信号实时发送到数据处理单元U2。数据处理单元U2对接收到的振动信号进行再采样,经过振动激励等效方法计算得到振动监测部位配合端的等效激励信号,并实时传送到有限元仿真单元U3。有限元仿真单元U3将振动监测部位配合端的等效激励信号施加于有限元仿真模型的对应位置上,进行仿真,仿真得到振动监测部位的应变图像和各个有限元点的应变数据,并将振动监测部位的应变图像和各个有限元点的应变数据以及等效振动监测部位配合端激励压缩打包,实时传送给监测分析单元U4。监测分析单元U4将收到的打包信号解压缩,将其中的应变图像进行展示,并对配合端的等效激励信号和有限元点的应变数据进行处理分析来实现对振动监测部位振动造成的机械健康状态的诊断。
在上述风机振动状态监测分析系统的基础上,本实施例还提供了一种风机振动状态检测方法,其对应的流程图如图2所示,具体包括以下步骤:
(1)采用n个振动传感器对风机振动监测部位的振动信号进行采集。
(2)按照设定采样周期对n个振动传感器在风机振动监测部位采集到的振动信号进行再采样,根据采样数据,采用振动激励等效方法计算风机振动监测部位的配合端的等效激励信号。
其中,振动激励等效方法的基本思想是:对于监测区域的振动信号,不论其振动源位置,均可视为是通过在监测区域的配合端施加激励引起的,即研究振动对监测区域影响时,可将振动等效为配合端的激励。通过采用振动激励等效方法,不需要采用复杂的算法对振动源进行定位,采用振动激励等效方法计算得到振动监测部位配合端的等效激励信号的步骤如下:
1)对振动传感器采集的振动信号进行再采样,设定采样周期为T采样,T<T采样<T+T′,T为振动信号在振动监测部位相距最远的两个振动传感器间的传播时间,T′为振动信号在风机振动监测部位的不同配合端和距离该配合端最近的振动传感器之间的传播时间中的最小值。
其中,如图3所示,振动信号在振动监测部位相距最远的两个振动传感器间的传播时间T为振动信号在振动传感器1和振动传感器3之间的传播时间。振动信号在风机振动监测部位的不同配合端和距离该配合端最近的振动传感器之间的传播时间中的最小值T′是指:若配合端1和振动传感器1之间的距离L1<配合端2和振动传感器3之间的距离L4,则T′为振动信号在配合端1和振动传感器1之间的传播时间;若配合端1和振动传感器1之间的距离L1>配合端2和振动传感器3之间的距离L4,则T′为振动信号在配合端2和振动传感器3之间的传播时间。通过对设定采样周期T采样进行设置,可以保证同一个振动信号在采样周期内只能在风机振动监测部位的振动传感器间传播一次。
2)判断每个设定采样周期内采集到的每个振动传感器的振动信号个数,此时分为三种情况:
情况1:
步骤A:若每个采样周期内采集到每个振动传感器的1个振动信号,即振动传感器1、2……n-1、n采集到的振动信号分别为X1(t)、X2(t)……Xn-1(t)、Xn(t),说明如果此设定采样周期内存在振动,振动是从其中一个配合端开始传输的。那么,采用互相关去噪方法对X1(t)、X2(t)……Xn-1(t)、Xn(t)进行去噪处理,提取振动成分信号S1(t),S2(t-t1,2)……Sn-1(t-t1,n-1),Sn(t-t1,n),振动信号与对应振动成分信号的关系为:
X1(t)=S1(t)+Z1(t)
X2(t)=S2(t-t1,2)+Z2(t)=λ1S1(t-t1,2)+Z2(t)
…
Xn-1(t)=Sn-1(t-t1,n-1)+Zn-1(t)=λn-2S1(t-t1,n-1)+Zn-1(t)
Xn(t)=Sn(t-t1,,n)+Zn(t)=λn-1S1(t-t1,n)+Zn(t)
其中,Z1(t)、Z2(t)……Zn-1(t)、Zn(t)为外界随机干扰信号,λ1、λ2……λn-2、λn-1为衰减因数,t1,2、t1,3……t1,n-1、t1,n为振动信号到达振动传感器1和其余传感器的时延差。
需要说明的是,在步骤A中,本实施例是采用互相关的方式对振动信号进行去噪处理,以提取中振动成分信号,也就是有效振动信号,对应的具体过程属于现有技术,此处不再赘述。当然,作为其他的实施方式,也可以采用现有技术中的中值滤波、均值滤波等方法来对振动信号进行滤波。
步骤B:根据提取的振动成分信号,得到振动成分信号的幅值V1、V2……Vn-1、Vn,对振动成分信号的幅值V1、V2……Vn-1、Vn以及对应振动传感器采集到振动信号的时间t1、t2……tn-1、tn进行判断。
步骤C:若V1>V2>……>Vn-1>Vn且t1<t2<……<tn-1<tn,则根据振动成分信号的幅值V1、V2和振动传感器1、2的距离L2,计算振动信号单位距离的衰减率β1,对应的计算公式为:
需要说明的是,在计算振动信号单位距离的衰减率β1时,利用振动成分信号的幅值V1、V2和对应振动传感器1、2的距离L2的目的是为了提高衰减率β1的准确性。作为其他的实施方式,可以使用振动成分信号的幅值V1、V2……Vn-1、Vn中任意两个振动成分信号的幅值以及对应两个振动传感器之间的距离,计算振动信号单位距离的衰减率。
步骤D:根据振动成分信号的幅值V1、振动信号单位距离的衰减率β1以及振动传感器1与其最近的配合端之间的距离L1,计算该最近配合端的等效激励信号V01,计算公式为:
需要说明的是,在计算等效激励信号V01时,使用振动成分信号的幅值V1以及振动传感器1与其最近的配合端之间的距离L1是为了提高等效激励信号V01的准确性。作为其他的实施方式,可以选择振动成分信号的幅值V1、V2……Vn-1、Vn中任意一个振动成分信号的幅值以及该选择的振动成分信号的幅值对应的振动传感器与目标配合端之间的距离,计算等效激励信号,这里的目标配合端就是指距离最大的一个振动成分信号的幅值对应的振动传感器最近的配合端。
另外,在情况1中,若V1<V2<……<Vn-1<Vn且t1>t2>……>tn-1>tn,则根据振动成分信号的幅值Vn-1、Vn和振动传感器n-1、n的距离Ln,计算振动信号单位距离的衰减率β2,对应的计算公式为:
根据振动成分信号的幅值Vn、振动信号单位距离的衰减率β2以及振动传感器n与其最近的配合端之间的距离Ln+1,计算该最近配合端的等效激励信号V02,计算公式为:
在情况1中,若不满足V1>V2>……>Vn-1>Vn且t1<t2<……<tn-1<tn,也不满足V1<V2<……<Vn-1<Vn且t1>t2>……>tn-1>tn,说明此采样周期内混有与振动信号同频的干扰信号,则丢弃数据。
情况2:
步骤a:若每个设定采样周期内采集到每个振动传感器的2个振动信号,即振动传感器1、2……n-1、n采集到的一组振动信号分别为X1,1(t)、X2,1(t)……Xn-1,1(t)、Xn,1(t),振动传感器1、2……n-1、n采集到的另外一组振动信号分别为X1,2(t)、X2,2(t)……Xn-1,2(t)、Xn,2(t),说明如果此设定采样周期内存在振动,振动是从两个配合端分别传输的,采用互相关去噪方法分别对X1,1(t)、X2,1(t)……Xn-1,1(t)、Xn,1(t)和X1,2(t)、X2,2(t)……Xn-1,2(t)、Xn,2(t)进行去噪处理,提取振动成分信号S1,1(t)、S2,1(t-t1,2)……Sn,1(t-t1,n)和S1,2(t)、S2,2(t+t1,2)……Sn,2(t+t1,n),t1,2、t1,3……t1,n-1、t1,n为振动信号到达振动传感器1和其余传感器的时延差。
需要说明的是,在步骤a中,本实施例也是采用互相关的方式对振动信号进行去噪处理,以提取中振动成分信号。作为其他的实施方式,同样也可以采用现有技术中的中值滤波、均值滤波等方法来对振动信号进行滤波。
步骤b:根据提取的振动成分信号,得到振动成分信号的幅值V1,1、V2,1……Vn-1,1、Vn,1和V1,2、V2,2……Vn-1,2、Vn,2,对振动成分信号的幅值V1,1、V2,1……Vn-1,1、Vn,1和V1,2、V2,2……Vn-1,2、Vn,2以及对应振动传感器采集到振动信号的时间t1,1、t2,1……tn-1,1、tn,1和t1,2、t2,2……tn-1,2、tn,2进行判断。
步骤c:若V1,1>V2,1>……>Vn-1,1>Vn,1,t1,1<t2,1<……<tn-1,1<tn,1且V1,2<V2,2<……<Vn-1,2<Vn,2,t1,2>t2,2>……>tn-1,2>tn,2,则根据振动成分信号的幅值V1,1、V2,1和振动传感器1、2的距离L2,计算振动信号单位距离的衰减率β1,根据振动成分信号的幅值Vn,2、Vn-1,2和振动传感器n-1、n的距离Ln,计算振动信号单位距离的衰减率β2。
其中,与情况1中的步骤C类似,在计算振动信号单位距离的衰减率β1和振动信号单位距离的衰减率β2时,作为其他的实施方式,也可以采用对应所有振动成分信号的幅值中的任意两个幅值以及这两个幅值对应的振动传感器的距离进行实现,此处不再赘述。
步骤d:根据振动成分信号的幅值V1,1、振动信号单位距离的衰减率β1以及振动传感器1与其最近的配合端之间的距离L1,计算该最近配合端的等效激励信号V01;根据振动成分信号的幅值Vn,2、振动信号单位距离的衰减率β2以及振动传感器n与其最近的配合端之间的距离Ln+1,计算该最近配合端的等效激励信号V02,对应的计算公式可以参考情况1。
其中,与情况1中的步骤D类似,在计算等效激励信号V01和等效激励信号V02时,作为其他的实施方式,也可以采用对应所有振动成分信号的幅值中的任意一个幅值以及这幅值对应的振动传感器与目标配合端的距离进行实现,此处不再赘述。
另外,在情况2中,若V1,1<V2,1<……<Vn-1,1<Vn,1,t1,1>t2,1>……>tn-1,1>tn,1且V1,2>V2,2>……>Vn-1,2>Vn,2,t1,2<t2,2<……<tn-1,2<tn,2,则根据振动成分信号的幅值Vn,1、Vn-1,1和振动传感器n-1、n的距离Ln,计算振动信号单位距离的衰减率β2,根据振动成分信号的幅值V1,2、V2,2和振动传感器1、2的距离L2,计算振动信号单位距离的衰减率β1。根据振动成分信号的幅值Vn,1、振动信号单位距离的衰减率β2以及振动传感器n与其最近的配合端之间的距离Ln+1,计算该最近配合端的等效激励信号V02。根据振动成分信号的幅值V1,2、振动信号单位距离的衰减率β1以及振动传感器1与其最近的配合端之间的距离L1,计算该最近配合端的等效激励信号V01,对应的计算公式可以参考情况1。
在情况2中,若不满足V1,1>V2,1>……>Vn-1,1>Vn,1,t1,1<t2,1<……<tn-1,1<tn,1且V1,2<V2,2<……<Vn-1,2<Vn,2,t1,2>t2,2>……>tn-1,2>tn,2,也不满足V1,1<V2,1<……<Vn-1,1<Vn,1,t1,1>t2,1>……>tn-1,1>tn,1且V1,2>V2,2>……>Vn-1,2>Vn,2,t1,2<t2,2<……<tn-1,2<tn,2,说明此设定采样周期内混有与振动信号同频的干扰信号,则丢弃数据。
情况3:不满足每个设定采样周期内采集到每个振动传感器的1个振动信号,也不满足每个设定采样周期内采集到每个振动传感器的2个振动信号,则说明此采样周期内存在纯干扰信号,丢弃数据。
为了使上述步骤1)和2)中的内容更加清楚明了,以图3中的布置有3个振动传感器的振动监测部位为例进行说明,对应的步骤流程图如图4所示:
按照设定采样周期对振动传感器采集到的振动信号进行再采样,根据每个设定采样周期内采集到的每个振动传感器的振动信号个数,分为三种情况:
情况1:如果一个设定采样周期内分别得到三个振动传感器的一个信号X1(t)、X2(t)、X3(t),则说明如果此设定采样周期内存在振动,振动是从配合端1或配合端2开始传输的。采用互相关去噪方法对X1(t)、X2(t)、X3(t)进行去噪处理后,提取振动成分信号S1(t)和S2(t-t1,2)、S3(t-t1,3),计算公式为:
X1(t)=S1(t)+Z1(t)
X2(t)=S2(t-t1,2)+Z2(t)=λ1S1(t-t1,2)+Z2(t)
X3(t)=S3(t-t1,3)+Z3(t)=λ2S1(t-t1,3)+Z3(t)
其中,Z1(t)、Z2(t)、Z3(t)为外界随机干扰信号,λ1为振动信号从振动传感器1到振动传感器2的衰减因数,λ2为振动信号从振动传感器1到振动传感器3的衰减因数,t1,2为振动信号到达振动传感器1和振动传感器2的时延差,t1,3为振动信号到达振动传感器1和振动传感器3的时延差。λ1、λ2可取大于或小于1的数;t1,2、t1,3可取正、负值。
根据提取的振动成分信号,得到振动成分信号的幅值V1、V2、V3,判断振动成分信号的幅值V1、V2、V3和振动传感器采集到振动信号的时间t1、t2、t3的关系:
如果V1>V2>V3且t1<t2<t3,则根据振动成分信号的幅值V1、V2和振动传感器1、2的距离L2计算振动信号单位距离的衰减率β1。根据振动成分信号的幅值V1、振动信号单位距离的衰减率β1和振动传感器1与配合端1之间的距离L1计算配合端1的等效激励信号V01,计算公式为:
如果V1<V2<V3且t1>t2>t3,则根据振动成分信号的幅值V2、V3和振动传感器2、3的距离L3计算振动信号单位距离的衰减率β2。根据振动成分信号的幅值V3、振动信号单位距离的衰减率β2和振动传感器3与配合端2之间的距离L4计算配合端2的等效激励信号V02,计算公式为:
另外,若不满足V1>V2>V3且t1<t2<t3,也不满足V1<V2<V3且t1>t2>t3,则说明此设定采样周期内混有与振动信号同频的干扰信号,丢弃数据。
情况2:如果一个设定采样周期内分别得到三个振动传感器的两个信号X1,1(t)、X2,1(t)、X3,1(t)和X1,2(t)、X2,2(t)、X3,2(t),则说明如果此设定采样周期内存在振动,振动是从两个配合端,即配合端1和配合端2分别传输的。采用互相关去噪方法对X1,1(t)、X2,1(t)、X3,1(t)进行去噪处理,得到振动成分信号S1,1(t)、S2,1(t-t1,2)、S3,1(t-t1,3),采用互相关去噪方法对X1,2(t)、X2,2(t)、X3,2(t)进行去噪处理,得到振动成分信号S1,2(t)、S2,2(t+t1,2)、S3,2(t+t1,3)。其中,t1,2为振动信号到达振动传感器1和振动传感器2的时延差,t1,3为振动信号到达振动传感器1和振动传感器3的时延差。
根据提取的振动成分信号,得到振动成分信号的幅值V1,1、V2,1、V3,1和V1,2、V2,2、V3,2,判断振动成分信号的幅值V1,1、V2,1、V3,1和V1,2、V2,2、V3,2以及振动传感器采集到振动信号的时间t1,1、t2,1、t3,1和t1,2、t2,2、t3,2的关系:
如果V1,1>V2,1>V3,1,t1,1<t2,1<t3,1且V1,2<V2,2<V3,2,t1,2>t2,2>t3,2,则根据振动成分信号的幅值V1,1、V2,1和振动传感器1、2的距离L2计算振动信号单位距离的衰减率β1,根据振动成分信号的幅值V3,2、V2,2和振动传感器2、3的距离L3计算振动信号单位距离的衰减率β2。根据振动成分信号的幅值V1,1、振动信号单位距离的衰减率β1和振动传感器1与配合端1之间的距离L1计算配合端1的等效激励信号V01;根据振动成分信号的幅值V3,2、振动信号单位距离的衰减率β2和振动传感器3与配合端2之间的距离L4计算配合端2的等效激励信号V02。
如果V1,1<V2,1<V3,1,t1,1>t2,1>t3,1且V1,2>V2,2>V3,2,t1,2<t2,2<t3,2,则根据振动成分信号的幅值V3,1、V2,1和振动传感器2、3的距离L3计算振动信号单位距离的衰减率β2,根据振动成分信号的幅值V1,2、V2,2和振动传感器1、2的距离L2计算振动信号单位距离的衰减率β1。根据振动成分信号的幅值V3,1、振动信号单位距离的衰减率β2和振动传感器3与配合端2之间的距离L4计算配合端2的等效激励信号V02;根据振动成分信号的幅值V1,2、振动信号单位距离的衰减率β1和振动传感器1与配合端1之间的距离L1计算配合端1的等效激励信号V01。
(3)将振动监测部位配合端的等效激励信号施加于有限元仿真模型的对应位置,仿真得到振动监测部位的应变图像和各个有限元点的应变数据。
其中,有限元仿真单元U3在有限元仿真软件中根据风机的实际情况,建立风机振动监测部位和配合部件的有限元模型。其中,现有的有限元仿真软件有很多,在本实施例中,所采用的有限元仿真软件为COMSOLMultiphysics。另外,有限元仿真单元U3在有限元仿真软件中建立风机振动监测部位和配合部件的有限元模型的具体过程属于现有技术,此处不再赘述。
然后将数据处理单元U2传送的经过振动激励等效方法计算得到的配合端的等效激励信号施加到有限元模型的对应位置上,进行应变仿真。将仿真得到的振动监测部位应变图像、各个有限元点的应变数据和数据处理单元U2传送的振动监测部位配合端的等效激励信号打包压缩,实时传送给监测分析单元U4。例如,图5给出了某次仿真得到的振动监测部位应变图像,该振动监测部位图像包括图5的左侧图形和右侧的条状图,条状图是对颜色表示应变大小的说明,相当于图例。
(4)对振动监测部位的应变图像进行展示,并对振动监测部位配合端的等效激励信号和有限元点的应变数据进行分析,得到振动监测部位对应的健康状态。
具体的,监测分析单元U4对有限元仿真单元U3传送过来的打包信号解压缩,将其中的应变图像进行展示,并对配合端的等效激励信号和有限元点的应变数据进行处理分析来实现对振动监测部位振动造成的机械健康状态的诊断。其中,对振动监测部位配合端的等效激励信号和有限元点的应变数据进行分析的步骤包括:
4.1)设定阈值,筛选出超过设定阈值的振动监测部位的配合端的等效激励信号。
其中,风机不同监测部位对应设定的阈值不同,需要根据现场试验和经验以及需要进行确定。
4.2)判断筛选出的振动监测部位的配合端的等效激励信号是否为周期性信号,此时分为三种情形:
第一种情形,若筛选出的振动监测部位的配合端的等效激励信号为两个,且两个配合端的等效激励信号均为周期性信号,则说明该振动监测部位所属零件存在故障,建议对此部件进行检查。
第二种情形:若筛选出的振动监测部位的配合端的等效激励信号中的其中一个为周期性信号(包括:筛选出的振动监测部位的配合端的等效激励信号只有一个且为周期性信号;筛选出的振动监测部位的配合端的等效激励信号有两个且只有一个为周期性信号),则判断是否满足与周期性配合端的等效激励信号所属配合端直接相连的部件的两个配合端的等效激励信号是周期性信号,若满足,则与周期性配合端的等效激励信号所属配合端直接相连的部件存在故障;否则,周期性配合端的等效激励信号所属配合端和与其直接相连的部件的配合区域存在故障。
第三种情形:若筛选出的振动监测部位的配合端的等效激励信号均不为周期性信号(包括:筛选出的振动监测部位的配合端的等效激励信号只有一个且不为周期性信号;筛选出的振动监测部位的配合端的等效激励信号有两个且均不为周期性信号),则将振动监测部位各有限元点仿真应变大小ε仿i与对应振动监测部位各点可承受额定应变ε额i进行比较:如果ε仿i>ε额i,则说明振动导致的该振动监测部位的健康状态为告警状态;如果ε仿i>η*ε额i,0<η<1,取η=0.8,则振动导致的该振动监测部位的健康状态为预警状态;否则,振动导致的该振动监测部位的健康状态为正常状态。
其中,振动监测部位各有限元点仿真应变大小ε仿i即为仿真得到的有限元点的应变数据,包含在有限元仿真单元U3传送的压缩包中。
例如,在图5中,振动监测部位为B,振动状态监测分析的步骤流程图如图6所示,此时,若筛选出的振动监测部位B的两个配合端的等效激励信号均为周期性信号,也就是等效配合端B1和B2的激励信号均为周期性信号,则说明振动监测部位B所属部件存在故障,建议对此部件进行检查。若振动监测部位B的两个配合端B1和B2的等效激励信号只有一个为周期性信号,则继续判断与周期性等效激励信号所属等效配合端B1或B2直接相连的部件A或C的两个配合端的等效激励信号是否是周期信号,如果是,则说明此直接相连的部件A或C存在故障,建议对部件进行检查;如果否,则说明配合区域,即A与B或B与C的配合区域存在故障,建议进行检查。
上述的风机振动状态检测方法基于振动传感器的采集信号,结合有限元分析,充分考虑了检测部件各部位的不同可承受能力,可以对监测部件整体的状态及状态变化进行直观的展示,并可以对监测部件的健康状态进行预警,解决了监测部件不全面,状态评估滞后性和展示不直观的缺陷。通过采用振动激励等效方法,不需要对振动源进行定位,算法简单,在多部件复杂配合的风机装置振动影响分析中尤为有效,并去掉了环境干扰信号的影响,使监测分析更准确。
另外,需要说明的是,上述的风机振动状态检测方法不仅可以应用在风机的故障检测过程中,同样可以应用在其他设备的故障检测过程中。
风机振动状态检测装置实施例:
本实施例提供了一种风机振动状态检测装置,包括处理器和存储器,该处理器用于处理存储在存储器中的指令,以实现一种风机振动状态检测方法,由于该风机振动状态检测方法已经在上述的风机振动状态检测方法实施例中进行了详细介绍,此处不再赘述。
风机振动状态配合端的等效激励确定方法实施例:
本实施例提供了一种风机振动状态配合端的等效激励确定方法,用于确定风机振动监测部位两个配合端的等效激励信号。当然,该方法不仅可以适用于风机振动监测部位的配合端等效激励信号确定,还可以适用于其他设备的振动监测部位的配合端等效激励信号确定。并且,计算出的配合端的等效激励信号不仅可以应用于监测部件的故障检测过程中,还有其他的用途。由于该风机振动状态配合端的等效激励确定方法已经在上述的风机振动状态检测方法实施例中进行了详细介绍,此处不再赘述。
风机振动状态检测系统实施例:
本实施例提供了一种风机振动状态检测系统,包括处理器、存储器以及设置在风机振动监测部位的振动传感器,处理器连接存储器和振动传感器,用于实现一种风机振动状态检测方法。由于该风机振动状态检测系统的工作过程已经在上述的风机振动状态检测方法实施例中进行了详细介绍,此处不再赘述。
最后应当说明的是,以上实施例仅用于说明本发明的技术方案而非对其保护范围的限制,尽管参照上述实施例对本申请进行了详细的说明,所属领域的普通技术人员应当理解,本领域技术人员阅读本申请后依然可对申请的具体实施方式进行种种变更、修改或者等同替换,但这些变更、修改或者等同替换,均在本发明的权利要求保护范围之内。
Claims (7)
1.一种风机振动状态检测方法,其特征在于,步骤如下:
(1)按照设定采样周期对n个振动传感器在风机振动监测部位采集到的振动信号进行再采样;
(2)根据得到的风机振动监测部位的振动信号,计算振动成分信号;
(3)根据振动成分信号的幅值以及振动传感器采集到振动信号的时间,计算风机振动监测部位配合端的等效激励信号;
(4)根据风机振动监测部位配合端的等效激励信号的周期性,确定风机振动监测部位的故障区域;
步骤(3)中计算风机振动监测部位配合端的等效激励信号的步骤为:
若n个振动传感器在每个设定采样周期内均采集到1个振动信号X1(t)、X2(t)……Xn-1(t)、Xn(t),并且满足V1>V2>……>Vn-1>Vn且t1<t2<……<tn-1<tn或满足V1<V2<……<Vn-1<Vn且t1>t2>……>tn-1>tn,其中,Vi为第i个振动传感器采集到的振动信号计算出的振动成分信号的幅值,ti为第i个振动传感器采集到振动信号的时间,则根据其中两个振动成分信号的幅值以及对应两个振动传感器之间的距离,计算振动信号单位距离的衰减率;
根据选择的其中一个振动成分信号的幅值、振动信号单位距离的衰减率以及该选择的其中一个振动成分信号的幅值对应的振动传感器与目标配合端之间的距离,计算该目标配合端的等效激励信号,目标配合端为距离最大的一个振动成分信号的幅值对应的振动传感器最近的配合端。
2.根据权利要求1所述的风机振动状态检测方法,其特征在于,计算风机振动监测部位配合端的等效激励信号的步骤为:
若n个振动传感器在每个设定采样周期内均采集到2个振动信号X1,1(t)、X2,1(t)……Xn-1,1(t)、Xn,1(t)和X1,2(t)、X2,2(t)……Xn-1,2(t)、Xn,2(t),并且满足V1,1>V2,1>……>Vn-1,1>Vn,1,t1,1<t2,1<……<tn-1,1<tn,1且V1,2<V2,2<……<Vn-1,2<Vn,2,t1,2>t2,2>……>tn-1,2>tn,2或满足V1,1<V2,1<……<Vn-1,1<Vn,1,t1,1>t2,1>……>tn-1,1>tn,1且V1,2>V2,2>……>Vn-1,2>Vn,2,t1,2<t2,2<……<tn-1,2<tn,2,其中,Vi,1为第i个振动传感器采集到的第一振动信号计算出的振动成分信号的幅值,ti,1为第i个振动传感器采集到第一振动信号的时间,Vi,2为第i个振动传感器采集到的第二振动信号计算出的振动成分信号的幅值,ti,2为第i个振动传感器采集到第二振动信号的时间,则根据各振动成分信号的幅值V1,1、V2,1……Vn-1,1、Vn,1中的两个振动成分信号的幅值以及对应两个振动传感器之间的距离,计算第一振动信号单位距离的衰减率,根据各振动成分信号的幅值V1,2、V2,2……Vn-1,2、Vn,2中的两个振动成分信号的幅值以及对应两个振动传感器之间的距离,计算第二振动信号单位距离的衰减率;
根据各振动成分信号的幅值V1,1、V2,1……Vn-1,1、Vn,1中选择的一个振动成分信号的幅值、第一振动信号单位距离的衰减率以及该各振动成分信号的幅值V1,1、V2,1……Vn-1,1、Vn,1中选择的一个振动成分信号的幅值对应的振动传感器与第一目标配合端之间的距离,计算该第一目标配合端的等效激励信号,第一目标配合端为距离各振动成分信号的幅值V1,1、V2,1……Vn-1,1、Vn,1中最大的一个振动成分信号的幅值对应的振动传感器最近的配合端;根据各振动成分信号的幅值V1,2、V2,2……Vn-1,2、Vn,2中选择的一个振动成分信号的幅值、第二振动信号单位距离的衰减率以及该各振动成分信号的幅值V1,2、V2,2……Vn-1,2、Vn,2中选择的一个振动成分信号的幅值对应的振动传感器与第二目标配合端之间的距离,计算该第二目标配合端的等效激励信号,第二目标配合端为距离各振动成分信号的幅值V1,2、V2,2……Vn-1,2、Vn,2中最大的一个振动成分信号的幅值对应的振动传感器最近的配合端。
3.根据权利要求2所述的风机振动状态检测方法,其特征在于,若振动监测部位的两个配合端的等效激励信号的幅值均大于设定阈值,且两个配合端的等效激励信号均为周期性信号,则判定该振动监测部位所属部件存在故障。
4.根据权利要求1或2所述的风机振动状态检测方法,其特征在于,若振动监测部位的至少一个配合端的等效激励信号的幅值大于设定阈值,且在大于设定阈值的等效激励信号中只有第一配合端的等效激励信号为周期性信号,则判断是否满足与该第一配合端直接相连的部件的两个等效激励信号是周期性信号;若满足,则与该第一配合端直接相连的部件存在故障,否则,该第一配合端和与其直接相连的部件的配合区域存在故障。
5.根据权利要求1或2所述的风机振动状态检测方法,其特征在于,还包括:将振动监测部位的配合端的等效激励信号施加于有限元仿真模型的对应位置,仿真得到振动监测部位的应变图像;若振动监测部位的配合端的等效激励信号的幅值大于设定阈值,且大于设定阈值的等效激励信号均不为周期性信号,则根据仿真得到振动监测部位各有限元点仿真应变与振动监测部位各点可承受额定应变的关系确定振动监测部位的健康状态。
6.根据权利要求1-2中任一项所述的风机振动状态检测方法,其特征在于,振动信号在振动监测部位相距最远的两个振动传感器间的传播时间<所述设定采样周期<振动信号在振动监测部位相距最远的两个振动传感器间的传播时间+振动信号在风机振动监测部位的不同配合端和距离该配合端最近的振动传感器之间的传播时间中的最小值。
7.一种风机振动状态检测装置,其特征在于,包括处理器和存储器以及设置在风机振动监测部位的振动传感器,所述处理器用于处理存储在存储器中的指令,以实现如权利要求1-6中任一项所述的风机振动状态检测方法;处理器连接存储器和振动传感器,用于:按照设定采样周期接收n个振动传感器发送过来的在风机振动监测部位采集到的振动信号,并将其存储在存储器中;根据存储器中存储的振动信号,计算振动成分信号;根据振动成分信号的幅值以及振动传感器采集到振动信号的时间,计算风机振动监测部位配合端的等效激励信号;根据风机振动监测部位配合端的等效激励信号的周期性,确定风机振动监测部位的故障区域。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911067941.6A CN110793615B (zh) | 2019-11-04 | 2019-11-04 | 风机振动状态检测方法、装置及配合端等效激励确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911067941.6A CN110793615B (zh) | 2019-11-04 | 2019-11-04 | 风机振动状态检测方法、装置及配合端等效激励确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110793615A CN110793615A (zh) | 2020-02-14 |
CN110793615B true CN110793615B (zh) | 2021-09-14 |
Family
ID=69442650
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911067941.6A Active CN110793615B (zh) | 2019-11-04 | 2019-11-04 | 风机振动状态检测方法、装置及配合端等效激励确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110793615B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101545459A (zh) * | 2008-03-25 | 2009-09-30 | 唐德尧 | 风力发电机塔架运行状态监控装置 |
CN103049670A (zh) * | 2013-01-14 | 2013-04-17 | 哈尔滨工程大学 | 管道激励源识别及其振动响应预测方法 |
CN104048825A (zh) * | 2014-05-28 | 2014-09-17 | 北京航空航天大学 | 一种多传感器融合的风电机组齿轮箱故障定位方法 |
US8990031B2 (en) * | 2009-12-10 | 2015-03-24 | Industrial Technology Research Institute | Torsional resonance frequency measuring device and method |
CN104484827A (zh) * | 2014-09-29 | 2015-04-01 | 许继集团有限公司 | 一种风力发电机组的故障频率成分提取方法 |
CN104915498A (zh) * | 2015-06-08 | 2015-09-16 | 广东工业大学 | 基于模型识别与等效简化的高速平台运动参数自整定方法 |
CN110186679A (zh) * | 2019-05-13 | 2019-08-30 | 温州大学 | 盾构主轴轴承的诊断方法 |
-
2019
- 2019-11-04 CN CN201911067941.6A patent/CN110793615B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101545459A (zh) * | 2008-03-25 | 2009-09-30 | 唐德尧 | 风力发电机塔架运行状态监控装置 |
US8990031B2 (en) * | 2009-12-10 | 2015-03-24 | Industrial Technology Research Institute | Torsional resonance frequency measuring device and method |
CN103049670A (zh) * | 2013-01-14 | 2013-04-17 | 哈尔滨工程大学 | 管道激励源识别及其振动响应预测方法 |
CN104048825A (zh) * | 2014-05-28 | 2014-09-17 | 北京航空航天大学 | 一种多传感器融合的风电机组齿轮箱故障定位方法 |
CN104484827A (zh) * | 2014-09-29 | 2015-04-01 | 许继集团有限公司 | 一种风力发电机组的故障频率成分提取方法 |
CN104915498A (zh) * | 2015-06-08 | 2015-09-16 | 广东工业大学 | 基于模型识别与等效简化的高速平台运动参数自整定方法 |
CN110186679A (zh) * | 2019-05-13 | 2019-08-30 | 温州大学 | 盾构主轴轴承的诊断方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110793615A (zh) | 2020-02-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104048825B (zh) | 一种多传感器融合的风电机组齿轮箱故障定位方法 | |
EP2476901A1 (en) | Method and monitoring apparatus for automated surveillance of a wind turbine and a method for creating a linear model | |
CN112834193B (zh) | 一种基于三维图的运营桥梁振动和健康状态异常预警方法 | |
JP2004523735A (ja) | 器械を取り付けられた機器およびプロセスのための推論信号生成装置 | |
EP2345894A2 (en) | Trending of vibration data taking into account torque effect | |
Samali et al. | Location and severity identification of notch-type damage in a two-storey steel framed structure utilising frequency response functions and artificial neural network | |
KR101129466B1 (ko) | 웨이블렛 스레스홀딩 함수와 4차 모멘트를 이용한 회전기계의 상태 진단 방법 | |
CN113569445A (zh) | 一种基于数字孪生技术的钢结构健康监测系统及方法 | |
CN101504324A (zh) | 拉索运营状态智能判别方法及系统 | |
CN108444662A (zh) | 一种基于日温度效应的桥梁损伤在线监测方法 | |
JP2012018045A (ja) | センサ異常診断装置及びセンサシステム | |
CN113467433A (zh) | 一种机械设备故障检测的方法及装置 | |
CN109470972B (zh) | 一种基于振动图像的变压器铁心松动诊断方法 | |
CN111311872A (zh) | 一种船体结构应力长期监测报警系统 | |
CN110990978A (zh) | 一种螺栓状态监测方法和装置 | |
CN114154779A (zh) | 一种基于emd-som的健康评估与故障诊断方法及系统 | |
CN104870952A (zh) | 土木结构的完整性 | |
CN112594356B (zh) | 一种高铁与地铁齿轮箱运行监测和故障诊断系统 | |
CN114739667A (zh) | 一种多模态信息融合轴承润滑状态监测装置及方法 | |
Liu et al. | Research on online monitoring technology for transmission tower bolt looseness | |
KR101046748B1 (ko) | 발전기 고정자 권선의 전용 진동 진단 방법 및 분석 시스템 | |
CN110793615B (zh) | 风机振动状态检测方法、装置及配合端等效激励确定方法 | |
CN102072926A (zh) | 一种发动机机体疲劳裂纹诊断的方法 | |
CN113485898A (zh) | 振动测点显示方法、装置、设备及存储介质 | |
CN110598680A (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 |