Effect of section geometry on development of shrinkage-induced deformations in box girder bridges

Non-uniform shrinkage strains can lead to significant additional deflections in large box girder bridges, leading to serviceability problems. This article examines experimentally and analytically the effect of different cross-section geometries on the shrinkage camber of bridge box girders. Small-scale beams were tested to determine the development of shrinkage strains across the beams of depth. Parameters investigated include cross section thickness, drying conditions, and type of concrete mix. Based on the experimental results, inverse analysis is utilised to obtain a surface factor and a hydro-shrinkage coefficient. In this study, such vales are used to determine, for the first time, shrinkage-induced bending deformations of long-span bridges using a hydro-mechanical approach. The results are then used to examine numerically the effect of different section geometries on the development of shrinkage camber. It is shown that the analytical predictions match the experimental results with an accuracy of 85%. A further parametric study is carried out to investigate the effects of specimen geometry and ambient relative humidity. The hydro-mechanical approach is further validated using shrinkage field data from the 230 m two-span box girder Yiju River Bridge (China). The approach proposed in this study is expected to contribute towards improving the predictions of the long term behaviour of box girder bridges and towards better bridge management.


Introduction
Box girders are widely used in the design and construction of prestressed concrete (PC) bridges due to their structural efficiency and cost-effectiveness. However, in recent years numerous of such bridges experienced excessive vertical deflections at mid-span, affecting significantly serviceability and aesthetics [1]. Previous research has attributed these large deflections to additional curvatures produced by non-uniform shrinkage strain across the depth of concrete girders [2,3]. Changes in the sectional shape and size are also expected to contribute to the development of shrinkage camber as the drying rate is not uniform over the cross section. To date, however, the shrinkage camber caused by sectional shape has not been examined experimentally.
The process of drying shrinkage is caused by diffusion of moisture from the inside to the outer surface of concrete when this is exposed to a drying environment. The most effective way to mitigate the effect of shrinkage deformations is to prevent the loss of moisture and to protect the concrete from rapid drying [4]. However, the complete elimination of internal deformations is practically impossible. As a result, it is important to develop practical ways to account for shrinkage strains during construction as well as for structural assessment required for bridge management.
Current design guidelines include approximate analytical models to predict the long-term shrinkage behavior of concrete (e.g. ACI 209 [5] and Model Code 2010 (MC 2010) [6]). These models account for factors such as type of concrete, ambient relative humidity (RH), as well as cross section size and shape. In these models, the theoretical cross section thickness is defined by an average drying path length representing the entire cross section [3], which reflects the expectation that uniform drying shrinkage produces only axial shortening deformation but no deflection. While the above models are widely used in bridge design/analysis, they were proved unsuitable for calculating deflections of long span bridges with box girders [7,8].
Different studies have investigated non-uniform shrinkage effects in box girder sections. Using the KB bridge (in Palau) as case study, Křístek [8] estimated the different shrinkage and creep effects of box girder bridges by assigning different shrinkage properties to each part of a whole box section, using different thicknesses for the web, top and bottom slabs. Křístek concluded that this approach can produce more accurate predictions than using the whole box section for calculating long-term deflections and stress redistributions. A more accurate 3D solid model of the KB bridge identified that the non-uniform distribution of drying shrinkage within the box girder section (due to the variation in thickness among different parts of the section) was one of the causes for excessive vertical deflections [2]. Malm and Sundquist [9] examined the shrinkage behaviour of the hollow box-girder Gröndal bridge (Sweden) using finite element (FE) analysis. The results showed that, after 2 years of service, nonuniform drying shrinkage accounted for 18% of the total bridge deflection. However, the above studies ignored the fact that shrinkage and creep depend on the section thickness, which varies across the cross section and along the bridge. Therefore, to accurately simulate the behaviour of an entire bridge, a large number of different geometries and shrinkage and creep models are needed, which makes the analysis computationally demanding and time-consuming.
To address the above drawbacks, more precise (yet complex) moisture diffusion analyses have been used to evaluate differential drying shrinkage deformations. For instance, numerical methods based on finite differences have been adopted to analyse the effect of internal RH changes over the cross section depth, along with a shrinkage coefficient that relates the internal RH and shrinkage deformations [10][11][12][13]. Multi-field models that consider thermal-diffusion coupling effects were also used to assess the thermalshrinkage stress of large concrete structures [14][15][16][17]. Theoretically, the above mentioned approach based on hydro-mechanical analysis can be used to predict the time dependent deformations of box girders with complex cross sections. However, the input parameters in the diffusion model (e.g. diffusion coefficient, surface factor and hydro-shrinkage coefficient) affect the accuracy of the predictions. The moisture diffusion coefficient depends on the internal RH of concrete and reflects the moisture transfer during drying. Existing guidelines (e.g. Model Code 90 [18] and MC2010 [6]) suggest analytical expressions to compute the moisture diffusion coefficient. The surface factor (also called transport coefficient or convection coefficient) represents the velocity of internal moisture that evaporates through the concrete surface to the environment. The surface factor depends on the w/c ratio and the moisture gradient and, consequently, its value should be determined on a case by case basis. Previous research also found that surface factor values had a large scatter [19,20]. The hydro-shrinkage coefficient (also called contraction coefficient and shrinkage coefficient) is the coupling condition between free shrinkage strain and varying moisture. Typical values of hydro-shrinkage coefficient range from 1.2 to 3.0 mm/m [21,22]. The surface factor and the hydro-shrinkage coefficient can be calibrated using a combination of tests results and parametric studies.
This article examines experimental and analytically the effect of cross-section geometry on the development of shrinkage deformation of box girders of bridges. Section 2 of this study discusses methods for calculating shrinkage camber. Subsequently, Sect. 3 presents and discusses results of an experimental programme aimed to obtain surface factors and hydroshrinkage coefficients from beam specimens. Section 4 proposes to use, for the first time, a hydromechanical coupling model to simulate the shrinkageinduced bending deformation of the tested beams. The model is further validated by predicting the deflections of an actual case study bridge (Sect. 5). Concluding remarks of this study are given in Sect. 6. This article aims to contribute towards the development of more accurate models to assess the overall deflection of box girder bridges so as to improve the performance, serviceability and economic viability of structures.
2 Calculation of shrinkage camber and moisture diffusion

