[go: up one dir, main page]

CN101025430A - Cage type asynchronous motor rotor strip-broken failure detecting method - Google Patents

Cage type asynchronous motor rotor strip-broken failure detecting method Download PDF

Info

Publication number
CN101025430A
CN101025430A CN 200710061634 CN200710061634A CN101025430A CN 101025430 A CN101025430 A CN 101025430A CN 200710061634 CN200710061634 CN 200710061634 CN 200710061634 A CN200710061634 A CN 200710061634A CN 101025430 A CN101025430 A CN 101025430A
Authority
CN
China
Prior art keywords
fault
ratio
rotor
component
frequency
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.)
Pending
Application number
CN 200710061634
Other languages
Chinese (zh)
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.)
North China Electric Power University
Original Assignee
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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN 200710061634 priority Critical patent/CN101025430A/en
Publication of CN101025430A publication Critical patent/CN101025430A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Tests Of Circuit Breakers, Generators, And Electric Motors (AREA)

Abstract

一种笼型异步电动机转子断条故障检测方法及装置,属检测技术领域,用于解决在线检测转子断条初发故障的问题。其技术方案是:它通过对采集的定子电流瞬时信号is做连续细化傅里叶变换,得到其基波即参考信号uS,再根据参考信号uS及其频率f1对定子电流瞬时信号is做自适应滤波,然后对滤波输出信号eT做连续细化傅里叶变换,确定当前(1±2s)f1边频分量与基波分量幅值之比并把它作为故障特征,最后根据检测阈值确定故障指数,依故障指数判断是否存在转子断条故障。本发明可以高灵敏度、高可靠性地在线检测异步电动机初发性转子断条故障。

A cage-type asynchronous motor rotor broken bar fault detection method and device, which belong to the detection technology field, are used to solve the problem of online detection of rotor broken bar initial faults. Its technical solution is: it obtains its fundamental wave, that is, the reference signal u S by continuously refining the Fourier transform of the collected stator current instantaneous signal i s , and then calculates the stator current instantaneous signal u S and its frequency f 1 The signal i s is adaptively filtered, and then the filtered output signal e T is continuously refined by Fourier transform to determine the ratio of the current (1±2s) f 1 side frequency component to the amplitude of the fundamental component and take it as a fault feature , and finally determine the fault index according to the detection threshold, and judge whether there is a broken rotor bar fault according to the fault index. The invention can detect the incipient broken rotor bar fault of the asynchronous motor on-line with high sensitivity and high reliability.

Description

