CN100468263C - Continuous miner remote real-time failure forecast and diagnosis method and device - Google Patents

Continuous miner remote real-time failure forecast and diagnosis method and device Download PDF

Info

Publication number
CN100468263C
CN100468263C CNB2007100127040A CN200710012704A CN100468263C CN 100468263 C CN100468263 C CN 100468263C CN B2007100127040 A CNB2007100127040 A CN B2007100127040A CN 200710012704 A CN200710012704 A CN 200710012704A CN 100468263 C CN100468263 C CN 100468263C
Authority
CN
China
Prior art keywords
centerdot
fault
detecting device
output
data
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.)
Expired - Fee Related
Application number
CNB2007100127040A
Other languages
Chinese (zh)
Other versions
CN101126929A (en
Inventor
张颖伟
刘建昌
高会军
孙得维
王滢
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northeastern University China
Original Assignee
Northeastern University China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northeastern University China filed Critical Northeastern University China
Priority to CNB2007100127040A priority Critical patent/CN100468263C/en
Publication of CN101126929A publication Critical patent/CN101126929A/en
Application granted granted Critical
Publication of CN100468263C publication Critical patent/CN100468263C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Landscapes

  • Testing And Monitoring For Control Systems (AREA)

Abstract

The utility model relates to a method and a device for long-range real-time failure prediction and diagnosis of continuous miner, wherein the method comprises following steps: collecting data, processing data, forming a detector, learning and training the immune neural network, predicting and diagnosing failure and displaying the result. While the device comprises a host computer, a data communication device, a lower computer and a local sensing transmitter, wherein the local sensing transmitter comprises an acceleration transducer, a vibration pickup, a flow sensor, a pressure sensor and a temperature sensor; the data communication device is connected with the host computer and the lower computer respectively and the lower computer is connected with the local sensing transmitter. The local sensing transmitter device sends the signals collected by the sensors to the lower computer, while the collected signals undergo A/D transformation in the lower computer at certain time and then are sent to the host computer via the data communication device in order to predict and diagnose the failure. The utility model has the advantages of high accuracy and small sample size, and can not only predict and diagnose the failure, but also absorb and utilize new technology.

Description

Continuous miner remote real-time failure forecast and diagnostic method and device
Technical field
The invention belongs to failure prediction and diagnostic techniques field, particularly a kind of continuous miner remote real-time failure forecast and diagnostic method and device.
Background technology
Machine is adopted by cutting unit, tractive unit, motor and four most of compositions of auxiliary equipment by the most widely used company at present.Cutting unit adopts gear drive, and tractive unit adopts hydrostatic transmission.Most of faults of Lian Caiji are mechanical part fault and hydraulic system fault.
Comprise that the diagnosis of (1) temperature monitoring, (2) utilization analytic ferrographic method adopt the machine fault and carry out monitoring, diagnosing, (3) and connect and adopt the machine fault diagnosis expert system and (4) adopt the machine fault diagnosis system based on the company of artificial neural network connecting and connect the machine of adopting method for diagnosing faults commonly used at present.
Above-mentioned (1), (2) traditional diagnosis method do not have initiatively the ability of " thinking study ", it can only make judgement according to the field data that is input to the diagnostic system the inside, if the data of input are incomplete or conflicting, it will can not get conclusion even get the wrong sow by the ear, and the traditional diagnosis method poor efficiency of relatively requiring great effort.Same system is adopted (1), (2) two kinds of conventional fault diagnosis method, may obtain the diagnosis of mutual contradiction sometimes, at this moment still can only lean on experienced maintenance personal manually to diagnose.Connect and to adopt the machine fault diagnosis expert system there are the following problems: 1. " bottleneck " of knowledge acquisition.The management of present most coal mine machinery equipments, maintenance work are very unsound, lack various raw readings and theoretical property and sum up; Fault history, diagnosis and the service experience etc. of equipment do not form system documentation.2. the fragility of system, promptly knowledge and experience is not comprehensive, and it is just powerless to run into the problem that did not solve.Have only superficiality and empirical knowledge, lack essence and rational knowledge, ignored the deep layer effects such as understanding to knowledge, in case the situation that heuristic rule is not considered occurs, system performance sharply descends.3. knowledge base management difficulty.When knowledge base was too huge, because expert system adopts serial mode, inference method was simple and control strategy is dumb, caused knowledge shot array problems such as search speed descends, rule contradicts each other.Because system does not possess association and self-learning capability, to the suitable trouble of renewal, modification and maintenance of knowledge base.4. inferential capability is poor, and intelligence is low.5. versatility is poor, and particular device needs specific expert system.(defective of Artificial Neural Networks-ANN) is: 1. Learning Algorithms is a nonlinear optimal problem to artificial neural network, has the local minimum problem inevitably.Although 2. many improvement error back propagations (BP) algorithm is arranged, yet its studying convergence speed is still very slow.3. network must be through more time-consuming off-line training, and will obtain comprehensively and abundant training sample is difficult.
Summary of the invention
At above-mentioned problems of the prior art, the present invention combines immune algorithm with neural network, and a kind of continuous miner remote real-time failure forecast and diagnostic method and device based on immune neural net is provided.
Apparatus of the present invention comprise that host computer, data communication apparatus, slave computer, on-the-spot sensing become and send device.On-the-spot sensing becomes send device to comprise acceleration, vibration, flow, pressure, temperature sensor.Slave computer comprises one or more Programmable Logic Controllers.What data communication apparatus adopted is the CAN bus.Wherein data communication apparatus links to each other with host computer, slave computer respectively, and slave computer send device to link to each other with on-the-spot sensing change.Adopt the sensing change of the on-the-spot erecting stage of machine in company and send device, on-the-spot sensing becomes the sensor that send in the device signal of gathering is delivered to slave computer, slave computer regularly is sent to host computer with acquired signal by the CAN bus, and host computer carries out analyzing and processing to data again, carries out failure prediction and diagnosis then.
The each several part function of apparatus of the present invention:
(1) device is sent in on-the-spot sensing change: comprise acceleration, vibration, flow, pressure, temperature sensor, be responsible for parameter acquisition and transmission;
(2) slave computer: the signal A/D of being responsible for gathering is changed, and signal is sent to host computer;
(3) data communication apparatus: the signal that slave computer is sent on the CAN bus is sent to host computer through the USB-CAN module converts;
(4) host computer: receive the slave computer data, carry out analyzing and processing, provide the failure condition of Lian Caiji by the inventive method;
The inventive method mainly contains following characteristics:
(1) the immune system Negative Selection Algorithm is improved, and sure selection algorithm and improved Negative Selection Algorithm have been carried out effective combination, improved the efficient that detecting device generates by eliminating redundant detecting device;
(2) simultaneously by the adjustment of matching threshold, the detecting device that can obtain promptly will take place He break down;
(3) the fault sample of the detecting device that obtains, make the both measurable fault that will take place of immune neural net, the diagnosable again fault that has taken place as neural network training.
Attainable failure prediction of the present invention and diagnosis comprise adopts machine reduction-gear housing failure prediction and diagnosis, main electrical fault prediction and diagnosis, Main Hydraulic Pump failure prediction and diagnosis and hydraulic system fault prediction and diagnosis to connecting.
Failure prediction of the present invention and diagnostic method may further comprise the steps:
The collection of step 1, data
Gather acceleration, vibration, flow, pressure, temperature data; And through the A/D conversion;
The processing of step 2, data
Because the nerve of BP neural network adopts the Sigmoid function, if crossing conference, the input quantity absolute value make neuron output saturated, make the weights adjustment enter the flat region of error curved surface then, will make this moment error not too responsive to the variation of weights.So, before neural network training, need make normalized to the input data, the scope that makes it between [0,1], the present invention adopt suc as formula:
x i ′ = x i - x min x max - x min ( HI - LO ) + LO - - - ( 1 )
Method is carried out normalized to data.X wherein iBe i data of measuring before the normalization, x MaxBe { x i, i=1,2 ..., the maximal value among the n; x MinBe { x i, i=1,2 ..., the minimum value among the n, x i' be the data after the normalization, HI=0.9, LO=0.1, here add that two of HI, LO mainly are in order to eliminate input complete 0 in data and complete 1, occur 0 or at 1 o'clock because can find to input or output from weights are adjusted formula, can make to adjust and stagnate, error no longer reduces.In addition in normalized, because of the data differences of considering various patterns little, all on same quantity rank, so in normalized, minimum and maximum value is the value in all sample datas.If the data differences in each fault mode is bigger, basic on different data ranks, so in normalized, minimum and maximum value is got the value of sample data in the various fault modes respectively, rather than the value in all sample datas.
Step 3, generation detecting device
The most original Negative Selection fault detection method has a lot of defectives: at first be that negative selection algorithm is not considered the mutual coupling between the detecting device, may cause some detecting devices to flock together, can not cover all nonego spaces, perhaps need very big detectors set.Therefore, when generating detecting device, the detecting device with the detecting device generation coupling that has generated is removed, generate the detecting device collection affirming that selection algorithm and improved Negative Selection Algorithm effectively combine to calculate again;
(1) definition oneself string S be contain limited character, length is the set of the character string of L.Oneself string is that the information data of L standard sign equipment state of the data pattern during real coding equipment operate as normal links together as the normal mode string, and chooses h and go here and there as the normal mode string R 0 ′ = [ X 1 , X 2 , · · · , X h ] , h = 1,2 , · · · . For certain fault mode i, i=1,2 ..., k (k is the kind of fault) links together the information data of the characterization device fault signature of the L of the fault mode of a real coding standard (but binding isotherm derivation) as the fault mode string
Figure C200710012704D00082
(being the detecting device individuality) chooses the h string as the fault mode string R i ′ = [ X hi + 1 , X hi + 2 , · · · , X hi + h ] . Obtain oneself's string S = ( R 0 ′ , R 1 ′ , · · · , R k ′ ) .
(2) set matching threshold z.The present invention obtains each detecting device collection, neural network training then by getting different matching threshold z respectively.By the graph of errors of neural metwork training process with select nominal situation at random and the fault diagnosis situation of the data detection immune neural net that various fault mode obtains, determine the size of more excellent matching threshold z then.Distance Matching principle that the present invention adopts European (Euclidean) characterizes the affinity size, as the formula (2), reaches necessary requirement (the more little then affinity of Euclidean distance is big more) if Euclidean distance d, shows affinity less than preset threshold z.
d = &Sigma; i = 1 L ( ab i - ag i ) 2 < z - - - ( 2 )
In the formula: d is an Euclidean distance; Ab tCoordinate for antibody; Ag iCoordinate for antigen; Z is matching threshold; L is the length of string.
(3) generate detecting device collection R iGenerate detecting device collection D at random, to each detecting device among the D, with the normal mode string among the S
Figure C200710012704D0008151547QIETU
Carry out affinity according to formula (2) and calculate (the detected set D that antibody in the formula and antigen are respectively oneself string S and generate at random),, then put into detecting device collection R if d less than preset threshold z, illustrates that affinity is enough big 0In, if can not satisfy condition, continue and first kind of fault mode string
Figure C200710012704D0008162158QIETU
Carry out affinity and calculate, if d, then puts into detecting device collection R less than preset threshold z 1In, if can not mate, will continue and following a kind of fault mode string
Figure C200710012704D0008162209QIETU
Carry out affinity and calculate, if d less than preset threshold z, puts into detecting device collection R so 2In, if can not mate, continue circulation, up to mating with last a kind of fault mode.Repeat this process till reaching needed detector number M.
(4) generation forecast detecting device collection.Getting matching threshold is z+a, calculates detecting device and first kind of fault mode string among the D
Figure C200710012704D0008162221QIETU
Affinity, if d less than preset threshold z+a and more than or equal to z, then puts into the detecting device collection
Figure C200710012704D00086
In, if can not mate, will continue and following a kind of fault mode string
Figure C200710012704D00087
Carry out affinity and calculate, if d less than preset threshold z+a and more than or equal to z, puts into the detecting device collection so
Figure C200710012704D00088
In, if can not mate, continue circulation, up to mating, obtain the detecting device collection with last a kind of fault mode R a = ( R 1 a , R 2 a , &CenterDot; &CenterDot; &CenterDot; , R k a ) . Repeat whole process till reaching needed detector number M.In like manner, get matching threshold and be respectively z+2a, z+3a, a ∈ R, the detecting device collection of each fault mode of different matching threshold correspondences is respectively R 2 a = ( R 1 2 a , R 2 2 a , &CenterDot; &CenterDot; &CenterDot; , R k 2 a ) With R 3 a = ( R 1 3 a , R 2 3 a , &CenterDot; &CenterDot; &CenterDot; , R k 3 a ) . The affinity of different matching threshold detecting devices is calculated by following European formula control
z &le; d = &Sigma; i = 1 L ( ab i - ag i ) 2 < z + a - - - ( 3 )
z + a &le; d = &Sigma; i = 1 L ( ab i - ag i ) 2 < z + 2 a - - - ( 4 )
z + 2 a &le; d = &Sigma; i = 1 L ( ab i - ag i ) 2 < z + 3 a - - - ( 5 )
In the formula: a is the regulated quantity of matching threshold; Z is matching threshold; D is an Euclidean distance.
The study of step 4, immune neural net and training
With the detecting device collection R that trains 0Normal sample as neural network training; With the detecting device collection R that trains i(i=1,2 ..., k), R a, R 2aAnd R 3aAs the fault sample of neural network training, determine that output 1 is fault, 0 is normal.Because some signal can change before fault took place, what answer is contingent fault sample in the future so matching threshold is the detecting device set pair of z+a, z+2a, z+3a, the fail result of the output of its neural network is indicating the fault that equipment promptly will take place, and the confidence level that will break down is represented in each matching threshold.
The neural network that the present invention adopts is three layers of BP network model.The BP network model is exactly the training algorithm that adopts error back propagation, i.e. the feedforward network be made up of non-linear unit of BP network.Can adopt a kind of learning algorithm of supervising formula, its basic thought is to make the error sum of squares of network output layer reach minimum by the weights of regulating networks and threshold value, just makes output valve as far as possible near expectation value.Neural network learning, diagnosis and prediction detailed process are as follows in the inventive method:
(1) determining of network structure and learning parameter: determine that according to input vector the input layer of three-layer neural network counts n, determines that according to the fault mode kind output layer node counts q, rule of thumb and by considering the study number of times and the discrimination of network behind the training study, comprehensively relatively the back determines that hidden layer node counts p and weight convergence factor ζ and error convergence factor delta.
(2) network status initialization: give each connection weight W Ij, V JtAnd threshold value θ j, γ tGive the random value between [1 ,+1].
(3) provide the neural network learning sample right: selected characteristic sample data vector A k = [ a 1 k , a 2 k , &CenterDot; &CenterDot; &CenterDot; , a n k ] And object vector Y k = [ y 1 k , y 2 k , &CenterDot; &CenterDot; &CenterDot; , y q k ] Offer neural network.
(4) begin study, each sample carried out following operation:
1. each neuronic output of computational grid hidden layer and output layer
Hidden layer output b j = f ( &Sigma; i = 1 n W ij &CenterDot; a i - &theta; j ) ( j = 1,2 , &CenterDot; &CenterDot; &CenterDot; , p ) - - - ( 7 )
Output layer output c t = f ( &Sigma; j = 1 p V jt &CenterDot; b j - &gamma; t ) ( t = 1,2 , &CenterDot; &CenterDot; &CenterDot; , q ) - - - ( 8 )
W in the formula IjThe connection weight of-input layer i unit to hidden layer j unit;
θ jThe threshold value of-hidden layer j unit;
V Jt-hidden layer j unit to output layer t unit connection weight;
γ t-output layer t cell threshode;
F-S type activation function f ( x ) = 1 1 + exp ( - x ) ;
2. calculation training error
The hidden layer training error e j k = [ &Sigma; t = 1 q d t k &CenterDot; V jt ] &CenterDot; b j &CenterDot; ( 1 - b j ) ( j = 1,2 , &CenterDot; &CenterDot; &CenterDot; , p ) - - - ( 9 )
The output layer training error d t k = ( y t k - c t ) &CenterDot; c t &CenterDot; ( 1 - c t ) ( t = 1,2 , &CenterDot; &CenterDot; &CenterDot; , q ) - - - ( 10 )
In the formula
Figure C200710012704D00104
-desired output;
c t-actual output;
b jThe output of-hidden layer j unit;
3. revise weights
Input layer is to the new connection weight between the hidden layer:
W ij ( N + 1 ) = W ij ( N ) + &beta; &CenterDot; e j k &CenterDot; a i k - - - ( 11 )
Hidden layer is to the new connection weight between the output layer:
V jt ( N + 1 ) = V jt ( N ) + &alpha; &CenterDot; d t k &CenterDot; b j - - - ( 12 )
In the formula
Figure C200710012704D00107
Respectively from formula (9) and (10);
N-study number of times;
α-learning coefficient, 0<α<1;
β-learning coefficient, 0<β<1;
(5) whether satisfy | b j(N+1)-b j(N) | 2<ζ, carried out for (6) step if satisfy, otherwise turned back to for (4) step.
(6) whether satisfy | b j(N)-Y k|<δ goes on foot if satisfy then carry out (7), otherwise turned back to for (2) step.
(7) weights and the threshold data that obtains deposited in the knowledge base of being set up.
Step 5, failure prediction and fault diagnosis
The neural network that trains is used for the failure prediction and the diagnosis of equipment, can judges according to output vector which kind of fault has taken place.For example: setting desired output and each modes relationships is: normal (0,0,0,0,0,0); Fault 1 (1,0,0,0,0,0); Fault 2 (0,1,0,0,0,0); Fault 3 (0,0,1,0,0,0); Prediction fault 1 (0,0,0,1,0,0) (confidence level 80%); Prediction fault 2 (0,0,0,0,1,0) (confidence level 80%); Prediction fault 3 (0,0,0,0,0,1) (confidence level 80%).Detect vector data to one group and input to the immune neural net system, after the prediction and diagnosis through immune neural net, which kind of pattern is actual output vector just can be defined as during near which kind of pattern desired output vector.When output vector is (0.2,0.1,0.2,0.1,0.1 , 0.8) time, with the prediction fault 3 (0,0,0,0,0, 1) more approaching, so determine fault for this reason.For the fault mode that increases newly, the process that generates detecting device by improved Negative Selection Algorithm obtains the detecting device collection, and neural network training gets final product then, uses very convenient.
Step 6, demonstration output result
The present invention is applied to Negative Selection Algorithm in the middle of the acquisition of neural network training fault sample according to immune characteristics, has improved the reliability of fault sample.Because some signal can change before fault took place, what choose fault sample correspondence that different matching threshold obtains is contingent fault sample in the future, Xun Lian neural network thus, the fail result of its output is indicating the fault that equipment promptly will take place, so the present invention can to or occurent equipment failure diagnose, can predict following imminent fault again.
The present invention---show as based on the remote real-time failure forecast of the Lian Caiji of immune neural net and the superiority of diagnostic method and device:
(1) is fit to the occasion that fault sample lacks,, comes neural network training according to the immune algorithm sample set that just can extensively be out of order as long as know a small amount of typical sample of fault; For the equipment of non-fault sample,, come neural network training according to the immune algorithm sample set that also can extensively be out of order if can reason out a spot of fault sample according to theory.
(2) for imminent fault, be difficult to extract feature, this method can be extracted the feature that is about to break down easily according to the adjustment of matching threshold.The neural network of training, both diagnosable fault, measurable again imminent fault.
Description of drawings
Fig. 1 is apparatus of the present invention structural representation;
Fig. 2 is example middle controller and controller CANOPEN wiring diagram;
Fig. 3 is example middle controller and host computer wiring diagram;
Fig. 4 chooses synoptic diagram for cut speed reduction unit measuring point;
Fig. 5 is for generating the process flow diagram of detecting device by improved Negative Selection Algorithm;
Fig. 6 adopts machine failure prediction and diagnosis immune neural net structure for connecting in the example;
Fig. 7 is an immune neural net training module process flow diagram in the example;
Fig. 8 reverse transmittance nerve network training process error curve diagram;
Fig. 9 immune neural net training process error curve diagram;
The verify error curve map of Figure 10 reverse transmittance nerve network;
The verify error curve map of Figure 11 immune neural net;
Figure 12 is a fault display interface in the example;
Figure 13 is a primary control program process flow diagram in the example;
Figure 14 is the output display interface of the immune neural net of cut reduction-gear housing;
Wherein: 1-measuring point, 1,2-measuring point, 2,3-measuring point, 3,4-measuring point, 4,5-motor, 6-Machine guarding, 7-Machine guarding, 8-motor, 9-left cutting mechanism, 10-reduction-gear housing, 11-connecting shaft, 12-right cutting mechanism.
Embodiment
Be used for connecting to adopt the failure prediction of machine cut reduction-gear housing and be diagnosed as example further specifying with the present invention below.
Apparatus of the present invention comprise that host computer, slave computer, data communication apparatus, on-the-spot sensing become and send device that wherein on-the-spot sensing becomes send device to comprise acceleration, vibration, flow, pressure, temperature sensor.As shown in Figure 1.
Slave computer adopts the third generation CAN controller of Finland Epec Oy company exploitation in this example---and EPECOY.E3002024 (being called for short EPEC2024), EPEC2024 is embedded CAN bus communication module carries CANopen and CAN2.0 bus interface.Only need the pairing stitch of CAN bus access controller is provided with its node number, just can realize the mutual communication of controller on the CAN bus, as shown in Figure 2.
Host computer links to each other by the embedded CAN bus communication module of USB interface and EPEC2024, CAN_GND, the CAN_L of the CAN bus communication module that EPEC2024 is embedded, CAN_H pin link to each other with XM4.3, XM4.6, the XM4.2 pin of controller respectively, as shown in Figure 3.Hop adopts twisted-pair feeder to realize host computer and the long-range carrier communication of slave computer, and communications protocol adopts the senior communications protocol of CANOPEN.
On-the-spot sensing becomes the signal that send device to adopt EPEC2024 controller pick-up transducers to be collected, and the transmission storage space of the data storage of being gathered to the EPEC2024 controller, by program data are sent on the CAN bus, transmission interface receives and sends to the information on the CAN bus and send into host computer.
Cutting mechanism is one of vitals of Lian Caiji, and it is made up of 2 motor, two cover Machine guardings, speed reduction unit, left and right sides cutting drum, cutter chain and a cutting arm etc.Two motors split in the cutting arm left and right sides with vertical machine y direction symmetry, send power to speed reduction unit by Machine guarding separately, and the shaft of speed reduction unit drives left and right cutting drum and the broken coal seam of cutter chain simultaneously.Because motor is axially vertical with the big direction of principal axis of intermediate calender rolls, so speed reduction unit adopted one-level spur gear, one-level cone gear (commutation) and three grades of decelerations of primary planet pinion, and reduction case is T-shaped.Cutting mechanism acts directly on the coal petrography, be cantilever position during cut, the stressed complexity of cylinder, the cut load change is bigger, cause mechanical disorder easily, in the course of the work, speed reduction unit bears high varying load, becomes high pulling torque, rate of breakdown is higher, most of faults of Lian Caiji appear at cutting unit speed reduction unit this part.In this example to five kinds of diagnosing malfunctions of distortion of its broken teeth, flank of tooth fatigue, deformation of tooth surface, tooth surface abrasion, axle.
Failure prediction and diagnostic method are carried out according to the following steps in this example:
The collection of step 1, data
When carrying out data acquisition, reasonably selecting measuring point is online prerequisite of carrying out failure prediction and diagnosis.According to the needs of analyzing and the design feature of equipment, according to the structure of reduction-gear housing and the selection principle of vibration signal measuring point, the measuring point arrangement as shown in Figure 4, measuring point 1 (3) and measuring point 2 (4) differ 180 °, 4 measuring points all are placed on the output shaft end cap.The frequency of 1~5f scope is divided into one section every 1f, measures extreme value in each segment, these 16 extreme values are joined together as the description of this data vector as the representative of this section.Adopt the HG-3518 data acquisition unit to carry out signal data acquisition, the instrument sampling number selects N=1024, maximum analysis frequency f=2000Hz.
The processing of step 2, data
Obtained 900 groups of data through experiment, before neural network training, the input data have been made normalized, the scope that makes it is between [0,1], and after employing was carried out normalized suc as formula (1) method to data, it was as shown in table 1 to have chosen 18 groups of data wherein.
The training sample data of table 1 neural network
Figure C200710012704D00131
Figure C200710012704D00141
Step 3, generation detecting device
(1) the self-string length of definition is the limited string assemble S of L=16.Choose three strings as the normal mode string R 0 &prime; = [ X 1 , X 2 , X 3 ] ; For certain fault mode i, (i=1,2 ..., 5), choose three strings as every kind of fault mode string R i &prime; = [ X 3 i + 1 , X 3 i + 2 , X 3 i + 3 ] .
(2) generate detecting device collection R iSet matching threshold z=0.500 earlier, generate detecting device collection D at random, for each detecting device among the D, with the normal mode among the S
Figure C200710012704D0014163205QIETU
Carry out affinity and calculate, if d less than preset threshold 0.500, then puts into detecting device collection R suc as formula (2) 0In, if can not satisfy condition, continue and broken teeth fault mode string
Figure C200710012704D0014152256QIETU
Carry out affinity and calculate, if d, then puts into detecting device collection R less than preset threshold 0.500 1In, if can not mate, will continue and following a kind of flank of tooth wear-out failure pattern string
Figure C200710012704D0014163218QIETU
Carry out affinity and calculate, if d less than preset threshold 0.500, puts into detecting device collection R so 2In, if can not mate, continue circulation, up to mating, obtain detecting device collection R with last a kind of bending shaft fault mode 5Repeat this process till every kind of pattern reaches needed detector number 400.So just generated the normal mode of operation R of cut reduction-gear housing 0, and the detecting device collection of each fault mode: broken teeth R 1, the tired R of the flank of tooth 2, deformation of tooth surface R 3, tooth surface abrasion R 4With bending shaft R 5
Ideal matching threshold value z=0.420 has been determined in the training of the immune neural net by three kinds of matching thresholds, then generation forecast detecting device collection.The setting matching threshold is z+a=0.420+0.020=0.440, promptly gets a=0.020.(the detecting device collection is R for detecting device among the calculating D and the data pattern among the S 0) carry out affinity calculating suc as formula (2), if d less than preset threshold 0.500, then deletes,, continue and broken teeth fault mode string if can not satisfy condition
Figure C200710012704D0014163436QIETU
Carry out affinity and calculate, if d less than preset threshold z+a=0.440 and more than or equal to z=0.420, then puts into the detecting device collection
Figure C200710012704D0014163443QIETU
In, if can not mate, will continue and following a kind of fault mode string
Figure C200710012704D0014163452QIETU
Carry out affinity and calculate, if d less than preset threshold z+a=0.440 and more than or equal to z=0.420, puts into the detecting device collection so
Figure C200710012704D0014163501QIETU
In, if can not mate, continue circulation, up to mating, obtain the detecting device collection with last a kind of fault mode
Figure C200710012704D00144
Repeat whole process till reaching needed detector number M=400.In like manner, get matching threshold and be respectively z+2a, z+3a (a=0.020), the detecting device collection of each fault mode of different matching threshold correspondences is respectively
Figure C200710012704D00145
With
Figure C200710012704D00146
The affinity of different matching threshold detecting devices is calculated by formula (3), (4), (5) and is calculated.
Be the detecting device R that z=0.420 and a=0.020 generate with matching threshold like this 0, (R 1, R 2..., R 5),
Figure C200710012704D00151
Sample training neural network as neural network.The immune neural net of training both can carry out failure prediction, can carry out fault diagnosis again.
The study of step 4, immune neural net and training
According to input vector determine input layer count n=16, according to fault mode kind output layer node count q=21, the hidden layer node number is a key of finding the solution network problem, it directly has influence on the recognition capability of network model to test sample book, rule of thumb and by the node that this immune neural net behind the training study is chosen hidden layer counts p=15.The final network structure of determining is 16-15-21, as shown in Figure 6.Selected weight convergence factor ζ=0.001 and error convergence factor delta=0.1, training objective 0.001, maximum training batches 2000 times.
The corresponding relation of table 2 cut reduction-gear housing typical fault and failure symptom (different threshold value)
Figure C200710012704D00152
Figure C200710012704D0016152619QIETU
Determine the network model and the training data of neural network this moment, then determined the encoding relation of the ideal output of pattern to be discriminated and neural network.Provided the corresponding relation of cut reduction-gear housing typical fault to be discriminated with failure symptom (different threshold value) in table 2, it is that equipment operation is normal that output is at 0 o'clock.
Comparison example
(1) contrast of immune neural net and BP neural metwork training graph of errors
Train the BP neural network to obtain the graph of errors of BP neural metwork training process shown in Figure 8 with the data that 900 groups of experiments obtain, as can be seen from the figure after having passed through 1398 training, just reached the target error of setting.
During matching threshold z=0.420 above, the detecting device collection R of acquisition 0(R 1, R 2..., R 5) train immune neural net as the training sample of immune neural net.Obtain the graph of errors of immune neural net training process shown in Figure 9, as can be seen from the figure after having passed through 1058 training, just reached the target error limit of setting.
By in the training error curve of Fig. 8 and Fig. 9 as can be seen this immune neural net only need 1058 training just can reach the training objective error of former BP neural network through training just to reach for 1398 times.So generally speaking immune neural net can be saved time of 1/4th than former BP neural network, so immune neural net error convergence speed is very fast.
(2) fault diagnosis examples comparative
With 60 groups of check data randomly drawing the BP neural network that trains is tested, in table 3~table 8, provided the actual output diagnostic result of the BP neural network after the check.With the error (∞-norm as shown in Equation 14) of ideal output with the ∞-norm mark BP neural network of actual output result's difference, calculate detection error as shown in table 9, for the ease of observing, this paper is according to the error result of the table 9 verify error curve of BP neural network that drawn accordingly, as shown in figure 10.
‖x—y‖ =max(|x 1—y 1|,|x 2—y 2|,…,|x 21—y 21|) (14)
Wherein<x 1, x 2..., x 21Be the real output value of BP neural network,<y 1, y 2..., y 21It is the idea output of BP neural network.
Table 3 nominal situation is the diagnostic result of corresponding BP neural network down
Nominal situation Broken teeth Flank of tooth fatigue Deformation of tooth surface Tooth surface abrasion Bending shaft
0.894 0.027 0.135 -0.391 0.033 0.073
0.652 0.087 -0.217 0.136 0.249 -0.276
0.617 -0.156 0.107 0.312 -0.061 0.327
0.583 0.269 0.266 0.290 0.224 0.055
0.877 0.156 0.301 -0.213 0.167 0.195
0.686 -0.030 0.177 0.292 -0.324 -0.122
0.597 0.195 0.314 0.136 0.164 0.192
0.751 0.354 -0.281 0.272 0.246 0.255
0.746 0.183 0.467 0.330 0.360 0.143
0.906 0.158 0.144 -0.281 0.278 0.185
The output of the BP neural network of table 4 broken teeth correspondence
Nominal situation Broken teeth Flank of tooth fatigue Deformation of tooth surface Tooth surface abrasion Bending shaft
-0.035 0.797 0.437 -0.212 0.215 -0.186
0.306 0.852 -0.213 0.192 0.374 0.360
-0.074 0.601 0.314 0.266 -0.136 0.316
0.133 0.751 0.285 0.348 0.171 0.500
0.115 0.615 -0.126 0.319 0.415 0.370
0.093 0.582 -0.475 0.309 0.414 0.345
0.079 0.911 0.193 0.242 0.243 0.381
0.145 0.687 0.329 0.271 0.181 0.265
-0.044 0.817 0.387 0.166 0.442 0.434
0.124 0.922 0.187 0.239 0.203 0.199
The output of the BP neural network of the tired correspondence of table 5 flank of tooth
Nominal situation Broken teeth Flank of tooth fatigue Deformation of tooth surface Tooth surface abrasion Bending shaft
0.103 0.384 0.684 0.243 0.275 0.168
0.082 0.262 0.710 0.304 0.304 -0.081
-0.085 0.319 0.744 0.263 0.231 0.125
0.176 0.341 0.572 0.397 0.460 0.285
0.184 0.218 0.843 0.179 -0.121 0.260
0.128 0.320 0.915 -0.148 0.254 0.211
0.038 -0.166 0.738 0.391 0.183 0.181
0.242 0.208 0.605 0.322 0.411 -0.028
0.145 0.210 0.757 0.231 0.196 0.222
-0.024 0.179 0.597 0.342 0.109 0.247
The output of the BP neural network of table 6 deformation of tooth surface correspondence
Nominal situation Broken teeth Flank of tooth fatigue Deformation of tooth surface Tooth surface abrasion Bending shaft
-0.101 0.180 0.262 0.546 0.475 0.168
0.067 -0.213 0.161 0.794 0.221 0.256
-0.230 0.314 0.341 0.651 0.381 -0.279
0.129 0.236 0.189 0.814 0.214 0.186
0.061 0.434 0.143 0.731 0.115 0.333
0.021 0.112 0.268 0.911 0.196 -0.162
0.142 0.255 0.312 0.607 0.360 0.169
0.282 0.296 0.243 0.813 0.247 0.219
0.110 0.220 -0.190 0.741 0.295 0.126
-0.012 0.190 0.178 0.605 0.301 0.243
The output of the BP neural network of table 7 tooth surface abrasion correspondence
Nominal situation Broken teeth Flank of tooth fatigue Deformation of tooth surface Tooth surface abrasion Bending shaft
-0.001 0.265 0.376 0.423 0.827 0.336
0.014 0.362 -0.188 -0.173 0.912 -0.188
0.115 0.381 0.266 0.380 0.682 0.396
0.170 0.307 0.439 0.321 0.643 0.338
0.093 -0.175 0.215 0.291 0.702 0.211
-0.128 0.371 0.272 0.263 0.592 0.228
0.092 0.270 -0.320 0.216 0.709 -0.124
0.167 0.269 0.148 -0.312 0.617 0.301
0.317 0.247 0.335 0.456 0.728 0.163
0.113 0.323 0.241 0.385 0.723 0.262
The output of the BP neural network of table 8 bending shaft correspondence
Nominal situation Broken teeth Flank of tooth fatigue Deformation of tooth surface Tooth surface abrasion Bending shaft
-0.103 -0.104 -0.213 0.140 0.257 0.901
0.092 0.170 0.326 0.296 -0.223 0.719
-0.102 0.152 0.381 0.117 0.283 0.882
0.107 0.144 0.315 0.378 -0.176 0.724
0.044 0.219 0.175 0.145 0.283 0.584
0.053 -0.355 0.142 0.402 0.262 0.827
-0.030 0.298 0.292 0.246 -0.181 0.728
0.124 0.192 -0.125 -0.299 0.302 0.620
0.092 -0.126 0.295 0.354 0.317 0.597
0.081 0.308 0.204 0.271 0.114 0.742
The detection deviation of table 9 BP neural network diagnostic result
Figure C200710012704D00191
Maximum 0.475 in 60 verify errors as can be seen from Table 9, minimum then is 0.186.Can find out after this network is through training it is to satisfy the fault diagnosis of even adopting machine cut reduction-gear housing substantially to require intuitively from the verify error curve.But the diagnostic result deviation is also bigger, and is most of between 0.350~0.460, forms erroneous judgement easily according to the output of BP neural network.
Then the immune neural net that trains is tested, in table 10~table 15, provided the actual output diagnostic result of the immune neural net after the check with 60 groups of check data randomly drawing.With ideal output and the actual error of the ∞ of difference-norm labelled immune neural network as a result of exporting, draw detection error result shown in table 16 at last, for the ease of observing, this paper is according to drawn the accordingly verify error curve of immune neural net of the error result of table 16, as shown in figure 11.
The output of the immune neural net of table 10 nominal situation correspondence
Nominal situation Broken teeth Flank of tooth fatigue Deformation of tooth surface Tooth surface abrasion Bending shaft
0.853 0.184 -0.014 0.209 0.122 0.205
0.981 -0.088 0.189 0.135 -0.139 0.191
0.862 0.123 -0.109 -0.077 0.051 -0.166
0.986 -0.103 0.034 0.105 0.152 0.073
0.862 0.086 0.103 0.132 0.038 0.152
0.845 0.171 0.097 0.067 0.102 0.201
0.953 0.050 0.035 0.086 0.115 0.143
0.852 0.049 0.100 -0.212 0.215 0.169
0.794 0.207 0.072 0.122 0.121 -0.148
0.902 0.142 0.104 0.201 0.077 0.056
The output of the immune neural net of table 11 broken teeth correspondence
Nominal situation Broken teeth Flank of tooth fatigue Deformation of tooth surface Tooth surface abrasion Bending shaft
0.082 0.998 -0.131 0.039 -0.123 0.064
-0.093 0.906 0.045 -0.184 0.202 -0.141
0.176 0.876 0.091 0.204 -0.114 0.013
-0.161 0.891 0.158 -0.155 0.226 0.171
0.050 0.798 -0.117 0.063 0.201 -0.149
0.107 0.848 0.130 0.180 0.150 0.099
-0.124 0.923 -0.204 0.191 -0.138 0.072
0.185 0.884 -0.082 -0.097 0.209 -0.176
-0.219 0.785 0.158 -0.076 -0.169 0.225
0.110 0.958 0.183 0.034 0.098 0.169
The output of the immune neural net of the tired correspondence of table 12 flank of tooth
Nominal situation Broken teeth Flank of tooth fatigue Deformation of tooth surface Tooth surface abrasion Bending shaft
0.089 0.213 0.984 0.143 0.073 0.175
0.123 -0.072 0.893 -0.117 0.148 -0.185
0.159 0.113 0.839 0.026 -0.092 0.099
-0.105 0.156 0.877 -0.109 0.212 0.224
0.046 -0.125 0.789 0.018 0.032 0.236
0.267 0.087 0.885 0.056 0.139 0.229
0.228 0.172 0.966 -0.118 -0.137 -0.194
-0.134 0.243 0.887 0.239 0.144 0.156
0.121 0.244 0.977 0.132 -0.024 0.117
0.165 0.153 0.874 0.218 0.123 0.237
The output of the immune neural net of table 13 deformation of tooth surface correspondence
Nominal situation Broken teeth Flank of tooth fatigue Deformation of tooth surface Tooth surface abrasion Bending shaft
0.170 0.072 0.104 0.994 0.141 -0.138
-0.222 -0.157 0.129 0.869 -0.211 0.078
0.081 0.312 -0.055 0.792 0.127 -0.165
-0.089 -0.225 0.193 0.896 -0.160 0.088
0.157 0.242 0.128 0.867 0.215 0.151
0.166 0.185 -0.186 0.982 0.174 0.155
-0.076 -0.220 0.234 0.885 0.184 -0.217
0.216 0.133 -0.283 0.876 -0.088 0.133
-0.098 -0.233 0.219 0.915 0.152 0.038
0.154 0.131 0.205 0.853 0.210 0.141
The output of the immune neural net of table 14 tooth surface abrasion correspondence
Nominal situation Broken teeth Flank of tooth fatigue Deformation of tooth surface Tooth surface abrasion decreases Bending shaft
0.249 0.116 0.229 0.172 0.925 0.287
-0.072 -0.098 -0.209 -0.169 0.878 0.201
0.040 0.087 0.130 0.240 0.876 0.193
-0.092 -0.015 -0.181 0.051 0.955 -0.056
0.106 0.214 0.133 0.247 0.797 0.073
-0.090 0.125 0.217 0.212 0.851 0.221
0.137 0.215 0.163 0.213 0.918 0.145
0.128 -0.244 0.230 0.176 0.877 -0.058
0.223 0.173 -0.228 0.218 0.941 0.290
-0.085 0.065 0.158 -0.174 0.901 0.159
The output of the immune neural net of table 15 bending shaft correspondence
Nominal situation Broken teeth Flank of tooth fatigue Deformation of tooth surface Tooth surface abrasion Bending shaft
0.082 0.209 0.285 0.072 0.115 0.937
-0.150 -0.178 -0.172 0.158 0.220 0.873
0.058 0.137 0.179 -0.223 -0.119 0.905
0.063 0.254 0.231 0.192 0.212 0.886
0.185 0.042 0.134 0.241 0.094 0.793
0.187 0.108 0.162 0.196 0.092 0.963
0.285 0.126 -0.253 0.257 0.109 0.875
0.223 -0.119 0.062 -0.241 -0.165 0.929
-0.139 0.097 0.166 0.184 0.211 0.903
0.173 0.158 0.210 0.264 0.062 0.859
The detection deviation of table 16 immune neural net diagnostic result
Figure C200710012704D00221
Maximum in 60 verify errors as can be seen from Table 16 also have only 0.290, and minimum then has only 0.131.Be to satisfy to connect the fault diagnosis of adopting machine cut reduction-gear housing after the immune neural net process is trained as can be seen by verify error.Can compare error-checking size very intuitively by Figure 10 and Figure 11.And the error convergence speed of comparing immune neural net with the BP neural network is very fast, and output bias is less, and accuracy is higher.
It is faster than BP neural network speed of convergence to draw the immune neural net speed of convergence according to the comparison of training error curve; The diagnostic result that can draw immune neural net according to the comparison of verify error curve is more accurate than the diagnostic result of BP neural network.So the method for diagnosing faults ratio based on immune neural net that invention proposes is more superior based on the BP Neural Network Method for Fault Diagnosis.
Step 5, failure prediction and fault diagnosis
In order better to verify accuracy and the practicality that immune neural net carries out failure prediction and diagnosis, the application practice by the eastern coal mining of god scene has fully shown accuracy of the present invention and practicality.
In the eastern coal mining production run of god, on Lian Caiji, used the inventive method and carried out failure prediction and diagnosis.The confidence level that can judge the fault that tooth surface abrasion takes place from the output of immune neural net on April 16 is 85%, because for suitable stage reason, the cut reduction-gear housing is not safeguarded this moment, still allows it proceed running.And, because construction environment is abominable, noise is bigger, to compare with other reduction-gear housing, this moment, imperceptible reduction-gear housing from the sense organ had fairly obvious variation.Find that reduction-gear housing sends big abnormal sound but spent two; body temperature is above 70 ℃; and can judge that the confidence level of the fault that flank of tooth fatigue takes place is 70% in the output of apparatus of the present invention at this moment, the confidence level that the fault of tooth surface abrasion takes place is 100%, should in time carry out shutdown maintenance.Find through tearing the reduction-gear housing maintenance open: occur a spot of spot corrosion, peeling phenomenon on the flank of tooth of one-level spur gear, and the primary planet pinion flank of tooth has a large amount of cuts, because this accident treatment is in time appropriate, avoided complete machine to be in and stagnated even state of paralysis, guaranteed working continuously of production.
The present invention---based on continuous miner remote real-time failure forecast and the diagnostic method and the device of immune neural net, can compare remote live judgement fast and accurately to connecting the operational situation of adopting machine cut reduction-gear housing, avoid the blindness of maintenance, fully guaranteed the normal operation of Lian Caiji.
Step 6, output is display interface as a result, as shown in figure 14
The major failure of cut reduction-gear housing has been listed in output as a result in the display interface: broken teeth, flank of tooth fatigue, deformation of tooth surface, tooth surface abrasion and bending shaft, correspondingly listed failure cause, output result, confidence level (for different threshold values) and historical record etc.It was fault greater than 0.65 o'clock that native system is set immune neural net output, whether the output result only shows fault, there is the result of fault output to adopt the pink colour font simultaneously, be convenient to observe the output result like this, click the historical record button, can see detailed historical record (date of breaking down, time, output result etc.).
For the software development part, select suitable developing instrument can reduce the difficulty of system development, accelerate development progress.The design of this monitoring software adopts Visual C++6.0 as developing instrument based on WindowsXP operating system.The design philosophy of this monitoring software is that the failure prediction and the diagnostic method that adopt the present invention to propose are monitored underground equipment, mainly comprises the generation of data acquisition, detecting device, the study of neural network and the realization of parts such as training, failure prediction and fault diagnosis, fault log management.
The main control process flow of host computer realization failure prediction and diagnosis is carried out according to the following steps, as shown in figure 13:
Step 1, beginning;
The setting of step 2, baud rate and controller equiment number;
Step 3, initialization USB-CAN modular converter;
Step 4, data acquisition;
The processing of step 5, data;
The generation of step 6, detecting device;
Step 7, neural network learning and training;
Whether step 8, failure judgement prediction and diagnosis have fault; If, execution in step 9; If not, return step 4;
The type of step 9, diagnostic process failure judgement;
Step 10, fault prompting are reported to the police, and send failure message to the down-hole, return step 4.
The element that Lian Caiji often breaks down also has main motor, hydraulic power unit.The fault of motor mainly shows as bolt looseness, bearing fault and resonance that mechanical vibration cause, make rotor eccentricity sweep that thorax, back end cover for motor eccentric wear, water cooling are bad to cause motor temperature rising etc.Generally speaking, when being out of order, hydraulic system all at first expects hydraulic pump, because pump is the power source of hydraulic system, hydraulic pump is hydraulic system " heart ", it converts mechanical energy to the hydraulic pressure energy, provide the fluid of enough flows and pressure to system, finish the needs of exercises to satisfy machinery, the action of all elements all will have relation with pump.The fault of in advance considering hydraulic pump is necessary, because can regularly trace fault rootstock by logic like this, helps fixing a breakdown.The hydraulic pump most common failure has: (1) not oil extraction; (2) vibration and noise; (3) pressure anomaly; (4) underfed or shakiness; (5) abnormal heating etc.
The present invention---can also predict and diagnose the fault of main motor, hydraulic power unit based on the method and apparatus of the continuous miner remote real-time failure forecast of immune neural net and diagnosis, can both well carry out remote live judgement fast and accurately to corresponding component, increased the Homes Using TV of Lian Caiji, greatly improve the reliability of Lian Caiji, guaranteed the efficient high yield of its workplace.

Claims (3)

1, a kind of continuous miner remote real-time failure forecast and diagnostic method, this method depends on host computer, slave computer, data communication apparatus and the change of on-the-spot sensing and send device, embedded the functional software of realizing failure prediction and diagnosis in the host computer, it is characterized in that this method may further comprise the steps:
Step 1, image data
Gather acceleration, vibration, flow, pressure, temperature data, and through the A/D conversion;
Step 2, deal with data
The input data are made normalized, and the scope that makes it adopts following formula between [0,1]:
x i &prime; = x i - x min x max - x min ( HI - LO ) + LO - - - ( 1 )
X wherein iBe i data of measuring before the normalization, x MaxBe { x i, i=1,2 ..., the maximal value among the n; x MinBe { x i, i=1,2 ..., the minimum value among the n, Be the data after the normalization, HI=0.9, LO=0.1;
Step 3, generation detecting device carry out according to the following steps
(1) definition oneself string S is for containing limited character, and length is the set of the character string of L;
Oneself string is that the information data of L standard sign equipment state of the data pattern during real coding equipment operate as normal links together as the normal mode string, and chooses h and go here and there as the normal mode string R 0 &prime; = [ X 1 , X 2 , &CenterDot; &CenterDot; &CenterDot; , X h ] , H=1,2, For certain fault mode i, i=1,2 ..., k links together the information data of the characterization device fault signature of the L of the fault mode of a real coding standard as the fault mode string
Figure C200710012704C00024
, choose the h string as the fault mode string R i &prime; = [ X hi + 1 , X hi + 2 , &CenterDot; &CenterDot; &CenterDot; , X hi + h ] , Obtain oneself's string S = ( R 0 &prime; , R 1 &prime; , &CenterDot; &CenterDot; &CenterDot; , R k &prime; ) ;
(2) set matching threshold z
Obtain each detecting device collection, neural network training then by getting different matching threshold z respectively; By the graph of errors of neural metwork training process with select nominal situation at random and the fault diagnosis situation of the data detection immune neural net that various fault mode obtains, determine more excellent matching threshold z then; Adopt the Euclidean distance matching principle to characterize the affinity size, as the formula (2), if Euclidean distance d, shows that affinity reaches necessary requirement less than preset threshold z;
d = &Sigma; i = 1 L ( ab i - ag i ) 2 < z - - - ( 2 )
In the formula: d is an Euclidean distance; Ab iCoordinate for antibody; Ag iCoordinate for antigen; Z is matching threshold; L is the length of string;
(3) generate detecting device collection R i
Generate detecting device collection D at random, to each detecting device among the D, with the normal mode string among the S
Figure C200710012704C00028
Carry out affinity according to formula (2) and calculate, wherein antibody in the formula and antigen are respectively the detected set D that the oneself goes here and there S and generates at random; If d less than preset threshold z, illustrates that affinity is enough big, then put into detecting device collection R 0In, if can not satisfy condition, continue and first kind of fault mode string
Figure C200710012704C00031
Carry out affinity and calculate, if d, then puts into detecting device collection R less than preset threshold z 1In, if can not mate, will continue and following a kind of fault mode string Carry out affinity and calculate, if d less than preset threshold z, puts into detecting device collection R so 2In, if can not mate, continue circulation, up to mating with last a kind of fault mode; Repeat this process till reaching needed detector number M;
(4) generation forecast detecting device collection
Getting matching threshold is z+a, calculates detecting device and first kind of fault mode string among the D
Figure C200710012704C00033
Affinity, if d less than preset threshold z+a and more than or equal to z, then puts into the detecting device collection
Figure C200710012704C00034
In, if can not mate, will continue and following a kind of fault mode string
Figure C200710012704C00035
Carry out affinity and calculate, if d less than preset threshold z+a and more than or equal to z, puts into the detecting device collection so
Figure C200710012704C00036
In, if can not mate, continue circulation, up to mating, obtain the detecting device collection with last a kind of fault mode R a = ( R 1 a , R 2 a , &CenterDot; &CenterDot; &CenterDot; , R k a ) ; Repeat whole process till reaching needed detector number M; In like manner, get matching threshold and be respectively z+2a, z+3a, a ∈ R, the detecting device collection of each fault mode of different matching threshold correspondences is respectively R 2 a = ( R 1 2 a , R 2 2 a , &CenterDot; &CenterDot; &CenterDot; , R k 2 a ) With R 3 a = ( R 1 3 a , R 2 3 a , &CenterDot; &CenterDot; &CenterDot; , R k 3 a ) ; The affinity of different matching threshold detecting devices is calculated by following European formula control:
z &le; d = &Sigma; i = 1 L ( ab i - ag i ) 2 < z + a - - - ( 3 )
z + a &le; d = &Sigma; i = 1 L ( ab i - ag i ) 2 < z + 2 a - - - ( 4 )
z + 2 a &le; d = &Sigma; i = 1 L ( ab i - ag i ) 2 < z + 3 a - - - ( 5 )
In the formula: a is the regulated quantity of matching threshold; Z is matching threshold; D is an Euclidean distance;
The study of step 4, immune neural net and training;
With the detecting device collection R that trains 0Normal sample as neural network training; With the detecting device collection R that trains 1, R a, R 2aAnd R 3aAs the fault sample of neural network training, i=1 wherein, 2 ..., k determines that output 1 is fault, 0 is normal; Because some signal can change before fault took place, what answer is contingent fault sample in the future so matching threshold is the detecting device set pair of z+a, z+2a, z+3a, the fail result of the output of its neural network is indicating the fault that equipment promptly will take place, and the confidence level that will break down is represented in each matching threshold;
Step 5, failure prediction and fault diagnosis;
The neural network that trains is used for the failure prediction and the diagnosis of equipment, judges according to output vector which kind of fault has taken place;
Step 6, demonstration output result.
2, continuous miner remote real-time failure forecast according to claim 1 and diagnostic method is characterized in that the study of the immune neural net in the described step 4 and training detailed process are as follows:
(1) determining of network structure and learning parameter: determine that according to input vector the input layer of three-layer neural network counts n, determines that according to the fault mode kind output layer node counts q, rule of thumb and by considering the study number of times and the discrimination of network behind the training study, comprehensively relatively the back determines that hidden layer node counts p and weight convergence factor ζ and error convergence factor delta;
(2) network status initialization: give each connection weight W Ij, V JtAnd threshold value θ j, γ tGive the random value between [1 ,+1];
(3) provide the neural network learning sample right: selected characteristic sample data vector A k = [ a 1 k , a 2 k , &CenterDot; &CenterDot; &CenterDot; , a n k ] And object vector Y k = [ y 1 k , y 2 k , &CenterDot; &CenterDot; &CenterDot; , y q k ] Offer neural network;
(4) begin study, each sample carried out following operation:
1. each neuronic output of computational grid hidden layer and output layer
Hidden layer output b j = f ( &Sigma; i = 1 n W ij &CenterDot; a i - &theta; j ) - - - ( 7 )
Output layer output c t = f ( &Sigma; j = 1 p V jt &CenterDot; b j - &gamma; t ) - - - ( 8 )
J=1 in the formula, 2 ..., p;
t=1,2,…,q;
W IjThe connection weight of-input layer i unit to hidden layer j unit;
θ jThe threshold value of-hidden layer j unit;
V Jt-hidden layer j unit to output layer t unit connection weight;
γ t-output layer t cell threshode;
F-S type activation function f ( x ) = 1 1 + exp ( - x ) ;
2. calculation training error
The hidden layer training error e j k = [ &Sigma; i = 1 q d t k &CenterDot; V jt ] &CenterDot; b j &CenterDot; ( 1 - b j ) - - - ( 9 )
The output layer training error d t k = ( y t k - c t ) &CenterDot; c t &CenterDot; ( 1 - c t ) - - - ( 10 )
In the formula
Figure C200710012704C00048
-desired output;
c t-actual output;
b jThe output of-hidden layer j unit;
3. revise weights
Input layer is to the new connection weight between the hidden layer:
W ij ( N + 1 ) = W ij ( N ) + &beta; &CenterDot; e j k &CenterDot; a i k - - - ( 11 )
Hidden layer is to the new connection weight between the output layer:
V jt ( N + 1 ) = V jt ( N ) + &alpha; &CenterDot; d t k &CenterDot; b j - - - ( 12 )
In the formula
Figure C200710012704C00053
Figure C200710012704C00054
Respectively from formula (9) and (10);
N-study number of times;
α-learning coefficient, 0<α<1;
β-learning coefficient, 0<β<1;
(5) whether satisfy | b j(N+1)-b j(N) |<ζ, carried out for (6) step if satisfy, otherwise turned back to for (4) step;
(6) whether satisfy | b j(N)-Y k|<δ goes on foot if satisfy then carry out (7), otherwise turned back to for (2) step;
(7) weights and the threshold data that obtains deposited in the knowledge base of being set up.
3, realize the device of described continuous miner remote real-time failure forecast of claim 1 and diagnostic method, it is characterized in that this device comprises host computer, data communication apparatus, slave computer, on-the-spot sensing becomes send device, on-the-spot sensing becomes send device to comprise acceleration, vibration, flow, pressure, temperature sensor, slave computer comprises one or more Programmable Logic Controllers, what data communication apparatus adopted is the CAN bus, wherein data communication apparatus respectively with host computer, slave computer links to each other, slave computer send device to link to each other with on-the-spot sensing change, adopt the sensing change of the on-the-spot erecting stage of machine in company and send device, on-the-spot sensing becomes the sensor that send in the device signal of gathering is delivered to slave computer, slave computer regularly is sent to host computer with acquired signal by the CAN bus, host computer carries out analyzing and processing to data again, carries out failure prediction and diagnosis then.
CNB2007100127040A 2007-09-05 2007-09-05 Continuous miner remote real-time failure forecast and diagnosis method and device Expired - Fee Related CN100468263C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2007100127040A CN100468263C (en) 2007-09-05 2007-09-05 Continuous miner remote real-time failure forecast and diagnosis method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2007100127040A CN100468263C (en) 2007-09-05 2007-09-05 Continuous miner remote real-time failure forecast and diagnosis method and device

Publications (2)

Publication Number Publication Date
CN101126929A CN101126929A (en) 2008-02-20
CN100468263C true CN100468263C (en) 2009-03-11

Family

ID=39094985

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2007100127040A Expired - Fee Related CN100468263C (en) 2007-09-05 2007-09-05 Continuous miner remote real-time failure forecast and diagnosis method and device

Country Status (1)

Country Link
CN (1) CN100468263C (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI574137B (en) * 2015-05-21 2017-03-11 Mitsubishi Electric Corp Remote monitoring devices, remote monitoring and maintenance systems, remote monitoring methods and remote monitoring program products
CN111650919A (en) * 2020-05-14 2020-09-11 中铁第四勘察设计院集团有限公司 Multi-dimensional monitoring escalator fault prediction and health management method and system

Families Citing this family (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8095261B2 (en) * 2009-03-05 2012-01-10 GM Global Technology Operations LLC Aggregated information fusion for enhanced diagnostics, prognostics and maintenance practices of vehicles
CN101697078B (en) * 2009-10-22 2011-11-16 北京航空航天大学 On-line fault diagnosis device and on-line fault diagnosis method of simple-type environmental control system
CN101699359B (en) * 2009-10-28 2012-01-04 上海理工大学 Method for visualizing fault state monitoring
US8676432B2 (en) 2010-01-13 2014-03-18 GM Global Technology Operations LLC Fault prediction framework using temporal data mining
CN102243497B (en) * 2011-07-25 2013-10-02 江苏吉美思物联网产业股份有限公司 Networking technology-based remote intelligent analysis service system used for engineering machinery
CN102998996A (en) * 2012-12-07 2013-03-27 陕西千山航空电子有限责任公司 Airborne real-time fault diagnosis method
CN103245373B (en) * 2013-04-09 2017-02-08 哈尔滨工程大学 Method for diagnosing faults of sensor of underwater robot
CN103745005A (en) * 2014-01-24 2014-04-23 清华大学 Dynamic artificial immunity fault diagnosis method based on vaccine transplantation
CN103914735B (en) * 2014-04-17 2017-03-29 北京泰乐德信息技术有限公司 A kind of fault recognition method and system based on Neural Network Self-learning
CN104020709B (en) * 2014-04-24 2017-08-25 燕山大学 A kind of dynamic control demo system of oil refining process displaying
CN104035433B (en) * 2014-05-30 2017-01-11 北京必可测科技股份有限公司 Remote precise diagnosis and maintenance system and method of industrial devices
JP6148316B2 (en) * 2015-07-31 2017-06-14 ファナック株式会社 Machine learning method and machine learning device for learning failure conditions, and failure prediction device and failure prediction system provided with the machine learning device
DE102016008987B4 (en) 2015-07-31 2021-09-16 Fanuc Corporation Machine learning method and machine learning apparatus for learning failure conditions, and failure prediction apparatus and failure prediction system including the machine learning apparatus
CN106598791B (en) * 2016-09-12 2020-08-21 湖南微软创新中心有限公司 Industrial equipment fault preventive identification method based on machine learning
CN106502181A (en) * 2016-12-20 2017-03-15 合肥中辰轻工机械有限公司 A kind of intelligence control system for bottle washer safety detection
CN107122538B (en) * 2017-04-24 2020-11-17 中国电力科学研究院 Temperature rise optimization method and system for extra-high voltage direct current converter station joint terminal
CN107421763B (en) * 2017-08-01 2019-11-15 上海庆科信息技术有限公司 A kind of equipment fault detection method and device
CN107643712A (en) * 2017-09-25 2018-01-30 中国航空工业集团公司西安飞机设计研究所 A kind of Multi-sensor intelligent monitoring device
US11237550B2 (en) * 2018-03-28 2022-02-01 Honeywell International Inc. Ultrasonic flow meter prognostics with near real-time condition based uncertainty analysis
CN110610226A (en) * 2018-06-14 2019-12-24 北京德知航创科技有限责任公司 Generator fault prediction method and device
CN109211390B (en) * 2018-08-24 2020-12-15 国网山东省电力公司电力科学研究院 Safety test and evaluation method for vibration and strength of power transmission tower
TWI673723B (en) * 2018-09-20 2019-10-01 帆宣系統科技股份有限公司 Intelligent pre-diagnosis and health management system and method
CN111061148A (en) * 2018-10-17 2020-04-24 帆宣系统科技股份有限公司 Intelligent pre-diagnosis and health management system and method
TWI667660B (en) * 2018-11-20 2019-08-01 帆宣系統科技股份有限公司 Intelligent pre-diagnosis and health management system modeling method and computer program product
CN111539546A (en) * 2019-02-01 2020-08-14 帆宣系统科技股份有限公司 Modeling method of intelligent pre-diagnosis and health management system and computer program product thereof
CN110017989B (en) * 2019-05-17 2021-06-08 山东科技大学 Method for diagnosing bearing fault of wind turbine
CN110457776B (en) * 2019-07-22 2022-06-14 电子科技大学 Test strategy rapid generation method based on fault decision network
CN110477846A (en) * 2019-08-23 2019-11-22 西安交通大学医学院第二附属医院 A kind of digestive endoscopy off-note real-time mark method and system
CN111105074A (en) * 2019-11-22 2020-05-05 江南机电设计研究所 Fault prediction method based on improved deep belief learning
CN112417954B (en) * 2020-10-13 2022-12-16 同济大学 Bearing fault mode diagnosis method and system for small sample data set

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
基于神经网络和专家系统的故障诊断. 徐常胜,周兆英,刘思行,张耀清.控制与决策,第10卷第4期. 1995 *
基于神经网络模式识别的板形模糊控制器. 刘建昌,王柱.东北大学学报(自然科学版),第26卷第8期. 2005 *
基于联合神经网络的冗余传感器故障诊断和信号重构. 张颖伟,王福利.信息与控制,第32卷第2期. 2003 *
采煤机故障诊断的神经网络方法. 张寒松,贾瑞清,王廷军,杨振鹏.矿山机械,第2004年卷第11期. 2004 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI574137B (en) * 2015-05-21 2017-03-11 Mitsubishi Electric Corp Remote monitoring devices, remote monitoring and maintenance systems, remote monitoring methods and remote monitoring program products
CN111650919A (en) * 2020-05-14 2020-09-11 中铁第四勘察设计院集团有限公司 Multi-dimensional monitoring escalator fault prediction and health management method and system

Also Published As

Publication number Publication date
CN101126929A (en) 2008-02-20

Similar Documents

Publication Publication Date Title
CN100468263C (en) Continuous miner remote real-time failure forecast and diagnosis method and device
CN100572810C (en) A kind of sucker-rod pumping machine oil pumping control method and system thereof
CN108051211B (en) A kind of wind generator set main shaft holds temperature pre-warning diagnostic method
EP3221579B1 (en) Wind turbine condition monitoring method and system
Garga et al. Hybrid reasoning for prognostic learning in CBM systems
CN1854706B (en) Methods and systems for diagnosing machinery
CN100470416C (en) Power plant thermal equipment intelligent state diagnosing and analyzing system
CN111043023B (en) Fracturing pump on-line monitoring and fault diagnosis system
CN106761668A (en) Oil well failure intelligent analysis decision system and method
CN106409120A (en) Machine learning method and machine learning device, and fault prediction device and fault prediction system
CN108256556A (en) Wind-driven generator group wheel box method for diagnosing faults based on depth belief network
CN104748962B (en) Planetary gear box intelligent diagnosis method based on stacking automatic encoding machine
CN105987822A (en) Method and system for predicting equipment failure
CN102741676A (en) Operation machine monitoring diagnosis device
CN103180793A (en) Anomaly diagnostic device and industrial machine
CN103814335A (en) Shovel, shovel management device, and shovel management method
NZ590408A (en) An online diagnostic method and system for a geothermal generation facility
CN116228186A (en) Ship cabin intelligent operation and maintenance system based on human engineering
CN116186946B (en) Hydraulic system fault diagnosis method and system based on diagnosis model
Ghulanavar et al. Diagnosis of faulty gears by modified AlexNet and improved grasshopper optimization algorithm (IGOA)
US20150275659A1 (en) Telemetry Diagnostics
Peng et al. Knowledge-based prognostics and health management of a pumping system under the linguistic decision-making context
CN205607641U (en) Online failure diagnosis device of simulation coal -winning machine rocking arm transmission system
CN110058576A (en) Equipment fault prognostics and health management method based on big data
CN116628544A (en) Fault diagnosis method based on hydraulic winch state on-line monitoring

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20090311

Termination date: 20091009