Abstract
In recent decades, carbon fiber reinforced plastics (CFRP) have been widely used to repair and maintain concrete structures around the world. Since the parameter uncertainties of load and resistance are very important for the reliability assessment of RC bridge strengthened by CFRP, this paper presents a method to estimate the reliability and residual life of RC bridges strengthened by CFRP. In the proposed method, uncertainties of material properties, geometry parameters, load model, and time-dependent resistance model are taken into account. The proposed method combines the inverse reliability method and the calculation method of load and resistance of RC bridge strengthened by CFRP and is illustrated by an example RC bridge strengthened by CFRP during the service stage. The results indicate that the proposed approach can provide valid information regarding parameter uncertainties for the reliability of RC bridge strengthened by CFRP during the service stage. Additionally, the effects of parameter uncertainty of the reliability and residual life of RC bridge strengthened by CFRP during the service stage are analyzed and discussed. The proposed method is more robust and reliable than the traditional method.
1. Introduction
High-performance materials, such as carbon fiber reinforced plastics (CFRP) and ultrahigh-performance concrete [1], have extremely high strength-to-weight ratios and are suitable for efficient repair of deficient structural members. These materials mentioned above can be utilized for rapid replacement and require less cost compared to traditional methods. The use of the CFRP to repair and strengthen reinforced concrete (RC) structures has been well applied.
With the increase of traffic volume and the deterioration of atmospheric environment, bridges often fail to achieve the expected service life. Accordingly, maintenance and reinforcement should be conducted regularly to meet the requirements of loading capacity and long service. In recent years, CFRP have been commonly used to strengthen concrete bridges. Gajdosova and Juraj [2] presented an investigation into the performance of slender rectangular reinforced concrete columns strengthened with CFRP in two manners. Dai et al. [3] presented the results of an experimental study on the behavior of concrete confined by CFRP jackets with a large rupture strain. Karimi et al. [4] conducted a composite column consisting of steel, concrete, and CFRP and assessed it through experimental testing and analytical modeling. Sun et al. [5] studied the seismic behavior of hollow bridge columns retrofitted with CFRP. Delgado et al. [6] presented the experimental analysis of retrofit solutions using CFRP sheets along the piers’ entire height to prevent shear failure. Gu et al. [7] presented the results of a research program that evaluated the confinement effectiveness of the type and the amount of CFRP used to retrofit circular concrete columns. Kumar et al. [8] and Baah [9] investigated the performance of CFRP to strengthen the concrete columns. Grace et al. [10, 11] applied the CFRP to reinforced T-beam and Composite Posttensioning System. Yang et al. [12, 13] studied the behaviors of CFRP in application of bridges. It can be seen that CFRP have been widely used in structural engineering but were only limited to deterministic analysis.
However, there are currently few studies concerning the change of reliability before and after reinforcement, as well as whether the service life requirements could be satisfied after reinforcement. Therefore, it is necessary to evaluate the reliability before and after reinforcement of RC bridges during the service stage. In recent years, some scholars have studied the reliability of CFRP reinforced bridges in service. Sun et al. [14] investigated a number of CFRP reinforcement tests and obtained uncertain parameters of calculation model of reinforced concrete components, and they proposed that aging and degradation of bonding surface of CFRP cannot be considered after reinforcement. Du et al. [15] established the limit state equation of a bridge strengthened with CFRP through many tests, calculated the reliability index, and verified the feasibility of the bridge reinforcement design. Forouzannia et al. [16] investigated the calibration of flexural resistance factors for load and resistance factor design of concrete bridge girders prestressed with CFRP. However, limited studies have been carried out to investigate the influence of parameter uncertainty on reliability index, and the time-dependent parameters have not been fully considered in the reliability models mentioned above.
The current bridge design code is developing in the direction of performance-based orientation, and the determination of design parameter is to meet the requirement of the prescribed reliability of bridge structures, especially the residual life of bridges during the service stage. The structure of this paper is as follows. Firstly, time-dependent resistance model of bridges strengthened with CFRP is developed. Secondly, the load model of bridges strengthened with CFRP is proposed. After that, the calculation method to assess the reliability of bridges strengthened with CFRP is introduced. Finally, with an example bridge, the time-dependent reliability index before and after reinforcement, and the residual life of RC girder bridge reinforced by CFRP are obtained, and the service life of the bridge after reinforcement is predicted using the inverse reliability method. Additionally, the variation coefficients of dead load, live load, and cross-sectional area of corroded steel bars are analyzed and discussed.
2. Time-Dependent Resistance Model of Bridges Strengthened with CFRP
Time-dependent resistance model of bridges strengthened with CFRP was developed based on initial rust time model of steel bars, corrosion rate model of steel bars, corrosion depth model of steel bars when the protective layer cracks, cracking time of the protective layer of concrete, time-dependent section area model of corroded steel bars, time-dependent model of yield strength and corrosion rate of corroded steel bars, time-dependent model of concrete strength, and cooperative working coefficient model of concrete and steel bar, which have been studied, but according to author’s knowledge, the reliability index and residual life estimation using the proposed inverse reliability method have not been investigated in the previous references. In calculating reliability index and residual life estimation using inverse reliability, the models of mean and standard deviation of time-dependent resistance after reinforcement are very important. The details of the time-dependent resistance model of bridges strengthened with CFRP are carefully described below.
2.1. Initial Rust Time Model of Steel Bars
The initial rust time of steel bars of bridges after service has been studied by many scholars. Niu and Cheng [17] established the carbonation model of concrete structures. When the atmosphere carbonates to the surface of steel bar (the carbonation depth is equal to the protective thickness of concrete layer C), the t is the time when steel bar begins to rust. Peng et al. [18] improved the prediction model of carbonization depth through taking multiple factors such as water cement ratio and calcium oxide content into account. However, the initial rust time of steel bars of real bridges calculated from the above model is mostly over 80 years, which was not consistent with the engineering practice.
Reference [19] introduced the time model of initial rust of steel bars as follows:where is the initial rust time of steel model (year); C is the protective thickness of concrete layer (cm); and is diffusion coefficient influenced by many factors such as concrete compactness and mix ratio. A range of mean values between 0.32 and 2.58 is given by [20] (); is the error function, ; is the chloride ion concentration on concrete surface (%); and is the critical chloride ion concentration, and a range of mean values (%) between 0.027% and 0.045% is given by [21].
2.2. Corrosion Rate Model of Steel Bars
Because the corrosion rate of steel bar is different before and after expansion cracking, the model of corrosion rate of steel bar is given in [22].(1)Before cracking of protective layer ():(2)After cracking of protective layer ():where is the cracking time of protective layer (year); is the corrosion rate of steel bars before cracking of protective layer (cm/year); is the corrosion rate of steel bars after cracking of protective layer (cm/year); is the correction factor of position of steel bars; 1.0 is for the steel bars in the middle and 1.6 is for those in the corner; is the environmental correction factor; 3.0–4.0 is for the outdoor humid environment; 2.5–3.5 is for the outdoor dry environment; T is the environmental temperature (°C); RH is the ambient humidity (%); C is the concrete cover thickness (mm); and is cubic compressive strength (MPa).
2.3. Corrosion Depth Model of Steel Bars When Protective Layer Cracks
The corrosion depth of steel bars when the protective layer cracks [23] is as follows:where is the corrosion depth of steel bars when the protective layer cracks (mm); is the uncertainty coefficient of calculation model; is the position influence coefficient of steel bars, 1.0 at corner and 1.35 at noncorner; C is the thickness of protective layer (mm); and D is the diameter of steel bar (mm).
2.4. The Cracking Time of Protective Layer of Concrete
The cracking time of protective layer of concrete is the initial rust time of steel bar plus the time from initial rust to cracking of protective layer. From (1), (2), and (4), we can indicate that the cracking time (year) of protective layer is
2.5. Time-Dependent Section Area Model of Corroded Steel Bars
It can be derived from the above analysis as follows:where is the initial section area of steel bar (mm2) and is the time-dependent section area (mm2).
2.6. Time-Dependent Model of Yield Strength and Corrosion Rate of Corroded Steel Bars
A yield strength model of steel bar was established based on abundant experimental studies as follows:where is the time-dependent yield strength of steel bars (MPa), is the corrosion rate of steel bars (%), and is the initial yield strength of steel bars before corrosion (MPa).
The model of steel corrosion rate can be derived from (6):
2.7. Time-Dependent Model of Concrete Strength
According to [24], the time-dependent model of concrete in atmospheric environment is as follows:where is the mean value of concrete strength measured in 28 days (MPa); is the standard deviation of concrete strength measured in 28 days (MPa); and t is the time (year).
2.8. Cooperative Working Coefficient Model of Concrete and Steel Bar
Reference [25] introduced the expression of the cooperative working coefficient of concrete and steel bar based on abundant research results and a series of tests. The expression is as follows:where is the steel corrosion rate (%).
2.9. Models of Mean and Standard Deviation of Time-Dependent Resistance after Reinforcement
When the component is the first type of T-section, the structure diagram is shown in Figure 1.