Shrinkage camber
Non-uniform shrinkage over the beam depth can produce both shortening deformation and bending deformations along the beam. By assuming that every small element in the cross section can deform freely, the free shrinkage strain e sh at an element located at (x, y) from the onset of shrinkage strain can be calculated as follows: where b is the hydro-shrinkage coefficient; and DH is the change of pore humidity at (x, y). If a plane strain distribution along the beam depth is assumed, the shrinkage curvature w sh can be obtained: where t s and t are the time at the onset of shrinkage strain and the analysis time, respectively; y À y c is the distance from the element to the neutral axis of the whole section; I is the moment of inertia of the whole section; and A is the area of the whole section. Likewise, the longitudinal average strain e sh is: where all the variables are as defined before. It should be noted that the internal reinforcement restrains the development of shrinkage-induced deformations, and that this can lead to stress and strain redistributions over the cross section's height. However, in current design practice, the restrain of reinforcement is usually neglected in typical prestressed bridges with relative low reinforcement ratios. As a result, such restraining effect are neglected in subsequent analysis of the case study bridge presented later in this article.

Moisture diffusion over cross section
The time-dependent shrinkage of concrete consists of four main sources, including drying, autogenous, plastic and carbonation shrinkage. Water in concrete is consumed during hydration, and therefore the internal moisture decreases and affects the humidity diffusion process. As a result, autogenous shrinkage is significant in high-strength concrete during the early ages. According to Fick's second law, the diffusion equation with the effect of self-desiccation [23] can be expressed as: where H is the pore humidity; D(H) is the moisture diffusion coefficient; oH a ot is the change rate of internal humidity due to self-desiccation; H 0 ðx; yÞ is the initial condition function which denotes the initial internal humidity distribution at an element located at (x, y); f is the surface factor; H en is the environmental humidity; and H s is the relative humidity on the exposed face S.
Equation (4) can be solved using finite difference methods that solve time-dependent nonlinear diffusion problems [24]. Accordingly, the FE governing equation in nonlinear diffusion analysis is [25]: where D ½ is the element humidity diffusion matrix; C ½ is the element humidity capacity matrix; and R ½ is the element boundary humidity flux. By using Euler backward implicit time integration and full Newton-Raphson iteration methods, the nonlinear transient equation can be rewritten as: where Dt is the time step; i is the iteration number; and DH f g is the element humidity variation matrix from t to Dt. According to the Euler backward difference method, i tþDt _ H È É can be expressed as: The hydro-shrinkage coefficient b links the moisture diffusion model and the structural shrinkage model described in Sect. 2.1. The time-dependent internal humidity distribution can be determined using Eqs. (5)-(7), whereas the shrinkage strain and curvature can be calculated using Eqs. (1)- (3). In this study, ADINA Ò [26] software is used to simulate the drying process in concrete. The decrease of internal humidity due to self-desiccation is modelled by the timedependent uniform drying process, as shown later in Sect. 4.1.
In summary, to determine the development of shrinkage-induced deformations using the hydro-mechanical approach described in the previous sections, the following main assumptions and calculation steps are necessary: 1. Calibrate the surface factor f and the hydroshrinkage coefficient b using inverse analyses; 2. Calculate the time-dependent internal humidity distribution assuming that only drying and selfdesiccation occur; and 3. Assume a plane strain distribution along the beam depth and determine both shortening deformation and bending deformations along the beam using section analysis.
3 Experimental programme

