Abstract

Concrete creep plays a significant role in the long-term performance of the prestressed concrete structure. However, most of the existing prediction models cannot accurately reflect the in-site concrete creep in a bridge construction environment. To improve the prediction accuracy of creep effects in concrete structures, an innovative creep analysis method is developed in this study. Parameters in the creep model in fib MC 2010 have been calibrated with respect to the long-term loading test results of the prestressed concrete beam. The measured strains of concrete and the midspan deflections of the test beam are compared with the predicted results using the creep model in fib MC 2010. It indicates that the results predicted by the calibrated creep model are in good agreement with the test results. However, the results predicted by the creep model in fib MC 2010 significantly deviate from the test results. This proposed creep analysis method can provide a new thought to improve the predicted effect of the creep effects on creep-sensitive structures.

1. Introduction

Concrete creep plays a significant role in the long-term performance of concrete structures [15]. Since the discovery of concrete creep in the early 20th century, relevant scholars have studied it for more than a hundred years. Yet, despite the long history of research, the understanding of concrete creep phenomenon is still far less than complete [69]. Except for B3 and B4 creep models developed by Bažant based on microprestress-solidification theory [10, 11], creep models developed by other scholars and design codes are empirical models based on laboratory tests [1218]. Creep in concrete is a highly complex phenomenon, which depends on a wide range of factors. The long-term deflections predicted by the creep models in different design codes have a great difference [1921]. Bažant et al. [22] analyzed the long-term deflection of Koror-Babeldaob (KB) Bridge in Palau, which collapsed in 1996, using ACI 209, CEB MC 1990, JSCE, GL2000, and B3 creep models. Granata et al. [23] investigated the effects of time-dependent phenomena on prestressed concrete girder bridges and made a compared analysis on the creep behaviors for actual cases of bridges by using the creep prediction models in Eurocode 2, ACI 209, B3, GL2000, and AASHTO 2007. Pan et al. [24] studied the creep strain of the high-strength concrete used in the continuous rigid frame of Sutong Bridge and compared it with the strains estimated by ACI 209, JTG D62-2004, B3, and GL2000 models. Vandewalle [25] investigated the effects of environmental conditions on creep compliance function experimentally and predicted the strain by the models in Eurocode 2 and B3. Yang et al. [26] carried out the experimental analysis of 32 m precast PC box girder and compared it with the theoretical results calculated by JTG 2004, fib MC 2010, and GL2000. These comparisons show the significant differences between the predictions supplied by these models especially for delayed deflections at very long times.

So far, many large-span bridges suffer from excessive long-term deflection [22, 27, 28], and the causes of this phenomenon are complicated, among which concrete creep is one of the main factors [29]. Accurate prediction of the long-term behavior of the long-span prestressed concrete bridges, reasonable creep, and shrinkage model is indispensable [30]. Due to the complexity of the actual environment condition, it is difficult to simulate accurately. For large and creep-sensitive structures, like long-span bridges, in-site creep test is an alternative way to predict creep in concrete for actual structures accurately. Furthermore, relevant design codes [17, 31] recommend this method to determine creep deformation for the actual structures. However, most concrete creep tests are carried out at the standard environmental condition, which is significantly different from the actual environment. The most tough thing is that the concrete creep model based on laboratory test is not reliable in practical engineering. Therefore, it is urgent to establish a creep analysis method suitable for the structures in site conditions.

In this study, long-term creep tests of six T-section beams with different prestressing levels under nonstandard environmental condition are carried out to simulate the in-site concrete creep in the bridge construction environment. Based on the test results, an innovative creep analysis method is developed and the parameters in the creep model in fib MC 2010 are calibrated by the least-squares optimization program. The measured strains of concrete and the midspan deflections of the test beams are compared with the predicted results using the creep model in fib MC 2010 and the calibrated creep model proposed in this paper. The objective of this study is to determine the parameters of the creep model with respect to the concrete creep test for a specific engineering project and propose the concrete creep model that meets the actual in-site environmental conditions. It should be noted that the creep model provided here can only be used to calculate the creep of concrete structures under the conditions in this study. The proposed creep analysis method can provide experience for the analysis of structural creep effect.

2. Experimental Program

2.1. Details of Specimens

Grade R52.5 Portland cement, gravel with continuous grading size of 5–10 mm and river sand, is used in this test, and 12% fly ash and 4% slag powder are added to the concrete mixture. Table 1 shows the concrete mixture proportions. When casting test beams, twelve 150 mm × 150 mm × 150 mm cube specimens and six 150 mm × 150 mm × 300 mm prism specimens were also cast, which were used to test the compressive strength and elastic modulus of concrete, respectively. Among them, six cube specimens and three prismatic specimens cured in the same conditions as the T-section beam.

