CN108709635B - Method and device for determining main frequency component information of rotary mechanical vibration signal - Google Patents
Method and device for determining main frequency component information of rotary mechanical vibration signal Download PDFInfo
- Publication number
- CN108709635B CN108709635B CN201810938946.0A CN201810938946A CN108709635B CN 108709635 B CN108709635 B CN 108709635B CN 201810938946 A CN201810938946 A CN 201810938946A CN 108709635 B CN108709635 B CN 108709635B
- Authority
- CN
- China
- Prior art keywords
- frequency component
- determining
- amplitude spectrum
- main frequency
- component information
- 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
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000001228 spectrum Methods 0.000 claims abstract description 108
- 238000005070 sampling Methods 0.000 claims abstract description 45
- 238000012545 processing Methods 0.000 claims abstract description 22
- 238000004590 computer program Methods 0.000 claims description 16
- 230000003595 spectral effect Effects 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 abstract description 10
- 238000010586 diagram Methods 0.000 description 13
- 230000006870 function Effects 0.000 description 12
- 238000010183 spectrum analysis Methods 0.000 description 5
- 210000001072 colon Anatomy 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention provides a method and a device for determining main frequency component information of a vibration signal of a rotating machine, wherein the method comprises the following steps: acquiring a rotating mechanical vibration time domain waveform sampling signal with preset duration and sampling frequency; carrying out fast Fourier transform on the time domain waveform sampling signal to obtain a complex sequence, and obtaining an amplitude spectrum according to the complex sequence; determining the main frequency component information of the vibration signal of the rotating machinery according to the preset frequency distinguishing interval between the main frequency components and the amplitude spectrum; the main frequency component information includes: the frequency and amplitude of the dominant frequency component. According to the technical scheme, the accuracy of determining the main frequency component information of the vibration signal of the rotary machine is improved, fault reasons and fault location analysis can be performed more accurately, technical support is provided for timely processing faults and recovering equipment operation as early as possible, and efficient operation of the equipment is guaranteed.
Description
Technical Field
The invention relates to the technical field of data processing, in particular to a method and a device for determining main frequency component information of a vibration signal of a rotating machine.
Background
A rotary machine is a machine that mainly depends on a rotary motion to perform a specific function. Typical rotating machines are steam turbines, gas turbines, centrifugal and axial compressors, fans, pumps, water turbines, electrical generators, aircraft engines, and the like. As soon as the rotary machine starts to operate, vibration is inevitably generated. If the vibration is too large, great harm is caused, and the working performance of the equipment is reduced or the equipment fails to work. In addition, vibration can cause accelerated wear, fatigue, cracking, etc. of certain components under additional loading, which can affect life or cause equipment failure. Therefore, it is necessary to test or monitor the vibration of the rotating machine. The vibration condition of the equipment can be quantitatively known through the vibration signal, and the reason of vibration generation can be known through frequency domain analysis of the vibration signal. Especially, when the equipment has faults, the vibration signals are tested and the time domain and the frequency domain are analyzed, so that the fault reason and the fault location can be obtained, and the technical guarantee is provided for timely processing the faults and recovering the operation as early as possible.
The spectral analysis is to represent the amplitude, phase or energy transformation of a signal by a frequency axis, and further analyze the frequency characteristics. After obtaining the time-domain waveform data of the vibration mode signal, it is necessary to perform analysis in the frequency domain to know the frequency-domain information of the vibration signal. This can be used to analyze the origin of the rotary machine vibration and cause analysis and fault location in the event of a fault.
The spectral analysis is implemented based on fast fourier, with amplitude, phase and power spectra. Amplitude spectra are mostly adopted in the vibration analysis of the rotating machinery. The amplitude spectrum is the frequency as independent variable and the amplitude of each frequency component composing the signal as dependent variable, and such frequency function is called amplitude spectrum, which characterizes the distribution of the amplitude of the signal along with the frequency as shown in fig. 1. The existing scheme has the problem of low accuracy in determining the main frequency component information of the vibration signal of the rotary machine.
In view of the above problems, no effective solution has been proposed.
Disclosure of Invention
The embodiment of the invention provides a method for determining main frequency component information of a vibration signal of a rotating machine, which is used for improving the accuracy of determining the main frequency component information of the vibration signal of the rotating machine and comprises the following steps:
acquiring a rotating mechanical vibration time domain waveform sampling signal with preset duration and sampling frequency;
carrying out fast Fourier transform on the time domain waveform sampling signal to obtain a complex sequence, and obtaining an amplitude spectrum according to the complex sequence;
determining the main frequency component information of the vibration signal of the rotating machinery according to the preset frequency distinguishing interval between the main frequency components and the amplitude spectrum; the main frequency component information includes: the frequency and amplitude of the dominant frequency component.
The embodiment of the invention also provides a device for determining the main frequency component information of the vibration signal of the rotating machinery, which is used for improving the accuracy of determining the main frequency component information of the vibration signal of the rotating machinery, and the device comprises:
the acquisition unit is used for acquiring a rotary mechanical vibration time domain waveform sampling signal with preset duration and sampling frequency;
the amplitude spectrum determining unit is used for carrying out fast Fourier transform on the time domain waveform sampling signal to obtain a complex sequence and obtaining an amplitude spectrum according to the complex sequence;
the main frequency component information determining unit is used for determining main frequency component information of the vibration signal of the rotating machinery according to preset frequency distinguishing intervals among main frequency components and the amplitude spectrum; the main frequency component information includes: the frequency and amplitude of the dominant frequency component.
The embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor realizes the method for determining the main frequency component information of the rotating mechanical vibration signal when executing the computer program.
Embodiments of the present invention also provide a computer-readable storage medium storing a computer program for executing the method for determining information of a main frequency component of a vibration signal of a rotating machine.
The technical scheme provided by the embodiment of the invention is as follows: firstly, acquiring a sampling signal and a sampling frequency of a rotary mechanical vibration time domain waveform with preset duration; secondly, performing fast Fourier transform on the time domain waveform sampling signal to obtain a complex sequence, and obtaining an amplitude spectrum according to the complex sequence; finally, determining the main frequency component information of the vibration signal of the rotating machinery according to the preset frequency distinguishing interval between the main frequency components and the amplitude spectrum; the main frequency component information includes: the frequency and amplitude of the primary frequency component; the method and the device have the advantages that the main frequency component information of the rotary mechanical vibration signal is determined according to the preset frequency distinguishing interval between the main frequency components, the repeated determination of the same main frequency component is avoided, the determined main frequency component corresponds to the same frequency component, and the second main frequency information and the third main frequency information which need to be determined are not included, all useful main frequency component information is accurately determined, so that the accuracy of determining the main frequency component information of the rotary mechanical vibration signal is improved, the fault reason and fault location analysis can be conveniently and accurately performed, the technical guarantee is provided for timely processing faults and recovering the equipment operation as early as possible, and the efficient operation of the equipment is ensured.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a schematic diagram of a frequency domain amplitude spectrum of a vibration signal according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a frequency domain amplitude spectrum of a vibration signal obtained in the prior art;
FIG. 3 is a flow chart illustrating a method for determining principal frequency component information of a vibration signal of a rotating machine according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of a time domain waveform of a vibration signal according to an embodiment of the present invention;
FIG. 5 is a frequency domain amplitude spectrogram (5 th order) of the main frequency component corresponding to FIG. 4 obtained by applying the method for determining the main frequency component information of the vibration signal of the rotating machine according to the embodiment of the present invention;
FIG. 6 is a frequency domain amplitude spectrogram (10 th order) of the main frequency component corresponding to FIG. 4 obtained by applying the method for determining the main frequency component information of the vibration signal of the rotating machine according to the embodiment of the present invention;
fig. 7 is a schematic structural diagram of an apparatus for determining information of main frequency components of a vibration signal of a rotating machine according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Before the embodiments of the present invention are described, the technical names related to the embodiments of the present invention are described first:
1. a rotating machine: a rotary machine is a machine that mainly depends on a rotary motion to perform a specific function. Typical rotating machines are steam turbines, gas turbines, centrifugal and axial compressors, fans, pumps, water turbines, electrical generators, aircraft engines, and the like.
2. Vibration: vibration is the reciprocating motion of an object (or a portion of an object) near an equilibrium position (the position of the object when stationary).
3. And (3) spectrum analysis: spectral analysis refers to a method of transforming a time domain signal into a frequency domain and analyzing the signal. The purpose of the spectrum analysis is to decompose a complex time history waveform into a plurality of single harmonic components through Fourier transform to research so as to obtain the frequency structure of a signal and information of each harmonic and phase.
The inventor finds that: in the prior art, a Fast Fourier Transform (FFT) is used to transform a vibration signal from a time domain to a frequency domain to obtain a complex sequence of the FFT, and further obtain an amplitude spectrum. And then sorting the amplitude spectrum from large to small, further calculating the frequency of the top 3-order main frequency component, and displaying the frequency on the amplitude spectrum. The disadvantages of the prior art are as follows:
(1) the phase magnitude corresponding to the dominant frequency component cannot be determined.
(2) The determination result cannot be automatically output to an Excel document or a txt document, and the main frequency component information also needs to be manually input in subsequent data processing and report compiling.
(3) It sometimes happens that the same main frequency component is repeatedly determined multiple times, thereby reducing useful determination information, as shown in fig. 2, the statistical 3 main frequency components are too little different, which all correspond to the same frequency component, while the information of the second main frequency and the third main frequency to be determined is not statistical (in fig. 2, f1 represents the first main frequency, 100.029 after the first colon is frequency size, 0.72 after the second colon represents amplitude size, the explanation about f2 line and f3 line is similar to that of f1, and in addition, the abscissa represents frequency and the ordinate represents amplitude in fig. 2).
In consideration of the technical problems, the inventor provides a scheme for determining the main frequency component information of the vibration signal of the rotary machine, which can automatically and accurately determine the main frequency component information (including frequency, amplitude and phase) of the vibration signal according to actual needs, and display and output the main frequency component information in an amplitude spectrogram to an Excel document or a txt document. The following describes a detailed description of the scheme for determining the information of the main frequency components of the vibration signal of the rotating machine.
Fig. 3 is a schematic flow chart of a method for determining information of main frequency components of a vibration signal of a rotating machine according to an embodiment of the present invention, as shown in fig. 3, the method includes the following steps:
step 101: acquiring a rotating mechanical vibration time domain waveform sampling signal with preset duration and sampling frequency;
step 102: carrying out fast Fourier transform on the time domain waveform sampling signal to obtain a complex sequence, and obtaining an amplitude spectrum according to the complex sequence;
step 103: determining the main frequency component information of the vibration signal of the rotating machinery according to the preset frequency distinguishing interval between the main frequency components and the amplitude spectrum; the main frequency component information includes: the frequency and amplitude of the dominant frequency component.
The technical scheme provided by the embodiment of the invention determines the main frequency component information of the rotary mechanical vibration signal according to the preset frequency distinguishing interval between the main frequency components, thereby avoiding the situation shown in figure 2: the same main frequency component is repeatedly determined, the determined main frequency component corresponds to the same frequency component, and the information of the second main frequency and the third main frequency which need to be determined is not included, so that all useful main frequency component information is accurately determined, the accuracy of determining the main frequency component information of the vibration signal of the rotating machine is improved, the fault reason and fault location analysis can be more accurately performed, the technical guarantee is provided for timely processing the fault and recovering the equipment operation as soon as possible, and the high-efficiency operation of the equipment is ensured.
The following describes each step of the embodiment of the present invention in detail.
In step 101, obtaining a sampling signal and a sampling frequency of a time-domain waveform of vibration of a rotating machine for a preset duration may include: sampling data (which can comprise time domain waveform sampling signals and sampling frequency) of integral multiple of the rotation period is obtained, so that the accuracy of determining the main frequency component information of the rotary mechanical vibration signal is improved.
In step 102, performing fast fourier transform on the time-domain waveform sampling signal to obtain a complex sequence, and obtaining an amplitude spectrum according to the complex sequence may include: and after the mean value of the time domain waveform sampling signal is subtracted, the fast Fourier transform is carried out to obtain a complex sequence, and then an amplitude spectrum is calculated, so that the accuracy of determining the main frequency component information of the rotary mechanical vibration signal is improved.
In the above step 103, the primary frequency discrimination interval is set in advance. In addition, besides setting the main frequency distinguishing interval, the method further includes a step of setting the order of the main frequency component, the number of the specific setting order may be flexibly set according to the actual working requirement, for example, a 5-order frequency domain amplitude spectrogram is shown in fig. 5, and in order to visually display information of more main frequency components, the statistical order of the main frequency component may be adjusted to 10-order, so as to obtain a 10-order frequency domain amplitude spectrogram as shown in fig. 6.
In step 103, determining the information of the main frequency components of the vibration signal of the rotating machine according to the preset frequency discrimination interval between the main frequency components and the amplitude spectrum may include:
determining the resolution of a frequency spectral line of the amplitude spectrum according to the array length and the sampling frequency corresponding to the time-domain waveform sampling signal;
determining a reference position interval in the amplitude spectrum array according to the resolution of the frequency spectral line of the amplitude spectrum and a preset frequency distinguishing interval between main frequency components;
and determining whether the frequency components to be confirmed are counted as the main frequency components according to the reference position interval in the amplitude spectrum array, the position serial number calculated in the amplitude spectrum array corresponding to the frequency components to be confirmed and the position serial number corresponding to the former confirmed main frequency component in the amplitude spectrum.
In an embodiment, determining whether the frequency component to be confirmed is counted as the main frequency component according to the reference position interval in the amplitude spectrum array, the position serial number calculated in the amplitude spectrum array corresponding to the frequency component to be confirmed, and the position serial number corresponding to the previously confirmed main frequency component in the amplitude spectrum may include:
a step of determining a first main frequency component; the step of determining the first dominant frequency component may comprise: acquiring the maximum amplitude and the corresponding frequency component in the amplitude spectrum array, and calculating the corresponding position serial number of the frequency component in the amplitude spectrum array;
a step of determining other main frequency components; the step of determining other major frequency components comprises: assigning the numerical values of the determined main frequency components on the corresponding position serial numbers in the amplitude spectrum array and all the position serial numbers in the interval of the left and right adjacent reference positions to be zero to obtain an updated amplitude spectrum array; acquiring the maximum amplitude and the corresponding frequency component in the updated amplitude spectrum array, calculating the corresponding position serial number of the frequency component in the updated amplitude spectrum array, and counting the frequency component as a main frequency component;
and continuing to perform the step of determining other primary frequency components based on the updated array of magnitude spectra until all primary frequency component information is found.
In specific implementation, according to the technical scheme provided in the foregoing embodiment, frequency distinguishing intervals between the main frequency components are preset, a reference position interval in the amplitude spectrum array is determined, and then, according to the reference position interval in the amplitude spectrum array and a position serial number calculated in the amplitude spectrum array corresponding to the frequency component to be confirmed (current frequency component), it is determined whether the frequency component to be confirmed is counted as the main frequency component, so that a situation shown in fig. 2 is avoided: the same main frequency component is repeatedly determined, the determined main frequency component corresponds to the same frequency component, and the information of the second main frequency and the third main frequency which need to be determined is not included, so that all useful main frequency component information is accurately determined, the accuracy of determining the main frequency component information of the vibration signal of the rotating machine is improved, the fault reason and fault location analysis can be more accurately performed, the technical guarantee is provided for timely processing the fault and recovering the equipment operation as soon as possible, and the high-efficiency operation of the equipment is ensured.
As the inventor finds that: in the existing scheme, sometimes an amplitude spectrogram of light is not enough, and the frequency, the amplitude and the corresponding phase size of main frequency components in the amplitude spectrogram need to be known specifically. However, the prior art cannot determine the phase magnitude corresponding to the dominant frequency component. Therefore, in consideration of this technical problem, the inventors propose the following embodiments.
In one embodiment, the primary frequency component information further includes: the phase magnitude of the primary frequency component;
the method for determining the main frequency component information of the rotary mechanical vibration signal further comprises the following steps: and determining the phase size of the main frequency component of the vibration signal of the rotating machine according to the preset frequency distinguishing interval between the main frequency components, the amplitude spectrum and the complex sequence.
In specific implementation, the method for determining the phase magnitude may be: when the current frequency component is determined to be a frequency component that can be counted according to the frequency differentiation interval (i.e. the current frequency component does not belong to the repeatedly counted frequencies f2 and f3 as shown in fig. 2), the position number corresponding to the current frequency component is found out from the complex sequence, and the phase magnitude of the current frequency component can be determined according to the real part and the imaginary part corresponding to the position number in the complex sequence (see the description of the following embodiments for details). The technical scheme provided by the embodiment of the invention determines the phase size of the main frequency component by distinguishing the interval of the frequency and the amplitude of the main frequency component, can more accurately perform fault reason and fault location analysis by combining the phase size of the main frequency component on the basis of the frequency and the amplitude of the main frequency component, provides technical guarantee for timely processing the fault and recovering the operation as early as possible, and ensures the high-efficiency operation of the equipment.
As the inventor finds that: in the prior art, the following data processing and report compiling need to manually input the main frequency component information, so that the inventor proposes the following implementation scheme in consideration of the technical problem.
In one embodiment, the method for determining the information of the main frequency components of the vibration signal of the rotating machine further comprises: and processing the main frequency component information and outputting the main frequency component information to a document. Wherein the document may be: excel documents or txt files.
In specific implementation, after the main frequency component information of the vibration signal of the rotary machine is determined according to the frequency distinguishing interval, the main frequency component information can be visually displayed on a frequency domain amplitude spectrogram, and the determination result of the main frequency component information is output to an Excel document or a txt document, so that the step of manually inputting the main frequency component information in subsequent data processing and report compiling is omitted, the labor cost is saved, and the working efficiency is improved.
In addition, the inventors have also found that: the specific implementation process also relates to the condition of main frequency component information statistics of multi-working-condition and multi-channel data. In view of this technical problem, the inventors propose the following embodiments.
In one embodiment, the method for determining the information of the main frequency components of the vibration signal of the rotating machine further comprises:
executing a batch processing program to determine main frequency component information corresponding to each working condition and channel for the main frequency component information of the multi-working condition and/or multi-channel rotating machinery vibration signal data;
and processing the main frequency component information corresponding to each working condition and channel, and outputting the main frequency component information to a document.
During specific implementation, when the main frequency component information of a plurality of working conditions and a plurality of channels is counted, a batch processing program is executed, the main frequency components of different working conditions and channels are determined respectively, then the statistical result is output to an Excel document or a txt file, a large amount of manual statistical input work is saved, wrong statistical results caused by human errors are avoided, the accuracy of determining the main frequency component information of the vibration signal of the rotary machine is improved, and meanwhile, the efficiency of determining the main frequency component information of the vibration signal of the rotary machine is improved.
The following is a complete flow chart illustrating how the present invention may be implemented.
(1) Obtaining a sampling signal x (t) (generally an integral multiple of a rotation period) of a time domain waveform of vibration of the rotating machinery for a certain time length, and assuming that the sampling frequency of the signal is fs。
(2) The separation interval JG between the main frequency components is set (can be preset according to actual needs, and has high flexibility).
(3) And setting the main frequency component order JS to be automatically counted, (which can be preset according to actual needs and has high flexibility).
(4) And subtracting the self average value from the time domain waveform sampling signal x (t), then performing Fast Fourier Transform (FFT) to obtain a complex sequence FFTX, and further calculating to obtain an amplitude spectrum SAmp.
(5) In specific implementation, the following procedure is executed to obtain the position number in the FFTX array (sequence) corresponding to the JS-order primary frequency component.
L, calculating the array length corresponding to the time domain waveform sampling signal x (t), wherein L represents the array length corresponding to the time domain waveform sampling signal x (t), and length (x (t)) represents the function of calculating the array length corresponding to x (t);
Deltf＝fsl', amplitude spectrum frequency line resolution, wherein Deltf represents frequency line resolution;
calculating to obtain corresponding sequence number position intervals (reference position intervals) in the amplitude spectrum array according to the frequency distinguishing intervals; wherein G represents the corresponding sequence number position interval in the amplitude spectrum array; round (JG/Deltf) represents a function for solving the position interval of the corresponding sequence number in the amplitude spectrum array;
for i＝1:JS
[ dx (i), xh (i) ] ═ max (samp); the max function returns the maximum amplitude of the array and the corresponding position serial number; wherein i represents the serial number of the main frequency component, dx (i) represents the amplitude of the ith main frequency component, and xh (i) represents the position serial number corresponding to the ith main frequency component; max (samp) represents a function for finding the maximum amplitude value and the corresponding position number;
if Xh(i)-G<1；
SAmp(1:Xh(i)+G)＝0；
elseif Xh(i)+G>length(SAmp)；
SAmp(Xh(i)-G:length(SAmp))＝0；
else；
amp (xh (i) -G: xh (i) + G) ═ 0; assigning a value of the maximum value to be zero at a distance of G positions from both sides; this step effectively avoids the repeated statistics of the same dominant frequency component as shown in fig. 2, thus reducing the number of useful statistics.
end; the calculation is circulated until the positions corresponding to the JS main frequency components are obtained;
end。
(6) the Dx array is the amplitude of JS main frequency components, the Xh array is the sequence number position corresponding to JS main frequency components, and the following procedures are executed to obtain the size and the phase of the main frequency components:
For i＝1:JS
mainf (i) ═ (xh (i) -1) × Deltf; ' Mainf is a JS array of principal frequency components;
phase (i) 180 × angle (FFTX (xh (i))/pi; phase is the corresponding Phase of JS main frequency components;
End。
(7) and displaying the statistical result of the main frequency components on the amplitude spectrogram, and outputting the statistical result to an Excel document or a txt file.
(8) And for the main frequency component information statistics of the multi-working-condition and multi-channel data, executing a batch processing program to automatically count the main frequency component information, and outputting the result to an Excel document or a txt file.
The following is a detailed example to illustrate how the invention may be practiced.
And selecting a vibration signal of a certain rotating machine for example demonstration.
(1) The rotating speed of the rotating mechanical unit is 200 rpm, the rotating period is 0.3 second, data of 6 seconds are selected, and time domain waveforms are shown in fig. 4. In fig. 4, the abscissa represents the sampling time and the ordinate represents the amplitude.
(2) Setting the interval of the main frequency component partition to be 1Hz and the statistical order of the main frequency component to be 5.
(3) The frequency domain amplitude spectrogram obtained by the method for automatically determining the main frequency component information of the vibration signal of the rotary machine is shown in fig. 5.
(4) In order to visually display more information of the main frequency components, the statistical order of the main frequency components may be adjusted to 10, and the obtained frequency domain amplitude spectrogram is shown in fig. 6. In fig. 5 and 6, the abscissa represents frequency and the ordinate represents amplitude.
Based on the same inventive concept, the embodiment of the present invention further provides an apparatus for determining the main frequency component information of the vibration signal of the rotating machine, as described in the following embodiments. Since the principle of the apparatus for determining the primary frequency component information of the vibration signal of the rotary machine for solving the problem is similar to the method for determining the primary frequency component information of the vibration signal of the rotary machine, the implementation of the apparatus for determining the primary frequency component information of the vibration signal of the rotary machine can be referred to the implementation of the method for determining the primary frequency component information of the vibration signal of the rotary machine, and repeated details are omitted. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware, or a combination of software and hardware is also possible and contemplated.
Fig. 7 is a schematic structural diagram of an apparatus for determining information of main frequency components of a vibration signal of a rotating machine according to an embodiment of the present invention, as shown in fig. 7, the apparatus includes:
the acquisition unit 02 is used for acquiring a sampling signal and a sampling frequency of a rotary mechanical vibration time domain waveform with preset duration;
an amplitude spectrum determining unit 04, configured to perform fast fourier transform on the time domain waveform sampling signal to obtain a complex sequence, and obtain an amplitude spectrum according to the complex sequence;
a main frequency component information determining unit 06 for determining main frequency component information of the vibration signal of the rotary machine based on a frequency discrimination interval between preset main frequency components and the amplitude spectrum; the main frequency component information includes: the frequency and amplitude of the dominant frequency component.
In one example, the primary frequency component information may further include: the phase magnitude of the primary frequency component;
the apparatus for determining information of the main frequency component of the vibration signal of the rotating machine may further include: and the phase size determining unit is used for determining the phase size of the main frequency component of the rotary mechanical vibration signal according to the preset frequency distinguishing interval between the main frequency components, the amplitude spectrum and the complex sequence.
In one example, the primary frequency component information determining unit may be specifically configured to:
determining the resolution of a frequency spectral line of the amplitude spectrum according to the array length and the sampling frequency corresponding to the time-domain waveform sampling signal;
determining a reference position interval in the amplitude spectrum array according to the resolution of the frequency spectral line of the amplitude spectrum and a preset frequency distinguishing interval between main frequency components;
and determining whether the frequency components to be confirmed are counted as the main frequency components according to the reference position interval in the amplitude spectrum array, the position serial number calculated in the amplitude spectrum array corresponding to the frequency components to be confirmed and the position serial number corresponding to the former confirmed main frequency component in the amplitude spectrum.
In an embodiment, determining whether the frequency component to be confirmed is counted as the main frequency component according to the reference position interval in the amplitude spectrum array, the position serial number calculated in the amplitude spectrum array corresponding to the frequency component to be confirmed, and the position serial number corresponding to the previously confirmed main frequency component in the amplitude spectrum may include:
a step of determining a first main frequency component; the step of determining a first dominant frequency component comprises: acquiring the maximum amplitude and the corresponding frequency component in the amplitude spectrum array, and calculating the corresponding position serial number of the frequency component in the amplitude spectrum array;
a step of determining other main frequency components; the step of determining other major frequency components comprises: assigning the numerical values of the determined main frequency components on the corresponding position serial numbers in the amplitude spectrum array and all the position serial numbers in the interval of the left and right adjacent reference positions to be zero to obtain an updated amplitude spectrum array; acquiring the maximum amplitude and the corresponding frequency component in the updated amplitude spectrum array, calculating the corresponding position serial number of the frequency component in the updated amplitude spectrum array, and counting the frequency component as a main frequency component;
and continuing to perform the step of determining other primary frequency components based on the updated array of magnitude spectra until all primary frequency component information is found.
In one example, the apparatus for determining information of the main frequency component of the vibration signal of the rotating machine may further include: and the output unit is used for processing the main frequency component information and outputting the main frequency component information to a document.
In one example, the apparatus for determining information of the main frequency component of the vibration signal of the rotating machine may further include: a batch processing unit to:
executing a batch processing program to determine main frequency component information corresponding to each working condition and channel for the main frequency component information of the multi-working condition and/or multi-channel rotating machinery vibration signal data;
and processing the main frequency component information corresponding to each working condition and channel, and outputting the main frequency component information to a document.
The embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor realizes the method for determining the main frequency component information of the rotating mechanical vibration signal when executing the computer program.
Embodiments of the present invention also provide a computer-readable storage medium storing a computer program for executing the method for determining information of a main frequency component of a vibration signal of a rotating machine.
The technical scheme provided by the implementation of the invention can achieve the following beneficial technical effects:
(1) the interval and the order can be distinguished according to the set frequency, the main frequency component information can be automatically and accurately counted, the frequency domain amplitude spectrogram is visually displayed, and the result is output to an Excel document or a txt file, so that the accuracy of determining the main frequency component information of the vibration signal of the rotary machine is improved, the fault reason and the fault location analysis can be more accurately performed, the technical guarantee is provided for timely processing the fault and early recovering the equipment operation, and the high-efficiency operation of the equipment is ensured.
(2) When the main frequency component information of a plurality of working conditions and a plurality of channels is counted, the batch processing program is executed, and the counting result is output to the Excel document or the txt file, so that a large amount of manual counting and inputting work is saved, wrong counting results caused by human errors are avoided, and the working quality and the working efficiency are improved.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above description is only a preferred embodiment of the present invention, and is not intended to limit the present invention, and various modifications and changes may be made to the embodiment of the present invention by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (8)
1. A method of determining information about the dominant frequency components of a rotating machine vibration signal, comprising:
acquiring a rotating mechanical vibration time domain waveform sampling signal with preset duration and sampling frequency;
carrying out fast Fourier transform on the time domain waveform sampling signal to obtain a complex sequence, and obtaining an amplitude spectrum according to the complex sequence;
determining the main frequency component information of the vibration signal of the rotating machinery according to the preset frequency distinguishing interval between the main frequency components and the amplitude spectrum; the main frequency component information includes: the frequency and amplitude of the primary frequency component;
determining the main frequency component information of the vibration signal of the rotating machinery according to the preset frequency distinguishing interval between the main frequency components and the amplitude spectrum, wherein the method comprises the following steps:
determining the resolution of a frequency spectral line of the amplitude spectrum according to the array length and the sampling frequency corresponding to the time-domain waveform sampling signal;
determining a reference position interval in the amplitude spectrum array according to the resolution of the frequency spectral line of the amplitude spectrum and a preset frequency distinguishing interval between main frequency components;
determining whether the frequency components to be confirmed are counted as the main frequency components according to the reference position interval in the amplitude spectrum array, the position serial number calculated in the amplitude spectrum array corresponding to the frequency components to be confirmed and the position serial number corresponding to the former confirmed main frequency component in the amplitude spectrum;
determining whether the frequency components to be confirmed are counted as the main frequency components according to the reference position interval in the amplitude spectrum array, the position serial number calculated in the amplitude spectrum array corresponding to the frequency components to be confirmed and the position serial number corresponding to the former confirmed main frequency component in the amplitude spectrum, wherein the step comprises the following steps:
a step of determining a first main frequency component; the step of determining a first dominant frequency component comprises: acquiring the maximum amplitude and the corresponding frequency component in the amplitude spectrum array, and calculating the corresponding position serial number of the frequency component in the amplitude spectrum array;
a step of determining other main frequency components; the step of determining other major frequency components comprises: assigning the numerical values of the determined main frequency components on the corresponding position serial numbers in the amplitude spectrum array and all the position serial numbers in the interval of the left and right adjacent reference positions to be zero to obtain an updated amplitude spectrum array; acquiring the maximum amplitude and the corresponding frequency component in the updated amplitude spectrum array, calculating the corresponding position serial number of the frequency component in the updated amplitude spectrum array, and counting the frequency component as a main frequency component;
and continuing to perform the step of determining other primary frequency components based on the updated array of magnitude spectra until all primary frequency component information is found.
2. The method of determining principal frequency component information of a rotating machine vibration signal as claimed in claim 1, wherein said principal frequency component information further comprises: the phase magnitude of the primary frequency component;
the method for determining the main frequency component information of the rotary mechanical vibration signal further comprises the following steps: and determining the phase size of the main frequency component of the vibration signal of the rotating machine according to the preset frequency distinguishing interval between the main frequency components, the amplitude spectrum and the complex sequence.
3. The method of determining principal frequency component information of a rotating machine vibration signal as set forth in claim 1, further including: and processing the main frequency component information and outputting the main frequency component information to a document.
4. The method of determining principal frequency component information of a rotating machine vibration signal as set forth in claim 1, further including:
executing a batch processing program to determine main frequency component information corresponding to each working condition and channel for the main frequency component information of the multi-working condition and/or multi-channel rotating machinery vibration signal data;
and processing the main frequency component information corresponding to each working condition and channel, and outputting the main frequency component information to a document.
5. An apparatus for determining information about the dominant frequency components of a rotating machine vibration signal, comprising:
the acquisition unit is used for acquiring a rotary mechanical vibration time domain waveform sampling signal with preset duration and sampling frequency;
the amplitude spectrum determining unit is used for carrying out fast Fourier transform on the time domain waveform sampling signal to obtain a complex sequence and obtaining an amplitude spectrum according to the complex sequence;
the main frequency component information determining unit is used for determining main frequency component information of the vibration signal of the rotating machinery according to preset frequency distinguishing intervals among main frequency components and the amplitude spectrum; the main frequency component information includes: the frequency and amplitude of the primary frequency component;
the primary frequency component information determining unit is specifically configured to:
determining the resolution of a frequency spectral line of the amplitude spectrum according to the array length and the sampling frequency corresponding to the time-domain waveform sampling signal;
determining a reference position interval in the amplitude spectrum array according to the resolution of the frequency spectral line of the amplitude spectrum and a preset frequency distinguishing interval between main frequency components;
determining whether the frequency components to be confirmed are counted as the main frequency components according to the reference position interval in the amplitude spectrum array, the position serial number calculated in the amplitude spectrum array corresponding to the frequency components to be confirmed and the position serial number corresponding to the former confirmed main frequency component in the amplitude spectrum;
determining whether the frequency components to be confirmed are counted as the main frequency components according to the reference position interval in the amplitude spectrum array, the position serial number calculated in the amplitude spectrum array corresponding to the frequency components to be confirmed and the position serial number corresponding to the former confirmed main frequency component in the amplitude spectrum, wherein the step comprises the following steps:
a step of determining a first main frequency component; the step of determining a first dominant frequency component comprises: acquiring the maximum amplitude and the corresponding frequency component in the amplitude spectrum array, and calculating the corresponding position serial number of the frequency component in the amplitude spectrum array;
a step of determining other main frequency components; the step of determining other major frequency components comprises: assigning the numerical values of the determined main frequency components on the corresponding position serial numbers in the amplitude spectrum array and all the position serial numbers in the interval of the left and right adjacent reference positions to be zero to obtain an updated amplitude spectrum array; acquiring the maximum amplitude and the corresponding frequency component in the updated amplitude spectrum array, calculating the corresponding position serial number of the frequency component in the updated amplitude spectrum array, and counting the frequency component as a main frequency component;
and continuing to perform the step of determining other primary frequency components based on the updated array of magnitude spectra until all primary frequency component information is found.
6. The apparatus for determining principal frequency component information of a rotating machine vibration signal according to claim 5, wherein said principal frequency component information further includes: the phase magnitude of the primary frequency component;
the device for determining the main frequency component information of the rotary mechanical vibration signal further comprises: and the phase size determining unit is used for determining the phase size of the main frequency component of the rotary mechanical vibration signal according to the preset frequency distinguishing interval between the main frequency components, the amplitude spectrum and the complex sequence.
7. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1 to 4 when executing the computer program.
8. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810938946.0A CN108709635B (en) | 2018-08-17 | 2018-08-17 | Method and device for determining main frequency component information of rotary mechanical vibration signal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810938946.0A CN108709635B (en) | 2018-08-17 | 2018-08-17 | Method and device for determining main frequency component information of rotary mechanical vibration signal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108709635A CN108709635A (en) | 2018-10-26 |
| CN108709635B true CN108709635B (en) | 2020-07-28 |
Family
ID=63873624
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810938946.0A Active CN108709635B (en) | 2018-08-17 | 2018-08-17 | Method and device for determining main frequency component information of rotary mechanical vibration signal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108709635B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109657971B (en) * | 2018-12-14 | 2021-05-28 | 浙江博远电子科技有限公司 | Data processing method and device based on cloud platform |
| CN111380607A (en) * | 2018-12-28 | 2020-07-07 | 财团法人工业技术研究院 | Motor vibration or noise frequency detection device and method thereof |
| EP3712577B1 (en) * | 2019-03-22 | 2023-07-26 | ABB Schweiz AG | Apparatus for equipment monitoring |
| CN110161363B (en) * | 2019-06-20 | 2021-04-06 | 广东石油化工学院 | Transformer running state vibration and sound detection method and system based on main frequency characteristic quantity |
| CN112525334B (en) * | 2020-11-18 | 2023-06-02 | 西安因联信息科技有限公司 | Dynamic equipment vibration multistable detection method |
| CN113188644A (en) * | 2021-03-31 | 2021-07-30 | 大唐东营发电有限公司 | Harmonic vibration amplitude calculation method for rotary mechanical vibration test system |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005214777A (en) * | 2004-01-29 | 2005-08-11 | Nsk Ltd | Frequency spectral analyzer |
| CN102095885A (en) * | 2010-10-28 | 2011-06-15 | 华南理工大学 | Vibration measuring method of engine speed by applying four-point energy centrobaric correction method |
| JP2011127527A (en) * | 2009-12-18 | 2011-06-30 | Tokyo Electric Power Co Inc:The | Vibration monitoring device for rotor |
| CN103234627A (en) * | 2013-04-17 | 2013-08-07 | 国家电网公司 | Complete alternation synchronous sampling and analyzing method for rotating machinery vibration signals |
| CN103323274A (en) * | 2013-05-24 | 2013-09-25 | 上海交通大学 | Rotating machinery condition monitoring and fault diagnosing system and method |
| CN105222950A (en) * | 2015-09-18 | 2016-01-06 | 中国水利水电科学研究院 | A kind of differentiating method of true and false frequency multiplication pressure fluctuation signal |
| CN104034412B (en) * | 2014-06-24 | 2016-06-08 | 西安交通大学 | A kind of rotary machine fault characteristic extraction method based on fractional order principle of holography |
| CN106092524A (en) * | 2016-05-13 | 2016-11-09 | 长兴昇阳科技有限公司 | A kind of method using vibration signal accurately to extract tach signal |
| CN107677362A (en) * | 2017-09-06 | 2018-02-09 | 浙江大学 | A kind of improved rotating machinery modulating frequency extracting method based on real cepstrum analysis |
| CN108151869A (en) * | 2017-11-27 | 2018-06-12 | 广州航新航空科技股份有限公司 | A kind of mechanical oscillation characteristic index extracting method, system and device |
| JP2018117431A (en) * | 2017-01-17 | 2018-07-26 | 株式会社トーエネック | Short-circuit diagnosis device for rotating machine and short-circuit diagnosis method for rotating machine |
-
2018
- 2018-08-17 CN CN201810938946.0A patent/CN108709635B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005214777A (en) * | 2004-01-29 | 2005-08-11 | Nsk Ltd | Frequency spectral analyzer |
| JP2011127527A (en) * | 2009-12-18 | 2011-06-30 | Tokyo Electric Power Co Inc:The | Vibration monitoring device for rotor |
| CN102095885A (en) * | 2010-10-28 | 2011-06-15 | 华南理工大学 | Vibration measuring method of engine speed by applying four-point energy centrobaric correction method |
| CN103234627A (en) * | 2013-04-17 | 2013-08-07 | 国家电网公司 | Complete alternation synchronous sampling and analyzing method for rotating machinery vibration signals |
| CN103323274A (en) * | 2013-05-24 | 2013-09-25 | 上海交通大学 | Rotating machinery condition monitoring and fault diagnosing system and method |
| CN104034412B (en) * | 2014-06-24 | 2016-06-08 | 西安交通大学 | A kind of rotary machine fault characteristic extraction method based on fractional order principle of holography |
| CN105222950A (en) * | 2015-09-18 | 2016-01-06 | 中国水利水电科学研究院 | A kind of differentiating method of true and false frequency multiplication pressure fluctuation signal |
| CN106092524A (en) * | 2016-05-13 | 2016-11-09 | 长兴昇阳科技有限公司 | A kind of method using vibration signal accurately to extract tach signal |
| JP2018117431A (en) * | 2017-01-17 | 2018-07-26 | 株式会社トーエネック | Short-circuit diagnosis device for rotating machine and short-circuit diagnosis method for rotating machine |
| CN107677362A (en) * | 2017-09-06 | 2018-02-09 | 浙江大学 | A kind of improved rotating machinery modulating frequency extracting method based on real cepstrum analysis |
| CN108151869A (en) * | 2017-11-27 | 2018-06-12 | 广州航新航空科技股份有限公司 | A kind of mechanical oscillation characteristic index extracting method, system and device |
Non-Patent Citations (2)
| Title |
|---|
| 《Vibration signal analysis using parameterized time–frequency method for features extraction of varying-speed rotary machinery》;Y.Yang;《Journal of Sound and Vibration》;20150120;全文 * |
| 《旋转机械变速工况下阶次全息方法研究》;汪华平;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20140315;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108709635A (en) | 2018-10-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108709635B (en) | Method and device for determining main frequency component information of rotary mechanical vibration signal | |
| US10393621B2 (en) | Method for assessing the condition of rotating machinery connected to an electric motor | |
| US5744723A (en) | Method for determining rotational speed from machine vibration data | |
| US8370109B2 (en) | Machine vibration baseline synthesizer | |
| JP6945371B2 (en) | Diagnostic equipment for rotating machine systems, power converters, rotating machine systems, and diagnostic methods for rotating machine systems. | |
| RU2300116C2 (en) | Mode of diagnostics od electrical engines of alternating current and of mechanical arrangements involved with them | |
| EP2574947A1 (en) | A method of determining stationary signals for the diagnostics of an electromechanical system | |
| KR101498527B1 (en) | Diagnosis system using vibration frequency analysis program for rotation equipment of power plant | |
| CN111401584A (en) | Automatic diagnosis method and system for equipment fault | |
| CN103941722B (en) | By component feature frequency multiplication amplitude Data Trend Monitor and the method for diagnostic device fault | |
| CN109765055A (en) | Rolling bearing fault testing method and system based on EWT, spectrum virtual value and KNN | |
| JP2022144024A (en) | Apparatus and program for detecting/predicting failure | |
| WO2014123443A1 (en) | Method and device for vibration diagnosis and forecasting sudden engine failure | |
| Rodriguez-Donate et al. | Wavelet-based general methodology for multiple fault detection on induction motors at the startup vibration transient | |
| MX2011006462A (en) | Rotating machine shaft signal monitoring method and system. | |
| CN117093945B (en) | Online monitoring and fault analysis system for rotary machine state | |
| CN112345827A (en) | Graphical differentiation of spectral frequency families | |
| CN109357751B (en) | Power transformer winding looseness defect detection system | |
| JP2018128284A (en) | State analysis device, display method, and program | |
| Cabal-Yepez et al. | FPGA-based online induction motor multiple-fault detection with fused FFT and wavelet analysis | |
| CN108520093B (en) | Mechanical equipment fault diagnosis method and device based on knowledge base | |
| CN102879084B (en) | On-line low-frequency vibration unstable state early warning method for steam turbine generator unit | |
| EP0105600A2 (en) | Apparatus for selecting a signal representing vibrations | |
| CN107436244B (en) | Equipment fault alarm method based on frequency segmentation vibration data acquisition | |
| JP2004020484A (en) | Abnormality monitoring device and program for monitoring abnormality |
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 |