1. Introduction
Fatigue of steel bridges has always been one of the major concerns of bridge engineers. An increasing number of steel bridges are suffering from fatigue cracks under the action of cyclic heavy trucks, which may threaten the safety of bridges. Therefore, it is urgent and significant to research the fatigue problem of steel bridges under overloaded truck traffic and take reasonable measures for vehicle weight limits.
The traditional idea for vehicle weight limits is checking and comparing the value of the stress induced by the overloaded truck with that induced by the design load [
1]. Many countries or regions have issued corresponding regulations to check whether a truck is permitted to cross a bridge. For example, the Bridge Formula B is widely applied in the United States to judge whether a heavy vehicle can cross a bridge [
2]. The Bridge Formula B is given as follows:
where
W denotes the permissible weight in newtons of any set of multi-axles;
N represents the number of axles in the selected axle set; and
B denotes the length between the two outermost axles in the selected axle set. However, some researchers criticized that the bridge damage caused by repeated truck loadings is not taken into account when the Bridge Formula B is adopted for permitting check [
3,
4]. Deng and Yan [
4] proposed a new weight-limiting method that can consider the effect of fatigue damage under the cyclic loadings of trucks based on the theory of the S-N curve and Palmgren–Miner’s rule. Furthermore, Wang et al. [
5] established a random traffic load model to consider the simultaneous existence of multiple trucks, and proposed a new method for determining the truck weight limit according to fatigue reliability analyses. However, these studies were performed on the basis of the theory of the S-N curve and Palmgren–Miner’s rule, in which the effect of original cracks on the fatigue life of steel bridges can hardly be taken into account [
6,
7].
The linear elastic fracture mechanics (LEFM) is an important branch of fracture mechanics, and its feasibility and reliability in fatigue analysis of steel bridges with cracks have been fully verified [
8]. The location and the size of structural cracks can be considered when the theory of LEFM is adopted for analyzing the remaining fatigue life of steel bridges [
6], which contributes to a more accurate fatigue analysis of in-service steel bridges with cracks. Due to the strong uncertainty and randomness during the process of crack growth, many scholars conducted probabilistic analyses on the fatigue problem of steel bridges based on the LEFM. Pipinato et al. [
9] used the probabilistic methods and LEFM to calculate the fatigue reliability of a highway steel girder bridge under the dual effects of earthquake loads and traffic loads. Guo and Chen [
10] investigated the growth of fatigue cracks of an old steel box girder bridge, in which fatigue reliability analyses on the basis of the LEFM were performed according to the long-term stress monitoring data. Leander and Al-Emrani [
11] proposed a probabilistic method for the fatigue evaluation of steel bridges on the basis of the LEFM, and investigated the impact of different modeling options on the reliability index of details. Although a large number of studies have been conducted to analyze crack growth in steel bridges based on the LEFM, few studies have been conducted to investigate the effects of traffic load conditions on the fatigue reliability of steel bridges with cracks based on the theory of LEFM, especially for the case of implementing truck weight limits. Therefore, it is necessary to provide an efficient method that can quickly analyze the fatigue reliability of steel bridges with cracks under different traffic load conditions, and provide valuable suggestions for setting reasonable truck load limits accordingly.
In recent years, many new artificial intelligence-based methods have been developed and applied to bridge health monitoring and fatigue analysis [
12,
13,
14,
15]. Guo et al. [
12] used Kohonen neural network and long short-term memory neural network to evaluate the bridge health state. Huynh [
13] proposed a new autonomous vision-based bolt-looseness detection method for splice bolted connections based on the Faster RCNN. Lu et al. [
14] proposed a machine learning-based framework to calculate the fatigue reliability of steel bridge decks. Yan et al. [
15] proposed a probabilistic machine learning method for fatigue analysis of steel bridges under the action of overloaded trucks.
In the present study, a new framework is proposed for fatigue reliability analysis of steel bridges with cracks, which can be used to quickly calculate the fatigue failure probability of steel bridges under truck weight limits according to different traffic load parameters. A typical 30.4-m-span steel-concrete composite girder bridge with a crack at the bottom flange of the girder was utilized as an instance to demonstrate the proposed framework. According to the transportation survey data from Wisconsin of the United States, a numerical model of random traffic load under load-limiting conditions is first established. Then, the limit state function of fatigue failure is obtained using the theory of LEFM when taking into account the crack detection error and the uncertainty process of crack growth. Finally, the ANN-MCS method is introduced to analyze fatigue reliability of the crack under consideration. Effects of four key parameters, including the average daily truck traffic, the gross vehicle weight limit, the violation rate, and the crack size, on the time-varying fatigue reliability of the steel bridge are quantitatively analyzed. Based on the analysis, a new method was proposed to set the gross vehicle weight limit for in-service steel bridges with cracks.
4. Method of LEFM
The theory of LEFM was demonstrated to perform well for fatigue analysis of steel bridges with cracks [
8]. The Paris Law, which is one of the most representative rule based on the theory of LEFM, describes the relationship between the rate of growth of a crack and stress-intensity factor (
K) induced by the cyclic loads, as show in Equation (2).
where
a denotes the crack size;
N denotes the number of cycles;
C and
m are constants related to material properties, and
m = 3 for carbon steels;
K denotes the stress-intensity factor, which describes the magnitude of the stress and strain values near the crack tip; and Δ
K is the stress-intensity factor range. The crack size of the fatigue detail under consideration can be calculated by solving the differential equation.
The stress intensity factor
K is influenced by many factors, including the crack shape, the crack size, and the component size. According to the research by Newman and Raju [
22], the
K for a semi-elliptical surface crack under the action of the remote tensile stress can be calculated to be:
where
σ denotes the remote uniform-tension stress. A reasonable estimation of
Q can be obtained in Equation (4) when
a/
c ≤ 1.
The stress-intensity boundary-correction factor in Equation (3) for surface cracks was taken as
where
In these equations,
a and
c represent the depth and the half-length of the surface crack;
ϕ represents the parametric angle of the ellipse, adopting
π/2 for the point at the maximum depth; and
w and
b represent the thickness and the half-width of the bottom flange, taking a value of 22.2 mm and 177.8 mm, respectively. The aspect ratio of the crack (
a/
c) was assumed to be 0.62 in the present study [
11]. It is worth noting that, as the purpose of the study mainly focuses on investigating the growth of the crack along the thickness direction of the bottom plate, the crack size in the present study denotes the crack depth of
a shown in
Figure 4.
A set of stress cycles were extracted from the stress history caused by traffic loads using the rainflow-counting algorithm. As it is complicated to use a set of stress ranges with variable amplitudes to calculate the crack growth rate, the stress ranges were equivalently converted into a constant stress range (
Sre) to simplify the calculation process [
7,
8]. The equivalent stress range (
Sre) is computed as:
where
is the frequency of the
ith stress range
Sr,i. Therefore, according to Equation (3), the stress-intensity factor range was calculated as:
7. Results and Discussion
Monte Carlo Simulation (MCS) is a commonly used method for structural reliability analysis, whose reliability can be guaranteed if the amount of random samples is sufficient [
28]. A sampling number of 10
6 for MCS, which is sufficient for reliability analysis, is adopted in this study [
29]. The target value of the reliability index is commonly taken as 2.0 to 3.5 for steel bridges [
30], and a reliability index of 2.0 is equivalent to
Pf = 2.3%. In the present research, the target probability of fatigue failure was set to 2.3% [
14,
15]. Through trail calculations, it was found that there was a huge gap in the computational time between the conventional Monte Carlo simulation and the ANN-MCS method for reliability analysis. For example, assuming
t = 10 and ADTT = 3000, it takes about 8.3 h and 0.8 s to perform a reliability analysis by the conventional Monte Carlo simulation and the ANN-MCS method, respectively. Therefore, applying the ANN-MCS method to fatigue reliability analysis has a significant advantage in computational efficiency while ensuring accuracy. In the subsequent part, the effects of four key parameters, namely, the ADTT, GVWL, VR, and
adtc, on the failure probability of the considered fatigue detail are investigated.
7.1. Effects of the Average Number of Daily Truck Traffic (ADTT)
The failure probability (
Pf) of the fatigue detail under consideration was calculated under different ADTTs as the remaining service time (
t) increases, as shown in
Figure 17, in which the ADTTs ranges from 1500 to 5000 while the values of GVWL, VR, PFL, and
adtc were adopted to be 600 kN, 0.1, 0.15, and 8 mm, respectively. It is seen from
Figure 17 that the time required for the failure probability of the fatigue detail increasing to over 2.3% reduces significantly as the ADTT increases. Specifically, the required time decreases from over 20 years to around 7 years as the ADTT increases from 1500 to 5000.
7.2. Effects of the Gross Vehicle Weight Limit (GVWL)
The failure probability (
Pf) of the fatigue detail under consideration was calculated under different GVWLs as the remaining service time increases, as shown in
Figure 18, in which the GVWL ranges from 300 kN to 500 kN while the values of ADTT, VR, PFL, and
adtc were adopted to be 3000, 0.1, 0.15, and 8 mm, respectively. It should be noted that no load limit was also considered and illustrated in
Figure 18 for comparison. It can be observed from
Figure 18 that the implementation of truck weight limits can effectively extend the remaining service time of the bridge. Specifically, even if the GVW is limited to 600 kN, the remaining service time with a 2.3% failure probability can also be effectively increased from 9 years to 12 years compared to no weight limit. Furthermore, if the GVWL is set to 350 kN, the service time can be significantly extended to 18 years. In addition, it is interesting to find that when the GVWL is set to 450, 500, and 600 kN, respectively, there is not much difference in the remaining service time, which may be due to the fact that the proportion of trucks with the GVW ranging from 450 to 600 kN is small.
7.3. Effects of the Violation Rate (VR)
With the GVWL preset to 600 kN and 300 kN, the effects of VR on the fatigue failure probability of the considered fatigue detail are shown in
Figure 19 and
Figure 20, respectively. It is found from
Figure 19 and
Figure 20 that the remaining service time with a failure probability of 2.3% sharply reduces from 27 years to 14.5 years when the VR increases from 0.1 to 0.5 if the GVWL is taken a value of 300 kN. However, the remaining service time with a failure probability of 2.3% only decreases from 12 years to 10.3 years when the VR increases from 0.1 to 0.5 if the GVWL is adopted to be 600 kN. This indicates that the role played by the VR is closely related to the value of GVWL, and high VR values lead to significant increases in the fatigue failure probability when the GVWL is set low.
7.4. Effects of the Detected Crack Size (adtc)
In the present research, the actual crack size was assumed to follow a lognormal distribution with
adtc as the mean value and 0.25 as the CoV. When ADTT = 3000, GVWL = 600 kN, VR = 0.1, and PFL = 0.15, the failure probability was calculated under different
adtc values, as exhibited in
Figure 21. It is found from
Figure 21 that the detected crack size
adtc has a considerable impact on the failure probability of the bridge. If the
adtc is less than or equal to 2 mm, the expected remaining service time with a 2.3% failure probability is greater than 50 years, which indicates that the fatigue failure is unlikely to happen in the short term. However, if the
adtc exceeds 7 mm, the estimated remaining service time with a 2.3% failure probability is less than 16 years. Furthermore, while the
adtc is greater than or equal to 9 mm, the predicted remaining service life with a 2.3% failure probability is less than 10 years, and corresponding measures should be taken to control the increase in the crack size, such as rehabilitation, traffic control, and truck weight limit enforcement.
7.5. Truck Weight Limits
A method of setting the gross vehicle weight limit for steel bridges with cracks is proposed in this section. When the target failure probability of fatigue detail is set to 2.3%, the recommended value of the GVWL can be determined based on the four parameters, including the average daily truck traffic (ADTT), violation rate (VR), detected crack size (
adtc), and remaining service time (
t). The relationship between the remaining fatigue life (
t) and the parameters under consideration was obtained by the ANN-MCS method. For example, assuming VR = 0.1, ADTT = 3000, PFL = 0.15, and
Pf = 2.3%, the relationship between
t and
adtc under different values of GVWL are presented as the solid lines in
Figure 22. Based on
Figure 22, a reasonable GVWL can be determined based on the detected crack size and the required remaining fatigue life. First, find the detected crack size from the abscissa, such as 8.5 mm, and make a vertical line marked as Line①, as shown by the dash-dotted line in
Figure 22. Then, find the expected remaining service time (
t) from the ordinate and draw a horizontal line, and then find the intersection point with the Line① made in the first step. Finally, find the nearest solid line right above the intersection and the GVWL value corresponding to this line is the accepted limit for gross vehicle weight.
For example, if adtc is 8.5 mm and the expected remaining service time (t) is 20 years, draw a horizontal Line②, and it can be seen that the solid line representing GVWL = 300 kN is just above the intersection of Line② and Line①. Therefore, the gross vehicle weight limit could set to 300 kN. Similarly, if the expected remaining service time is 15 years, draw the horizontal Line③ and it can be found that the solid line denoting GVWL = 350 kN is just above the intersection of Line③ and Line①. Therefore, the accepted gross vehicle weight limit, in this case, is 350 kN.
8. Summary and Conclusions
Very few previous studies investigated the influences of the gross vehicle weight limit (GVWL) on the fatigue failure probability of in-service steel bridges with cracks. In this study, a new framework was proposed for fatigue reliability analysis of steel bridges with cracks, in which the computational efficiency in calculating the fatigue failure probability was significantly improved through introducing the ANN. The framework proposed in the present study can help predict the fatigue failure probability of steel bridges with cracks based on four key parameters, including the average daily truck traffic, gross vehicle weight limit, violation rate, and crack size. In addition, with the goal of controlling fatigue failure probability to be less than the target value of 2.3%, a new method was proposed to set the gross vehicle weight limit for in-service steel bridges with cracks based on the detected crack size and the desired remaining service time.
A random traffic load model under weight limits was established based on the transportation survey data from Wisconsin of the United States. A typical 30.4-m-span steel-concrete composite girder bridge with a crack at the bottom flange of the girder was utilized as an example to illustrate the feasibility of the suggested framework for fatigue reliability analysis. Results show that when using the Monte Carlo method for fatigue reliability analysis of steel bridges with cracks, the conventional process of calculating the limit state function is time-consuming. As a substitution, the proposed ANN-MCS method shows great performance in terms of high accuracy, high efficiency, and good robustness in calculating the limit state function. Through parametric analysis, it is observed that both the flow-controlling measures and gross vehicle weight limits are effective for prolonging the remaining fatigue life of steel bridges with cracks. Strict implementation of the gross vehicle weight limit is essential because high violation rates can lead to significant increases in the probability of fatigue failure. In addition, the probability of fatigue failure is particularly sensitive to the detected crack size and when the size is greater than 9 mm, the predicted remaining service life with a failure probability of 2.3% is less than 10 years and measures should be taken in time to control the increase in the crack size.
It should be noted that only the vehicle-induced fatigue damage was taken into consideration in the present study, and the influences of many other factors, including corrosion and seismic loads, on the fatigue damage were ignored and ought to be given attention in future research.