Six post-tensioned prestressed concrete T-section beams identified as PC1∼PC6 with a length of 4 m were produced. The details of the test prestressed concrete beam are shown in Figure 1(a). The prestressing tendon is made of a single strand, with Young’s modulus of 1.95 × 105 MPa, cross-sectional area of 139 mm2, and ultimate tensile stress of 1860 MPa. To prevent the local crushing failure of concrete, steel anchorage plate and reinforced net are provided. Details of anchorage are shown in Figure 1(b).

Prestressed concrete T-section beams are divided into two groups according to the concrete age at tensioning. Table 2 lists details of stress in prestressing tendons of each test beam, in which fpt is the ultimate tensile stress of prestressing tendons. The embedded plastic bellow is grouted immediately after tensioning. The grouting material is cement mortar with a water-cement ratio of 0.4 and an expansive agent ratio of 10%.

2.2. Creep Test

A four-point bending strategy is adopted in this experiment. Figure 2(a) shows the prestressed concrete T-section beam in the test. Figure 2(b) shows the arrangement of measuring points (MA and MB). The embedded vibrating wire strain gauges (VWSG) are used to measure the strain at each measuring point and the dial indicators are used to measure the deformation development at midspan of the beam. Figures 3(a) to 3(f) illustrate the loading histories of the six beams. It should be noted that loads of PC1, PC2, and PC3 exhibited large fluctuations due to the use of the self-locking hydraulic jacks at the initial stage. Therefore, the self-locking hydraulic jacks were replaced by the mechanical jacks in the subsequent creep tests.

2.3. Test Results and Discussion

The compressive strength test and modulus of elasticity test of concrete were carried out at the age of 28 days. The cubic compressive strength of standard curing specimen is 46.86 MPa and the elastic modulus of prism specimens is 33.62 × 104 MPa. The compressive strength of specimens in the same curing environment as the T-section beam is 42.5 MPa and the elastic modulus is 32.47 × 104 MPa. It can be seen that the difference of the curing environment will affect the compressive strength and modulus of elasticity of concrete.

Figures 4(a) and 4(b) demonstrate the variation of the average daily temperature and the average surrounding relative humidity of the measurement points during the test. It can be seen from Figure 4 that there are two types of fluctuation for temperature and relative humidity; i.e., the small local fluctuation reflects the difference between the daily average temperature and relative humidity, and the large range represents the seasonal variation.

Based on the stress and strain in concrete, the test can be divided into three stages, as shown in Figure 5. From casting concrete to tensioning (Stage I), the strains in this stage are mainly shrinkage strain of concrete. Considering the restraining effect of non-prestressed steel on concrete shrinkage, it will produce tension stress in concrete and compressive stress in non-prestressed steel. From tensioning to apply the external loads (Stage II), the stress in concrete and non-prestressed steel is mainly caused by the prestress force. Due to the time-dependent strain of concrete, the stress and strain in concrete, non-prestressed steel and prestressing tendons also change with time. With the increase of shrinkage and creep strain of concrete, the stress in prestressing tendons and concrete decreases, and the creep rate of concrete shows a decreasing trend. From the application of the external load to the observation time (Stage III), there is a sudden change of concrete strain due to the external load, which becomes the dominant factor of the concrete stress in this stage. The strain and stress in non-prestressed steel and prestressing tendons still change with the strain of concrete. Table 3 lists the details of the strain at the measuring point for PC1 to PC6. Here, t0 and tF are the time of prestressing and applying the external load, respectively. εmc,A(t0) and εmc,B(t0) denote the strains of Points MA and MB at the time of t0 under the prestress, respectively. εmc,A(t0F) and εmc,B(t0F) denote the strains of Points MA and MB at the time of tF under the prestress, respectively. εmc,A(tF) and εmc,B(tF) represent the strains of Points MA and MB at the time of tF under the prestress and external load, respectively.

Figure 6 presents the curves of the total strain of concrete (temperature strain, shrinkage strain, and creep strain) versus time measured at MA and MB points at the midspan of the specimens PC1-PC6. From Figure 6, it can be seen that the total strain of concrete measured at MA and MB is closely related to the temperature variation trend, as indicated in Figure 4(a). During the 200–300 d, 500–680 d, 850–1050 d, and 1250–1350 d, the temperature tends to rise, and the total strain of concrete tends to increase rapidly.

3. Creep Analysis Method

3.1. Calculation of the Strain

After anchorage, the internal forces of concrete, non-prestressed steel, and prestressing tendons are under equilibrium condition (static equilibrium); that is,where t0 is the concrete age in days at tensioning; t is the concrete age in days at the moment considered; is the tensile force in concrete at time t; and are the tensile force and the stress in non-prestressed steel, respectively; and are the tensile force and the stress in prestressing tendons, respectively; and and are the cross-sectional area of prestressing tendons and non-prestressed steel, respectively.