Test specimens and exposure conditions
Four non-symmetric I cross section beams (group I) were produced and tested to examine the camber due to non-uniform shrinkage. The beams were cast sideways to eliminate differences in concrete properties along the height of the beam. The simplified scaled-down geometry of the short I beams replicated approximately half of a box girder section, as shown in Fig. 1. The beams had a total depth of 300 mm, a total length of 500 mm and a web thickness of 50 mm, as shown in Fig. 1. The thicknesses of the top and bottom slabs (flanges) were 50 and 100 mm, respectively. For comparison, four control concrete prisms (specimens B and BS) of 100 9 100 9 450 mm were also cast.
No reinforcement was used in the beams and prisms. The I beam specimens were subjected to different exposure conditions to promote concrete drying shrinkage. To amplify the effect of differential shrinkage, the surfaces of the web and bottom slab of two specimens (IP-1 and IP-2) were epoxy coated to prevent drying, whereas two beam counterparts (I-1 and I-2) were fully exposed to the environment. Two prisms (BS) were fully epoxy coated on all surfaces, whereas two prisms (B) had only one surface exposed to the environment. Table 1 summarises the main characteristics of the tested specimens and their corresponding the mean concrete compressive strength f c .

Material properties
The I beams and prisms were cast using two batches of normal-strength concrete (mixes C1 and C2) typical of bridge construction in China. To investigate the effect of water-cement (w/c) ratio on drying and autogenous shrinkage, w/c ratios of 0.37 and 0.40 were selected for mixes C1 and C2, respectively. The following proportions were used for mix C1: cement = 460 kg/m 3 , gravel = 1270 kg/m 3 , sand = 520 kg/m 3 , whereas for mix C2: cement = 410 kg/m 3 , gravel = 1275 kg/m 3 , sand = 545 kg/m 3 . For each batch, f c was obtained from concrete cubes tested at 28 days according to the Chinese JTG D60 bridge design code [27]. All specimens were moist-cured for 3 days following casting, after which the epoxy coating was applied onto the faces of the selected specimens. During testing, all specimens were kept in standard laboratory conditions. It should be noted that the tested beams have the same exposure conditions on the inner and outer faces of the simulated cross section. This is because, in actual bridge construction practice, the inner parts of the box sections are also exposed to the environment during the construction period of the cantilever segments (which can last for several months or years).
The I beams were rotated 90°and simply supported on two steel rods to ensure that only free deformations occurred during the test. To measure the development of strain during drying, two vibrating wire strain gauges were installed along the centerline of the beams: one at the top slab and one beneath the bottom slab, as shown in Fig. 1. The strain gauges had a measurement range of ±1500 le and an accuracy of 1 le. The prisms were also supported on two steel rods and monitored by vibrating wire strain gauges. Measurements were read every day during the first week of testing. The beams and prims were monitored for a period of 115 days. Figure 2a-b show the evolution of shrinkage strains measured at the top (Top test data points) and bottom (Bot test data points) of the un-sealed specimens I-1 and I-2. As expected, the shrinkage strains were different at these locations, thus confirming that both axial shortening and shrinkage bending occurred during the tests. The difference between the measured data from top and bottom strain gauges (Top-Bot test) indicates the magnitude of the camber deformation, which results in additional bending deformations. The results indicate that during the first 58 days, the shrinkage strain on the top developed faster than that on the bottom, while the bending deformation increased. However, after 58 days, the camber decreased after reaching a peak value. These results were expected since, at the initial stages of the test, the development of drying shrinkage at the (thinner) top slab was faster than at the bottom slab. Consequently, an increasing upward camber deformation occurred over time but only up to a peak point. After this, drying shrinkage at the top slab developed more slowly than at the bottom slab because the majority of internal moisture in the former has been consumed by convection. The Top-Bot curve decreased with time, thus leading to a downward camber deformation. Figure 2c-d also show the evolution of shrinkage strains at the top and bottom of the sealed specimens IP-1 and IP-2, respectively. As expected, the difference between shrinkage strains on the top and bottom slabs is up to 55% higher than those of the corresponding un-sealed specimens (shown in Fig. 2a-b). However, shrinkage strains were still measured on the sealed surface of the specimens. This can be attributed to autogenous shrinkage but also to restraining effects of the sealed flange.

