CN108956117B - The minimizing technology of electric and magnetic oscillation component, Diagnosis of Rotating Machinery method and device - Google Patents
The minimizing technology of electric and magnetic oscillation component, Diagnosis of Rotating Machinery method and device Download PDFInfo
- Publication number
- CN108956117B CN108956117B CN201810787683.8A CN201810787683A CN108956117B CN 108956117 B CN108956117 B CN 108956117B CN 201810787683 A CN201810787683 A CN 201810787683A CN 108956117 B CN108956117 B CN 108956117B
- Authority
- CN
- China
- Prior art keywords
- frequency
- wave crest
- periphery
- rotating machinery
- frequency spectrum
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
-
- 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
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/34—Testing dynamo-electric machines
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
Problem of the present invention is that accurately removing the frequency component of electric and magnetic oscillation caused by inverter power supply from the signal of vibrating sensor for being installed on the rotating machinery driven as inverter power supply.For this purpose, carrying out Fourier transformation to the time of vibration waveform got by vibrating sensor, frequency spectrum is calculated.The frequency of the maximum wave crest on the periphery of the integral multiple of the carrier frequency of inverter power supply is set as reference frequency.The auto-correlation function of the frequency spectrum on calculating benchmark frequency periphery.Seek the interval of the wave crest of the auto-correlation function on reference frequency periphery.It will be as the wave crest existing for the peak separation of wave crest is extracted as object wave crest at equal intervals existing for the reference frequency periphery.Reduce the level of the frequency component of object wave crest in frequency spectrum.
Description
Technical field
The present invention relates to the minimizing technology of electric and magnetic oscillation component, Diagnosis of Rotating Machinery method and Diagnosis of Rotating Machinery dresses
It sets.
Background technique
The equipment in factory etc. stops in order to prevent, is determined the vibration of rotating machinery in the past to monitor abnormal set
Standby diagnosis (for example, referring to non-patent literature 1).
In recent years, it is used more and more by the motor of Driven by inverter.The electronic equipment that inverter mode drives
There is the advantage that following, that is, only by change setting (frequency modulating signal of inverter power supply), just can simply change
Variable speed is operated.But it is generated caused by carrier frequency when running well as the motor of Driven by inverter
Electric and magnetic oscillation, this in vibration diagnosis become noise and interfere vibration anomaly monitoring.Therefore, it is supervised in the vibration of rotating machinery
In view apparatus, following method is devised, that is, from the signal of the vibrating sensor measured, remove electricity caused by inverter
The frequency component of magnetic vibration, Lai Jinhang vibration diagnosis.
In patent document 1, a kind of minimizing technology of electric and magnetic oscillation component including following processing is disclosed.
Detect the given frequency on the basis of the frequency for the integral multiple for becoming the carrier frequency determined by inverter power supply
Multiple peak values in rate range.
, as benchmark peak value, to seek the reference peak and each peak value with the immediate peak value of the integral multiple of carrier frequency
Whole frequency intervals.
Reference frequency interval is determined according to frequency interval, extracts the peak for the integral multiple that frequency interval is benchmark frequency interval
Value is as removal object peak value.
The method of patent document 1 is that upward convex point is defined as wave crest and finds wave crest first, and from the crest location
Method as peak separation is sought, although ideal spectrum waveform is effectively, since actual data are with certain
Obtained from the sample rate of time carries out A/D transformation to the signal of vibrating sensor, and frequency spectrum is also discrete data, therefore by
There is error in each crest location and added up in FFT leaks (leakage) the problems such as, error is possible to become larger.In addition,
The method of patent document 1 it is also possible to by simple noise, be not that wave crest caused by inverter electric and magnetic oscillation is also detected as
Except object wave crest.
Patent document 2 discloses a kind of method of the frequency component of electric and magnetic oscillation caused by determining inverter.But it is right
For this method, due to the carrier frequency and modulating frequency according to inverter power supply, electromagnetic vibration caused by inverter is calculated
Dynamic frequency, and removal can not then be determined so if modulating frequency, that is, motor revolving speed is unknown by being set as removal object wave crest
Object wave crest.
Citation
Patent document
Patent document 1:JP speciallys permit No. 5565120
Patent document 2:JP special open 2016-116251
Non-patent literature
Non-patent literature 1: Ji Mingzhu on well, " live query To answer え Ru real trample vibratory drilling method To I Ru device diagnostic (return
Answer the device diagnostic based on actual vibration method of scene enquirement) ", Japanese plant protects association, in September, 1998
Summary of the invention
(subject to be solved by the invention)
Problem of the present invention is that from the signal for the vibrating sensor for being installed on the rotating machinery driven by inverter power supply
In, accurately remove the frequency component of electric and magnetic oscillation caused by inverter power supply.In addition, problem of the present invention is that, base
Diagnosis of Rotating Machinery is accurately proceed in the removal of such electric and magnetic oscillation component.
(means for solving the problems)
The 1st aspect of the present invention provides a kind of minimizing technology of electric and magnetic oscillation component, in the removal of the electric and magnetic oscillation component
In method, the vibration of the rotating machinery is obtained by being installed on by the vibrating sensor of the rotating machinery of inverter power supply driving
Dynamic time waveform carries out Fourier transformation to time of vibration waveform, frequency spectrum is calculated, by the inverter power supply in the frequency spectrum
The frequency of maximum wave crest on periphery of integral multiple of carrier frequency be set as reference frequency, calculate the institute on the reference frequency periphery
The auto-correlation function for stating frequency spectrum, the interval of the wave crest of the auto-correlation function by seeking the reference frequency periphery, thus
The peak separation of wave crest at equal intervals is sought existing for the reference frequency periphery, is extracted before and after the reference frequency every institute
Wave crest existing for peak separation is stated as object wave crest, and reduces the electricity of the frequency component of object wave crest described in the frequency spectrum
It is flat.The level of the object wave crest is e.g. reduced to wave crest by the reduction of the level of the frequency component of the object wave crest
The level of peak foot.
The 2nd aspect of the present invention provides a kind of Diagnosis of Rotating Machinery method, in the Diagnosis of Rotating Machinery method, to logical
The frequency spectrum after the removal of electric and magnetic oscillation component is carried out inverse Fourier transform by the minimizing technology for crossing the electric and magnetic oscillation component to be come
Time of vibration waveform is calculated, the time of vibration waveform as obtained from the inverse Fourier transform is based on, determines rotating machinery
State.
The 3rd aspect of the present invention provides a kind of Diagnosis of Rotating Machinery device, has: vibrating sensor, be installed on by
The rotating machinery of inverter power supply driving；Fourier transformation portion, the rotation to being got by the vibrating sensor
The time of vibration waveform of favourable turn tool carries out Fourier transformation to calculate frequency spectrum；Auto-correlation function calculation part, to reference frequency week
The auto-correlation function of the frequency spectrum on side is calculated, wherein the reference frequency is the inverter power supply in the frequency spectrum
Carrier frequency integral multiple periphery maximum wave crest frequency；Peak separation test section, by seeking the benchmark frequency
The interval of the wave crest of the auto-correlation function on rate periphery, to seek existing for the reference frequency periphery wave crest at equal intervals
Peak separation；Object wave crest test section extracts the wave crest existing for the peak separation before and after the reference frequency
As object wave crest；Level reduction portion reduces the level of the frequency component of object wave crest described in the frequency spectrum；Fourier
Inverse transformation portion, the frequency spectrum after being reduced to the level of the frequency component for making the object wave crest by level reduction portion into
Row inverse Fourier transform；And determination unit, it is based on the time of vibration waveform as obtained from the inverse Fourier transform, is determined
The state of rotating machinery.The reduction of the level of the frequency component of the object wave crest in level reduction portion is, for example, by institute
State object wave crest level be reduced to wave crest peak foot level.
(invention effect)
The minimizing technology of related electric and magnetic oscillation component according to the present invention can be driven from being installed on by inverter power supply
Rotating machinery vibrating sensor signal in accurately remove the frequency point of electric and magnetic oscillation caused by inverter power supply
Amount.In addition, Diagnosis of Rotating Machinery method involved according to the present invention and Diagnosis of Rotating Machinery device, by accurately going
Except electric and magnetic oscillation component caused by inverter power supply, so as to realize high-precision Diagnosis of Rotating Machinery.
Detailed description of the invention
Fig. 1 is the structure chart of Diagnosis of Rotating Machinery device involved in embodiments of the present invention.
Fig. 2 is the flow chart of the processing for illustrating to be executed by Diagnosis of Rotating Machinery device.
Fig. 3 A is the curve graph for showing acceleration time waveform.
Fig. 3 B is the curve graph for showing the frequency spectrum as obtained from the FFT of acceleration time waveform.
Fig. 3 C is the curve graph for showing the frequency spectrum near reference frequency.
Fig. 4 is the curve graph of the calculated auto-correlation function of frequency spectrum shown near reference frequency.
Fig. 5 is the curve graph of frequency spectrum near reference frequency that is showing extraction for illustrating object wave crest.
Fig. 6 is the curve graph for showing the frequency spectrum near the reference frequency after wave crest component is eliminated.
Acceleration obtained from Fig. 7 is the inverse FFT of the frequency spectrum near the reference frequency shown after being eliminated as wave crest component
The curve graph of time waveform.
Fig. 8 is the schematic diagram eliminated for illustrating wave crest component.
(symbol description)
1 Diagnosis of Rotating Machinery device
2 inverter power supplies
3 motor
4 piezoelectric acceleration sensors
5 processing units
6 pretreatment portions
7 storage units
8 operational parts
9 input units
10 output sections
11 amplifiers
12 bandpass filters
13 A/D converters
21 Fast Fourier Transform (FFT) portions (Fourier transformation portion)
22 auto-correlation function calculation parts
23 peak separation test sections
24 object wave crest test sections
25 wave crest component elimination portions (level reduction portion)
26 inverse fast Fourier transform portions (inverse Fourier transform portion)
27 determination units.
Specific embodiment
The embodiment of invention described below includes following method, that is, removal is installed on to be driven by inverter power supply
Rotating machinery vibrating sensor signal in include inverter power supply caused by electric and magnetic oscillation component.In this method
In, in order to determine the frequency component of electric and magnetic oscillation caused by inverter power supply, even if not extracting whole waves of certain frequency range
Peak separation is simultaneously sought and not with modulating frequency from its positional relationship for input in peak, also determines the object reduced as level
Wave crest interval.That is, the peak search caused by inverter power supply, this method has feature below.
Seek equally spaced multiple waves on the carrier frequency periphery as the feature of frequency component caused by inverter
The interval at peak (line frequency spectrum).
Modulating frequency is not used in the interval calculation of the wave crest.
Fig. 1 shows Diagnosis of Rotating Machinery device 1 involved in embodiments of the present invention.Become in the present embodiment and examines
The rotating machinery of disconnected object is the motor 3 driven by inverter power supply 2.
Diagnosis of Rotating Machinery device 1 in present embodiment has: piezoelectric acceleration sensor (vibrating sensor) 4,
It is installed on the bearing portion of motor 3；And processing unit 5, the signal from piezoelectric acceleration sensor 4 is handled.
Processing unit 5 carries out necessary pretreated pretreatment to the output from piezoelectric acceleration sensor 4 in addition to having
Except portion 6, it is also equipped with storage unit 7, operational part 8, input unit 9 and output section 10.Processing unit 5 can be by addition to cpu
Hardware and software mounted therein also comprising storage device as RAM, ROM construct.
It is additional for making an uproar after the output from piezoelectric acceleration sensor 4 is amplified by amplifier 11 in pretreatment portion 6
The bandpass filter 12 of sound removal, and then A/D (analog/digital) transformation is carried out by A/D converter 13.By these, treated
The acceleration time waveform of the motor 3 obtained by piezoelectric acceleration sensor 4 is stored in storage unit 7.
Operational part 8 removes the electric and magnetic oscillation from inverter power supply 2 from the acceleration time waveform for be stored in storage unit 7
Component, and the judgement of the state according to the acceleration time waveform operating motor 3 for having been removed electric and magnetic oscillation component.Determine
As a result for example it is output to the output section 10 as display.
Operational part 8 in present embodiment has: Fast Fourier Transform (FFT) portion (Fourier transformation portion) 21, auto-correlation function
Calculation part 22, peak separation test section 23, object wave crest extraction unit 24, wave crest component elimination portion (level reduction portion) 25, quickly
Inverse Fourier transform portion (inverse Fourier transform portion) 26 and determination unit 27.
The summary of the processing executed by operational part 8 is shown in the flow chart (step S1~S13) of Fig. 2.Fast Fourier
Transformation component 21 executes step S1.Auto-correlation function calculation part 22 executes step S3.Peak separation test section 23 executes step S4.It is right
As wave crest extraction unit 24 executes step S5, S8.Wave crest component elimination portion 25 executes step S6, S9.Inverse fast Fourier transform portion
26 execute step S12.Determination unit 27 executes step S13.
Hereinafter, illustrating the processing executed by operational part 8 referring to Fig. 2.In the following description, as needed, join together
According to Fig. 3 A to Fig. 7.Fig. 3 A is an example of the acceleration time waveform got by piezoelectric acceleration sensor 4, inverter power supply 2
Carrier frequency be 12kHz, the case where modulating frequency (output frequency) of inverter power supply 2 is 20Hz.Fig. 3 B is to pass through to Fig. 7
To various waveforms obtained from the processing of the acceleration time waveform of Fig. 3 A in operational part 8.In the following description, exist by
Fig. 3 A to Fig. 7 is collectively referred to as the case where " reference example ".
Firstly, in step sl, carrying out Fast Fourier Transform (FFT) to acceleration time waveform, calculating frequency spectrum.Fig. 3 B is logical
Cross frequency spectrum obtained from the Fast Fourier Transform (FFT) of the acceleration time waveform of Fig. 3 A, as carrier frequency 12kHz and its 2
There is line frequency spectrum in the frequency band periphery of 24kHz again.Fig. 3 B is the figure for being nearby exaggerated the 12kHz of the frequency spectrum of Fig. 3 A,
There is the wave crest at 2 times of the interval 40Hz of multiple modulating frequency 20Hz as inverter power supply 2.In addition, with inverter electricity
There is the case where wave crest in the interval of the integral multiple of the modulating frequency in source, such as records in patent document 2.
Then, in step s 2, reference frequency fc is determined.Here, so-called reference frequency fc, is the inverter electricity in frequency spectrum
The frequency of the maximum wave crest on the periphery (for example, range of 12kHz ± 0.2kHz) of the integral multiple of the carrier frequency in source 2.Due to carrying
Wave frequency rate is known to the specification from inverter power supply 2, therefore rough frequency can be specified by user, can use input
Portion 9 is inputted.It can also be shown using frequency spectrum as image in output section 10, by user based on the display come designated carrier frequency
The frequency of the maximum wave crest on the periphery of the integral multiple of rate.In reference example, the maximum wave crest on specified carrier frequency periphery
Frequency be accurately 11984.25Hz, be set to reference frequency fc.
In step s3, auto-correlation function is calculated for frequency spectrum.Auto-correlation function is generally mostly used for discovery time waveform
Periodicity, but herein for frequency spectrum calculate auto-correlation.If calculating the auto-correlation function of frequency spectrum, the value of auto-correlation function is pressed
According to the equally spaced wave crest in frequency spectrum each peak separation and become sharp wave crest, sought by program processing between wave crest
Every becoming easy.
The frequency spectrum of the computing object of auto-correlation function need not be set as entire scope, as long as by the significant reference frequency of its wave crest
The peripheral extent of fc is set as object.In the case where reference example, due to the modulating frequency (output frequency) of inverter power supply 2
Set maximum value as 60Hz, therefore the interval of desirable wave crest is up to 120Hz.It also takes out at this moment and calculates benchmark frequency
The range of the front and back 600Hz of rate fc.In addition, front and back of the range of the lag of auto-correlation function similarly with reference frequency fc
600Hz is suitable, and is set as 200 point part (frequency resolutions of the Fast Fourier Transform (FFT) in reference example as index (index)
For 3.05Hz).
Then, in step s 4, according to the auto-correlation function obtained by step S3, seek depositing on the periphery of reference frequency fc
Wave crest at equal intervals peak separation P.
Fig. 4 be to the auto-correlation function on the reference frequency fc frequency spectrum periphery of the frequency spectrum (Fig. 3 B, Fig. 3 C) in reference example just
The figure drawn of lag side.If successively seeking the index (lag) of the position of the wave crest (upward convex) of the auto-correlation function
Difference, then become 13,13,13,13,14,13 ....In reference example, mode 13 peak separation P has been set as.Wave crest
Interval P is also possible to the average value of the index difference of the position of the wave crest of auto-correlation function.In reference example, due in quick Fu
The frequency resolution of leaf transformation is 3.05Hz, therefore peak separation frequency becomes the interval 3.05 × 13=39.65Hz, is able to confirm that
40Hz with 2 times of the modulating frequency 20Hz of inverter is substantially uniform.
Then, in step s 5, for frequency spectrum, each peak separation P is extracted in the frequency side bigger than reference frequency fc
The wave crest at place is as object wave crest.Since peak separation P is discrete value (index unit) and with ± 1 error, only
Will make index advanced peak separation P (in reference example for 13) ± 1 in the range of there are wave crests, then become object wave crest.
Then, execute in frequency spectrum the level of the frequency component of object wave crest reduction, i.e., execution object wave crest point and its
The elimination (elimination of wave crest component) of the frequency component of the point of front and back.
Referring to Fig. 8, illustrate an example that wave crest component is eliminated.Fig. 8 schematically shows a part of frequency spectrum.In addition, in Fig. 8
In, symbol XkIndicate object wave crest.Firstly, seeking from object wave crest XkConsecutive points Xk+1To Xk+nN (be in this embodiment 5
It is a) average value A and standard deviation.If Xk+1It is greater than standard deviation with the absolute value of the difference of average value A, then next to Xk+2
To Xk+n+1It similarly calculated, determined, repeatedly the processing, in the time with the difference of average value lower than the point appearance of standard deviation
The point is set as 1 point of the peak foot of wave crest by point.For the consecutive points X of object wave crestk-1To Xk-nAlso it is similarly handled, is determined
Another extracts the point of the peak foot of wave crest.In fig. 8, wave crest Xk+2、Xk-2For peak pin point.Then, to the point (wave crest of two peak feet
Xk+2、Xk-2) linear interpolation is carried out, eliminate the frequency component between two peak feet.By (taking absolute value for original Fourier transformation result
Plural number before) ingredient of real number and imaginary number makes its reduction multiplied by the ratio reduced by interpolation.Symbol L indicates to carry out
The straight line of interpolation.
Step S5, S6 repeats specified number.That is, executing step S5, S6 for whole object wave crests.
Then, about frequency spectrum, the side smaller than reference frequency fc for frequency executes processing identical with step S5, S6,
Execute the extraction of object wave crest and the elimination (step S8~S10) of wave crest component.
If being set with other reference frequencies fc, the processing (step of step S2~S11 is executed to reference frequency fc
S11)。
In Fig. 5, circle mark is added to the object wave crest in frequency spectrum to show.
Fig. 6 is to reduce extracted wave crest, carries out linear interpolation to two peak feet, and the frequency component between two peak feet is disappeared
Frequency spectrum after removing.Original frequency spectrum is also shown with dotted line lightly.More than, the removal of electric and magnetic oscillation component is completed.
Then, in step s 12, the frequency spectrum (being Fig. 6 in reference example) after eliminating to wave crest component carries out in quick Fu
Leaf inverse transformation is back to acceleration time waveform.Fig. 7 shows the removal of the electric and magnetic oscillation component in reference example treated acceleration
Waveform, and light the acceleration time waveform shown before removal processing.If comparing the acceleration amplitude of removal before and after the processing
RMS value, then remove before processing as 2.48m/s2, and be 0.64m/s after removal processing2, being able to confirm that has reduced 74%.
Then, in step s 13, using electric and magnetic oscillation component removal treated acceleration time waveform (in reference example
In be Fig. 7), to execute judgement relevant to the vibrational state of motor 3.It is a variety of known to the example of such judgement, such as
It is judged as in the case where acceleration amplitude has been more than preset threshold value and abnormal vibrations has occurred.It can also will determine result
It is output to output section 10.Electric and magnetic oscillation component removal treated acceleration time waveform can also be by various filters
Reason, is used in simple diagnosis, accurate diagnosis (such diagnosis is recorded in non-patent literature 1 etc.).
It according to the present embodiment, can be with inverter in the vibration diagnosis of the rotating machinery driven by inverter power supply 2
The revolving speed of modulating frequency, that is, motor 3 of power supply 2 independently accurately remove piezoelectric acceleration sensor 4 signal in include
Inverter power supply 2 caused by electric and magnetic oscillation component.As a result, it is possible to which the benchmark of previous vibration diagnosis is directly applied to
In addition to the data of electric and magnetic oscillation component, it is able to carry out more accurate device diagnostic.
In embodiments, it is illustrated by taking piezoelectric acceleration sensor 4 as an example, but power type speed is passed
Sensor, non-contact displacement transducer etc. also can be using the present invention.In embodiments, with the output signal of vibrating sensor
Be be illustrated in case where vibration acceleration but it is also possible to be the speed integrated to vibration acceleration or
The displacement of further progress integral.Sensor-disposed portion position be set as the bearing portion of motor but it is also possible to be motor with
Outer rotating machinery, sensor-disposed portion position are also possible to the firm housing section of rotating machinery.
Claims (5)
1. a kind of minimizing technology of electric and magnetic oscillation component, wherein
The vibration of the rotating machinery is obtained by being installed on by the vibrating sensor of the rotating machinery of inverter power supply driving
Time waveform,
Fourier transformation is carried out to time of vibration waveform, calculates frequency spectrum,
The frequency of the maximum wave crest on the periphery of the integral multiple of the carrier frequency of the inverter power supply in the frequency spectrum is set as
Reference frequency,
The auto-correlation function of the frequency spectrum on the reference frequency periphery is calculated,
By seeking the interval of the wave crest of the auto-correlation function on the reference frequency periphery, to seek in the benchmark frequency
The peak separation of wave crest, the peak separation are the auto-correlation letters on the reference frequency periphery at equal intervals existing for rate periphery
The mode or average value of the difference of the position of several wave crests,
The wave crest existing for the peak separation is extracted before and after the reference frequency as object wave crest,
And reduce the level of the frequency component of object wave crest described in the frequency spectrum.
2. the minimizing technology of electric and magnetic oscillation component according to claim 1, wherein
The reduction of the level of the frequency component of the object wave crest is the peak foot that the level of the object wave crest is reduced to wave crest
Level.
3. a kind of Diagnosis of Rotating Machinery method, wherein
The frequency spectrum after removing electric and magnetic oscillation component to the minimizing technology of the electric and magnetic oscillation component by claims 1 or 2
It carries out inverse Fourier transform and calculates time of vibration waveform,
Based on the time of vibration waveform as obtained from the inverse Fourier transform, the state of rotating machinery is determined.
4. a kind of Diagnosis of Rotating Machinery device, has:
Vibrating sensor is installed on the rotating machinery driven by inverter power supply；
Fourier transformation portion, to the time of vibration waveform of the rotating machinery got by the vibrating sensor into
Row Fourier transformation calculates frequency spectrum；
Auto-correlation function calculation part calculates the auto-correlation function of the frequency spectrum on reference frequency periphery, wherein described
Reference frequency is the frequency of the maximum wave crest on the periphery of the integral multiple of the carrier frequency of the inverter power supply in the frequency spectrum；
Peak separation test section, the interval of the wave crest of the auto-correlation function by seeking the reference frequency periphery, from
And the peak separation of wave crest at equal intervals is sought existing for the reference frequency periphery, the peak separation is the reference frequency
The mode or average value of the difference of the position of the wave crest of the auto-correlation function on periphery；
Object wave crest test section extracts before and after the reference frequency wave crest existing for the peak separation as object
Wave crest；
Level reduction portion reduces the level of the frequency component of object wave crest described in the frequency spectrum；
Inverse Fourier transform portion, after being reduced to the level for the frequency component for making the object wave crest by level reduction portion
The frequency spectrum carries out inverse Fourier transform；With
Determination unit is based on the time of vibration waveform as obtained from the inverse Fourier transform, determines the state of rotating machinery.
5. Diagnosis of Rotating Machinery device according to claim 4, wherein
The reduction of the level of the frequency component of the object wave crest in level reduction portion is by the electricity of the object wave crest
Pancake as low as the peak foot of wave crest level.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2017229230A JP6420885B1 (en) | 2017-11-29 | 2017-11-29 | Method for removing electromagnetic vibration component, diagnostic method for rotating machine, and diagnostic device for rotating machine |
| JP2017-229230 | 2017-11-29 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108956117A CN108956117A (en) | 2018-12-07 |
| CN108956117B true CN108956117B (en) | 2019-11-08 |
Family
ID=64098847
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810787683.8A Active CN108956117B (en) | 2017-11-29 | 2018-07-18 | The minimizing technology of electric and magnetic oscillation component, Diagnosis of Rotating Machinery method and device |
Country Status (3)
| Country | Link |
|---|---|
| JP (1) | JP6420885B1 (en) |
| KR (1) | KR101966270B1 (en) |
| CN (1) | CN108956117B (en) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2020085603A (en) * | 2018-11-21 | 2020-06-04 | 東芝産業機器システム株式会社 | Method for measuring specific frequency |
| DE102019127211A1 (en) * | 2019-03-05 | 2020-09-10 | Computational Systems, Inc. | System for separating periodic amplitude peaks from non-periodic amplitude peaks in machine vibration data |
| US20210033450A1 (en) * | 2019-07-29 | 2021-02-04 | Computational Systems, Inc. | Vibration Waveform DC Disturbance Removal |
| KR20210017577A (en) | 2019-08-09 | 2021-02-17 | 주식회사 엘지화학 | Quantitatively diagnosis method for manufacturing facility |
| CN110530507B (en) * | 2019-08-29 | 2021-10-15 | 郑州大学 | Edge computing method, medium and system for rotating equipment monitoring |
| CN111964583B (en) * | 2020-07-08 | 2022-05-27 | 瑞声新能源发展(常州)有限公司科教城分公司 | Motor vibration displacement estimation method, device and medium |
| CN116490760A (en) * | 2020-11-16 | 2023-07-25 | 三菱电机株式会社 | Diagnostic device for motor |
| JP7611711B2 (en) | 2021-01-13 | 2025-01-10 | 株式会社Kelk | Analysis device, analysis system, and analysis method |
| JP7489927B2 (en) * | 2021-01-27 | 2024-05-24 | 公益財団法人鉄道総合技術研究所 | Diagnostic device and diagnostic method |
| CN114355187A (en) * | 2021-12-15 | 2022-04-15 | 武汉船用电力推进装置研究所（中国船舶重工集团公司第七一二研究所） | A method of extracting and analyzing motor vibration force wave test |
| CN116930853B (en) * | 2023-07-17 | 2024-12-24 | 海南电网有限责任公司 | A meter box aging detection method and device based on vibration analysis |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1457424A (en) * | 2001-02-28 | 2003-11-19 | 古河电气工业株式会社 | Method of measuring optical fiber drawing tensile force |
| JP5565120B2 (en) * | 2010-06-09 | 2014-08-06 | 富士電機株式会社 | High-frequency electromagnetic vibration component removal method and high-frequency electromagnetic vibration component removal device, rolling bearing diagnosis method and bearing diagnosis device for a rotating machine |
| CN104459676A (en) * | 2014-11-05 | 2015-03-25 | 上海大学 | System and method for simultaneously measuring lengths of two optical fibers |
| CN104599682A (en) * | 2015-01-13 | 2015-05-06 | 清华大学 | Method for extracting pitch period of telephone wire quality voice |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS587167B2 (en) | 1978-11-13 | 1983-02-08 | 沖電気工業株式会社 | Level measurement method |
| JPH03291539A (en) * | 1990-04-09 | 1991-12-20 | Toshiba Corp | Detecting method for abnormality of roller bearing of electric motor |
| US8544331B2 (en) * | 2008-12-04 | 2013-10-01 | University Of Ottawa | Parameter independent detection of rotating machinery faults |
| JP5738711B2 (en) * | 2011-07-29 | 2015-06-24 | 株式会社東芝 | Rotating machine state monitoring device, rotating machine state monitoring method, and rotating machine state monitoring program |
| JP2015126360A (en) * | 2013-12-26 | 2015-07-06 | 株式会社豊田中央研究所 | Signal processing apparatus and signal processing method |
| JP6160519B2 (en) * | 2014-03-07 | 2017-07-12 | 株式会社Ｊｖｃケンウッド | Noise reduction device |
| JP6450575B2 (en) * | 2014-12-10 | 2019-01-09 | 旭化成エンジニアリング株式会社 | Inverter noise elimination method and diagnostic method for equipment including inverter |
| JP6044647B2 (en) * | 2015-01-13 | 2016-12-14 | 株式会社明電舎 | Dynamometer control device and inertia moment estimation method using the same |
| KR101832244B1 (en) * | 2016-11-11 | 2018-02-27 | 한국과학기술연구원 | Bearing Test apparatus for testing behavior of the bearing |
-
2017
- 2017-11-29 JP JP2017229230A patent/JP6420885B1/en active Active
-
2018
- 2018-07-18 CN CN201810787683.8A patent/CN108956117B/en active Active
- 2018-09-19 KR KR1020180112090A patent/KR101966270B1/en active IP Right Grant
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1457424A (en) * | 2001-02-28 | 2003-11-19 | 古河电气工业株式会社 | Method of measuring optical fiber drawing tensile force |
| JP5565120B2 (en) * | 2010-06-09 | 2014-08-06 | 富士電機株式会社 | High-frequency electromagnetic vibration component removal method and high-frequency electromagnetic vibration component removal device, rolling bearing diagnosis method and bearing diagnosis device for a rotating machine |
| CN104459676A (en) * | 2014-11-05 | 2015-03-25 | 上海大学 | System and method for simultaneously measuring lengths of two optical fibers |
| CN104599682A (en) * | 2015-01-13 | 2015-05-06 | 清华大学 | Method for extracting pitch period of telephone wire quality voice |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2019100761A (en) | 2019-06-24 |
| JP6420885B1 (en) | 2018-11-07 |
| KR101966270B1 (en) | 2019-04-05 |
| CN108956117A (en) | 2018-12-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108956117B (en) | The minimizing technology of electric and magnetic oscillation component, Diagnosis of Rotating Machinery method and device | |
| EP2005125B1 (en) | A method and a system for monitoring the condition and operation of periodically moving objects | |
| CN102087138A (en) | System and method for vibration analysis and phase analysis of vibration waveforms using dynamic statistical averaging of tachometer data to accurately calculate rotational speed | |
| US10311703B1 (en) | Detection of spikes and faults in vibration trend data | |
| CN111256993A (en) | Method and system for diagnosing fault type of main bearing of wind turbine generator | |
| US10197594B2 (en) | Tachometer signal jitter reduction system and method | |
| JP6550476B2 (en) | Method for analyzing a signal and device for carrying out the method | |
| JP2021076533A (en) | Device and method for detecting rubbing of rotary machine | |
| JP6430234B2 (en) | Vibration analysis apparatus and program for rotating machine | |
| CN115496102A (en) | Wind turbine generator blade fault diagnosis method and device, equipment and storage medium | |
| JP7083293B2 (en) | Status monitoring method and status monitoring device | |
| WO2020202567A1 (en) | Vibro-acoustic analysis method and device and equipment-abnormality-location estimation method and device | |
| JP7394031B2 (en) | Abnormality detection device and abnormality detection method for rolling bearings | |
| CN105300688A (en) | RMS-based self-adaptive quick evaluating method for rotating speed of gearbox | |
| CN114061947A (en) | Fault diagnosis method and system for gearbox variable speed based on sparse time-frequency analysis | |
| JP2013160749A (en) | Facility diagnostic method and facility diagnostic device of rotary machine | |
| JP6305303B2 (en) | Vibration diagnosis apparatus, method and program for rotating equipment | |
| Henry et al. | Prism signal processing for machine condition monitoring I: Design and simulation | |
| CN110537082B (en) | Vibration detection device and abnormality determination system | |
| JP6036178B2 (en) | Respiratory sound analyzer, respiratory sound analysis method and respiratory sound analysis program | |
| JP6922065B1 (en) | Rotation speed detection system and rotation speed detection device for rotating equipment | |
| Klein et al. | Methods for diagnostics of bearings in non-stationary environments | |
| CN110866519B (en) | Rolling bearing fault diagnosis method based on Fourier decomposition and multiscale permutation entropy deviation value | |
| JP4209793B2 (en) | Abnormality diagnosis method based on acoustic signal and program used for executing the method | |
| Minemura et al. | Acoustic feature representation based on timbre for fault detection of rotary machines |
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 |