CN106646222B - A kind of three-phase cage asynchronous motor rotor broken splits the diagnostic method of conducting bar number - Google Patents
A kind of three-phase cage asynchronous motor rotor broken splits the diagnostic method of conducting bar number Download PDFInfo
- Publication number
- CN106646222B CN106646222B CN201610843046.9A CN201610843046A CN106646222B CN 106646222 B CN106646222 B CN 106646222B CN 201610843046 A CN201610843046 A CN 201610843046A CN 106646222 B CN106646222 B CN 106646222B
- Authority
- CN
- China
- Prior art keywords
- stator
- signal
- rotor
- phase
- conducting bar
- 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
- 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
- G01R31/346—Testing of armature or field windings
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02K—DYNAMO-ELECTRIC MACHINES
- H02K15/00—Methods or apparatus specially adapted for manufacturing, assembling, maintaining or repairing of dynamo-electric machines
- H02K15/02—Methods or apparatus specially adapted for manufacturing, assembling, maintaining or repairing of dynamo-electric machines of stator or rotor bodies
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Tests Of Circuit Breakers, Generators, And Electric Motors (AREA)
Abstract
The invention discloses the diagnostic methods that a kind of three-phase cage asynchronous motor rotor broken splits conducting bar number.According to the threephase stator voltage and current instantaneous signal of three-phase cage asynchronous motor, stator voltage virtual value, stator current virtual value, active power, apparent energy and power factor are obtained;Meanwhile by doing FFT spectrum analysis to electric current instantaneous signal, f in stator current is obtained1The virtual value I of frequency componentS;Meanwhile to electric current instantaneous signal i1Hilbert transform is carried out, and filters out the DC component of Hilbert modulus, and carry out FFT spectrum analysis, obtains the signal 2sf to be analyzed1The amplitude of frequency component, and the cage type asynchronous motor Hilbert modulus 2sf derived1Numerical value corresponding relationship between the amplitude and rotor rupture conducting bar number of frequency component diagnoses rotor rupture conducting bar number.The defect that MHMSA method is unable to quantitative Diagnosis rotor rupture conducting bar number is compensated for, for discovery in time and excludes rotor bar breaking fault, it is ensured that the safe and stable operation of motor creates advantage.
Description
Technical field
The present invention relates to technical field of motors, and in particular to a kind of three-phase cage asynchronous motor rotor broken splits conducting bar number
Diagnostic method.
Background technique
In three-phase cage asynchronous motor operational process, rotor bar is by radial electromagnetic force, electric rotating magnetic force, centrifugation
The effect of the alternate stresses such as power, thermal flexure amount of deflection power, rotor manufacturing defect, may cause broken bar fault, such failure hair in addition
Raw probability is about 10%.
Rotor broken bar is typical gradual failure, initial stage usual 1 bar failure, then with neighbouring its of fracture conducting bar
His conducting bar continues to be broken, and motor output drop is finally caused even to shut down.Therefore, it is necessary to implement rotor bar breaking fault detection and and
When to fracture conducting bar number make diagnosis.
After rotor bar breaking fault occurs for cage type asynchronous motor, will occur (1 ± 2s) f in its stator current1Frequency
Extra current component (s is revolutional slip, f1For frequency of supply), it can be as rotor bar breaking fault feature.And stator is electric
Stream signal is easy to acquire, therefore the stator current signal frequency spectrum analysis method based on Fast Fourier Transform (FFT) (FFT) is answered extensively
It is detected for rotor bar breaking fault.Such methods are commonly referred to as MCSA (Motor Current Signal Analysis) method.
Initially, MCSA directly carries out FFT spectrum analysis to stable state stator current signal, according in spectrogram whether there is (1
±2s)f1Frequency component judges rotor, and whether there is or not broken strips, and according to (1-2s) f in stator current1The virtual value I of frequency componentLWith f1
The virtual value I of frequency componentSDiagnosis is made to rotor rupture conducting bar number, as shown in formula (1).
Wherein, Nb, N be respectively rotor rupture conducting bar number and conducting bar sum.
But formula (1) is and the rotor bar breaking fault premised on ignoring rotor speed fluctuation (i.e. hypothesis invariablenes turning speed)
Necessarily rotor speed is caused to fluctuate, therefore formula (1) is that perfect condition does not tally with the actual situation, diagnoses rotor rupture accordingly and lead
Number is short of accuracy.
For this purpose, MCSA transfers to diagnose rotor rupture conducting bar number using formula (2).
Wherein, IRFor (1+2s) f in stator current1The virtual value of frequency component.Formula (2) is actually based on formula (1)
, but count and the fluctuation of speed, thus it is more in line with actual conditions.
Obviously, when MCSA diagnoses rotor rupture conducting bar number, no matter formula (1) or formula (1) are used, with the proviso that: it is necessary
Stator current (1-2s) f is clearlyed distinguish by FFT spectrum analysis1Frequency component, (1+2s) f1Frequency component and f1Frequency point
Amount.But when broken strip slight due to rotor, (1 ± 2s) f1The virtual value of component is relative to f1Component is very small, and asynchronous motor is transported
Revolutional slip s very little (in general, s < 0.05) when row, (1 ± 2s) f1With f1Numerical value is close, if directly doing FFT spectrum analysis,
Then (1 ± 2s) f1Component may be by f1The leakage of component is flooded.In this case, (1 ± 2s) f1The virtual value of frequency component
ILWith IRIt can not know, then formula (1) can not be applied with formula (2), and diagnosis rotor rupture conducting bar number is not known where to begin certainly.This
It is the shortcoming of MCSA.
To overcome stator current (1 ± 2s) f in MCSA method1Component is easy by f1The defect that the leakage of component is flooded, hair
Exhibition, which is formd, analyzes (Motor Hilbert Modulus Signal Analysis, MHMSA) based on Hilbert modulus signal
Cage type asynchronous motor rotor strip-broken failure detecting method, core is: in cage type asynchronous motor rotor normal condition
Under, Hilbert modulus is only a DC component;Once rotor bar breaking fault occurs, will occur in its Hilbert modulus
One additional 2sf1The wave component of frequency.Therefore, FFT spectrum point is done to the Hilbert modulus signal of cage type asynchronous motor
It analyses and investigates the 2sf1The detection of rotor bar breaking fault can be realized in the wave component of frequency.Great mass of data shows party's Faxian
Write be better than MCSA, can it is sensitive, reliably detect rotor bar breaking fault.
But for now, MHMSA method not yet solves the problems, such as the diagnosis of rotor rupture conducting bar number.Specifically,
Existing MHMSA method cannot quantitatively determine the number of rotor rupture conducting bar, and reason is: not yet discovery three-phase cage is asynchronous
Motor Hilbert modulus 2sf1Numerical value corresponding relationship between the amplitude and rotor rupture conducting bar number of frequency component.
Therefore, during application MHMSA method, how to diagnose three-phase cage asynchronous motor rotor broken and split conducting bar number
Mesh just becomes the problem faced at present in relation to technical staff.
Summary of the invention
In view of this, the present invention provides the diagnosis sides that a kind of three-phase cage asynchronous motor rotor broken splits conducting bar number
Method can find and exclude in time rotor bar breaking fault, it is ensured that the safe and stable operation of motor.
A kind of three-phase cage asynchronous motor rotor broken splits the diagnostic method of conducting bar number, specifically comprises the following steps:
Step 1: acquiring three-phase cage asynchronous motor threephase stator voltage transient signal and three by the sample frequency of setting
Phase stator current instantaneous signal;
Step 2: obtain m phase stator voltage virtual value, m indicate separate, m=1,2,3;
Step 3: obtaining the arithmetic mean of instantaneous value of threephase stator voltage effective value and as stator voltage virtual value;
Step 4: obtaining the virtual value of m phase stator current;
Step 5: obtaining the arithmetic mean of instantaneous value of threephase stator current effective value and as stator current virtual value;
Step 6: obtaining instantaneous active function according to threephase stator voltage transient signal and threephase stator electric current instantaneous signal
Rate;
Step 7: the method by calculating average value, obtains the DC component mean (p) of instantaneous active power p, i.e., it is active
Power;
Step 8: obtaining apparent energy according to stator voltage virtual value and stator current virtual value;Active power and view exist
Power ratio, as power factor
Step 9: doing FFT spectrum analysis to the electric current instantaneous signal of any one phase stator in threephase stator, stator is obtained
F in electric current1The virtual value I of frequency componentS;Wherein, f1For frequency of supply;
Step 10: carrying out Hilbert transform to electric current instantaneous signal;
Step 11: obtaining Hilbert modulus Mod;
Step 12: obtaining the DC component of Hilbert modulus Mod, i.e. average value A0;
Step 13: filtering out the DC component A in Hilbert modulus Mod0, obtain signal M to be analyzed;
Step 14: carrying out FFT spectrum analysis to signal M to be analyzed, its 2sf is obtained1The amplitude A of frequency component1;Wherein,
S is revolutional slip;
Step 15: according toIt is disconnected to obtain diagnosis rotor
Split conducting bar number Nb, wherein N is rotor bar sum.
Preferably, the threephase stator voltage transient signal and threephase stator electric current of the three-phase cage asynchronous motor are instantaneous
Signal is synchronous acquisition.
Preferably, the sample frequency of each signal is 1000Hz, when sampling a length of 10s, i.e. sampling number is 10000 points.
The utility model has the advantages that
The present invention is according to three-phase cage asynchronous motor Hilbert modulus 2sf1The amplitude of frequency component judges that rotor is disconnected
Number, compensates for the defect that traditional MHMSA method is unable to quantitative Diagnosis rotor rupture conducting bar number, for discovery in time and excludes
Rotor bar breaking fault, it is ensured that the safe and stable operation of motor creates advantage.
Detailed description of the invention
Fig. 1 is the electrical schematic diagram of present invention acquisition three-phase cage asynchronous motor stator voltage, current signal;
Fig. 2 is Y100L-2 type three-phase cage asynchronous motor (3kW, 380V) full load, in 1 broken strip of rotor
Hilbert modulus FFT spectrum;
Fig. 3 is Y100L-2 type three-phase cage asynchronous motor (3kW, 380V) full load, in 1 broken strip of rotor
Stator current FFT spectrum.
Specific embodiment
The present invention will now be described in detail with reference to the accompanying drawings and examples.
The present invention provides the diagnostic methods that a kind of three-phase cage asynchronous motor rotor broken splits conducting bar number, wherein each
Meaning of parameters are as follows:
S is revolutional slip;f1For frequency of supply;FFT is Fast Fourier Transform (FFT);The full name in English of MCSA is Motor
Current Signal Analysis;ILFor (1-2s) f in stator current1The virtual value of frequency component;ISFor in stator current
f1The virtual value of frequency component;NbFor rotor rupture conducting bar number;N is rotor bar sum;IRFor (1+2s) f in stator current1
The virtual value of frequency component;The full name in English of MHMSA is Motor Hilbert Modulus Signal Analysis;vmFor
M phase stator voltage instantaneous signal;imFor m phase stator current instantaneous signal;M be separate, m=1,2,3;VmFor m phase stator
Voltage effective value;vmIt (k) is k-th of sampled value of m phase stator voltage instantaneous signal;N is sampling number;VRMSFor stator electricity
It is pressed with valid value;ImFor m phase stator current virtual value;imIt (k) is k-th of sampled value of m phase stator current instantaneous signal;IRMS
For stator current virtual value;P is instantaneous active power;For power factor;For power-factor angle;For the 1st phase stator
Electric current instantaneous signal i1Hilbert transform;T is the time;τ is delay;Mod is Hilbert modulus;A0For Hilbert mould
Measure the average value of Mod, that is, DC component;Mean () represents operation of averaging;M is that Hilbert modulus Mod filters out direct current
Component A0Signal to be analyzed afterwards;A1Represent the 2sf of Hilbert modulus Mod1The amplitude of frequency component;iL' represent rotor broken bar
Stator current caused by failure (1-2s) f1Frequency component;IL' it is stator current (1-2s) f caused by rotor bar breaking fault1
The virtual value of frequency component;αLFor stator current caused by rotor bar breaking fault (1-2s) f1The initial phase angle of frequency component;E is
Natural constant;J is imaginary unit;Δ T is pulsating torque;P is number of pole-pairs;Ψ is the virtual value of fundamental flux;αΨFor fundamental wave magnetic
Logical initial phase angle;iL" for stator current (1-2s) f caused by torque pulsation1Frequency component;iRFor caused by torque pulsation
Stator current (1+2s) f1Frequency component;IL" for stator current (1-2s) f caused by torque pulsation1Frequency component it is effective
Value;IRFor stator current caused by torque pulsation (1+2s) f1The virtual value of frequency component;ZSFor the modulus value of stator impedance;αS
For stator impedance angle;J, rotary inertia;MsFor the 2sf of Hilbert modulus Mod1Frequency component;γ is phase angle;PT is voltage
Mutual inductor;CT is current transformer.
Step 1: acquiring three-phase cage asynchronous motor threephase stator voltage transient signal v by certain frequencymIt is fixed with three-phase
Electron current instantaneous signal im, wherein m indicate separate, m=1,2,3;
Step 2: obtaining m phase stator voltage v according to formula (3)mVirtual value Vm;
Wherein vm(k) k-th of sampled value of m phase stator voltage is represented, n is sampling number;
Step 3: obtaining the arithmetic mean of instantaneous value of threephase stator voltage effective value according to formula (4), having as stator voltage
Valid value VRMS;
Step 4: obtaining m phase stator current i according to formula (5)mVirtual value Im;
Wherein im(k) k-th of sampled value of m phase stator current is represented, n is sampling number;
Step 5: obtaining the arithmetic mean of instantaneous value of threephase stator current effective value according to formula (6), having as stator current
Valid value IRMS;
Step 6: according to threephase stator voltage transient signal vmWith threephase stator electric current instantaneous signal im, obtain instantaneous active
Power p;
Step 7: obtaining the DC component mean (p) of instantaneous active power p, i.e. active power using formula (8);
Step 8: obtaining power factor according to formula (9)
Wherein, 3VRMSIRMSFor apparent energy.
Step 9: to the electric current instantaneous signal i of any one phase stator in threephase stator1FFT spectrum analysis is done, is determined
F in electron current1The virtual value I of frequency componentS;Wherein, f1For frequency of supply;
Step 10: to electric current instantaneous signal i1Carry out Hilbert transform
Wherein, t is the time;τ is delay;
Step 11: obtaining Hilbert modulus Mod according to formula (11):
Step 12: obtaining DC component (i.e. average value) A of Hilbert modulus Mod0;
Step 13: filtering out the DC component A in Hilbert modulus Mod0, obtain signal M to be analyzed;
M=Mod-A0 (13)
Step 14: carrying out FFT spectrum analysis to signal M to be analyzed, its 2s f is obtained1The amplitude A of frequency component1;Its
In, s is revolutional slip;
Step 15: obtaining diagnosis rotor rupture conducting bar number according to formula (14);
Above-mentioned three-phase cage asynchronous motor rotor broken splits the diagnostic method of conducting bar number, the cage type asynchronous motor
Threephase stator voltage transient signal vmWith threephase stator electric current instantaneous signal imBe synchronous acquisition (m indicate separate, m=1,2,
3), the sample frequency of each signal is 1000Hz, when sampling a length of 10s, i.e. sampling number is 10000 points.
Wherein, the present invention is deduced Hilbert modulus 2sf1The amplitude of frequency component and rotor rupture conducting bar number it
Between numerical value corresponding relationship, and propose a kind of three-phase cage asynchronous motor rotor broken as diagnosis criterion and split conducting bar number
Quantitative Diagnosis method concrete foundation and derivation process are as follows:
It is well known that cage type asynchronous motor rotor strip-broken failure will lead to appearance (1-2s) f in its stator current1Frequency
Additional components iL', it is as follows:
Wherein, IL' indicate stator current (1-2s) f caused by rotor bar breaking fault1The virtual value of frequency component, αLTable
Show stator current caused by rotor bar breaking fault (1-2s) f1The initial phase angle of frequency component;T indicates the time, is variable, and e is certainly
Right constant, j are imaginary unit.
The additional components and fundamental flux effect will generate a pulsating torque Δ T, as follows:
Δ T=-3PI 'LΨsin[2π(2sf1)t+αL-αΨ] (16)
Wherein, P indicates number of pole-pairs;Ψ indicates the virtual value of fundamental flux, αΨIndicate the initial phase angle of fundamental flux.
The pulsating torque necessarily causes rotor speed with 2sf1Frequency fluctuation.This will form phase-modulation to fundamental flux,
And (1-2s) f is incuded in the stator windings1、(1+2s)f1The additional electromotive force of frequency.Wherein, frequency is (1-2s) f1Induction
Electromotive force generates (1-2s) f1The stator current additional components of frequency, can be regarded as to iL' reaction, be denoted as iL", show
In formula (17);And frequency is (1+2s) f1Induced electromotive force generate (1+2s) f1The stator current additional components of frequency, are remembered
For iR, it is shown in formula (18).
In formula (17), formula (18), IL" indicate stator current (1-2s) f caused by torque pulsation1Frequency component iL" have
Valid value, IRIndicate stator current (1+2s) f caused by torque pulsation1Frequency component iRVirtual value;f1For frequency of supply, αSTable
Show the impedance angle of stator impedance, ZSIndicate the modulus value of stator impedance;J indicates rotary inertia, and s is revolutional slip.
Therefore, in rotor bar breaking fault, the 1st phase stator current instantaneous signal i1It can indicate are as follows:
According to formula (19): in rotor bar breaking fault, practical stator current includes two (1-2s) f1Frequency
Component --- iL' and iL″.Therefore, FFT spectrum analysis is carried out to stator current signal and obtains (1-2s) f1 frequency component
Virtual value ILIt actually should be iL' and iL" resultant effect.From formula (15), formula (17): iL' and iL" phase difference be αS+
Pi/2, and the α for asynchronous motorS≈ pi/2, this means that iL' and iL" phase difference be approximately π, that is, the two is anti-
Phase.Therefore, IL=IL′-IL", that is:
IL'=IL+IL″ (20)
In formula (20), ILIt " can not be obtained and carrying out FFT spectrum analysis to stator current signal.But according to
Known to formula (17), formula (18): IL″≈IR(because s < < 1), and IRFor stator current (1+2s) f1The virtual value of frequency component is
It can be analyzed and to be obtained by stator current signal FFT spectrum.Therefore, formula (21) is set up:
I′L≈IL+IR (21)
Further it is necessary to note: formula (1) is to ignore before rotor speed fluctuation (i.e. hypothesis invariablenes turning speed or J=∞) is
It mentions, at this time IL"=0, IR=0, therefore IL′≈IL.This means that: I in formula (1)LReal meaning be actually IL′.It is based on
This because, and consider certainly exist rotor speed fluctuation, convolution (1), formula (21) can derive formula (2).It is above-mentioned to be
MCSA (Motor Current Signal Analysis) diagnoses the theoretical basis of rotor rupture conducting bar number according to formula (2).
Accordingly, it is considered to arrive αS≈π/2、IL″≈IR, formula (19) can simplify as formula (22) form:
According to formula (22), the 1st phase stator current instantaneous signal i can be derived1Hilbert transform, such as formula (23) institute
Show:
Wherein, τ is delay;Hilbert modulus Mod is further derived, can be obtained:
Wherein, ISFor f in stator current1The virtual value of frequency component;According to formula (24) it is found that the Hilbert modulus packet
Include DC component, 2s f1、4s f1Frequency component, it is extremely complex.In rotor bar breaking fault, cage type asynchronous motor is uncommon
To occur a 2sf in your Bert modulus1The wave component of frequency.Therefore, FFT spectrum analysis is done simultaneously to Hilbert modulus signal
Investigate the 2sf1The wave component of frequency can carry out rotor bar breaking fault detection.Here it is MHMSA (Motor Hilbert
Modulus Signal Analysis) method theoretical basis.As previously mentioned, in place of MCSA Shortcomings --- (1 ± 2s) f1
Component may be by f1The leakage of component is flooded.But in MHMSA method, the DC component in Hilbert modulus signal is to hold
Easily filter out, therefore above-mentioned " flooding " problem be not present, therefore this method is significantly better than MCSA, can it is sensitive, reliably detect
Rotor bar breaking fault.
But for now, MHMSA method not yet solves the problems, such as the diagnosis of rotor rupture conducting bar number.The present invention is then
It successfully solves the above problem, is deduced three-phase cage asynchronous motor Hilbert modulus 2sf1The amplitude of frequency component with
Numerical value corresponding relationship between rotor rupture conducting bar number, and a kind of three-phase cage asynchronous electric is proposed as diagnosis criterion
The quantitative Diagnosis method of machine rotor fracture conducting bar number.
For real electrical machinery, there are following numerical relations: IS>>IL、IS>>IR.Therefore, convolution (24) is it is found that the Xi Er
The main frequency content of Bert modulus is DC component and 2sf1Frequency component, the two are shown in formula (25), formula (26).
Ms=4IS(IL′-IR)cos[2π(2sf1)t+αL]+4ISIRcos[2π(2sf1)t+αL-2αΨ] (26)
Wherein, A0For DC component (i.e. average value) M of Hilbert modulus ModsIndicate Hilbert modulus Mod's
2sf1Frequency component.
For formula (26), Hilbert modulus 2sf is further derived1The expression formula of frequency component can obtain:
In formula (27), γ indicates a phase angle,The Hilbert modulus can be obtained according to formula (27)
2sf1The amplitude A of frequency component1, as shown in formula (28):
Formula (25) are substituted into formula (28), can be obtained:
In addition, for asynchronous motorTherefore formula (29) can indicate are as follows:
In formula (30),Indicate power-factor angle,As power factor.
It can be obtained according to formula (30):
It can be obtained according to formula (2), formula (21) and formula (31):
Wherein, NbFor rotor rupture conducting bar number, N is rotor bar sum;Above formula is Hilbert modulus 2sf1Frequency
Numerical value corresponding relationship between the amplitude and rotor rupture conducting bar number of component, in this, as diagnosis criterion i.e. it is proposed that for simultaneously
The quantitative Diagnosis method of rotor rupture conducting bar number suitable for MHMSA method.
Embodiment:
The present invention measures stator voltage, electric current instantaneous signal using circuit shown in Fig. 1, and the circuit is by voltage changer, electricity
Current-to-voltage converter, data acquisition card and portable computer composition.Three voltage changers are connected to asynchronous electric respectively
Machine stator three-phase windings, signal output end are connected to the 1st, 2,3 analog signal input channels of data acquisition card respectively.Described three
Current-to-voltage converter is connected to asynchronous motor stator three-phase windings respectively, and signal output end is connected to data acquisition card respectively
The 4th, 5,6 analog signal input channels.The output port of the data acquisition card connects the USB port of portable computer.Signal is adopted
Truck uses auspicious rich China's RBH8351 type data acquisition card, and the model of portable computer is Thinkpad X100e.Data acquisition card
It is integrated with the circuits such as low-pass filter, signal acquisition holding, analog/digital conversion.Stator voltage, electric current instantaneous signal send to signal and adopt
Truck, data acquisition card are connected to portable computer by USB interface.Portable computer controls data acquisition card with appropriate frequency
Stator voltage, electric current instantaneous signal are sampled, and is stored in hard disk, then stator voltage, current signal are carried out by portable computer
Processing, and rotor rupture conducting bar number is diagnosed, steps are as follows:
A. threephase stator voltage transient signal v is measuredmWith threephase stator electric current instantaneous signal im(m indicate separate, m=1,2,
3)。
Three no-load voltage ratios are used to measure threephase stator voltage transient signal for the voltage changer of 220V/2.5V.For high pressure
Motor measures threephase stator voltage transient signal in voltage transformer pt secondary side;For low voltage motor, directly in motor wiring
Threephase stator voltage transient signal is measured at terminal.
Three no-load voltage ratios are used to measure threephase stator electric current instantaneous signal for the current-to-voltage converter of 10A/1.0V.For big
Middle size motor measures threephase stator electric current instantaneous signal in Current Transmit secondary side;For micro-machine, directly in motor
Threephase stator electric current instantaneous signal is measured at connecting terminal.
Above-mentioned threephase stator voltage transient signal and threephase stator electric current instantaneous signal are that synchronized sampling measures, Mei Yixin
Number sample frequency be 1000Hz, sampling when a length of 10s, i.e. sampling number is 10000 points.
B. m phase stator voltage v is calculatedmVirtual valueWherein vm(k) it represents its k-th
Sampled value.
C. the arithmetic mean of instantaneous value of threephase stator voltage effective value is calculatedAnd using it as stator voltage virtual value
VRMS。
D. m phase stator current i is calculatedmVirtual valueWherein im(k) it is fixed to represent m phase
K-th of sampled value of electron current.
E. the arithmetic mean of instantaneous value of threephase stator current effective value is calculatedAnd using it as stator current virtual value
IRMS。
F. instantaneous active power p is calculated:
G. DC component (that is, average value) mean (p), mean () for calculating instantaneous active power p indicate fortune of averaging
It calculates.
H. power factor is calculatedAnd it is applied to formula (32).
I. to the 1st phase stator current instantaneous signal i1FFT spectrum analysis is done, its f is obtained1The virtual value I of frequency componentS, and
It is applied to formula (32).
J. the 1st phase stator current instantaneous signal i is calculated1Hilbert transform
K. Hilbert modulus Mod is calculated,
L. DC component (that is, average value) A of Hilbert modulus Mod is calculated0=mean (Mod), and be applied to
Formula (32).
M. the DC component A in Hilbert modulus Mod is filtered out0, obtain signal M, M=Mod-A to be analyzed0。
N. FFT spectrum analysis is carried out to signal M to be analyzed, obtains its 2sf1The amplitude A of frequency component1, and be applied to
Formula (32).
O. rotor rupture conducting bar number is diagnosed according to formula (32), it may be assumed that
Wherein, NbFor rotor rupture conducting bar number;N is rotor bar sum.
Rotor rupture conducting bar is carried out to a Y100L-2 type three-phase cage asynchronous motor (3kW, 380V) using this method
Number diagnosis, effect are satisfactory.The rotor bar sum of the motor is N=20.
Fig. 2, Fig. 3 respectively indicate Hilbert modulus FFT spectrum of the motor full load in 1 broken strip of rotor with
Stator current FFT spectrum, specific data are shown in table 1, table 2.
According to fig. 2 and combine table 1, it is known that: in 1 broken strip of rotor, diagnostic method of the present invention for MHMSA
Diagnostic result is 0.80 ≈, 1 broken strip, and with the substantially disorderly conjunction of actual conditions, this diagnostic result is effective.
According to Fig. 3 and combine table 2, it is known that: in 1 broken strip of rotor, existing diagnostic method for MCSA is examined
Disconnected result is 0.72 ≈, 1 broken strip.
Contrast table 1, table 2, it is known that: the diagnostic result of the diagnostic method of the present invention for MHMSA with it is existing for MCSA
The diagnostic result of diagnostic method is approximate consistent.
This indicates that three-phase cage asynchronous motor rotor broken proposed by the invention splits the quantitative Diagnosis side of conducting bar number
Method be it is effective, based on diagnosis criterion --- three-phase cage asynchronous motor Hilbert modulus 2sf1Frequency component
Numerical value corresponding relationship between amplitude and rotor rupture conducting bar number, i.e. formula (32) are also correct.
In summary, there are offices for the MHMSA method of current three-phase cage asynchronous motor rotor strip-broken fault detection
Limit --- the number of rotor rupture conducting bar cannot be quantitatively determined.The present invention is deduced three-phase cage asynchronous motor Martin Hilb
Special modulus 2sf1Numerical value corresponding relationship between the amplitude and rotor rupture conducting bar number of frequency component, and as diagnosis criterion
Propose a kind of quantitative Diagnosis method that three-phase cage asynchronous motor rotor broken splits conducting bar number.This method is directed to and is suitable for
The MHMSA method of three-phase cage asynchronous motor rotor strip-broken fault detection compensates for it and is unable to quantitative Diagnosis rotor rupture and leads
The defect of number, has engineering value and application prospect.
1 rotor rupture conducting bar number diagnostic data of table (the method for the present invention is used for MHMSA)
2 rotor rupture conducting bar number diagnostic data of table (existing method is used for MCSA)
It is broken conducting bar number | IS(A) | IL(A) | IR(A) | Diagnostic result | Diagnosis basis |
1 | 5.56 | 0.0952 | 0.1057 | 0.72≈1 | Formula (2), formula (9) |
In conclusion the above is merely preferred embodiments of the present invention, being not intended to limit the scope of the present invention.
All within the spirits and principles of the present invention, any modification, equivalent replacement, improvement and so on should be included in of the invention
Within protection scope.
Claims (3)
1. the diagnostic method that a kind of three-phase cage asynchronous motor rotor broken splits conducting bar number, which is characterized in that specifically include as
Lower step:
Step 1: the sample frequency acquisition three-phase cage asynchronous motor threephase stator voltage transient signal and three-phase by setting are fixed
Electron current instantaneous signal;
Step 2: obtain m phase stator voltage virtual value, m indicate separate, m=1,2,3;
Step 3: obtaining the arithmetic mean of instantaneous value of threephase stator voltage effective value and as stator voltage virtual value;
Step 4: obtaining the virtual value of m phase stator current;
Step 5: obtaining the arithmetic mean of instantaneous value of threephase stator current effective value and as stator current virtual value;
Step 6: obtaining instantaneous active power according to threephase stator voltage transient signal and threephase stator electric current instantaneous signal;
Step 7: the method by calculating average value, obtains the DC component mean (p) of instantaneous active power p, i.e. wattful power
Rate;
Step 8: obtaining apparent energy according to stator voltage virtual value and stator current virtual value;Active power and apparent energy
The ratio between, as power factor
Step 9: doing FFT spectrum analysis to the electric current instantaneous signal of any one phase stator in threephase stator, stator current is obtained
Middle f1The virtual value I of frequency componentS;Wherein, f1For frequency of supply;
Step 10: carrying out Hilbert transform to electric current instantaneous signal;
Step 11: obtaining Hilbert modulus Mod;
Step 12: obtaining the DC component of Hilbert modulus Mod, i.e. average value A0;
Step 13: filtering out the DC component in Hilbert modulus Mod, signal M to be analyzed is obtained;
Step 14: carrying out FFT spectrum analysis to signal M to be analyzed, its 2sf is obtained1The amplitude A of frequency component1;Wherein, s is
Revolutional slip;
Step 15: according toObtain diagnosis rotor rupture conducting bar
Number Nb, wherein N is rotor bar sum.
2. a kind of three-phase cage asynchronous motor rotor broken as described in claim 1 splits the diagnostic method of conducting bar number, special
Sign is that the threephase stator voltage transient signal of the three-phase cage asynchronous motor is synchronous with threephase stator electric current instantaneous signal
Acquisition.
3. a kind of three-phase cage asynchronous motor rotor broken as described in claim 1 splits the diagnostic method of conducting bar number, special
Sign is, the sample frequency of each signal is 1000Hz, when sampling a length of 10s, i.e. sampling number is 10000 points.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610843046.9A CN106646222B (en) | 2016-09-22 | 2016-09-22 | A kind of three-phase cage asynchronous motor rotor broken splits the diagnostic method of conducting bar number |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610843046.9A CN106646222B (en) | 2016-09-22 | 2016-09-22 | A kind of three-phase cage asynchronous motor rotor broken splits the diagnostic method of conducting bar number |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106646222A CN106646222A (en) | 2017-05-10 |
CN106646222B true CN106646222B (en) | 2019-06-21 |
Family
ID=58853013
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610843046.9A Active CN106646222B (en) | 2016-09-22 | 2016-09-22 | A kind of three-phase cage asynchronous motor rotor broken splits the diagnostic method of conducting bar number |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106646222B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113947125A (en) * | 2021-07-05 | 2022-01-18 | 华北电力大学(保定) | High-precision diagnosis method for rotor broken number of fifteen-phase asynchronous motor based on ESPRIT-PSA and LGBM |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102944842A (en) * | 2012-11-30 | 2013-02-27 | 华北电力大学(保定) | Detecting method for rotor broken bar fault of cage-type asynchronous motor |
CN103018669A (en) * | 2012-11-30 | 2013-04-03 | 华北电力大学(保定) | Method for detecting failure of rotor broken bar of cage type asynchronous motor at high reliability |
CN103235260A (en) * | 2013-04-15 | 2013-08-07 | 哈尔滨工业大学 | Submersible motor rotor broken bar fault recognition method based on HHT (Hilbert-Huang transform) |
US9433365B1 (en) * | 2012-07-26 | 2016-09-06 | National Yang-Ming University | System and method for indentifying rotors in fractionated signals in persistent atrial fibrillation ablation |
-
2016
- 2016-09-22 CN CN201610843046.9A patent/CN106646222B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9433365B1 (en) * | 2012-07-26 | 2016-09-06 | National Yang-Ming University | System and method for indentifying rotors in fractionated signals in persistent atrial fibrillation ablation |
CN102944842A (en) * | 2012-11-30 | 2013-02-27 | 华北电力大学(保定) | Detecting method for rotor broken bar fault of cage-type asynchronous motor |
CN103018669A (en) * | 2012-11-30 | 2013-04-03 | 华北电力大学(保定) | Method for detecting failure of rotor broken bar of cage type asynchronous motor at high reliability |
CN103235260A (en) * | 2013-04-15 | 2013-08-07 | 哈尔滨工业大学 | Submersible motor rotor broken bar fault recognition method based on HHT (Hilbert-Huang transform) |
Non-Patent Citations (3)
Title |
---|
A New Criterion for Rotor Broken Bar Fault Diagnosis in Line-start and Inverter-fed Induction Motors using Hilbert-Huang Transform;Jawad Faiz 等;《2012 IEEE International Conference on Power Electronics, Drives and Energy Systems》;20121219;1-6 |
基于Hilbert模量频谱分析的异步电机转子断条故障研究;马宏忠 等;《电机与控制学报》;20090531;第13卷(第3期);371-376 |
基于复调制滤波与ESPRIT的异步电动机转子断条故障检测新方法;许伯强 等;《电测与仪表》;20140610;第51卷(第11期);76-80 |
Also Published As
Publication number | Publication date |
---|---|
CN106646222A (en) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103454585B (en) | A kind of permagnetic synchronous motor loss of excitation method for diagnosing faults based on dead electricity residual voltage | |
CN106199424B (en) | Permanent magnet synchronous motor turn-to-turn short circuit fault diagnosis method | |
CN109991539A (en) | Multiphase corner connection brushless exciter rotating diode one-phase open circuit detection method and system | |
CN101025430A (en) | Cage type asynchronous motor rotor strip-broken failure detecting method | |
CN108196154B (en) | Fault detection and fault positioning method for rotary rectifier of aviation three-stage synchronous motor | |
CN103018669B (en) | Method for detecting failure of rotor broken bar of cage type asynchronous motor at high reliability | |
CN109738780A (en) | One tube open circuit detection method of multiphase corner connection brushless exciter rotating diode and system | |
Kim et al. | High-frequency signal injection-based rotor bar fault detection of inverter-fed induction motors with closed rotor slots | |
CN106646223B (en) | A kind of three-phase cage asynchronous motor rotor broken splits the diagnostic method of conducting bar number | |
Alwodai et al. | Modulation signal bispectrum analysis of motor current signals for stator fault diagnosis | |
Wang et al. | Induction motor broken rotor bar fault diagnosis based on third-order energy operator demodulated current signal | |
Grantham et al. | Rapid parameter determination for induction motor analysis and control | |
CN107017808B (en) | The continuous estimation method of synchronous motor rotor position based on pulsation exciting current response | |
CN102081129B (en) | Method for diagnosing short circuit fault of stator winding of switched reluctance motor | |
CN106646222B (en) | A kind of three-phase cage asynchronous motor rotor broken splits the diagnostic method of conducting bar number | |
CN106443446B (en) | A kind of method for diagnosing faults suitable for permanent magnet machine rotor | |
CN105607000A (en) | Method for measuring locked-rotor current of three-phase asynchronous motor under no-load condition | |
CN106646224A (en) | Three-phase asynchronous motor load fluctuation detection method | |
Pietrzak et al. | Stator phase current STFT analysis for the PMSM stator winding fault diagnosis | |
Zhu et al. | Practical aspects of broken rotor bars detection in PWM voltage-source-inverter-fed squirrel-cage induction motors | |
Drif et al. | The instantaneous power factor approach for rotor cage faults diagnosis in three-phase induction motors | |
Halder et al. | Stator inter-turn fault diagnosis in permanent magnet synchronous motor | |
Ayhan et al. | Statistical analysis on a case study of load effect on PSD technique for induction motor broken rotor bar fault detection | |
Babaa et al. | Combined electrical faults detection and diagnosis using current signature analysis | |
Liu et al. | Stator inter-turn fault detection for the converter-fed induction motor based on the adjacent-current phase-shift |
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 |