Numerical analysis
Based on the experimental results from Sect. 3, this section presents an inverse analysis performed to obtain the surface factor and the hydro-shrinkage coefficient. A parametric analysis is then carried out to examine the influence of specimen geometry and ambient RH on the shrinkage behaviour of the tested beams. It should be mentioned that solar radiation can significantly affect the long term drying process in actual concrete bridges. However, to simulate actual environmental conditions of a bridge, a more complex hydro-thermo-mechanical model would be needed to couple the effect of moisture, temperature and deformation fields, thus resulting in very demanding computational analyses. As a result, current guidelines use ''average'' RH and environmental temperature, which are also adopted in the calculations of this study. To determine suitable parameters for a hydro-mechanical model, a moisture diffusion analysis is first performed. This study adopts the approach suggested by the recent MC2010 [6], where the moisture diffusion coefficient D(H) is calculated as: where H is the relative pore humidity; A surface factor f and a hydro-shrinkage coefficient are required for the moisture diffusion analysis.
Previous literature indicates that the surface factor f depends on the water-cement ratio and the moisture gradient, and a wide range of values were suggested [19,28]. Likewise, the hydro-shrinkage coefficient was found to range between 1.1 9 10 -3 and 3.0 9 10 -3 [21,22].
In addition, free shrinkage strain is difficult to measure directly due to internal restraining effects. The FEA package ADINA Ò [26] was used to perform inverse analyses and estimate the surface factor and the hydro-shrinkage coefficient [19]. The results indicate that the best match is obtained for f values of 0.3 and 0.5 mm/day for mixes C1 and C2, respectively (see Fig. 3). Likewise, the use of a hydro-shrinkage coefficient values of 0.0014-0.0015 led to good predictions for both concrete mixes. Consequently, these values are used in subsequent analyses.
At an early age, the internal moisture in concrete is consumed due to self-desiccation, thus leading to     autogenous shrinkage. Since the experiments were conducted in a relatively short time, the effect of concrete self-desiccation needs to be accounted for in the analysis. As water is consumed by hydration, the internal RH decreases and interacts with the humidity diffusion process. The autogenous shrinkage depends on the water to binder ratio and on temperature. Whilst some predictive models have been previously proposed to compute autogenous shrinkage [23,29], this study adopts the approach proposed in MC2010. Accordingly, the basic shrinkage e cbs is calculated as: where e cbs0 is the basic notional coefficient (which depends on the cement type and concrete compressive strength); and b bs ðtÞ is a time development function for basic shrinkage strain. The time development function for autogenous strain b bs (t) can be computed as follows: By differentiating Eq. (10) with respect to t, Eq. (9) can be rewritten as: Thus the change of internal humidity due to selfdesiccation can be calculated as: where all the variables are as defined above.
To validate the autogenous shrinkage model, FE analyses were carried out to predict the behavior of prisms with sealed and unsealed surfaces. The surface factor f was set to '0' in the analysis to simulate the effect of sealed surfaces, thus assuming that no moisture evaporated from the surface of the specimens and that only autogenous shrinkage occurred. Conversely, on the unsealed surfaces, f was set to 0.3 and 0.5 mm/day for mixes C1 and C2, respectively. Figure 4 shows that the FE models predict the results of the tested prisms with good accuracy, in particular for the sealed specimens. The differences between the numerical predictions and the experimental results can be attributed to (a) the use of 'average' environmental conditions and temperature in the FE modelling, and (b) the intrinsic inhomogeneity of concrete structure, which led to different conditions along the prisms.

