Gear Health Monitoring and Prognosis
|
|
- Randolph Williams
- 5 years ago
- Views:
Transcription
1 Gear Health Monitoring and Prognosis Matej Gas perin, Pavle Bos koski, -Dani Juiric ic Department of Systems and Control Joz ef Stefan Institute Ljubljana, Slovenia Abstract Many industrial set-ups nowadays rely on different types of gears to transmit rotational torque. A gear usually operates in a closed environment, immersed into lubricant (oil). Therefore, direct monitoring of its condition is not possible. We have constructed a pilot test bed in order to develop and test different methods for gear health monitoring and predicting the time of failure. In this paper we present a novel approach to the estimation of gear failure, based on measured vibration signals. From these signals features are extracted first and then the evolution over time is monitored. Expectation-Maximization algorithm is used to estimate the parameters of a discrete linear state-space model that is used to model the features time series. Time to reach safety alarm threshold is determined by making the prediction using the estimated linear model. Keywords - Condition monitoring; Gear systems; Signal processing; Expectation-Maximization I. I NTRODUCTION An open problem in industry is that - according to the estimates - over 95% of installed drives belong to the older generation with no embedded diagnostic functionality. This means that most of them are poorly monitored or even not monitored at all. As many of such machines will still be operating for some time, an upgrade with low-cost intelligent condition monitoring module would significantly improve surveillance, reduce maintenance costs, decrease failure rates and increase equipment availability. Indeed, timely localization of the root cause would allow for more efficient proactive maintenance as well as less costly and better planed fault accommodation. The goal of this work is to construct a reliable machine fault prognostic system to forecast damage propagation in rotary machinery. For this purpose, an experimental test bed has been constructed, which consists of an motor-generator pair, connected with single stage gearbox (c.f. Figure 1). Gear health can be assessed by processing vibration signals by means of Hilbert transform (HT) [1], [2]. Components from Hilbert spectrum can be used as fault indicators (or features). Feature us a real number calculated every ten minutes from five seconds long signal acquisition session. The sequence of features can be viewed as a time series i.e. a realization os the stochastic process. The objective of this work is to extend feature extraction with the prediction of the feature time series in order to estimate the time of occurrence of the safety alarm. There are several different techniques available for timeseries prediction [3], among them neural networks and neurofuzzy systems. In this paper we assume that the time series can be viewed as an output of a second order linear, discrete time, stochastic dynamic system. Because of the progress of faults in gears, the systems parameters are expected to change in time. For online estimation of model parameters an Expectation-Maximization (EM) algorithm is used. The EM algorithm is a well established method, which arose in the mathematical statistics community [4], but has found wide engineering applications in different areas. One of the applications is the use of the algorithm for estimating statespace model parameters in the presence of hidden variables (unmeasured states) [5]. The outline of the paper is as follows. Chapter 2 will present the experimental setup used for this study, experimental protocol and feature extraction. Chapter 3 will describe the stochastic time series properties and prediction along with the EM algorithm for model parameter estimation. Results of prediction using the data from the test bed are presented in Chapter 4. II. T HE EXPERIMENTAL S ET- UP The experimental test bed consists of a motor-generator pair with a single stage gearbox (Fig. 1). The motor is a standard DC motor powered through a Simoreg DC drive. A generator is being used as a break. The generated power is being fed back in the system, thus achieving the breaking force. Fig. 1. The test bed The test bed is equipped with 8 accelerometers. The mounting position and sensitivity axis of each accelerometer are shown in Figure 2.
2 Fig. 2. Vibration sensors placement scheme A. The Experimental protocol The test run was done with a constant torque of 82.5 Nm and constant speed of 99rpm. This speed of 99rpm generates GMF f gm = 396Hz, rotational speed of input shaft f i = 16.5Hz, and rotational speed of output shaft f o = 24.75Hz. The signals were sampled with sampling frequency f s = 8kHz. Each acquisition session lasted for 5 seconds. The acquisition was repeated every 1 minutes as illustrated in Figure 3. (a) Output gear (b) Input gear Fig. 4. Output and input gear at the end of the experiment (notice heavily pitted teeth) FE 1 y 1 VIBRATION SENSORS FE 2 y 2 TEST RIG FE N y N Fig. 3. The concept of signal acquisition Fig. 5. Feature extraction procedure In order to speed up the experiment, the contact surface between the gears was decreased to 1/3 of the original surface. In this manner the fault evolution horizon is made shorter. The displacement is shown in Fig. 4. The overall experiment run took 65 hours. At the end both gears were heavily damaged i.e. on both gears spalling can be seen on all teeth, which progressed even into a plastic deformation of some of them,as shown in Fig. 4. B. Feature extraction Signals from all 8 sensors were acquired simultaneously. These signals acquired at each acquisition session, were analyzed using envelope analysis [2]. From each sensor s i {s 1,, s 8 }, at each acquisition session k, a feature vector [y si,1(k), y si,2(k),, y si,m(k)] was derived, where m is the total number of extracted features. Each element of the feature set represents the value of the amplitude of specific spectral component from the envelope spectrum for the particular sensor s i. III. STOCHASTIC TIME SERIES PREDICTION A. State space model of time series State-space representation is a very general model, that can describe a whole set of different models. In our case we assume that condition of the machine is a dynamic process influenced by random tribological inputs which occur due to impact between moving surfaces. Condition can be viewed as a random process, which can be described by a state space model: x t+1 = f(x t,u,w t,θ), y t = g(x t,u,e t,θ) (1) where y t is a measured data, u measured system input x t an unmeasured system state, w t and e t denote process and measurement noise, respectively. Θ denotes a vector of static model parameters. For practical use, the expression (1) can be in many cases simplified and transformed into linear form. w t and e t are
3 assumed to be mutually independent white noise processes with variance Q and R respectively. Additionally, where no distinct measurable inputs are present, u t is considered to be. The resulting model is described as x t+1 = Ax t + w t, y t = Cx t + e t, (2) We assume that the system starts with initial vector x with mean µ and covariance matrix Σ. The observed system output data, indexed by time (y t ) represent the time series we wish to analyze or predict. B. EM algorithm for dynamic state-space system estimation Expectation-Maximization is applied as an iterative method to estimate a vector of unknown parameters (Θ), given measurement data (Y = [y 1, y 2,..., y n ]). In other words, we wish to find the set of parameters Θ, such that the data likelihood function (p(y Θ)) is a maximum. This estimate of Θ is known as maximum likelihood (ML) estimate. Usually, the log-likelihood function of the parameters is defined as, L(Θ) = lnp(y X,Θ), (3) where X = [x 1,x 2,...,x n ]. Since ln(x) is a strictly increasing function, the value of Θ that maximizes p(y X, Θ) also maximizes L(Θ). { Θ k+1 = arg max EX Y,Θk {lnp(y,x Θ)} } (4) Θ EM algorithm is an iterative procedure for maximizing L(Θ), meaning that after k th iteration, we obtain the estimate for Θ, denoted by Θ k. The advantage of the algorithm is, that it can also operate when system states (x) are not known. Because the output is dependant on the unobserved system states, direct maximization is not possible. The EM algorithm alternates between two steps, first maximizing the likelihood function with respect to the system states (E-step) and than with respect to the parameters (M-step). 1) E-step: Given an estimete of the parameter values (Θ k = {A, Q, C, R, x, Q }), Rauch-Tung-Striebel (RTS) smoother provides an optimum estimate of the unobserved state sequence (x t ) of a state space model (2). In dynamical systems with hidden states the E-step corresponds directly to solving the smoothing problem. In other words, for any time t we would like to compute p(x t y 1:T ) (5) where y 1:T = [y 1, y 2,..., y T ], T > t. The RTS smoother procedure is as follows: First, Kalman filter is applied to the observable data in a forward manner (t = 1, 2,..., n), starting with initial estimate of the state mean and variance (x,σ ). TABLE I RTS SMOOTHER ALGORITHM Forward filter Initialization x = x P = P Computation: For t = 1, 2,..., n x t+1 t = Ax t t P t+1 t = AP t t A T + Q K t = P t+1 t C T (CP t+1 t C T + R) 1 x t+1 t+1 = x t+1 t + K t(y t+1 Cx t+1 t ) P t+1 t+1 = P t+1 t K tcp t+1 t Backward filter Computation: For T = n, n 1,...,1 J t = P t t A T P 1 t+1 t x t T = x t t + J t(x t+1 T x t+1 t ) P t T = P t t J t(p t+1 T P t+1 t )J T t Recursive smoother is applied in a backward manner (t = n, n 1,...,1), taking the filtered state estimate at time n as initial condition. Summary of the algorithm is given in Table I. 2) M-step: As stated above, the vector of unknown system parameters in the case of state space model is given as Θ = {A, Q, C, R, x, Σ } (6) If the system states can be observed in addition to system outputs, the joint pdf can be written as p(y Θ) = p(y X, Θ)p(X Θ) Assuming Gaussian distribution of initial state and ignoring the constants, the complete data log-likelihood can be written as 2 lnl(θ) = ln Σ + (x µ ) Σ 1 (x µ ) + n ln Q + (x t Ax t 1 ) Q 1 (x t Ax t 1 ) + n ln R + (y t Cx t ) R 1 (y t Cx t ) (7) where x N (µ,σ ). Taking the expected value of the expression with respect to the current parameter estimate (Θ k ) and complete observed data Y n l(θ Θ k ) = E { 2 lnl(θ) Y n,θ k } (8) Using the results from RTS smoother equations: (x t x t ) = xn t xn t + P n t (9) (x t x t 1 ) = xn t xn t 1 + P n t,t 1 (1) (x t 1 x t 1) = x n t 1x n t 1 + P n t 1 (11) (x t ) = x n t (12)
4 and taking expectation over the expression (7) yields where l(θ Θ k ) = = ln Σ + tr [ Σ 1 (P n + (x n µ )(x n µ ) ) ] + n ln Q + tr [ Q 1 {S 11 S 1 A AS 1 + AS A } ] + n ln R + tr [ R 1 {S 22 S 2 C CS 2 + CS 11 C } ] (13) S 11 = S 1 = S = S 22 = S 2 = x n t x n t + P n t (14) x n t x n t 1 + P n t,t 1 (15) x n t 1 xnt 1 + P n t 1 (16) n E(y t y t ) = (y t y t ) (17) y t x n t (18) Or goal is to find the maximum of the function l(θ Θ n ) with respect to parameters A,C,Q,R,x,Σ. Calculating derivatives of the Eq. (13) with respect to all parameters we obtain the following result A = S 1 S 1 (19) Q = n 1 ( S 11 S 1 S 1 ) S 1 (2) C = S 2 S 1 11 (21) R = n 1 ( S 22 S 2 S 1 ) 11 S 2 (22) x = x n (23) Q = P n (24) Both E and M steps are iterated until the increase in the likelihood at the current time step compared to the previous one is greater than the selected threshold. 3) Time series prediction: With known model parameters, predicting the future values of the time series is straightforward. We start at the time of the last measure, with the filter estimate of the state vector distribution (x n n, P n n ), where n is the time index of the last measured sample. At every future time step t > n, the mean and variance of the output vector y t can be calculated analytically. IV. EXPERIMENTAL RESULTS The algorithm has been used to predict first passage time in the setup as described in section 2. Each feature is represented by a time series of 39 samples. The goal is, to use the algorithm for online prognosis, which is done in the following way. At time t, we estimate the parameters of the underlying stochastic model (Eq. 2) based on data window y (t Nw+1):t, where Nw is window length. Then we use the model to predict future behavior, that is time T when y(t) for the first time crosses the alarm value y i.e. y(t) y. The idea behind the approach is that as the condition of the machine will change over time, so will the values of the model parameters. Model parameters will determine the trend in the feature value, while noise covariance parameters stand for the varying noise component (which increases in the times close to failure). A. Model structure The underlying model is assumed to be of the following form: x 1 (t + 1) = x 2 (t + 1) = a 11 x 1 (t) + a 21 x 2 (t) + w 1 (t) a 22 x 2 (t) + w 2 (t) y(t) = c 1 x 1 (t) + e(t) (25) or in matrix state-space form x t+1 = Ax t + w t, B. Online tracking of model parameters y t = Cx t + e t, (26) Whenever a new measurement is obtained, the algorithm estimates model parameters using the last 1 samples. Initial parameters of the model are set to the following values: [ ] 1 1 A =, 1 [ ] 1 Q =, 1 C = [ 1 ] R = [ 1 ] [ ] x = [ ] 1 Σ =, 1 (27) Convergence of the EM procedure is achieved when the relative increase of the log likelihood function is less than 1 4. Figure 6 shows the measured time series, first eigenvalue of the estimated matrix A and diagonal elements of the estimated noise covariance matrix Q. Up to 4 hours of operation, it can be seen that the feature values are relatively small and there is no trend present. Therefore the estimated system is stable (eigenvalue is less than 1). However, when the feature value starts to increase, the system eigenvalue becomes greater than one.
5 Feature values time series Feature value First eigenvalue of matrix A Diagonal elements of Q 1 x Predicted time [h] (a) Feature measured time 5 4 Fig. 6. Values of the estimated parameters The increase in system noise parameters corresponds to stronger noise component present in the signal, which can be clearly seen from time series (between 4 and 65 hours). Fig. 7. (b) Prediction Results of MC analysis using eight sensor data feature C. Prediction of the first passage time The goal of our prediction is to determine the time, at which the feature will exceed a certain value, we approach in the following way. At a selected time (t), model parameters are estimated using windowed data (y (t Nw+1):t ). Based on the estimated model, a Monte-Carlo simulation is preformed. For every run of MC simulation a time when the predicted feature value becomes greater than a certain threshold is determined. The procedure is repeated at different times and each for each of them the simulation is repeated 1 times. The large number of iteration is necessary due to relatively large value of noise covariances which cause a high dispersion of the results. The final result of the Monte-Carlo simulation analysis is a distribution of the times, at which feature value became greater than the selected threshold for an array of different times. From the full time series, an actual time at which this occurred at the actual machine has been determined. The described results for a feature representing 8 th sensor data (c.f. 2) of the vibrations are shown in Figure 7. It can be seen, that the first estimates of the approximate time are made around 2 hours before the feature actually achieves that value. The estimation variance gradually decreases and, in this case, 15 hours before the variance settles inside a ±5 hours of the true value. To evaluate the performance of our algorithm, we have compared its results to the ones obtained using a simple linear regression. for this, we used a time series consisting of third sensor feature values. The results are shown in Figure 8, where black boxplots represent the results from EM based Feature value 15 1 Predicted time [h] 5 Feature value Safety alarm value (a) Feature (b) Prediction Fig. 8. Comparison of the MC simulation results from third sensor data for EM based alrogithm (black) and linear regression (grey) algorithm and grey boxplots from linear regression. Both times, the simulation has been repeated 1 times. It can be seen, that EM based algorithm provides a more accurate prediction of the first passage time. This results confirm our hypothesis, that the complex dynamics of the feature value, coupled with strong noise component require a more advanced models in order to predict the future time
6 series evolution. V. CONCLUSION It has been shown that second order dynamic linear stochastic state space model is sufficient to model the dynamic behavior of the vibration feature value in the presented setup. Model parameters have been estimated using Expectation-Maximization algorithm, which iteratively maximizes the model likelihood function with respect to hidden states and to unknown parameters. The results show, that the model obtained in this way can effectively predict the future behavior of the feature and can therefore be used to predict the time of safety alarm. From current experiments, we conclude that using this method, the accurate prediction can be made 15 to 2 hours in advance. This offers the machine operators or maintenance a reasonable amount of time to replace the gear system without causing unnecessary production downtime. In this stage, our algorithm has only been tested in the case of constant load. In real application, it is common that the load as well as rounds per minute are changing in time. Next step in the development of the procedure is to modify the algorithm so that it will include these as a measured system input. REFERENCES [1] R. Rubini and U. Meneghetti, Application of the envelope and wavelet transform and analyses for the diagnosis of incipient faults in ball bearings, Mechanical Systems and Signal Processing, vol. 15, pp , 21. [2] D. Ho and R. B. Randall, Optimisation of bearing diagnostic techniques using simulated and actual bearing fault signals, Mechanical Systems and Signal Processing, vol. 14, pp , 2. [3] W. Q. Wang, M. F. Golnaragh, and F. Ismail, Prognosis of machine health condition using neuro-fuzzy systems, Mechanical Systems and Signal Processing, vol. 18, pp , 23. [4] A. Dempster, N. Lard, and D. Rubin, Maximum likelihood from incomplete data via the em algorithm, Journal of the Royal Statistical Society, Series B, vol. 39, pp. 1 38, [5] S. Gibson and B. Ninness, Robust maximum-likelihood estimation of multivariable dynamic systems, Automatica, vol. 41, pp , 25.
Prognosis of gear health using stochastic dynamical models with online parameter estimation
Prognosis of gear health using stochastic dynamical models with online parameter estimation Matej Gašperin 1, Pavle Boškoski 1, and Dani Juričič 1 1 Jožef Stefan Institute, Ljubljana, Slovenia matej.gasperin@ijs.si
More informationModel-Based Prognostics Under Non-stationary Operating Conditions
Model-Based Prognostics Under Non-stationary Operating Conditions Matej Gašperin, Pavle Boškoski, -Dani Juričić Joˇzef Stefan Institute, Ljubljana, Slovenia matej.gasperin@ijs.si, pavle.boskoski@ijs.si,
More informationABSTRACT INTRODUCTION
ABSTRACT Presented in this paper is an approach to fault diagnosis based on a unifying review of linear Gaussian models. The unifying review draws together different algorithms such as PCA, factor analysis,
More informationFactor Analysis and Kalman Filtering (11/2/04)
CS281A/Stat241A: Statistical Learning Theory Factor Analysis and Kalman Filtering (11/2/04) Lecturer: Michael I. Jordan Scribes: Byung-Gon Chun and Sunghoon Kim 1 Factor Analysis Factor analysis is used
More information2198. Milling cutter condition reliability prediction based on state space model
2198. Milling cutter condition reliability prediction based on state space model Hongkun Li 1, Shuai Zhou 2, Honglong Kan 3, Ming Cong 4 School of Mechanical Engineering, Dalian University of Technology,
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond January 18, 2017 Contents 1 Batch and Recursive Estimation 2 Towards Bayesian Filtering 3 Kalman Filter and Bayesian Filtering and Smoothing
More informationMonitoring Radar Mechanical Drive Systems FAA's choice of monitoring solutions
Monitoring Radar Mechanical Drive Systems FAA's choice of monitoring solutions Mission-critical Ask the public to name mission-critical systems, and air traffic control radar will be at the top of the
More informationEM-algorithm for Training of State-space Models with Application to Time Series Prediction
EM-algorithm for Training of State-space Models with Application to Time Series Prediction Elia Liitiäinen, Nima Reyhani and Amaury Lendasse Helsinki University of Technology - Neural Networks Research
More informationBearing fault diagnosis based on EMD-KPCA and ELM
Bearing fault diagnosis based on EMD-KPCA and ELM Zihan Chen, Hang Yuan 2 School of Reliability and Systems Engineering, Beihang University, Beijing 9, China Science and Technology on Reliability & Environmental
More informationIntelligent Fault Classification of Rolling Bearing at Variable Speed Based on Reconstructed Phase Space
Journal of Robotics, Networking and Artificial Life, Vol., No. (June 24), 97-2 Intelligent Fault Classification of Rolling Bearing at Variable Speed Based on Reconstructed Phase Space Weigang Wen School
More informationAbstract. 1. Introduction
Abstract Repairable system reliability: recent developments in CBM optimization A.K.S. Jardine, D. Banjevic, N. Montgomery, A. Pak Department of Mechanical and Industrial Engineering, University of Toronto,
More informationA CUSUM approach for online change-point detection on curve sequences
ESANN 22 proceedings, European Symposium on Artificial Neural Networks, Computational Intelligence and Machine Learning. Bruges Belgium, 25-27 April 22, i6doc.com publ., ISBN 978-2-8749-49-. Available
More informationMachine Learning 4771
Machine Learning 4771 Instructor: ony Jebara Kalman Filtering Linear Dynamical Systems and Kalman Filtering Structure from Motion Linear Dynamical Systems Audio: x=pitch y=acoustic waveform Vision: x=object
More informationChapter 3. Experimentation and Data Acquisition
48 Chapter 3 Experimentation and Data Acquisition In order to achieve the objectives set by the present investigation as mentioned in the Section 2.5, an experimental set-up has been fabricated by mounting
More informationThe Kalman Filter ImPr Talk
The Kalman Filter ImPr Talk Ged Ridgway Centre for Medical Image Computing November, 2006 Outline What is the Kalman Filter? State Space Models Kalman Filter Overview Bayesian Updating of Estimates Kalman
More informationComputational and Experimental Approach for Fault Detection of Gears
Columbia International Publishing Journal of Vibration Analysis, Measurement, and Control (2014) Vol. 2 No. 1 pp. 16-29 doi:10.7726/jvamc.2014.1002 Research Article Computational and Experimental Approach
More informationParameter Estimation in a Moving Horizon Perspective
Parameter Estimation in a Moving Horizon Perspective State and Parameter Estimation in Dynamical Systems Reglerteknik, ISY, Linköpings Universitet State and Parameter Estimation in Dynamical Systems OUTLINE
More informationLecture 6: April 19, 2002
EE596 Pat. Recog. II: Introduction to Graphical Models Spring 2002 Lecturer: Jeff Bilmes Lecture 6: April 19, 2002 University of Washington Dept. of Electrical Engineering Scribe: Huaning Niu,Özgür Çetin
More informationMultilevel Analysis of Continuous AE from Helicopter Gearbox
Multilevel Analysis of Continuous AE from Helicopter Gearbox Milan CHLADA*, Zdenek PREVOROVSKY, Jan HERMANEK, Josef KROFTA Impact and Waves in Solids, Institute of Thermomechanics AS CR, v. v. i.; Prague,
More informationComputer Intensive Methods in Mathematical Statistics
Computer Intensive Methods in Mathematical Statistics Department of mathematics johawes@kth.se Lecture 16 Advanced topics in computational statistics 18 May 2017 Computer Intensive Methods (1) Plan of
More informationDefense Technical Information Center Compilation Part Notice
UNCLASSIFIED Defense Technical Information Center Compilation Part Notice ADP013496 TITLE: Methods to Estimate Machine Remaining Useful Life Using Artificial Neural Networks DISTRIBUTION: Approved for
More informationFE670 Algorithmic Trading Strategies. Stevens Institute of Technology
FE670 Algorithmic Trading Strategies Lecture 3. Factor Models and Their Estimation Steve Yang Stevens Institute of Technology 09/12/2012 Outline 1 The Notion of Factors 2 Factor Analysis via Maximum Likelihood
More informationExpectation Maximization
Expectation Maximization Seungjin Choi Department of Computer Science and Engineering Pohang University of Science and Technology 77 Cheongam-ro, Nam-gu, Pohang 37673, Korea seungjin@postech.ac.kr 1 /
More informationMarkov Chain Monte Carlo Methods for Stochastic Optimization
Markov Chain Monte Carlo Methods for Stochastic Optimization John R. Birge The University of Chicago Booth School of Business Joint work with Nicholas Polson, Chicago Booth. JRBirge U of Toronto, MIE,
More informationParticle Filtering Approaches for Dynamic Stochastic Optimization
Particle Filtering Approaches for Dynamic Stochastic Optimization John R. Birge The University of Chicago Booth School of Business Joint work with Nicholas Polson, Chicago Booth. JRBirge I-Sim Workshop,
More informationCHAPTER 6 FAULT DIAGNOSIS OF UNBALANCED CNC MACHINE SPINDLE USING VIBRATION SIGNATURES-A CASE STUDY
81 CHAPTER 6 FAULT DIAGNOSIS OF UNBALANCED CNC MACHINE SPINDLE USING VIBRATION SIGNATURES-A CASE STUDY 6.1 INTRODUCTION For obtaining products of good quality in the manufacturing industry, it is absolutely
More informationVibration Signals Analysis and Condition Monitoring of Centrifugal Pump
Technical Journal of Engineering and Applied Sciences Available online at www.tjeas.com 2013 TJEAS Journal-2013-3-12/1081-1085 ISSN 2051-0853 2013 TJEAS Vibration Signals Analysis and Condition Monitoring
More informationCHAPTER 4 FAULT DIAGNOSIS OF BEARINGS DUE TO SHAFT RUB
53 CHAPTER 4 FAULT DIAGNOSIS OF BEARINGS DUE TO SHAFT RUB 4.1 PHENOMENON OF SHAFT RUB Unwanted contact between the rotating and stationary parts of a rotating machine is more commonly referred to as rub.
More informationLinear Dynamical Systems
Linear Dynamical Systems Sargur N. srihari@cedar.buffalo.edu Machine Learning Course: http://www.cedar.buffalo.edu/~srihari/cse574/index.html Two Models Described by Same Graph Latent variables Observations
More informationImpeller Fault Detection for a Centrifugal Pump Using Principal Component Analysis of Time Domain Vibration Features
Impeller Fault Detection for a Centrifugal Pump Using Principal Component Analysis of Time Domain Vibration Features Berli Kamiel 1,2, Gareth Forbes 2, Rodney Entwistle 2, Ilyas Mazhar 2 and Ian Howard
More informationMarkov Chain Monte Carlo Methods for Stochastic
Markov Chain Monte Carlo Methods for Stochastic Optimization i John R. Birge The University of Chicago Booth School of Business Joint work with Nicholas Polson, Chicago Booth. JRBirge U Florida, Nov 2013
More informationUnsupervised Learning Methods
Structural Health Monitoring Using Statistical Pattern Recognition Unsupervised Learning Methods Keith Worden and Graeme Manson Presented by Keith Worden The Structural Health Monitoring Process 1. Operational
More informationCS281A/Stat241A Lecture 17
CS281A/Stat241A Lecture 17 p. 1/4 CS281A/Stat241A Lecture 17 Factor Analysis and State Space Models Peter Bartlett CS281A/Stat241A Lecture 17 p. 2/4 Key ideas of this lecture Factor Analysis. Recall: Gaussian
More informationReliability Theory of Dynamic Loaded Structures (cont.) Calculation of Out-Crossing Frequencies Approximations to the Failure Probability.
Outline of Reliability Theory of Dynamic Loaded Structures (cont.) Calculation of Out-Crossing Frequencies Approximations to the Failure Probability. Poisson Approximation. Upper Bound Solution. Approximation
More informationA Data-driven Approach for Remaining Useful Life Prediction of Critical Components
GT S3 : Sûreté, Surveillance, Supervision Meeting GdR Modélisation, Analyse et Conduite des Systèmes Dynamiques (MACS) January 28 th, 2014 A Data-driven Approach for Remaining Useful Life Prediction of
More informationDynamic Approaches: The Hidden Markov Model
Dynamic Approaches: The Hidden Markov Model Davide Bacciu Dipartimento di Informatica Università di Pisa bacciu@di.unipi.it Machine Learning: Neural Networks and Advanced Models (AA2) Inference as Message
More information2256. Application of empirical mode decomposition and Euclidean distance technique for feature selection and fault diagnosis of planetary gearbox
56. Application of empirical mode decomposition and Euclidean distance technique for feature selection and fault diagnosis of planetary gearbox Haiping Li, Jianmin Zhao, Jian Liu 3, Xianglong Ni 4,, 4
More informationLecture 6: Bayesian Inference in SDE Models
Lecture 6: Bayesian Inference in SDE Models Bayesian Filtering and Smoothing Point of View Simo Särkkä Aalto University Simo Särkkä (Aalto) Lecture 6: Bayesian Inference in SDEs 1 / 45 Contents 1 SDEs
More informationA Study of Covariances within Basic and Extended Kalman Filters
A Study of Covariances within Basic and Extended Kalman Filters David Wheeler Kyle Ingersoll December 2, 2013 Abstract This paper explores the role of covariance in the context of Kalman filters. The underlying
More informationHidden Markov models for time series of counts with excess zeros
Hidden Markov models for time series of counts with excess zeros Madalina Olteanu and James Ridgway University Paris 1 Pantheon-Sorbonne - SAMM, EA4543 90 Rue de Tolbiac, 75013 Paris - France Abstract.
More informationImproving Electromotor Process in Water Pump by Using Power Spectral Density, Time Signal and Fault Probability Distribution Function
Improving Electromotor Process in Water Pump by Using Power Spectral Density, Time Signal and Fault Probability Distribution Function Hojjat Ahmadi, Zeinab Khaksar Department of Agricultural Machinery
More informationX t = a t + r t, (7.1)
Chapter 7 State Space Models 71 Introduction State Space models, developed over the past 10 20 years, are alternative models for time series They include both the ARIMA models of Chapters 3 6 and the Classical
More informationLatent Variable Models and EM algorithm
Latent Variable Models and EM algorithm SC4/SM4 Data Mining and Machine Learning, Hilary Term 2017 Dino Sejdinovic 3.1 Clustering and Mixture Modelling K-means and hierarchical clustering are non-probabilistic
More informationState Feedback Controller for Position Control of a Flexible Link
Laboratory 12 Control Systems Laboratory ECE3557 Laboratory 12 State Feedback Controller for Position Control of a Flexible Link 12.1 Objective The objective of this laboratory is to design a full state
More information1 Kalman Filter Introduction
1 Kalman Filter Introduction You should first read Chapter 1 of Stochastic models, estimation, and control: Volume 1 by Peter S. Maybec (available here). 1.1 Explanation of Equations (1-3) and (1-4) Equation
More informationVibration Analysis of Shaft Misalignment and Diagnosis Method of Structure Faults for Rotating Machinery
Available online at www.ijpe-online.com Vol. 3, No. 4, July 07, pp. 337-347 DOI: 0.3940/ijpe.7.04.p.337347 Vibration Analysis of Shaft Misalignment and Diagnosis Method of Structure Faults for Rotating
More informationVibroacoustic source separation using an improved cyclic Wiener filter
Vibroacoustic source separation using an improved cyclic Wiener filter Konstantinos C. Gryllias 1, Jérôme Antoni 1 and Mario Eltabach 2 1 LVA, ISA-Lyon, University of Lyon, F-69621, France 2 CETIM, F-60300,
More informationChapter 4: Factor Analysis
Chapter 4: Factor Analysis In many studies, we may not be able to measure directly the variables of interest. We can merely collect data on other variables which may be related to the variables of interest.
More informationFAULT DIAGNOSIS IN ROTATING MACHINE USING FULL SPECTRUM OF VIBRATION AND FUZZY LOGIC
Journal of Engineering Science and Technology Vol. 12, No. 11 (2017) 2952-2964 School of Engineering, Taylor s University FAULT DIAGNOSIS IN ROTATING MACHINE USING FULL SPECTRUM OF VIBRATION AND FUZZY
More informationConstrained State Estimation Using the Unscented Kalman Filter
16th Mediterranean Conference on Control and Automation Congress Centre, Ajaccio, France June 25-27, 28 Constrained State Estimation Using the Unscented Kalman Filter Rambabu Kandepu, Lars Imsland and
More informationMachine learning - HT Maximum Likelihood
Machine learning - HT 2016 3. Maximum Likelihood Varun Kanade University of Oxford January 27, 2016 Outline Probabilistic Framework Formulate linear regression in the language of probability Introduce
More informationSensor Tasking and Control
Sensor Tasking and Control Sensing Networking Leonidas Guibas Stanford University Computation CS428 Sensor systems are about sensing, after all... System State Continuous and Discrete Variables The quantities
More informationParticle Filtering for Data-Driven Simulation and Optimization
Particle Filtering for Data-Driven Simulation and Optimization John R. Birge The University of Chicago Booth School of Business Includes joint work with Nicholas Polson. JRBirge INFORMS Phoenix, October
More informationDetection of bearing faults in high speed rotor systems
Detection of bearing faults in high speed rotor systems Jens Strackeljan, Stefan Goreczka, Tahsin Doguer Otto-von-Guericke-Universität Magdeburg, Fakultät für Maschinenbau Institut für Mechanik, Universitätsplatz
More informationCS534 Machine Learning - Spring Final Exam
CS534 Machine Learning - Spring 2013 Final Exam Name: You have 110 minutes. There are 6 questions (8 pages including cover page). If you get stuck on one question, move on to others and come back to the
More informationOptimization of the detection of train wheel defects. SNCF Innovation and Research Department Paris, FRANCE 1
Optimization of the detection of train wheel defects 1 R. Ziani SNCF Innovation and Research Department Paris, FRANCE 1 Abstract This paper describes how statistical models and learning algorithms could
More informationApplication of Vector Time Series Modeling and T-squared Control Chart to Detect Early Gearbox Deterioration
International Journal of Performability Engineering, Vol.10, No. 1, January 2014, pp.105-114. RAMS Consultants Printed in India Application of Vector Time Series Modeling and T-squared Control Chart to
More informationA methodology for fault detection in rolling element bearings using singular spectrum analysis
A methodology for fault detection in rolling element bearings using singular spectrum analysis Hussein Al Bugharbee,1, and Irina Trendafilova 2 1 Department of Mechanical engineering, the University of
More informationWHEELSET BEARING VIBRATION ANALYSIS BASED ON NONLINEAR DYNAMICAL METHOD
15 th November 212. Vol. 45 No.1 25-212 JATIT & LLS. All rights reserved. WHEELSET BEARING VIBRATION ANALYSIS BASED ON NONLINEAR DYNAMICAL METHOD 1,2 ZHAO ZHIHONG, 2 LIU YONGQIANG 1 School of Computing
More informationSmoothers: Types and Benchmarks
Smoothers: Types and Benchmarks Patrick N. Raanes Oxford University, NERSC 8th International EnKF Workshop May 27, 2013 Chris Farmer, Irene Moroz Laurent Bertino NERSC Geir Evensen Abstract Talk builds
More informationHMM and IOHMM Modeling of EEG Rhythms for Asynchronous BCI Systems
HMM and IOHMM Modeling of EEG Rhythms for Asynchronous BCI Systems Silvia Chiappa and Samy Bengio {chiappa,bengio}@idiap.ch IDIAP, P.O. Box 592, CH-1920 Martigny, Switzerland Abstract. We compare the use
More informationUnsupervised Learning with Permuted Data
Unsupervised Learning with Permuted Data Sergey Kirshner skirshne@ics.uci.edu Sridevi Parise sparise@ics.uci.edu Padhraic Smyth smyth@ics.uci.edu School of Information and Computer Science, University
More informationoutput dimension input dimension Gaussian evidence Gaussian Gaussian evidence evidence from t +1 inputs and outputs at time t x t+2 x t-1 x t+1
To appear in M. S. Kearns, S. A. Solla, D. A. Cohn, (eds.) Advances in Neural Information Processing Systems. Cambridge, MA: MIT Press, 999. Learning Nonlinear Dynamical Systems using an EM Algorithm Zoubin
More informationNote Set 5: Hidden Markov Models
Note Set 5: Hidden Markov Models Probabilistic Learning: Theory and Algorithms, CS 274A, Winter 2016 1 Hidden Markov Models (HMMs) 1.1 Introduction Consider observed data vectors x t that are d-dimensional
More information2262. Remaining life prediction of rolling bearing based on PCA and improved logistic regression model
2262. Remaining life prediction of rolling bearing based on PCA and improved logistic regression model Fengtao Wang 1, Bei Wang 2, Bosen Dun 3, Xutao Chen 4, Dawen Yan 5, Hong Zhu 6 1, 2, 3, 4 School of
More informationAn EM algorithm for Gaussian Markov Random Fields
An EM algorithm for Gaussian Markov Random Fields Will Penny, Wellcome Department of Imaging Neuroscience, University College, London WC1N 3BG. wpenny@fil.ion.ucl.ac.uk October 28, 2002 Abstract Lavine
More informationA New Approach to Tune the Vold-Kalman Estimator for Order Tracking
A New Approach to Tune the Vold-Kalman Estimator for Order Tracking Amadou Assoumane, Julien Roussel, Edgard Sekko and Cécile Capdessus Abstract In the purpose to diagnose rotating machines using vibration
More informationCrib Sheet : Linear Kalman Smoothing
Crib Sheet : Linear Kalman Smoothing Gabriel A. Terejanu Department of Computer Science and Engineering University at Buffalo, Buffalo, NY 14260 terejanu@buffalo.edu 1 Introduction Smoothing can be separated
More information1 What is a hidden Markov model?
1 What is a hidden Markov model? Consider a Markov chain {X k }, where k is a non-negative integer. Suppose {X k } embedded in signals corrupted by some noise. Indeed, {X k } is hidden due to noise and
More informationDual Estimation and the Unscented Transformation
Dual Estimation and the Unscented Transformation Eric A. Wan ericwan@ece.ogi.edu Rudolph van der Merwe rudmerwe@ece.ogi.edu Alex T. Nelson atnelson@ece.ogi.edu Oregon Graduate Institute of Science & Technology
More informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Dan Cawford Manfred Opper John Shawe-Taylor May, 2006 1 Introduction Some of the most complex models routinely run
More informationBearings-Only Tracking in Modified Polar Coordinate System : Initialization of the Particle Filter and Posterior Cramér-Rao Bound
Bearings-Only Tracking in Modified Polar Coordinate System : Initialization of the Particle Filter and Posterior Cramér-Rao Bound Thomas Brehard (IRISA/CNRS), Jean-Pierre Le Cadre (IRISA/CNRS). Journée
More informationAlgorithmisches Lernen/Machine Learning
Algorithmisches Lernen/Machine Learning Part 1: Stefan Wermter Introduction Connectionist Learning (e.g. Neural Networks) Decision-Trees, Genetic Algorithms Part 2: Norman Hendrich Support-Vector Machines
More informationIn-Flight Engine Diagnostics and Prognostics Using A Stochastic-Neuro-Fuzzy Inference System
In-Flight Engine Diagnostics and Prognostics Using A Stochastic-Neuro-Fuzzy Inference System Dan M. Ghiocel & Joshua Altmann STI Technologies, Rochester, New York, USA Keywords: reliability, stochastic
More informationLaboratory 11 Control Systems Laboratory ECE3557. State Feedback Controller for Position Control of a Flexible Joint
Laboratory 11 State Feedback Controller for Position Control of a Flexible Joint 11.1 Objective The objective of this laboratory is to design a full state feedback controller for endpoint position control
More informationPractical Bayesian Optimization of Machine Learning. Learning Algorithms
Practical Bayesian Optimization of Machine Learning Algorithms CS 294 University of California, Berkeley Tuesday, April 20, 2016 Motivation Machine Learning Algorithms (MLA s) have hyperparameters that
More informationSystem Identification, Lecture 4
System Identification, Lecture 4 Kristiaan Pelckmans (IT/UU, 2338) Course code: 1RT880, Report code: 61800 - Spring 2016 F, FRI Uppsala University, Information Technology 13 April 2016 SI-2016 K. Pelckmans
More informationPartially Observable Markov Decision Processes (POMDPs)
Partially Observable Markov Decision Processes (POMDPs) Sachin Patil Guest Lecture: CS287 Advanced Robotics Slides adapted from Pieter Abbeel, Alex Lee Outline Introduction to POMDPs Locally Optimal Solutions
More informationA nonlinear dynamic vibration model of defective bearings: The importance of modelling the finite size of rolling elements
A nonlinear dynamic vibration model of defective bearings: The importance of modelling the finite size of rolling elements Alireza Moazenahmadi, Dick Petersen and Carl Howard School of Mechanical Engineering,
More informationWe Prediction of Geological Characteristic Using Gaussian Mixture Model
We-07-06 Prediction of Geological Characteristic Using Gaussian Mixture Model L. Li* (BGP,CNPC), Z.H. Wan (BGP,CNPC), S.F. Zhan (BGP,CNPC), C.F. Tao (BGP,CNPC) & X.H. Ran (BGP,CNPC) SUMMARY The multi-attribute
More informationTowards a Bayesian model for Cyber Security
Towards a Bayesian model for Cyber Security Mark Briers (mbriers@turing.ac.uk) Joint work with Henry Clausen and Prof. Niall Adams (Imperial College London) 27 September 2017 The Alan Turing Institute
More informationGear Fault Identification by using Vibration Analysis
International Journal of Current Engineering and Technology E-ISSN 2277 4106, P-ISSN 2347 5161 2016 INPRESSCO, All Rights Reserved Available at http://inpressco.com/category/ijcet Research Article Ajanalkar
More informationStatistical Estimation
Statistical Estimation Use data and a model. The plug-in estimators are based on the simple principle of applying the defining functional to the ECDF. Other methods of estimation: minimize residuals from
More informationLecture 7: Optimal Smoothing
Department of Biomedical Engineering and Computational Science Aalto University March 17, 2011 Contents 1 What is Optimal Smoothing? 2 Bayesian Optimal Smoothing Equations 3 Rauch-Tung-Striebel Smoother
More informationProbabilistic Graphical Models
Probabilistic Graphical Models Brown University CSCI 2950-P, Spring 2013 Prof. Erik Sudderth Lecture 13: Learning in Gaussian Graphical Models, Non-Gaussian Inference, Monte Carlo Methods Some figures
More informationObservation and analysis of the vibration and displacement signature of defective bearings due to various speeds and loads
Observation and analysis of the vibration and displacement signature of defective bearings due to various speeds and loads Alireza-Moazen ahmadi a) Carl Howard b) Department of Mechanical Engineering,
More informationAn Acoustic Emission Approach to Assess Remaining Useful Life of Aging Structures under Fatigue Loading PI: Mohammad Modarres
An Acoustic Emission Approach to Assess Remaining Useful Life of Aging Structures under Fatigue Loading PI: Mohammad Modarres Outline Objective Motivation, acoustic emission (AE) background, approaches
More informationExpectation propagation for signal detection in flat-fading channels
Expectation propagation for signal detection in flat-fading channels Yuan Qi MIT Media Lab Cambridge, MA, 02139 USA yuanqi@media.mit.edu Thomas Minka CMU Statistics Department Pittsburgh, PA 15213 USA
More informationSwitching Regime Estimation
Switching Regime Estimation Series de Tiempo BIrkbeck March 2013 Martin Sola (FE) Markov Switching models 01/13 1 / 52 The economy (the time series) often behaves very different in periods such as booms
More informationON NUMERICAL ANALYSIS AND EXPERIMENT VERIFICATION OF CHARACTERISTIC FREQUENCY OF ANGULAR CONTACT BALL-BEARING IN HIGH SPEED SPINDLE SYSTEM
ON NUMERICAL ANALYSIS AND EXPERIMENT VERIFICATION OF CHARACTERISTIC FREQUENCY OF ANGULAR CONTACT BALL-BEARING IN HIGH SPEED SPINDLE SYSTEM Tian-Yau Wu and Chun-Che Sun Department of Mechanical Engineering,
More informationLecture 4: Hidden Markov Models: An Introduction to Dynamic Decision Making. November 11, 2010
Hidden Lecture 4: Hidden : An Introduction to Dynamic Decision Making November 11, 2010 Special Meeting 1/26 Markov Model Hidden When a dynamical system is probabilistic it may be determined by the transition
More informationUsing Operating Deflection Shapes to Detect Misalignment in Rotating Equipment
Using Operating Deflection Shapes to Detect Misalignment in Rotating Equipment Surendra N. Ganeriwala (Suri) & Zhuang Li Mark H. Richardson Spectra Quest, Inc Vibrant Technology, Inc 8205 Hermitage Road
More informationDynamic System Identification using HDMR-Bayesian Technique
Dynamic System Identification using HDMR-Bayesian Technique *Shereena O A 1) and Dr. B N Rao 2) 1), 2) Department of Civil Engineering, IIT Madras, Chennai 600036, Tamil Nadu, India 1) ce14d020@smail.iitm.ac.in
More informationLecture 8: Bayesian Estimation of Parameters in State Space Models
in State Space Models March 30, 2016 Contents 1 Bayesian estimation of parameters in state space models 2 Computational methods for parameter estimation 3 Practical parameter estimation in state space
More informationDynamic Tests on Ring Shear Apparatus
, July 1-3, 2015, London, U.K. Dynamic Tests on Ring Shear Apparatus G. Di Massa Member IAENG, S. Pagano, M. Ramondini Abstract Ring shear apparatus are used to determine the ultimate shear strength of
More informationGaussian with mean ( µ ) and standard deviation ( σ)
Slide from Pieter Abbeel Gaussian with mean ( µ ) and standard deviation ( σ) 10/6/16 CSE-571: Robotics X ~ N( µ, σ ) Y ~ N( aµ + b, a σ ) Y = ax + b + + + + 1 1 1 1 1 1 1 1 1 1, ~ ) ( ) ( ), ( ~ ), (
More informationROLLER BEARING FAILURES IN REDUCTION GEAR CAUSED BY INADEQUATE DAMPING BY ELASTIC COUPLINGS FOR LOW ORDER EXCITATIONS
ROLLER BEARIG FAILURES I REDUCTIO GEAR CAUSED BY IADEQUATE DAMPIG BY ELASTIC COUPLIGS FOR LOW ORDER EXCITATIOS ~by Herbert Roeser, Trans Marine Propulsion Systems, Inc. Seattle Flexible couplings provide
More informationEarly Fault Detection and Optimal Maintenance Control for Partially Observable Systems Subject to Vibration Monitoring
Early Fault Detection and Optimal Maintenance Control for Partially Observable Systems Subject to Vibration Monitoring by Chen Lin A thesis submitted in conformity with the requirements for the degree
More informationIEOR E4570: Machine Learning for OR&FE Spring 2015 c 2015 by Martin Haugh. The EM Algorithm
IEOR E4570: Machine Learning for OR&FE Spring 205 c 205 by Martin Haugh The EM Algorithm The EM algorithm is used for obtaining maximum likelihood estimates of parameters when some of the data is missing.
More information: Probabilistic Engineering Analysis and Design Professor Youn, Byeng Dong
CHAPTER 9. HEALTH DIAGNOSTICS AND PROGNOSTICS 9.1 Introduction Last several decades, tremendous advance has been made on the physics-based analysis and design under uncertainties. However, it is still
More information