A kind of rotating machinery Incipient Fault Diagnosis method optimizing multi-kernel support vector machine based on genetic Annealing
Technical field
The present invention relates to the integral intelligent diagnostic method of a kind of rotating machinery Incipient Fault Diagnosis, particularly relate to a kind of rotating machinery initial failure integral intelligent diagnostic method optimizing multi-kernel support vector machine based on genetic Annealing.
Background technology
The initial failure of rotating machinery is the subject matter threatening rotating machinery safe and reliable operation; if able to the initial failure of rotating machinery is carried out reliable diagnosis; grasp the running status of rotating machinery in advance; the initial failure gradual change gradually that can avoid rotating machinery develops into typical fault, ultimately results in the malignant events such as rotating machinery shutdown suddenly and occurs. In a broad sense, rotating machinery Incipient Fault Diagnosis is for the rotating machinery being currently in operation, the running status that identification rotating machinery is current, finds out its source of trouble and trouble location, analyzes its fault cause, and proposes rational solution. Its research contents relates generally to mechanical kinetics, materialogy, plant equipment measuring technology, Engineering Signal treatment technology, computer technology etc., its main method is based on the theory of mechanical kinetics and materialogy, analyze the failure mechanism of rotating machinery, sensor is adopted to measure the operating state signal that rotating machinery is current, signal is analyzed process, and construct intellectualized algorithm the running status that rotating machinery is current is carried out diagnosis identification.
Tradition is in rotating machinery Incipient Fault Diagnosis process, and technical staff is first with the failure mechanism of the knowledge analysis rotating machinery of mechanical kinetics and materialogy, and utilizes various kinds of sensors to gather the operating state signal of rotating machinery, it is analyzed processing to all kinds of operating state signals obtained subsequently, extracts the principal character index of status signal, artificial intelligence, computer technology is finally utilized to set up fault diagnosis model, the value size according to characteristic index, infer the running status that rotating machinery is current. for rotating machinery, vibration signal is the most direct to the reflection of its running status, therefore, rotating machinery is usually employing vibration signal be analyzed, but the vibration signal intensity of reflection rotating machinery initial failure is relatively small, although fault features information can be collected by modern high-precision sensor, but it is affected by the interference of the environmental factorss such as signal route of transmission, propagation medium, is difficult to accurately extract the vibration performance of rotating machinery initial failure. in the building process of fault diagnosis model, what application was maximum at present is two classification method, k nearest neighbor sorting algorithm, artificial neural network (ArtificialNeuralNetwork, and support vector machine (SupportVectorMachine ANN), SVM) etc., two classification method, k nearest neighbor sorting algorithm is simple, but reliability is low, the reliability of ANN makes moderate progress, but study easily occurred or owes study and be absorbed in local minimum, SVM performance among this is relatively optimum, but parameter selects and the structure of kernel function also has certain randomness, therefore rationally reliable initial failure intelligent diagnostics model also not a duck soup is set up.
Present invention is generally directed to the retrofit that rotating machinery carries out based on the intelligent diagnostics algorithm of support vector machine. Carry out in failure diagnostic process in conventionally employed support vector machine, it is usually the direct support vector machine that the fault features index of acquisition is input to and carries out classification identification, structure for support vector machine, what generally adopt is monokaryon support vector machine, and the parameter for support vector machine generally adopts trial and error procedure or simple optimized algorithm to obtain. In the present invention, structure for support vector machine, on the basis of conventional single-core SVM, a kind of new machine learning model formed by being weighted all kinds of individual core merging, i.e. multinuclear hybrid supporting vector machine (Multi-kernelSupportVectorMachine, MSVM), MSVM inherits extensive Generalization Ability and the self-learning capability of conventional single-core SVM. Rotate in the diagnosis process of machinery initial failure utilizing MSVM, the kernel functional parameter of each individual core and the rationally selection of weight factor thereof will directly affect the identification performance of MSVM, by complex inheritance algorithm (GeneticAlgorithm, it is called for short GA) and simulated annealing (SimulatedAnnealing, it is called for short SA) respective advantage, construct genetic annealing algorithms (GeneticSimulatedAnnealingAlgorithm, it is called for short GA-SA), and then achieve the best to MSVM control parameter preferentially. The rotating machinery Incipient Fault Diagnosis method optimizing multi-kernel support vector machine based on genetic Annealing being consequently formed, the method has good capacity of resisting disturbance, and extensive Generalization Ability is strong, and the initial failure intelligent diagnostics for rotating machinery provides a kind of new approach.
Summary of the invention
Large rotating machinery comprises numerous transmission parts, including each type gear, bearing, axle system etc., once there is typical fault, not only results in the shutdown suddenly of rotating machinery, more seriously the personal safety of line operator can be produced great threat. The typical fault of rotating machinery is a process incremental, that progressively develop. If the commitment that can produce at rotating machinery fault, accomplish discovery as early as possible, as early as possible diagnosis, repair as early as possible, then being avoided that risk factor is accumulative transfinites. It can thus be seen that the Precise Diagnosis of the initial failure of rotating machinery is crucial. And the fault features of rotating machinery is very faint, both difficulty was effectively caught, again difficult effectively identification. Faint in order to solve rotating machinery fault features, the non-linear problem being difficult to identification between fault signature and fault mode, a kind of rotating machinery Incipient Fault Diagnosis method optimizing multi-kernel support vector machine based on genetic Annealing is proposed, be the technical scheme is that and adopt time domain, frequency domain, time-frequency domain signal processing method to extract the time domain in vibration signal, frequency domain, time and frequency domain characteristics index, this time domain, frequency domain, time and frequency domain characteristics index comprehensive are got up, constructs hybrid domain feature set; Adopt gaussian kernel and polynomial kernel to build multi-kernel support vector machine, and complex inheritance algorithm and simulated annealing build genetic annealing algorithms; Optimize all kinds of parameters of multi-kernel support vector machine with genetic annealing algorithms, obtain all kinds of optimized parameters of multi-kernel support vector machine; The multi-kernel support vector machine that the hybrid domain feature set of acquisition is input to genetic annealing algorithms optimization rotates the accurate recognition of machinery initial failure pattern. Specifically comprise the following steps that
Step one, gathers the vibrational state information in rotating machinery running by ICP acceleration transducer, subsequently this vibration signal is carried out noise reduction, filter preprocessing.
Step 2, at rotating machinery each vibration signal pretreated, calculating its time domain, frequency domain, time and frequency domain characteristics index value, constructs hybrid domain feature set.
Step 3, based on gaussian kernel support vector machine and Polynomial SVM, by Gauss karyomerite KRBFWith multinomial overall situation core KpolyCombine, construct multi-kernel support vector machine:
Kmix(xi,xj)=�� KRBF(xi,xj)+(1-��)Kpoly(xi,xj)
In formula: xiAnd xjFor the characteristic vector of the input space, the kernel functional parameter of �� and d respectively gaussian kernel function and Polynomial kernel function, �� (0 < �� < 1) is weight factor.
Step 4, heredity-annealingization algorithm optimization multi-kernel support vector machine, form chromosomal real number gene vectors X with the weight factor �� of multinuclear SVM, penalty parameter c, kernel functional parameter �� and d:
X=[��, c, ��, d]
Utilize multinuclear SVM to the classification accuracy f of hybrid domain feature set training sample to evaluate each chromosomal quality, adopt Metropolis acceptance criterion to determine new chromosome pij' whether enter population of future generation, in Metropolis acceptance criterion, adopt dynamic temperature-regulating strategy, if ti+1=(1-1/i)m��ti, wherein tiMoving back temperature parameter when evolving for i & lt, m �� Z is for controlling the decrease speed of temperature, in order to prevent the optimal solution of current population from losing in the next generation, adopts elite retention strategy, compares the optimal solution x in the fitness f of new explanation x and population*Fitness f*If f is better than f*, then x is replaced respectively with x and f*And f*��
Step 5, is input to the multi-kernel support vector machine after optimization by hybrid domain feature set, utilizes optimum multi-kernel support vector machine automatically to pick out the fault mode of rotating machinery initial failure.
In described step 2, the vibration characteristic signals that rotating machinery produces in running is the Main Basis of reflection running state of rotating machine, by calculating the time domain of vibration signal, frequency domain character index value, can the appearance of quantitatively characterizing rotating machinery initial failure, time domain, frequency domain character index are a lot, at this data analysis by experiment, have chosen and bearing initial failure is characterized obvious 6 temporal signatures indexs and 2 frequency domain character indexs; Described temporal signatures index includes root amplitude, root-mean-square value, variance, average amplitude, kurtosis and maximum; Described frequency domain character index includes means frequency and tune spread;
In order to obtain more fault characteristic information, reflect the initial failure state of rotating machinery comprehensively, also adopt integrated empirical mode decomposition (ensembleempiricalmodedecomposition, EEMD) extracting the time and frequency domain characteristics index in original vibration signal, vibration signal decomposes, through EEMD, the final result obtained and is:
Wherein: cjT () is original vibration signal x (t) limited the Intrinsic Mode component (Intrinsicmodefunction, IMF) obtained from high frequency to low frequency adaptive decomposition, r (t) is the residual components after decomposing. When rotating machinery breaks down, the frequency band distribution of vibration signal and spectrum energy all can change, and at this, first the vibration signal of rotating machinery carried out EEMD decomposition, obtain n IMF component; Secondly the energy T={E of each IMF component is calculated1,E2,��,En, namely
The IMF component energy feature obtained is normalized, i.e. T '={ E1/E,E2/E,��,En/ E}, whereinBy being analyzed above it can be seen that contain 6 temporal signatures, 2 frequency domain characters and n IMF component energy feature in the hybrid domain feature set of rotating machinery Incipient Fault Diagnosis, wherein n takes the integer between 5��10.
Adopting karyomerite function and overall situation Kernel multi-kernel support vector machine in described step 3, gaussian kernel function is typical karyomerite function, describes as follows:
KRBF(xi,xj)=exp (-| | xi-xj||2/��2)
Polynomial kernel function is typical overall situation kernel function, describes as follows:
Karyomerite function learning ability is strong, and generalization ability is weak; Overall situation kernel function generalization ability is strong, learning capacity is weak, traditional SVM adopts single core function, can effectively solve the problem that the classification identification problem of simple data, but the rotating machinery Incipient Fault Diagnosis problem for multiple complex data sources but has some limitations, in order to make SVM have better learning capacity and Generalization Capability, it is proposed to by karyomerite KRBFWith overall situation core KpolyCombine, construct multi-kernel support vector machine: Kmix(xi,xj), multi-kernel support vector machine has had the advantage of each individual core concurrently, has better classification identification performance.
Optimizing multi-kernel support vector machine detailed process based on heredity-annealing algorithm in step 4 is, is first initialize to control parameter, and arranging population scale is S, and independent evolutionary generation is N, and maximum evolutionary generation is M, and the crossover probability of population is p1, mutation probability is p2, initial annealing temperature is T0, overall situation evolutionary generation variable l1=0, randomly generate initial population;
According to the MSVM test result to training sample, evaluate the fitness of each individuality of current population, based on elite retention strategy, currently most individuality is deposited in memory apparatus; If optimum individual meets the condition of convergence in current population, then evolutionary process terminates, and returns global optimum individual, otherwise initializes Local Evolution algebraic variable l2=0;
For determining annealing temperature TiUnder population, in order to produce a new population, it is necessary to randomly select individual x from current populationi,xj, by crossover probability p1Carry out the real-valued operation that intersects, produce two new individual xi',x'j, calculate fitness f (xi'),f(x'j), select according to Metropolis acceptance criterion new individual, update memory apparatus simultaneously; Individuality after intersecting is pressed mutation probability p2Carry out real-valued mutation operation, calculate the fitness that variation is individual, select according to Metropolis acceptance criterion new individual, update memory apparatus simultaneously; l1=l1+ 1, l2=l2+ 1, if l2< N, then continue to produce new population;
If l1< M, is now not reaching to maximum iteration time, then revise the annealing temperature T of populationi, start iteration again, otherwise compare current population optimal solution and the value in memory apparatus, export optimal solution.
Optimize the characteristic vector of multi-kernel support vector machine as genetic Annealing after extracting the hybrid domain feature set normalized characterizing rotating machinery initial failure state; Utilize genetic Annealing to optimize multi-kernel support vector machine training sample is trained, obtain the optimal control parameter vector of multi-kernel support vector machine; Test sample is carried out identification diagnosis by the multi-kernel support vector machine Reconstruction initial failure intelligent diagnostics model finally utilizing optimum, exports final diagnostic result.
The present invention, based on the rolling bearing performance degeneration fuzzy granulation Forecasting Methodology of comentropy, adopts technique scheme can obtain following technique effect:
1) integrated application time domain, frequency domain, time-frequency domain signal processing method, never ipsilateral characterizes the time domain of running state of rotating machine, frequency domain, time-frequency domain statistical nature, it is possible to more accurately, more fully get the initial failure Weak characteristic of rotating machinery.
2) multi-kernel support vector machine is applied to the Incipient Fault Diagnosis of rotating machinery as the Incipient Fault Diagnosis model of a kind of superior performance, and the various coupling faults of rotating machinery can be made more identification conclusion accurately by it.
3) genetic annealing algorithms overcomes the premature problem of traditional genetic algorithm, solves the select permeability of multi-kernel support vector machine initial parameter, has higher versatility, robustness and degree of accuracy.
In sum, the present invention optimizes the rotating machinery Incipient Fault Diagnosis method of multi-kernel support vector machine based on genetic Annealing, analyse in depth the Major Difficulties in traditional rotating machinery Incipient Fault Diagnosis method, combine time domain, frequency domain, time-frequency domain signal processing method fault signature in early days extracts, multi-kernel support vector machine intelligent fault diagnosis in early days and genetic annealing algorithms advantage in multiparameter parallel optimization etc., define the integrated rotating machinery Incipient Fault Diagnosis pattern built by " the multi-kernel support vector machine fault diagnosis that fault features is extracted �� optimized ", maintain again the high accuracy in engineer applied simultaneously, high efficiency, the requirements such as replicability, a kind of new approach is provided for rotating machinery initial failure intelligent diagnostics.
Accompanying drawing explanation
Accompanying drawing 1 is the operation principle schematic diagram that in step 4 of the present invention, genetic Annealing optimizes multi-kernel support vector machine.
Accompanying drawing 2 is the operation principle schematic diagram that the present invention optimizes the rotating machinery Incipient Fault Diagnosis method of multi-kernel support vector machine based on genetic Annealing.
Detailed description of the invention
Below in conjunction with accompanying drawing 2, the embodiment of the rotating machinery initial failure intelligent diagnostics that the present invention optimizes multi-kernel support vector machine based on genetic Annealing is elaborated. The main purpose of the present embodiment is to be extracted by the Weak characteristic of initial failure in characteristic of rotating machines vibration signal by time domain, frequency domain, time-frequency domain signal processing method, multi-kernel support vector machine model is constructed based on Radial basis kernel function and Polynomial kernel function, genetic annealing algorithms is built based on genetic algorithm and simulated annealing, optimize multi-kernel support vector machine with genetic annealing algorithms, and then realize the Intelligence Diagnosis of rotating machinery initial failure. Embodiment comprises the following specific steps that:
Step one, gathers the vibrational state information in rotating machinery running by ICP acceleration transducer, subsequently this vibration signal is carried out noise reduction, filter preprocessing.
Described ICP sensor is the piezoelectric acceleration sensor of built-in miniature IC amplifier, traditional piezoelectric acceleration sensor and charge amplifier are integrated in one by it, directly can connect with record, display and collecting device, simplify test system, improve measuring accuracy and reliability. It is widely used in the fields such as nuclear explosion, Aero-Space, railway, bridge, building, car and boat, machinery, water conservancy, electric power, oil, chemical industry, environmental protection, earthquake.
Step 2, to at rotating machinery each vibration signal pretreated, calculate its temporal signatures index (subordinate list 1), frequency domain character index (subordinate list 2), time and frequency domain characteristics index respectively, wherein time and frequency domain characteristics index adopts integrated empirical mode decomposition technology (EEMD) to build by calculating the energy of mode function (IMF) component in each layer, gets up to build hybrid domain feature set by temporal signatures index (subordinate list 1), frequency domain character index (subordinate list 2), time and frequency domain characteristics index comprehensive;
Further, the vibration characteristic signals that rotating machinery produces in running is the Main Basis of reflection running state of rotating machine, and vibration signal can produce different characteristic parameter values to characterize the wave character of vibration signal at the statistical analysis of time domain and frequency domain. Temporal signatures value therein is the original foundation of rotary machinery fault diagnosis; Frequency domain character value is by carrying out spectrum transformation to original vibration signal, it is possible to extract more characteristic information. Therefore, by calculating the time domain of vibration signal, frequency domain character index value, it is possible to the appearance of quantitatively characterizing rotating machinery initial failure. Time domain, frequency domain character index are a lot, at this data analysis by experiment, have chosen and rotating machinery initial failure is characterized obvious 6 temporal signatures indexs and 2 frequency domain character indexs, as shown in Table 1 and Table 2, in practical engineering application, it is possible to add time domain, frequency domain character index flexibly according to practical situation.
Table 1 temporal signatures index
Table 2 frequency domain character index
Temporal signatures index in table 1 mainly reflects vibration signal amplitude energy size and seasonal effect in time series distribution situation; Frequency domain character index in table 2 mainly reflects the change of vibration signal main band position and the dispersion intensity of frequency spectrum. In order to obtain more fault characteristic information, reflect the initial failure state of rotating machinery comprehensively, integrated empirical mode decomposition (ensembleempiricalmodedecomposition, EEMD) is also adopted to extract the time and frequency domain characteristics index in original vibration signal.
Further, described integrated empirical mode decomposition (EEMD) is a kind of signal decomposition method based on data-driven, it is at empirical mode decomposition (empiricalmodedecomposition, EMD) on basis, by superposition white Gaussian noise, to reduce the modal overlap situation in signal decomposition process, the decomposition step of EEMD is as follows:
(1) in primary signal x (t), superposition average is 0, and standard deviation is the white Gaussian noise n of constantiT (), obtains signal x to be decomposedi(t):
xi(t)=x (t)+ni(t)
Wherein: xiT () is the signal after adding i & lt white Gaussian noise. The white noise n addediT () is too big or too little all can not effectively solve mode mixing problem compared with original vibration signal x (t), the too little change that will not cause EMD decomposition result of white noise ratio, then can make to decompose the harmonic components permeating some falsenesses in each IMF component obtained too greatly, under normal circumstances, white noise takes 0.1��0.4 times that standard deviation is primary signal standard deviation.
(2) to xiT () carries out empirical mode decomposition, obtain several IMF components | cij(t) and remainder ri(t), wherein cijT () decomposes, after representing that i & lt adds white Gaussian noise, the i-th IMF obtained.
(3) repeating M step (1), step (2), obtained IMF component carries out ensemble average, eliminate and repeatedly add the white noise impact on decomposition result, obtaining the EEMD IMF result finally decomposed is:
In formula: perform the EMD number of times M decomposed more big, be more conducive to the raising of Decomposition Accuracy under normal circumstances, but when M big to a certain extent time, the improvement of result is also extremely limited, at this through rotating machinery vibrating status signal is analysed in depth, take M=100.
(4) final result that vibration signal obtains through EEMD decomposition is:
Wherein: cjT () is original vibration signal x (t) limited the Intrinsic Mode component (Intrinsicmodefunction, IMF) obtained from high frequency to low frequency adaptive decomposition, r (t) is the residual components after decomposing.
When initial failure occurs in rotating machinery, the frequency band distribution of vibration signal and spectrum energy all can change, and at this, on the EEMD basis decomposed, first the vibration signal of rotating machinery carried out EEMD decomposition, obtain n IMF component; Secondly the energy T={E of each IMF component is calculated1,E2,��,En, namely
The IMF component energy feature obtained is normalized, i.e. T '={ E1/E,E2/E,��,En/ E}, whereinBy being analyzed above it can be seen that contain 6 temporal signatures, 2 frequency domain characters and n IMF component energy feature in the hybrid domain feature set of rotating machinery Incipient Fault Diagnosis, wherein n generally takes the integer between 5��10.
Step 3, based on gaussian kernel support vector machine and Polynomial SVM, by Gauss karyomerite KRBFWith multinomial overall situation core KpolyCombine, construct multi-kernel support vector machine:
Kmix(xi,xj)=�� KRBF(xi,xj)+(1-��)Kpoly(xi,xj)
In formula: xiAnd xjFor the characteristic vector of the input space, the kernel functional parameter of �� and d respectively gaussian kernel function and Polynomial kernel function, �� (0 < �� < 1) is weight factor.
Further, utilizing SVM to carry out in fault identification process, the selection of kernel function is most important, the discriminant function that different kernel functions is corresponding different, and then directly affects the identification precision of SVM. The kernel function of SVM mainly includes local kernels and global kernels, and gaussian kernel function is typical karyomerite function, describes as follows:
KRBF(xi,xj)=exp (-| | xi-xj||2/��2)
Polynomial kernel function is typical overall situation kernel function, describes as follows:
Karyomerite function learning ability is strong, and generalization ability is weak; Overall situation kernel function generalization ability is strong, and learning capacity is weak. Traditional SVM adopts single core function, it is possible to effectively solve the classification identification problem of simple data, but the rotating machinery Incipient Fault Diagnosis problem for multiple complex data sources but has some limitations. In order to make SVM have better learning capacity and Generalization Capability, it is proposed to by karyomerite KRBFWith overall situation core KpolyCombine, construct mixed nucleus SVM:
Kmix(xi,xj)=�� KRBF(xi,xj)+(1-��)Kpoly(xi,xj)
In formula: xiAnd xjFor the characteristic vector of the input space, the kernel functional parameter of �� and d respectively gaussian kernel function and Polynomial kernel function, �� (0 < �� < 1) is for regulating parameter. Mixed nucleus SVM has had the advantage of each individual core concurrently, has better classification identification performance.
Step 4, genetic Annealing algorithm optimization multi-kernel support vector machine, form chromosomal real number gene vectors X with the weight factor �� of multinuclear SVM, penalty parameter c, kernel functional parameter �� and d:
X=[��, c, ��, d]
Utilize multinuclear SVM to the classification accuracy f of hybrid domain feature set training sample to evaluate each chromosomal quality, adopt Metropolis acceptance criterion to determine new chromosome pij' whether enter population of future generation, in Metropolis acceptance criterion, adopt dynamic temperature-regulating strategy, if ti+1=(1-1/i)m��ti, wherein tiMoving back temperature parameter when evolving for i & lt, m �� Z is for controlling the decrease speed of temperature, in order to prevent the optimal solution of current population from losing in the next generation, adopts elite retention strategy, compares the optimal solution x in the fitness f of new explanation x and population*Fitness f*If f is better than f*, then x is replaced respectively with x and f*And f*��
Further, MSVM parameter algorithm flow process is optimized as shown in Figure 1 based on genetic annealing algorithms:
Specifically comprise the following steps that
(1) control parameter is initialized: arranging population scale is S, and independent evolutionary generation is N, and maximum evolutionary generation is M, and the crossover probability of population is p1, mutation probability is p2, initial annealing temperature is T0, overall situation evolutionary generation variable l1=0, randomly generate initial population;
(2) according to the MSVM test result to training sample, evaluate the fitness of each individuality of current population, based on elite retention strategy, currently most individuality is deposited in memory apparatus;
(3) if optimum individual meets the condition of convergence in current population, then evolutionary process terminates, and returns global optimum individual, otherwise initializes Local Evolution algebraic variable l2=0;
(4) for determining annealing temperature TiUnder population implement the following, to produce a new population:
I. from current population, randomly select individual xi,xj, by crossover probability p1Carry out the real-valued operation that intersects, produce two new individual xi',x'j, calculate fitness f (xi'),f(x'j), select according to Metropolis acceptance criterion new individual, update memory apparatus simultaneously;
Ii. the individuality after intersecting is pressed mutation probability p2Carry out real-valued mutation operation, calculate the fitness that variation is individual, select according to Metropolis acceptance criterion new individual, update memory apparatus simultaneously.
(5)l1=l1+ 1, l2=l2+ 1, if l2< N, then go to step (4), otherwise continue;
(6) if l1< M then revises the annealing temperature T of populationi, go to step (2), otherwise compare current population optimal solution and the value in memory apparatus, export optimal solution.
So far, the overall process of rotating machinery Incipient Fault Diagnosis is namely completed.
The present invention is not restricted to the described embodiments, and every utilizes principles of the invention and mode, through the technical scheme that conversion and replacement are formed, all in protection scope of the present invention.