CN106646222A - Method for diagnosing number of broken conducting bars of rotor of three-phase squirrel-cage asynchronous motor - Google Patents

Method for diagnosing number of broken conducting bars of rotor of three-phase squirrel-cage asynchronous motor Download PDF

Info

Publication number
CN106646222A
CN106646222A CN201610843046.9A CN201610843046A CN106646222A CN 106646222 A CN106646222 A CN 106646222A CN 201610843046 A CN201610843046 A CN 201610843046A CN 106646222 A CN106646222 A CN 106646222A
Authority
CN
China
Prior art keywords
rotor
stator
signal
phase
asynchronous motor
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.)
Granted
Application number
CN201610843046.9A
Other languages
Chinese (zh)
Other versions
CN106646222B (en
Inventor
许伯强
关涛
孙丽玲
王艳武
回志澎
寿海明
冀欣
冀路明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
PLA NAVY EQUIPMENT RESEARCH INSTITUTE SHIP ARGUMENT INSTITUTE
North China Electric Power University
Original Assignee
PLA NAVY EQUIPMENT RESEARCH INSTITUTE SHIP ARGUMENT INSTITUTE
North China Electric Power University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by PLA NAVY EQUIPMENT RESEARCH INSTITUTE SHIP ARGUMENT INSTITUTE, North China Electric Power University filed Critical PLA NAVY EQUIPMENT RESEARCH INSTITUTE SHIP ARGUMENT INSTITUTE
Priority to CN201610843046.9A priority Critical patent/CN106646222B/en
Publication of CN106646222A publication Critical patent/CN106646222A/en
Application granted granted Critical
Publication of CN106646222B publication Critical patent/CN106646222B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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 a method for diagnosing the number of broken conducting bars of the rotor of a three-phase squirrel-cage asynchronous motor. The method comprises the following steps: getting the stator voltage effective value, stator current effective value, active power, apparent power and power factor according to three-phase stator voltage and current transient signals of a three-phase squirrel-cage asynchronous motor; making an FFT spectrum analysis of the current transient signal to get the effective value IS of the frequency component f1 in the stator current; carrying out Hilbert transform on a current transient signal i1, filtering the DC component of the Hilbert modulus, and making an FFT spectrum analysis to get the amplitude of the frequency component 2sf1 of the to-be-analyzed signal; and getting the number of broken conducting bars of the rotor through diagnosis based on the derived corresponding relationship between the amplitude of the frequency component 2sf1 of the Hilbert modulus of the squirrel-cage asynchronous motor and the number of broken conducting bars of the rotor. The method makes up for the defect of the MHMSA method which cannot diagnose the number of broken conducting bars of a rotor quantitatively, and creates favorable conditions for discovering and eliminating a rotor bar breakage fault in time and ensuring secure and stable operation of motors.

Description

A kind of three-phase cage asynchronous motor rotor broken splits the diagnostic method of sliver number
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 sliver number Diagnostic method.
Background technology
In three-phase cage asynchronous motor running, rotor bar is subject to radial electromagnetic force, electric rotating magnetic force, centrifugation The effect of the alternate stresses such as power, thermal flexure amount of deflection power, in addition rotor manufacturing defect, may cause broken bar fault, this kind of failure to be sent out Raw probability is about 10%.
Rotor broken bar is typical gradual failure, initial stage usual 1 bar failure, then neighbouring with fracture sliver its His sliver continues to rupture, and finally causes motor output drop even to be shut down.Therefore, it is necessary to implement rotor bar breaking fault detection and and When to rupture sliver number make diagnosis.
Cage type asynchronous motor occurs after rotor bar breaking fault, and (1 ± 2s) f is will appear from its stator current1Frequency Extra current component (s is revolutional slip, f1For frequency of supply), can be as rotor bar breaking fault feature.And stator is electric Stream signal is easy to collection, therefore the stator current signal frequency spectrum analysis method based on Fast Fourier Transform (FFT) (FFT) extensively should For rotor bar breaking fault detection.This kind of method is commonly referred to 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 whether there is broken strip, 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 sliver number, as shown in formula (1).
Wherein, Nb, N be respectively rotor rupture sliver number and sliver sum.
But, formula (1) be to ignore rotor speed fluctuation (assuming invariablenes turning speed) premised on, and rotor bar breaking fault Rotor speed is necessarily caused to fluctuate, therefore formula (1) is that perfect condition does not tally with the actual situation, and rotor rupture is diagnosed accordingly and is led Bar number is short of accuracy.
For this purpose, MCSA transfers to diagnose rotor rupture sliver 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 more conform to actual conditions.
Obviously, during MCSA diagnosis rotor rupture sliver number, no matter using formula (1) or formula (1), 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 during broken strip slight due to rotor, (1 ± 2s) f1The virtual value of component is relative to f1Component is very little, and asynchronous motor is transported Revolutional slip s very littles are (in general, s during row<0.05), (1 ± 2s) f1With f1Numerical value is close to, 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 IRCannot know, then formula (1) cannot be applied with formula (2), diagnosis rotor rupture sliver number is not known where to begin certainly.This It is the weak point of MCSA.
Stator current (1 ± 2s) f in overcome MCSA methods1Component is easily by f1The defect that the leakage of component is flooded, sends out Exhibition is defined based on Hilbert modulus signal analysis (Motor Hilbert Modulus Signal Analysis, MHMSA) Cage type asynchronous motor rotor strip-broken failure detecting method, its core is：In cage type asynchronous motor rotor normal condition Under, its Hilbert modulus is only a DC component；Once generation rotor bar breaking fault, will appear from 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 Analyse and investigate the 2sf1The wave component of frequency, you can realize the detection of rotor bar breaking fault.Great mass of data shows, party's Faxian Write be better than MCSA, can with it is sensitive, reliably detect rotor bar breaking fault.
But, for now, MHMSA methods not yet solve the problems, such as the diagnosis of rotor rupture sliver number.Specifically, Existing MHMSA methods can not quantitatively determine the number of rotor rupture sliver, and reason is：Not yet find that three-phase cage is asynchronous Motor Hilbert modulus 2sf1Numerical value corresponding relation between the amplitude and rotor rupture sliver number of frequency component.
Therefore, in using MHMSA procedures, how to diagnose three-phase cage asynchronous motor rotor broken and split sliver number Mesh, just becomes the difficult problem faced at present about technical staff.
The content of the invention
In view of this, the invention provides a kind of three-phase cage asynchronous motor rotor broken splits the diagnosis side of sliver number Method, can in time find and exclude 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 sliver number, specifically includes following steps：
Step one, by setting sample frequency collection three-phase cage asynchronous motor threephase stator voltage transient signal and three Phase stator current instantaneous signal；
Step 2, the virtual value for obtaining m phase stator voltages, m represents separate, m=1,2,3；
Step 3, obtain threephase stator voltage effective value arithmetic mean of instantaneous value and as stator voltage virtual value；
Step 4, the virtual value for obtaining m phase stator currents；
Step 5, obtain threephase stator current effective value arithmetic mean of instantaneous value and as stator current virtual value；
Step 6, according to threephase stator voltage transient signal and threephase stator electric current instantaneous signal, obtain instantaneous active work( Rate；
Step 7, the method by calculating mean value, obtain DC component mean (p) of instantaneous active power p, i.e., active Power；
Step 8, according to stator voltage virtual value and stator current virtual value, obtain apparent energy；Active power and apparent Power ratio, as power factor
Step 9, to threephase stator in any one phase stator electric current instantaneous signal do FFT spectrum analysis, obtain stator F in electric current1The virtual value I of frequency componentS；Wherein, f1For frequency of supply；
Step 10, Hilbert transform is carried out to electric current instantaneous signal；
Step 11, acquisition Hilbert modulus Mod；
Step 12, the DC component for obtaining Hilbert modulus Mod, i.e. mean value A0；
Step 13, DC component A filtered in Hilbert modulus Mod0, obtain signal M to be analyzed；
Step 14, treat signal Analysis M and carry out FFT spectrum analysis, obtain its 2sf1Amplitude A of frequency component1；Wherein, S is revolutional slip；
Step 15, basisObtain diagnosis rotor rupture Sliver number Nb, wherein, N is rotor bar sum.
It is preferred that the threephase stator voltage transient signal of the three-phase cage asynchronous motor is instantaneous with threephase stator electric current Signal is synchronous acquisition.
It is preferred that the sample frequency of each signal is 1000Hz, a length of 10s, i.e. sampling number are 10000 points during sampling.
Beneficial effect：
The present invention is according to three-phase cage asynchronous motor Hilbert modulus 2sf1The amplitude of frequency component judges that rotor breaks Bar number, compensate for the defect that traditional MHMSA methods are unable to quantitative Diagnosis rotor rupture sliver number, be to find in time and exclude Rotor bar breaking fault, it is ensured that the safe and stable operation of motor creates advantage.
Description of the drawings
Fig. 1 is present invention collection three-phase cage asynchronous motor stator voltage, the electrical schematic diagram of 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
Develop simultaneously below in conjunction with the accompanying drawings embodiment, describes the present invention.
The invention provides a kind of three-phase cage asynchronous motor rotor broken splits the diagnostic method of sliver number, wherein, respectively Meaning of parameters is：
S is revolutional slip；f1For frequency of supply；FFT is Fast Fourier Transform (FFT)；The English full name 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 sliver number；N is rotor bar sum；IRFor (1+2s) f in stator current1 The virtual value of frequency component；The English full name of MHMSA is Motor Hilbert Modulus Signal Analysis；vmFor M phase stator voltage instantaneous signals；imFor m phase stator current instantaneous signals；M be separate, m=1,2,3；VmFor m phase stators Voltage effective value；vmK () is k-th sampled value of m phase stator voltage instantaneous signals；N is sampling number；VRMSFor stator electricity It is pressed with valid value；ImFor m phase stator current virtual values；imK () is k-th sampled value of m phase stator current instantaneous signals；IRMS For stator current virtual value；P is instantaneous active power；For power factor；For power-factor angle；For the 1st phase stator electricity Stream instantaneous signal i1Hilbert transform；T is the time；τ is time delay；Mod is Hilbert modulus；A0For Hilbert modulus The mean value of Mod, that is, DC component；Mean () represents computing of averaging；M filters direct current point for Hilbert modulus Mod Amount A0Signal to be analyzed afterwards；A1Represent the 2sf of Hilbert modulus Mod1The amplitude of frequency component；iL' represent rotor broken bar event Stator current (1-2s) f caused by barrier institute1Frequency component；IL' for rotor bar breaking fault institute caused by stator current (1-2s) f1Frequently The virtual value of rate component；αLFor stator current (1-2s) f caused by rotor bar breaking fault institute1The initial phase angle of frequency component；E is certainly Right constant；J is imaginary unit；Δ T is pulsating torque；P is number of pole-pairs；Ψ is the virtual value of fundamental flux；αΨFor fundamental flux Initial phase angle；iL" for stator current (1-2s) f caused by torque pulsation institute1Frequency component；iRTo determine caused by torque pulsation institute Electron current (1+2s) f1Frequency component；IL" for stator current (1-2s) f caused by torque pulsation institute1The virtual value of frequency component； IRFor stator current (1+2s) f caused by torque pulsation institute1The virtual value of frequency component；ZSFor the modulus value of stator impedance；αSIt is fixed Sub- impedance angle；J, rotary inertia；MsFor the 2sf of Hilbert modulus Mod1Frequency component；γ is phase angle；PT is mutual induction of voltage Device；CT is current transformer.
Step one, by certain frequency gather three-phase cage asynchronous motor threephase stator voltage transient signal vmIt is fixed with three-phase Electron current instantaneous signal im, wherein, m represents separate, m=1,2,3；
Step 2, according to formula (3), obtain m phase stator voltages vmVirtual value Vm；
Wherein vmK () represents k-th sampled value of m phase stator voltages, n is sampling number；
Step 3, according to formula (4), obtain the arithmetic mean of instantaneous value of threephase stator voltage effective value, have as stator voltage Valid value VRMS；
Step 4, according to formula (5), obtain m phase stator current imVirtual value Im；
Wherein imK () represents k-th sampled value of m phase stator currents, n is sampling number；
Step 5, according to formula (6), obtain the arithmetic mean of instantaneous value of threephase stator current effective value, have 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, using formula (8), obtain DC component mean (p) of instantaneous active power p, i.e. active power；
Step 8, according to formula (9), obtain power factor
Wherein, 3VRMSIRMSFor apparent energy.
Step 9, to threephase stator in any one phase stator electric current instantaneous signal i1FFT 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 time delay；
Step 11, according to formula (11), obtain Hilbert modulus Mod：
Step 12, DC component (i.e. mean value) A for obtaining Hilbert modulus Mod0；
Step 13, DC component A filtered in Hilbert modulus Mod0, obtain signal M to be analyzed；
M=Mod-A0 (13)
Step 14, treat signal Analysis M and carry out FFT spectrum analysis, obtain its 2s f1Amplitude A of frequency component1；Its In, s is revolutional slip；
Step 15, according to formula (14), obtain diagnosis rotor rupture sliver number；
Above-mentioned three-phase cage asynchronous motor rotor broken splits the diagnostic method of sliver number, the cage type asynchronous motor Threephase stator voltage transient signal vmWith threephase stator electric current instantaneous signal imBe synchronous acquisition (m represents separate, m=1,2, 3), the sample frequency of each signal is 1000Hz, and a length of 10s, i.e. sampling number are 10000 points during sampling.
Wherein, the present invention is deduced Hilbert modulus 2sf1The amplitude of frequency component and rotor rupture sliver number it Between numerical value corresponding relation, and propose a kind of three-phase cage asynchronous motor rotor broken as diagnosis criterion and split sliver number Quantitative Diagnosis method concrete foundation and derivation be：
It is well known that cage type asynchronous motor rotor strip-broken failure will cause occur (1-2s) f in its stator current1Frequency Additional components iL', it is as follows：
Wherein, IL' represent stator current (1-2s) f caused by rotor bar breaking fault institute1The virtual value of frequency component, αLTable Show stator current (1-2s) f caused by rotor bar breaking fault institute1The initial phase angle of frequency component；T represents the time, is variable, and e is certainly Right constant, j is imaginary unit.
The additional components will produce a pulsating torque Δ T with fundamental flux effect, as follows：
Δ T=-3PI 'LΨsin[2π(2sf1)t+αL-αΨ] (16)
Wherein, P represents number of pole-pairs；Ψ represents the virtual value of fundamental flux, αΨRepresent 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 sensed in the stator windings1、(1+2s)f1The additional electromotive force of frequency.Wherein, frequency is (1-2s) f1Sensing Electromotive force produces (1-2s) f1The stator current additional components of frequency, can be regarded as to iL' reaction, be designated as iL", show In formula (17)；And frequency is (1+2s) f1Induced electromotive force produce (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" represent stator current (1-2s) f caused by torque pulsation institute1Frequency component iL" have Valid value, IRRepresent stator current (1+2s) f caused by torque pulsation institute1Frequency component iRVirtual value；f1For frequency of supply, αSTable Show the impedance angle of stator impedance, ZSRepresent the modulus value of stator impedance；J represents rotary inertia, and s is revolutional slip.
Therefore, in the case of rotor bar breaking fault, the 1st phase stator current instantaneous signal i1Can be expressed as：
Understood according to formula (19)：In the case of rotor bar breaking fault, stator current is actual to include 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 components Virtual value ILActually 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/2s, 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), IL" cannot be obtained by carrying out FFT spectrum analysis to stator current signal.But according to Formula (17), formula (18) understand：IL″≈IR(because s<<1), IRFor stator current (1+2s) f1The virtual value of frequency component, be Can be analyzed by stator current signal FFT spectrum and be obtained.Therefore, formula (21) is set up：
I′L≈IL+IR (21)
Further it is necessary to note：Formula (1) is to ignore rotor speed fluctuation (assuming invariablenes turning speed, or J=∞) for front Carry, now IL"=0, IR=0, therefore IL′≈IL.This means：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 foundation of rotor rupture sliver number according to formula (2).
Accordingly, it is considered to arrive αS≈π/2、IL″≈IR, formula (19) can be reduced to 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 time delay；Hilbert modulus Mod is further derived, can be obtained：
Wherein, ISFor f in stator current1The virtual value of frequency component；Understood according to formula (24), the Hilbert modulus bag Include DC component, 2s f1、4s f1Frequency component, it is extremely complex.In the case of rotor bar breaking fault, cage type asynchronous motor is wished A 2sf is will appear from your Bert modulus1The wave component of frequency.Therefore, FFT spectrum is done to Hilbert modulus signal to analyze simultaneously Investigate the 2sf1The wave component of frequency, you can carry out rotor bar breaking fault detection.Here it is MHMSA (Motor Hilbert Modulus Signal Analysis) method theoretical foundation.As it was previously stated, MCSA Shortcomings part --- (1 ± 2s) f1 Component may be by f1The leakage of component is flooded.But in MHMSA methods, the DC component in Hilbert modulus signal is to hold Easily filter, therefore not there are problems that above-mentioned " flooding ", therefore the method is significantly better than MCSA, can with it is sensitive, reliably detect Rotor bar breaking fault.
But, for now, MHMSA methods not yet solve the problems, such as the diagnosis of rotor rupture sliver number.It is of the invention then The problems referred to above are successfully solved, three-phase cage asynchronous motor Hilbert modulus 2sf is deduced1The amplitude of frequency component with Numerical value corresponding relation between rotor rupture sliver number, and propose a kind of three-phase cage asynchronous electric as diagnosis criterion The quantitative Diagnosis method of machine rotor fracture sliver number.
For real electrical machinery, there is following numerical relation：IS>>IL、IS>>IR.Therefore, knowable to convolution (24), the Xi Er The main frequency content of Bert modulus is DC component and 2sf1Frequency component, the two is 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. mean value) M of Hilbert modulus ModsRepresent Hilbert modulus Mod 2sf1Frequency component.
For formula (26), Hilbert modulus 2sf is further derived1The expression formula of frequency component, can obtain：
In formula (27), γ represents a phase angle, Can Hilbert modulus 2sf according to formula (27)1Amplitude A of frequency component1, as shown in formula (28)：
Formula (25) is substituted into formula (28), can be obtained：
In addition, for asynchronous motorTherefore formula (29) can be expressed as：
In formula (30),Represent power-factor angle,As power factor.
Can be obtained according to formula (30)：
Can be obtained according to formula (2), formula (21) and formula (31)：
Wherein, NbFor rotor rupture sliver number, N is rotor bar sum；Above formula is Hilbert modulus 2sf1Frequency Numerical value corresponding relation between the amplitude and rotor rupture sliver number of component, can be proposed for simultaneously in this, as diagnosis criterion Suitable for the quantitative Diagnosis method of the rotor rupture sliver number of MHMSA methods.
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 pocket computer composition.Three voltage changers are connected to respectively asynchronous electric Machine stator three-phase windings, its signal output part is connected to respectively the 1st, 2,3 analog signal input channels of data acquisition card.Described three Current-to-voltage converter is connected to respectively asynchronous motor stator three-phase windings, and its signal output part is connected to respectively data acquisition card The 4th, 5,6 analog signal input channels.The output port of the data acquisition card connects the USB port of pocket computer.Signal is adopted Truck adopts auspicious rich magnificent RBH8351 types data acquisition card, and the model of pocket computer is Thinkpad X100e.Data acquisition card It is integrated with the circuits such as low pass filter, signals collecting holding, analog/digital conversion.Stator voltage, electric current instantaneous signal are delivered to signal and are adopted Truck, data acquisition card is connected to pocket computer by USB interface.Pocket computer control signal capture card is with appropriate frequency Sampling stator voltage, electric current instantaneous signal, and hard disk is stored in, then stator voltage, current signal are carried out by pocket computer Process, and diagnose rotor rupture sliver number, step is as follows：
A. threephase stator voltage transient signal v is measuredmWith threephase stator electric current instantaneous signal im(m represents separate, m=1,2, 3)。
Adopt three no-load voltage ratios that threephase stator voltage transient signal is measured for the voltage changer of 220V/2.5V.For high pressure Motor, in voltage transformer pt secondary side threephase stator voltage transient signal is measured；For low voltage motor, directly in motor wiring Threephase stator voltage transient signal is measured at terminal.
Adopt three no-load voltage ratios that threephase stator electric current instantaneous signal is measured for the current-to-voltage converter of 10A/1.0V.For big Middle size motor, in Current Transmit secondary side threephase stator electric current instantaneous signal is measured；For micro-machine, directly in motor Threephase stator electric current instantaneous signal is measured at binding post.
Above-mentioned threephase stator voltage transient signal and threephase stator electric current instantaneous signal are that synchronized sampling is measured, Mei Yixin Number sample frequency be 1000Hz, sampling when a length of 10s, i.e. sampling number be 10000 points.
B. m phase stator voltages v are calculatedmVirtual valueWherein vmK () 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 are calculatedmVirtual valueWherein imK () represents m phases fixed K-th 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, mean value) mean (p) of instantaneous active power p is calculated, mean () represents fortune of averaging Calculate.
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, mean value) A of Hilbert modulus Mod is calculated0=mean (Mod), and be applied to Formula (32).
M. DC component A in Hilbert modulus Mod is filtered0, obtain signal M, M=Mod-A to be analyzed0。
N. treating signal Analysis M carries out FFT spectrum analysis, obtains its 2sf1Amplitude A of frequency component1, and be applied to Formula (32).
O. rotor rupture sliver number is diagnosed according to formula (32), i.e.,：
Wherein, NbFor rotor rupture sliver number；N is rotor bar sum.
Rotor rupture sliver is carried out to Y100L-2 types three-phase cage asynchronous motor (3kW, 380V) using the method Number is diagnosed, and effect is satisfactory.The rotor bar sum of the motor is N=20.
Fig. 2, Fig. 3 represent respectively Hilbert modulus FFT spectrum of the motor full load in the case of 1 broken strip of rotor with Stator current FFT spectrum, concrete data are shown in table 1, table 2.
According to Fig. 2 and with reference to table 1, it is known that：In the case of 1 broken strip of rotor, for the diagnostic method of the present invention of 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 with reference to table 2, it is known that：In the case of 1 broken strip of rotor, the existing diagnostic method for MCSA is examined Disconnected result is 0.72 ≈, 1 broken strip.
Contrast table 1, table 2, it is known that：For MHMSA diagnostic method of the present invention diagnostic result with it is existing for MCSA The diagnostic result of diagnostic method is approximate consistent.
This is indicated that：Three-phase cage asynchronous motor rotor broken proposed by the invention splits the quantitative Diagnosis side of sliver number Method is effective, the diagnosis criterion of its foundation --- three-phase cage asynchronous motor Hilbert modulus 2sf1Frequency component Numerical value corresponding relation between amplitude and rotor rupture sliver number, i.e. formula (32) are also correct.
In summary, the MHMSA method presence office of current three-phase cage asynchronous motor rotor strip-broken fault detect Limit --- can not quantitatively determine the number of rotor rupture sliver.The present invention is deduced three-phase cage asynchronous motor Martin Hilb Special modulus 2sf1Numerical value corresponding relation between the amplitude and rotor rupture sliver number of frequency component, and as diagnosis criterion Propose a kind of quantitative Diagnosis method that three-phase cage asynchronous motor rotor broken splits sliver number.The method is directed to and is applied to The MHMSA methods of three-phase cage asynchronous motor rotor strip-broken fault detect, compensate for it and are unable to quantitative Diagnosis rotor rupture leading The defect of bar number, possesses construction value and application prospect.
Rotor rupture sliver number diagnostic data (the inventive method, for MHMSA) of table 1
Rotor rupture sliver number diagnostic data (existing method, for MCSA) of table 2
Fracture sliver 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 sum, presently preferred embodiments of the present invention is these are only, is not intended to limit protection scope of the present invention. All any modification, equivalent substitution and improvements within the spirit and principles in the present invention, made etc., should be included in the present invention's Within protection domain.

