Time-domain Markov chain Monte Carlo–based Bayesian damage detection of ballasted tracks using nonlinear ballast stiffness model

2020 ◽  
pp. 147592172096695
Author(s):  
Heung-Fai Lam ◽  
Mujib Olamide Adeagbo ◽  
Yeong-Bin Yang

This article reports the development of a methodology for detecting ballast damage under a sleeper based on measured sleeper vibration following the Bayesian statistical system identification framework. To ensure the methodology is applicable under large amplitude vibration of the sleeper (e.g. under trainload), the nonlinear stress–strain behavior of railway ballast is considered. This, on one hand, significantly reduces the problem of modeling error, but, on the other hand, increases the number of uncertain model parameters. The uncertainty associated with the identified model parameters of the rail–sleeper–ballast system may be very high. To overcome this difficulty, the Markov chain Monte Carlo–based Bayesian model updating is adopted in the proposed methodology for the approximation of the posterior probability density function of uncertain model parameters. Owing to the nonlinear behavior of the system, the model updating is performed in the time domain instead of the modal domain. The applicability of the proposed damage detection methodology was first verified numerically using simulated impact hammer test data in two damaged cases perturbed with Gaussian white noise. Second, impact hammer tests of in situ sleepers in the full-scale in-door ballasted track test panel were carried out to collect data for the experimental verification of the proposed methodology. Artificial ballast damage was simulated under the target concrete sleeper by replacing normal-sized ballast particles (∼60 mm) by small-sized ballast particles (∼15 mm). The proposed methodology successfully identified the location and severity of ballast damage under the sleeper. From the calculated posterior marginal probability density functions of model parameters, one can quantify the uncertainties associated with the damage detection results. The proposed methodology is an essential step in the development of a long-term railway track health monitoring system utilizing train-induced vibration.

2017 ◽  
Vol 17 (3) ◽  
pp. 706-724 ◽  
Author(s):  
Heung F Lam ◽  
Jia H Yang ◽  
Qin Hu ◽  
Ching T Ng

This article reports the development of a Bayesian method for assessing the damage status of railway ballast under a concrete sleeper based on vibration data of the in situ sleeper. One of the important contributions of the proposed method is to describe the variation of stiffness distribution of ballast using Lagrange polynomial, for which the order of the polynomial is decided by the Bayesian approach. The probability of various orders of polynomial conditional on a given set of measured vibration data is calculated. The order of polynomial with the highest probability is selected as the most plausible order and used for updating the ballast stiffness distribution. Due to the uncertain nature of railway ballast, the corresponding model updating problem is usually unidentifiable. To ensure the applicability of the proposed method even in unidentifiable cases, a computational efficient Markov chain Monte Carlo–based Bayesian method was employed in the proposed method for generating a set of samples in the important region of parameter space to approximate the posterior (updated) probability density function of ballast stiffness. The proposed ballast damage detection method was verified with roving hammer test data from a segment of full-scale ballasted track. The experimental verification results positively show the potential of the proposed method in ballast damage detection.


2008 ◽  
Vol 10 (2) ◽  
pp. 153-162 ◽  
Author(s):  
B. G. Ruessink

When a numerical model is to be used as a practical tool, its parameters should preferably be stable and consistent, that is, possess a small uncertainty and be time-invariant. Using data and predictions of alongshore mean currents flowing on a beach as a case study, this paper illustrates how parameter stability and consistency can be assessed using Markov chain Monte Carlo. Within a single calibration run, Markov chain Monte Carlo estimates the parameter posterior probability density function, its mode being the best-fit parameter set. Parameter stability is investigated by stepwise adding new data to a calibration run, while consistency is examined by calibrating the model on different datasets of equal length. The results for the present case study indicate that various tidal cycles with strong (say, >0.5 m/s) currents are required to obtain stable parameter estimates, and that the best-fit model parameters and the underlying posterior distribution are strongly time-varying. This inconsistent parameter behavior may reflect unresolved variability of the processes represented by the parameters, or may represent compensational behavior for temporal violations in specific model assumptions.


2017 ◽  
Vol 14 (18) ◽  
pp. 4295-4314 ◽  
Author(s):  
Dan Lu ◽  
Daniel Ricciuto ◽  
Anthony Walker ◽  
Cosmin Safta ◽  
William Munger

Abstract. Calibration of terrestrial ecosystem models is important but challenging. Bayesian inference implemented by Markov chain Monte Carlo (MCMC) sampling provides a comprehensive framework to estimate model parameters and associated uncertainties using their posterior distributions. The effectiveness and efficiency of the method strongly depend on the MCMC algorithm used. In this work, a differential evolution adaptive Metropolis (DREAM) algorithm is used to estimate posterior distributions of 21 parameters for the data assimilation linked ecosystem carbon (DALEC) model using 14 years of daily net ecosystem exchange data collected at the Harvard Forest Environmental Measurement Site eddy-flux tower. The calibration of DREAM results in a better model fit and predictive performance compared to the popular adaptive Metropolis (AM) scheme. Moreover, DREAM indicates that two parameters controlling autumn phenology have multiple modes in their posterior distributions while AM only identifies one mode. The application suggests that DREAM is very suitable to calibrate complex terrestrial ecosystem models, where the uncertain parameter size is usually large and existence of local optima is always a concern. In addition, this effort justifies the assumptions of the error model used in Bayesian calibration according to the residual analysis. The result indicates that a heteroscedastic, correlated, Gaussian error model is appropriate for the problem, and the consequent constructed likelihood function can alleviate the underestimation of parameter uncertainty that is usually caused by using uncorrelated error models.


2002 ◽  
Vol 6 (5) ◽  
pp. 883-898 ◽  
Author(s):  
K. Engeland ◽  
L. Gottschalk

Abstract. This study evaluates the applicability of the distributed, process-oriented Ecomag model for prediction of daily streamflow in ungauged basins. The Ecomag model is applied as a regional model to nine catchments in the NOPEX area, using Bayesian statistics to estimate the posterior distribution of the model parameters conditioned on the observed streamflow. The distribution is calculated by Markov Chain Monte Carlo (MCMC) analysis. The Bayesian method requires formulation of a likelihood function for the parameters and three alternative formulations are used. The first is a subjectively chosen objective function that describes the goodness of fit between the simulated and observed streamflow, as defined in the GLUE framework. The second and third formulations are more statistically correct likelihood models that describe the simulation errors. The full statistical likelihood model describes the simulation errors as an AR(1) process, whereas the simple model excludes the auto-regressive part. The statistical parameters depend on the catchments and the hydrological processes and the statistical and the hydrological parameters are estimated simultaneously. The results show that the simple likelihood model gives the most robust parameter estimates. The simulation error may be explained to a large extent by the catchment characteristics and climatic conditions, so it is possible to transfer knowledge about them to ungauged catchments. The statistical models for the simulation errors indicate that structural errors in the model are more important than parameter uncertainties. Keywords: regional hydrological model, model uncertainty, Bayesian analysis, Markov Chain Monte Carlo analysis


Sign in / Sign up

Export Citation Format

Share Document