Abstract
The remaining service life of a bridge is an important criterion for the bridge master to make decisions of maintenance and rehabilitation. In the assessment of fatigue problems of steel and composite bridges, a crucial step is to determine the fatigue stress spectra, which are usually obtained through field measurements. In this study, the standard fatigue vehicle spectrum of Ting Kau Bridge in Hong Kong is obtained from the dynamic weigh-in-motion system that monitors the traffic flow and volume and axle and gross vehicle weights. An advanced probabilistic loading model is proposed to simulate the vehicles running along the bridge based on measurements from the weigh-in-motion system considering not only the diurnal variations of traffic flow within different time intervals but also the rate of change of the annual traffic. The statistical fatigue reliability assessment at critical locations can then be carried out to estimate the remaining service life of the bridge.
Keywords
Introduction
The ever-evolving developments in design and construction technologies have made spans of over 1000 m possible for cable-stayed bridges. Steel is normally used in the decks either fully or partly to reduce the self-weight for those with long spans. In the design and construction of cable-stayed bridges, the design life is often 100 years or more. Over the service life, the effects of environmental corrosion, material ageing, long-term load effect, fatigue and so on will inevitably lead to accumulation of damage and possible reduction of strength. In extreme cases, catastrophic failures could happen. Fatigue has been one of the most critical forms of damage for cable-stayed bridges with steel or composite decks, and assessing the remaining fatigue lives of steel bridges has been the primary focus in recent years (Liu et al., 2010). Fatigue of steel girders is caused by cumulative damage under cyclic actions of various kinds of vehicle loads with different weights and physical dimensions.
In the assessment of fatigue problems, a crucial step is to determine the fatigue stress spectra. Usually, the fatigue stress spectra of bridges can be obtained through field measurements such as from structural health monitoring (SHM) systems which are able to record the responses for fatigue assessment and performance prediction (Elkordy et al., 1994). However, SHM systems are always expensive, and there are a lot of structural components on which it is difficult to install strain gauges (Guo and Chen, 2013). As a result, some other effective tools are required to obtain the fatigue stress spectra of the details that are hardly accessible. In the method of Liu and Mahadevan (2007) for predicting stochastic fatigue life under variable amplitude loading, the S-N curve is stochastic and the fatigue damage accumulation rule is nonlinear, making it computationally intensive for long-span cable-stayed bridges.
A statistical load model was proposed by Guo and Chen (2013) recently so that the uncertain information of vehicle types, number of axles, axle weights, axle spacing, transverse position of vehicle and so on could be properly accounted for. As a sequel to this statistical load model, a statistical approach for modelling traffic load is developed based on the probabilistic traffic characteristics extracted from the measured data of weigh-in-motion (WIM) system installed on Ting Kau Bridge in Hong Kong. Lim (1992) noted the diurnal variations in data measured by WIM systems and the correlation between the average steering axle mass for different hours of the day. Vanderstaay et al. (2006) revealed the characteristics that the mass magnitude depended strongly on the axle group type and that the mass during the night was greater than the mass during the day generally. Righiniotis and Chryssanthopoulos (2003) and Righiniotis (2006) pointed out that traffic growth should be taken into account in fatigue life prediction because of the effects of vehicle loads. The variation of annual traffic has been studied for Runyang Bridge in China, giving rise to a traffic growth model to fit the growth trend (Guo et al., 2008). The traffic load model proposed in this study attempts to reflect the diurnal variations of traffic flow within different time intervals and the rate of change of the annual traffic on a bridge (Zhang et al., 2014). In conjunction with this statistical traffic load model, the stress time-histories at critical locations of the bridge are obtained, and the probability density curves of the fatigue life can be constructed for further analyses.
WIM system
Ting Kau Bridge is a cable-stayed bridge in Hong Kong with a total length of 1177 m completed in 1998. As shown in Figure 1, the Ting Kau main span and Tsing Yi main span are 448 and 475 m, respectively, while the two side spans are both 127 m. The bridge connects Ting Kau at the north and Tsing Yi Island at the south over Rambler Channel. The WIM system collects information of vehicle speed, heading direction, lane used, axle weight, axle spacing between neighbouring axles and so on. The system is capable of dynamic measurement of vehicular axle loads and speed with a programmed vehicle classification system. The WIM system is located at the Tsing Yi end of Ting Kau Bridge on both carriageways as shown in Figure 2.

