Next Article in Journal
Direction of Arrival Estimation for MIMO Radar via Unitary Nuclear Norm Minimization
Previous Article in Journal
Improving Passive Time Reversal Underwater Acoustic Communications Using Subarray Processing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Probabilistic Fatigue Life Updating for Railway Bridges Based on Local Inspection and Repair

1
School of Urban and Environmental Engineering, Ulsan National Institute of Science and Technology (UNIST), Ulsan 44919, Korea
2
Fire Research Institute, Korea Institute of Civil Engineering and Building Technology, Hwaseong 18544, Korea
3
Department of Transportation and Logistics Engineering, Hanyang University ERICA Campus, Ansan 15588, Korea
4
Department of Mechanical and Aerospace Engineering, Trine University, Angola, IN 46703, USA
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(4), 936; https://doi.org/10.3390/s17040936
Submission received: 21 February 2017 / Revised: 12 April 2017 / Accepted: 19 April 2017 / Published: 24 April 2017
(This article belongs to the Section Physical Sensors)

Abstract

:
Railway bridges are exposed to repeated train loads, which may cause fatigue failure. As critical links in a transportation network, railway bridges are expected to survive for a target period of time, but sometimes they fail earlier than expected. To guarantee the target bridge life, bridge maintenance activities such as local inspection and repair should be undertaken properly. However, this is a challenging task because there are various sources of uncertainty associated with aging bridges, train loads, environmental conditions, and maintenance work. Therefore, to perform optimal risk-based maintenance of railway bridges, it is essential to estimate the probabilistic fatigue life of a railway bridge and update the life information based on the results of local inspections and repair. Recently, a system reliability approach was proposed to evaluate the fatigue failure risk of structural systems and update the prior risk information in various inspection scenarios. However, this approach can handle only a constant-amplitude load and has limitations in considering a cyclic load with varying amplitude levels, which is the major loading pattern generated by train traffic. In addition, it is not feasible to update the prior risk information after bridges are repaired. In this research, the system reliability approach is further developed so that it can handle a varying-amplitude load and update the system-level risk of fatigue failure for railway bridges after inspection and repair. The proposed method is applied to a numerical example of an in-service railway bridge, and the effects of inspection and repair on the probabilistic fatigue life are discussed.

1. Introduction

Railway bridges, which play an important role in the daily lives of numerous people, are known to be prone to the risk of fatigue failure. Many structural systems are subjected to the risk of fatigue-induced failure caused by repeated loading over their life cycle, which is a particularly critical safety concern for railway bridges. Several methods have been proposed for the fatigue damage mitigation of structures such as vibration control (e.g., tuned mass dampers) [1], local stress reduction (e.g., crack-arrest holes) [2], and compressive stress introduction (e.g., peening) [3]. To make appropriate decisions about effective bridge maintenance and repair, it is required to predict the fatigue life accurately. However, fatigue life prediction is a challenging task because such a prediction for a bridge requires considering results of inspection and repair, which are affected by various sources of uncertainty, including material properties, anticipated vehicle loads, and environmental conditions.
Two approaches are commonly employed for the fatigue damage evaluation and life prediction of bridge structures. The first approach is the traditional stress–life (S–N) curve method, in which the relationship between the constant-amplitude stress range, S, and the number of cycles to failure, N, is determined by performing appropriate fatigue experiments and described by an S–N curve. In this approach, the Palmgren–Miner linear damage hypothesis, also called Miner’s rule [4], extends this approach to varying-amplitude loadings. The second method is the fracture mechanics approach, which is dominantly dedicated to exploring the features and disciplines of crack initiation and growth while considering the stress field at the crack tip.
In general, the two approaches are applied sequentially, with the S–N curve method being used at the bridge design stage or for the preliminary evaluation of fatigue life and the fracture mechanics approach being used for more refined crack-based assessment of the remaining fatigue life or for effective decision-making on inspection and maintenance strategies [5]. The S–N curve is obtained entirely from experiments and has limitations in estimating the probabilistic bridge life based on the current condition of a bridge. Meanwhile, the fracture mechanics approach is usually applied to predict the propagation life from an initial crack or defect. The method based on linear elastic fracture mechanics (LEFM) relates the growth of an initial crack of size to the number of fatigue cycles. Many investigations have been performed on bridge fatigue condition assessment using the fracture mechanics approach. Fisher [6] illustrated more than 25 case studies on fatigue crack phenomena in steel bridges using the fracture mechanics approach and other theories. Agerskov and Nielsen [7] investigated the fatigue damage accumulation in steel bridges under random loadings and determined the fatigue life of welded joints in steel highway bridges by carrying out fracture mechanics analysis. By using an LEFM model to predict crack growth, MacDougall et al. [8] quantified the differences between the fatigue lives of a short-span bridge and a medium-span bridge under successive passages of either a steel-sprung vehicle or an air-sprung vehicle. Using experimental data obtained from the structural components on the Kinuura Bridge and theoretical predictions based on LEFM, Xiao et al. [9] found that the lack of penetration zones of 2–3 mm resulted in the low fatigue strength of butt-welded joints.
Structural health monitoring (SHM) systems have often been used to probabilistically evaluate the fatigue life of a bridge based on the measurement of its current condition [10]. LEFM approaches are more suitable for fatigue life evaluation than S–N curve approaches because LEFM can describe fatigue crack propagation quantitatively in terms of the crack length, whereas the S–N curve can only identify whether the structure will fail or survive. Recently, long-term SHM of bridges has been one of the major topics of interest for researchers and engineers in the fields of civil, mechanical, materials, and computer science engineering [11]. The design and implementation of such an SHM system involve the integration of analytical skills and instrumentation technologies with the knowledge and experience of bridge design, construction, inspection, maintenance, and management. An online SHM system can provide reliable information pertaining to the integrity, durability, and reliability of bridges. The information can then be incorporated into a bridge management and maintenance system for optimizing maintenance actions and improving design standards, specifications, codes, and guidelines. The SHM system is, in fact, an augmentation but not a substitute for systems used in current practice in bridge maintenance and management, not only through the use of advanced technologies in sensing, data acquisition, computing, communication, and data and information management but also through the effective integration of these technologies into an intelligent system. In addition, online SHM strategies have been applied to detect structural damage by system identification and model updating [12,13,14]. An accurate estimation of the actual situation pertaining to the fatigue life and remaining life of critical components is an important task of the SHM system; this task can be accomplished by using the continuously measured data of dynamic strain from the long-term SHM system. The fatigue performance of steel bridges depends on a number of factors such as material characteristics, stress history, and environmental conditions, and all these factors exhibit uncertainty and randomness during the service life of the bridge. On the other hand, when field measurement data are used for fatigue condition assessment, the uncertainties related to the data and the inaccuracies resulting from the use of data-processing techniques persist and are hardly avoidable. In view of these facts, it is more appropriate to conduct fatigue life assessment in a probabilistic manner rather than in a deterministic manner.
However, the results of bridge inspection and repair, which is one of the important SHM datum, are not used in the abovementioned studies. Lee and Song [15] proposed a new system reliability approach for updating the structural fatigue life by using formulations for fatigue life evaluation. In their research, however, the proposed method was applied to a structure where the loading condition was simplified as a constant amplitude. In addition, it is not feasible to update the prior risk information after bridges are repaired.
In this paper, a new method is proposed for the effective fatigue life evaluation of railway bridges by addressing the following challenges: (i) the formulations proposed by Lee and Song [15] are further developed so that the proposed method can handle varying-amplitude loads; and (ii) the proposed method enables updating the system-level risk of fatigue failure for railway bridges after inspection and repair. To demonstrate the method, it is applied to a numerical example of an in-service railway bridge. SHM data for the railway bridge subjected to actual train loading are introduced and analyzed by using rainflow counting, and the fatigue life of the bridge is evaluated and updated based on various results of inspection and repair.

