CN110793615B - Fan vibration state detection method and device and matching end equivalent excitation determination method - Google Patents
Fan vibration state detection method and device and matching end equivalent excitation determination method 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
Abstract
The invention relates to a method and a device for detecting the vibration state of a fan and a method for determining equivalent excitation of a matching end, wherein the method comprises the following steps: sampling vibration signals collected by the n vibration sensors at the vibration monitoring part of the fan according to a set sampling period; calculating a vibration component signal according to a vibration signal of the fan vibration monitoring part obtained by sampling; calculating an equivalent excitation signal at the matching end of the vibration monitoring part of the fan according to the amplitude of the vibration component signal and the time for the vibration sensor to acquire the vibration signal; and determining the fault area of the fan vibration monitoring part according to the periodicity of the equivalent excitation signal at the matching end of the fan vibration monitoring part. The invention determines the fault condition of the fan monitoring component by calculating the equivalent excitation signal of the matching end, the calculation method is simple and easy to implement, and the fault area can be accurately determined.
Description
Technical Field
The invention relates to a method and a device for detecting a vibration state of a fan and a method for determining equivalent excitation of a matching end, and belongs to the technical field of wind power generation.
Background
As a clean and renewable new energy source, wind power is rapidly developed, and especially, offshore wind power development has become a global hotspot and frontier. The wind driven generator (called a fan for short) is used as a key device of wind power generation, and the wind driven generator inevitably breaks down due to complex operating environment, so that the generating efficiency is seriously influenced if the wind driven generator is abnormally stopped. Therefore, it is necessary and urgent to monitor the critical operating signals of the wind turbine, analyze and diagnose the fault of the wind turbine, and determine the fault area. The faults of the fan mainly include electrical faults and mechanical faults, the monitoring, analyzing and diagnosing system related to the electrical signals of the fan is mature at present, and the monitoring, analyzing and fault diagnosing technology related to the mechanical signals of the fan is relatively deficient. Vibration-induced faults are most common in mechanical faults of wind turbines, and therefore it is necessary and urgent to monitor and analyze the mechanical health of the wind turbine in response to the effects of vibration of critical components of the wind turbine.
In the prior art, the monitoring and analysis of the mechanical health state of the vibration influence of the fan component are realized by arranging a vibration sensor at a vibration monitoring part and transmitting the acquired vibration data to a master control PLC (programmable logic controller) and a server connected with the PLC through a cable or an optical fiber. In the prior art, the vibration monitoring and analysis of the fan part is carried out on vibration signals acquired by a vibration sensor, and the overall influence of vibration on a monitoring part is not involved. In fact, the influence of the vibration on the monitoring component is different due to different bearing capacities of different parts, and therefore, the finally determined fault area of the fan monitoring part is inaccurate only according to the collected partial vibration signals. Because the health of the fan components is evaluated based on the vibration signals collected by the sensors, the evaluation hysteresis is caused, and the efficiency of the fan power generation is seriously influenced when or after the fault occurs. In addition, interference signals are attached to vibration signals acquired by the vibration sensor, so that a large error exists in analysis based on the acquired signals, and the judgment of a fault area is also influenced. In addition, the health state of the fan monitoring part is displayed in a data and signal mode in the prior art, and the displaying is not visual. In addition, most of the existing analysis on the influence caused by vibration adopts a vibration positioning method, but the method is obviously not suitable for a device with complex matching of components, such as a fan, and the vibration positioning method is generally complex in calculation process.
Disclosure of Invention
The invention aims to provide a fan vibration state detection method, a fan vibration state detection device and a matching end equivalent excitation determination method, which are used for solving the problem that the determination of a fan detection component fault area is inaccurate in the prior art.
In order to solve the technical problem, the invention provides a method for detecting the vibration state of a fan, which comprises the following steps:
resampling the vibration signals collected by the n vibration sensors at the vibration monitoring part of the fan according to a set sampling period;
calculating a vibration component signal according to the obtained vibration signal of the vibration monitoring part of the fan;
calculating an equivalent excitation signal at the matching end of the vibration monitoring part of the fan according to the amplitude of the vibration component signal and the time for the vibration sensor to acquire the vibration signal;
and determining the fault area of the fan vibration monitoring part according to the periodicity of the equivalent excitation signal at the matching end of the fan vibration monitoring part.
The invention has the beneficial effects that: the vibration signals of the fan vibration monitoring part are all considered to be caused by applying excitation on the matching end of the monitoring part, so that the equivalent excitation signal of the matching end of the vibration monitoring part is calculated according to the collected vibration signals, and the position of the vibration monitoring part with a fault can be obtained according to the periodicity of the equivalent excitation signal of the matching end, so that the fault range can be determined. The method determines the fault condition of the fan monitoring component by calculating the equivalent excitation signal of the matching end, is simple and feasible, can accurately determine the fault area, and is particularly effective in the analysis of the vibration influence of the fan device with complex matching of multiple components.
Further, in order to reliably obtain the equivalent excitation signal so as to accurately obtain the health state of the fan vibration monitoring part, the step of calculating the equivalent excitation signal at the matching end of the fan vibration monitoring part is as follows:
if the n vibration sensors acquire 1 vibration signal X in each set sampling period1(t)、X2(t)……Xn-1(t)、Xn(t) and satisfy V1>V2>……>Vn-1>VnAnd t is1<t2<……<tn-1<tnOr satisfy V1<V2<……<Vn-1<VnAnd t is1>t2>……>tn-1>tnWherein V isiAmplitude, t, of vibration component signal calculated for vibration signal collected by ith vibration sensoriFor the ith vibration sensor, vibration is collectedCalculating the attenuation rate of the unit distance of the vibration signal according to the amplitude of two vibration component signals and the distance between the two corresponding vibration sensors;
and calculating the equivalent excitation signal of the target matching end according to the selected amplitude of one vibration component signal, the attenuation rate of the unit distance of the vibration signal and the distance between the vibration sensor corresponding to the selected amplitude of one vibration component signal and the target matching end, wherein the target matching end is the closest matching end of the vibration sensor corresponding to the amplitude of the vibration component signal with the largest distance.
Further, in order to reliably obtain the equivalent excitation signal so as to accurately obtain the health state of the fan vibration monitoring part, the step of calculating the equivalent excitation signal at the matching end of the fan vibration monitoring part is as follows:
if the n vibration sensors acquire 2 vibration signals X in each set sampling period1,1(t)、X2,1(t)……Xn-1,1(t)、Xn,1(n) and X1,2(t)、X2,2(t)……Xn-1,2(t)、Xn,2(t) and satisfy V1,1>V2,1>……>Vn-1,1>Vn,1，t1,1<t2,1<……<tn-1,1<tn,1And V is1,2<V2,2<……<Vn-1,2<Vn,2，t1,2>t2,2>……>tn-1,2>tn,2Or satisfy V1,1<V2,1<……<Vn-1,1<Vn,1，t1,1>t2,1>……>tn-1,1>tn,1And V is1,2>V2,2>……>Vn-1,2>Vn,2，t1,2<t2,2<……<tn-1,2<tn,2Wherein V isi,1Amplitude, t, of vibration component signal calculated for first vibration signal collected by ith vibration sensori,1Time, V, of acquisition of first vibration signal for ith vibration sensori,2Collected for the ith vibration sensorAmplitude, t, of vibration component signal calculated from two vibration signalsi,2The time when the second vibration signal is acquired by the ith vibration sensor is determined according to the amplitude V of each vibration component signal1,1、V2,1……Vn-1,1、Vn,1Calculates the attenuation rate of the first vibration signal per unit distance based on the amplitudes V of the vibration component signals and the distances between the two vibration sensors1,2、V2,2……Vn-1,2、Vn,2Calculating the attenuation rate of the unit distance of the second vibration signal according to the amplitude of the two vibration component signals and the distance between the two corresponding vibration sensors;
according to the amplitude V of each vibration component signal1,1、V2,1……Vn-1,1、Vn,1The amplitude of the selected one vibration component signal, the attenuation rate of the first vibration signal per unit distance, and the amplitude V of each vibration component signal1,1、V2,1……Vn-1,1、Vn,1Calculating the equivalent excitation signal of the first target mating end corresponding to the distance between the vibration sensor and the first target mating end corresponding to the amplitude of the selected vibration component signal, wherein the first target mating end is the amplitude V of each vibration component signal1,1、V2,1……Vn-1,1、Vn,1The closest mating end of the vibration sensor corresponding to the amplitude of the largest vibration component signal; according to the amplitude V of each vibration component signal1,2、V2,2……Vn-1,2、Vn,2The amplitude of the selected one vibration component signal, the attenuation rate of the second vibration signal per unit distance, and the amplitude V of each vibration component signal1,2、V2,2……Vn-1,2、Vn,2Calculating the equivalent excitation signal of the second target matching end by the distance between the vibration sensor corresponding to the amplitude of the selected vibration component signal and the second target matching end, wherein the second target matching end is the amplitude V of each vibration component signal1,2、V2,2……Vn-1,2、Vn,2Vibration corresponding to the amplitude of the largest vibration component signalThe nearest mating end of the motion sensor.
Further, in order to reliably determine a fault area affecting the health state of the vibration monitoring portion, if the amplitudes of the equivalent excitation signals at the two mating ends of the vibration monitoring portion are both greater than a set threshold value and the equivalent excitation signals at the two mating ends are both periodic signals, it is determined that a component to which the vibration monitoring portion belongs has a fault.
Further, in order to reliably determine a fault area affecting the health state of the vibration monitoring portion, if the amplitude of the equivalent excitation signal of at least one of the mating ends of the vibration monitoring portion is greater than a set threshold, and only the equivalent excitation signal of the first mating end is a periodic signal among the equivalent excitation signals greater than the set threshold, it is determined whether two equivalent excitation signals of a component directly connected to the first mating end are satisfied as periodic signals; if so, a component directly connected to the first mating end fails, otherwise, the first mating end and the mating area of the component directly connected thereto fail.
Furthermore, in order to facilitate visually seeing the vibration influence condition of the vibration monitoring part, the method further comprises the following steps: applying the equivalent excitation signal of the matching end of the vibration monitoring part to the corresponding position of the finite element simulation model, and simulating to obtain a strain image of the vibration monitoring part; and if the amplitude of the equivalent excitation signal at the matching end of the vibration monitoring part is greater than a set threshold value and the equivalent excitation signals greater than the set threshold value are not periodic signals, determining the health state of the vibration monitoring part according to the relationship between the simulation strain of each finite element point of the vibration monitoring part and the rated strain bearable of each point of the vibration monitoring part obtained by simulation.
Further, in order to ensure that the same vibration signal can only be transmitted once between the sensors of the vibration monitoring portion in the sampling period, the transmission time of the vibration signal between the two vibration sensors at the vibration monitoring portion with the farthest distance is less than the set sampling period < the transmission time of the vibration signal between the two vibration sensors at the vibration monitoring portion with the farthest distance + the minimum value of the transmission time of the vibration signal between the different mating ends of the fan vibration monitoring portion and the vibration sensor closest to the mating end.
In order to solve the technical problem, the invention further provides a fan vibration state detection device, which comprises a processor and a memory, wherein the processor is used for processing the instructions stored in the memory so as to realize the fan vibration state detection method.
The invention has the beneficial effects that: the vibration signals of the fan vibration monitoring part are all considered to be caused by applying excitation on the matching end of the monitoring part, so that the equivalent excitation signal of the matching end of the vibration monitoring part is calculated according to the collected vibration signals, and the position of the vibration monitoring part with a fault can be obtained according to the periodicity of the equivalent excitation signal of the matching end, so that the fault range can be determined. The method determines the fault condition of the fan monitoring component by calculating the equivalent excitation signal of the matching end, is simple and feasible, can accurately determine the fault area, and is particularly effective in the analysis of the vibration influence of the fan device with complex matching of multiple components.
In order to solve the technical problem, the invention also provides a method for determining the equivalent excitation of the matching end of the fan in the vibration state, which comprises the following steps:
resampling the vibration signals collected by the n vibration sensors at the vibration monitoring part of the fan according to a set sampling period;
calculating a vibration component signal according to the obtained vibration signal of the vibration monitoring part of the fan;
and calculating an equivalent excitation signal at the matching end of the vibration monitoring part of the fan according to the amplitude of the vibration component signal and the time for acquiring the vibration signal by the vibration sensor.
The invention has the beneficial effects that: the vibration signals of the fan vibration monitoring part are regarded as being caused by applying excitation to the matching end of the vibration monitoring part, the equivalent excitation signal of the matching end of the vibration monitoring part can be calculated according to the collected vibration signals, and a plurality of functions can be realized by utilizing the equivalent excitation signal, for example, the fault area of the fan monitoring part can be accurately determined by utilizing the periodicity of the equivalent excitation signal.
In order to solve the technical problem, the invention further provides a fan vibration state detection system, which comprises a processor, a memory and a vibration sensor arranged at a fan vibration monitoring part, wherein the processor is connected with the memory and the vibration sensor and is used for: receiving vibration signals which are sent by n vibration sensors and collected at a fan vibration monitoring part according to a set sampling period, and storing the vibration signals in a memory; calculating a vibration component signal according to the vibration signal stored in the memory; calculating an equivalent excitation signal at the matching end of the vibration monitoring part of the fan according to the amplitude of the vibration component signal and the time for the vibration sensor to acquire the vibration signal; and determining the fault area of the fan vibration monitoring part according to the periodicity of the equivalent excitation signal at the matching end of the fan vibration monitoring part.
The invention has the beneficial effects that: the vibration signals of the fan vibration monitoring part are all considered to be caused by applying excitation on the matching end of the monitoring part, so that the equivalent excitation signal of the matching end of the vibration monitoring part is calculated according to the collected vibration signals, and the position of the vibration monitoring part with a fault can be obtained according to the periodicity of the equivalent excitation signal of the matching end, so that the fault range can be determined. The method determines the fault condition of the fan monitoring component by calculating the equivalent excitation signal of the matching end, is simple and feasible, can accurately determine the fault area, and is particularly effective in the analysis of the vibration influence of the fan device with complex matching of multiple components.
Drawings
FIG. 1 is a schematic view of the overall structure of a fan vibration state monitoring and analyzing system according to the present invention;
FIG. 2 is a flow chart of a method for detecting a vibration state of a fan according to the present invention;
FIG. 3 is a schematic view of the vibration sensor arrangement of the vibration monitoring section of the present invention;
FIG. 4 is a flow chart diagram of a mating end equivalent excitation determination method of the present invention;
FIG. 5 is a schematic diagram of a finite element simulation of a vibration monitoring portion according to the present invention;
FIG. 6 is a flow chart of the vibration status monitoring and analyzing method according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail with reference to the accompanying drawings and specific embodiments.
The embodiment of the fan vibration state detection method comprises the following steps:
as shown in fig. 1, the present embodiment provides a fan vibration state monitoring and analyzing system, which includes: the vibration sensing unit U1, the data processing unit U2, the finite element simulation unit U3 and the monitoring analysis unit U4 are sequentially connected in a communication mode. The vibration sensing unit U1 is a plurality of vibration sensors disposed non-uniformly at the vibration monitoring position of the fan, 3 are taken as an example, and the arrangement is shown in fig. 3, wherein the distance L2 between the vibration sensors 1 and 2> the distance L3 between the vibration sensors 2 and 3. The working principle of the fan vibration state monitoring and analyzing system based on finite element simulation analysis is as follows:
the vibration sensing unit U1 sends the acquired vibration signals to the data processing unit U2 in real time. The data processing unit U2 resamples the received vibration signals, calculates the equivalent excitation signals of the matching end of the vibration monitoring part by a vibration excitation equivalent method, and transmits the equivalent excitation signals to the finite element simulation unit U3 in real time. The finite element simulation unit U3 applies the equivalent excitation signal at the matching end of the vibration monitoring part to the corresponding position of the finite element simulation model for simulation, obtains the strain image of the vibration monitoring part and the strain data of each finite element point through simulation, and performs excitation, compression and packaging on the strain image of the vibration monitoring part, the strain data of each finite element point and the matching end of the equivalent vibration monitoring part, and transmits the strain images, the strain data of each finite element point and the matching end of the equivalent vibration monitoring part to the monitoring and analyzing unit U4 in real time. The monitoring and analyzing unit U4 decompresses the received packaged signal, displays the strain image therein, and processes and analyzes the equivalent excitation signal at the mating end and the strain data at the finite element point to realize the diagnosis of the mechanical health state caused by the vibration of the vibration monitoring portion.
On the basis of the system for monitoring and analyzing the vibration state of the fan, the embodiment further provides a method for detecting the vibration state of the fan, and the corresponding flow chart is shown in fig. 2, and the method specifically comprises the following steps:
(1) and acquiring vibration signals of the vibration monitoring part of the fan by adopting n vibration sensors.
(2) And resampling the vibration signals collected by the n vibration sensors at the fan vibration monitoring part according to a set sampling period, and calculating the equivalent excitation signal of the matching end of the fan vibration monitoring part by adopting a vibration excitation equivalent method according to sampling data.
The basic idea of the vibration excitation equivalent method is as follows: the vibration signal of the monitoring area, regardless of the position of the vibration source, can be considered to be caused by applying excitation at the mating end of the monitoring area, i.e. when the influence of the vibration on the monitoring area is studied, the vibration can be equivalent to the excitation of the mating end. By adopting the vibration excitation equivalent method, the vibration source does not need to be positioned by adopting a complex algorithm, and the step of calculating the equivalent excitation signal of the matching end of the vibration monitoring part by adopting the vibration excitation equivalent method comprises the following steps:
1) resampling the vibration signal collected by the vibration sensor, and setting the sampling period as TSampling，T＜TSamplingT + T 'is the propagation time of the vibration signal between two vibration sensors which are farthest away from the vibration monitoring part, and T' is the minimum value of the propagation time of the vibration signal between different matching ends of the fan vibration monitoring part and the vibration sensor which is closest to the matching end.
As shown in fig. 3, the propagation time T of the vibration signal between the two vibration sensors most distant from the vibration monitoring portion is the propagation time of the vibration signal between the vibration sensor 1 and the vibration sensor 3. The minimum value T' in the propagation time of the vibration signal between different matching ends of the fan vibration monitoring part and the vibration sensor closest to the matching end is as follows: if the distance L1 between the mating end 1 and the vibration sensor 1 is<The distance L4 between the mating end 2 and the vibration sensor 3, then T' is the propagation time of the vibration signal between the mating end 1 and the vibration sensor 1; if the distance L1 between the mating end 1 and the vibration sensor 1 is>The distance L4 between the mating end 2 and the vibration sensor 3, T' is the vibration signal at the mating end 2 and the vibration transmissionThe propagation time between sensors 3. By setting a sampling period TSamplingThe vibration sensor is arranged, so that the same vibration signal can only be transmitted once between the vibration sensors at the vibration monitoring part of the fan in the sampling period.
2) The number of the vibration signals of each vibration sensor collected in each set sampling period is judged, and at the moment, the three conditions are divided into three conditions:
case 1:
step A: if 1 vibration signal of each vibration sensor is acquired in each sampling period, namely the vibration signals acquired by the vibration sensors 1, 2 … … n-1 and n are respectively X1(t)、X2(t)……Xn-1(t)、Xn(t) indicating that if there is vibration within the set sampling period, the vibration is transmitted from one of the mating ends. Then, the cross-correlation denoising method is adopted to X1(t)、X2(t)……Xn-1(t)、Xn(t) denoising and extracting vibration component signal S1(t)，S2(t-t1,2)……Sn-1(t-t1,n-1)，Sn(t-t1,n) The relationship between the vibration signal and the corresponding vibration component signal is:
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)
wherein Z is1(t)、Z2(t)……Zn-1(t)、Zn(t) is an external random interference signal, λ1、λ2……λn-2、λn-1As attenuation factor, t1,2、t1,3……t1,n-1、t1,nIs the time delay difference of the vibration signal arriving at the vibration sensor 1 and the rest of the sensors.
It should be noted that, in step a, the embodiment performs denoising processing on the vibration signal in a cross-correlation manner to extract a medium vibration component signal, that is, an effective vibration signal, and a corresponding specific process belongs to the prior art and is not described herein again. Of course, as another embodiment, the vibration signal may be filtered by a method such as median filtering or mean filtering in the related art.
And B: obtaining the amplitude V of the vibration component signal according to the extracted vibration component signal1、V2……Vn-1、VnFor amplitude V of vibration component signal1、V2……Vn-1、VnAnd the time t corresponding to the vibration signal collected by the vibration sensor1、t2……tn-1、tnAnd (6) judging.
And C: if V1>V2>……>Vn-1>VnAnd t is1<t2<……<tn-1<tnAccording to the amplitude V of the vibration component signal1、V2And distance L2 between the vibration sensors 1 and 2, and calculating attenuation rate beta of vibration signal per unit distance1The corresponding calculation formula is:
in addition, the attenuation coefficient β per unit distance of the vibration signal is calculated1Using the amplitude V of the vibration component signal1、V2The distance L2 between the vibration sensor 1 and the vibration sensor 2 is set to increase the damping rate β1The accuracy of (2). As another embodiment, the amplitude V of the vibration component signal may be used1、V2……Vn-1、VnThe amplitude of any two vibration component signals and the amplitude between corresponding two vibration sensorsAnd calculating the attenuation rate of the vibration signal per unit distance.
Step D: according to the amplitude V of the vibration component signal1Damping rate beta of vibration signal per unit distance1And the distance L1 between the vibration sensor 1 and its nearest mating end, the equivalent excitation signal V of which is calculated01The calculation formula is as follows:
it should be noted that the equivalent excitation signal V is calculated01Using the amplitude V of the vibration component signal1And the distance L1 between the vibration sensor 1 and its nearest mating end is to increase the equivalent excitation signal V01The accuracy of (2). As another example, the amplitude V of the vibration component signal may be selected1、V2……Vn-1、VnCalculating an equivalent excitation signal according to the amplitude of any vibration component signal and the distance between the vibration sensor corresponding to the amplitude of the selected vibration component signal and the target matching end, wherein the target matching end is the closest matching end of the vibration sensor corresponding to the amplitude of the vibration component signal with the largest distance.
In addition, in case 1, if V1<V2<……<Vn-1<VnAnd t is1>t2>……>tn-1>tnAccording to the amplitude V of the vibration component signaln-1、VnAnd the distance Ln between the vibration sensor n-1 and the vibration sensor n, and calculating the attenuation rate beta of the vibration signal per unit distance2The corresponding calculation formula is:
according to the amplitude V of the vibration component signalnDamping rate beta of vibration signal per unit distance2And the distance Ln + between the vibration sensor n and its nearest mating end1, calculating the equivalent excitation signal V of the nearest matching end02The calculation formula is as follows:
in case 1, if V is not satisfied1>V2>……>Vn-1>VnAnd t is1<t2<……<tn-1<tnNor satisfy V1<V2<……<Vn-1<VnAnd t is1>t2>……>tn-1>tnIf the interference signal with the same frequency as the vibration signal is mixed in the sampling period, the data is discarded.
Case 2:
step a: if 2 vibration signals of each vibration sensor are acquired in each set sampling period, namely a group of vibration signals acquired by the vibration sensors 1, 2 … … n-1 and n are respectively X1,1(t)、X2,1(t)……Xn-1,1(t)、Xn,1(t), the other group of vibration signals collected by the vibration sensors 1, 2 … … n-1 and n are respectively X1,2(t)、X2,2(t)……Xn-1,2(t)、Xn,2(t) if there is vibration in the set sampling period, the vibration is transmitted from two matching ends respectively, and the X is subjected to cross-correlation denoising respectively1,1(t)、X2,1(t)……Xn-1,1(t)、Xn,1(t) and X1,2(t)、X2,2(t)……Xn-1,2(t)、Xn,2(t) denoising and extracting vibration component signal S1,1(t)、S2,1(t-t1,2)……Sn,1(t-t1,n) And S1,2(t)、S2,2(t+t1,2)……Sn,2(t+t1,n)，t1,2、t1,3……t1,n-1、t1,nIs the time delay difference of the vibration signal arriving at the vibration sensor 1 and the rest of the sensors.
In step a, the embodiment also performs denoising processing on the vibration signal in a cross-correlation manner to extract a medium vibration component signal. As another embodiment, the vibration signal may be filtered by a method such as median filtering or mean filtering in the related art.
Step b: obtaining the amplitude V of the vibration component signal according to the extracted vibration component signal1,1、V2,1……Vn-1,1、Vn,1And V1,2、V2,2……Vn-1,2、Vn,2For amplitude V of vibration component signal1,1、V2,1……Vn-1,1、Vn,1And V1,2、V2,2……Vn-1,2、Vn,2And the time t corresponding to the vibration signal collected by the vibration sensor1,1、t2,1……tn-1,1、tn,1And t1,2、t2,2……tn-1,2、tn,2And (6) judging.
Step c: if V1,1>V2,1>……>Vn-1,1>Vn,1，t1,1<t2,1<……<tn-1,1<tn,1And V is1,2<V2,2<……<Vn-1,2<Vn,2，t1,2>t2,2>……>tn-1,2>tn,2According to the amplitude V of the vibration component signal1,1、V2,1And distance L2 between the vibration sensors 1 and 2, and calculating attenuation rate beta of vibration signal per unit distance1According to the amplitude V of the vibration component signaln,2、Vn-1,2And the distance Ln between the vibration sensor n-1 and the vibration sensor n, and calculating the attenuation rate beta of the vibration signal per unit distance2。
Wherein, similarly to step C in case 1, the attenuation rate β per unit distance of the vibration signal is calculated1Attenuation rate beta of unit distance from vibration signal2In other embodiments, any two amplitudes corresponding to the amplitudes of all the vibration component signals and the distance between the vibration sensors corresponding to the two amplitudes may be used, and details are not described here.
Step d: according to the amplitude V of the vibration component signal1,1Damping rate beta of vibration signal per unit distance1And the distance L1 between the vibration sensor 1 and its nearest mating end, the equivalent excitation signal V of which is calculated01(ii) a According to the amplitude V of the vibration component signaln,2Damping rate beta of vibration signal per unit distance2And the distance Ln +1 between the vibration sensor n and the nearest mating end thereof, calculating the equivalent excitation signal V of the nearest mating end02The corresponding calculation formula can be referred to case 1.
Wherein, similarly to step D in case 1, the equivalent excitation signal V is calculated01And an equivalent excitation signal V02In other embodiments, any amplitude corresponding to the amplitudes of all vibration component signals and the distance between the vibration sensor and the target matching end corresponding to the amplitude may also be used, and details are not repeated here.
In addition, in case 2, if V1,1<V2,1<……<Vn-1,1<Vn,1，t1,1>t2,1>……>tn-1,1>tn,1And V is1,2>V2,2>……>Vn-1,2>Vn,2，t1,2<t2,2<……<tn-1,2<tn,2According to the amplitude V of the vibration component signaln,1、Vn-1,1And the distance Ln between the vibration sensor n-1 and the vibration sensor n, and calculating the attenuation rate beta of the vibration signal per unit distance2According to the amplitude V of the vibration component signal1,2、V2,2And distance L2 between the vibration sensors 1 and 2, and calculating attenuation rate beta of vibration signal per unit distance1. According to the amplitude V of the vibration component signaln,1Damping rate beta of vibration signal per unit distance2And the distance Ln +1 between the vibration sensor n and the nearest mating end thereof, calculating the equivalent excitation signal V of the nearest mating end02. According to the amplitude V of the vibration component signal1,2Damping rate beta of vibration signal per unit distance1And the distance L1 between the vibration sensor 1 and its nearest mating end, which is calculatedTerminal equivalent excitation signal V01The corresponding calculation formula can be referred to case 1.
In case 2, if V is not satisfied1,1>V2,1>……>Vn-1,1>Vn,1，t1,1<t2,1<……<tn-1,1<tn,1And V is1,2<V2,2<……<Vn-1,2<Vn,2，t1,2>t2,2>……>tn-1,2>tn,2Nor satisfy V1,1<V2,1<……<Vn-1,1<Vn,1，t1,1>t2,1>……>tn-1,1>tn,1And V is1,2>V2,2>……>Vn-1,2>Vn,2，t1,2<t2,2<……<tn-1,2<tn,2If the set sampling period is mixed with an interference signal having the same frequency as the vibration signal, the data is discarded.
Case 3: if the condition that 1 vibration signal of each vibration sensor is collected in each set sampling period is not met, and 2 vibration signals of each vibration sensor are not collected in each set sampling period, the condition that a pure interference signal exists in the sampling period is indicated, and data are discarded.
In order to make the content in the above steps 1) and 2) more clear, the vibration monitoring portion with 3 vibration sensors disposed in fig. 3 is taken as an example for explanation, and the flow chart of the corresponding steps is shown in fig. 4:
the vibration signal that gathers vibration sensor is resampled according to setting for the sampling cycle, and the vibration signal number of every vibration sensor who gathers in every setting for the sampling cycle is divided into three kinds of condition:
case 1: if a signal X of three vibration sensors is obtained in a set sampling period1(t)、X2(t)、X3(t) it means that if there is vibration within this set sampling period, the vibration is transmitted from the mating end 1 or the mating end 2. Using cross-correlation denoising method to X1(t)、X2(t)、X3(t) extracting vibration component signal S after denoising treatment1(t) and S2(t-t1,2)、S3(t-t1,3) The calculation formula is as follows:
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)
wherein Z is1(t)、Z2(t)、Z3(t) is an external random interference signal, λ1For damping the vibration signal from the vibration sensor 1 to the vibration sensor 2, lambda2Damping factor, t, for vibration signals from vibration sensor 1 to vibration sensor 31,2For time delay differences of the arrival of the vibration signals at the vibration sensor 1 and the vibration sensor 2, t1,3Is the time delay difference of the vibration signal arriving at the vibration sensor 1 and the vibration sensor 3. Lambda [ alpha ]1、λ2May be a number greater or less than 1; t is t1,2、t1,3Can take positive and negative values.
Obtaining the amplitude V of the vibration component signal according to the extracted vibration component signal1、V2、V3Determining the amplitude V of the vibration component signal1、V2、V3And the time t when the vibration sensor acquires the vibration signal1、t2、t3The relationship of (1):
if V1>V2>V3And t is1<t2<t3According to the amplitude V of the vibration component signal1、V2Calculating the attenuation rate beta of the vibration signal per unit distance from the distance L2 of the vibration sensors 1 and 21. According to the amplitude V of the vibration component signal1Damping rate beta of vibration signal per unit distance1And the distance L1 between the vibration sensor 1 and the mating end 1, calculating the equivalent excitation signal V of the mating end 101The calculation formula is as follows:
if V1<V2<V3And t is1>t2>t3According to the amplitude V of the vibration component signal2、V3Calculating the attenuation rate beta of the vibration signal per unit distance from the distance L3 of the vibration sensors 2 and 32. According to the amplitude V of the vibration component signal3Damping rate beta of vibration signal per unit distance2And the distance L4 between the vibration sensor 3 and the mating end 2, calculating the equivalent excitation signal V of the mating end 202The calculation formula is as follows:
if V is not satisfied1>V2>V3And t is1<t2<t3Nor satisfy V1<V2<V3And t is1>t2>t3If so, it means that the interference signal with the same frequency as the vibration signal is mixed in the set sampling period, and the data is discarded.
Case 2: if two signals X of three vibration sensors are respectively obtained in a set sampling period1,1(t)、X2,1(t)、X3,1(t) and X1,2(t)、X2,2(t)、X3,2(t), it means that if there is vibration in this set sampling period, the vibration is transmitted from both mating ends, i.e., the mating end 1 and the mating end 2, respectively. Using cross-correlation denoising method to X1,1(t)、X2,1(t)、X3,1(t) carrying out denoising treatment to obtain a vibration component signal S1,1(t)、S2,1(t-t1,2)、S3,1(t-t1,3) Using cross-correlation denoising method to X1,2(t)、X2,2(t)、X3,2(t) carrying out denoising treatment to obtain a vibration component signal S1,2(t)、S2,2(t+t1,2)、S3,2(t+t1,3). Wherein, t1,2For time delay differences of the arrival of the vibration signals at the vibration sensor 1 and the vibration sensor 2, t1,3Is the time delay difference of the vibration signal arriving at the vibration sensor 1 and the vibration sensor 3.
Obtaining the amplitude V of the vibration component signal according to the extracted vibration component signal1,1、V2,1、V3,1And V1,2、V2,2、V3,2Determining the amplitude V of the vibration component signal1,1、V2,1、V3,1And V1,2、V2,2、V3,2And the time t when the vibration sensor acquires the vibration signal1,1、t2,1、t3,1And t1,2、t2,2、t3,2The relationship of (1):
if V1,1>V2,1>V3,1，t1,1<t2,1<t3,1And V is1,2<V2,2<V3,2，t1,2>t2,2>t3,2According to the amplitude V of the vibration component signal1,1、V2,1Calculating the attenuation rate beta of the vibration signal per unit distance from the distance L2 of the vibration sensors 1 and 21According to the amplitude V of the vibration component signal3,2、V2,2Calculating the attenuation rate beta of the vibration signal per unit distance from the distance L3 of the vibration sensors 2 and 32. According to the amplitude V of the vibration component signal1,1Damping rate beta of vibration signal per unit distance1And the distance L1 between the vibration sensor 1 and the mating end 1, calculating the equivalent excitation signal V of the mating end 101(ii) a According to the amplitude V of the vibration component signal3,2Damping rate beta of vibration signal per unit distance2And the distance L4 between the vibration sensor 3 and the mating end 2, calculating the equivalent excitation signal V of the mating end 202。
If V1,1<V2,1<V3,1，t1,1>t2,1>t3,1And V is1,2>V2,2>V3,2，t1,2<t2,2<t3,2According to the amplitude V of the vibration component signal3,1、V2,1Calculating the attenuation rate beta of the vibration signal per unit distance from the distance L3 of the vibration sensors 2 and 32According to the amplitude V of the vibration component signal1,2、V2,2Calculating the attenuation rate beta of the vibration signal per unit distance from the distance L2 of the vibration sensors 1 and 21. According to the amplitude V of the vibration component signal3,1Damping rate beta of vibration signal per unit distance2And the distance L4 between the vibration sensor 3 and the mating end 2, calculating the equivalent excitation signal V of the mating end 202(ii) a According to the amplitude V of the vibration component signal1,2Damping rate beta of vibration signal per unit distance1And the distance L1 between the vibration sensor 1 and the mating end 1, calculating the equivalent excitation signal V of the mating end 101。
(3) And applying the equivalent excitation signal at the matching end of the vibration monitoring part to the corresponding position of the finite element simulation model, and simulating to obtain a strain image of the vibration monitoring part and strain data of each finite element point.
The finite element simulation unit U3 establishes a finite element model of the fan vibration monitoring part and the matching parts in finite element simulation software according to the actual condition of the fan. There are many existing finite element simulation software, and in this embodiment, the adopted finite element simulation software is comsolmutithphysics. In addition, the specific process of the finite element simulation unit U3 establishing a finite element model of the fan vibration monitoring part and the mating component in the finite element simulation software belongs to the prior art, and is not described herein again.
And then applying an equivalent excitation signal of the matching end, which is obtained by calculation through a vibration excitation equivalent method and transmitted by the data processing unit U2, to a corresponding position of the finite element model to perform strain simulation. And packaging and compressing the simulated strain image of the vibration monitoring part, strain data of each finite element point and an equivalent excitation signal at the matching end of the vibration monitoring part transmitted by the data processing unit U2, and transmitting the equivalent excitation signal to the monitoring analysis unit U4 in real time. For example, fig. 5 shows a strain image of a vibration monitoring region obtained by a certain simulation, the strain image of the vibration monitoring region includes a left-side graph and a right-side bar graph of fig. 5, and the bar graphs are illustrations of colors representing strain magnitudes and correspond to legends.
(4) And displaying the strain image of the vibration monitoring part, and analyzing the equivalent excitation signal at the matching end of the vibration monitoring part and the strain data of the finite element point to obtain the corresponding health state of the vibration monitoring part.
Specifically, the monitoring and analyzing unit U4 decompresses the packed signal transmitted by the finite element simulation unit U3, displays the strain image therein, and processes and analyzes the equivalent excitation signal at the matching end and the strain data at the finite element point to diagnose the mechanical health state caused by the vibration of the vibration monitoring portion. The method comprises the following steps of analyzing equivalent excitation signals at the matching end of the vibration monitoring part and strain data of finite element points, wherein the steps comprise:
4.1) setting a threshold value, and screening out equivalent excitation signals of the matching end of the vibration monitoring part exceeding the set threshold value.
The corresponding set thresholds of different monitoring parts of the fan are different, and need to be determined according to field tests, experience and needs.
4.2) judging whether the equivalent excitation signal of the matching end of the screened vibration monitoring part is a periodic signal, wherein the method is divided into three situations:
in the first case, if two equivalent excitation signals are selected at the coupling end of the screened vibration monitoring part, and the equivalent excitation signals at the two coupling ends are both periodic signals, it is said that the part to which the vibration monitoring part belongs has a fault, and it is proposed to inspect the part.
The second case: if one of the equivalent excitation signals of the matched end of the screened vibration monitoring part is a periodic signal (comprising that only one equivalent excitation signal of the matched end of the screened vibration monitoring part is a periodic signal, and only one equivalent excitation signal of the matched end of the screened vibration monitoring part is a periodic signal), judging whether the equivalent excitation signals of the two matched ends of the component directly connected with the matched end to which the equivalent excitation signal of the periodic matched end belongs are periodic signals, and if so, judging that the component directly connected with the matched end to which the equivalent excitation signal of the periodic matched end belongs has a fault; otherwise, the matching region of the matching end of the equivalent excitation signal of the periodic matching end and the part directly connected with the matching end has a fault.
The third situation: if the equivalent excitation signals of the matching end of the screened vibration monitoring part are not periodic signals (including that only one equivalent excitation signal of the matching end of the screened vibration monitoring part is not periodic signal, and two equivalent excitation signals of the matching end of the screened vibration monitoring part are not periodic signals), the simulation strain magnitude epsilon of each finite element point of the vibration monitoring part is determinedImitation iCan bear rated strain epsilon with each point of the corresponding vibration monitoring partForehead iAnd (3) comparison: if epsilonImitation i>εForehead iIf so, the health state of the vibration monitoring part caused by the vibration is indicated as an alarm state; if epsilonImitation i>η*εForehead i，0<η<1, if eta is 0.8, the health state of the vibration monitoring part caused by vibration is an early warning state; otherwise, the health state of the vibration monitoring part caused by the vibration is a normal state.
Wherein, the simulation strain magnitude epsilon of each finite element point of the vibration monitoring partImitation iThe strain data of the finite element points obtained by the simulation is contained in the compression packet transmitted by the finite element simulation unit U3.
For example, in fig. 5, the vibration monitoring portion B is shown, and the flow chart of the steps of the vibration state monitoring analysis is shown in fig. 6, at this time, if the equivalent excitation signals at the two mating ends of the screened vibration monitoring portion B are both periodic signals, that is, the excitation signals at the equivalent mating ends B1 and B2 are both periodic signals, it is indicated that the component to which the vibration monitoring portion B belongs has a fault, and it is recommended to check the component. If only one of the equivalent excitation signals of the two mating ends B1 and B2 of the vibration monitoring part B is a periodic signal, continuously judging whether the equivalent excitation signals of the two mating ends of the component A or C directly connected with the equivalent mating end B1 or B2 to which the periodic equivalent excitation signal belongs are periodic signals, if so, indicating that the directly connected component A or C has a fault, and recommending that the component is checked; if not, the matching area, namely the matching area of A and B or B and C has a fault, is suggested to be checked.
The fan vibration state detection method is based on the collected signals of the vibration sensor, combines finite element analysis, fully considers different bearing capacities of all parts of the detection component, can visually display the overall state and state change of the monitoring component, can early warn the health state of the monitoring component, and overcomes the defects of incomplete monitoring component, state evaluation hysteresis and non-visual display. By adopting the vibration excitation equivalent method, the vibration source does not need to be positioned, the algorithm is simple, the method is particularly effective in the vibration influence analysis of the fan device with complex and matched components, the influence of environmental interference signals is removed, and the monitoring analysis is more accurate.
In addition, it should be noted that the fan vibration state detection method can be applied not only to the fan fault detection process, but also to the fault detection processes of other devices.
Fan vibration state detection device embodiment:
the embodiment provides a fan vibration state detection device, which includes a processor and a memory, where the processor is configured to process an instruction stored in the memory to implement a fan vibration state detection method, and the fan vibration state detection method is described in detail in the above embodiment of the fan vibration state detection method, and is not described here again.
The embodiment of the equivalent excitation determining method of the matching end of the fan vibration state comprises the following steps:
the embodiment provides an equivalent excitation determining method for a matching end of a fan vibration state, which is used for determining equivalent excitation signals of two matching ends of a fan vibration monitoring part. Of course, the method can be suitable for determining equivalent excitation signals of the matching ends of the fan vibration monitoring parts, and can also be suitable for determining equivalent excitation signals of the matching ends of the vibration monitoring parts of other equipment. In addition, the calculated equivalent excitation signal of the matching end can be applied to the fault detection process of the monitoring component and also has other purposes. Since the method for determining the equivalent excitation at the matching end of the fan vibration state has been described in detail in the embodiment of the fan vibration state detection method, the details are not described herein.
Fan vibration state detection system embodiment:
the embodiment provides a fan vibration state detection system which comprises a processor, a memory and a vibration sensor arranged at a fan vibration monitoring part, wherein the processor is connected with the memory and the vibration sensor and is used for realizing a fan vibration state detection method. Since the working process of the fan vibration state detection system has been described in detail in the above embodiment of the fan vibration state detection method, further description is omitted here.
Finally, it should be noted that the above embodiments are only used for illustrating the technical solutions of the present invention and not for limiting the protection scope thereof, and although the present application is described in detail with reference to the above embodiments, those skilled in the art should understand that after reading the present application, various changes, modifications or equivalents of the embodiments of the present application can be made, and these changes, modifications or equivalents are within the protection scope of the claims of the present invention.
Claims (7)
1. A method for detecting the vibration state of a fan is characterized by comprising the following steps:
(1) resampling the vibration signals collected by the n vibration sensors at the vibration monitoring part of the fan according to a set sampling period;
(2) calculating a vibration component signal according to the obtained vibration signal of the vibration monitoring part of the fan;
(3) calculating an equivalent excitation signal at the matching end of the vibration monitoring part of the fan according to the amplitude of the vibration component signal and the time for the vibration sensor to acquire the vibration signal;
(4) determining a fault area of the fan vibration monitoring part according to the periodicity of the equivalent excitation signal at the matching end of the fan vibration monitoring part;
the step of calculating the equivalent excitation signal of the matching end of the fan vibration monitoring part in the step (3) comprises the following steps:
if the n vibration sensors acquire 1 vibration signal X in each set sampling period1(t)、X2(t)……Xn-1(t)、Xn(t) and satisfy V1>V2>……>Vn-1>VnAnd t is1<t2<……<tn-1<tnOr satisfy V1<V2<……<Vn-1<VnAnd t is1>t2>……>tn-1>tnWherein V isiAmplitude, t, of vibration component signal calculated for vibration signal collected by ith vibration sensoriCalculating the attenuation rate of the unit distance of the vibration signal according to the amplitude of two vibration component signals and the distance between the two corresponding vibration sensors when the ith vibration sensor acquires the vibration signal;
and calculating the equivalent excitation signal of the target matching end according to the selected amplitude of one vibration component signal, the attenuation rate of the unit distance of the vibration signal and the distance between the vibration sensor corresponding to the selected amplitude of one vibration component signal and the target matching end, wherein the target matching end is the closest matching end of the vibration sensor corresponding to the amplitude of the vibration component signal with the largest distance.
2. The fan vibration state detection method according to claim 1, wherein the step of calculating the equivalent excitation signal at the coupling end of the fan vibration monitoring part comprises:
if the n vibration sensors acquire 2 vibration signals X in each set sampling period1,1(t)、X2,1(t)……Xn-1,1(t)、Xn,1(t) and X1,2(t)、X2,2(t)……Xn-1,2(t)、Xn,2(t) and satisfy V1,1>V2,1>……>Vn-1,1>Vn,1，t1,1<t2,1<……<tn-1,1<tn,1And V is1,2<V2,2<……<Vn-1,2<Vn,2，t1,2>t2,2>……>tn-1,2>tn,2Or satisfy V1,1<V2,1<……<Vn-1,1<Vn,1，t1,1>t2,1>……>tn-1,1>tn,1And V is1,2>V2,2>……>Vn-1,2>Vn,2，t1,2<t2,2<……<tn-1,2<tn,2Wherein V isi,1Amplitude, t, of vibration component signal calculated for first vibration signal collected by ith vibration sensori,1Time, V, of acquisition of first vibration signal for ith vibration sensori,2Amplitude, t, of vibration component signal calculated for second vibration signal collected by ith vibration sensori,2The time when the second vibration signal is acquired by the ith vibration sensor is determined according to the amplitude V of each vibration component signal1,1、V2,1……Vn-1,1、Vn,1Calculates the attenuation rate of the first vibration signal per unit distance based on the amplitudes V of the vibration component signals and the distances between the two vibration sensors1,2、V2,2……Vn-1,2、Vn,2Calculating the attenuation rate of the unit distance of the second vibration signal according to the amplitude of the two vibration component signals and the distance between the two corresponding vibration sensors;
according to the amplitude V of each vibration component signal1,1、V2,1……Vn-1,1、Vn,1The amplitude of the selected one vibration component signal, the attenuation rate of the first vibration signal per unit distance, and the amplitude V of each vibration component signal1,1、V2,1……Vn-1,1、Vn,1Calculating the equivalent excitation signal of the first target mating end corresponding to the distance between the vibration sensor and the first target mating end corresponding to the amplitude of the selected vibration component signal, wherein the first target mating end is the amplitude V of each vibration component signal1,1、V2,1……Vn-1,1、Vn,1Of the largest one of the vibration component signalsThe closest matching end of the vibration sensor corresponding to the amplitude; according to the amplitude V of each vibration component signal1,2、V2,2……Vn-1,2、Vn,2The amplitude of the selected one vibration component signal, the attenuation rate of the second vibration signal per unit distance, and the amplitude V of each vibration component signal1,2、V2,2……Vn-1,2、Vn,2Calculating the equivalent excitation signal of the second target matching end by the distance between the vibration sensor corresponding to the amplitude of the selected vibration component signal and the second target matching end, wherein the second target matching end is the amplitude V of each vibration component signal1,2、V2,2……Vn-1,2、Vn,2The amplitude of the largest one of the vibration component signals corresponds to the closest mating end of the vibration sensor.
3. The fan vibration state detection method according to claim 2, wherein if the amplitudes of the equivalent excitation signals at the two mating ends of the vibration monitoring portion are both greater than a set threshold value, and the equivalent excitation signals at the two mating ends are both periodic signals, it is determined that a component to which the vibration monitoring portion belongs has a fault.
4. The fan vibration state detection method according to claim 1 or 2, characterized in that if the amplitude of the equivalent excitation signal of at least one of the mating ends of the vibration monitoring portion is greater than a set threshold, and only the equivalent excitation signal of the first mating end is a periodic signal among the equivalent excitation signals greater than the set threshold, it is determined whether two equivalent excitation signals of a component directly connected to the first mating end are satisfied as periodic signals; if so, a component directly connected to the first mating end fails, otherwise, the first mating end and the mating area of the component directly connected thereto fail.
5. The fan vibration state detection method according to claim 1 or 2, further comprising: applying the equivalent excitation signal of the matching end of the vibration monitoring part to the corresponding position of the finite element simulation model, and simulating to obtain a strain image of the vibration monitoring part; and if the amplitude of the equivalent excitation signal at the matching end of the vibration monitoring part is greater than a set threshold value and the equivalent excitation signals greater than the set threshold value are not periodic signals, determining the health state of the vibration monitoring part according to the relationship between the simulation strain of each finite element point of the vibration monitoring part and the rated strain bearable of each point of the vibration monitoring part obtained by simulation.
6. The fan vibration state detection method according to any one of claims 1 to 2, wherein a propagation time of the vibration signal between two vibration sensors most distant from the vibration monitoring portion < the set sampling period < a propagation time of the vibration signal between two vibration sensors most distant from the vibration monitoring portion + a minimum value of a propagation time of the vibration signal between different mating ends of the fan vibration monitoring portion and a vibration sensor closest to the mating end.
7. A fan vibration state detection device, which is characterized by comprising a processor, a memory and a vibration sensor arranged at a fan vibration monitoring part, wherein the processor is used for processing instructions stored in the memory to realize the fan vibration state detection method according to any one of claims 1 to 6; the processor is connected with the memory and the vibration sensor and is used for: receiving vibration signals which are sent by n vibration sensors and collected at a fan vibration monitoring part according to a set sampling period, and storing the vibration signals in a memory; calculating a vibration component signal according to the vibration signal stored in the memory; calculating an equivalent excitation signal at the matching end of the vibration monitoring part of the fan according to the amplitude of the vibration component signal and the time for the vibration sensor to acquire the vibration signal; and determining the fault area of the fan vibration monitoring part according to the periodicity of the equivalent excitation signal at the matching end of the fan vibration monitoring part.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911067941.6A CN110793615B (en) | 2019-11-04 | 2019-11-04 | Fan vibration state detection method and device and matching end equivalent excitation determination method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911067941.6A CN110793615B (en) | 2019-11-04 | 2019-11-04 | Fan vibration state detection method and device and matching end equivalent excitation determination method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110793615A CN110793615A (en) | 2020-02-14 |
| CN110793615B true CN110793615B (en) | 2021-09-14 |
Family
ID=69442650
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911067941.6A Active CN110793615B (en) | 2019-11-04 | 2019-11-04 | Fan vibration state detection method and device and matching end equivalent excitation determination method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110793615B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101545459A (en) * | 2008-03-25 | 2009-09-30 | 唐德尧 | Device for monitoring running state of wind generator tower |
| CN103049670A (en) * | 2013-01-14 | 2013-04-17 | 哈尔滨工程大学 | Pipe excitation source identification and prediction method of vibration response to pipe excitation source |
| CN104048825A (en) * | 2014-05-28 | 2014-09-17 | 北京航空航天大学 | Multi-sensor fusion type fault locating method for gearbox of wind turbine generator |
| US8990031B2 (en) * | 2009-12-10 | 2015-03-24 | Industrial Technology Research Institute | Torsional resonance frequency measuring device and method |
| CN104484827A (en) * | 2014-09-29 | 2015-04-01 | 许继集团有限公司 | Wind power generator set failure frequency component extraction method |
| CN104915498A (en) * | 2015-06-08 | 2015-09-16 | 广东工业大学 | Model identification and equivalent simplification based high-speed platform motion parameter self-adjusting method |
| CN110186679A (en) * | 2019-05-13 | 2019-08-30 | 温州大学 | The diagnostic method of shield main shaft bearing |
-
2019
- 2019-11-04 CN CN201911067941.6A patent/CN110793615B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101545459A (en) * | 2008-03-25 | 2009-09-30 | 唐德尧 | Device for monitoring running state of wind generator tower |
| US8990031B2 (en) * | 2009-12-10 | 2015-03-24 | Industrial Technology Research Institute | Torsional resonance frequency measuring device and method |
| CN103049670A (en) * | 2013-01-14 | 2013-04-17 | 哈尔滨工程大学 | Pipe excitation source identification and prediction method of vibration response to pipe excitation source |
| CN104048825A (en) * | 2014-05-28 | 2014-09-17 | 北京航空航天大学 | Multi-sensor fusion type fault locating method for gearbox of wind turbine generator |
| CN104484827A (en) * | 2014-09-29 | 2015-04-01 | 许继集团有限公司 | Wind power generator set failure frequency component extraction method |
| CN104915498A (en) * | 2015-06-08 | 2015-09-16 | 广东工业大学 | Model identification and equivalent simplification based high-speed platform motion parameter self-adjusting method |
| CN110186679A (en) * | 2019-05-13 | 2019-08-30 | 温州大学 | The diagnostic method of shield main shaft bearing |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110793615A (en) | 2020-02-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3993825B2 (en) | Inference signal generator for instrumented equipment and processes | |
| CN104048825B (en) | A kind of gearbox of wind turbine Fault Locating Method of Multi-sensor Fusion | |
| CN103471841B (en) | A kind of rotating machinery vibrating failure diagnosis method | |
| EP2476901A1 (en) | Method and monitoring apparatus for automated surveillance of a wind turbine and a method for creating a linear model | |
| Delvecchio et al. | Advanced signal processing tools for the vibratory surveillance of assembly faults in diesel engine cold tests | |
| CN112834193B (en) | Operation bridge vibration and health state abnormity early warning method based on three-dimensional graph | |
| CN108444662A (en) | A kind of bridge damnification on-line monitoring method based on degree/day effect | |
| CN102901560A (en) | Safe comprehensive monitoring system for structure of offshore jacket platform | |
| 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 (en) | Method for condition monitoring of rotating machines via wavelet thresholding function and 4th-order moment | |
| EP2345894A2 (en) | Trending of vibration data taking into account torque effect | |
| EP3480455B1 (en) | Wind turbine monitoring device, wind turbine monitoring method, wind turbine monitoring program, and storage medium | |
| CN113467433A (en) | Method and device for detecting faults of mechanical equipment | |
| CN113569445A (en) | Steel structure health monitoring system and method based on digital twinning technology | |
| JP2012018045A (en) | Sensor failure diagnosis device and sensor system | |
| CN109470972B (en) | Transformer core looseness diagnosis method based on vibration image | |
| CN110990978A (en) | Bolt state monitoring method and device | |
| CN111311872A (en) | Long-term monitoring and alarming system for stress of hull structure | |
| CN102072926B (en) | Method for diagnosing body fatigue crack of motor | |
| KR101046748B1 (en) | Dedicated Vibration Diagnosis Method and Analysis System for Generator Stator Windings | |
| CN110793615B (en) | Fan vibration state detection method and device and matching end equivalent excitation determination method | |
| CN110598680A (en) | Method and system for evaluating health state of mechanical equipment and readable storage medium | |
| CN110411730A (en) | A kind of rotating equipment failures judgment method, system and readable storage medium storing program for executing | |
| EP1686443A1 (en) | Methods, systems, and computer program products for implementing condition monitoring activities | |
| CN114715354A (en) | Health management device for marine equipment and equipment state detection and fault diagnosis method |
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 |