CN101520651A - Analysis method for reliability of numerical control equipment based on hidden Markov chain - Google Patents
Analysis method for reliability of numerical control equipment based on hidden Markov chain Download PDFInfo
- Publication number
- CN101520651A CN101520651A CN200910060949A CN200910060949A CN101520651A CN 101520651 A CN101520651 A CN 101520651A CN 200910060949 A CN200910060949 A CN 200910060949A CN 200910060949 A CN200910060949 A CN 200910060949A CN 101520651 A CN101520651 A CN 101520651A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- numerical control
- control equipment
- reliability
- performance characteristic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Feedback Control In General (AREA)
Abstract
The invention provides an analysis method for reliability of numerical control equipment based on a hidden Markov chain. The method particularly comprises the following steps: 1, monitoring dynamic performance signals of the numerical control equipment, and abstracting the performance characteristic parameter values showing the change of the reliability of the numerical control equipment; 2, constructing a predictive model of the performance characteristic parameter values; 3, using the predictive model to predict the performance characteristic parameter values within the time needed for vectorization, and adopting the source coding method to vectorize the predictive values of the performance characteristic parameter values; 4, adopting the discrete hidden Markov chain model to identify the state transition probability of the numerical control equipment; and 5, utilizing the Chapman-Kolmogorov differential equation to establish a relational expression of the operational state and the state transition probability so as to deduce the probability of the numerical control equipment in different operating states, namely obtaining the reliability of the numerical control equipment. The method can accurately analyze, evaluate and predict the change of the reliability of the numerical control equipment before the numerical control equipment goes wrong, thereby avoiding the fault of the numerical control equipment and improving the operational reliability of the numerical control equipment.
Description
Technical field
The present invention belongs to the numerical control equipment technical field, specifically is a kind of reliability of numerical control equipment variation analysis method based on discrete Hidden Markov chain model.
Background technology
Reliability is a kind of inherent attribute of numerical control equipment, is meant that numerical control equipment carries out it in assignment of mission and condition thereof in the next stipulated time and require the ability of the various key functions of task.It has reflected the mass property in the numerical control equipment use, time factor and environmental factor when having considered the numerical control equipment operation.
At present, the general thinking of reliability of numerical control equipment analysis is by gathering and analyzing numerically controlled equipment composition parts and the design of machine system and the historical failure and the lifetime data of manufacturing information, field failure data and like product, start with from the Analysis on Fault Diagnosis of numerical control equipment, according to numerical control equipment work and functional characteristics, trouble-shooting pattern and reason, assessment reliability of numerical control equipment level, and the reliability innovative approach is proposed.Wherein, mostly the approach of field failure data acquisition is to follow the tracks of for a long time tens numerical control equipments, and fault data derives from the field notes and the maintenance report of numerical control equipment.This shows, the same model numerical control equipment that existing analysis method for reliability of numerical control equipment needs and works long hours, quantity is abundant, and a large amount of faults and lifetime data will be arranged; By the analysis of failure data, and based on certain statistics judgment criterion, the failure probability distributions of identification equipment, and then adopt mathematical statistics method to infer reliability of numerical control equipment, find out the weak link in the numerical control equipment use.
Yet this subject matter that faces based on the reliability of numerical control equipment research of fault data has: on the one hand, it needs a large amount of fault and lifetime data, to find the life-span distribution curve of numerical control equipment.May obtain the data of a large amount of numerical control equipment fault and lifetime data, particularly machine system in the actual engineering hardly.In addition, numerical control equipment breaks down through long-time running, and this moment, equipment may be in end of lifetime, caused the reliability result of gained to lack actual directive significance.On the other hand, it has only considered the normal and fault two states of numerical control equipment, does not consider performance state variation feature in the numerical control equipment use.For example, numerical control equipment breaks down and goes to toward having experienced a series of deterioration state, is embodied in the output performance parameter and progressively departs from its setting, until final generation disabler.This causes the user can't understand the real-time change conditions of reliability of numerical control equipment.
Therefore, press for a kind of new method that is used for the reliability of numerical control equipment field of design and development, it does not need to wait for that numerical control equipment breaks down through long-play, just can analyze, assess and predict the reliability of numerical control equipment change conditions exactly, thereby guides user is taked rational preventive measure in advance, prevents that fault from taking place.
Summary of the invention
The objective of the invention is to be at the deficiencies in the prior art, a kind of reliability of numerical control equipment variation analysis method based on discrete Hidden Markov chain model is provided, before breaking down, numerical control equipment analyzes, assess and predict the reliability of numerical control equipment change conditions exactly, find out weak link in the numerical control equipment use etc., avoid the generation of numerical control equipment fault, improve the operational reliability of numerical control equipment.
A kind of analysis method for reliability of numerical control equipment based on hidden Markov chain may further comprise the steps:
Dynamic property signal when (1) the monitoring numerical control equipment moves extracts the performance characteristic parameter value that more than one reflection reliability of numerical control equipment changes at more than one time point;
(2), make up the forecast model of performance characteristic parameter value according to the performance characteristic parameter value that extracts;
(3) use forecast model prediction T
0Each performance characteristic parameter value of individual time point carries out vector quantization to these predicted values, obtains the output observation sequence of performance characteristic parameter;
(4), adopt discrete Hidden Markov chain model to determine the running status transition probability matrix of numerical control equipment according to the output observation sequence;
(5) utilize Qie Puman-Andrei Kolmogorov differential equation to set up the relational expression of running status and status change probability matrix, infer that numerical control equipment is in the probability under the different running statuses.
Described step (2) adopts non-stationary autoregression integration running mean algorithm to set up the forecast model of j performance characteristic ginseng value
Wherein, t represents the predicted time point, and p is the autoregression order, and q is the running mean order, φ
zBe running mean coefficient, measuring error
Described step (3) is specially: the output observation sequence of performance characteristic parameter is expressed as
J=1 ..., m, m are the performance characteristic number of parameters, t represents the predicted time point, wherein,
i
0=2 ..., N-1, t
1∈ (1 ..., T
0),
Represent that j performance characteristic parameter is at t
1The predicted value of time point,
i
1=1 ..., N-1,
Represent that j performance characteristic parameter is in i
1The value of defining of zone of dispersion.
Described step (4) is specially: make numerical control equipment comprise 0~k+1 kind running status, the output observation sequence and the predefined original state transition probability matrix A that obtain according to step (3)
0, initial observation value probability matrix B
0With initial probability distribution vector π
0, utilize Bao Mu-Wei Erqi (Baum-Welch) algorithm iteration calculate numerical control equipment constantly status change model λ of t=(π), π represents the probability distribution vector for A, B,
The status change probability matrix
a
Cd{ c} is in state { probability of d}, 0≤a constantly at t+1 for numerical control equipment is in state constantly at t
Cd≤ 1, and
The observed reading probability matrix
b
CjFor numerical control equipment is in state constantly at t { observed reading appears during c}
Probability, 0≤b
Cj≤ 1, and
The beneficial effect that the present invention has is embodied in:
Numerical control equipment is similar to human body, and its function is decline gradually along with time lengthening, and the function that is embodied in numerical control equipment and performance be deterioration progressively, causes the no plan downtime rate increase of numerical control equipment generation, has reduced the numerical control equipment operational reliability.Patent of the present invention is monitored, is analyzed and predict by performance degradation and reliability change to numerical control equipment, can help enterprise to find the potential failure mode of numerical control equipment, formulate rational maintenance schedule and sudden inefficacy takes place to avoid numerical control equipment, thereby realize the reliability growth of numerical control equipment, this all has great importance for improving numerical control equipment utilization factor, minimizing numerical control equipment maintenance cost, prolongation numerical control equipment serviceable life and structure numerical control equipment health status monitoring system etc.
Description of drawings
Fig. 1 general principles figure of the present invention.
Fig. 2 vector quantization figure of the present invention.
Fig. 3 status change topology diagram of the present invention.
Fig. 4 reliability change of the present invention result of calculation figure.
Embodiment
Below in conjunction with drawings and Examples technical solutions according to the invention are described further.
With reference to Fig. 1, a kind of reliability of numerical control equipment change computing method based on discrete Hidden Markov chain model mainly comprise the performance characteristic parameter monitoring and the modules such as extraction, performance degradation model, the quantification of performance characteristic parameter vector, status change model and Calculation of Reliability of numerical control equipment.
1, performance characteristic parameter monitoring and extraction.
The performance characteristic parameter monitoring is by being installed on the numerical control equipment or on every side testing tools such as the sensor various dynamic property signals (as vibration signal, motor message, hydraulic pressure signal etc.) when coming the testing equipment operation, extracting characteristic parameter through signal Processing and analysis with extracting.The performance characteristic parameter of being extracted is relevant with the type and the disposal route thereof of signal, but it should follow following principle generally: (1) performance characteristic parameter should have clear physical meaning, (2) the performance characteristic parameter is easy to measure and is relatively stable, the key parameter that (3) performance characteristic parameter changes for the reflection reliability of numerical control equipment etc.For example,, and adopt time domain method to analyze the time history of vibration signal if gather the vibration signal of numerical control equipment with acceleration transducer, then the performance characteristic parameter that can extract have that peak-peak, root mean square refer to, peak factor and the nargin factor etc.
According to the testing tool of prior layout, put and measure respectively the l in the numerical control equipment use preset time, obtains the dynamic property signal.The signal of gathering is handled and is analyzed, extract the performance characteristic parameter that can reflect timely and accurately that the numerical control equipment performance state changes, promptly D={ (T, Y) | (t
i, y
Ij); I=1,2 ..., l; J=1,2 ..., m}, wherein t
iBe i Measuring Time point, y
IjFor at t
iJ the performance characteristic parameter that the moment measures.
2, set up the performance degradation model of numerical control equipment.
The discrete-time series { (t of check feature characteristic parameter
i, y
Ij); I=1,2 ..., l; J=1,2 ..., the stationarity of m}.If non-stationary then needs to adopt { the y of d jump point-score with non-stationary
IjTime series is converted into { w stably
IjTime series, promptly
Adopt the nARIMA algorithm to { w
IjTime series carries out modeling.Because j performance characteristic parameter w
IjEach value w in the not only preceding p of value step with it
(i-1) j, w
(i-2) j..., w
(i-p) jRelevant, and the random perturbation value ε in q preceding with it step
(i-1) j, ε
(i-2) j..., ε
(i-q) j, then obtain p rank autoregression q rank moving average model:
Wherein, θ
zBe autoregressive coefficient, φ
zBe the running mean coefficient,
Be measuring error,
Be variance.Introduce backward shift operator B, then formula (1) is converted to formula (2):
Wherein, the order p of model and q and unknown parameter θ
1..., θ
p, φ
1..., φ
qAnd σ
εDirectly have influence on model accuracy.Adopt minimal information criterion (AIC criterion) to determine model order, adopt nonlinear least square method to come the estimation model unknown parameter.In case determine model order and model parameter, then obtained non-stationary nARIMA (p, d, q) performance degradation model:
Wherein, φ
zBe running mean coefficient, measuring error
P is the autoregression order, and the representation model predicted value is relevant with preceding p rank performance parameter value, and q is the running mean order, and expression is carried out the correction on q rank to the model predication value error.
3, the performance characteristic parameter vector quantizes
In case set up the performance degradation model, then can predict the performance characteristic parameter situation of change of numerical control equipment in following a period of time.For example, { y
Tj, j=1 ..., m} represents one group of t moment performance characteristic parameter that do as one likes energy degradation model forecasting institute gets.These performance characteristic parameters have reflected the implicit state variation feature of numerical control equipment from the side.
Yet because numerical control equipment itself is very complicated, the performance characteristic parameter that is observed in its use can not be corresponding simply one by one with state.Therefore, the present invention uses discrete Hidden Markov chain model, performance state variation characteristic and reliability change conditions thereof when setting up numerical control equipment status change model and coming analyzing numerically controlled equipment operation.
Because the input quantity of status change model is necessary for limited discrete value, and the performance characteristic parameter that is obtained is a continuous real number value, adopts Lloyd (Laue moral) algorithm that the performance characteristic parameter is carried out vector quantization for this reason and handle to form discrete code set.
With reference to Fig. 2, the vector quantization process of performance characteristic parameter is: according to subregion vector (partition) and codebook vectors (codebook), vector quantization is divided into N zone in proper order with the characteristic parameter span, discrete observation value of each zone mapping.The required predicted time of the vector quantization T that counts
0Generally get three to six unit interval points, at this T
0In the individual predicted time point section first is to T
0Time point, j characteristic ginseng value is at each regional index value
Be defined as:
i
0=2,…,N-1,t
1∈(1,…,T
0);
Wherein, subregion vector (partition) is defined by the performance failure canonical matrix H of equipment, that is:
The numerical control equipment that forecasting institute is got is at 1~T
0Performance characteristic parameter set { y constantly
TjCarry out vector quantization, obtain the output observation sequence of performance characteristic parameter:
Wherein, the required predicted time point T of vector quantization
0Generally get the 3rd to the 6th unit interval point.
4, status change model
With reference to Fig. 3, set up the status change model of equipment.Running status transition probability is meant that the equipment running status is constant or turns to the probability of poorer running status.
In the status change topology diagram, state space S={0,1,2 ..., k, k+1} are numerical control equipment performance running status collection.Running status comprises normal condition, deterioration state and failure state.Circle is represented latent state among Fig. 3,0} represents normal condition, and 1,2 ..., k} represents k deterioration state, { k+1} represents failure state; Directed arc is called the transition arc, the transfer of expression state.Under the situation of not considering to keep in repair, numerical control equipment performance degradation process is irreversible; Each state can only be to the state transitions of its higher numbering in right side, and each state also can shift to self simultaneously.Therefore, make numerical control equipment be expressed as at t status change probability matrix A constantly:
Wherein, a
Cd{ c} is in state { probability of d}, 0≤a constantly at t+1 for numerical control equipment is in state constantly at t
Cd≤ 1, and
The state do as one likes of numerical control equipment can the characteristic parameter observed reading be come perception.According to numerical control equipment at T0 output observation sequence constantly
Make numerical control equipment be expressed as at t observed reading probability matrix B constantly:
Wherein, b
CjFor numerical control equipment is in state constantly at t { observed reading appears during c}
Probability, 0≤b
Cj≤ 1, and
Under the starting condition, numerical control equipment is in normal operating condition, and then the probability distribution of original state vector is π
0=(1,0 ..., 0), status change probability matrix A
0With observed reading probability matrix B
0Take random device to choose, thereby obtain the status change model λ of numerical control equipment under the starting condition
0=(A
0, B
0, π.)。
After running status collection, status change probability matrix A and observed reading probability matrix B definition, adopt discrete hidden Markov chain model solution status change probability matrix A and observed reading probability matrix B, be specially: with the performance characteristic parameter at 1~T
0Output observation sequence constantly
Be input to original state transition model λ
0=(A
0, B
0, π.) in, utilize Baum-Welch (Bao Mu-Wei Erqi) algorithm carries out iterative computation, makes model parameter progressively tend to more reasonably than the figure of merit, thus obtain numerical control equipment constantly status change model λ of t=(A, B, π).In order to verify the rationality of this status change model, adopt Forward-Backward (front and back to) probability of algorithm computation observation sequence O under given λ, i.e. P (O| λ).If P (O| λ) surpasses expectation value 0.8, then think status change model λ=(A, B are feasible π).
5, Calculation of Reliability
In case determine numerical control equipment t constantly status change model λ=(A, B π), then calculate the reliability change conditions of numerical control equipment.Detailed process is: make P
c(t)=P (q
t=c) the expression numerical control equipment is in the probability of c state constantly at t.According to Qie Puman-Andrei Kolmogorov differential equation, have:
P′(t)=P(t)·A (4)
Wherein, P (t)=(P
0(t), P
1(t) ..., P
k(t), P
(k+1)(t)) be state vector, P ' is the single order differential state vector of P (t) (t), and A is a t status change matrix constantly.Formula (4) is carried out Laplace (Laplce) conversion:
If numerical control equipment is in normal condition usually in starting condition, then have:
P(0)=(P
0(0),P
1(0),…,P
k(0),P
(k+1)(0))=(1,0,…,0,0)
Simultaneously, in t moment status change matrix A substitution formula (5), then calculate:
P(s)=(P
0(s),P
1(s),…,P
k(s),P
(k+1)(s))
Then, P (s) is carried out the Laplace inverse transformation, then obtain numerical control equipment is in different conditions constantly at t probability P (t)=(P
0(t), P
1(t) ..., P
k(t), P
(k+1)(t)), calculate numerical control equipment t reliability R (t)=1-P constantly
(k+1)Reliability index such as (t).And then, the reliability change conditions in the analyzing numerically controlled equipment use, the weak link of discovery numerical control equipment.
Embodiment:
This embodiment has provided the specific implementation process of the present invention in engineering practice, simultaneous verification should the invention validity.
With reference to Fig. 4, be research object with Ou Tai OTM-650 CNC milling machine, by the circular motion precision deterioration that simulating assembly causes because of backlass, infer reliability of numerical control equipment change rule.Test unit comprises Ou Tai OTM-650 CNC milling machine, Central China " century star " digital control system HMC-21M, Heidenhain KGM182 plane grating measurement mechanism and ordinary PC etc.Wherein, plane grating KGM182 is a kind of high-accuracy testing tool that is used for non-contact measurement equipment two dimensional motion contour accuracy.Making worktable by NC programming is that 50mm and speed of feed are that 2000mm/min does counterclockwise/circular motion clockwise with the radius in XOY plane, adopts the circular motion track of plane grating on-machine measurement equipment.Allow numerical control equipment move 480 hours, every interval sampling in about 20 hours 1 time obtains 24 groups of circular traces.According to the Circular Test of NC Machine Tools test stone that ISO230-4:1996 provides, 24 groups of circular traces of check and analysis to measure gained calculate corresponding circular motion accuracy characteristic amounts such as circle hysteresis, circle deviation and radius deviation, and are as shown in table 1.
Table 1 circular motion accuracy characteristic amount
Because every kind of characteristic quantity only has 24 measured values, in order to improve match and precision of prediction, adopt interpolation 96 data of sampling generation again according to measured value, 120 circular motion accuracy characteristic amounts are arranged altogether.For every kind of characteristic quantity, choose preceding 100 data as the model training collection, back 20 data are as the model checking collection.Find circular motion accuracy characteristic amount right and wrong stably through check, adopt first difference method that they are converted into time series { w stably for this reason
Ij.Lagging behind with circle is example, determines the order of performance degradation model according to AIC criterion, when (p, when q) being (5,2), the information function of model reaches minimum value 1.7939.And then, adopt nonlinear least square method to determine the model parameter of nARIMA (5,1,2), that is:
Thereby foundation circle hysteresis degradation model is:
y
ij=0.6847y
(i-1)j-0.8273y
(i-2)j+0.7267y
(i-3)j-0.0877y
(i-4)j-0.1747y
(i-5)j+…
+ε
ij-0.0498ε
(i-1)j-0.9502ε
(i-2)j;i=5,…,100,j=4
Utilize back 20 data to verify this model, the FPE of model (final prediction deviation) is 1.5459e-5.Reach a conclusion thus: this model satisfies requirement of actual application.
The process of hysteresis degradation model is justified in similar foundation, can set up the performance degradation model of following radius deviation and circle deviation respectively:
(1) radius deviation (F
Max) degradation model is:
y
ij=y
(i-1)j+0.9165y
(i-2)j-1.43y
(i-3)j+0.5338y
(i-4)j+0.0033y
(i-5)j+…
+ε
ij-0.0048ε
(i-1)j-0.9952ε
(i-2)j;i=5,…,100,j=1
(2) radius deviation (F
Min) degradation model is:
y
ij=1.267y
(i-1)j-0.1327y
(i-2)j-0.723y
(i-3)j+0.4271y
(i-4)j-0.1187y
(i-5)j+0.2803y
(i-6)j+…
+0.2942y
(i-7)j-0.1943y
(i-8)j+0.1601y
(i-9)j-0.1513y
(i-10)j+0.4735y
(i-11)j-0.166y
(i-12)j+…
-0.4402y
(i-13)j+ε
ij-0.7509ε
(i-1)j+0.2458ε
(i-2)j;i=13,…,100,j=2
(3) justifying deviation (G) degradation model is:
y
ij=0.5102y
(i-1)j+0.5102y
(i-2)j+ε
ij-0.0043ε
(i-1)j+0.207ε
(i-2)j;i=2,…,100,j=3
According to above-mentioned performance degradation model, the processing technology of combining with digital control equipment requires and uses rules, with the numerical control equipment failure state be divided into normally, slight degradation, seriously deterioration and inefficacy one of four states, and definite circular motion precision failure criteria, as shown in table 2.Then, adopt the Lloyd algorithm that the performance characteristic parameter of numerical control equipment is carried out vector quantization, obtain numerical control equipment t circular motion precision output observation sequence 0 constantly.
Table 2 circular motion precision failure criteria
Because numerical control equipment is in normal operating condition under starting condition, the probability distribution of original state vector π then
0=[1,0,0,0]
T, status change probability matrix A
0With observed reading probability matrix B
0Take random device to choose respectively, have:
Thereby set up the status change model λ under the original state
0=(A
0, B
0, π
0).With circular motion precision output observation sequence O substitution λ
0=(A
0, B
0, π
0) in, utilize the Baum-Welch algorithm, come the physical training condition transition model by iterative computation.It is 100 that total training iterations is set, and along with the increase of iterations, model output logarithm probable value progressively increases, enter the convergence state of saturation through model after 10 iteration, the output probability value almost no longer increases, and obtains t status change model λ=(A constantly, B, π), that is:
Simultaneously, adopt the Forward-Backward algorithm to calculate t circular motion constantly precision output observation sequence O in that λ=(π) the probable value P under (O| λ) is 80.84% for A, B.Reach a conclusion thus: the status change model λ of foundation=(A, B π) can satisfy application request.
At last, with λ=(π) π in and A are updated in the formula 5 for A, B, find the solution this formula, and obtain t state vector P (t)=(P of numerical control equipment constantly through the Laplace inverse transformation
0(t), P
1(t), P
2(t), P
3(t)), wherein:
P
0(t)=exp(0.9443t)
P
2(t)=13.77exp(0.9443t)-13.93exp(0.9397t)+0.1633exp(0.6142t)
P
3(t)=-8.28exp(0.9443t)+7.74exp(0.9397t)-0.0141exp(0.6142t)+0.5556exp(t)
Then numerical control equipment is R (t)=1-P in t fiduciary level constantly
3(t), with reference to Fig. 4.As shown in Figure 4, fiduciary level had been lower than 0.8 after numerical control equipment moved about 300 hours, and numerical control equipment is in serious deterioration state, and the probability that breaks down is very high.
Claims (4)
1, a kind of analysis method for reliability of numerical control equipment based on hidden Markov chain may further comprise the steps:
Dynamic property signal when (1) the monitoring numerical control equipment moves extracts the performance characteristic parameter value that more than one reflection reliability of numerical control equipment changes at more than one time point;
(2), make up the forecast model of performance characteristic parameter value according to the performance characteristic parameter value that extracts;
(3) use forecast model prediction T
0Each performance characteristic parameter value of individual time point carries out vector quantization to these predicted values, obtains the output observation sequence of performance characteristic parameter;
(4), adopt discrete Hidden Markov chain model to determine the running status transition probability matrix of numerical control equipment according to the output observation sequence;
(5) utilize Qie Puman-Andrei Kolmogorov differential equation to set up the relational expression of running status and status change probability matrix, infer that numerical control equipment is in the probability under the different running statuses.
2, a kind of analysis method for reliability of numerical control equipment based on hidden Markov chain according to claim 1 is characterized in that, described step (2) adopts non-stationary autoregression integration running mean algorithm to set up the forecast model of j performance characteristic ginseng value
Wherein, t represents the predicted time point, and p is the autoregression order, and q is the running mean order, and φ z is the running mean coefficient, measuring error
3, a kind of analysis method for reliability of numerical control equipment based on hidden Markov chain according to claim 1 is characterized in that, described step (3) is specially: the output observation sequence of performance characteristic parameter is expressed as
M is the performance characteristic number of parameters, and t represents the predicted time point, wherein,
4, a kind of analysis method for reliability of numerical control equipment according to claim 1 based on hidden Markov chain, it is characterized in that, described step (4) is specially: make numerical control equipment comprise 0~k+1 kind running status, the output observation sequence and the predefined original state transition probability matrix A that obtain according to step (3)
0, initial observation value probability matrix B
0With initial probability distribution vector π
0, utilize Bao Mu-Wei Erqi (Baum-Welch) algorithm iteration calculate numerical control equipment constantly status change model λ of t=(π), π represents the probability distribution vector for A, B,
The status change probability matrix
a
Cd{ c} is in state { probability of d}, 0≤a constantly at t+1 for numerical control equipment is in state constantly at t
Cd≤ 1, and
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910060949XA CN101520651B (en) | 2009-03-03 | 2009-03-03 | Analysis method for reliability of numerical control equipment based on hidden Markov chain |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910060949XA CN101520651B (en) | 2009-03-03 | 2009-03-03 | Analysis method for reliability of numerical control equipment based on hidden Markov chain |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101520651A true CN101520651A (en) | 2009-09-02 |
CN101520651B CN101520651B (en) | 2011-06-29 |
Family
ID=41081286
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910060949XA Expired - Fee Related CN101520651B (en) | 2009-03-03 | 2009-03-03 | Analysis method for reliability of numerical control equipment based on hidden Markov chain |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101520651B (en) |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799674A (en) * | 2010-03-28 | 2010-08-11 | 华中科技大学 | Method for analyzing service state of numerical control equipment |
CN102254184A (en) * | 2011-07-18 | 2011-11-23 | 华中科技大学 | Method for fusing multi-physical-domain feature information |
CN102411339A (en) * | 2011-11-30 | 2012-04-11 | 华中科技大学 | Method for evaluating performance reliability of numerical control equipment |
CN102496028A (en) * | 2011-11-14 | 2012-06-13 | 华中科技大学 | Breakdown maintenance and fault analysis method for complicated equipment |
CN102968569A (en) * | 2012-11-30 | 2013-03-13 | 西南大学 | Reliability assessment method for safety instrument system based on Markov model and D-S evidence theory |
CN103034170A (en) * | 2012-11-27 | 2013-04-10 | 华中科技大学 | Numerical control machine tool machining performance prediction method based on intervals |
CN104950709A (en) * | 2015-06-10 | 2015-09-30 | 周海昇 | Device efficiency analyzing system based on operating data characteristic recognition |
CN105607578A (en) * | 2014-11-20 | 2016-05-25 | 中国科学院沈阳计算技术研究所有限公司 | Method and device for numerical control system hardware security protection based on Markov determination method |
CN106095655A (en) * | 2016-05-31 | 2016-11-09 | 北京蓝海讯通科技股份有限公司 | A kind of method for detecting abnormality, application and monitoring device |
CN106484983A (en) * | 2016-09-29 | 2017-03-08 | 中国电力科学研究院 | A kind of appraisal procedure of power system relay protection device running status and device |
CN106803101A (en) * | 2016-12-30 | 2017-06-06 | 北京交通大学 | Odometer method for diagnosing faults based on HMM |
CN107562041A (en) * | 2017-09-22 | 2018-01-09 | 广东工业大学 | Goat fault early warning method, device, equipment and computer-readable recording medium |
CN108122007A (en) * | 2017-12-20 | 2018-06-05 | 西安航空学院 | A kind of foundation and its method of discrimination based on human factor forced outage model |
CN108764496A (en) * | 2018-05-16 | 2018-11-06 | 重庆大学 | Fault of numerical control machine tool repairs and reliability analysis system |
CN110757216A (en) * | 2019-11-04 | 2020-02-07 | 吉林大学 | Tool magazine manipulator reliability test method based on half Markov process |
CN111765075A (en) * | 2020-05-20 | 2020-10-13 | 天津市天锻压力机有限公司 | Hydraulic forging press pump source fault prediction method and system |
CN112685933A (en) * | 2020-12-24 | 2021-04-20 | 中国人民解放军海军工程大学 | Method for predicting residual service life of roller screw pair |
CN113051700A (en) * | 2019-12-27 | 2021-06-29 | 北京国双科技有限公司 | Equipment reliability monitoring method and device |
CN114693143A (en) * | 2022-04-07 | 2022-07-01 | 北京工业大学 | Health state evaluation method, system, equipment and medium of numerical control machine tool |
CN115577990A (en) * | 2022-12-08 | 2023-01-06 | 中国人民解放军国防科技大学 | Method for establishing equipment system task reliability model |
-
2009
- 2009-03-03 CN CN200910060949XA patent/CN101520651B/en not_active Expired - Fee Related
Cited By (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799674A (en) * | 2010-03-28 | 2010-08-11 | 华中科技大学 | Method for analyzing service state of numerical control equipment |
CN102254184A (en) * | 2011-07-18 | 2011-11-23 | 华中科技大学 | Method for fusing multi-physical-domain feature information |
CN102496028A (en) * | 2011-11-14 | 2012-06-13 | 华中科技大学 | Breakdown maintenance and fault analysis method for complicated equipment |
CN102496028B (en) * | 2011-11-14 | 2013-03-20 | 华中科技大学 | Breakdown maintenance and fault analysis method for complicated equipment |
CN102411339A (en) * | 2011-11-30 | 2012-04-11 | 华中科技大学 | Method for evaluating performance reliability of numerical control equipment |
CN103034170A (en) * | 2012-11-27 | 2013-04-10 | 华中科技大学 | Numerical control machine tool machining performance prediction method based on intervals |
CN103034170B (en) * | 2012-11-27 | 2014-10-29 | 华中科技大学 | Numerical control machine tool machining performance prediction method based on intervals |
CN102968569A (en) * | 2012-11-30 | 2013-03-13 | 西南大学 | Reliability assessment method for safety instrument system based on Markov model and D-S evidence theory |
CN102968569B (en) * | 2012-11-30 | 2016-04-06 | 西南大学 | Based on the safety instrumented systems reliability estimation method of Markov model and D-S evidence theory |
CN105607578A (en) * | 2014-11-20 | 2016-05-25 | 中国科学院沈阳计算技术研究所有限公司 | Method and device for numerical control system hardware security protection based on Markov determination method |
CN104950709A (en) * | 2015-06-10 | 2015-09-30 | 周海昇 | Device efficiency analyzing system based on operating data characteristic recognition |
CN106095655A (en) * | 2016-05-31 | 2016-11-09 | 北京蓝海讯通科技股份有限公司 | A kind of method for detecting abnormality, application and monitoring device |
CN106095655B (en) * | 2016-05-31 | 2018-06-12 | 北京蓝海讯通科技股份有限公司 | A kind of method for detecting abnormality, application and monitoring device |
CN106484983A (en) * | 2016-09-29 | 2017-03-08 | 中国电力科学研究院 | A kind of appraisal procedure of power system relay protection device running status and device |
CN106803101A (en) * | 2016-12-30 | 2017-06-06 | 北京交通大学 | Odometer method for diagnosing faults based on HMM |
CN106803101B (en) * | 2016-12-30 | 2019-11-22 | 北京交通大学 | Odometer method for diagnosing faults based on Hidden Markov Model |
CN107562041A (en) * | 2017-09-22 | 2018-01-09 | 广东工业大学 | Goat fault early warning method, device, equipment and computer-readable recording medium |
CN108122007A (en) * | 2017-12-20 | 2018-06-05 | 西安航空学院 | A kind of foundation and its method of discrimination based on human factor forced outage model |
CN108122007B (en) * | 2017-12-20 | 2022-04-01 | 西安航空学院 | Establishment and discrimination method of forced shutdown model based on human factors |
CN108764496A (en) * | 2018-05-16 | 2018-11-06 | 重庆大学 | Fault of numerical control machine tool repairs and reliability analysis system |
CN110757216A (en) * | 2019-11-04 | 2020-02-07 | 吉林大学 | Tool magazine manipulator reliability test method based on half Markov process |
CN113051700A (en) * | 2019-12-27 | 2021-06-29 | 北京国双科技有限公司 | Equipment reliability monitoring method and device |
CN113051700B (en) * | 2019-12-27 | 2023-12-29 | 北京国双科技有限公司 | Equipment reliability monitoring method and device |
CN111765075A (en) * | 2020-05-20 | 2020-10-13 | 天津市天锻压力机有限公司 | Hydraulic forging press pump source fault prediction method and system |
CN112685933A (en) * | 2020-12-24 | 2021-04-20 | 中国人民解放军海军工程大学 | Method for predicting residual service life of roller screw pair |
CN112685933B (en) * | 2020-12-24 | 2022-10-18 | 中国人民解放军海军工程大学 | Method for predicting residual service life of roller screw pair |
CN114693143A (en) * | 2022-04-07 | 2022-07-01 | 北京工业大学 | Health state evaluation method, system, equipment and medium of numerical control machine tool |
CN115577990A (en) * | 2022-12-08 | 2023-01-06 | 中国人民解放军国防科技大学 | Method for establishing equipment system task reliability model |
Also Published As
Publication number | Publication date |
---|---|
CN101520651B (en) | 2011-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101520651B (en) | Analysis method for reliability of numerical control equipment based on hidden Markov chain | |
CN110298455B (en) | Mechanical equipment fault intelligent early warning method based on multivariate estimation prediction | |
Luo et al. | Early fault detection of machine tools based on deep learning and dynamic identification | |
CN110303380B (en) | Method for predicting residual life of cutter of numerical control machine tool | |
CN109840362B (en) | Multi-objective optimization-based integrated just-in-time learning industrial process soft measurement modeling method | |
CN101799674B (en) | Method for analyzing service state of numerical control equipment | |
CN101797704B (en) | Method for thermal deformation error compensation of digital control gear hobbing machine | |
Lindemann et al. | Anomaly detection and prediction in discrete manufacturing based on cooperative LSTM networks | |
CN101870076B (en) | Method for predicting service life of guide pair of numerical control machine on basis of performance degradation model | |
CN104714537B (en) | A kind of failure prediction method based on the relative mutation analysis of joint and autoregression model | |
CN109583520B (en) | State evaluation method of cloud model and genetic algorithm optimization support vector machine | |
CN102937534B (en) | Method for predicting fault of electromechanical device based on combined prediction model | |
CN102062832B (en) | Power system online prewarning method based on microdisturbance signal low-frequency oscillation mode identification | |
CN104390657A (en) | Generator set operating parameter measuring sensor fault diagnosis method and system | |
CN110262450B (en) | Fault prediction method for cooperative analysis of multiple fault characteristics of steam turbine | |
CN113239132B (en) | Online out-of-tolerance identification method for voltage transformer | |
CN116380445B (en) | Equipment state diagnosis method and related device based on vibration waveform | |
CN105243393A (en) | Characteristic-based fault forecasting method for complex electromechanical system | |
CN105868164A (en) | Soft measurement modeling method based on monitored linear dynamic system model | |
CN111860839A (en) | Shore bridge fault monitoring method based on multi-signal fusion and Adam optimization algorithm | |
CN102539192A (en) | Failure prediction method based on ICA reconstruction | |
CN103245759A (en) | Product quality monitoring method based on autoregression total projection to latent structures (T-PLS) | |
Li et al. | Fault forecasting of a machining center tool magazine based on health assessment | |
CN103472732A (en) | Improved multivariate controller performance monitoring method based on Mahalanobis distance | |
CN114800036B (en) | Equipment health state assessment method |
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: 20110629 Termination date: 20120303 |