Claims (3)

1. a kind of three-phase cage asynchronous motor rotor broken splits the diagnostic method of sliver number, it is characterised in that specifically include as Lower step：
Step one, determine by sample frequency collection three-phase cage asynchronous motor threephase stator voltage transient signal and the three-phase of setting Electron current instantaneous signal；
Step 2, the virtual value for obtaining m phase stator voltages, m represents separate, m=1,2,3；
Step 3, obtain threephase stator voltage effective value arithmetic mean of instantaneous value and as stator voltage virtual value；
Step 4, the virtual value for obtaining m phase stator currents；
Step 5, obtain threephase stator current effective value arithmetic mean of instantaneous value and as stator current virtual value；
Step 6, according to threephase stator voltage transient signal and threephase stator electric current instantaneous signal, obtain instantaneous active power；
Step 7, the method by calculating mean value, obtain DC component mean (p) of instantaneous active power p, i.e. wattful power Rate；
Step 8, according to stator voltage virtual value and stator current virtual value, obtain apparent energy；Active power and apparent energy Ratio, as power factor
Step 9, to threephase stator in any one phase stator electric current instantaneous signal do FFT spectrum analysis, obtain stator current Middle f1The virtual value I of frequency componentS；Wherein, f1For frequency of supply；
Step 10, Hilbert transform is carried out to electric current instantaneous signal；
Step 11, acquisition Hilbert modulus Mod；
Step 12, the DC component for obtaining Hilbert modulus Mod, i.e. mean value A0；
Step 13, DC component A filtered in Hilbert modulus Mod0, obtain signal M to be analyzed；
Step 14, treat signal Analysis M and carry out FFT spectrum analysis, obtain its 2sf1Amplitude A of frequency component1；Wherein, s is Revolutional slip；
Step 15, basisObtain diagnosis rotor rupture sliver Number Nb, wherein, N is rotor bar sum.
2. a kind of three-phase cage asynchronous motor rotor broken as claimed in claim 1 splits the diagnostic method of sliver number, and it is special Levying is, the threephase stator voltage transient signal of the three-phase cage asynchronous motor is synchronous with threephase stator electric current instantaneous signal Collection.
3. a kind of three-phase cage asynchronous motor rotor broken as claimed in claim 1 splits the diagnostic method of sliver number, and it is special Levying is, the sample frequency of each signal is 1000Hz, and a length of 10s, i.e. sampling number are 10000 points during sampling.
CN201610843046.9A 2016-09-22 2016-09-22 A kind of three-phase cage asynchronous motor rotor broken splits the diagnostic method of conducting bar number Active CN106646222B (en)

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 true CN106646222A (en) 2017-05-10
CN106646222B 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)

