CN104215905A - Motor fault diagnosis method based on Mahalanobis-Taguchi system and Box-Cox transformation - Google Patents
Motor fault diagnosis method based on Mahalanobis-Taguchi system and Box-Cox transformation Download PDFInfo
- Publication number
- CN104215905A CN104215905A CN201410452287.1A CN201410452287A CN104215905A CN 104215905 A CN104215905 A CN 104215905A CN 201410452287 A CN201410452287 A CN 201410452287A CN 104215905 A CN104215905 A CN 104215905A
- Authority
- CN
- China
- Prior art keywords
- motor
- lambda
- box
- characteristic parameter
- mahalanobis distance
- 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
Links
- 238000003745 diagnosis Methods 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000009466 transformation Effects 0.000 title abstract 5
- 238000013461 design Methods 0.000 claims abstract description 6
- 238000006243 chemical reaction Methods 0.000 claims description 18
- 230000003862 health status Effects 0.000 claims description 14
- 238000012360 testing method Methods 0.000 claims description 10
- 238000007476 Maximum Likelihood Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 230000001131 transforming effect Effects 0.000 abstract 1
- 238000012549 training Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 2
- 241001391944 Commicarpus scandens Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000013501 data transformation Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Tests Of Circuit Breakers, Generators, And Electric Motors (AREA)
Abstract
The invention discloses a motor fault diagnosis method based on a Mahalanobis-Taguchi system and Box-Cox transformation. The motor fault diagnosis method includes the steps of firstly, acquiring vibration signals when a motor is in normal operation and in fault; secondly, computing time-domain and time-frequency-domain characteristic parameters of the vibration signals; thirdly, applying a robust design concept in the Mahalanobis-Taguchi system to the characteristic parameters and selecting useful characteristic parameters for computing Mahalanobis distances when the motor is operated in different health states; fourthly, transforming the Mahalanobis distances obtained by computation in the step three to Gaussian distribution data by means of Box-Cox transformation, and determining a Mahalanobis distance range when the motor is operated in different health states by the aid of the Gaussian distribution property and inverse Box-Cox transformation; fifthly, subjecting the motor in an unknown operation state to signal acquisition, useful characteristic parameter computation and Mahalanobis distance computation, comparing the Mahalanobis distance range determined in the step four and judging the health condition of the tested motor, so that motor fault diagnosis is achieved. The motor fault diagnosis method based on the Mahalanobis-Taguchi system and Box-Cox transformation is capable of achieving diagnosis for electrical faults of the motor effectively.
Description
Technical field
The present invention relates to a kind of Method of Motor Fault Diagnosis.
Background technology
Along with the development of AC electric power systems, alternating current generator is used in commercial production widely.Because it uses frequent, load torque fluctuation, the reason such as overheated of transshipping, motor is easy to break down.Electrical fault not only can affect the production of industry, but also may cause the security incidents such as casualties.Therefore, for guaranteeing motor security of operation, reliable, the enforcement of Method of Motor Fault Diagnosis then seems particularly important.The fault type of motor can be divided into mechanical fault and electric fault two aspect.The diagnosis of motor electric fault under unbalanced input voltages that what the present invention was specifically related to is.Existing diagnostic method reasoning from logic is complicated, assesses the cost higher, is more difficultly actually applied.
Summary of the invention
In order to overcome the more difficult deficiency being applicable to motor electric fault diagnosis of existing method for diagnosing faults, the invention provides a kind of Method of Motor Fault Diagnosis converted based on horse field system and Box-Cox effectively realizing motor electric fault diagnosis.
The technical solution adopted for the present invention to solve the technical problems is:
Based on the Method of Motor Fault Diagnosis that horse field system and Box-Cox convert, it is characterized in that: said method comprising the steps of:
Step one, gather the vibration signal of motor when normally working and break down;
Step 2, the vibration signal obtained step one carry out the calculating of characteristic parameter, specifically: calculating vibration signal being carried out to time domain charactreristic parameter, simultaneously to vibration signal by wavelet package transforms, carry out the calculating of characteristic parameter on time-frequency domain;
Step 3, to the characteristic parameter calculated in step 2, characteristic parameter when normally working with motor is benchmark, uses the based Robust Design theory in the system of horse field to select useful characteristic parameter, and recalculates mahalanobis distance based on these characteristic parameters;
Step 4, to the non-negative that step 3 obtains, the mahalanobis distance of non-gaussian distribution, use Box-Cox conversion, convert the data of Gaussian distribution to, utilize the character of Gaussian distribution and inverse Box-Cox conversion, determine the mahalanobis distance scope of machine operation when different health status;
The one group of data x=[x be made up of mahalanobis distance
1, x
2..., x
n], N is the number of mahalanobis distance sample, i-th data x in x
ithe data y obtained after Box-Cox conversion
i(λ) calculate by following formula:
Wherein convert parameter lambda to be estimated by the maximum likelihood function of following formula:
Wherein, y=[y
1, y
2..., y
n],
before and after Box-Cox converts, the distribution situation of data can be undertaken judging or verifying by histogram or normal probability plot;
Then use the character of Gaussian distribution, obtain the upper and lower limit of the data of the rear Gaussian distributed of Box-Cox conversion, and determine the mahalanobis distance scope of machine operation under different health status by the inverse Box-Cox conversion of following formula:
Step 5, to the motor of unknown duty through vibration signals collecting, useful feature parameter calculates and mahalanobis distance calculates, the mahalanobis distance scope determined in contrast step 4, judges the health status of testing of electric motors, thus realizes the fault diagnosis of motor.
Further, in described step 3, use the based Robust Design theory in the system of horse field to select useful characteristic parameter, comprise the following steps:
Step 3.1, to set after a stack features standard parameter characteristic of correspondence vector as z
i, its dimension is n, the mahalanobis distance MD that these group data are corresponding
ifor
Wherein, C
-1for the inverse matrix of correlation matrix,
for character vector z
itransposition;
Step 3.2, use Taguchi's method are optimized characteristic parameter and choose;
Quantity n according to characteristic parameter in proper vector selects two-level orthogonal array L
p(t
c), wherein p is test number (TN), and t=2 is number of levels, and definition level '+' is " using this characteristic parameter ", and level '-' is " not using this characteristic parameter ", and c is columns, represents and allows to arrange maximum characteristic parameter numbers; According to 13 characteristic parameter selection orthogonal arrage L that step 2 builds
16(2
15); To a jth characteristic parameter, if
with
be respectively the signal to noise ratio (S/N ratio) average of this characteristic parameter under level '+' and level '-';
represent and adopt this characteristic parameter on average to detect effect to fault,
represent and do not adopt this characteristic parameter on average to detect effect to fault, therefore select satisfied
useful feature parameter.
In described step 2, vibration signal is carried out to the calculating of time domain charactreristic parameter, include valid value, peak value, kurtosis, measure of skewness, crest factor, nargin coefficient, impulse ratio, form factor and root mean square; Use wavelet package transforms that vibration signal is decomposed 2 layers simultaneously, calculate 4 characteristic parameters about energy in wavelet packet coefficient.These 13 parameters use alphabetical A respectively
1, A
2..., A
13represent.
Technical conceive of the present invention is: the diagnosis of motor electric fault under unbalanced input voltages that what the present invention was specifically related to is.Mahalanobis distance has been used to the health status of characterization device, but this parameter has the characteristic of non-negative and non-gaussian distribution, determines the mahalanobis distance scope of machine operation under different health status, usually more difficult.Converted by Box-Cox, and use the character of Gaussian distribution, in the training stage, construct the mahalanobis distance scope of machine operation under different health status; At test phase, calculate the mahalanobis distance of unknown health status testing of electric motors, by the contrast of mahalanobis distance scope built with the training stage, judge the health status of testing of electric motors, realize the fault diagnosis of motor.
Beneficial effect of the present invention is mainly manifested in: the diagnosis effectively realizing motor electric fault, reliability is good.
Accompanying drawing explanation
Fig. 1 is the process flow diagram of Method of Motor Fault Diagnosis based on horse field system and Box-Cox conversion of the present invention.
Fig. 2 is motor when normally working, the histogram of mahalanobis distance.
Fig. 3 is the mahalanobis distance of motor when normally working when adopting Box-Cox conversion, the relation curve of maximum likelihood function f (y, λ) and conversion parameter lambda.
Fig. 4 be the mahalanobis distance of motor when normally working after Box-Cox conversion, the normal probability plot of sample data.
Embodiment
Below in conjunction with accompanying drawing, the invention will be further described.
With reference to Fig. 1 ~ Fig. 4, a kind of Method of Motor Fault Diagnosis converted based on horse field system and Box-Cox, comprises the following steps:
Step one, in the training stage, gather the vibration signal of motor faulty motor when normally working and under unbalanced input voltages respectively, wherein fault-signal includes three from small to large, the signal that fault level increases successively, represents respectively by fault 1, fault 2 and fault 3.
Step 2, to the vibration signal collected, carry out calculation of characteristic parameters, build the characteristic data set of a higher-dimension.Specifically vibration signal is carried out to the calculating of time domain charactreristic parameter, include valid value, peak value, kurtosis, measure of skewness, crest factor, nargin coefficient, impulse ratio, form factor and root mean square; Use wavelet package transforms that vibration signal is decomposed 2 layers simultaneously, calculate 4 characteristic parameters about energy in wavelet packet coefficient.These 13 parameters use alphabetical A respectively
1, A
2..., A
13represent.
Step 3, calculate characteristic parameter in step 2, characteristic parameter when normally working with motor is benchmark, uses the based Robust Design theory in the system of horse field, and use two-level orthogonal array to test, result is as shown in table 1, extracts useful characteristic parameter A
1, A
2, A
4, A
5, A
9, A
10.The useful feature parameter extracted is used to recalculate mahalanobis distance.Mahalanobis distance when normally working for motor, its data are as shown in table 2, and the column distribution plan of data as shown in Figure 2, can see these data not Gaussian distributed.
Table 1
Table 2
Step 4, the mahalanobis distance that step 3 is obtained, use Box-Cox conversion, be the data of Gaussian distribution by the data transformation of non-gaussian distribution, then the character of Gaussian distribution is used, namely the Data distribution8 of 95.4% has in the scope within 2 standard deviations at distance average, and by inverse Box-Cox conversion, determine the mahalanobis distance scope of machine operation under different health status.When normally working for motor, find best conversion parameter lambda=0.39 by maximum likelihood function, as shown in Figure 3, it is-0.10 that the mahalanobis distance of non-gaussian distribution is transformed into a class mean, and variance is the data of 0.56 Gaussian distribution, as shown in Figure 4.Use the character of Gaussian distribution, namely about have the Data distribution8 of 95.4% in the scope of [-1.22,1.02], then by inverse Box-Cox conversion,
The mahalanobis distance scope obtained when motor normally works is [0.19,2.36].
With reference to same method, can determine to break down 1 respectively when motor, when fault 2 and fault 3, corresponding mahalanobis distance scope is [5.59,19.70] respectively, [25.32,51.35] and [62.95,101.28].
Step 5, at test phase, to the motor of unknown duty through vibration signals collecting, useful feature parameter calculates and mahalanobis distance calculates, and by from set up in training process, corresponding to the contrast of the mahalanobis distance scope of the different health status of motor, the health status of testing of electric motors is judged, thus realizes the fault diagnosis of motor.
Claims (3)
1., based on the Method of Motor Fault Diagnosis that horse field system and Box-Cox convert, it is characterized in that: said method comprising the steps of:
Step one, gather the vibration signal of motor when normally working and break down;
Step 2, the vibration signal obtained step one carry out the calculating of characteristic parameter, specifically: calculating vibration signal being carried out to time domain charactreristic parameter, simultaneously to vibration signal by wavelet package transforms, carry out the calculating of characteristic parameter on time-frequency domain;
Step 3, the characteristic parameter that step 2 is obtained, the characteristic parameter obtained when normally working with motor is benchmark, uses the based Robust Design theory in the system of horse field to select useful characteristic parameter, and recalculates mahalanobis distance based on these characteristic parameters;
Step 4, to the non-negative that step 3 obtains, the mahalanobis distance of non-gaussian distribution, use Box-Cox conversion, convert the data of Gaussian distribution to, utilize the character of Gaussian distribution and inverse Box-Cox conversion, determine the mahalanobis distance scope of machine operation when different health status;
Sample data x=[the x be made up of mahalanobis distance
1, x
2..., x
n], N is the number of mahalanobis distance sample, i-th data x in x
ithe data y obtained after Box-Cox conversion
i(λ) calculate by following formula:
Wherein convert parameter lambda to be estimated by the maximum likelihood function of following formula:
In formula, y=[y
1, y
2..., y
n],
before and after Box-Cox converts, the distribution situation of data can be undertaken judging or verifying by histogram or normal probability plot;
Then use the character of Gaussian distribution, obtain the upper and lower limit of the data of the rear Gaussian distributed of Box-Cox conversion, and determine the mahalanobis distance scope of machine operation under different health status by the inverse Box-Cox conversion of following formula:
Step 5, to the motor of unknown duty through vibration signals collecting, useful feature parameter calculates and mahalanobis distance calculates, the mahalanobis distance scope determined in contrast step 4, judges the carrying out of testing of electric motors health status, thus realizes the fault diagnosis of motor.
2. as claimed in claim 1 based on the Method of Motor Fault Diagnosis that horse field system and Box-Cox convert, it is characterized in that: in described step 3, use the based Robust Design theory in the system of horse field to select useful characteristic parameter, comprise the following steps:
Step 3.1, to set after a stack features standard parameter characteristic of correspondence vector as z
i, its dimension is n, the mahalanobis distance MD that these group data are corresponding
ifor
Wherein, C
-1for the inverse matrix of correlation matrix,
for character vector z
itransposition;
Step 3.2, use Taguchi's method are optimized characteristic parameter and choose;
Quantity n according to characteristic parameter in proper vector selects two-level orthogonal array L
p(t
c), wherein p is test number (TN), and t=2 is number of levels, and definition level '+' is " using this characteristic parameter ", and level '-' is " not using this characteristic parameter ", and c is columns, represents and allows to arrange maximum characteristic parameter numbers; According to 13 characteristic parameter selection orthogonal arrage L that step 2 builds
16(2
15); To a jth characteristic parameter, if
with
be respectively the signal to noise ratio (S/N ratio) average of this characteristic parameter under level '+' and level '-';
represent and adopt this characteristic parameter on average to detect effect to fault,
represent and do not adopt this characteristic parameter on average to detect effect to fault, therefore select satisfied
useful feature parameter.
3. as claimed in claim 1 or 2 based on the Method of Motor Fault Diagnosis that horse field system and Box-Cox convert, it is characterized in that: in described step 2, vibration signal is carried out to the calculating of time domain charactreristic parameter, include valid value, peak value, kurtosis, measure of skewness, crest factor, nargin coefficient, impulse ratio, form factor and root mean square; Use wavelet package transforms that vibration signal is decomposed 2 layers simultaneously, calculate 4 characteristic parameters about energy in wavelet packet coefficient.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410452287.1A CN104215905A (en) | 2014-09-05 | 2014-09-05 | Motor fault diagnosis method based on Mahalanobis-Taguchi system and Box-Cox transformation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410452287.1A CN104215905A (en) | 2014-09-05 | 2014-09-05 | Motor fault diagnosis method based on Mahalanobis-Taguchi system and Box-Cox transformation |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104215905A true CN104215905A (en) | 2014-12-17 |
Family
ID=52097606
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410452287.1A Pending CN104215905A (en) | 2014-09-05 | 2014-09-05 | Motor fault diagnosis method based on Mahalanobis-Taguchi system and Box-Cox transformation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104215905A (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105300692A (en) * | 2015-08-07 | 2016-02-03 | 浙江工业大学 | Bearing fault diagnosis and prediction method based on extended Kalman filtering algorithm |
CN106021719A (en) * | 2016-05-19 | 2016-10-12 | 浙江工业大学 | Unscented Kalman filter algorithm-based method for predicting residual life of bearing |
CN109477869A (en) * | 2016-07-25 | 2019-03-15 | 三菱电机株式会社 | The diagnostic device of motor |
CN110187206A (en) * | 2019-05-22 | 2019-08-30 | 中国人民解放军国防科技大学 | Fault detection method for suspension system in non-Gaussian process under complex working condition |
CN111553048A (en) * | 2020-03-23 | 2020-08-18 | 中国地质大学(武汉) | Method for predicting sintering process operation performance based on Gaussian process regression |
CN112101458A (en) * | 2020-09-16 | 2020-12-18 | 河海大学常州校区 | Taguchi function-signal-to-noise ratio-based characteristic measurement method and device |
CN112380932A (en) * | 2020-11-02 | 2021-02-19 | 上海三菱电梯有限公司 | Vibration signal characteristic value selection method and elevator health state evaluation or fault diagnosis method |
CN113609964A (en) * | 2021-08-03 | 2021-11-05 | 西安市双合软件技术有限公司 | Motor abnormal vibration early warning method and system |
-
2014
- 2014-09-05 CN CN201410452287.1A patent/CN104215905A/en active Pending
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105300692B (en) * | 2015-08-07 | 2017-09-05 | 浙江工业大学 | A kind of bearing failure diagnosis and Forecasting Methodology based on expanded Kalman filtration algorithm |
CN105300692A (en) * | 2015-08-07 | 2016-02-03 | 浙江工业大学 | Bearing fault diagnosis and prediction method based on extended Kalman filtering algorithm |
CN106021719A (en) * | 2016-05-19 | 2016-10-12 | 浙江工业大学 | Unscented Kalman filter algorithm-based method for predicting residual life of bearing |
CN109477869B (en) * | 2016-07-25 | 2021-02-26 | 三菱电机株式会社 | Diagnostic device for motor |
CN109477869A (en) * | 2016-07-25 | 2019-03-15 | 三菱电机株式会社 | The diagnostic device of motor |
CN110187206B (en) * | 2019-05-22 | 2021-11-02 | 中国人民解放军国防科技大学 | Fault detection method for suspension system in non-Gaussian process under complex working condition |
CN110187206A (en) * | 2019-05-22 | 2019-08-30 | 中国人民解放军国防科技大学 | Fault detection method for suspension system in non-Gaussian process under complex working condition |
CN111553048A (en) * | 2020-03-23 | 2020-08-18 | 中国地质大学(武汉) | Method for predicting sintering process operation performance based on Gaussian process regression |
CN111553048B (en) * | 2020-03-23 | 2023-09-22 | 中国地质大学(武汉) | Method for predicting operation performance of sintering process based on Gaussian process regression |
CN112101458A (en) * | 2020-09-16 | 2020-12-18 | 河海大学常州校区 | Taguchi function-signal-to-noise ratio-based characteristic measurement method and device |
CN112101458B (en) * | 2020-09-16 | 2024-04-19 | 河海大学常州校区 | Characteristic measurement method and device based on field function-signal-to-noise ratio |
CN112380932A (en) * | 2020-11-02 | 2021-02-19 | 上海三菱电梯有限公司 | Vibration signal characteristic value selection method and elevator health state evaluation or fault diagnosis method |
CN112380932B (en) * | 2020-11-02 | 2022-10-14 | 上海三菱电梯有限公司 | Vibration signal characteristic value selection method and elevator health state evaluation or fault diagnosis method |
CN113609964A (en) * | 2021-08-03 | 2021-11-05 | 西安市双合软件技术有限公司 | Motor abnormal vibration early warning method and system |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104215905A (en) | Motor fault diagnosis method based on Mahalanobis-Taguchi system and Box-Cox transformation | |
CN102288857B (en) | Fault arc identification and detection method and detection protection device | |
CN106443316B (en) | Multi-information detection method and device for deformation state of power transformer winding | |
CN106338385B (en) | A kind of rotary machinery fault diagnosis method based on unusual spectral factorization | |
CN101819253B (en) | Probabilistic neural network-based tolerance-circuit fault diagnosis method | |
CN103761372B (en) | A kind of multi-electrical level inverter Fault Diagnosis Strategy based on pivot analysis with more classification Method Using Relevance Vector Machines | |
CN105129109B (en) | Method for evaluating health of aircraft aileron actuator system based on multi-fractal theory and self-organizing map (SOM) network | |
CN102222151B (en) | Analog circuit fault prediction method based on ARMA (Autoregressive Moving Average) | |
CN108229544A (en) | A kind of NPC three-level photovoltaic inverters open fault diagnostic method | |
CN104933012B (en) | A kind of measured deviation On-line Fault recognition methods of traction substation transformer | |
CN103064008B (en) | A kind of Nolinear analog circuit soft fault diagnostic method based on Hilbert-Huang transform | |
CN104020401B (en) | The appraisal procedure of transformer insulated heat ageing state based on cloud models theory | |
CN103760464B (en) | Solve and the low current neutral grounding system fault route selecting method of SVM based on analysis diagram | |
CN103983452B (en) | Utilize the method that hybrid domain characteristic vector and grey correlation analysis carry out Fault Pattern Recognition to epicyclic gearbox | |
CN103995237A (en) | Satellite power supply system online fault diagnosis method | |
CN102636740A (en) | Method for predicting faults of power electronic circuit based on FRM-RVM (fuzzy rough membership-relevant vector machine) | |
CN103077402A (en) | Transformer partial-discharging mode recognition method based on singular value decomposition algorithm | |
CN105372531A (en) | Transformer insulation thermal aging parameter correlation calculation method based on Weibull distribution model | |
CN103745119A (en) | Oil-immersed transformer fault diagnosis method based on fault probability distribution model | |
CN109472288A (en) | A kind of pump-storage generator vibration hybrid feature extraction and classification method | |
CN103439631A (en) | Method and system for detecting corrosion state of grounding grid | |
CN101782625B (en) | Power electronic system fault diagnostic method based on Gradation-boosting algorithm | |
CN104793124A (en) | Switched circuit fault diagnosis method based on wavelet transformation and ICA (independent component analysis) feature extraction | |
CN105510729A (en) | Overheating fault diagnosis method of transformer | |
CN105182141A (en) | Transformer fault detection method and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20141217 |