FE modelling of hydro-mechanical behaviour
As discussed in Sect. 2, the calculation of shrinkage deformations requires: (1) analysing the time-dependent moisture distribution over the cross section of the element, and (2) determining the shrinkage camber and resultant deformation through section analysis. In this study, the hydro-mechanical behaviour of the I  beams was simulated using the diffusion model developed in Sects. 2 and 4, and a mechanical model (fibre-based section analysis [30]). These two models are linked by the hydro-shrinkage coefficient b (see Eqs. 1, 4 and 12). 2D diffusion elements were adopted to calculate the internal distribution of RH, whereas fibre elements were used to calculate time-dependent deformations. Both diffusion and mechanical models were discretised using quadrilateral elements with the same mesh size. The analyses to calculate the internal distribution of RH were performed using the general FEA package ADINA Ò [26], while the section analyses were carried out using routines programmed in Matlab Ò software [31]. Based on the RH measured in the laboratory, a constant ambient RH = 70% was used for the exposed surface in the diffusion model. The time-dependent internal RH due to self-desiccation was also considered in the FE modelling. The differences between the predictions and the test data can be attributed to the variation of environmental conditions (RH) and temperature during the tests, which cannot be captured by the FE models. This is confirmed by the smaller difference obtained between predictions and test results for the unsealed surface conditions (see Fig. 2b) when compared to partially sealed surface conditions (Fig. 2c). These results confirm that it is possible to develop appropriate numerical models to predict the shrinkage behaviour of the I beams with an acceptable accuracy. The results also indicate that the calibrated models can be used to perform further parametric analyses on the I beams. As expected, as the RH stabilises through the section, the difference in shrinkage between Top and Bottom slabs reduces, consequently reducing the shrinkage camber. The experimental and FE results presented in Fig. 2a-d confirm that the cross section shape can lead to shrinkage camber and that it influences the development of shrinkage-induced deformations at different depths of the specimen. The results also indicate that the numerical approach adopted in this study can predict the shrinkage camber and shrinkage strain of the tested specimens with good accuracy. However, due to the limited number of specimens tested in this study, the influence of different cross section geometry and ambient RH on shrinkage camber was not examined. These effects are examined in the following section, which presents the results of a parametric study that was carried out to analyse the shrinkage camber of I beams with different cross sections and ambient RH.

Influence of section depth
To determine how the section depth influences shrinkage camber (expressed as curvature), the height of the web was scaled up 2-5 times. This led to web depths of 150, 300, 450, 600 and 750 mm. The rest of the section dimensions were kept the same. The results in Fig. 5 show that the maximum value of shrinkage camber decreases as the height of the web increases. However, as the web thickness was the same, all peak values of camber are observed at around 30-40 days. Figure 6 compares results of shrinkage camber for different section sizes. In this case, the original section was scaled up geometrically by factors of 2-5 times, thus leading to sections of (flange width 9 total depth) 250 9 300, 500 9 600, 750 9 900, 1000 9 1200 and 1250 9 1500 mm. The results indicate that size affects both the maximum value of shrinkage camber and the time at which it is achieved. For example, the change in scale factor from 1 to 5  Figure 7 compares the response of the original model subjected to environmental RH varying from 40 to 80%. The results indicate that, as the RH increases, the maximum value of the shrinkage camber decreases from 0.95/mm (RH = 40%) to 0.32/mm (RH = 80%). It is also shown that the higher the RH, the earlier the camber reduces and tends to zero. Such trend is attributed to the internal moisture balancing with the external ambient RH.

