Next Article in Journal
Comparison of Different Sets of Features for Human Activity Recognition by Wearable Sensors
Next Article in Special Issue
Air Quality Monitoring for Vulnerable Groups in Residential Environments Using a Multiple Hazard Gas Detector
Previous Article in Journal
Qualitative Dynamics of Suspended Particulate Matter in the Changjiang Estuary from Geostationary Ocean Color Images: An Empirical, Regional Modeling Approach
Previous Article in Special Issue
An Ontology-Driven Approach for Integrating Intelligence to Manage Human and Ecological Health Risks in the Geospatial Sensor Web
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Model Order Identification for Cable Force Estimation Using a Markov Chain Monte Carlo-Based Bayesian Approach

1
College of Computer Science and Software Engineering, Shenzhen University, Shenzhen 51086, China
2
Institute of Intelligence Cyber Sensing System, Shenzhen Academy of Robotics, Shenzhen 518057, China [email protected] (S.Z.)
*
Authors to whom correspondence should be addressed.
Submission received: 26 September 2018 / Revised: 21 November 2018 / Accepted: 25 November 2018 / Published: 29 November 2018
(This article belongs to the Special Issue Computational Intelligence-Based Sensors)

Abstract

:
The tensile force on the hanger cables of a suspension bridge is an important indicator of the structural health of the bridge. Tensile force estimation methods based on the measured frequency of the hanger cable have been widely used. These methods empirically pre-determinate the corresponding model order of the measured frequency. However, because of the uncertain flexural rigidity, this empirical order determination method not only plays a limited role in high-order frequencies, but also hinders the online cable force estimation. Therefore, we propose a new method to automatically identify the corresponding model order of the measured frequency, which is based on a Markov chain Monte Carlo (MCMC)-based Bayesian approach. It solves the limitation of empirical determination in the case of large flexural rigidity. The tensile force and the flexural rigidity of cables can be calculated simultaneously using the proposed method. The feasibility of the proposed method is validated via a numerical study involving a finite element model that considers the flexural rigidity and via field application to a suspension bridge.

Graphical Abstract

1. Introduction

As one of the primary load-bearing components of a cable system bridge, such as a cable-stayed bridge or a suspension bridge, the cable force should be adjusted during construction and monitored when the bridge is in service. Therefore, it is essential that the cable force is accurately measured. At present, the vibration method, according to which cable forces are estimated from the measured vibration frequency, is the most widely used method [1,2,3] for measuring cable forces due to its simplicity, speed, and economy relative to the method according to which cable forces are estimated via magneto-elastic permeability [4] or from a cutter ring containing strain gauges.
Currently, vibration methods that ignore cable sag are widely used to determine the tensile force of the hanger cables of a suspension bridge because the sag of vertical cables is so small as to hardly influence the measured frequency. The cable tension can be estimated based on taut string theory [5,6,7] or by established formulas for specified cables [8,9].
However, methods based on a single cable model are prone to errors when flexural rigidity values are large. Therefore, inverse analysis and system identification methods to estimate the hanger cable force based on finite element models (FEMs) have recently been proposed. These methods match the theoretical frequency from the numerical model with the measured frequency from an existing hanger cable. Park [10] first built a hanger cable FEM to obtain the modal frequency in an ideal geometric condition and then used the back analysis technique to optimize the variables using measured frequencies. Papadimitriou [11] discussed two different categories of model classified by different end plate connections and boundary conditions, using the Bayesian approach to estimate the cable force. Xie [12] used a genetic algorithm to optimize this inverse problem.
All of the aforementioned cable force estimation methods require the model order to be determined empirically, the accuracy of which depends on expert understanding of the structure. Generally, for measured frequencies of a low order, such as the first three orders, it is typically easy to specify the order from the similarity between the measured and theoretical frequencies from taut string theory or the interval between measured frequencies. However, in the case of high flexural rigidity, the empirical method is limited because of the sensitive effect of flexural rigidity, and a rapid increase occurs in the ratio between the measured frequency and higher-order theoretical frequencies [5,13]. In addition, in some cases, the frequency spectrum of cable vibration may indicate that low regions have a very small response or are submerged by noise frequency related to the low order, and that high regions have distinct dominant frequencies related to the high order. Hence, empirical determination plays a limited role in the case of high flexural rigidity, which hinders real-time cable force estimation.
To address this problem, we propose a new method for automatic model order identification that simultaneously estimates the cable force and the flexural rigidity of the hanger cable. The new method regards the model order as an unknown parameter, along with the flexural rigidity and cable force. This approach involves both continuous and discrete variables, and a stochastic method is adopted that can well-represent this sort of mixed mode. Here, a popular stochastic method, the Bayesian approach [14], is used. The feasibility of the proposed method is validated via two cases: a numerical study involving an FEM considering the flexural rigidity and a field application to a suspension bridge.