2. Proposed Methodology

2.1. Limit-State Function Formulations for Fatigue Failure

In a structural reliability problem, the failure probability calculation requires constructing the so-called limit-state function, which represents the failure event of interest as an analytical function of random variables and deterministic parameters [10,15]. As mentioned in Section 1, many methods for the fatigue reliability analysis of bridges have been developed based on the S–N or LEFM approach. A widely used LEFM-based numerical model for crack propagation is the Paris equation [16]. Based on this equation, Lee and Cho [10] and Lee and Song [15] derived a series of formulations for estimating the structural risk of fatigue-induced failure. However, the formulations are based on constant-amplitude loads, and they cannot be applied to a railway bridge because they cannot consider cyclic loads with varying amplitudes, which is a realistic pattern generated by train traffic. Therefore, in this study, the limit-state functions proposed by Lee and Cho [10] and Lee and Song [15] are further developed for evaluating the fatigue life of a railway bridge with varying-amplitude loads.
First, consider the following Paris equation [16], which is a widely-known model of crack propagation:
d a d N = C ( Δ K ) m
where a denotes the crack length, N denotes the number of loading cycles, C and m represent the material properties, and ΔK is the range of the stress intensity factor. This range can be evaluated by using Newman’s approximation [17]:
Δ K = S Y ( a ) π a
where ΔS denotes the stress amplitude and Y(a) is the geometry function, which is a dimensionless parameter that depends on the geometry around the crack. The following equation can be obtained by substituting Equation (2) into Equation (1):
1 [ Y ( a ) π a ] m d a = C Δ S m d N
If the load amplitude is a constant, Equation (3) can be developed to an equation for the fatigue life evaluation of a structure, as described by Lee and Song [15]. However, when a train passes over a bridge, the stress experienced by a structural member is found to be fluctuating, as shown in Figure 1, and this makes it impossible to use the formulations devised for a constant amplitude.
The rainflow cycle counting algorithm [18] is widely used to handle such a varying-amplitude load and assess the fatigue life of structures, especially in mechanical engineering. Usually, the algorithm extracts cycles from the load, stress, or strain history obtained from measurement or simulation results. From the cycle counting, several cycles and half-cycles with different amplitude and mean values are obtained. This counting algorithm can be used with the Paris equation, which enables it to compute the expected fatigue life under random loading conditions in a theoretical manner [18].
The stress amplitude ΔSi and the corresponding number of cycles ni are used for the integration of Equation (3) from the initial condition to the current time point, thereby providing the following relationship between the current crack length and the time duration as follows:
a 0 a 1 [ Y ( a ) π a ] m d a = C ν 0 T i = 1 N a m p [ n i ( Δ S i ) m ]
where a0 is the initial crack length, ν0 is the frequency of train loading, T is the time duration, and Namp is the total number of stress amplitude values. If a crack failure is assumed to occur when the crack length exceeds the critical crack length ac, the time required for crack growth from a0 to ac, T, is expressed by:
T = 1 C ν 0 i = 1 N a m p [ n i ( Δ S i ) m ] a 0 a c 1 [ Y ( a ) π a ] m d a
The limit-state function for the failure of the j-th structural member or location (“member” is used hereafter) within a given time interval [0, Ts] is expressed by:
g ( X ) = T j 0 T s = 1 C ν 0 i = 1 N a m p [ n i ( Δ S i ) m ] a 0 a c 1 [ Y ( a ) π a ] m d a T s
where X denotes the vector of random variables and Tj0 is the time required for the fatigue failure of the j-th structural member. In structural reliability problems, g(X) ≤ 0 typically indicates the occurrence of a failure event.
Although several studies [10,15,19,20] for the probability estimation of fatigue failure events were conducted with the Paris equation in Equation (1), other crack propagation models such as the Walker equation [21] and the Forman equation [22] can also be introduced to model crack growth more accurately. Various crack growth models have been suggested to improve the Paris equation, and any model can replace the Paris equation in the proposed method, if it satisfies the following conditions: (1) the time duration of the crack growth can be analytically expressed in terms of random variables and deterministic variables, as done in Equations (1)–(5); and (2) the statistical properties (e.g., distribution type, mean, standard deviation, and correlation) of the random variables are available. The first condition can be met by several crack growth models. In this research, however, the Paris equation is selected because the statistical characteristics of the equation parameters were investigated sufficiently as discussed in the previous studies [10,15].
If a railway bridge is considered to be a structural system consisting of multiple structural members, the failure of the bridge system should be expressed by a Boolean function of the failure events of individual members; to obtain such a function, system reliability analysis is required. In this research, a bridge is considered to fail if any of its members fails because of fatigue, and the probability of this system event, Psys, is described as:
P s y s = P [ j = 1 N m e m ( T j 0 < T s ) ]
where Nmem is the total number of members in consideration.