Cage type asynchronous motor rotor strip-broken failure detecting method
Technical field
The present invention relates to a kind of method and device that can online detection cage type asynchronous motor rotor strip-broken fault, belong to the detection technique field.
Background technology
Cage type asynchronous motor is in operational process, rotor bar is subjected to the effect of alterante stresses such as radial electromagnetic force, rotating electromagnetic power, centrifugal force, thermal flexure amount of deflection power, the rotor manufacturing defect all may cause broken bar fault in addition, and this kind fault probability of happening is about 15%.
Rotor broken bar is typical gradual fault, common 1,2 bar failure of initial stage, then development decline even shutdown so that motor is exerted oneself gradually.Therefore, must implement the online detection of rotor bar breaking fault, particularly the online detection of the property sent out rotor bar breaking fault just.
After the cage type asynchronous motor generation rotor bar breaking fault, (1 ± 2s) f will appear in its stator current 1(s is a revolutional slip to the extra current component of frequency, f 1Be line frequency), this current component is called the side frequency component, can be used as the rotor bar breaking fault feature.And stator current signal is easy to gather, and therefore the stator current signal frequency spectrum analysis method based on Fourier transform is widely used in the rotor bar breaking fault detection.
Initial rotor strip-broken failure detecting method is that the stable state stator current signal is directly carried out spectrum analysis, according to whether there being (1-2s) f in the spectrogram 1Frequency component judges that rotor has or not disconnected bar.When slightly breaking bar owing to rotor, (1-2s) f 1The amplitude of component is with respect to f 1Component is very little, and asynchronous motor when operation revolutional slip s is very little, (1-2s) f 1With f 1These two frequency numerical value are approaching, if directly do the Fourier spectrum analysis, and (1-2s) f then 1Component may be by f 1The leakage of component is flooded.This is the weak point of the method.
In order to remedy the deficiency of the method, development has formed adaptive filter method and starting current Time-Varied Spectrum Analysis method, and the core of adaptive filter method is: at first adopt adaptive filter method to offset stator current f 1Frequency component is carried out spectrum analysis afterwards again, and this can give prominence to rotor bar breaking fault characteristic component---(1-2s) f in spectrogram 1Frequency component, thus the rotor bar breaking fault detection sensitivity significantly improved.
Fig. 2 is the theory diagram of adaptive filter method.Among Fig. 2, i SRepresent actual stator current signal, it comprises signal S to be extracted TWith noise n T, and u SIt is reference signal.Here, S TBe (1-2s) f in the stator current 1Frequency component, n TBe the f in the stator current 1Frequency component, and e TThen represent i sMake auto adapted filtering and handle resulting signal afterwards.If the response of sef-adapting filter is y T, obviously, e T=i S-y TAccording to e TSize, adjust the parameter of wave filter, appropriate change y by adaptive algorithm T, can make y TUnder the meaning of least mean-square error, offset n T, and e TTo under the meaning of least mean-square error, approach signal S to be extracted T
When adopting adaptive filter method, noise u SAdopt test circuit shown in Figure 3 to obtain.Obviously, resistance R 1On voltage signal be i among Fig. 2 s, and resistance R 2On voltage signal when net capacity is enough big, only contain f 1Frequency component can be used as noise u SAmong Fig. 3, resistance R 1Effect be that secondary side current signal with current transformer CT is converted into a voltage signal that amplitude is suitable, resistance R 2Be connected to the secondary side of voltage transformer pt.
The weak point of adaptive filter method is following two aspects.At first, this method needs test circuit shown in Figure 3 to obtain noise signal u S, hardware circuit is slightly complicated.Secondly and since motor itself intrinsic asymmetric, air gap eccentric centre, rotor misalignment and other factor, even asynchronous motor is in normal operating condition, also may comprise (1 ± 2s) f in its stator current 1And other frequency component, and for different asynchronous motors, (Fig. 4, Fig. 5 are respectively Y100L-2 type, the Y100L1-4 type rotor normal motor stator a phase current auto adapted filtering frequency spectrum under full load conditions to the situation complexity; Fig. 5, Fig. 6 represent the stator a phase current auto adapted filtering frequency spectrum of Y100L-2 type motor under fully loaded and sliver of rotor and two bar failure situations respectively).This easily obscures mutually with the rotor bar breaking fault feature, causes erroneous judgement, influences the fault detect reliability.Adaptive filter method can not be taken into account above-mentioned factor, and it is not enough that reliability is still disliked.
Starting current Time-Varied Spectrum Analysis method is gathered the stator current signal of cage type asynchronous motor in start-up course, and it is done the segmentation spectrum analysis, obtains time varying spectrum, judges that in view of the above rotor has or not disconnected bar.This method possesses distinct advantages: at first, s is bigger in start-up course, thus in spectrogram (1-2s) f 1Frequency component can be away from f 1Frequency component, this can reduce the requirement to the spectrum analysis frequency resolution; Secondly, (1-2s) f in start-up course 1Frequency component and f 1The ratio of the amplitude of the ratio of frequency component amplitude during much larger than steady-state operation, thereby make (1-2s) f 1Frequency component is easy to detect.
The weak point of starting current Time-Varied Spectrum Analysis method is: must wait for that electric motor starting can use.Similar with adaptive filter method, this method can not take into account motor itself intrinsic asymmetric, air gap eccentric centre, rotor misalignment and other factor, reliability still remains to be improved.
Summary of the invention
The object of the present invention is to provide a kind of can high sensitivity, the asynchronous motor rotor strip-broken failure detecting method and the device of the online detection cage type asynchronous motor rotor strip-broken fault in high reliability ground.
The alleged problem of the present invention realizes with following technical proposals:
A kind of cage type asynchronous motor rotor strip-broken failure detecting method the steps include: that it passes through stator current momentary signal i to gathering sDo continuous refinement Fourier transform, obtaining its first-harmonic is reference signal u S, again according to reference signal u SAnd frequency f 1To stator current momentary signal i SDo auto adapted filtering, then to the filtering output signal e TDo continuous refinement Fourier transform, determine current (1 ± 2s) f 1Side frequency component and the ratio of fundametal compoment amplitude and it as fault signature, determine fault index according to detection threshold at last, judge whether to exist rotor bar breaking fault according to fault index.
Above-mentioned cage type asynchronous motor rotor strip-broken failure detecting method, fault index calculates as follows:
A. measure a phase stator current momentary signal i s:
For high-voltage motor, adopt 1 current clamp to measure a phase stator current momentary signal at current transformer CT secondary side; For low voltage motor, adopt 1 current clamp directly to measure a phase stator current momentary signal at the connecting terminal of motor place;
B. to stator current momentary signal i sDo continuous refinement Fourier transform, determine the frequency f of its fundametal compoment 1, amplitude I mWith initial phase angle φ, form reference signal u S:
For sample frequency is f s, sampling number is the time series i (t of N k),
u S(k)=I msin(2πf 1kT S+φ+π)
Wherein, k=0,1,2, Λ, N-1, I m, f 1, φ determines by continuous refinement Fourier transform.
C. according to the fundametal compoment frequency f 1, reference signal u STo stator current momentary signal i sDo auto adapted filtering, offset its fundametal compoment, obtain the filtering output signal e T
D. to the filtering output signal e TDo continuous refinement Fourier transform, inquiry (1 ± 2s) f in its continuous refinement spectrogram 1The side frequency component information is determined current (1 ± 2s) f 1Side frequency component and f 1The ratio ratio of component amplitude (1-2s) f1+ ratio (1+2s) f1,
Wherein, ratio (1-2s) f1Be (1-2s) f 1Side frequency component and f 1The ratio of component amplitude, ratio (1+2s) f1Be (1+2s) f 1Side frequency component and f 1The ratio of component amplitude;
E. determine fault index:
Under the situation of not setting up normal motor sample reference paper as yet, according to conventional experience detection threshold is set and (is approximately
Figure A20071006163400071
Z rBe the rotor slot number), ratio (1-2s) f1+ ratio (1+2s) f1Be fault index with the ratio of detection threshold;
F. whether exist according to the fault index failure judgement:
Fault index numerical value<1, the expression motor is in health status, and its numerical value is littler, and health status is clearer and more definite; Fault index numerical value>1, the expression motor is in malfunction, and its numerical value is bigger, and malfunction is more serious.
Above-mentioned cage type asynchronous motor rotor strip-broken failure detecting method is for can be in the filtering output signal e TContinuous refinement spectrogram in accurately the inquiry (1 ± 2s) f 1The side frequency component information, should calculate revolutional slip s earlier:
s = 1 - P Z r ( f rsh f 1 ± v ) , v = 1,3,5 , Λ
Wherein, f RshFor rotor tooth slot harmonic component frequency, P are the motor number of pole-pairs, Z rBe the rotor slot number.
Then according to the fundametal compoment frequency f 1, revolutional slip s, in the filtering output signal e TContinuous refinement spectrogram in the inquiry (1 ± 2s) f 1The side frequency component information is determined current (1 ± 2s) f 1Side frequency component and f 1The ratio ratio of component amplitude (1-2s) f1+ ratio (1+2s) f1
Above-mentioned cage type asynchronous motor rotor strip-broken failure detecting method, for eliminate real electrical machinery itself the influence of intrinsic asymmetric, air gap eccentric centre, rotor misalignment and other factor, should be under the rotor normal condition, according to revolutional slip s and fault signature ratio (1-2s) f1+ ratio (1+2s) f1Concrete numerical value set up sample database, and the detection threshold of adjusting in view of the above:
If the current numerical value of revolutional slip between sample data revolutional slip upper and lower limit, then adopts the linear interpolation mode that detection threshold is set; Otherwise, determine immediate with it sample data revolutional slip, as detection threshold, and make safety factor be not less than 1 the fault signature numerical value of correspondence.
The present invention is by gathering asynchronous motor stator winding current signal, and data acquisition card is sent to computing machine with this signal, by computing machine current signal is handled, and judges whether to exist rotor bar breaking fault, and operating process is simple and convenient.Utilize stator current (1 ± 2s) f 1Frequency component is as fault signature, continuous refinement Fourier transform, auto adapted filtering, the estimation of rotor tooth slot harmonic revolutional slip, detection threshold automatic-adjusting technique are organically combined, when improving sensitivity, eliminated motor itself intrinsic asymmetric, air gap eccentric centre, rotor misalignment and other factor to extracting the influence of fault signature, effectively prevented erroneous judgement.But high sensitivity, the online detection asynchronous motor in high reliability ground be the property sent out rotor bar breaking fault just.
Description of drawings
The invention will be further described below in conjunction with accompanying drawing.
Fig. 1 is an electric theory diagram of the present invention;
Fig. 2 is the theory diagram of adaptive filter method;
Fig. 3 is the schematic diagram of adaptive filter method signal acquisition circuit;
Fig. 4 is the stator a phase current auto adapted filtering frequency spectrum of Y100L-2 type rotor normal motor under full load conditions;
Fig. 5 is the stator a phase current auto adapted filtering frequency spectrum of Y100L1-4 type rotor normal motor under full load conditions;
The stator a phase current auto adapted filtering frequency spectrum that Fig. 6 is a Y100L-2 type motor under fully loaded and 1 bar failure situation;
The stator a phase current auto adapted filtering frequency spectrum that Fig. 7 is a Y100L-2 type motor under fully loaded and 2 bar failure situations;
Fig. 8 is the normal sample notebook data of experiment motor;
Stator a phase current when Fig. 9 is experiment motor semi-load and a bar failure;
Stator a phase current frequency spectrum when Figure 10 is experiment motor semi-load and a bar failure.
Each label is among the figure: PT, voltage transformer (VT), CT, current transformer, M, motor, R1, R2, resistance.
The meaning of used each symbol: s, revolutional slip in the literary composition; f 1, line frequency (fundamental frequency); i s, stator current signal; S T, auto adapted filtering signal to be extracted; n T, noise signal; u S, (auto adapted filtering) reference signal; e T, the filtering output signal; y T, filter response; I m, the reference signal amplitude; φ, reference signal initial phase angle; f s, sample frequency; N, sampling number; I (t k), the current sample instantaneous value; t k, sampling instant; Ratio (1-2s) f1, (1-2s) f 1Side frequency component and f 1The ratio of component amplitude; Ratio (1+2s) f1, (1+2s) f 1Side frequency component and f 1The ratio of component amplitude; f Rsh, rotor tooth slot harmonic component frequency, P, motor number of pole-pairs; Z r, the rotor slot number; T s, sampling time interval; A (n), b (n), a (0) represent fourier coefficient; Δ f, frequency discrimination unit.
Embodiment
The present invention adopts circuit shown in Figure 1 to detect, this circuit is made up of current transformer CT, data acquisition card and computing machine, described current transformer is connected on the phase line of asynchronous motor stator winding, its signal output part connects the simulating signal input channel 5 (input terminal 5 and 17) of data acquisition card, and the output port of described data acquisition card connects the USB mouth of computing machine.Data acquisition card adopts auspicious rich magnificent RBH8321 type data acquisition card, and the model of computing machine is DELL M1210, data acquisition card is integrated circuit such as low-pass filter, signals collecting maintenance, mould/number conversion.The stator current momentary signal is delivered to data acquisition card, and data acquisition card is connected in portable computer by USB interface.Portable computing machine control signal capture card is with appropriate frequency sampling stator current momentary signal, and is stored in hard disk, by computing machine current signal handled again, judges whether to exist rotor bar breaking fault.This matched with devices software is based on Windows XP operating system and adopt the establishment of Visual C++ application development platform.
Use continuous refinement Fourier transform and auto-adaptive filtering technique and can guarantee to extract in high sensitivity motor stator electric current (1 ± 2s) f 1Frequency component; Use rotor tooth slot harmonic revolutional slip estimation technique and can judge correctly whether this component is really caused by rotor bar breaking fault; Application then can suitably be provided with detection threshold based on the detection threshold of sample learning from the strategy of adjusting, and avoids fault omission and erroneous judgement.
Use continuous refinement fourier transform method, can obtain the accurate and analytical expression of a certain main frequency component in the signal to be analyzed, i.e. frequency, amplitude and initial phase angle.
For sample frequency is f s, sampling number is the time series i (t of N k), discrete Fourier progression is:
a ( n ) = 2 N Σ k = 0 N - 1 i ( t k ) cos ( 2 πkn / N ) b ( n ) = 2 N Σ k = 0 N - 1 i ( t k ) sin ( 2 πkn / N ) a ( 0 ) = 1 N Σ k = 0 N - 1 i ( t k ) n = 0,1,2 , Λ , N - 1 - - - ( 1 )
Wherein, t k=kT s, T s=1/f s, k=0,1,2, Λ, N-1, a (n), b (n), a (0) represent fourier coefficient.
Fast Fourier Transform (FFT) is the special circumstances of above-mentioned discrete transform, i.e. N=2 m(m is a positive integer), at this moment, Fourier transform can adopt the recursion fast algorithm.This conversion, frequency discrimination unit are Δ f=f s/ N, N is inversely proportional to sampling number.Obviously,, must increase sampling number exponentially if wish to improve the frequency discrimination ability, sampling number one timing, the frequency discrimination ability can't further improve.
Time series i (t k) comprise signal 0 to f sInformation in/2 these frequency domains if spectrum curve is regarded as continuously, thinks that promptly the n in the formula (1) is a continuous real number that belongs to interval [0, N/2], and formula (1) can be rewritten an accepted way of doing sth (2).At this moment, the frequency discrimination ability no longer is subjected to the restriction of sampling number, and the value of frequency f is continuous.
a ( f ) = 2 N &Sigma; k = 0 N - 1 i ( t k ) cos ( 2 &pi;kf / f s ) b ( f ) = 2 N &Sigma; k = 0 N - 1 i ( t k ) sin ( 2 &pi;kf / f s ) - - - ( 2 ) 0 < f &le; f s / 2
When using continuous refinement Fourier transform, refinement scope, refinement density can be carried out step by step, to improve computing velocity.
The sensitivity and the reliability of the fault detect of refinement Fourier pair raising cage type asynchronous motor rotor strip-broken are significant continuously.At first, can accurately determine motor stator current first harmonics frequency f by continuous refinement Fourier transform 1Secondly, can determine motor stator current first harmonics component accurate and analytical expression, i.e. frequency f 1, amplitude I mWith initial phase angle φ, form reference signal u in view of the above S, handle stator current signal is done auto adapted filtering.
Reference signal u SDetermine according to formula (3).
u S ( k ) = I m sin ( 2 &pi; f 1 k T S + &phi; + &pi; ) k = 0,1,2 , &Lambda; , N - 1 - - - ( 3 )
Fig. 2 is an adaptive filtration theory, and its basic ideas are: adopt adaptive filter method to offset motor stator electric current f 1Component, outstanding (1-2s) f in spectrogram 1Component---rotor bar breaking fault feature, thus the sensitivity that rotor bar breaking fault detects significantly improved.Reference signal u SUsing continuous refinement Fourier transform determines.
In motor operation course, because rotor mmf teeth groove harmonic wave and first-harmonic air-gap flux reciprocation will comprise rotor tooth slot harmonic component in the stator current.According to its frequency f Rsh, motor number of pole-pairs P and rotor slot count ZT and can determine motor slip ratio by formula (4).
s = 1 - P Z r ( f rsh f 1 &PlusMinus; v ) , v = 1,3,5 , &Lambda; - - - ( 4 )
In engineering reality, often select (1 ± 2s) f 1Frequency component and f 1The ratio ratio of component amplitude (1-2s) f1, ratio (1+2s) f1As fault signature, by judging whether its numerical value surpasses a certain threshold value and realize that rotor bar breaking fault detects.
For take into account real electrical machinery itself intrinsic asymmetric, air gap eccentric centre, rotor misalignment and other factor, should adopt detection threshold based on sample learning from the strategy of adjusting to improve sensitivity and reliability.
Suppose that initial rotor is normal, according to revolutional slip s and fault signature ratio (1-2s) f1+ ratio (1+2s) f1Concrete numerical value set up sample database, this is because fault signature ratio (1-2s) fi+ ratio (1+2s) f1Basically only depend on revolutional slip s.
Fig. 8 represents the normal sample notebook data of Y100L-2 type experiment motor.Set up sample database, should pay attention to its simple and direct property, and, contain the revolutional slip normal fluctuation as far as possible in conjunction with the motor actual operating.
As soon as the normal motor sample database is set up, can carry out from adjusting detection threshold according to the current numerical value of revolutional slip, specific as follows: if the current numerical value of revolutional slip between sample data revolutional slip upper and lower limit, then adopts the linear interpolation mode that detection threshold is set; Otherwise, determine immediate with it sample data revolutional slip, as detection threshold, and make safety factor be not less than 1 the fault signature numerical value of correspondence.
From Fig. 6, Fig. 7 as can be seen, rotor bar breaking fault after taking place in real electrical machinery, and its stator current frequency spectrum often comprises a plurality of spectrums peak, and (1 ± 2s) f 1Side frequency component spectrum peak-to-peak value may not be maximum.This be attributable to real electrical machinery itself intrinsic asymmetric, air gap eccentric centre, rotor misalignment and other factor.These the spectrum peaks very easily with rotor bar breaking fault feature---(1 ± 2s) f 1Frequency component is obscured mutually, causes erroneous judgement, influences the reliability that rotor bar breaking fault detects.Therefore, when carrying out the rotor bar breaking fault detection, should predict f 1With the concrete numerical value of s, thereby in the stator current frequency spectrum, on purpose inquire about (1 ± 2s) f 1The side frequency component, and it is done quantitative test, to guarantee the broken bar fault detecting reliability.
According to Fig. 4,5 as can be known, owing to technology, manufacturing and reason is installed, intrinsic asymmetric, air gap eccentric centre, rotor misalignment, the other factors in addition to a certain degree of any real electrical machinery certainty itself, even motor is in normal operating condition, also may comprise (1 ± 2s) f in its stator current 1And other frequency component.And for different asynchronous motors, the situation complexity.
How suitably this be provided with detection threshold to take into account this problem of sensitivity and reliability simultaneously with regard to having proposed.If it is too high that detection threshold is provided with, then be unfavorable for improving sensitivity.For example, rotor broken bar detection method in the past is mostly with (1-2s) f 1Component and f 1Likening to of component amplitude is fault signature, and detection threshold generally is set at 1~2%, and be obviously relatively conservative for Y100L-2 type experiment motor, will cause the fault omission.Because during bar failure of this rotor, its stator current (1-2s) f 1Component and f 1The ratio of component amplitude only is 0.67%.On the other hand, low if detection threshold was provided with, so that be lower than motor stator current (1-2s) f when normal operation 1Component and f 1The ratio of component amplitude must cause the fault erroneous judgement, and reliability is not just known where to begin yet.
The analysis showed that more than the rotor bar breaking fault that carries out high sensitivity, high reliability detects, must at first clear and definite normal motor (1 ± 2s) f 1Side frequency component and f 1The ratio of component amplitude is provided with suitable detection threshold in view of the above, to avoid fault omission and erroneous judgement.
The asynchronous motor just property sent out rotor strip-broken failure detecting method can be realized above-mentioned purpose by using based on the detection threshold of sample learning from the strategy of adjusting.
Fig. 6, Fig. 7 show, behind the cage type asynchronous motor generation rotor bar breaking fault, and stator current (1-2s) f 1, (1+2s) f 1Frequency component is that coupling occurs, and for disconnected bar detection sensitivity and reliability consideration, should inquire about (1 ± 2s) f simultaneously 1Frequency component information.
Stator a phase current frequency spectrum when the stator a phase current when Fig. 9 represents Y100L-2 type experiment motor semi-load and a bar failure, Figure 10 are represented Y100L-2 type experiment motor semi-load and a bar failure, concrete data are shown in table 2.
Experimental result under table 2 motor semi-load and the bar failure situation
f 1(Hz) 50.01
s(%) 1.78
(1-2s)f 1(Hz) 48.27
(1+2s)f 1(Hz) 51.83
f 1Component amplitude (A) 3.7113
(1-2s)f 1Component amplitude (A) 0.0227
(1+2s)f 1Component amplitude (A) 0.0235
Fault signature amount ratio (1-2s)f1+ratio (1+2s)f1(%) 1.24
Do not set up the detection threshold (%) under the normal motor sample reference paper situation as yet 1.0
Set up the detection threshold (%) under the normal motor sample reference paper situation 0.40