2. Dynamics of a Cable with the Flexural Rigidity

By ignoring cable sag, taking flexural rigidity into account, and having simply supported ends, the dynamic equation of a cable can be described as Equation (1) via the dynamic displacement of gravity direction v ( x ,   t ) at time t and length x [15]:
N 2 v ( x , t ) x 2 E I 4 v ( x , t ) x 4 = m 2 v ( x , t ) t 2
where N is the cable force, E I is the flexural rigidity, and m is the mass per unit length of the cable. For different conditions, various deformation equations can be derived from Equation (1) as a partial differential equation. Among them, Shimada [7] proposed a linear relationship (2) between ( f k / k ) 2 and k 2 , where k is the model order and f k is the k-th nature frequency, which is widely used for cable force estimation.
( f k k ) 2 = ( 1 4 m L 2 ) N + ( k 2 π 2 4 m L 4 ) E I
To scale the flexural stiffness effect on cables, a dimensionless parameter ξ is introduced, expressed as
ξ = L N E I

3. Model Order Identification

In the proposed method, the model order is regarded as an unknown parameter in the cable vibration model. In addition, this vibration model has another two unknown parameters—cable force and flexural rigidity. The input is multiple measured frequencies extracted from the frequency spectrum. The range of model order is limited considering a completely slack cable. Finally, the three parameters are optimized through Bayesian estimation. However, an integral is required for each item of the Bayesian formula. We need to balance the number of integral points, as insufficiency will result in an obscure or erroneous output, and excessive data will severely prolong the computational time. The MCMC method can be used to complete this integral task. For a random sampling point (parameter vector), the MCMC method determines acceptance or rejection, and another point is resampled based on the previously accepted point. The Bayesian approach is a feature of the acceptance mechanism.

3.1. Frequency-Domain Analysis

Measured frequencies of the cable are determined by identifying the peaks from the power spectrum (PS) [1]. Researchers often use fast Fourier transform (FFT) in PS implementation. However, the calculation of Fourier spectra overlooks that a signal is a sample of a stochastic process; as such, much of the signal content is merely noise. Thus, it is a completely different task to obtain spectra from a theoretical signal versus a real cable vibration signal, the latter of which involves a stochastic process. Recently, output-only system identification for modal analysis, such as frequency domain decomposition (FDD) [16] and stochastic space identification (SSI) [17], has received considerable attention. Compared with the single FFT, these methods have the advantage of the PS being more distinguishable when the test structure is excited by ambient forces.
Here, FDD is adopted for frequency extraction. The power spectral density matrix of outputs of the structure can be defined as a matrix G y y ( j ω ) as discrete frequencies ω = ω i . Then, to obtain the measured frequencies, the power spectral density matrix in each frequency is decomposed by singular value decomposition, i.e., G y y ( j ω i ) =   U i S i U i H , where matrix U i = [ u i 1 , u i 2 , , u i m ] contains the singular values u i j , and matrix S i contains the singular values s i j . The dominant peaks of these singular values belong to the measured frequencies.

3.2. Model Order Range Limitation