2.2. Reliability Updating through Inspection and Repair Events

A new approach to update the probability of fatigue failure based on inspection results was recently proposed [15]. However, the formulations provided in the research were developed for aircraft and cannot be used to update the prior risk information after repairing events. In the present research, the previously proposed formulations are introduced and then further developed so that it can be applied to railway bridges after both of inspection and repairing events.
For an event of interest Ei, (e.g., a member failure in Equation (6) or a system failure in Equation (7)), its probability Pi can be updated through an inspection event as follows:
P i , u p = P ( E i | I E j )
where IEj denotes the inspection event for the j-th member. This formulation can be extended to utilize multiple available inspection results to update the probability Pi, i.e.,:
P i , u p = P ( E i | I E j I E k I E l )
Inspection results are classified into two types of events—equality and inequality events—depending on whether a crack is detected (and measured) or not. Considering that inspections are generally made at multiple locations, three possible combinations of inspection results exist: inequality, equality, and mixed cases [15,23].
First, an inequality event indicates that no crack is detected through inspection at a certain single location. There are two possible explanations as to why no crack is detected. First, a crack is too small to be detected. Second, although a relatively large crack actually exists on the inspected member, it may be missed because of human error or the limitations of the detecting device. In either case, this event can be formulated as:
I E j : g j , n o ( X ) = T I T j d ( X ) < 0 where T j d = 1 C ν 0 i = 1 N a m p [ n i ( Δ S i ) m ] a j 0 a d 1 [ Y ( a ) π a ] m d a
where Tjd denotes the time required for the crack to grow to a detectable crack size in the j-th member, ad is the detectable crack size during the inspection process, and TI denotes the time at which the member is inspected. The detectable crack size ad is related to a specific inspection method and is modeled as a random variable reflecting the actual probability of detection.
For a single inequality event, the conditional probability in Equation (8) can be calculated as:
P i , u p = P ( E i | I E j ) = P ( E i ( g j , n o ( X ) < 0 ) ) P ( g j , n o ( X ) < 0 )
Likewise, the conditional probability for multiple inequality events can be calculated as:
P i , u p = P ( E i | I E j I E k ) = P [ E i ( g j , n o ( X ) < 0 ) ( g k , n o ( X ) < 0 ) ] P [ ( g j , n o ( X ) < 0 ) ( g k , n o ( X ) < 0 ) ]
The updated probabilities can be obtained from Equations (11) and (12) by computing the probabilities in the numerator and the denominator by component and system reliability analysis.
Second, an equality event indicates that a crack is detected and measured. The event is formulated as:
I E j : g j , y e s ( X ) = T I T j m ( X ) = 0 where T j m = 1 C ν 0 i = 1 N a m p [ n i ( Δ S i ) m ] a j 0 a m + ε m 1 [ Y ( a ) π a ] m d a = 0
where Tjm denotes the time required for the crack to grow to the measured crack size in the j-th member, am is the measured crack size, and εm is the measuring error. In the equation, it should be noted that Tjm is equal to the inspection time TI, which makes the formulation different from the limit-state formulations in Equation (10). The conditional probability for the equality event is calculated as:
P i , u p = P ( E i | I E j ) = P ( E i ( g j , y e s ( X ) = 0 ) ) P ( g j , y e s ( X ) = 0 )
Compared to the terms in Equation (11), the probability terms in Equation (14) are more challenging to compute using structural reliability methods because the probabilities of both the numerator and the denominator are zero [24]. Through an analytical derivation presented in a previous study [15], Equation (14) is transformed to:
P i , u p = [ θ P [ E i ( g j , y e s + θ 0 ) ] θ P [ ( g j , y e s + θ 0 ) ] ] θ = 0
The updated probability in the equation can be calculated via numerical differentiation. Likewise, the conditional probability for multiple equality events is formulated as follows [25]:
P i , u p = P ( E i | I E j I E k ) = [ n θ 1 θ n P [ E i ( g j , y e s + θ 1 0 ) ( g k , y e s + θ n 0 ) ] n θ 1 θ n P [ ( g j , y e s + θ 1 0 ) ( g k , y e s + θ n 0 ) ] ] θ 1 = θ n = 0
where n is the number of observed equality events. The updated probability Pi,up in Equation (16) can also be calculated using the n-th-order numerical differentiation. However, the equation needs to be used with caution because such a high-order numerical differentiation can produce significant error unless the probability calculations in the numerator and the denominator are extremely accurate. Therefore, the equality case including multiple equality events is not discussed in this research.
Lastly, with regard to mixed cases, the simplest case involves a single inequality event and a single equality event. The updated probability is formulated as:
P i , u p = P [ E i ( g j , n o ( X ) < 0 ) ( g j , y e s ( X ) = 0 ) ] P [ ( g j , n o ( X ) < 0 ) ( g j , y e s ( X ) = 0 ) ]
This equation can be transformed to:
P i , u p = [ θ P [ E i ( g j , n o < 0 ) ( g j , y e s + θ 0 ) ] θ P [ ( g j , n o < 0 ) ( g j , y e s + θ 0 ) ] ] θ = 0
Finally, this formulation can be generalized for a mixed case involving multiple inequality and equality inspection events as:
P i , u p = { n θ 1 θ n P [ E i ( g j , n o < 0 ) ( g k , n o < 0 ) ( g l , y e s + θ 1 0 ) ( g l , y e s + θ n 0 ) ] / n θ 1 θ n P [ ( g j , n o < 0 ) ( g k , n o < 0 ) ( g l , y e s + θ 1 0 ) ( g l , y e s + θ n 0 ) } θ 1 = θ n = 0
It should be noted that Equation (19) is the combination of Equations (12) and (16). Because of the aforementioned difficulty in carrying out high-order numerical differentiation, the mixed case involving multiple equality events is not handled in this research.
Overall, system reliability updating after repair is similar to that through inspection. However, when a crack is detected, measured, and repaired, the crack size in the repaired member changes; therefore, structural analysis needs to be performed with a new initial crack length for the repaired member. When the repair is made at Tr, according to Equation (6), the limit-state function for the failure of the repaired member j within an inspection cycle [0, Ts] is described as:
g ( X ) = T j r + T r T s T j r = 1 C ν 0 i = 1 N a m p [ n i ( Δ S i ) m ] a j r a j c 1 [ Y ( a ) π a ] m d a
where ajr is the initial crack length in the repaired member. ajr depends on the quality of the repair methods, e.g., grinding and welding.
In summary, the fatigue failure events of a structural member and structural system are expressed by analytical functions of random variables and deterministic parameters shown in Equations (6) and (7), respectively. These equations are derived through the analytical derivation from Equation (1). Likewise, the events of no-crack detection, crack detection and size measurement, and repair are expressed by Equations (10), (13) and (20), respectively. With these equations, the original probability of a fatigue failure event which can be estimated from Equations (6) and (7) can be updated based on the conditional probability, as shown in Equations (9), (16) and (19).
A reliability analysis needs to be performed to calculate the probability of a fatigue failure event with these formulations, and numerous reliability analysis methods have been developed so far [26]. Among them, the first-order reliability method (FORM) [27] and the second-order reliability method (SORM) [28] are used to perform component reliability analysis. In this research, FERUM (Finite Element Reliability Using MATLAB) [29], which is a reliability software package developed at the University of California, Berkeley, was modified and then introduced to perform reliability analysis. In addition, for a structural system consisting of multiple structural members, a failure event needs to be expressed by a system event which requires system reliability analysis (SRA) [30,31]. In this study, it is assumed that the railway bridge system fails if any of its structural members fail, and the Matrix-based System Reliability (MSR) method [31] is employed to calculate the system probability.

3. Numerical Example

3.1. Example of Bridge Considered in This Study: Calumet Bridge

The proposed method of fatigue life updating is tested by applying it to a real bridge named Calumet Bridge located in Illinois, USA. The bridge is located over the Little Calumet River near Chicago, IL (Figure 2).
The following four bridges are shown in Figure 2: the test bridge is the second to left truss bridge (the bridge marked with an arrow), the west-most bridge is used for the Metra lines, another truss bridge carrying freight and passenger trains is located to the east of the test bridge, and the east-most bridge is closed for traffic. In Kim et al. and Moreu et al. [32,33], the structural condition of the Calumet Bridge was assessed based on in-service data obtained from an SHM system.
The bridge is approximately 95 m long, 21 m tall, and 10 m wide and is made of A36 American Society for Testing and Materials (ASTM) steel. The bridge opened for service in 1971, with an expected life of 100 years. The test bridge carries two tracks—the CN 1 track on the west side and the CN 2 track on the east side—both of which are open for south- and north-bound freight and passenger trains. Approximately ten freight trains and six passenger (Amtrak) trains run on either the CN 1 track or the CN 2 track on a daily basis. The following aspects of the bridge make it suitable and unique for use as the test bridge: (i) The test bridge is made of steel, which is the most common type of bridge construction material used in the U.S., but vulnerable to fatigue; (ii) Only the train loads are applied to the test bridge, whereas adjacent bridges are open to car traffic as well; (iii) A large amount of traffic comprising various types of trains on the daily basis allows rich types of excitation on the bridge. This bridge handles a large amount of traffic per day (approximately 20 trains, including both freight and Amtrak traffic).
Wireless smart sensor (WSS) networks were installed for a year on the bridge for creating an efficient SHM system with various types of sensors [32]. The wireless Imote2 platform was used in the research [34]. This platform can stack a sensor board via two external connectors, thus allowing the network to have a flexible design with any type of sensor board. The sensor boards used were (i) a normal accelerometer board; (ii) a high-sensitivity accelerometer board; and (iii) a strain sensor (denoted as SHM-S sensor).
Using the measured data from the WSS networks, the authors in Kim et al. [35] developed and calibrated a model that can efficiently estimate structural behavior subject to a moving train. The model was validated by comparing the measured axial strain from the network with the estimated strain from the model. The model is adopted in this paper, and the maximum axial strains estimated from the model, while a test car in Kim et al. [32] ran at 53.11 km/h is used in the paper. More details regarding the bridge including section properties and the model can be found in [32,35].

3.2. Monitoring Data and Rainflow Counting

Calumet Bridge consists of over hundreds of members, and it is not feasible to monitor all the members. Therefore, in this research, 10 global members in the calibrated model that show the highest levels of axial stresses under a test train [32] are selected for monitoring. The 10 members in the decreasing order of their maximum stresses are Members 13, 27, 8, 6, 7, 5, 10, 9, 17 and 4, whose locations are shown in Figure 3. Table 1 lists the maximum stresses of the 10 members; the stress levels of Members 13 and 27 are much higher than those of the other members. In addition, Figure 4 shows the stresses of Members 13 and 27 under a passing train.
Next, the rainflow counting algorithm is applied to the stress responses of the 10 selected members. The rainflow cycle counting algorithm is widely used for the fatigue life assessment of machine components or structures under varying-amplitude loading. Usually, the algorithm extracts cycles from the load, stress, or strain history obtained from measurement or simulation. In this research, a code for rainflow counting developed by Nieslony [18] is used and applied to the simulated results using the model developed by Kim et al. [35], and the number of amplitudes and the number of cycles are then obtained. The data are used to calculate i = 1 N a m p [ n i ( Δ S i ) m ] in Equation (4).

3.3. Random Variables and Deterministic Parameters

In this research, the statistical information of random variables is determined through a comprehensive literature review [10,15,19,36,37,38], and the mean, coefficient of variation (c.o.v.), and distribution type of the random variables in this numerical example are summarized in Table 2.
Although each of the random variables in this numerical example is assumed to follow the normal, lognormal, or exponential distribution, other distribution types such as Gumbel and Weibull distributions can also be introduced if they are more appropriate to represent the uncertainty of random variable. This is possible because the proposed method employs FERUM which enables the reliability with a total of 16 different probability distributions including normal, lognormal, exponential, Gumbel, and Weibull distributions.
In this example, only C is considered as a random variable and m is assumed to be a constant (3.344), because it was shown through a preliminary analysis that considering m as another random variable had a negligible effect on the failure probability result despite of a large computational cost.
Statistical independence is assumed among all the random variables, except the following variables: (1) values of the Paris equation parameter C of the 10 members (correlation coefficient: 0.6); (2) initial crack lengths (a0) in the 10 members (correlation coefficient: 0.6); (3) initial crack lengths in the repaired members (ar); (4) detectable crack sizes (ad); and (5) crack measurement errors (εm). The correlation coefficients in these cases are not known, so they are initially assumed to be 0.6, which indicates that the 10 members are manufactured by the same manufacturer and that their material properties are thus highly correlated.
In addition, the following deterministic parameters are used: L-bracket width (b) = 650 mm, critical crack length (ac) = 12.7 or 19.05 mm, and average train traffic (assumed based on Moreu et al. [33]) = 20/day. For the geometry function Y(a) in Equation (2), the following function for I-beams is introduced based on Tada et al. [36]:
Y ( a ) = 1.122 0.231 ( a b ) + 10.550 ( a b ) 2 21.710 ( a b ) 3 + 30.382 ( a b ) 4

4. Analysis Results

4.1. Fatigue Life Evaluation

Using the proposed method, the failure probabilities of the 10 selected members are calculated at various periods of the service time. Because the failure probabilities turned out to be very small, reliability index, a useful probability indicator which can be calculated by the following equation, is introduced:
β = Φ 1 ( 1 P f )
where β is the reliability index, Φ−1(·) is the inverse function of the standard normal CDF (cumulative distribution function), and the Pf is the failure probability.
Figure 5 shows the reliability indices of the 10 selected members at various periods of the service time. Overall, the reliability indices of the bridge decrease as the service life increases, which means that the probability of failure increases as the use of the bridge increases.
The American Association of State Highway and Transportation Officials (AASHTO) Bridge Design Code recommends a target reliability index of 3.5 (i.e., a failure probability of 2.33 × 10−4) with a service life of 75 years for steel and prestressed concrete components [10]. By finding the intersections between the reliability index curves and the target reliability index (i.e., black lines in Figure 5), the fatigue lives of all the 10 members are estimated. Accordingly, the fatigue lives of Members 13 and 27 are evaluated as 75.4 and 78.2 years, respectively, whereas the fatigue lives of the other members are estimated to be longer than 100 years. As mentioned in Section 3.2, this is mainly because the stress levels of Members 13 and 27 are much higher than those of the other members. Therefore, the analysis results of only Members 13 and 27 are presented hereafter.
In addition to the analysis results of the structural members, the fatigue life of the bridge system is evaluated. In the evaluation, as mentioned in Section 2.2, the system failure event is assumed to occur when any of the 10 members fails. Furthermore, the correlation coefficient between the values of the Paris equation parameter C and between the initial crack lengths a0 of the 10 members is assumed to be 0.6, which accounts for the high statistical dependence due to the assumption that all the members are produced by the same manufacturer. In addition, the fatigue life of the bridge system is evaluated using a range of correlation coefficient values (i.e., 0.0, 0.2, 0.4, 0.6 and 0.8). This evaluation is carried out to investigate the effect of the correlation coefficient as a parametric study, and the analysis results are shown in Figure 6.
The correlation coefficient has a negligible impact on the variations in the system bridge life against fatigue. This is because the system failure event is dominated by the failure events of Members 13 and 27, even though the failure event of the bridge system is assumed to occur if any of the 10 members fails. Therefore, in this example, the correlation coefficient is set to be 0.6 for the following analyses. With this correlation coefficient and the target reliability index of 3.5, the fatigue life of this bridge system is evaluated as 70.7 years.

4.2. Updated Fatigue Life Based on Local Inspection and Repair

To test the proposed method, it is applied to hypothetical scenarios of local inspection and repair, which are listed in Table 3.

4.2.1. Inequality Cases

Figure 7 shows the reliability index results of Members 13 and 27 and the bridge system for the inequality cases. First, the updated probabilities of different scenarios are compared with each other to investigate the impact of various inspection conditions on the reliability updating, such as the number of inspections, crack-detecting resolution, inspection interval, and measured crack length.
To show the reliability index updates more clearly, Figure 8 shows how the reliability indices of Member 13, Member 27, and the bridge system are updated in Scenarios 1–4. In Scenario 1, the updated reliability indices decrease for the two members and the system because no crack is detected even after 50 years. The updated indices decrease further in Scenario 2 because no crack is detected even when a better crack-detecting device is used (with a smaller mean of the detectable crack size ad). In Scenario 3, no crack is detected even though an inspection is made at a later time (i.e., 75 years) than in Scenario 1, which further increases the reliability indices. In Scenario 4, it is assumed that a crack is not observed at any of the members. Because there are 10 members of interest, the inspection event is described as the intersection of 10 inequality events, as shown in Equation (12). Compared to Scenario 1, nine additional signs indicating the better safety of the bridge are obtained in Scenario 4, which further increases the reliability indices, as shown in Figure 8. From Figure 8, the fatigue lives of Member 13, Member 27, and the bridge system in Scenarios 1–4 can be evaluated, and the results are shown in Table 4.

4.2.2. Equality Cases

Figure 9 shows the reliability index results for the equality cases, and the updated indices in Scenarios 5–7 are shown in Figure 10. From Scenarios 5 to 7, the measured crack size increases from 0.1 to 1.0 mm. It is clearly seen that the longer the crack size measured, the more likely the members and the system are to fail. In addition, the fatigue lives of Member 13, Member 27, and the bridge system in Scenarios 5–7 are shown in Table 5.

4.2.3. Mixed Cases

Figure 11 shows the results obtained by the proposed method for the mixed cases (i.e., Scenarios 8 and 9), and Figure 12 shows how the reliability indices are updated in these scenarios. For comparison purposes, the updated reliability indices in Scenario 7 are also presented in the figure. First, it is clearly seen that the updated failure probability becomes considerably lower in Scenario 8 than in Scenario 7 even though the same crack size is observed at the same inspection time point. This is due to the additional inequality events (i.e., no more cracks being detected at the other locations) observed in the mixed case for Scenario 8. However, after Members 13 and 27 are repaired, the updated reliability indices are increased. Furthermore, the fatigue lives of Member 13, Member 27, and the bridge system are shown in Table 6.
Lastly, it is noteworthy that the reliability updating could be performed in only a few minutes in each scenario (using a general personal computer with 3.60 GHz CPU and 8.00 GB RAM), which demonstrates the computational efficiency of the proposed method. MCS is not a feasible option because it would require a huge number of samples to obtain converged results. As discussed by Lee and Song [38], it would require 3 × 105 samples and a similar number of FE analyses to achieve 5% c.o.v. in the probability result, which is practically impossible even though each analysis takes only a few minutes. However, the proposed method makes system reliability updating feasible and provides considerable savings in terms of computational time.

5. Conclusions

This study presents a new system reliability approach to evaluate the fatigue failure risk of railway bridges and update their prior risk information in various scenarios of inspection and repair. The proposed method can handle a varying-amplitude load and can update the system-level risk of fatigue failure for truss bridges after inspection and repair. To demonstrate the proposed method, it is applied to a numerical example of an actual railway bridge. The analysis results show that the proposed method facilitates the updating of reliability estimates for railway bridges subjected to the risk of fatigue-induced failures. In addition, the effects of various inspection and repair scenarios on the reliability updating are investigated in the numerical example, such as the number of inspections, crack-detecting resolution, inspection interval, inspection location, measured crack length, and repair. The investigation results confirm that the proposed method successfully provides component and system reliability updates under various inspection and repair scenarios.

Acknowledgments

This study was supported by a grant (17SCIP-B066018-05) from the Smart Civil Infrastructure Research Program funded by the Ministry of Land, Infrastructure and Transport (MOLIT) of the Korean government and the Korea Agency for Infrastructure Technology Advancement (KAIA).

Author Contributions

Young-Joo Lee developed the reliability analysis program and performed the probabilistic analysis; Robin E. Kim developed the structural model of Calumet Bridge by calibrating it based on structural health monitoring data obtained from WSS networks. Wonho Suh analyzed the effect of train loads on Calumet Bridge and helped to perform reliability analysis considering the effect. Kiwon Park contributed to the analytical derivation of formulations and the literature review.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Domaneschi, M.; Martinelli, L.; Perotti, F. Wind and earthquake protection of cable-supported bridges. Bridge Eng. 2015, 169, 157–171. [Google Scholar] [CrossRef]
  2. Makabe, C.; Murdani, A.; Kuniyoshi, K.; Irei, Y.; Saimoto, A. Crack-growth arrest by redirecting crack growth by drilling stop holes and inserting pins into them. Eng. Fail. Anal. 2009, 16, 475–483. [Google Scholar] [CrossRef]
  3. Branco, C.M.; Infante, V.; Baptista, R. Fatigue behaviour of welded joints with cracks, repaired by hammer peening. Fatigue Fract. Eng. Mater. Struct. 2004, 27, 785–798. [Google Scholar] [CrossRef]
  4. Miner, M.A. Cumulative damage in fatigue. J. Appl. Mech. 1945, 12, 159–164. [Google Scholar]
  5. Chryssanthopoulos, M.K.; Righiniotis, T.D. Fatigue reliability of welded steel structures. J. Constr. Steel Res. 2006, 62, 1199–1209. [Google Scholar] [CrossRef]
  6. Fisher, J.W. Fatigue and Fracture in Steel Bridges: Case Studies; John Wiley & Sons: New York, NY, USA, 1984. [Google Scholar]
  7. Agerskov, H.; Nielsen, J.A. Fatigue in steel highway bridges under random loading. ASCE J. Struct. Eng. 1999, 125, 152–162. [Google Scholar] [CrossRef]
  8. MacDougall, C.; Green, M.F.; Shillinglaw, S. Fatigue damage of steel bridges due to dynamic vehicle loads. J. Bridg. Eng. 2006, 11, 320–328. [Google Scholar] [CrossRef]
  9. Xiao, Z.-G.; Yamada, K.; Inoue, J.; Yamaguchi, K. Fatigue cracks in longitudinal ribs of steel orthotropic deck. Int. J. Fatigue 2006, 28, 409–416. [Google Scholar] [CrossRef]
  10. Lee, Y.J.; Cho, S. SHM-based probabilistic fatigue life prediction for bridges based on FE model updating. Sensors 2016, 16, 317. [Google Scholar] [CrossRef] [PubMed]
  11. Nagayama, T.; Sim, S.H.; Miyamori, Y.; Spencer, B.F. Issues in structural health monitoring employing smart sensors. Smart Struct. Syst. 2007, 3, 299–320. [Google Scholar] [CrossRef]
  12. Kurt, M.; Eriten, M.; McFarland, D.M.; Bergman, L.A.; Vakakis, A.F. Methodology for model updating of mechanical components with local nonlinearities. J. Sound Vib. 2015, 357, 331–348. [Google Scholar] [CrossRef]
  13. Kurt, M.; Moore, K.J.; Eriten, M.; McFarland, D.M.; Bergman, L.A.; Vakakis, A.F. Nonlinear model updating applied to the IMAC XXXII Round Robin benchmark system. Mech. Syst. Signal Process. 2017, 88, 111–122. [Google Scholar] [CrossRef]
  14. Chen, H.; Kurt, M.; Lee, Y.S.; McFarland, D.M.; Bergman, L.A.; Vakakis, A.F. Experimental system identification of the dynamics of a vibro-impact beam with a view towards structural health monitoring and damage detection. Mech. Syst. Signal Process. 2014, 46, 91–113. [Google Scholar] [CrossRef]
  15. Lee, Y.J.; Song, J. System reliability updating of fatigue-induced sequential failures. J. Struct. Eng. 2014, 140, 04013074. [Google Scholar] [CrossRef]
  16. Paris, P.C.; Erdogan, F. An effective approximation to evaluate multinormal integrals. Struct. Saf. 1963, 20, 51–67. [Google Scholar]
  17. Newman, J.C.; Raju, I.S. An empirical stress intensity factor equation for the surface crack. Eng. Fract. Mech. 1981, 15, 185–192. [Google Scholar] [CrossRef]
  18. Nieslony, A. Determination of fragments of multiaxial service loading strongly influencing the fatigue of machine components. Mech. Syst. Signal Process. 2009, 23, 2712–2721. [Google Scholar] [CrossRef]
  19. Moan, T.; Song, R. Implications of inspection updating on system fatigue reliability of offshore structures. J. Offshore Mech. Arct. Eng. 2000, 122, 173–180. [Google Scholar] [CrossRef]
  20. Karamchandani, A.; Dalane, J.I.; Bjerager, P. Systems reliability approach to fatigue of structures. J. Struct. Eng. 1992, 118, 684–700. [Google Scholar] [CrossRef]
  21. Forman, R.G.; Kearney, V.E.; Engle, R.M. Numerical Analysis of Crack Propagation in Cyclic Loaded Structures. J. Basic Eng. 1967, 89, 459–464. [Google Scholar] [CrossRef]
  22. Walker, K. The effect of stress ratio during crack propagation and fatigue for 2024-T3 and 7075-T6 aluminum. In Effects of Environment and Complex Load History on Fatigue Life; ASTM International: Philadelphia, PA, USA, 1970. [Google Scholar]
  23. Jiao, G.; Moan, T. Methods of reliability model updating through additional events. Struct. Saf. 1990, 9, 139–153. [Google Scholar] [CrossRef]
  24. Straub, D. Reliability updating with equality information. Probab. Eng. Mech. 2011, 26, 254–258. [Google Scholar] [CrossRef]
  25. Madsen, H.O. Model updating in reliability theory. In Proceedings of the 5th International Conference on Applications of Statistics and Probability in Civil Engineering (ICASP5), Vancouver, BC, Canada, 25–29 May 1987. [Google Scholar]
  26. Hohenbichler, M.; Rackwitz, R. First-order concepts in system reliability. Struct. Saf. 1983, 1, 177–188. [Google Scholar] [CrossRef]
  27. Melchers, R.E. Structural Reliability Analysis and Prediction; John Wiley & Son Ltd.: New York, NY, USA, 1999. [Google Scholar]
  28. Der Kiureghian, A. First- and second-order reliability methods. In Engineering Design Reliability Handbook; CRC Press: Boca Raton, FL, USA, 2005; Chapter 14. [Google Scholar]
  29. Sudret, B.; Der Kiureghian, A. Stochastic Finite Element Methods and Reliability: A State-of-the-Art Report; Department of Civil and Environmental Engineering, University of California: Berkeley, CA, USA, 2000. [Google Scholar]
  30. Song, J.; Der Kiureghian, A. Bounds on system reliability by linear programming. J. Eng. Mech. 2003, 129, 627–636. [Google Scholar] [CrossRef]
  31. Song, J.; Kang, W.H. System reliability and sensitivity under statistical dependence by matrix-based system reliability method. Struct. Saf. 2009, 31, 148–156. [Google Scholar] [CrossRef]
  32. Kim, R.E.; Moreu, F.; Spencer, B.F. System identification of an in-service railroad bridge using wireless smart sensors. Smart Struct. Syst. 2015, 15, 683–698. [Google Scholar] [CrossRef]
  33. Moreu, F.; Kim, R.E.; Spencer, B.F. Railroad bridge monitoring using wireless smart sensors. Struct. Control Health Monit. 2016, 24, e1863. [Google Scholar] [CrossRef]
  34. Crossbow Technology Inc. Imote2—High-Performance Wireless Sensor Network Node. Available online: http://wsn.cse.wustl.edu/images/e/e3/Imote2_Datasheet.pdf (accessed on 17 January 2017).
  35. Kim, R.E.; Moreu, F.; Spencer, B.F. Hybrid model for railroad bridge dynamics. J. Struct. Eng. 2016, 142, 04016066. [Google Scholar] [CrossRef]
  36. Tada, H.; Paris, P.C.; Irwin, G.R. The Stress Analysis of Cracks Handbook; ASME: New York, NY, USA, 2000. [Google Scholar]
  37. Yazdani, N.; Albrecht, P. Crack growth rates of structural steel in air and aqueous environments. Eng. Fract. Mech. 1989, 32, 997–1007. [Google Scholar] [CrossRef]
  38. Lee, Y.J.; Song, J. Finite-element-based system reliability analysis of fatigue-induced sequential failures. Reliab. Eng. Syst. Saf. 2012, 108, 131–141. [Google Scholar] [CrossRef]
Figure 1. Generic stress response of a railway bridge under passing train.
Figure 1. Generic stress response of a railway bridge under passing train.
Sensors 17 00936 g001
Figure 2. Test bridge (i.e., Calumet Bridge) marked as CN1 and CN2 [32,33].
Figure 2. Test bridge (i.e., Calumet Bridge) marked as CN1 and CN2 [32,33].
Sensors 17 00936 g002
Figure 3. Locations of 10 members with the maximum stresses.
Figure 3. Locations of 10 members with the maximum stresses.
Sensors 17 00936 g003
Figure 4. Stress history of Members 13 and 27.
Figure 4. Stress history of Members 13 and 27.
Sensors 17 00936 g004
Figure 5. Reliability indices of 10 selected members obtained from proposed method. (a) Members 13, 27, 8, 6 and 7; (b) Members 5, 10, 9, 17 and 4.
Figure 5. Reliability indices of 10 selected members obtained from proposed method. (a) Members 13, 27, 8, 6 and 7; (b) Members 5, 10, 9, 17 and 4.
Sensors 17 00936 g005
Figure 6. Reliability indices of bridge system with varying correlation coefficients.
Figure 6. Reliability indices of bridge system with varying correlation coefficients.
Sensors 17 00936 g006
Figure 7. Reliability updating results by the proposed method for inequality cases (Scenarios 1–4 in Table 3). (a) Scenario 1; (b) Scenario 2; (c) Scenario 3; (d) Scenario 4.
Figure 7. Reliability updating results by the proposed method for inequality cases (Scenarios 1–4 in Table 3). (a) Scenario 1; (b) Scenario 2; (c) Scenario 3; (d) Scenario 4.
Sensors 17 00936 g007
Figure 8. Comparison of updated reliability indices for inequality cases (Scenarios 1–4 in Table 3). (a) Member 13; (b) Member 27; (c) System.
Figure 8. Comparison of updated reliability indices for inequality cases (Scenarios 1–4 in Table 3). (a) Member 13; (b) Member 27; (c) System.
Sensors 17 00936 g008
Figure 9. Reliability updating results by the proposed method for inequality cases (Scenarios 5–7 in Table 3). (a) Scenario 5; (b) Scenario 6; (c) Scenario 7.
Figure 9. Reliability updating results by the proposed method for inequality cases (Scenarios 5–7 in Table 3). (a) Scenario 5; (b) Scenario 6; (c) Scenario 7.
Sensors 17 00936 g009
Figure 10. Comparison of updated reliability indices for equality cases (Scenarios 5–7 in Table 3). (a) Member 13; (b) Member 27; (c) System.
Figure 10. Comparison of updated reliability indices for equality cases (Scenarios 5–7 in Table 3). (a) Member 13; (b) Member 27; (c) System.
Sensors 17 00936 g010
Figure 11. Reliability updating results by the proposed method for inequality cases (Scenarios 8 and 9 in Table 3). (a) Scenario 8; (b) Scenario 9.
Figure 11. Reliability updating results by the proposed method for inequality cases (Scenarios 8 and 9 in Table 3). (a) Scenario 8; (b) Scenario 9.
Sensors 17 00936 g011
Figure 12. Comparison of updated reliability indices for equality cases (Scenarios 8 and 9 in Table 3). (a) Member 13; (b) Member 27; (c) System.
Figure 12. Comparison of updated reliability indices for equality cases (Scenarios 8 and 9 in Table 3). (a) Member 13; (b) Member 27; (c) System.
Sensors 17 00936 g012aSensors 17 00936 g012b
Table 1. 10 members with the maximum stresses.
Table 1. 10 members with the maximum stresses.
Member No.13278675109174
Maximum stress (MPa)32.8232.4821.2620.7320.7320.2220.1619.2618.9518.48
Table 2. Statistical properties of random variables.
Table 2. Statistical properties of random variables.
Random Variable (RV)Meanc.o.v.Distribution TypeNumber of RVs
Paris law parameter (C)1.537 × 10−12 m/cycle/(MPa·mm)m0.2258Lognormal10
Initial crack length (a0)0.11 mm1.0Exponential10
Initial crack length in repaired member (ar)0.11 mm1.0Exponential# of repair events
Detectable crack size (ad)1.0 mm1.0Exponential# of inequality events
Crack measurement error (εm)00.1 *Normal# of equality events
Live load scale factor (S)10.1Lognormal1
*: standard deviation.
Table 3. Hypothetical scenarios of inspection and repair.
Table 3. Hypothetical scenarios of inspection and repair.
CaseScenario NumberScenario Description
Inequality1No crack is detected in Member 13 (TI = 50 years & mean of ad = 1.0 mm)
2No crack is detected in Member 13 (TI = 50 years & mean of ad = 0.5 mm)
3No crack is detected in Member 13 (TI = 75 years & mean of ad = 1.0 mm)
4No crack is detected anywhere (TI = 50 years & mean of ad = 1.0 mm)
Equality50.1 mm crack is found in Member 13 (TI = 50 years)
60.5 mm crack is found in Member 13 (TI = 50 years)
71.0 mm crack is found in Member 13 (TI = 50 years)
Mixed81.0 mm crack is found in Member 13, but nowhere else (TI = 50 years & mean of ad = 1.0 mm)
90.5 mm crack is found in Member 13, but nowhere else, and Members 13 and 27 are repaired (TI = 50 years & mean of ad = 1.0 mm)
Table 4. Fatigue lives for inequality cases (Scenarios 1–4 in Table 3).
Table 4. Fatigue lives for inequality cases (Scenarios 1–4 in Table 3).
ScenarioFatigue Life (Years)
Member 13Member 27System
Initial design75.478.270.7
Scenario 188.884.279.7
Scenario 297.988.185
Scenario 392.585.881.9
Scenario 494.288.884.2
Table 5. Fatigue lives for equality cases (Scenarios 5–7 in Table 3).
Table 5. Fatigue lives for equality cases (Scenarios 5–7 in Table 3).
ScenarioFatigue Life (Years)
Member 13Member 27System
Initial design75.478.270.7
Scenario 573.375.369
Scenario 663.670.762.5
Scenario 754.265.354.1
Table 6. Fatigue lives for mixed cases (Scenarios 8 and 9 in Table 3).
Table 6. Fatigue lives for mixed cases (Scenarios 8 and 9 in Table 3).
ScenarioFatigue Life (Years)
Member 13Member 27System
Initial design75.478.270.7
Scenario 754.265.354.1
Scenario 861.665.460.2
Scenario 9>100>100>100

Share and Cite

MDPI and ACS Style

Lee, Y.-J.; Kim, R.E.; Suh, W.; Park, K. Probabilistic Fatigue Life Updating for Railway Bridges Based on Local Inspection and Repair. Sensors 2017, 17, 936. https://doi.org/10.3390/s17040936

AMA Style

Lee Y-J, Kim RE, Suh W, Park K. Probabilistic Fatigue Life Updating for Railway Bridges Based on Local Inspection and Repair. Sensors. 2017; 17(4):936. https://doi.org/10.3390/s17040936

Chicago/Turabian Style

Lee, Young-Joo, Robin E. Kim, Wonho Suh, and Kiwon Park. 2017. "Probabilistic Fatigue Life Updating for Railway Bridges Based on Local Inspection and Repair" Sensors 17, no. 4: 936. https://doi.org/10.3390/s17040936

APA Style

Lee, Y. -J., Kim, R. E., Suh, W., & Park, K. (2017). Probabilistic Fatigue Life Updating for Railway Bridges Based on Local Inspection and Repair. Sensors, 17(4), 936. https://doi.org/10.3390/s17040936

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