CN110057581A - Rotary machinery fault diagnosis method based on interval type reliability rule-based reasoning - Google Patents
Rotary machinery fault diagnosis method based on interval type reliability rule-based reasoning Download PDFInfo
- Publication number
- CN110057581A CN110057581A CN201910354802.5A CN201910354802A CN110057581A CN 110057581 A CN110057581 A CN 110057581A CN 201910354802 A CN201910354802 A CN 201910354802A CN 110057581 A CN110057581 A CN 110057581A
- Authority
- CN
- China
- Prior art keywords
- rule
- fault
- interval type
- reliability
- interval
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/02—Gearings; Transmission mechanisms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/02—Gearings; Transmission mechanisms
- G01M13/028—Acoustic or vibration analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N5/00—Computing arrangements using knowledge-based models
- G06N5/04—Inference or reasoning models
- G06N5/048—Fuzzy inferencing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Software Systems (AREA)
- Theoretical Computer Science (AREA)
- Computational Linguistics (AREA)
- Fuzzy Systems (AREA)
- Artificial Intelligence (AREA)
- Automation & Control Theory (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Acoustics & Sound (AREA)
- Complex Calculations (AREA)
Abstract
The present invention relates to a kind of rotary machinery fault diagnosis methods based on interval type reliability rule-based reasoning.The present invention can classify to the fault signature data obtained under various faults mode, and construct the interval type reliability rule base of Fault characteristic parameters and fault type；It is online to obtain feature input parameter and reference values match degree, and computation rule activates weight；Interval type reliability is modified with activation weight, obtains new interval type reliability；It utilizesDempsterThe interval type reliability that these are activated by rule merges to obtain new interval type reliability, and according to the decision rule under the evidence of section, obtains fault type corresponding to the online fault signature.The present invention uses interval type reliability, the support that fault mode occurs for description fault characteristic signals, and the obtained failure result of decision contains more information capacities, is more favorable for policymaker and judges.
Description
Technical field
The present invention relates to a kind of rotary machinery fault diagnosis methods based on interval type reliability rule-based reasoning, belong to mechanical event
Barrier diagnosis and maintenance area.
Background technique
It is be unable to do without rotating machinery in modernization industry production, equipment is reduced by fault diagnosis technology to the greatest extent
Failure, reduce loss, this for reduce cost of equipment maintenance and promoted its working efficiency play the role of it is vital.Especially
It is to the rotating machinery of some large sizes, structure is complicated, and maintenance cost is high, long period, once serious failure occurs,
Easily cause serious harm.The generation for judging failure ahead of time using fault diagnosis technology thus, for the peace of lifting means
Full property has important practical significance.
Modernization industry production in be unable to do without mechanical equipment operation, reduced to the greatest extent before device fails because
Probability caused by equipment fault reduces loss, this is played for reducing cost of equipment maintenance and promoting its working efficiency to pass
Important role, the especially rotating machinery to some large sizes, structure is complicated, and maintenance cost is high, long period, once
Serious failure occurs, it is easy to cause serious harm.The generation of judgement failure ahead of time thus, for safeguarding the people
The security of the lives and property have great meaning.
Two critical issues are faced in the fault diagnosis of rotating machinery, first is that due to measurement environment in interference or
The real-time change of person's equipment itself operating status, so that the fault characteristic signals collected for troubleshooting usually all have
Stronger uncertainty also has led to the uncertainty for breakdown judge result.Interval type reliability, description event are used thus
The support that fault mode occurs for barrier characteristic signal, can more objectively describe the uncertainty of failure generation in this way；Two
It is when being differentiated using various characteristic signals to fault mode, due to probabilistic influence, what single feature provided sentences
Disconnected confidence results are necessarily inaccurate and unilateral, so needing using corresponding information fusion method to each single source feature
The section reliability that signal provides is merged, and the accuracy of failure decision is increased.
Summary of the invention
The purpose of the present invention is to propose to a kind of rotary machinery fault diagnosis method based on interval type reliability rule-based reasoning, from
Sampling analysis is carried out to corresponding fault signature data under different types of faults mode, and is converted into each fault type
The reliability of generation does the decision that is out of order according to the confidence interval obtained after fusion by the fusion of multi-source confidence interval.
Rotating machinery method for diagnosing faults proposed by the present invention based on section evidence fusion, including following step
It is rapid:
Rotary machinery fault diagnosis method based on interval type reliability rule-based reasoning, comprising the following steps:
(1) framework of identification Θ={ ζ of the rotary machinery fault diagnosis method based on interval type reliability rule-based reasoning is set1,
ζ2,ζ3, wherein ζ1Indicate the imbalance fault in rotating machinery, ζ2Indicate the unbalanced fault in rotating machinery, ζ3Indicate rotation
Pedestal looseness fault in machinery.
(2) set the revolving speed of rotating machinery as p, unit: rev/min (r/min), here p ∈ [1000r/min,
3000r/min], being mounted on the vibration acceleration sensor on the positions such as pedestal or bracket with Δ t ∈ [16s, 48s] is interval,
Time domain vibration acceleration signal can be acquired under various fault modes, obtained signal sequence ω (r), r=1,2 ..., n, n are
The number in sampling period.
(3) the time domain vibration signal sequence ω (r) obtained in step (2) is subjected to Fast Fourier Transform (FFT), is transformed to phase
Then the frequency-domain spectrum figure answered therefrom is chosen the amplitude of 1 times of fundamental frequency, 2 times of fundamental frequencies and 3 times of fundamental frequencies as Fault characteristic parameters, is divided
X is not denoted as it1(r)、x2(r)、x3(r)。
(4) interval type reliability rule base is established, Expressive Features parametric variable x is used for1、x2、x3It is non-between fault mode
Linear mapping relation, wherein kth rule RkRepresentation it is as follows:
In formula (1),For input characteristic parameter variable xi(i=1,2,3) input reference
Set, element therein meetqiIndicate the value number of character pair parameter reference values,WithFor x under various fault modesiMinimum and maximum value, βj,k(j=1,2,3, k=1,2 ..., L) it is to fault mode
ζj, the interval type reliability of j=1,2,3 imparting:
βj,k=[aj,bj] (2)
Wherein ajIndicate the minimum value of interval efficiency reliability, bjIndicate the maximum value of interval efficiency reliability, and 0≤aj≤bj
≤ 1, and meet constraint conditionAnd
(5) for the characteristic parameter sample x of r-th of sampling period acquisitioni(r) it brings into and defines kth rule in step (4)
In, obtain xiWith kth rule input referenceMatching degree:
1. whenOrWhen, xiForWithMatching degreeValue is 1, for other references
The matching degree value of value is 0.
2. whenWhen, xiForWithMatching degreeValue difference
Are as follows:
Matching degree value for other reference values is 0.
(6) according to formula (3), the activation weight of every rule in computation interval confidence rule base:
In formula (4),L is the regular number in rule base, θkIt is 0≤θ for initial rules weightk≤ 1, δi
For input characteristic parameter initial attribute weight, 0≤δi≤1。
(7) Dempster rule of combination is used, the rule being activated in the confidence rule base of section is merged to obtain every
The section certainty value of one fault mode, the specific steps are as follows:
The activation weight that (7-1) is obtained according to formula (4), to the interval type reliability β in kth rulej,kDiscount is carried out,
Interval type reliability after discount taken:
mj,k(ζj)=wkβj,k, j=1,2,3, k=1,2 ..., L
mΘ,k(Θ)=1-wk (5)
Wherein, mj,kFor j-th of interval type reliability after discount, mΘ,kIt is expressed as the interval type reliability of complete or collected works' element.
Interval type reliability (the m corresponding with the 1st, 2 rules of (7-2) for acquisition in step (7-1)j,1,mΘ,1) and
(mj,2,mΘ,2), they are merged using following Dempster rule of combination, obtains fused interval type reliability are as follows:
The combinatorial operation of the above interval type reliability can be realized by fmincon majorized function in MATLAB.
The result that (7-3) is obtained by formula (6)(mj,3,
mΘ,3) continue in step (7-2) Dempster rule of combination fusion, obtain new interval type reliabilityIt successively merges, final interval type can be obtained
Reliability Wherein in k expression system
Regular number.
(8) the interval type reliability judgement set type result of decision combined according to Dempster rule, meets following
Two conditions can determine that fault type is ζj, j=1,2,3:
(8-1)The left and right endpoint of confidence interval is respectively greater than the interval type reliability of other fault modes
Left and right endpoint.
(8-2)Right endpoint be respectively less than 0.3.
Rotary machinery fault diagnosis method proposed by the present invention based on interval type reliability rule-based reasoning, to different faults mould
The fault signature monitoring data obtained under formula are divided, and the mapping relations rule of Fault characteristic parameters and fault type is constructed
Library；It will acquire feature input parameter to be matched with reference value, and calculate activation weight；Discount is carried out to interval type reliability, is obtained
To new interval type reliability；It merges to obtain new interval type letter using the interval type reliability that these are activated by Dempster rule
Degree, and according to the decision rule under the evidence of section, obtain fault type corresponding to the online fault signature.
Beneficial effects of the present invention:
One, due to input information there is uncertainty, also have led to the uncertainty for breakdown judge result,
The present invention uses interval type reliability, the support that fault mode occurs for description fault characteristic signals, obtained failure thus
The result of decision contains more information capacities, is more favorable for policymaker and judges.
Two, when due to being differentiated using various characteristic signals, there is probabilistic influence, so single feature
The judgement provided can not be for the result is that reliable, so the method that the present invention uses the fusion of corresponding information, for single
Characteristic signal carry out reliability fusion, to increase the accuracy of decision.
Detailed description of the invention
Fig. 1 is the flow diagram of the method for the present invention.
Fig. 2 is rotor fault diagnosis system structure chart in the embodiment of the method for the present invention.
Specific embodiment
Rotary machinery fault diagnosis method based on interval type reliability rule-based reasoning, comprising the following steps:
(1) framework of identification Θ={ ζ of the rotary machinery fault diagnosis method based on interval type reliability rule-based reasoning is set1,
ζ2,ζ3, wherein ζ1Indicate the imbalance fault in rotating machinery, ζ2Indicate the unbalanced fault in rotating machinery, ζ3Indicate rotation
Pedestal looseness fault in machinery.
(2) set the revolving speed of rotating machinery as p, unit: rev/min (r/min), here p ∈ [1000r/min,
3000r/min], being mounted on the vibration acceleration sensor on the positions such as pedestal or bracket with Δ t ∈ [16s, 48s] is interval,
Time domain vibration acceleration signal can be acquired under various fault modes, obtained signal sequence ω (r), r=1,2 ..., n, n are
The number in sampling period.
(3) the time domain vibration signal sequence ω (r) obtained in step (2) is subjected to Fast Fourier Transform (FFT), is transformed to phase
Then the frequency-domain spectrum figure answered therefrom is chosen the amplitude of 1 times of fundamental frequency, 2 times of fundamental frequencies and 3 times of fundamental frequencies as Fault characteristic parameters, is divided
X is not denoted as it1(r)、x2(r)、x3(r)。
(4) interval type reliability rule base is established, Expressive Features parametric variable x is used for1、x2、x3It is non-between fault mode
Linear mapping relation, wherein kth rule RkRepresentation it is as follows:
In formula (1),For input characteristic parameter variable xi(i=1,2,3) input reference
Set, element therein meetqiIndicate the value number of character pair parameter reference values,WithFor x under various fault modesiMinimum and maximum value；βj,k(j=1,2,3, k=1,2 ..., L) it is to fault mode
ζj, the interval type reliability of j=1,2,3 imparting:
βj,k=[aj,bj] (2)
Wherein ajIndicate the minimum value of interval efficiency reliability, bjIndicate the maximum value of interval efficiency reliability, and 0≤aj≤bj
≤ 1, and meet constraint conditionAnd
For the ease of the understanding of interval type reliability rule base, illustrate here, it is assumed thatIn each reference value value
ForIn each reference value value beIn each reference value value beThe then raw L=64 rule of common property,
Set initial θk=1, δ1=δ2=δ3=1, part of rule format is as follows:
R1If: x1=0.05 and x2=0.035 and x3=0.28, then β1,1=[0.5632,0.9672], β2,1=
[0.1246,0.1365],β3,1=[0.1126,0.1365]
R2If: x1=0.05 and x2=0.125 and x3=0.42, then β1,2=[0.5631,0.9463], β2,2=
[0.1258,0.2367],β3,2=[0.1248,0.3365]
R3If: x1=0.05 and x2=0.25 and x3=0.67, then β1,3=[0.6321,0.9532], β2,3=
[0.1124,0.3236)],β3,3=[0.1214,0.2362]
......
R62If: x1=0.38 and x2=0.035 and x3=0.28, then β1,62=[0.1129,0.2563], β2,62=
[0.2145,0.4259],β3,62=[0.4569,0.9653]
R63If: x1=0.05 and x2=0.45 and x3=0.42, then β1,63=[0.1127,0.3896], β2,63=
[0.1116,0.3256],β3,63=[0.6563,0.9782]
R64If: x1=0.05 and x2=0.46 and x3=0.48, then β1,64=[0.3456,0.4569], β2,64=
[0.3213,0.4452],β3,64=[0.4563,0.9863]
(5) for the characteristic parameter sample x of r-th of sampling period acquisitioni(r) it brings into and defines kth rule in step (4)
In, obtain xiWith kth rule input referenceMatching degree:
1. whenOrWhen, xiForWithMatching degreeValue is 1, for other references
The matching degree value of value is 0.
2. whenWhen, xiForWithMatching degreeValue difference
Are as follows:
Matching degree value for other reference values is 0.
(6) according to formula (3), the activation weight of every rule in computation interval confidence rule base:
In formula (4),L is the regular number in rule base, θkIt is 0≤θ for initial rules weightk≤ 1,
δiFor input characteristic parameter initial attribute weight, 0≤δi≤1。
For ease of understanding, here for example, by taking the model in step (4) as an example, each input variable computation rule
Weight is activated, using the activation weight w of every rule in formula (3) and formula (4) computation interval type reliability algorithmk, it is assumed that
Input variable x1(r)=0.1663, x2(r)=0.1612, x3(r)=0.1221 regular R can, be activated21、R22、R25、R26、R37、
R38、R41、R42, can be calculated activation weight wkAre as follows: w1=0.075, w2=0.64, w3=0.0171, w4=0.16, w5=
0.0092、w6=0.0838, w7=0.011, w8=0.0035, remaining is 0, that is, has activated 8 in interval type reliability rule base
Rule.
The activation weight that (7-1) is obtained according to formula (4), to the interval type reliability β in kth rulej,kDiscount is carried out,
Interval type reliability after discount taken:
mj,k(ζj)=wkβj,k, j=1,2,3, k=1,2 ..., L
mΘ,k(Θ)=1-wk (5)
Wherein, mj,kFor j-th of interval type reliability after discount, mΘ,kThe interval type reliability L of complete or collected works' element is expressed as rule
The then number in library.
It for the ease of the understanding to formula (5), is illustrated here, it is assumed that w1=0.075, w2=0.64, w3=
0.0171, w4=0.16, w5=0.0092, w6=0.0838, w7=0.011, w8=0.0035；β1,1=[0.5632,
0.9672],
β2,1=[0.1246,0.1365], β3,1=[0.1126,0.1256], β1,2=[0.4263,0.9236], β2,2=
[0.1362,0.2423], β3,2=[0.1213,0.2352], β1,3=[0.5631,0.9463], β2,3=[0.1258,
0.2367],β3,3=[0.1248,0.3365], β1,4=[0.6321,0.9532], β2,4=[0.1124,0.3236], β3,4=
[0.1214,0.2362],β1,5=[0.1263,0.9542], β2,5=[0.1236,0.2354], β3,5=[0.1124,
0.2263],β1,6=[0.1213,0.2145], β2,6=[0.4569,0.9653], Β3,6=[0.1123,0.3563], β1,7=
[0.7563,0.9782],β2,7=[0.1563,0.2687], β3,7=[0.1129,0.2563], β1,8=[0.1012,
0.1818],β2,8=[0.5623,0.9016], β3,8=[0.1456,0.1569]；
It can be obtained according to formula (5):
mj,1(ζj)={ [0.0422,0.0725], [0.0093,0.0102], [0.0084,0.0094] }, mj,2(ζj)=
{[0.2728,0.5911],[0.0872,0.1551],[0.0776,0.1505]},mj,3(ζj)
={ [0.0096,0.0162], [0.0022,0.0040], [0.0021,0.0058] }, mj,4(ζj)=
{[0.1011,0.1525],[0.0180,0.0518],[0.0194,0.0378]},
mj,5(ζj)={ [0.0012,0.0088], [0.0011,0.0022], [0.0010,0.0021] }, mj,6(ζj)
={ [0.0102,0.0180], [0.0383,0.0809], [0.0094,0.0299] }, mj,7(ζj)
={ [0.0083,0.0108], [0.0017,0.0030], [0.0012,0.0028] }, mj,8(ζj)
={ [0.0004,0.0006], [0.0020,0.0032], [0.0005,0.0005] }；mΘ,1(Θ)=0.925,
mΘ,2(Θ)
=0.36, mΘ,3(Θ)=0.9829, mΘ,4(Θ)=0.84, mΘ,5(Θ)=0.9908, mΘ,6(Θ)=
0.9162,
mΘ,7(Θ)=0.989, mΘ,8(Θ)=0.9965；
Interval type reliability (the m corresponding with the 1st, 2 rules of (7-2) for acquisition in step (7-1)j,1,mΘ,1) and
(mj,2,mΘ,2), it is combined using following Dempster rule and merges them, obtain fused interval type reliability are as follows:
The combinatorial operation of the above interval type reliability can be realized by fmincon majorized function in MATLAB.
The result that (7-3) is obtained by formula (6)(mj,3,
mΘ,3) continue in step (7-2) Dempster rule of combination fusion, obtain new interval type reliabilityIt successively merges, final interval type can be obtained
Reliability Wherein in k expression system
Regular number.
In order to deepen to provide example here, it is assumed that m to the understanding of formula (6)j,1(ζj)={ [0.0205,0.0669],
[0.014,0.0199], [0.0092,0.0166] }, mθ,1(Θ)=[0.9298,0.9298], mj,2(ζj)={ [0.2061,
0.6330], [0.1040,0.1697], [0.1151,0.1698] }, mθ,2(Θ)=[0.3591,0.3591], mj,3(ζj)=
{ [0.0073,0.1065], [0.0033,0.0039], [0.0021,0.0049] }, mθ,3(Θ)=[0.9829,0.9829],
mj,4(ζj)={ [0.0420,0.1478], [0.0190,0.0383], [0.0213,0.0425] }, mθ,4(Θ)=[0.8438,
0.8438], mj,5(ζj)={ [0.0051,0.0090], [0.0014,0.0033], [0.0010,0.0027] }, mθ,5(Θ)=
[0.9908,0.9908], mj,6(ζj)={ [0.0388,0.0797], [0.0105,0.0199], [0.0096,0.0189] },
mθ,6(Θ)=[0.9162,0.9162], mj,7(ζj)={ [0.0014,0.0019], [0.0003,0.0007], [0.0003,
0.0007] }, mθ,7(Θ)=[0.9978,0.9978], mj,8(ζj)={ [0.0029,0.0052], [0.0118,0.0183],
[0.0030,0.0046] }, mθ,8(Θ)=[0.9796,0.9796]；
It can be obtained according to the Dempster rule of combination in step (7-2) 0.3437], then
By obtained result and (mj,3,mΘ,3) according to the Dempster rule of combination fusion in step (7-2), obtain final interval type
Reliability The wherein regular number in k expression system.
(8) it has been obtained in the fields such as target identification, system evaluation, fault diagnosis based on the decision-making technique of evidential reasoning wide
General application, the interval type reliability judgement set type result of decision combined according to Dempster rule, meets following two
A condition can be obtained interval type reliability and determine fault type ζj, j=1,2,3:
(8-1)The left and right endpoint of confidence interval is respectively greater than the interval type reliability of other fault modes
Left and right endpoint.
(8-2)Right endpoint be respectively less than 0.3.
Below in conjunction with attached drawing, the embodiment of the method for the present invention is discussed in detail:
The flow chart of the method for the present invention is as shown in Figure 1, core is: time-domain signal is obtained by Fast Fourier Transform (FFT)
To frequency-domain spectrum figure, different fundamental frequencies are chosen as fault signature and input parameter；Interval type reliability rule base is constructed to describe to input
The Nonlinear Mapping relationship of variable and fault signature mode；Discount is carried out to obtained interval type reliability and obtains new interval type letter
Degree；The result of decision of interval type reliability fusion is obtained using Dempster rule, accurately describes failure using decision rule
The type of generation；
1, framework of identification Θ={ ζ of the rotary machinery fault diagnosis method based on interval type reliability rule-based reasoning is set1,
ζ2,ζ3, wherein ζ1Indicate the imbalance fault in rotating machinery, ζ2Indicate the unbalanced fault in rotating machinery, ζ3Indicate rotation
Pedestal looseness fault in machinery.
2, method as shown in Figure 2 acquires data, sets the revolving speed of rotating machinery as 1500r/min, is mounted on base
The time interval Δ t=16s of vibration acceleration sensor on the positions such as seat or bracket, can acquire under three kinds of fault modes
Time domain vibration acceleration signal obtains signal sequence ω (r), r=1, and 2 ..., n, n are the number in sampling period.
3, the time domain vibration signal sequence ω (r) obtained in step (2) is subjected to Fast Fourier Transform (FFT), be transformed to corresponding
Frequency-domain spectrum figure, then therefrom choose 1 times of fundamental frequency, the amplitude of 2 times of fundamental frequencies and 3 times of fundamental frequencies is used as Fault characteristic parameters, distinguish
It is denoted as x1(r)、x2(r)、x3(r), r=1,2 ..., n are acquired n=40 times altogether here in order to specific description.
4, interval type reliability rule base is established, Expressive Features parametric variable x is used for1、x2、x3It is non-between fault mode
Linear mapping relation:
Choose the semantic values of input variable:
x1(r) fuzzy variable describes are as follows: minimum (PS1), very little (PM1), general (PN1), very big (PQ1),
x2(r) fuzzy variable describes are as follows: minimum (PS2), very little (PM2), general (PN2), very big (PQ2),
x3(r) fuzzy variable describes are as follows: minimum (PS3), very little (PM3), general (PN3), very big (PQ3), specific assignment
As shown in table 1- table 3:
1 x of table1(r) semantic values and reference value
2 x of table2(r) semantic values and reference value
3 x of table3(r) semantic values and reference value
Building interval type reliability rule base is defined according to 3 reference value of table 1- table, and meets constraint conditionAndExist respectivelyOne element of middle extraction is as x1(r),x2(r),x3(r) reference value, thus group
Composition rule, total to can produce L=4 × 4 × 4 rule rule bases building as shown in table 4:
4 interval type reliability rule base of table
5, in conjunction with given input x1、x2、x3Interval type reliability after obtaining discount with step (7), can be obtained final area
Between type reliability；
Firstly, the matching degree of input characteristic parameter and kth rule input reference is obtained according to formula (3), further according to formula
(4) the activation weight of every rule in the confidence rule base of section is obtained, the interval type letter after discount is then obtained by step (7)
Degree, finally, the rule of fusion activation, final interval type reliability can be obtained using Dempster rule of combination As shown in table 5- table 7.
5 x of table1(r) sample interval merges reliability
6 x of table2(r) sample interval merges reliability
7 x of table3(r) sample interval merges reliability
6, after obtaining table 5- table 7, according to the method for the present invention the step of (8) section fusion decision rule can determine whether it is true
Real fault type.
Claims (1)
1. the rotary machinery fault diagnosis method based on interval type reliability rule-based reasoning, it is characterised in that this method includes following step
It is rapid:
(1) framework of identification Θ={ ζ of the rotary machinery fault diagnosis method based on interval type reliability rule-based reasoning is set1,ζ2,
ζ3, wherein ζ1Indicate the imbalance fault in rotating machinery, ζ2Indicate the unbalanced fault in rotating machinery, ζ3Indicate whirler
Pedestal looseness fault in tool；
(2) revolving speed of rotating machinery is set as p, unit: rev/min, p ∈ [1000r/min, 3000r/min] here, peace
Vibration acceleration sensor on pedestal or backing positions is interval with Δ t ∈ [16s, 48s], under various fault modes
Time domain vibration acceleration signal can be acquired, is obtained signal sequence ω (r), r=1,2 ..., n, n are the number in sampling period；
(3) the time domain vibration signal sequence ω (r) obtained in step (2) is subjected to Fast Fourier Transform (FFT), be transformed to corresponding
Then frequency-domain spectrum figure is therefrom chosen the amplitude of 1 times of fundamental frequency, 2 times of fundamental frequencies and 3 times of fundamental frequencies as Fault characteristic parameters, is remembered respectively
For x1(r)、x2(r)、x3(r)；
(4) interval type reliability rule base is established, Expressive Features parametric variable x is used for1、x2、x3It is non-linear between fault mode
Mapping relations, wherein kth rule RkRepresentation it is as follows:
In formula (1),For input characteristic parameter variable xiInput reference set, element therein
MeetqiIndicate the value number of character pair parameter reference values,WithFor various fault modes
Lower xiMinimum and maximum value, βj,kFor to fault mode ζjThe interval type reliability of imparting:
βj,k=[aj,bj] (2)
Wherein ajIndicate the minimum value of interval efficiency reliability, bjIndicate the maximum value of interval efficiency reliability, and 0≤aj≤bj≤ 1, and
Meet constraint conditionAnd
(5) for the characteristic parameter sample x of r-th of sampling period acquisitioni(r) it brings into step (4) and defines in kth rule, obtain
Obtain xiWith kth rule input referenceMatching degree:
1. whenOrWhen, xiForWithMatching degreeValue is 1, for other reference values
Matching degree value is 0；
2. whenWhen, xiForWithMatching degreeValue be respectively as follows:
Matching degree value for other reference values is 0；
(6) according to formula (3), the activation weight of every rule in computation interval confidence rule base:
In formula (4),L is the regular number in rule base, θkIt is 0≤θ for initial rules weightk≤ 1, δiIt is defeated
Enter characteristic parameter initial attribute weight, 0≤δi≤1；
(7) Dempster rule of combination is used, the rule being activated in the confidence rule base of section is merged to obtain each event
The section certainty value of barrier mode, the specific steps are as follows:
The activation weight that (7-1) is obtained according to formula (4), to the interval type reliability β in kth rulej,kDiscount is carried out, is obtained
Interval type reliability after discount:
mj,k(ζj)=wkβj,k
mΘ,k(Θ)=1-wk (5)
Wherein, mj,kFor j-th of interval type reliability after discount, mΘ,kIt is expressed as the interval type reliability of complete or collected works' element；
Interval type reliability (the m corresponding with the 1st, 2 rules of (7-2) for acquisition in step (7-1)j,1,mΘ,1) and (mj,2,
mΘ,2), they are merged using following Dempster rule of combination, obtains fused interval type reliability are as follows:
The result that (7-3) is obtained by formula (6)(mj,3,mΘ,3)
The Dempster rule of combination fusion in step (7-2) is continued with, new interval type reliability is obtainedIt successively merges, final interval type can be obtained
Reliability Wherein in k expression system
Regular number；
(8) the interval type reliability judgement set type result of decision combined according to Dempster rule, meets following two
Condition can determine that fault type is ζj:
(8-1)The left and right endpoint of confidence interval is respectively greater than the left and right of the interval type reliability of other fault modes
Endpoint；
(8-2)Right endpoint be respectively less than 0.3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910354802.5A CN110057581B (en) | 2019-04-29 | 2019-04-29 | Rotary machine fault diagnosis method based on interval type credibility rule reasoning |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910354802.5A CN110057581B (en) | 2019-04-29 | 2019-04-29 | Rotary machine fault diagnosis method based on interval type credibility rule reasoning |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110057581A true CN110057581A (en) | 2019-07-26 |
| CN110057581B CN110057581B (en) | 2020-10-09 |
Family
ID=67321722
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910354802.5A Active CN110057581B (en) | 2019-04-29 | 2019-04-29 | Rotary machine fault diagnosis method based on interval type credibility rule reasoning |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110057581B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110850206A (en) * | 2019-11-13 | 2020-02-28 | 武汉理工大学 | Inverter capacitor aging fault diagnosis method based on confidence rule reasoning |
| CN112101797A (en) * | 2020-09-18 | 2020-12-18 | 中国人民解放***箭军工程大学 | Dynamic fault diagnosis method and system for complex industrial system |
| CN112197968A (en) * | 2020-09-29 | 2021-01-08 | 上海电力大学 | Confidence rule base-based intelligent diagnosis method for dynamic and static rub-impact faults of steam turbine rotor |
| CN112541155A (en) * | 2020-09-30 | 2021-03-23 | 上海电力大学 | Steam turbine rotor unbalance detection method based on confidence rule base |
| CN113569478A (en) * | 2021-07-27 | 2021-10-29 | 郑州铁路职业技术学院 | Rotary machine fault diagnosis method based on element influence degree and confidence rule base |
| CN117055537A (en) * | 2023-10-12 | 2023-11-14 | 深圳丰汇汽车电子有限公司 | Motorcycle electric control system fault diagnosis method based on Internet of things |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1920511A (en) * | 2006-08-01 | 2007-02-28 | 东北电力大学 | Fusion diagnosing method of centrifugal pump vibration accidents and vibration signals sampling device |
| JP2008197062A (en) * | 2007-02-15 | 2008-08-28 | Ricoh Co Ltd | Method of modeling transfer belt drive mechanism and method of supporting design of control system using the same |
| CN101441718A (en) * | 2008-12-19 | 2009-05-27 | 福建三元达通讯股份有限公司 | Sensor information fuse device and method |
| CN102928877A (en) * | 2012-11-14 | 2013-02-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Seismic property combination method based on Dempster/Shafe (D-S) evidence theory |
| CN106650122A (en) * | 2016-12-27 | 2017-05-10 | 宝鸡文理学院 | Equipment variable working condition operation risk evaluation method |
| CN108520266A (en) * | 2018-03-01 | 2018-09-11 | 西北工业大学 | A kind of Time Domain Fusion method for diagnosing faults based on DS evidence theories |
| CN108593286A (en) * | 2018-04-12 | 2018-09-28 | 中国神华能源股份有限公司 | The method for diagnosing faults of rotating machinery and the trouble-shooter of rotating machinery |
| CN109115491A (en) * | 2018-10-16 | 2019-01-01 | 杭州电子科技大学 | A kind of evidence fusion method of Electrical Propulsion Ship shafting propulsion system mechanical fault diagnosis |
-
2019
- 2019-04-29 CN CN201910354802.5A patent/CN110057581B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1920511A (en) * | 2006-08-01 | 2007-02-28 | 东北电力大学 | Fusion diagnosing method of centrifugal pump vibration accidents and vibration signals sampling device |
| JP2008197062A (en) * | 2007-02-15 | 2008-08-28 | Ricoh Co Ltd | Method of modeling transfer belt drive mechanism and method of supporting design of control system using the same |
| CN101441718A (en) * | 2008-12-19 | 2009-05-27 | 福建三元达通讯股份有限公司 | Sensor information fuse device and method |
| CN102928877A (en) * | 2012-11-14 | 2013-02-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Seismic property combination method based on Dempster/Shafe (D-S) evidence theory |
| CN106650122A (en) * | 2016-12-27 | 2017-05-10 | 宝鸡文理学院 | Equipment variable working condition operation risk evaluation method |
| CN108520266A (en) * | 2018-03-01 | 2018-09-11 | 西北工业大学 | A kind of Time Domain Fusion method for diagnosing faults based on DS evidence theories |
| CN108593286A (en) * | 2018-04-12 | 2018-09-28 | 中国神华能源股份有限公司 | The method for diagnosing faults of rotating machinery and the trouble-shooter of rotating machinery |
| CN109115491A (en) * | 2018-10-16 | 2019-01-01 | 杭州电子科技大学 | A kind of evidence fusion method of Electrical Propulsion Ship shafting propulsion system mechanical fault diagnosis |
Non-Patent Citations (5)
| Title |
|---|
| XIAOBIN XU: "Data classification using evidence reasoning rule", 《KNOWLE DGE-BASE D SYSTEMS》 * |
| 刘平: "基于区间值信度结构的故障诊断方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 徐晓滨: "基于证据推理规则的信息融合故障诊断方法", 《控制理论与应用》 * |
| 王晶晶: "基于证据理论与不精确概率的振动故障诊断", 《传感器与微***》 * |
| 靳留乾: "基于证据理论的不确定性推理方法及其应用", 《计算机工程与应用》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110850206A (en) * | 2019-11-13 | 2020-02-28 | 武汉理工大学 | Inverter capacitor aging fault diagnosis method based on confidence rule reasoning |
| CN112101797A (en) * | 2020-09-18 | 2020-12-18 | 中国人民解放***箭军工程大学 | Dynamic fault diagnosis method and system for complex industrial system |
| CN112101797B (en) * | 2020-09-18 | 2023-10-27 | 中国人民解放***箭军工程大学 | Dynamic fault diagnosis method and system for complex industrial system |
| CN112197968A (en) * | 2020-09-29 | 2021-01-08 | 上海电力大学 | Confidence rule base-based intelligent diagnosis method for dynamic and static rub-impact faults of steam turbine rotor |
| CN112541155A (en) * | 2020-09-30 | 2021-03-23 | 上海电力大学 | Steam turbine rotor unbalance detection method based on confidence rule base |
| CN113569478A (en) * | 2021-07-27 | 2021-10-29 | 郑州铁路职业技术学院 | Rotary machine fault diagnosis method based on element influence degree and confidence rule base |
| CN117055537A (en) * | 2023-10-12 | 2023-11-14 | 深圳丰汇汽车电子有限公司 | Motorcycle electric control system fault diagnosis method based on Internet of things |
| CN117055537B (en) * | 2023-10-12 | 2024-01-02 | 深圳丰汇汽车电子有限公司 | Motorcycle electric control system fault diagnosis method based on Internet of things |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110057581B (en) | 2020-10-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110057581A (en) | Rotary machinery fault diagnosis method based on interval type reliability rule-based reasoning | |
| CN104634571B (en) | A kind of Fault Diagnosis of Roller Bearings based on LCD MF | |
| Qian et al. | Edge computing: A promising framework for real-time fault diagnosis and dynamic control of rotating machines using multi-sensor data | |
| CN109506921B (en) | Fault diagnosis and early warning method for rotary machine | |
| CN106006344B (en) | Staircase On-line Fault early warning system and method for diagnosing faults | |
| CN106248356B (en) | A kind of rotary machinery fault diagnosis method based on kurtosis index | |
| CN106441896A (en) | Characteristic vector extraction method for rolling bearing fault mode identification and state monitoring | |
| CN105806613A (en) | Planetary gear case fault diagnosis method based on order complexity | |
| CN106769049A (en) | A kind of Fault Diagnosis of Roller Bearings based on Laplce's score value and SVMs | |
| CN103291600A (en) | Fault diagnosis method for hydraulic pump based on EMD-AR (empirical mode decomposition-auto-regressive) and MTS (mahalanobis taguchi system) | |
| CN105003453A (en) | Online monitoring and fault diagnosis system of mine fan | |
| CN105651504A (en) | Rotary machinery fault feature extraction method based on self-adaptive wavelet energy | |
| CN109115491A (en) | A kind of evidence fusion method of Electrical Propulsion Ship shafting propulsion system mechanical fault diagnosis | |
| CN204113701U (en) | A kind of mine fan on-line monitoring and fault diagnosis system | |
| CN105424366A (en) | Bearing fault diagnosis method based on EEMD adaptive denoising | |
| CN105910823A (en) | Rolling bearing fault diagnosis method | |
| CN104408322B (en) | Rotating mechanical device fault diagnosis method capable of synthesizing multisource fault probability likelihood credibility | |
| CN109765054A (en) | A kind of Fault Diagnosis of Roller Bearings | |
| CN105137324B (en) | A kind of more detection point failure component localization methods based on emulation disaggregated model | |
| CN105675274A (en) | Time-domain parameter and D-S evidence theory-based rotor running state monitoring method | |
| CN105571638A (en) | Machinery device fault combination prediction system and method | |
| CN102123177A (en) | Fault detection system for rotary machine based on network and on-line detection method thereof | |
| CN109976311A (en) | Single order fixed set point control system sensor fault diagnosis method and diagnostic system | |
| CN106203362B (en) | A kind of rotary machinery fault diagnosis method based on pulse index | |
| CN109031118A (en) | A kind of Method of Motor Fault Diagnosis, device and motor |
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 |