When ,where is the elastic modulus of CFRP (MPa); is the strain of CFRP; is the cross-sectional area of CFRP (cm2); is the critical height of CFRP when it reaches the allowable tensile strain and concurrent compression of concrete, ; is the ultimate tensile strain of concrete, taken as 0.0033; is the lag strain of CFRP when considering secondary loading, taken as 0 when the secondary loading is not considered; is the allowable tensile strain of CFRP, whose value should be less than or equal to the ultimate tensile strain of CFRP; is the ultimate tensile strain of CFRP, ; is the ultimate tensile strength of CFRP (MPa); is the number of layers of CFRP; is the thickness of CFRP (mm); is the relative limit height of compression zone of flexural components.
In addition,
From (10) and (11), the height of concrete compression zone x and the tensile strain of CFRP sheet can be obtained. Considering the uncertainties of material properties, geometric parameters, and calculation model, the mean value of resistance is expressed as follows:When and the CFRP sheet is larger than or equal to its allowable tensile strain, the moment of concrete in the compression zone is calculated as follows:where is the indefinite coefficient of calculation model; is the time-dependent strength of steel bar (MPa); is the time-dependent cross-sectional area of steel bar (mm2); is the bending compressive strength of time-dependent concrete (MPa).
According to (12) and (13), the expression of standard deviation can be obtained as follows:where is the partial derivative of to , and it is taken as mean value at each random variable; is the standard deviation of random variable .
3. Load Model
3.1. Dead Load Effect
From [26], the dead load does not change with time, following the normal distribution. The mean value, standard deviation, and variation coefficient of the dead load effect are as follows:where C is the coefficient of dead load effect and Gk is uniformly distributed loads. For simply supported girder bridges, .
3.2. Living Load Effect Model of Vehicle
Generally, a dimensionless parameter KSQ is used to express the ratio of the actual value to the calculated value according to the code standard of vehicle load, . Here, SQ represents the actual value and SQK represents the calculated value according to the code standard of vehicle load effect. The dimensionless parameter KSQ (shown in Table 1) of vehicle load under dense and general operating conditions can be obtained from [14].
When T = 100 years, the time-varying mean and standard deviation model of vehicle load effect are as follows:
4. Time-Varying Reliability Method
The structural time-varying function is , where R(t) represents the time-varying resistance that follows lognormal distribution, SG represents the bridge deadweight that follows normal distribution, and SQ(t) represents the live load of vehicle that generally follows extreme I distribution. When the function Z(t) is greater than 0, the structure is generally considered safe and reliable. When the function Z(t) is less than 0, the structure is in the state of failure and the failure probability of structure is the value of Z(t). When the function Z(t) is equal to 0, the structure is at the limit state of bearing capacity. The relationship between the failure probability of structures and the reliability index is expressed aswhere is the inverse function of normal distribution.
In this paper, the Monte Carlo method is used to iteratively calculate the reliability index through FERUM (Finite Element Reliability Using Matlab) software. The details of FERUM can be found in [27].
5. Life Estimation Based on Inverse Reliability Method
In this section, an inverse reliability method was proposed to estimate the residual life of bridges reinforced by CFRP in the service stage.
For a target reliability index , the problem of estimation life can be described as follows:where is the norm of a vector and is the transformation of the limit state function from the original space ( space) to the standard normal space ( space). is the vector of standard normal variables related to the design point, is the vector of basic random variables, and is the design parameter.
The inverse FORM algorithm can be solved by the following recurrence relationship:where is the vector of standard normal variables related to the design point at the th iteration; is the deterministic design parameter at the th iteration; is the vector of gradient operators with respect to ; is the vector of gradient operators with respect to; is the inner product of two vectors; and is the step size at th iteration, which is determined through a line search algorithm. The iterative algorithm proceeds iteratively until convergence is reached, i.e., whenwhere is the convergence tolerance, and to usually yields satisfactory estimates. In this paper, the inverse reliability method is used to iteratively estimate the residual life of RC beam bridges by CFRP through the program of FERUM.
The steps involved in the proposed method for estimating residual life of RC beam bridges reinforced by CFRP are as follows:(1)Select the initial values for the random variables and an unknown deterministic design parameter (based on engineering or statistical knowledge and/or assumptions). Mean values are often a reasonable choice for the random variables.(2)Initialize the iteration counter to .(3)Calculate .(4)Solve (22) and (23) to obtain and .(5)Substitute and into (20) and (21) to determine and .(6)Check the convergence criterion using (24); if unsatisfied, update the iteration counter and go to Step 3; otherwise, STOP (end of calculation).
6. Examples
The example bridge is located in China. It is a T-type reinforced concrete simple-supported beam bridge which has already had 50 years of service life [28]. The environment can be expressed as RH = 70% and temperature T = 20°C [29]. The length of the bridge is 100 m, and the deck is composed of the 7.0 m roadway and the 2 ∗ 0.15 m guardrail. There are five spans of the bridge with the 16.8 m span length. The section size of the bridge is shown in Figure 2. The concrete strength grade is C25, and the thickness of protective layer is 30 mm. The layout of the main tensioned bars is and in 3 rows. The design load of the bridge is Highway-II grade upgraded from the original Vehicle-15 grade due to the rapid development of road transport in recent years. Now the bottom of the beam is strengthened with 1.2 mm CFRP sheet, whose modulus of elasticity is 165 000 MPa and ultimate tensile strength is 2800 MPa.

(a)

(b)

(c)
Only the side beam was analyzed because it carries the most unfavourable load. Assuming that the performance of CFRP is unchanged and the bridge corrosion is delayed by 15 years after reinforcement. The mean value and standard deviation of time-varying resistance can be obtained based on the calculating steps (shown in Table 2). The results are consistent with [28].
The mean and standard deviation of side girder dead load are 566.1751 kN·m and 24.381 kN·m, respectively. The transverse distribution coefficients of the main girders were calculated, and the most unfavourable load of Highway-II grade was applied. Combining that with (17) and assuming that the bridge is in general operation state during 85% of the time and is in intensive operation state during the rest of the time, we can calculate the time-varying mean value and standard deviation of the vehicle load on the side girders (shown in Table 3). The results are consistent with [28].
The time-varying reliability index (shown in Table 4 and Figure 3) of bridges can be calculated by Monte Carlo method based on structural function. The time-varying failure probability of bridges is shown in Table 5. The calculation results are consistent with [28].

Figure 3 shows that after 50 years of service, the reliability index of the bridge is only 3.5764 and the failure probability is 1.7418 × 10−4, which cannot satisfy the requirement of the code [30], where the reliability index should be greater than 4.2 and the failure probability less than 1.3346 × 10−5. After reinforcement, the reliability index is rapidly increased to 4.611, and the failure probability is 2.0037 × 10−6. In the next 50 years, it can satisfy the requirements of reliability index and failure probability in the code. This shows that the reinforcement scheme is feasible. The results are consistent with [28].
In this paper, parameter sensitivity analysis of dead load and vehicle load was carried out. The influence of variation coefficient on bridge reliability index was mainly discussed. From Figures 4 and 5, we can see that the reliability index decreases with the increase of variation coefficient. However, the reduction of vehicle load reliability index is particularly obvious. It is confirmed that the main reason for degradation of bridge performance at present is the increase of traffic volume. The variation coefficient of bridge dead load has little effect on the reliability index.


In the above analysis, the results of reliability assessment were used in the forward first reliability method, in which the expectancy life estimation was calculated using the trial and error that was redundancy. Here, the expected life was estimated using the proposed method based on the inverse reliability approach. For example, a different target reliability index was recommended for the target reliability, and the residual life of the bridge can be seen in Table 6. From Table 6, it can be concluded that the residual life of the example bridge was close to the results of Table 4 because the target reliability was almost the same and the results are validated by forward reliability analysis. The calculation results are consistent with [28].
The results of residual life estimation of RC beam bridge reinforced by CFRP at the service stage from Table 6 indicate that, with increasing of the target reliability index, the remaining life estimation of RC beam bridge reinforced by CFRP at the service stage is expected to decrease. For more precise remaining life estimation of RC beam bridge reinforced by CFRP at the service stage, it is necessary to develop the analysis technique considering the effect of structural parameter randomness. This is very important for reasonable residual life estimation of RC beam bridges reinforced by CFRP at the service stage, which exhibit wide dispersion in structural parameters. This problem can be solved by the proposed method. The proposed method is more robust and reliable than the traditional method.
A parametric analysis was carried out to study the influence of the mean value of random variables and the coefficients of variation of random variables on the residual life estimation of RC beam bridge reinforced by CFRP at the service stage. For the sake of simplicity, only strengthening bridge at the target reliability 4.2 in the previous section is considered.
When keeping the coefficients of variation of all random variables unchanged and only varying the mean values under a ±15% of all random variables, the residual life estimations of RC beam bridge reinforced by CFRP at the service stage are compared in Figures 6–8. With the increase of the mean value of random variables of dead load and live load, the residual life estimation of RC beam bridge reinforced by CFRP at the service stage is expected to decrease. While an opposite influence can be observed for the random variables of CFRP. With the increase of the mean value of the random variables of CFRP, the residual life estimation of RC beam bridge reinforced by CFRP at the service stage is expected to increase.



By varying the coefficients of variation of random variables, the residual life of RC beam bridge reinforced by CFRP at the service stage is estimated using the proposed method mentioned previously. When keeping the mean values of all random variables unchanged and only varying the coefficients of variation under a ±15% of all random variables. The results are illustrated in Figures 9–11. It is observed that the residual life estimation of RC beam bridge reinforced by CFRP at the service stage decreases as the coefficients of variations of all random variables increase.



7. Conclusion
In this study, a new method is proposed to estimate the reliability of RC bridge strengthened by CFRP, and uncertainties of material properties, geometry parameters, load model, and time-dependent resistance model are taken into account. The proposed method combines the concepts of inverse reliability and the calculation method of load and resistance of RC bridge strengthened by CFRP. A RC bridge strengthened by CFRP is selected as an example to verify the proposed method, and the main conclusions and findings can be summarized as follows.(1)Based on the limit state function established in this paper, the Monte Carlo method is used to calculate the reliability index of the RC bridge strengthened by CFRP.(2)Taking an existing T-beam bridge as an example, the time-varying reliability before and after reinforcement is calculated. Results show that the reliability index decreases with aging. After strengthening the RC bridge using CFRP, the reliability index will increase.(3)The effects of various parameters on the reliability of RC bridge strengthened by CFRP during the service stage are analyzed. The results indicate that the variation coefficient of dead load, vehicle load, and cross-sectional area of corroded steel bars could affect the reliability of the RC bridge. Additionally, it is found that the variation coefficient of live load is more important than that of dead load and cross-sectional area of corroded steel bars.(4)For more precise residual life estimation of RC girder bridge reinforced by CFRP at the service stage, it is necessary to develop analysis technique considering the effect of structural parameter randomness. The proposed method is more robust and reliable than the traditional method.
Data Availability
The data supporting the conclusions of the study are included within the manuscript.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors would like to acknowledge the financial support from the Fundamental Research Funds for the Central Universities (Grant no. 300102210105), Natural Science Basic Resarch Program of Shaanxi (Program no. 2020JQ-377), Guangdong Provincial Transportation and Transportation Office Technology Funding Project (2016-02-016), and Shanxi Natural Science Basic Research Funding Project (201801D221217).