Schematic plan and elevation of Ting Kau Bridge: (a) plan showing girder labelling and (b) elevation.

Layout of weigh-in-motion dynamic weigh-bridge system: (a) subsystem for detection of over-weigh vehicles and (b) layout on lanes.
For data cleansing, a record of an individual vehicle moving along Ting Kau Bridge is considered unreliable if
Any one axle spacing or axle weight is zero or negative;
The axle spacing is greater than 16 m;
The vehicle speed is less than 50 km/h;
The gross vehicle weight (GVW) is greater than the maximum legislated weight for a vehicle in Hong Kong, which is 44 tons;
Any one axle weight is less than 0.5 ton with axle spacing larger than 2 m (Enright and O’Brien, 2011).
The data of WIM in 84 days in 24-h records are filtered. After filtering, 5% of vehicles are eliminated from the database.
Variation of traffic flow
Diurnal variation of traffic flow
In order to investigate the features of diurnal variations of WIM data for Ting Kau Bridge, the WIM data collected are divided into eight types according to the number of axles of vehicles moving along the bridge, namely, from 2 to 9 axles as shown in Figure 3. The numbers of the eight vehicle types are counted for each hour, whose results are shown in Figure 4, where n1 to n8 denote the first to the eighth vehicle types, respectively. There are two peaks of the numbers of operating vehicles in a day. The two periods around the peaks are called the morning rush hours and evening rush hours. To take into account the variation of daily transportation properly, the daily traffic in this study is divided into four time intervals for further analysis, including 07:00 to 11:00 (morning rush hours), 11:00 to 17:00, 17:00 to 21:00 (evening rush hours) and 21:00 to 07:00 on the next day. The proportion of vehicles of different types in 1 h is analysed and shown in Table 1 and Figure 5. Figure 5 shows that vehicles are predominantly of type 1 which has two axles, while vehicles of types 3–6 constitute a minor portion. Vehicles of types 7 and 8 are ignored in further analysis because of their extremely small proportion.

Vehicle type classification for Ting Kau Bridge: (a) vehicle type 1, (b) vehicle type 2, (c) vehicle type 3, (d) vehicle type 4, (e) vehicle type 5, (f) vehicle type 6, (g) vehicle type 7 and (h) vehicle type 8.

Diurnal variations of numbers of vehicles for Ting Kau Bridge: (a) n1, (b) n2–n5 and (c) n6–n8.
Proportion of vehicles of different types in 1 h for the four time intervals.

Schematic diagram of percentage of each vehicle type in 1 h for the four time intervals.
Variation of annual traffic
Ting Kau Bridge was opened to traffic in 1999 and the traffic has gained rapid growth in recent years as shown in Table 2. About 5% annual growth of traffic can be observed as shown in Figure 6. In other words, the traffic flow will be doubled in roughly 20 years unless there are other constraints such as traffic flow capacity. As the design life of bridges is 120 years as specified by British Standards (BS) 5400 (1980), this annual growth rate of traffic cannot be ignored; otherwise, any estimates of fatigue life will be grossly erroneous.
Annual traffic flow of Ting Kau Bridge.

Annual traffic growth of Ting Kau Bridge.
Statistical fatigue life estimation
Statistical load model
Considering the random variables of each vehicle type, the probability density functions (PDFs) of axle weights and axle spacing in four different time intervals of each day can be obtained. Some plots of the distribution considering the diurnal variation for different time intervals are shown in Figures 7 to 9. Furthermore, the parameters for all distributions are listed in Tables 3 and 4. The symbol ‘(1)’ in the first column is used to denote those for the first time interval and so on. It is found that most of the axle spacing of the six vehicle types follows the t location-scale distribution, while most of the axle weights follow the normal and lognormal distributions.

Distribution of axle weight for the first axle of vehicle type 2 in four time intervals: (a) 7:00–11:00, (b) 11:00–17:00, (c) 17:00–21:00 and (d) 21:00–7:00 next day.

Distribution of axle weight for the third axle of vehicle type 3 in four time intervals: (a) 7:00–11:00, (b) 11:00–17:00, (c) 17:00–21:00 and (d) 21:00–7:00 next day.

Distribution of axle spacing for the first axle of vehicle type 1 in four time intervals: (a) 7:00–11:00, (b) 11:00–17:00, (c) 17:00–21:00 and (d) 21:00–7:00 next day.
Parameters of distribution for random axle spacing of six vehicle types in the first time interval.
Parameters of distribution for random axle weight of six vehicle types in the first time interval.
When the statistical loading is applied for fatigue life estimation, the parameters of each individual vehicle and the arrangement of vehicles in each lane follow certain statistical distributions derived from the traffic data measured by the WIM system. Table 5 shows the proportion of traffic volume for different vehicle types and different lanes at different time intervals. The statistical traffic load model is employed to obtain the stress time-histories at critical locations of the bridge. Based on the simulated statistical stress time-histories, rain-flow counting (Downing, 1982) is used to calculate the number of statistical cycles. The probabilistic distributions of fatigue lives can be calculated by using the statistical distribution of number of cycles. The results of fatigue life obtained using the proposed method are therefore also uncertain. Therefore, Monte Carlo simulation can be used to simulate the large number of occurrences so that the variations can be revealed. Based on the principle of the mathematical statistics, the Monte Carlo method as a kind of numerical method adopting statistical sampling theory solves mathematical and physical problems approximately. Hence, it can also be considered as a statistical experiment or stochastic simulation. The mean stress effect may cause significant difference in fatigue life assessment because, with the revised stress ranges, the cycles to fatigue vary greatly. When considering the mean stress effect, the fatigue life can also be predicted (Zhang and Au, 2013).
Proportion of traffic in each lane.
Basic procedures of fatigue analysis integrating WIM data
The probabilistic analysis of the fatigue lives of Ting Kau Bridge can be performed by combining with the statistical traffic load model according to the following procedures as shown in Figure 10.
1. Establishment of statistical load model
The axle weights and axle spacing of each vehicle type for different time intervals are treated as random variables, whose probability density distributions can be obtained by curve fitting.
2. Calculation of influence line
A moving unit load is applied on different lanes (i.e. slow, middle and fast lanes) of the baseline finite element model. As a result, influence lines of stresses at various critical locations can be obtained for different structural components. A number of load steps corresponding to a static analysis are defined for the bridge.
3. Fatigue analysis (a) Determination of vehicle type
Table 5 shows the proportion of traffic volume for different vehicle types at different time intervals. Based on the information, vehicle samples with the amount equal to the number of vehicles crossing the bridge in each time interval can be simulated. The vehicle type for each sample can be determined. A random number Tvehicle uniformly distributed in the interval [0, 1] is generated first and then the vehicle type is determined by comparing the value of Tvehicle with the cumulative density of different kinds of vehicle. For example, for Tvehicle ≤ 0.9023 in the first time interval, loads of vehicle type 1 will be applied, while loads of vehicle type 2 are used for 0.9023 < Tvehicle ≤ 0.9534.
(b) Selection of lane information
Table 5 also shows the proportion of the traffic volume for different lanes at different time intervals. Vehicle samples travelling along different lanes during each time interval can be obtained. A random number Tlane uniformly distributed in the interval [0, 1] is generated first and then the lane is determined by comparing the value of Tlane with the cumulative density of different lanes. For example, if loads of vehicle type 1 are to be applied and Tlane ≤ 0.155, the simulated vehicle is applied on the first lane, while the vehicle is applied on the second lane for 0.155 < Tlane ≤ 0.4407.
(c) Generation of axle weights and spacing samples
Once the vehicle type and the lane information are known, the next target is to generate the samples of axle weights and spacing for each vehicle. The random number following normal PDF, lognormal PDF and t location distribution can be generated by the truncated Latin Hypercube sampling method (Helton and Davis, 2003), and the correlations between the axle weights and axle spacing can be neglected to simplify the analysis.
(d) Calculation of stress time-histories for each vehicle sample
According to the vehicle type, the lane, axle loads and the axle spacing of the sample determined in the above steps, loads of a vehicle are generated. The influence lines of the details concerned for the specified lane are used for determining the stress time-histories under the moving axles obtained. The results associated with the influence lines obtained from Step (2) are multiplied by the axle weights, and the results of all axles are superimposed according to the axle spacing, so as to obtain the stress time-histories.
(e) Assessment of fatigue life considering mean stress effect
According to the simulated stress time-histories, rain-flow counting is carried out in the MATLAB environment to calculate the number of cycles. The steps from (a) to (e) should be repeated for each sample of vehicle in different time intervals. The number of repetitions should be equal to the number of vehicles crossing the bridge. The daily accumulated number of cycles can be used to estimate the fatigue life for different critical locations considering the mean stress effect.
4. Statistical distribution of fatigue life
The Monte Carlo method is used here to repeat the procedures of fatigue analysis for a significant number of times so that a significant number of random fatigue lives can be obtained, from which the statistical distribution can be obtained using the fitting approach. A total of 100 repetitions are adopted here.
5. Statistical fatigue life considering annual growth of traffic
According to the stress time-histories simulated from the influence lines at various critical locations and the statistical traffic information, the traditional Monte Carlo method is used to repeat the procedures of fatigue analysis for a significant number of times so that the fatigue lives can be estimated taking into account the annual growth of traffic. A flat rate of increase is assumed for different types of vehicles. Combining with statistical distribution of number of cycles, the statistical distributions of fatigue lives for different critical locations can be obtained.

Flowchart of probabilistic approach of fatigue life assessment.
Results and discussions
For comparison, the measurements at the linear strain gauges of ‘SS-GLE-04’ and ‘SS-GLW-04’ are studied, which are deployed to record the strains at the lower flange of the two outer girders, that is, the first and the fourth main girders. In interpreting the labels, ‘E’ means the east side where the strain gauge is deployed on the first girder, ‘W’ means the west side where the strain gauge is deployed on the fourth girder and ‘L’ is the location as shown in Figure 1. The strain time-histories on 22 November 2007 are as shown in Figure 11. Using the rain-flow counting method, the fatigue lives obtained from data of strain gauges are 490 and 286 years, respectively.

Data segments of strain time-histories of (a) ‘SS-GLE-04’ and (b) ‘SS-GLW-04’ on 22 November 2007.
According to the process described above, the probability density curves of fatigue life without considering diurnal variations at locations ‘SS-GLE-04’ and ‘SS-GLW-04’ are presented in Figures 12 and 13 which roughly follow the normal distribution. With the diurnal variations considered, the fatigue lives at the same locations still follow the normal distribution as shown in Figures 14 and 15, but the values may be different. As shown in Figures 16 and 17, the fatigue lives considering the mean stress effect are almost one-third of those without considering the mean stress effect. The mean values of fatigue lives calculated using the probabilistic approach here are very close to those obtained from strain gauge measurements, which indicates reasonable accuracy of this method. Finally, the fatigue life variations in 5 years at locations of ‘SS-GLE-04’ and ‘SS-GLW-04’ considering mean stress effect and the hypothetical annual growth rate of traffic of 5% are shown in Figures 18 and 19, respectively. However, the results should be interpreted cautiously as the growth rate may be limited by the traffic flow capacity.

Probability density curve of fatigue life at location ‘SS-GLE-04’ without considering diurnal variations.

Probability density curve of fatigue life at location ‘SS-GLW-04’ without considering diurnal variations.

Probability density curve of fatigue life at location ‘SS-GLE-04’ considering diurnal variations but without considering mean stress effect.

Probability density curve of fatigue life at location ‘SS-GLW-04’ considering diurnal variations but without considering mean stress effect.

Probability density curve of fatigue life at location ‘SS-GLE-04’ considering diurnal variations and mean stress effect.

Probability density curve of fatigue life at location ‘SS-GLW-04’ considering diurnal variations and mean stress effect.

Fatigue life variation at location of ‘SS-GLE-04’ considering annual growth of traffic of 5%.

Fatigue life variation at location of ‘SS-GLW-04’ considering annual growth of traffic of 5%.
Conclusion
A statistical traffic load model is developed based on WIM data. The proposed advanced traffic load model is able to incorporate the variation of daily traffic flow within different time intervals and the trend of the annual traffic flow. Moreover, the uncertain information of vehicle type, number of axles, axle weight, axle spacing, transverse position of vehicle, vehicle speed and heading direction can be accounted for properly. Combining with the statistical traffic load model, a probabilistic approach is proposed to evaluate the variation of fatigue lives of Ting Kau Bridge with time of service. The fatigue lives at the monitored sections of the bridge are all within their designed limit of 120 years.
Footnotes
Acknowledgements
The authors gratefully acknowledge the Highways Department of the Hong Kong Government for assistance received in producing this article as well as permission of its publication.
Declaration of conflicting interests
The author(s) declared the following potential conflicts of interest with respect to the research, authorship and/or publication of this article: Any opinions expressed or conclusions reached in this article are entirely of the authors.
Funding
The author(s) received no financial support for the research, authorship and/or publication of this article.