Assuming that concrete is perfectly bonded with non-prestressed steel and prestressing tendons (deformation compatibility condition), thenwhere and are the change in strain of concrete at the centroid of non-prestressed steel and prestressing tendons, respectively.

The concrete stress is not constant due to the interaction of the concrete creep and the prestressing loss. On the one hand, the concrete creep results in the prestressing loss. On the other hand, the prestressing loss reduces the concrete stress, which in turn affects the development of concrete creep. For variable stresses condition, the principle of superposition is assumed to be valid. Therefore, the constitutive equation for concrete can be expressed as [12]where is the creep strain of concrete from time t0 to time t; and are the change in shrinkage strain and thermal strain of concrete, respectively; is the initial stress of concrete at time t0; is the change in temperature; and are the temperature of concrete based on the test data at time t and t0, respectively; and is the thermal expansion coefficient of concrete. for normal weight concrete suggested by fib MC 2010 is used in this study. is the creep compliance, which can be expressed as follows according to fib MC 2010 [12].where and are the elastic modulus of the concrete in MPa at the age of t0 and 28 days, respectively; is the creep coefficient of concrete, and t0 is the initial loading age of concrete in days. It should be noted that the relative humidity takes 60% as the annual average ambient relative humidity based on the measured data.

The stress of non-prestressed steel can be expressed aswhere is the strain in non-prestressed steel from t0 to t; is Young’s modulus of the non-prestressed steel.

Time-dependent stress of prestressing tendons, which includes stress changes caused by deformation and relaxation of prestressed tendons [31], can be evaluated aswhere is the strain in prestressing tendons from t0 to t; is the elastic modulus of the prestressing tendons; is the prestress loss due to relaxation of the prestressed tendons; and is the relaxation reduction factor. Generally, is in the range of 0.5∼0.9 [3235]. For simplicity, is adopted in this study according to the Eurocode [31]. is the absolute value of stress change in prestressing tendons caused by relaxation and the expression of can be found in fib MC 2010.

Combined with material mechanics and the plane-section assumption, the stress of concrete at the location interested can be expressed aswhere is the cross-sectional area of concrete; represents the moment of inertia of cross-section of concrete; is the bending moment caused by the self-weight of the beam and the external load. In stage II, , and in stage III, , in which is the bending moment induced by the external load. denotes the distance from the location interested to the centroid of the concrete section; and are the distance from the centroid of the non-prestressed steel and the prestressing tendons to the centroid of the concrete section, respectively.

According to equation (7), the concrete stress at the centroid of the non-prestressed steel and the prestressing tendons at time t can be expressed as follows:where is the distance from the measuring point to the centroid of the concrete section; and are the reinforcement ratio of non-prestressed steel and prestressing tendons, respectively; and , in which i and j represent the p, s and m, where the subscript pc, sc and mc represent the concrete at the centroid of prestressing tendons, the concrete at the centroid of non-prestressed steel, and the concrete at the measuring points, respectively. If eic and ejc are on the same side of the centroid of the section, take “+.” Otherwise, take “−.”

Figures 7(a) and 7(b) show the stress and strain distributions for the section of the test beam in Stage II and Stage III , respectively. In Figure 7, and are the distance from the measuring points MA and MB to the centroid of the concrete section, respectively. and are the stress in concrete at the measuring points MA and MB, respectively.

The instantaneous prestress force measured from the load sensors after anchorage is taken as the initial prestress force. The bending moment at time t0 is determined with respect to the self-weight of the concrete. According to the initial conditions, and equations (1)–(3), (5), and (8a), the tensile force in the non-prestressed steel at time t0 can be obtained as follows:

The step-by-step approach is used in the following analysis. The time interval is subdivided into a number of subintervals with variant time increment. The time interval Δti increases from 0.01 days to 10 days in this study (Δti = 0.01 d for t ≤ 7 d; Δti = 0.5 d for 7 d < t ≤ 30 d; Δti = 1 d for 30 d < t ≤ 365 d; Δti = 5 d for 365 d < t ≤ 700 d; Δti = 10 d for t > 700 d). The stress increment is assumed to be constant in a time interval Δti and applied at the beginning of the ith interval as can be seen in Figure 8. Therefore, equation (3) can be rewritten as [36]

It should be noted that the stress of concrete in Stage I is minimal. Therefore, it is not included in the strain analysis of the test beams.

3.2. Calibration of the Creep Model and Verification