Under the effect of flexural rigidity, the measured differs from the theoretical frequency obtained from taut string theory. This difference is significant if the flexural rigidity is high. As a result, it is challenging to confirm the order of measured frequencies. This section describes a method to limit the range of the model order and to provide optimal methods via the measured frequencies. An appropriate range of these orders must be provided. The measured frequency can be expressed further by the theoretical frequency value f ˙ from taut string theory given by
f k ¯ 2 = f ˙ k 2 ( 1 + C )
where C = ( k L ) 2 E I N . If the flexural rigidity is low, then it only minimally influences the measured frequency, i.e., the minimum value of the constant C is 0. Conversely, if the flexural rigidity is large, then the measured frequency is more strongly influenced by higher model orders for the same cable length. Consider a condition E I / N = β , where β represents a completely slack cable. Given a constraint k L / β , the range of C can be defined within 0   C 1 . Thus, the range of measured natural frequency is given as
f k ¯ 2 2 f ˙ j 2
where 2 f ˙ j 2 represents the nearest values to f k ¯ 2 . Because of the sequence, the measured frequencies should not share the same order; moreover, measured frequencies of high value should have a higher order than those of low value. Finally, for a set of measured frequencies with number s, the range of the order of the s-th measured frequency k s can be defined as
{ j + 1 k L / β k s < k s + 1     ( k Z )

3.3. Bayesian Parameter Estimation

After several measured frequencies are extracted from the PS, a parameterized model is built as follows:
f k s ¯ 2 k i 2 4 m L 2 N k i 4 π 2 4 m L 4 E I = ε i
where f k s ¯ is the s-th measured frequency, for which the range of the order k i is defined in the last section, and ε i is the i-th error value. The Euclidean norm of ε i is minimized to invert the measured frequency for the model parameters (model order, cable force, and flexural rigidity). This inversion problem can be solved using a deterministic method, such as the Gauss-Newton method, or using a stochastic method, such as the Bayesian approach. Because of the nonlinearity of the forward modeling and the general non-uniqueness and ill-posed nature of the inversion problem [18], many local optimal solutions may exist in this high-dimensional inversion problem. Relative to the deterministic method, the stochastic method provides a global solution and a superior approach to quantify uncertainty in the inversion problem. More importantly, the stochastic method has the capacity to represent mixed models with continuous and discrete variables. Therefore, a stochastic method—the Bayesian approach—is adopted here.
The Bayesian approach is based on Bayes’ theorem. This theorem states that, given the observed data, the probability of the model (also called the posterior probability function) depends on the likelihood function, the prior distribution, and the evidence. The Bayesian approach has been applied in many fields and has been proved to be a powerful approach to solving high-dimension inversion problems, such as finite-element-model updating [19,20,21,22] and artificial neural networks. The fundamental rule that governs the Bayesian approach is described as
P ( θ | D ) = P ( D | θ ) P ( θ ) P ( D )
where θ is a vector of the updating parameters; P ( θ ) is the probability distribution function of the updating parameters in the absence of any data, which is known as the prior distribution; D is the model output; P ( D ) is a normalization factor; and P ( D | θ ) is the likelihood probability distribution function. In addition, the quantity P ( θ | D ) is the posterior probability distribution function following a set of model outputs.
In this study, given the number of measured frequencies s, the vector of updating parameters θ is formed as θ = { k 1 , , k s , N , E I } . The prior distribution function P ( θ ) reflects prior information, which is defined as
P ( θ ) = i = 1 s + 2 P ( θ i )
In neural networks, the Gaussian prior has been successfully used to identify a large unknown parameter between two nodes, also known as weight; therefore, we assume this prior can successfully be used to identify a small number of parameters in this study. However, we find that the uniform distribution is more appropriate for the prior distribution of the model order, perhaps due to the integer property and sparse values of the model order. Studies have considered the solution for the prior distribution of integer values [15,23] in a linear model. In this study, the method described by [24] is used with the value range decided by Equation (9). The likelihood probability distribution function P ( D | θ ) is defined as the sum of squares of the error vector and is written as
P ( D | θ ) = 1 Z D e x p ( i = 1 s ε i 2 )
where Z D = e x p ( i = 1 n ε i 2 ) d θ . The posterior distribution is still proportional to the numerator of Equation (8) after ignoring the normalization factor P ( D ) .

3.4. Markov Chain Monte Carlo (MCMC) Method

It is apparent that both the prior distribution and the likelihood probability distribution are based on integral calculation. Due to the relatively high dimension of the updating parameter vector, it is challenging to solve the posterior distribution analytically. Thus, an MCMC method is used to avoid multi-dimensional integration by the sampling based on the posterior distribution.
MCMC methods belong to a class of algorithms for sampling from a probability distribution. These methods are based on building a Markov chain that has the desired distribution as its equilibrium distribution. After multiple iterations using the specific accept/rejection rule, the state of the chain is used as a sample of the desired distribution. Considering that the conditional distributions of some components in this study are nonstandard [25], the Metropolis-Hastings (M-H) algorithm is used as the sampler. Step-by-step instructions for the MH algorithm used in this study are presented in Figure 1.
Note that the prior distribution of the model order is the uniform distribution when the prior distribution of other updating parameters is the normal distribution. Because the acceptance ratio is related to the comparison between the posterior distribution of the proposal parameter vector and the next predicted parameter vector, considering the normalization factor is a constant, the normalization factor is ignored and does not require calculation.
In this study, we select a Gelman-Rubin diagnostic [26] to monitor the convergence of MCMCs. The Gelman-Rubin diagnostic is a general approach to monitoring the convergence of MCMC multi-chain outputs with starting values that are over-dispersed relative to the posterior distribution. The diagnostic assumes normality of the marginal posterior distributions. Convergence is assessed by comparing the estimated between-chains and within-chain variances for each model parameter. Large differences between these variances indicate non-convergence.

4. Feasibility Validation

4.1. Numerical Study

This section verifies our proposed method by judging whether the known parameters (model order, cable force, and flexural rigidity) are equal to the solutions for the frequencies from a simulated cable. Consider a cable without sag and with both ends hinged, as shown in Figure 2. The cable prototype is a hanger cable described previously [27]. Its mass per unit length is 13.6 kg/m and the diameter of the cross section is 56 mm. Four cables with different flexural rigidity values are used, with flexural rigidities of 63%, 158%, 317%, and 31% of the ideal value (Table 1). A linear FEM of the cable is programmed in MATLAB for eigenvalue analyses. The measured frequencies of the four cables are listed in Table 2. Four consecutive frequencies of each cable are selected to validate the proposed technique; these values correspond to the bold values in Table 2. The model order of the selected frequency is regarded as an unknown parameter. The four consecutive frequencies are matched with parameters of model order K1, K2, K3, and K4. Both the cable force and the flexural rigidity are initialized from 0.3, with an iteration number of 10,000.
Figure 3 shows the convergence of the model order of Cable 1. The convergence of the cable force and the flexural rigidity in the four cables are shown in Figure 4, demonstrating good agreement between the actual value under the ideal cable model and the ultimate converged value. Uncertainty in the estimation, expressed as the standard deviation (2σ), is denoted by the gray shaded area along the corresponding curves. With the cable force and flexural rigidity approaching the true values, the estimation uncertainty tends to decrease as the accuracy of the estimation improves. However, for Cable 4, when the cable force approaches the true value, the flexural rigidity wanders between 0.1 and 2.5. This wandering is caused by the low level of flexural rigidity. When flexural rigidity is low, the natural frequency is influenced only minimally, thus allowing the identified flexural rigidity to take on values over a wide range.
Because of the integer property, the curve of the chain values is far from smooth. The model order ultimately converges to the actual value. Comparing the convergence of cable force and flexural rigidity in Cable 1 shown in Figure 4(a1,a2), the model order rapidly converges to a stable condition when the cable force and flexural rigidity require additional iteration counts to reach the converging lines.

4.2. Field Application

The Huangpu suspension bridge, one of the 20 longest suspension bridges in the world, spans the main channel of the Pearl River, carrying the highway linking Caotang and Bicun in Guangzhou, China. It includes 144 hangers. To verify the proposed technique, we select hanger Cable N2 in the bridge as the study cable. Empirically estimating the model order is straightforward due to the low flexural rigidity of the cable. The properties of Cable N2 are shown in Table 3. Under an ambient vibration condition, we measure the acceleration time responses of the cable via the wireless sensor nodes with two-axis MEMS accelerometers attached along the longitudinal direction on the cable, as shown in Figure 5. The sensor nodes are wirelessly connected, have a sampling frequency of 100 Hz, and are synchronized via GPS timing. Figure 6a shows the acceleration responses of the cable under ambient vibration, and its PS is presented in Figure 6b.
FDD combined with Welch’s averaged and modified periodogram method is used to obtain the frequency spectrum. The window size for dividing the signal into sections is 1024, with an overlap number of 512. FDD can reduce the noise components in the frequency spectrum, as highlighted in Figure 7. The peak frequencies in the spectrum using FDD hold the same values as those in the single Fourier spectrum, whereas the sidelobes are decreased. Because of the large flexural rigidity, the interval between these frequencies is approximately regarded as an integral multiple of the first-order nature frequency. Thus, the model orders marked in Figure 6b are empirically determined to correspond to 16, 18, 20, 22, and 24. Next, these five frequencies are used to verify the proposed method, with orders 16, 18, 20, 22, and 24 associated with parameters K1, K2, K3, K4, and K5, respectively.
Figure 8 shows the convergence of cable force using the proposed method without predefining the model orders and the Gauss-Newton method with the predefined orders. Figure 9 shows the convergence of the five unknown model orders. The cable tensions estimated by taut string theory using peak frequencies of either weak or strong amplitude are represented by the green area, which is also based on the empirical predefined model order.
According to Figure 8 and Figure 9, the proposed method can effectively estimate the cable tension while automatically identifying the order of measured frequencies. The Gauss-Newton method requires the pre-determination of an accurate model order; otherwise, the estimated cable force has a large error. The proposed method has a large error in the cable force at the beginning of the iteration process that is caused by the large range of model order. After many iterations, the cable force converges to the same value as in the Gauss-Newton method. Histograms of the obtained samples of cable force and flexural rigidity are presented in Figure 10. The mean values of cable force and flexural rigidity are 0.258 MN and 0.065 MN·m2 when the standard deviations are 0.0031 MN and 0.0076 MN·m2, respectively.
The Gelman-Rubin diagnostic presents the scale reduction for each parameter by the shrink factor, as well as the auto-correlation. Figure 11 shows that the shrink factor of parameter N becomes stable after about 5000 iterations and wanders around 1 when the shrink factor of EI is less than 2. For both parameter N and EI, the plotted auto-correlation curve quickly decays to zero and then wanders around it. This indicates a good chain with a rapid rate and low degree of correlation between our samples, and also reveals that MCMC sampling has been done in an independent manner and the stationarity of the Markov chain has been attained. Based on the good performance for an iteration number of 10,000, we run MCMC for 10 and 20 times 10,000 again and find that the three runs have almost the same shape of histogram and convergence performance.

5. Conclusions

An automatic method for model order identification is proposed as an alternative to the existing method of cable force estimation, which uses empirical determination of the order of measured frequencies. The proposed method simultaneously calculates the cable force under the effect of flexural rigidity and includes the following features: (1) a measured frequency extraction procedure via Welch’s frequency spectrum, (2) a method to determine the range of the model orders of these frequencies, and (3) a Bayesian approach using these ranged orders as unknown parameters to determine both the cable force and the flexural rigidity. The Bayesian approach used here adopts an MCMC sampler, and the orders are given a uniform prior distribution. The feasibility of the proposed method is verified via a numerical study using an FEM model of cables and a field experiment, in which accelerometers are installed on hangers in a suspension bridge to collect the cable vibration data. The result of the field application shows that the proposed method can accurately calculate the model order and estimate the cable force, providing a value that is the same as the empirical value and the value obtained from the Gauss-Newton method. Note that the latter method must pre-define the model order, whereas the former method does not.
The proposed method can also work on more or less model order frequencies. In the field application, we use five model orders as unknown parameters to investigate the capacity of the number of model orders in our proposed method, and we repeat the method 50 times for different numbers of model order. The result is shown in Figure 12. The proposed method works well when the number of the model order is less than six. However, the performance becomes worse with a greater number of model orders. It is difficult to find an appropriate indicator to show the effect of different initializations on the identification. So, for cable force and flexural rigidity, we use 100 Gaussian random values in the same orders of magnitude of the predicted value. For instance, if the cable force is predicted to be 0.2 MN, the initialized value can be set from 0.1 to 0.99. Four model orders are chosen and the values are randomly taken from the defined limited formula. The accuracy rate of the proposed method for this random initialization is 94%.

Author Contributions

Conceptualization, G.Z. and Z.L.; Methodology, G.Z. and S.Z.; Software, S.Z. and Y.L.; Validation, J.H., S.Z. and Y.L.; Formal Analysis, S.Z.; Supervision, G.Z. and Z.L.

Funding

The authors gratefully acknowledge support from the National Natural Science Foundation of China (grant number: 61673278) and the Science and Technology Innovation Commission of Shenzhen Municipality (SZSTI) (grant numbers: JCYJ20150828104330541, JCYJ20170818105455447, JCYJ20170818105431734, and KQTD20140630154026047).

Acknowledgments

The authors would like to thank Meng Chen and Key Laboratory of Earthquake Monitoring and Disaster Reduction Technology of Guangdong Province, China for their help in preparing the hardware and environment for our experiments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Russell, J.C.; Lardner, T.J. Experimental Determination of Frequencies and Tension for Elastic Cables. J. Eng. Mech. 1998, 124, 1067–1072. [Google Scholar] [CrossRef]
  2. Mordini, A.; Savov, K.; Wenzel, H. Damage Detection on Stay Cables using an Open Source-based Framework for Finite Element Model Updating. Struct. Heal. Monit. 2008, 7, 91–102. [Google Scholar] [CrossRef]
  3. Ma, L. A highly precise frequency-based method for estimating the tension of an inclined cable with unknown boundary conditions. J. Sound Vib. 2017, 409, 65–80. [Google Scholar]
  4. Wang, M.L.; Chen, Z.L.; Koontz, S.S.; Lloyd, G.M. Magnetoelastic permeability measurement for stress monitoring in steel tendons and cables. Proc. SPIE. Int. Soc. Opt. Eng. 2000, 3995, 492–500. [Google Scholar]
  5. Casas, J.R. A combined method for measuring cable forces: The cable-stayed Alamillo Bridge, Spain. Struct. Eng. Int. 1994, 4, 235–240. [Google Scholar] [CrossRef]
  6. Smith, S.W.; Campbell, J.E. Testing and model verification of the Maysville Kentucky Bridge stay cables. In Proceedings of the XX International Modal Analysis Conference, Los Angeles, CA, USA, 4–7 February 2002; pp. 1050–1056. [Google Scholar]
  7. Shimada, T. A Study on the Maintenance and Management of the Tensile Measurement for the Cable of Bridge. Ph.D. Thesis, Kobe University, Kobe, Japan, 1995. [Google Scholar]
  8. Zui, H.; Shinke, T.; Namita, Y.H. Practical formulas for estimation of cable tension by vibration method. J. Struct. Eng. 1996, 122, 651–656. [Google Scholar] [CrossRef]
  9. Fang, Z.; Wang, J. Practical Formula for Cable Tension Estimation by Vibration Method. J. Bridg. Eng. 2012, 17, 161–164. [Google Scholar] [CrossRef]
  10. Park, D.; Kim, N. Back analysis technique for tensile force on hanger cables of a suspension bridge. J. Vib. Control 2014, 20, 761–772. [Google Scholar] [CrossRef]
  11. Papadimitriou, C.; Giakoumi, K.; Argyris, C.; Spyrou, L.A.; Panetsos, P. Bayesian estimation of tension in bridge hangers using modal frequency measurements. Struct. Monit. Maint. 2016, 3, 349–375. [Google Scholar] [CrossRef]
  12. Xie, X.; Li, X. Genetic Algorithm-based tension identification of hanger by solving inverse eigenvalue problem. Inverse Probl. Sci. Eng. 2014, 22, 966–987. [Google Scholar] [CrossRef]
  13. Mehrabi, A.B.; Tabatabai, H. Unified Finite Difference Formulation for Free Vibration. J. Struct. Eng. 1998, 124, 1313–1322. [Google Scholar] [CrossRef]
  14. Sansò, F.; Venuti, G. Integer variables estimation problems: The Bayesian approach. Ann. Geophys. 1997, 40, 1415–1431. [Google Scholar]
  15. Irvine, H.M. Cable Structures; The MIT Press: Cambridge, MA, USA, 1981. [Google Scholar]
  16. Brien, O.; Eugene, J.; Malekjafarian, A.; Obrien, E. Application of output-only modal method in monitoring of bridges using an instrumented vehicle. In Proceedings of the Civil Engineering Research in Ireland Conference, Belfast, UK, 28–29 August 2014. [Google Scholar]
  17. Ubertini, F.; Comanducci, G.; Cavalagli, N. Vibration-based structural health monitoring of a historic bell-tower using output-only measurements and multivariate statistical analysis. Struct. Heal. Monit. 2016, 15, 438–457. [Google Scholar] [CrossRef]
  18. Chen, J.; Kemna, A.; Hubbard, S.S. A comparison between Gauss-Newton and Markov chain Monte Carlo based methods for inverting spectral induced polarization data for Cole-Coleparameters. Geophysics 2008, 73, 247–259. [Google Scholar] [CrossRef]
  19. Beck, J.L. Bayesian system identification based on probability logic. Struct. Ctrl. Heal. Monit. 2010, 17, 825–847. [Google Scholar] [CrossRef]
  20. Erazo, K.; Nagarajaiah, S. An offline approach for output-only Bayesian identification of stochastic nonlinear systems using unscented Kalman filtering. J. Sound Vib. 2017, 397, 222–240. [Google Scholar] [CrossRef]
  21. Li, Z.; Feng, M.Q.; Luo, L.; Feng, D.; Xu, X. Statistical analysis of modal parameters of a suspension bridge based on Bayesian spectral density approach and shm data. Mech. Sys. Signal Proc. 2018, 98, 352–367. [Google Scholar] [CrossRef]
  22. Sun, H.; Büyüköztürk, O. Bayesian Model Updating Using Incomplete Modal Data without Mode Matching. Proc. SPIE 2016, 9805, 98050D. [Google Scholar]
  23. Barger, K.; Bunge, J. Objective Bayesian estimation for the number of species. Bayesian Anal. 2010, 5, 765–786. [Google Scholar] [CrossRef]
  24. Garcia, J.G.; Roncagliolo, P.A.; Carlos, H. A Bayesian technique for real and integer parameters estimation in linear models and its application to GNSS high precision positioning. IEEE Trans. Signal Process. 2016, 64, 923–933. [Google Scholar] [CrossRef]
  25. Liang, F.; Liu, C.; Carroll, R. Advanced Markov Chain Monte Carlo Methods: Learning from Past Samples; Wiley: New York, NY, USA, 2011; ISBN 1119956803. [Google Scholar]
  26. Gelman, A.; Rubin, D.B. Inference from iterative simulation using multiple sequences. Stat. Sci. 1992, 7, 457–511. [Google Scholar] [CrossRef]
  27. Zhan, S.; Chen, M.; Li, Z.; Sun, W.; Li, W.J.; Zhang, G. Large bridge cable tension estimation with cyber vibration sensing system. In Proceedings of the IEEE International Conference on Cyber Technology in Automation, Control, and Intelligent Systems, Chengdu, China, 19–22 June 2016; pp. 121–125. [Google Scholar]
Figure 1. M-H algorithm for MCMC sampling.
Figure 1. M-H algorithm for MCMC sampling.
Sensors 18 04187 g001
Figure 2. Cable model with hinged ends.
Figure 2. Cable model with hinged ends.
Sensors 18 04187 g002
Figure 3. Convergence of the model orders (a) K1, (b) K2, (c) K3, and (d) K4 of Cable 1. K1, K2, K3, and K4 are assumed as unknown parameters which are related to the four consecutive frequencies from the simulated cable.
Figure 3. Convergence of the model orders (a) K1, (b) K2, (c) K3, and (d) K4 of Cable 1. K1, K2, K3, and K4 are assumed as unknown parameters which are related to the four consecutive frequencies from the simulated cable.
Sensors 18 04187 g003
Figure 4. Convergence of the cable force and the flexural rigidity of (a) Cable 1, (b) Cable 2, (c) Cable 3, and (d) Cable 4.
Figure 4. Convergence of the cable force and the flexural rigidity of (a) Cable 1, (b) Cable 2, (c) Cable 3, and (d) Cable 4.
Sensors 18 04187 g004
Figure 5. Wireless sensor nodes with two-axis MEMS accelerometers attached on hanger cables.
Figure 5. Wireless sensor nodes with two-axis MEMS accelerometers attached on hanger cables.
Sensors 18 04187 g005
Figure 6. (a) Time history and (b) frequency spectrum of the acceleration response under ambient vibration.
Figure 6. (a) Time history and (b) frequency spectrum of the acceleration response under ambient vibration.
Sensors 18 04187 g006
Figure 7. Difference between the single Fourier spectrum and the spectrum obtained using Welch’s method: (a) comparison of the whole spectra and (b) magnified sections.
Figure 7. Difference between the single Fourier spectrum and the spectrum obtained using Welch’s method: (a) comparison of the whole spectra and (b) magnified sections.
Sensors 18 04187 g007
Figure 8. Cable force convergence in the proposed method without predefined model orders and the Gauss-Newton method with predefined orders.
Figure 8. Cable force convergence in the proposed method without predefined model orders and the Gauss-Newton method with predefined orders.
Sensors 18 04187 g008
Figure 9. Convergence of the model order.
Figure 9. Convergence of the model order.
Sensors 18 04187 g009
Figure 10. Histogram of the obtained samples of (a) flexural rigidity and (b) cable force.
Figure 10. Histogram of the obtained samples of (a) flexural rigidity and (b) cable force.
Sensors 18 04187 g010
Figure 11. Shrink factor of (a) cable force and (c) flexural rigidity; with ACF of (b) cable force and (d) flexural rigidity in the last 5000 iterations.
Figure 11. Shrink factor of (a) cable force and (c) flexural rigidity; with ACF of (b) cable force and (d) flexural rigidity in the last 5000 iterations.
Sensors 18 04187 g011
Figure 12. The accuracy performance for different numbers of model order.
Figure 12. The accuracy performance for different numbers of model order.
Sensors 18 04187 g012
Table 1. Structural properties of four cables with no sag.
Table 1. Structural properties of four cables with no sag.
Cable No.ξN (MN)EI (×106 N·m2)L (m)
120.00.80.210
240.00.320.550
370.70.51100
4198.00.3920.1100
Table 2. Measured frequencies of the simulated cables.
Table 2. Measured frequencies of the simulated cables.
Cable No.Model Order
23456789
124.238.354.5 (K1)73.4 (K2)95.2 (K3)120 (K4)148.4180
234.5 (K1)6.1 (K2)7.8 (K3)9.7 (K4)11.713.816.1
31.82.83.74.75.7 (K1)6.7 (K2)7.8 (K3)8.8 (K4)
41.6 (K1)2.4 (K2)3.2 (K3)4.1 (K4)4.95.76.57.3
Table 3. Properties of the selected hanger cable.
Table 3. Properties of the selected hanger cable.
Length (m)Diameter (mm)Elastic Modulus (Pa)Theoretical Sectional Area (mm2)Mass per Unit Length (kg/m)
110561.1 × 1011149013.6

Share and Cite

MDPI and ACS Style

Zhan, S.; Li, Z.; Hu, J.; Liang, Y.; Zhang, G. Model Order Identification for Cable Force Estimation Using a Markov Chain Monte Carlo-Based Bayesian Approach. Sensors 2018, 18, 4187. https://0-doi-org.brum.beds.ac.uk/10.3390/s18124187

AMA Style

Zhan S, Li Z, Hu J, Liang Y, Zhang G. Model Order Identification for Cable Force Estimation Using a Markov Chain Monte Carlo-Based Bayesian Approach. Sensors. 2018; 18(12):4187. https://0-doi-org.brum.beds.ac.uk/10.3390/s18124187

Chicago/Turabian Style

Zhan, Shaodong, Zhi Li, Jianmin Hu, Yiping Liang, and Guanglie Zhang. 2018. "Model Order Identification for Cable Force Estimation Using a Markov Chain Monte Carlo-Based Bayesian Approach" Sensors 18, no. 12: 4187. https://0-doi-org.brum.beds.ac.uk/10.3390/s18124187

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop