Railway ballast damage detection by Markov chain Monte Carlo-based Bayesian method

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.