In general, the creep models given in the codes are based on a large number of laboratory tests, which generally reflects the trends of concrete creep. However, there are significant differences between the creep strain calculated by the creep model in the code and the measured strain in actual structures. The main reason is that the environmental conditions on the site are quite different from the environmental conditions in the laboratory. For large-span prestressed concrete bridges, which are sensitive to concrete creep, it is necessary to establish a creep analysis method suitable for site conditions. Therefore, an innovative creep analysis method is proposed in this study based on member tests.

The creep coefficient in the proposed creep analysis methodology can be expressed aswhere αc1, αc2, αc3, αc4, αc5, αc6, and αc7 are parameters.

According to the principle of least squares, the estimators of parameters given in equation (11) are the solution of the following optimization problem:where is the strain of concrete measured from the test at time ti; is the calculated strain of concrete at time ti; n is the total number of measured values.

To quantify the difference between the calculated and the measured strains more clearly, the errors are normalized and the error coefficient is defined as follows in this paper:where is the average strain of concrete at all measuring points.

The flow chart for determining parameters of the creep model is shown in Figure 9.

The optimization program is developed according to the flowchart mentioned previously. Table 4 lists the parameters obtained from the optimization program and the corresponding error coefficient. It is generally accepted that the smaller the , the smaller the error between the calculated value and the measured value, and vice versa. It can be seen from Table 4 that the creep method presented in this paper can significantly reduce the error coefficient. Figure 10 shows the comparison between the measured concrete strains and the concrete strains predicted by the creep model in fib MC 2010 and the calibrated creep model. It can be seen that the concrete strains predicted by the calibrated creep model are in good agreement with the test strains. However, the concrete strains calculated with respect to the creep models provided in fib MC 2010 are smaller than the measured values. This means that the creep models given by fib MC 2010 do not necessarily accurately predict concrete strain for the specific conditions. It should be noted that the thermal strains of concrete have been subtracted from the total strains to compare clearly. It can be seen from Figure 10 that, after subtracting the thermal strains, the flat branches and ascending branches of concrete strain evolution curves as shown in Figure 6 vanish.

4. Analysis of the Deflection

Based on the creep analysis methodology mentioned previously, the deflections of the test beam at the midspan are calculated and compared with the test results below. The relationship between the deflection and curvature can be described aswhere is the radius of curvature for the section x at time t; x is the distance from the section interested to the support; and are the concrete strain for the section x at the centroid of the prestress tendons and the non-prestressed steel, respectively.

The beam is divided into (m − 1) segments, each with a length of Δx along the length direction, as illustrated in Figure 11. The deflection of the k − 1th, kth, and k + 1th segments at time t is , and , respectively. The radius of curvature of the kth segment at time t is . Figure 12 shows the geometric relationship between deflection and radius of curvature and the deflection angle . Considering that the radius of curvature is much larger than Δx, the relationship between the deflection angle and the deflection and the relationship between the deflection angle and the radius of curvature can be approximately expressed as

According to equations (15a)–(15c), the following expression can be obtained:

The linear system of equations can be expressed as

According to the boundary conditions and equations (14) and (17), the deflection of the midspan can be obtained.

Figure 13 demonstrates a comparison between the predicted midspan deflections and the measured deflections. From Figure 13, it can be seen that the midspan deflections calculated by the calibrated creep model are in good agreement with the measured deflections; however, the creep model in fib MC 2010 underestimates the midspan deflections, which indicates the superiority of the proposed method.

5. Conclusions and Remarks

This paper deals with the realistic creep analysis for practical prestressed concrete structures. Long-term creep tests of six T-section prestressed concrete beams with different prestressing levels are carried out. The strains of concrete and the midspan deflections of the test beam are monitored and analyzed. To improve the prediction accuracy of creep effects in concrete structures, an innovative creep analysis method on the basis of the prestressed concrete test results is proposed. The analysis in this study is based on the creep model recommended in fib MC 2010. This paper can be considered to have successfully achieved the following objectives:(1)Developed an innovative creep analysis method based on the superposition principle of creep and the creep model in fib MC 2010.(2)Calibrated the corresponding parameters of the creep model on the basis of the optimization program and the prestressed concrete test results.(3)Successfully predicted the actual concrete strain evolution and verified the innovative creep analysis method by the deflection of the midspan of test beams.

The proposed method can provide a reference for creep prediction of prestressed concrete structures, which are sensitive to concrete creep and need a creep model suitable for the structures in site conditions. It must be pointed out that the calibrated creep model can only be used to calculate the creep of a certain concrete under the specific field conditions and cannot be used for other situations. In addition, the concrete shrinkage model used in long-span prestressed concrete structures should also be modified according to the field measured data.

Data Availability

Some or all data, models, or code generated or used during the study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors gratefully acknowledge the financial support from the National Basic Research Program of China (973 Program) under the grant of 2015CB057703.