Claims (4)

1、一种笼型异步电动机转子断条故障检测方法,其特征是,它通过对采集的定子电流瞬时信号is做连续细化傅里叶变换,得到其基波参考信号uS,再根据参考信号uS及其频率f1对定子电流瞬时信号is做自适应滤波,然后对滤波输出信号eT做连续细化傅里叶变换,确定当前(1±2s)f1边频分量与基波分量幅值之比并把它作为故障特征,最后根据检测阈值确定故障指数,依故障指数判断转子断条故障。1. A method for detecting broken rotor bars of cage-type asynchronous motors, which is characterized in that it obtains its fundamental wave reference signal u S by performing continuous refinement Fourier transform on the collected stator current instantaneous signal i s , and then according to The reference signal u S and its frequency f 1 do adaptive filtering on the stator current instantaneous signal i s , and then perform continuous refinement Fourier transform on the filtered output signal e T to determine the current (1±2s) f 1 side frequency component and The ratio of the fundamental wave component amplitude is used as the fault feature, and finally the fault index is determined according to the detection threshold, and the rotor broken bar fault is judged according to the fault index. 2、根据权利要求1所述笼型异步电动机转子断条故障检测方法,其特征是,故障指数按如下步骤计算:2. According to claim 1, the method for detecting broken rotor bars of cage-type asynchronous motors is characterized in that the fault index is calculated according to the following steps: a.测取一相定子电流瞬时信号isa. Measure the instantaneous signal i s of one-phase stator current: b.对定子电流瞬时信号is做连续细化傅里叶变换,确定其基波分量的频率f1、幅值Im和初相角φ,形成参考信号uSb. Continuously refine the Fourier transform of the instantaneous signal i s of the stator current to determine the frequency f 1 , amplitude I m and initial phase angle φ of the fundamental component to form a reference signal u S : 对于采样频率为fs,采样点数为N的时间序列i(tk),For the time series i(t k ) whose sampling frequency is f s and the number of sampling points is N, uS(K)=Im sin(2πf1kTS+φ+π)u S (K)=I m sin(2πf 1 kT S +φ+π) 其中,k=0,1,2,Λ,N-1,Im、f1、φ通过连续细化傅里叶变换确定。Wherein, k=0, 1, 2, Λ, N-1, Im , f 1 , φ are determined by continuous refinement Fourier transform. c.根据基波分量频率f1、参考信号uS对定子电流瞬时信号is做自适应滤波,抵消其基波分量,得到滤波输出信号eTc. According to the frequency f 1 of the fundamental wave component and the reference signal u S , perform adaptive filtering on the instantaneous signal is of the stator current to offset its fundamental wave component, and obtain the filtered output signal e T ; d.对滤波输出信号eT做连续细化傅里叶变换,在其连续细化频谱图中查询(1±2s)f1边频分量信息,确定当前(1±2s)f1边频分量与f1分量幅值之比ratio(1-2s)f1+ratio(1+2s)f1d. Perform continuous refinement Fourier transform on the filtered output signal e T , query (1±2s) f 1 side frequency component information in its continuous refinement spectrogram, and determine the current (1±2s) f 1 side frequency component Ratio (1-2s)f1 +ratio (1+2s)f1 with the amplitude of the f 1 component, 其中,ratio(1-2s)f1是(1-2s)f1边频分量与f1分量幅值之比,ratio(1+2s)f1是(1+2s)f1边频分量与f1分量幅值之比;Among them, ratio (1-2s) f1 is the ratio of (1-2s) f 1 side frequency component to f 1 component amplitude, ratio (1+2s) f1 is (1+2s) f 1 side frequency component and f 1 ratio of component magnitudes; e.确定故障指数:e. Determine the failure index: 在尚未建立正常电机样本参考文件的情况下,根据常规经验设置检测阈值,近似为
Figure A2007100616340002C1
其中,Zr为转子槽数,ratio(1-2s)f1+ratio(1+2s)f1与检测阈值的比值即为故障指数;
In the case that no normal motor sample reference files have been established, the detection threshold is set according to conventional experience, which is approximately
Figure A2007100616340002C1
Among them, Zr is the number of rotor slots, and the ratio of ratio (1-2s)f1 +ratio (1+2s)f1 to the detection threshold is the failure index;
f.根据故障指数判断故障:f. Judging the fault according to the fault index: 若故障指数数值<1,表示电机处于健康状态,且其数值愈小,健康状态愈明确;故障指数数值>1,表示电机处于故障状态,且其数值愈大,故障状态愈严重。If the fault index value<1, it means that the motor is in a healthy state, and the smaller the value, the clearer the health state; if the fault index value>1, it means that the motor is in a fault state, and the larger the value, the more serious the fault state.
3、根据权利要求1或2所述笼型异步电动机转子断条故障检测方法,其特征是,确定转差率s:3. According to claim 1 or 2, the method for detecting broken rotor bars of cage-type asynchronous motors is characterized by determining the slip s: sthe s == 11 -- PP ZZ rr (( ff rshrsh ff 11 &PlusMinus;&PlusMinus; vv )) ,, vv == 1,3,51,3,5 ,, &Lambda;&Lambda; 其中,frsh为转子齿槽谐波分量频率、P为电机极对数,Zr为转子槽数,Among them, f rsh is the harmonic frequency of the rotor cogging, P is the number of pole pairs of the motor, Z r is the number of rotor slots, 然后根据基波分量频率f1、转差率s,在滤波输出信号eT的连续细化频谱图中查询(1±2s)f1边频分量信息,确定当前(1±2s)f1边频分量与f1分量幅值之比ratio(1-2s)f1+ratio(1+2s)f1Then, according to the fundamental component frequency f 1 and the slip rate s, query the (1±2s) f 1 edge frequency component information in the continuous refinement spectrogram of the filtered output signal e T to determine the current (1±2s) f 1 edge Ratio (1-2s)f1 +ratio (1+2s)f1 of frequency component and f 1 component amplitude. 4、根据权利要求3所述笼型异步电动机转子断条故障检测方法,其特征是,在电机转子正常情况下,根据转差率s与故障特征ratio(1-2s)f1+ratio(1+2s)f1的具体数值建立样本数据库,并据此整定检测阈值:4. According to claim 3, the method for detecting the broken bars of the cage-type asynchronous motor rotor is characterized in that, under the normal condition of the motor rotor, according to the slip s and the fault characteristic ratio (1-2s)f1 +ratio (1+ 2s) Establish a sample database based on the specific value of f1 , and set the detection threshold accordingly: 若转差率当前数值介于样本数据转差率上、下限之间,则采用线性内插方式设置检测阈值;否则,确定与之最接近的样本数据转差率,将对应的故障特征数值作为检测阈值,并使可靠系数不小于1。If the current value of the slip is between the upper and lower limits of the sample data slip, use linear interpolation to set the detection threshold; otherwise, determine the closest sample data slip, and use the corresponding fault characteristic value as Check the threshold and make the reliability factor not less than 1.
CN 200710061634 2007-03-28 2007-03-28 Cage type asynchronous motor rotor strip-broken failure detecting method Pending CN101025430A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 200710061634 CN101025430A (en) 2007-03-28 2007-03-28 Cage type asynchronous motor rotor strip-broken failure detecting method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 200710061634 CN101025430A (en) 2007-03-28 2007-03-28 Cage type asynchronous motor rotor strip-broken failure detecting method