Description of the bridge
The Yiju River Bridge at Chongqing, China (see Fig. 8) is selected as a case study to validate the hydromechanical FE approach proposed in this study. The two-span prestressed concrete box girder bridge was segmentally constructed over the Yiju River in 2016.   The RH inside the top slab of the box girder section was measured using a digital sensor with RH and temperature accuracies of ±2% and ±0.3°C, respectively. Three digital RH sensors were installed on the top slab of the mid-span section ''B'' at depths of 0, 30 and 100 mm from the top slab surface. The sensors were fixed to steel rebars, which in turn were welded to the girder's reinforcing cage. PVC tubes of 30 mm diameter protected the sensors during concrete casting. The data was collected using DC powered equipment.

FE validation and results
A total of 26 sectional models were used to simulate the 52 cantilever segments of the bridge in ADINA Ò software. Due to symmetry, only half of the bridge's cross section was modelled to reduce computational time. In the diffusion analysis, an average measured environment humidity H e = 0.70 was applied on the box girder surfaces as the convection boundary. The elastic modulus (E = 3.45 9 10 4 MPa) was obtained from tests on cubes, whereas the surface factor (f = 0.0005 m/day) and the hydro-shrinkage coefficient (b = 1.58 9 10 -3 ) were determined from a parametric analysis based on MC2010 and measured data.
The digital sensor on the box girder surface recorded ambient RH ranging from 51 to 100%, while the ambient temperature varied from 20 to 40°C during the measuring period (27 days). The short measuring period was due to the practical and financial constraints of the project. Figure 9a shows that the internal RH at depths of 30 and 100 mm decrease over time due to self-desiccation and drying of concrete. The measured internal RH at 30 mm depth decreases from 100 to 87% after 27 days. As expected, the RH at 30 mm is lower than that at 100 mm depth. The field data also show that the internal RH is significantly affected by the variation of ambient RH and temperature, as shown by the large scatter of data at 30 mm depth. Figure 9a also shows that the sectional diffusion analysis over the cross section of the bridge captures well the trend of the internal RH at both depths. Since the FE models use an average ambient RH, the calculated time-dependent RH curves are much smoother than the measured data. Figure 9b compares the calculated distribution of the internal RH at ages of 30, 2000 and 10,000 days (approximately 0, 5 and 27 years, respectively). It is shown that, within the first 30 days, the drying depth from the exposed surface is 50-100 mm and the internal RH at the center of the slab and web varies less than 1%. However, after 2000 days, the drying depth increases and some parts of the box section (e.g. thin top slabs) dry at a faster rate than the web. The internal RH of the thicker sections of the top and bottom slabs remains relatively high (up to 95%) even after 10,000 days, thus leading to a low shrinkage strain.
To study the long term shrinkage camber and shrinkage strain, the Yiju River bridge was divided into 26 sections in ADINA Ò . The sections had the same drying boundary conditions used in the previous analysis. The restrain of the bridge supports was considered by assigning suitable structural boundaries in the analysis of time-dependent deflection. As this study examines long periods of time, an average RH of 70% and an average temperature of 20°C were adopted in the FE models, which are deemed representative of future conditions at the Yiju River bridge. Dead load was neglected in the analysis. Figure 10a compares the development of shrinkage camber over time at sections ''A'', ''B'' and ''C'', which are representative of the bridge's sections at the pier, midspan and closure segments. The results confirm that shrinkage camber occurs at all sections due to differential drying across the box section, and that the magnitude of shrinkage camber depends on the size of the section. For example, the bottom slab thickness of section ''A'' is 1.92 m, or 2.5 times the slab thickness of section ''B'', and the maximum shrinkage camber at ''A'' is 6 9 10 -6 /m, i.e. 1.5 times the value at ''B'' (4 9 10 -6 /m). Figure 10a also shows that the  Fig. 10b show that the calculated shrinkage strain closely matches the MC 2010 predictions (Fig. 10b) when a value b ¼ 1:58 Â 10 À3 is used. Moreover, the relative error of shrinkage strain predictions and the MC 2010 models is around ±5%, which is considered sufficiently accurate. As expected, the results show that the time dependent axial shrinkage strain is significantly affected by the nominal thickness of the box girder section. The shrinkage strain at ''C'' calculated by the proposed approach is 381 le at 10,000 days, which is 1.68 times the value at ''A''. Figure 11 shows the calculated box girder deflection due to shrinkage camber at two points 'S1' and 'S2', which are approximately located at the free end and mid-span of the cantilever girder. As expected, the varying thickness of the web and bottom slab produces significant upward vertical displacements along the cantilever girder. This deflection increases over time up to 11,000 days, when it reaches maximum values of 3.0 cm and 1.8 cm at 'S1' and 'S2', respectively. After this peak point, the deflection decreases since the  majority of internal moisture in the top slab has been consumed by convection, the shrinkage at the top slab developed more slowly than at the bottom slab. According to the design, creep and prestressing losses will produce a 4.3 cm downward deflection in the bridge after 27 years of service. However, around 70% of such deflection can be potentially compensated by the shrinkage camber of the box girder section. It should be also noted that the influence of creep was disregarded in the analyses. Accordingly, current research by the authors is investigating the influence of creep on shrinkage-induced deformations, and their findings will be published in a forthcoming paper.
Overall, the results of this study indicate that the long term behaviour of long-span box girder bridges is affected by shrinkage camber, and therefore such additional curvatures need to be included in the analysis and design. Based on the results of this study, it can be concluded that the use of data from simple tests, along with inverse FE analyses implementing the adopted hydro-mechanical FE approach, can be used to simulate the shrinkage behaviour of concrete girders exposed to different ambient conditions. Due to the limited test data used for the calibration of the FE models adopted in this study, future research should focus on the determination of suitable surface factors and hydro-shrinkage coefficient values.