Cited By (1)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
JAWAD FAIZ 等: "A New Criterion for Rotor Broken Bar Fault Diagnosis in Line-start and Inverter-fed Induction Motors using Hilbert-Huang Transform", 《2012 IEEE INTERNATIONAL CONFERENCE ON POWER ELECTRONICS, DRIVES AND ENERGY SYSTEMS》 *
许伯强 等: "基于复调制滤波与ESPRIT的异步电动机转子断条故障检测新方法", 《电测与仪表》 *
马宏忠 等: "基于Hilbert模量频谱分析的异步电机转子断条故障研究", 《电机与控制学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
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

Also Published As

Publication number Publication date
CN106646222B (en) 2019-06-21

Similar Documents

Publication Publication Date Title
CN101025430A (en) Cage type asynchronous motor rotor strip-broken failure detecting method
CN103454585B (en) A kind of permagnetic synchronous motor loss of excitation method for diagnosing faults based on dead electricity residual voltage
CN103018669B (en) Method for detecting failure of rotor broken bar of cage type asynchronous motor at high reliability
CN103969579B (en) Permanent-magnetic synchronous motor stator unbalanced fault diagnostic method and system
Wolbank et al. Monitoring of rotor-bar defects in inverter-fed induction machines at zero load and speed
CN101363901B (en) Method for detecting early failure of generator by enhancing transformations by electrical current characteristic
CN102135600A (en) Fault detection device and method of asynchronous motor
CN109991539A (en) Multiphase corner connection brushless exciter rotating diode one-phase open circuit detection method and system
CN109738780A (en) One tube open circuit detection method of multiphase corner connection brushless exciter rotating diode and system
CN101672894A (en) Motor rotor turn-to-turn short circuit diagnosis method based on high-order harmonic shaft voltage signals
CN106199424A (en) A kind of permagnetic synchronous motor interturn short-circuit failure diagnosing method
CN104111391B (en) Three-phase instantaneous power based motorized spindle fault monitoring and diagnosing method
Kim et al. High-frequency signal injection-based rotor bar fault detection of inverter-fed induction motors with closed rotor slots
CN106526479A (en) Grid-harmonic-characteristic-based fault diagnosis method for electric system of electric shovel
CN103744023A (en) Double-feed wind power generator stator winding asymmetric fault detection method
CN106646223A (en) Diagnosis method for number of three-phase cage asynchronous motor rotor broken guiding bars
Wang et al. Induction motor broken rotor bar fault diagnosis based on third-order energy operator demodulated current signal
CN107017808B (en) The continuous estimation method of synchronous motor rotor position based on pulsation exciting current response
CN106646222A (en) Method for diagnosing number of broken conducting bars of rotor of three-phase squirrel-cage asynchronous motor
Pietrzak et al. Stator phase current STFT analysis for the PMSM stator winding fault diagnosis
Mariun et al. An experimental study of induction motor current signature analysis techniques for incipient broken rotor bar detection
CN106646224A (en) Three-phase asynchronous motor load fluctuation detection method
CN107091988B (en) A kind of three-phase motor current signal instantaneous frequency and instantaneous envelope extracting method
CN105607000A (en) Method for measuring locked-rotor current of three-phase asynchronous motor under no-load condition
CN108445340A (en) The detection method of five-phase PMSM inverter open fault

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