Publications (1)

Publication Number Publication Date
CN101025430A true CN101025430A (en) 2007-08-29

Family

ID=38743910

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 200710061634 Pending CN101025430A (en) 2007-03-28 2007-03-28 Cage type asynchronous motor rotor strip-broken failure detecting method

Country Status (1)

Country Link
CN (1) CN101025430A (en)

Cited By (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102135600A (en) * 2011-01-28 2011-07-27 中国人民解放军海军潜艇学院 Fault detection device and method of asynchronous motor
CN102279343A (en) * 2011-07-23 2011-12-14 华北电力大学(保定) Online detection method of broken rotor bar fault of cage-type asynchronous motor based on ESPRIT-SAA
CN102279341A (en) * 2011-07-23 2011-12-14 华北电力大学(保定) Cage asynchronous motor rotor broken-bar fault detection method based on electronic stability program rotation invariant technology (ESPRIT) and pattern search algorithm (PSA)
CN102279364A (en) * 2011-07-23 2011-12-14 华北电力大学(保定) Method for detecting broken bar fault of cage type asynchronous motor rotor based on multiple signal classification (MUSIC) and pattern search algorithm (PSA)
CN102279342A (en) * 2011-07-23 2011-12-14 华北电力大学(保定) Cage asynchronous motor rotor broken-bar fault detection method based on multi-signal classification technology (MUSIC) and simulation annealing algorithm (SAA)
CN102707232A (en) * 2012-06-01 2012-10-03 深圳市海亿达能源科技股份有限公司 Online monitoring device of state of motor device and monitoring method thereof
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
CN103033749A (en) * 2012-12-31 2013-04-10 南京工程学院 Digital detection system for detecting rotor broken bar failure of asynchronous motor and detection method of digital detection system
CN103185862A (en) * 2011-12-31 2013-07-03 上海汽车集团股份有限公司 Method and device for diagnosis of direct current motor breakdown in automobile starter
CN103261904A (en) * 2010-12-21 2013-08-21 西门子公司 Monitoring and fault diagnosis of an electric machine
CN103328995A (en) * 2011-02-03 2013-09-25 罗伯特·博世有限公司 Method and device for detecting a malfunction of an electric machine
CN103827683A (en) * 2011-09-30 2014-05-28 Abb技术有限公司 A method of determining stationary signals for the diagnostics of an electromechanical system
CN104765003A (en) * 2015-02-11 2015-07-08 湘潭大学 Asynchronous motor rotor bar break fault diagnosis method for engineering machine internet of things
CN105862335A (en) * 2016-04-06 2016-08-17 珠海格力电器股份有限公司 Method and device for determining eccentric load of washing machine
CN106680716A (en) * 2017-01-20 2017-05-17 浙江大学 Fault diagnosis method for the permanent magnet motor bearing based on the sensor-less
CN106842023A (en) * 2017-01-22 2017-06-13 浙江大学 The method for diagnosing faults of electric rotating machine
CN106872893A (en) * 2017-01-24 2017-06-20 杨明扬 A kind of motor rotor broken-bar fault processing method and system
CN107091985A (en) * 2017-04-21 2017-08-25 武汉科技大学 Asynchronous Motor Rotor-Bar Fault diagnostic system and method
CN107544025A (en) * 2017-08-30 2018-01-05 马鞍山马钢华阳设备诊断工程有限公司 A kind of Asynchronous Motor Rotor-Bar Fault determination methods of composite electrical signal and vibration signal
CN107589373A (en) * 2017-08-30 2018-01-16 马鞍山马钢华阳设备诊断工程有限公司 A kind of Asynchronous Motor Rotor-Bar Fault determination methods
CN107748885A (en) * 2017-06-30 2018-03-02 南京理工大学 A kind of method that ambiguous characters are identified
CN109031118A (en) * 2018-06-29 2018-12-18 珠海格力电器股份有限公司 Motor fault diagnosis method and device and motor
CN110531261A (en) * 2019-08-27 2019-12-03 国网江苏省电力有限公司检修分公司 A kind of non-intrusion type initial failure of asynchronous electromotor detection method and system
CN111398813A (en) * 2020-03-30 2020-07-10 西安交通大学 A method for diagnosing broken rotor bars of squirrel-cage three-phase asynchronous motors
CN111458640A (en) * 2020-04-22 2020-07-28 西安交通大学 Three-phase current-based alternating current asynchronous motor rotor broken bar fault diagnosis method
CN111537881A (en) * 2020-05-26 2020-08-14 华润智慧能源有限公司 Fault diagnosis method, device and equipment for asynchronous motor and readable storage medium
CN112254910A (en) * 2020-10-10 2021-01-22 华中科技大学 A method and system for fault detection of motor rotor broken bars based on fusion correlation spectrum
CN113703371A (en) * 2021-09-01 2021-11-26 山西华控伟业科技有限公司 Device and method for detecting equipment fault
CN114441169A (en) * 2021-12-02 2022-05-06 北京科锐特科技有限公司 Method for detecting outer raceway fault of rolling bearing of three-phase asynchronous motor based on multiple modulation components of rotation frequency in stator current signal
CN115586345A (en) * 2022-12-09 2023-01-10 深圳市双合电气股份有限公司 Motor self-adaptive rotating speed estimation method and system based on stator current
RU2791428C1 (en) * 2022-10-17 2023-03-07 федеральное государственное бюджетное образовательное учреждение высшего образования "Ивановский государственный энергетический университет имени В.И. Ленина" (ИГЭУ) Method for detecting broken rods in a short-closed rotor winding of asynchronous electric motor
CN115955161A (en) * 2023-03-15 2023-04-11 清华大学无锡应用技术研究院 Adaptive asynchronous induction motor slip estimation method, device, equipment and medium
CN118501695A (en) * 2024-07-19 2024-08-16 国网山东省电力公司莱芜供电公司 Air gap eccentricity fault detection method for variable frequency induction motor based on frequency adaptive FRFT filtering

Cited By (50)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103261904A (en) * 2010-12-21 2013-08-21 西门子公司 Monitoring and fault diagnosis of an electric machine
CN103261904B (en) * 2010-12-21 2016-05-04 西门子公司 The monitoring of electrically powered machine and fault diagnosis
CN102135600A (en) * 2011-01-28 2011-07-27 中国人民解放军海军潜艇学院 Fault detection device and method of asynchronous motor
CN103328995A (en) * 2011-02-03 2013-09-25 罗伯特·博世有限公司 Method and device for detecting a malfunction of an electric machine
CN103328995B (en) * 2011-02-03 2015-09-16 罗伯特·博世有限公司 For identifying the method and apparatus of the fault characteristic of motor
US9213065B2 (en) 2011-02-03 2015-12-15 Robert Bosch Gmbh Method and device for detecting a malfunction of an electric machine through the modulation of variables influencing the torque of the machine
CN102279343A (en) * 2011-07-23 2011-12-14 华北电力大学(保定) Online detection method of broken rotor bar fault of cage-type asynchronous motor based on ESPRIT-SAA
CN102279341A (en) * 2011-07-23 2011-12-14 华北电力大学(保定) Cage asynchronous motor rotor broken-bar fault detection method based on electronic stability program rotation invariant technology (ESPRIT) and pattern search algorithm (PSA)
CN102279364A (en) * 2011-07-23 2011-12-14 华北电力大学(保定) Method for detecting broken bar fault of cage type asynchronous motor rotor based on multiple signal classification (MUSIC) and pattern search algorithm (PSA)
CN102279342A (en) * 2011-07-23 2011-12-14 华北电力大学(保定) Cage asynchronous motor rotor broken-bar fault detection method based on multi-signal classification technology (MUSIC) and simulation annealing algorithm (SAA)
CN102279343B (en) * 2011-07-23 2013-08-07 华北电力大学(保定) Online detection method of broken rotor bar fault of cage-type asynchronous motor based on ESPRIT-SAA
CN103827683A (en) * 2011-09-30 2014-05-28 Abb技术有限公司 A method of determining stationary signals for the diagnostics of an electromechanical system
CN103827683B (en) * 2011-09-30 2017-09-12 Abb瑞士股份有限公司 It is determined that the method for the stationary signal diagnosed for Mechatronic Systems
CN103185862B (en) * 2011-12-31 2016-02-17 上海汽车集团股份有限公司 The diagnostic method of the direct current generator fault in automobile starter and device
CN103185862A (en) * 2011-12-31 2013-07-03 上海汽车集团股份有限公司 Method and device for diagnosis of direct current motor breakdown in automobile starter
CN102707232A (en) * 2012-06-01 2012-10-03 深圳市海亿达能源科技股份有限公司 Online monitoring device of state of motor device and monitoring method thereof
CN102707232B (en) * 2012-06-01 2015-10-07 深圳市海亿达能源科技股份有限公司 Motor apparatus state on_line monitoring device and monitoring method thereof
CN103018669A (en) * 2012-11-30 2013-04-03 华北电力大学(保定) Method for detecting failure of rotor broken bar of cage type asynchronous motor at high reliability
CN102944842B (en) * 2012-11-30 2015-11-18 华北电力大学(保定) A kind of cage type asynchronous motor rotor strip-broken failure detecting method
CN103018669B (en) * 2012-11-30 2015-03-25 华北电力大学(保定) Method for detecting failure of rotor broken bar of cage type asynchronous motor at high reliability
CN102944842A (en) * 2012-11-30 2013-02-27 华北电力大学(保定) Detecting method for rotor broken bar fault of cage-type asynchronous motor
CN103033749A (en) * 2012-12-31 2013-04-10 南京工程学院 Digital detection system for detecting rotor broken bar failure of asynchronous motor and detection method of digital detection system
CN104765003A (en) * 2015-02-11 2015-07-08 湘潭大学 Asynchronous motor rotor bar break fault diagnosis method for engineering machine internet of things
CN105862335A (en) * 2016-04-06 2016-08-17 珠海格力电器股份有限公司 Method and device for determining eccentric load of washing machine
CN105862335B (en) * 2016-04-06 2018-06-01 珠海格力电器股份有限公司 Method and device for determining eccentric load of washing machine
CN106680716A (en) * 2017-01-20 2017-05-17 浙江大学 Fault diagnosis method for the permanent magnet motor bearing based on the sensor-less
CN106680716B (en) * 2017-01-20 2019-01-08 浙江大学 A method of the magneto bearing failure diagnosis based on position-sensor-free
CN106842023A (en) * 2017-01-22 2017-06-13 浙江大学 The method for diagnosing faults of electric rotating machine
CN106842023B (en) * 2017-01-22 2019-05-21 浙江大学 The method for diagnosing faults of rotating electric machine
CN106872893A (en) * 2017-01-24 2017-06-20 杨明扬 A kind of motor rotor broken-bar fault processing method and system
CN106872893B (en) * 2017-01-24 2019-03-22 杨明扬 A kind of motor rotor broken-bar fault processing method and system
CN107091985A (en) * 2017-04-21 2017-08-25 武汉科技大学 Asynchronous Motor Rotor-Bar Fault diagnostic system and method
CN107748885B (en) * 2017-06-30 2021-02-12 南京理工大学 Method for recognizing fuzzy character
CN107748885A (en) * 2017-06-30 2018-03-02 南京理工大学 A kind of method that ambiguous characters are identified
CN107589373A (en) * 2017-08-30 2018-01-16 马鞍山马钢华阳设备诊断工程有限公司 A kind of Asynchronous Motor Rotor-Bar Fault determination methods
CN107544025A (en) * 2017-08-30 2018-01-05 马鞍山马钢华阳设备诊断工程有限公司 A kind of Asynchronous Motor Rotor-Bar Fault determination methods of composite electrical signal and vibration signal
CN109031118A (en) * 2018-06-29 2018-12-18 珠海格力电器股份有限公司 Motor fault diagnosis method and device and motor
CN110531261A (en) * 2019-08-27 2019-12-03 国网江苏省电力有限公司检修分公司 A kind of non-intrusion type initial failure of asynchronous electromotor detection method and system
CN111398813A (en) * 2020-03-30 2020-07-10 西安交通大学 A method for diagnosing broken rotor bars of squirrel-cage three-phase asynchronous motors
CN111458640A (en) * 2020-04-22 2020-07-28 西安交通大学 Three-phase current-based alternating current asynchronous motor rotor broken bar fault diagnosis method
CN111537881A (en) * 2020-05-26 2020-08-14 华润智慧能源有限公司 Fault diagnosis method, device and equipment for asynchronous motor and readable storage medium
CN112254910A (en) * 2020-10-10 2021-01-22 华中科技大学 A method and system for fault detection of motor rotor broken bars based on fusion correlation spectrum
CN112254910B (en) * 2020-10-10 2022-05-17 华中科技大学 A method and system for fault detection of motor rotor broken bars based on fusion correlation spectrum
CN113703371A (en) * 2021-09-01 2021-11-26 山西华控伟业科技有限公司 Device and method for detecting equipment fault
CN114441169A (en) * 2021-12-02 2022-05-06 北京科锐特科技有限公司 Method for detecting outer raceway fault of rolling bearing of three-phase asynchronous motor based on multiple modulation components of rotation frequency in stator current signal
CN114441169B (en) * 2021-12-02 2024-02-09 北京科锐特科技有限公司 Method for detecting faults of outer raceway of rolling bearing of three-phase asynchronous motor
RU2791428C1 (en) * 2022-10-17 2023-03-07 федеральное государственное бюджетное образовательное учреждение высшего образования "Ивановский государственный энергетический университет имени В.И. Ленина" (ИГЭУ) Method for detecting broken rods in a short-closed rotor winding of asynchronous electric motor
CN115586345A (en) * 2022-12-09 2023-01-10 深圳市双合电气股份有限公司 Motor self-adaptive rotating speed estimation method and system based on stator current
CN115955161A (en) * 2023-03-15 2023-04-11 清华大学无锡应用技术研究院 Adaptive asynchronous induction motor slip estimation method, device, equipment and medium
CN118501695A (en) * 2024-07-19 2024-08-16 国网山东省电力公司莱芜供电公司 Air gap eccentricity fault detection method for variable frequency induction motor based on frequency adaptive FRFT filtering

Similar Documents

Publication Publication Date Title
CN101025430A (en) Cage type asynchronous motor rotor strip-broken failure detecting method
CN102707232B (en) Motor apparatus state on_line monitoring device and monitoring method thereof
Zhang et al. A novel detection method of motor broken rotor bars based on wavelet ridge
CN100543440C (en) Failure testing method of asynchronous motor bearing
Filippetti et al. Condition monitoring and diagnosis of rotor faults in induction machines: State of art and future perspectives
CN103926533B (en) Permagnetic synchronous motor loss of excitation on-line fault diagnosis method and system
CN100561246C (en) Method for Automatically Identifying Turn-to-Turn Short-Circuit Faults of Generator Rotor Windings under Load Conditions
US11722079B2 (en) Detection system and method for rotor dynamic turn-to-turn short circuit fault of synchronous generator
CN102135600A (en) Fault detection device and method of asynchronous motor
CN103344368B (en) Based on the squirrel cage asynchronous motor efficiency on-line monitoring method can surveying electric parameters
CN101101319A (en) Generator rotor interturn short-circuit state detection device
CN103235260A (en) Submersible motor rotor broken bar fault recognition method based on HHT (Hilbert-Huang transform)
CN101025434A (en) Asynchronous motor stator winding inter-turn short circuit failure on-line detecting method and device
CN101363901B (en) Method for Early Fault Detection of Motor Using Current Signature Enhancement Transformation
CN102495368A (en) Non-invasive online detection method and system for rotor broken-bar fault of three-phase cage type asynchronous motor
CN104111391B (en) Three-phase instantaneous power based motorized spindle fault monitoring and diagnosing method
CN112254910B (en) A method and system for fault detection of motor rotor broken bars based on fusion correlation spectrum
CN106597277B (en) Frequency control Rotor Bars Breaking Fault in Squirrel-cage Asynchronous Motors dynamic diagnosis method and device
CN106772028B (en) Short Circuit Between Generator Rotor Windings diagnostic device and method based on excitation step method
CN111537881A (en) Fault diagnosis method, device and equipment for asynchronous motor and readable storage medium
Zhang et al. Diagnosis of mechanical unbalance fault in permanent magnet synchronous machine drives
CN106646223B (en) A kind of three-phase cage asynchronous motor rotor broken splits the diagnostic method of conducting bar number
Pilloni et al. Fault detection in induction motors
CN106842023A (en) The method for diagnosing faults of electric rotating machine
CN201177650Y (en) Cage induction motor with detection coil and detection device for broken rotor bar

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Open date: 20070829