Conclusions
This article examined experimentally and analytically the effect of different geometries on the development of shrinkage camber in bridge box girders. Small-scale concrete specimens were subjected to different exposure conditions to measure the development of shrinkage. Based on the experimental results, a parametric analysis was performed to obtain the surface factor and the hydro-shrinkage coefficient. The test results were then used to examine analytically the effect of geometry and ambient relative humidity on the development of shrinkage deformation of the specimens using a hydro-mechanical approach. The proposed approach was then verified using the Yiju River Bridge (China) as a case study. Based on the experimental and analytical results presented in this article, the following conclusions are drawn:

Deflection(cm)
Time(days) S1 S2 Fig. 11 Calculated deflection of the cantilever girder of Yiju River Bridge 1. Test results from concrete beams with I cross section indicate that axial shortening (due to uniform shrinkage) and shrinkage camber (due to different cross section geometry) occurred simultaneously. The shrinkage camber increased initially up to a peak point, after which it decreased.
The results also indicate that the shrinkage camber can recover when the RH of the testing environment increases. 2. Based on test results and FE inverse analyses, surface factors f = 0.3-0.5 mm/day and a hydroshrinkage coefficients b = 0.0014-0.0015 were derived for typical concrete mixes used in bridge design/construction in China. The use of these values in a hydro-mechanical approach proved very effective at simulating the experimental shrinkage-induced bending deformations of the I beams subjected to different exposure conditions. 3. Results from parametric FE analyses confirm that section size and ambient humidity are the main factors affecting the development of shrinkage deformations in concrete girders. Moreover, the section size affects both the maximum value of shrinkage camber and the time at which it is achieved. However, only the thickness of the cross section components had a major influence on the response after the maximum camber was achieved. 4. The proposed hydro-mechanical FE approach provides a reasonable estimate of the RH of sections of the Yiju River Bridge. Based on these estimates, it is shown that in the long term (27 years), approximately 70% of the vertical deflection can be potentially compensated by the shrinkage camber of the box girder section (excluding creep effects). 5. The results of this study confirm that the long term behaviour of long-span box girder bridges is affected by shrinkage camber, and therefore such additional curvatures need to be included during initial analysis and subsequent assessment.