Research Article  Open Access
Shiwei Li, Yongqing Yang, Qianhui Pu, Wangqing Wen, Aiguo Yan, "ThermoHydroMechanical Combined Effect Analysis Model for EarlyAge Concrete Bridges and Its Application", Advances in Civil Engineering, vol. 2020, Article ID 8864109, 15 pages, 2020. https://doi.org/10.1155/2020/8864109
ThermoHydroMechanical Combined Effect Analysis Model for EarlyAge Concrete Bridges and Its Application
Abstract
The early cracking of concrete beam bridges remains a concern in civil engineering. An analytical model considering the combined effect of thermohydromechanical processes forms the basis for assessing the cracking risk of girders during construction. Based on the equivalent hydration theory, the temperature and moisture conduction processes and the evolution of the mechanical properties of concrete were modeled as a function of the equivalent age. A coupling model for the temperature and moisture fields was established, and a theoretical framework for analyzing the thermohydromechanical combined effect was presented. Based on this, a numerical analysis method was proposed and implemented into ABAQUS; the results were validated with some typical tests. Finally, a longspan prestressed concrete (PC) box girder bridge with balanced cantilever construction was taken as an example, and the causes of web cracking and its impact degree were analyzed. The results show that the rate of moisture conduction is significantly lower than the rate of temperature conduction; even for thinwalled components, there exists a significant humidity gradient on the surface layer. The humidityinduced shrinkage and restraint of the precast members are the main causes of web cracking.
1. Introduction
During the setting and hardening period, concrete is subjected to the combined effect of hydration heat, ambient temperature and humidity, shrinkage and creep, and selfweight, that is, the combined effect of thermohydromechanical processes [1 ,2]. This leads to elastic and inelastic deformations. If these deformations are constrained by internal or external conditions, an additional stress will be developed in the concrete. Moreover, the mechanical properties of concrete materials vary with time. When the combined effect exceeds the resistance of the materials, irreversible damage occurs [3, 4].
When the damage due to the combined effect of thermohydromechanical processes accumulates to a certain extent, earlyage cracks are induced in the concrete. The earlyage cracking of concrete bridges is a universal phenomenon. Subramaniam [5] investigated the cracking of concrete bridge decks in New York and found that a large number of cracks occurred within 48 h of concrete pouring. Li [6] investigated eight concrete girder bridges under construction in Northeast China and found that the main girder concrete cracked to different degrees during the hardening process. Any irreversible damage to concrete bridges in the process of completion can significantly affect the construction quality and structural durability.
The focus in practical applications is on the prevention of nonstructural cracks. However, the preventive measures are ineffective because of the poor understanding of the cracking mechanism. The cracking phenomenon remains to be quite frequent during construction [7]. The role of single factors in the formation of cracks is often emphasized. Li [6] focused on the mechanism of cracking induced by shrinkage. Huang et al. [8] attributed the cracking of concrete girders to excessive early hydration heat. For concrete bridges, the concrete mix ratio is varied, the ambient environment is complex and dynamic, and the temperature and humidity conduction processes in concrete have coupling characteristics [9, 10]. Therefore, it is necessary to establish a suitable method to analyze the combined effect of the thermohydromechanical processes.
To analyze this combined effect during concrete hardening, much research has been devoted to this subject in recent decades. Ulm and Coussy [11, 12] proposed a thermochemomechanical model for the prediction of deformation and cracking accounting for strength growth through chemoplastic coupling within the theory of elastoplasticity, but the effects of creep and relative humidity were not considered. Cervera et al. [13, 14] formulated a coupled thermochemomechanical model for the behavior of concrete at early ages, but the effect of the relative humidity is not considered. Lackner and Mang [15] proposed a chemomechanical model for earlyage cracking of concrete in which the effect of relative humidity is also neglected. Gawin et al. [16, 17] developed a hydrothermalchemomechanical model of concrete at early ages which comprehensively accounts for the aging and hydrothermal coupling but lacks a comprehensive constitutive law for earlyage concrete. Di Luzio and Cusatis [18] proposed the solidificationmicroprestressmicroplane (SMM) model describing the mechanical behavior of concrete at early age under variable hydrothermal conditions. Cui et al. [19] formulated a sequential coupling analysis method for concrete members; however, those models cannot be applied to the sequential casting of multisegment concrete bridges, and they do not consider the restraint of steel bars on shrinkage deformation and its contribution to crack resistance.
In this study, to accurately analyze the thermohydromechanical combined effect in the construction process of concrete bridges, the temperature and moisture conduction processes and the evolution of the mechanical properties of concrete were modeled as a function of the equivalent age based on the equivalent hydration theory. A coupling model for the temperature and moisture fields was established, and a theoretical framework for the combined effect analysis was presented. The corresponding numerical analysis method was proposed and implemented into ABAQUS; the results were validated by conducting some typical tests. Finally, a longspan prestressed concrete (PC) box girder bridge with balanced cantilever construction was taken as an example, and the causes of web cracking and its impact degree were analyzed.
2. Coupling Analysis Theory of Temperature and Moisture Fields
2.1. Theory of Temperature Conduction
Due to the combined effect of hydration heat and ambient environment (atmospheric temperature and solar radiation), concrete remains in a dynamic thermal equilibrium with the external environment. During concrete hardening, the temperature field of concrete is mainly driven by hydration heat. When the hydration reaction is largely completed, the temperature field is mainly driven by the ambient temperature. Because of the influence of temperature conduction rate, the response of the internal temperature of concrete lags behind the change in the ambient temperature.
2.1.1. Hydration Heat
The active chemical components in the cement clinker react with the mixing water, releasing a lot of heat. The total heat release (H_{total}) considering the influence of clinker composition can be calculated as [20]where, C_{3}S, C_{2}S, C_{3}A, and C_{4}AF are the mass ratios of each chemical component; free CaO is the mass ratio of free calcium oxide in the cement; slag indicates blast furnace slag; FA is the mass ratio of fly ash; and p_{CaO} is the proportion ratio of CaO in FA.
When the active chemical components in cement clinker cannot be accurately obtained, H_{total} can be calculated as [21]where W is the amount of cement; F is the amount of admixture; k is the coefficient of reduction and equals 0.25 for fly ash; Q_{0} is a constant for heat of hydration of cement and equals 330 kJ/m^{3} for silicate cement 425# or 525#.
The hydration heat release (H(t)) at different times can be expressed as a function of the equivalent age (t_{e}), as shown inwhere α(t_{e}) is a function of the hydration degree and is related to t_{e}. It can be expressed aswhere η is the hydration time parameter and β is the hydration rate parameter, both of which are related to the clinker composition, cement fineness, and admixture; α_{u} is the ultimate degree of hydration; η, β, and α_{u} are obtained bywhere Blaine is the fineness coefficient of cement (m^{2}·kg^{−1}); FA_CaO is the mass ratio of CaO to FA; is the watercement ratio; HI is the hydraulic index and can be expressed aswhere CaO, MgO, Al_{2}O_{3}, and SiO_{2} are the mass ratios of each chemical component with respect to the admixtures.
t_{e} considering the temperature history is given bywhere E_{A} is the activation energy; for T ≥ 20°C, E_{A} = 33.5 kJ·mol^{−1}, and T < 20°C, E_{A} = 33.5 + 1.47 × (293 − T) (kJ·mol^{−1}); R is the universal gas constant; and T_{ref} is the reference temperature, which is typically 20°C.
2.1.2. Temperature Conduction
According to the principle of thermal equilibrium, the heat absorbed due to an increase in the structural temperature is equal to the sum of internal hydration heat and external inflow heat [22]. According to Fourier heat conduction theory, the differential equation for concrete temperature conduction can be expressed as [20]where k is the current thermal conductivity; k_{u} is the thermal conductivity of hardened concrete; c is the specific heat of concrete; h is the relative humidity; ξ is the influence coefficient of humidity on temperature conduction; ρ is the density of concrete; W_{c}, W_{a}, and are the amounts of cement, aggregate, and water, respectively, by weight; c_{c}, c_{a}, and are the specific heats of aggregate, cement, and water, respectively; c_{cef} is the fictitious specific heat of the hydrated part of cement and is equal to (8.4 T_{c} + 339) kJ·kg^{−1}.
The thermal conductivity k and specific heat of concrete c are related to the temperature of concrete, while the temperature of concrete is related to its equivalent age t_{e}. The equivalent age t_{e} comprehensively reflects the influence of the hydration process and the environmental temperature of concrete. Both k and c are described as functions of t_{e}, as shown in equation (8). Thus, the key performance parameters of the concrete temperature field evolve with the hydration process.
2.1.3. Initial and Boundary Conditions
The molding temperature of concrete is generally considered the initial condition of temperature conduction.
Based on the different contact conditions between the concrete surface and surrounding media, the boundary conditions can be divided into four types [21]. The heat transfer between concrete and air is considered based on the third type; this boundary condition can consider the influence of the type and thickness of the mold on the heat transfer of concrete, as shown in equation (12) [20]:where χ_{s} is the equivalent surface conduction of the concrete surface; χ is the surface conductance; and h_{i} and λ_{i} are the thickness and conductivity of the ith layer, respectively.
The heat transfer between the successive segments is considered based on the fourth type of the boundary conditions.
2.2. Theory of Moisture Conduction
Evaporation to the environment and cement hydration reaction are the two main factors leading to the reduction in the moisture content of concrete [23].
2.2.1. Humidity Reduction due to Cement Hydration
Zhang et al. [23] proposed a twostage calculation model for humidity reduction due to cement hydration based on a humidity conduction test. The model is expressed inwhere H_{su} is the relative humidity (RH) considering selfdesiccation at ultimate degree hydration; α_{c} is the critical hydration degree at which the humidity inside the concrete starts to decrease from 100% level; and κ is the drying rate parameter.
2.2.2. Moisture Conduction
According to the second Fick’s law, the differential equation for concrete moisture conduction can be expressed as [23]where H_{s} is the humidity reduction due to cement hydration and D(h) is the moisture transfer coefﬁcient. Experimental results have shown that D(h) increases rapidly with an increase in the humidity and that the temperature increases the humidity conduction rate, as shown in equation (15) [24]:where D_{1} is the maximum D(h) for h = 100%, as shown in equation (16); ω = D_{0}/D_{1}, D_{0} is minimum D(h) for h = 0%; h_{c} is the relative pore humidity at 0.5 D_{1}; n is an exponent; f(t_{e}) is the temperature influence function, as shown in equation (17).where D_{1.0} = 1 × 10^{−8} m^{2}/s and f_{cm} is the mean compressive strength of concrete expressed in MPa.
where U_{aT} and U_{ar} are the apparent activation energies at temperature T and at the reference temperature T_{ref}, respectively.
2.2.3. Initial and Boundary Conditions
The initial condition for moisture conduction in poured concrete is generally 100%.
The humidity transfer rate at the concrete surface is related to ambient humidity, wind speed, and temperature. It can be calculated using equation (18) [25]:where f_{boundary} is a proportionality factor and h_{surf} and h_{en} are the concrete surface humidity and environmental humidity, respectively.
3. Analysis of Stress Field
In addition to the corresponding selfstresses due to the nonuniform temperature and humidity gradients, an additional tensile stress will occur when the shrinkage deformation due to the temperature and humidity fields is constrained. Moreover, creep will cause stress relaxation. Remarkably, when concrete is subjected to the combined effect of thermohydromechanical processes, its mechanical properties evolve with time.
3.1. Mechanical Parameters of EarlyAge Concrete
The elastic modulus E, tensile strengths f_{t}, and compressive strength f_{c} of concrete varying with the equivalent age can be expressed as [26]where E_{28}, f_{t28}, and f_{c28} are the elastic modulus (GPa), tensile strength (MPa), and compressive strength (MPa) of concrete at 28 d, respectively; t_{0} is the time at which the concrete strength begins to develop (in h); and the parameters t_{0}, s, n_{E}, and n_{t} can be taken as 10 h, 0.211, 0.339, and 0.621, respectively, when there is no reliable test data.
3.2. Relationship between Temperature Field Change and Temperature Deformation
Temperature deformation can be expressed as a product of the temperature deformation coefficient and the temperature difference, as shown inwhere α_{T} is the temperature deformation coefficient, which is in the range of 1.0–1.2 × 10^{−5}/^{°}C, and ΔT is the temperature difference (in °C).
3.3. Relationship between Humidity Field Change and Shrinkage Deformation
When the relative humidity in the atmosphere is lower than that in concrete, the moisture in the concrete will diffuse to the environment and cause drying shrinkage. The hydration reaction of cement consumes free water in the capillaries of concrete, thus decreasing the interior relative humidity and causing an autogenous shrinkage of concrete [27]. The relationship between moistureinduced shrinkage and interior RH from the time of concrete setting can be expressed as [27]where ε_{W} is the shrinkage strain, ε_{c} is the shrinkage during the humidity developing stage I (RH = 100%), which is principally governed by w/c, and α_{h} is the shrinkage strain due to unit humidity reduction.
3.4. Tensile Creep of EarlyAge Concrete
The double power law (DPL) creep model proposed by Bazant and Panula [28] is selected because it has shown good accuracy in many earlyage tensile creep tests [29, 30]. The model is expressed inwhere J(t, τ) is the creep compliance function at time t, τ is the loading age, E(τ) is the asymptotic elastic modulus, and q, d, and p are the material parameters, calibrated based on experimental concrete creep tests.
3.5. Strength Criterion with Age Evolution
In this study, the Lubliner criterion [31], which has been widely used for concrete in complex stress states, was chosen, as shown inwhere < > are Macaulay brackets and < x > is taken as ; I1 is the first stress invariant; J_{2} represents the second deviatoric stress invariant; is the maximum effective stress eigenvalue; α and β are undetermined parameters determined by the uniaxial tensile strength and compressive strength, respectively, and can be taken as 0.12 and 7.68, respectively; γ is equal to 3 in the low confining pressure stress state.
4. Numerical Analysis Method of ThermoHydroMechanical Combined Effect
To realize the analysis of the thermohydromechanical combined effect, the temperature and humidity fields of concrete are determined first, and the additional deformation is then applied to the concrete structure. However, because of the coupling characteristics of the temperature and moisture conduction, iteration at each time step will increase the computing costs significantly. Therefore, the whole iteration method is used. This method does not significantly increase the calculation cost and shows sufficient calculation accuracy. The processes can be described as follows.
First, the initial temperature field is analyzed under the specified humidity environment. According to the environmental characteristics and mix proportion of concrete, the initial and boundary conditions of the temperature field analysis and the relevant basic parameters (c and k) are determined. The temperature field changes at each time step are calculated through the differential equation of the concrete temperature conduction, and the temperature change data at all the nodes in the time history are stored.
The initial humidity field is then analyzed. In addition to determining the initial conditions, boundary conditions, and related calculation parameters (D(h)), the results of the temperature field analysis are imported to calculate the humidity conductivity and solve the humidity conduction differential equation. Similarly, the time history humidity data at all the nodes need to be stored.
Finally, the calculation results of the humidity field are imported to the process of temperature field analysis, and a secondary temperature field analysis is carried out.
The above process is a global iteration, and the desired accuracy is achieved after several iterations.
When the temperature and humidity fields of concrete are determined, the stress field can be analyzed. The temperature and humidity data at each time step are used to calculate the additional deformation in terms of the temperature expansion coefficient and humidity deformation coefficient and the mechanical parameters, which vary with the effective age, are also adjusted.
The finite element software ABAQUS and its secondary development platform were selected for this complex process. The core programming steps are described as follows.
Because the differential equations of the temperature and moisture conduction (equations (7) and (13)) have a consistent form, they can be solved using subroutine UMATHT. The boundary conditions of the temperature and humidity conduction are defined by subroutine FILM, which can simulate the change in the temperature and humidity environment with time. The subroutine UEXPAN is used to extract the temperature and humidity data at each time step and calculate the corresponding increment strain. Moreover, the subroutine USDFLD is used to consider the influence of creep. In addition, the evolution of the mechanical parameters is simulated by defining field variables (FIELD). Figure 1 shows the specific process.
5. Model Verification
5.1. Temperature Field Simulation of a Prefabricated Box Girder by Electric Heating
A prefabricated concrete box girder bridge of the ChongqingGuizhou Railway was constructed in winter. To accelerate construction, measures of electric heating and insulation layer were employed. To understand the actual effect of these measures, Yao et al. [32] conducted a temperature test.
The length of the test beam was 26.58 m. After 6 h of casting, the inner form was removed, and four electric heaters with a power of 2000 W were evenly installed along the span direction. At the beamend of the box girder, the retaining wall with masonry was overlapped by a polyethylene foam layer. After 24 h of casting, the outer form was replaced by the polyethylene foam layer. The electric heaters were removed after working for three days; the insulation layer remained overlapped to avoid fast cooling. Table 1 lists the mix proportion of concrete. Figures 2 and 3 show the heating device box girder section and temperature measurement points.

The subroutine UMATHT was used for the temperature field analysis. The concrete was modeled using DC3D8. Based on the field measurements, the hydration heat function can be calculated using (26), and the unit of t_{e} is h:
The concrete molding temperature was 20°C. The ambient temperature, which changed only slightly, was approximately taken as 5°C. Based on the measured data, the electric heating inside the box girder is considered by constant thermal radiation (32°C). The equivalent surface conductance of the outer surface of concrete between concrete and thermal insulation layers is χ_{s} = 2 W/(m^{2}°C) and that of the flange end is χ_{s} = 0.5 W/(m^{2}·°C). The equivalent surface conductance between concrete and air is χ_{s} = 15 W/(m^{2}·°C). The heat conductivity k_{u} is calculated using equation (9) and taken as 2.618 W/(m^{2}·°C). The specific heat capacity c is calculated using equation (10) and taken as 913.22 J/(kg·°C). The density of concrete is 2497 kg·m^{−3}.
The results of the temperature measurement at A3 and A4 (as shown in Figure 3) are compared with the simulation results shown in Figure 4; the results are in good agreement, regardless of the peak temperature or the change trend. This demonstrates that the proposed method can effectively simulate the temperature field of concrete bridges.
5.2. Simulation of Moisture Conduction Process and Shrinkage Deformation
Zhang et al. [33] conducted a moisture conduction test on earlyage concrete under two drying conditions: only surface sealing (specimen I) and surface sealing followed by exposure to the environment (specimen II). Table 2 lists the concrete mix proportion. The specimens were cuboidshaped and 400 mm × 100 mm × 60 mm in size. A layer of 1 mm thick polytetrafluoroethylene (PTFE) sheet was laid on the bottom of the specimen to reduce the friction resistance.

Figure 5 shows the test apparatus. The digital temperature and humidity sensors, which are located at the geometric center of the specimen, were used. The free deformation of the concrete was measured using LVDT.
When the density and specific heat of concrete are set to 1, the differential equation of the moisture conduction has the same mathematical form as that of the temperature conduction. Therefore, the moisture field analysis method is consistent with the temperature field analysis method. The concrete was also modeled using DC3D8 with 2300 elements. The ambient temperature and humidity vary only slightly and are taken as 23°C and 16%, respectively. Table 3 lists the other relevant calculation parameters. The deformation due to the change in the humidity is calculated using equation (23); ε_{c} is 210 × 10^{−6}; and α_{h} is 15 × 10^{−6}/%.

Figures 6 and 7 show the humidity change and free deformation of the two specimens. The simulation and actual results are mostly consistent. Figure 8 shows the development of the humidity gradient of specimen I; Figure 9 shows a comparison between the temperature and humidity conduction processes. The results show that the humidity conduction rate is much lower than the temperature conduction rate; even for the thinwalled components, there remains a significant humidity gradient beneath the surface.
5.3. Simulation of Tensile Creep in EarlyAge Concrete
To understand the role of tensile creep in earlyage concrete, Khan et al. [34] evaluated the creep deformation of concrete dogbone specimens under a sustained tensile stress. Figure 10 shows the specimen size and loading apparatus. Table 4 lists the mix proportion of concrete; Table 5 lists the main mechanical properties.
(a)
(b)


The loading ages for the dogbone creep specimens were 1 and 2 d. The loading stress level was 50% of the measured tensile strength at the loading age, that is, 0.7 MPa on day 1 of aging and 1.05 MPa on day 2 of aging.
Two user subroutines (USDFLD and UEXPAN) were required to implement the creep analysis through the ABAQUS secondary development platform. The DPL creep function was chosen for this analysis, as shown in equation (24); q, d, and p are 0.59, 0.26, and 0.41, respectively. Figure 11 shows the simulation and test results of the creep coefficient. The simulated values were closer to the test results.
6. EarlyAge Cracking Analysis of a PC Box Girder Bridge with Cantilever Casting
6.1. Engineering Background
A PC box girder bridge with a span of 110 + 190+100 m was segmentally constructed in 2009. The web thickness varies from 1 m at the 0# block to 0.7 m at the midspan. Visible cracks occurred in 2# block during construction. Wang and Fang evaluated the change in the web temperature of 2# block during hardening in situ [35].
Table 6 lists the mix proportion of the main girder of C60. Figure 12 shows the evolution of the tensile strength and elastic modulus. The diameters of the vertical and transverse steel bars are Φ20 and Φ12, respectively. Figure 13 shows the temperature measuring points in the 2# block web.

6.2. Calculation Model
The shrinkage deformation of 2# block is constrained by 1# block. Therefore, to accurately analyze the response of 2# block under the combined action of multifields, the influence of 1# block should be considered. A 1/2 model of 1# block and 2# block with solid elements was constructed considering the symmetry of the structure. The constraint effects on the free deformation and contribution to concrete tension of the steel bar are considered, as shown in Figure 13. Similarly, the subroutines UMATHT and FILM are used to analyze the temperature and humidity conduction processes, and the temperature load and creep correction are realized using the subroutines UEXPAN and USDFLD.
For the temperature and humidity field analysis, the DC3D8 element is used to model the concrete. For the stress field analysis, the C3D8 concrete element and T3D4 rebar element are used to model the concrete girder and steel bar, respectively. The foregoing analysis results show that a significant humidity gradient is formed beneath the surface. Therefore, the thickness of the concrete element near the surface should not be too large, and the mesh should be consistent due to the continuity of the thermohydromechanical analysis. There are 193596 concrete elements and 19800 rebar elements, as shown in Figure 14. The active and inactive element methods are used to simulate the construction process. The changes in the ambient temperature and humidity are simulated using the subroutine FILM.
(a)
(b)
The hydration heat function is calculated using (27). The specific heat capacity of the concrete with the change in the effective age is calculated using (10). For the temperature boundary conditions, the ambient temperature is fitted with the sum of several trigonometric functions to the measured value, as shown in (28). When n is taken as 7, the prediction accuracy is good, as shown in Figure 15. k_{u} and χ_{s} are taken as 2.6 and 15 W/(m^{2}·°C), respectively. The density of concrete is 2459.7 kg·m^{−3}.
The humidity reduction due to cement hydration reaction is calculated using equation (12). The humidity conductivity coefficient is calculated using the CEBFIP (2010) formula considering temperature correction, as shown in equation (15). For the humidity boundary conditions, the main girder concrete should be kept wet during hardening. Therefore, the external humidity environment of the concrete is affected by the level of artificial intervention. To explore the influence of different curing effects on the humidity field, three calculation conditions (60, 80, and 100%) are defined based on the varying environmental humidity. According to the research results obtained by Oliveira et al. [25], the surface factor f_{boundary} is approximately linear with D_{1}, as shown in equation (29):
For the stress field analysis, the evolutions of the tensile strength and elastic modulus can be calculated using equations (19) and (20). The stress relaxation due to the creep is calculated using (24). q, d, and p are taken as 0.6, 0.2, and 0.2, respectively. The Lubliner criterion [31] is chosen to judge whether the concrete is cracked, as shown in (25).
6.3. Computational Results
6.3.1. Analysis of Temperature Field Calculation Results
Figure 15 shows a comparison between the measured data and the simulation results at measuring points 1 and 3. The simulated and measured results are consistent, indicating that the temperature field of the bridge is reflected in the simulation.
6.3.2. Analysis of Moisture Field Calculation Results
Figure 16 shows the humidity simulation results at points 1 and 3 under various calculation conditions. The results show that the humidity at the web center (point 3) is constant, mainly because the humidity conduction speed is too low, and the decrease in the humidity is mainly due to the cement hydration reaction. However, the humidity evolution on the web surface (point 1) is different, and the decrease in the humidity is mainly due to the humidity exchange with the environment. The lower the ambient humidity is, the higher the humidity gradient along the thickness of the web increases with time.
6.3.3. Analysis of Stress Field Calculation Results
Figure 17 shows the web cracking simulation results under three calculation conditions; the web cracks are colored in red. The results show the following:(1)The cracks first appear at the interface, mainly because the free deformation of 2# block is restrained by 1# block.(2)The degree of web cracking under different humidity environments is different, indicating that the external humidity environment has a significant impact on the mechanical behavior of the structure.
(a)
(b)
(c)
Based on whether the creep effect is considered or not, Figures 18 and 19 show the calculated values of the principal tensile stress at points 1 and 3 (in Figure 13). The results show that the creep can cause stress relaxation, which improves the crack resistance of the web. When considering the stress relaxation due to the creep and the crack resistance of the web reinforcement, the thermal moisture coupling effect does not cause any cracks at the center of the web (Figure 17).
Figure 20 shows the calculated values of the stress due to temperature and humidity at measuring points 1 and 3.
For the temperature stress (in Figure 20),
(1)The peak value is reached in approximately 33 h, because the temperature difference between the inside and outside of the web reaches the extreme value. However, the elastic modulus of the concrete is relatively low (approximately 14 GPa), and the surface tensile stress due to the temperature gradient is approximately 0.9 MPa;(2)With the completion of hydration (at 95 h), the temperature field of the web is driven by the ambient temperature, and the temperature stress fluctuates with the change in the ambient temperature.For humidity stress, (1)Before the removal of the mold, the cement hydration degree does not reach the critical hydration degree α_{c}, and the internal humidity of the web remains in the first stage. In this stage, the humidity change law is not influenced by the external environment humidity; therefore, the humidity stress under the different calculation conditions has the same evolution law;(2)After removing the mold (at 72 h), the surface layer of the web exchanges humidity with the environment, and the humidity gradient on the internal and external surfaces of the web increases with the age. When the environmental humidity is lower, the moisture conduction on the surface layer of the web is faster; however, the humidity at the center of the web is the same. Therefore, the humidity gradient along the thickness direction is higher, and the nonlinear distribution of the humidity gradient will make the web surface experience tension, and the center of the web is affected. With the decrease in the external humidity, the surface tensile stress increases, and the central tensile stress decreases. The calculation results, shown in Figure 19, confirm this rule;(3)When measures to preserve moisture are ineffective, the stress due to the humidity field will be significantly greater than that due to the temperature field.
7. Conclusion
A new method is proposed to analyze the thermohydromechanical combined effect during the construction of concrete beam bridges. A comparison between the calculation and test results shows that the proposed method can effectively simulate the combined effect of multifields under complex construction environments.
The humidity conduction rate is much lower than humidity conduction rate; even for thinwalled components, there remains a significant humidity gradient on the surface layer. When measures to avoid moisture are ineffective, the stress due to the humidity field will be significantly greater than that due to the temperature field; therefore, more attention should be paid to the stages of design and construction maintenance.
For webs and other thinwalled components, the combined effect of humidityinduced shrinkage and the constraint on the early casted segment lead to cracks in the web.
With regard to the elastoplastic mechanics of concrete, the concrete in the crack area will produce not only creep deformation but also plastic deformation. It is necessary to establish a viscoplastic constitutive relationship that evolves with the concrete age; this requires further study.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
This work was funded by China Postdoctoral Science Foundation (2020M672460) and the Major Scientific Research Projects of China Railway Siyuan Survey and Design Group Company Limited (2019K130).
References
 S. W. Li, “Study on timedepending structural behavior of longspan cantilever casting PC box girder bridges based on thermo–hydro–mechanical combined effects analysis,” Southwest Jiaotong University, Chengdu, China, 2019, Doctoral thesis. View at: Google Scholar
 S. W. Li, Q. Y. Yang, Q. H. Qian, D. Yang, B. L. Sun, and X. B. Li, “Threedimensional nonlinear creep and shrinkage effects of a longspan prestressed concrete box girder bridge,” Structural Concrete, vol. 20, no. 2, pp. 638–649, 2019. View at: Publisher Site  Google Scholar
 A. Elsafty and A. AbdelMohti, “Investigation of likelihood of cracking in reinforced concrete bridge decks,” International Journal of Concrete Structures and Materials, vol. 7, no. 1, pp. 79–93, 2013. View at: Publisher Site  Google Scholar
 V. Gribniak, G. Kaklauskas, D. Čygas, D. Bačinskas, R. Kupliauskas, and A. Sokolov, “Investigation of concrete cracking effect in deck slab of continuous bridges,” The Baltic Journal of Road and Bridge Engineering, vol. 5, no. 2, pp. 83–88, 2010. View at: Publisher Site  Google Scholar
 K. V. L. Subramaniam, “Identification of earlyage cracking in concrete bridge decks,” Journal of Performance of Constructed Facilities, vol. 30, no. 6, pp. 1–9, 2016. View at: Publisher Site  Google Scholar
 G. D. Li, “Mechanism and control for earlyage shrinkage crack of HPC used for bridge engineering,” Harbin Institute of Technology, Harbin, China, 2014, Doctoral thesis. View at: Google Scholar
 T. M. Wang, Control of Cracking in Engineering Structure, China Architecture & Building Press, Beijing, China, 2nd edition, 2017.
 Y. Huang, G. Liu, S. Huang, R. Rao, and C. Hu, “Experimental and finite element investigations on the temperature field of a massive bridge pier caused by the hydration heat of concrete,” Construction and Building Materials, vol. 192, pp. 240–252, 2018. View at: Publisher Site  Google Scholar
 X. Hu, C. Shi, Z. Shi, B. Tong, and D. Wang, “Early age shrinkage and heat of hydration of cementfly ashslag ternary blends,” Construction and Building Materials, vol. 153, pp. 857–865, 2017. View at: Publisher Site  Google Scholar
 S. E Giuseppe, “Thermohygrochemomechanical model of concrete at early ages and its extension to tumor growth numerical analysis,” Università Degli Studi Di Padova, Cachan, France, 2013, Doctoral thesis. View at: Google Scholar
 F.J. Ulm and O. Coussy, “Modeling of thermochemomechanical couplings of concrete at early ages,” Journal of Engineering Mechanics, vol. 121, no. 7, pp. 785–794, 1995. View at: Publisher Site  Google Scholar
 F.J. Ulm and O. Coussy, “Strength growth as chemoplastic hardening in early age concrete,” Journal of Engineering Mechanics, vol. 122, no. 12, pp. 1123–1132, 1996. View at: Publisher Site  Google Scholar
 M. Cervera, J. Oliver, and T. Prato, “Thermochemomechanical model for concrete. I: hydration and aging,” Journal of Engineering Mechanics, vol. 125, no. 9, pp. 1018–1027, 1999. View at: Publisher Site  Google Scholar
 M. Cervera, J. Oliver, and T. Prato, “Thermochemomechanical model for concrete. II: damage and creep,” Journal of Engineering Mechanics, vol. 125, no. 9, pp. 1028–1039, 1999. View at: Publisher Site  Google Scholar
 R. Lackner and H. A. Mang, “Chemoplastic material model for the simulation of earlyage cracking: from the constitutive law to numerical analyses of massive concrete structures,” Cement and Concrete Composites, vol. 26, no. 5, pp. 551–562, 2004. View at: Publisher Site  Google Scholar
 D. Gawin, B. A. F. Pesavento, and B. A. Schrefler, “Hygrothermochemomechanical modelling of concrete at early ages and beyond. Part I: hydration and hygrothermal phenomena,” International Journal for Numerical Methods in Engineering, vol. 67, no. 3, pp. 299–331, 2006. View at: Publisher Site  Google Scholar
 D. Gawin, B. A. F. Pesavento, and B. A. Schrefler, “Hygrothermochemomechanical modelling of concrete at early ages and beyond. Part II: shrinkage and creep of concrete,” International Journal for Numerical Methods in Engineering, vol. 67, no. 3, pp. 332–363, 2006. View at: Publisher Site  Google Scholar
 G. Di Luzio and G. Cusatis, “Solidificationmicroprestressmicroplane (SMM) theory for concrete at early age: theory, validation and application,” International Journal of Solids and Structures, vol. 50, no. 6, pp. 957–975, 2013. View at: Publisher Site  Google Scholar
 W. Cui, W. Chen, and N. Wang, “Thermo–hydro–mechanical coupling analysis of earlyage concrete with behavioral changes considered and its application,” China Civil Engineering Journal, vol. 48, no. 2, pp. 44–53, 2015. View at: Google Scholar
 Z. Ge, Predicting Temperature and Strength Development of the Field Concrete, Iowa State University, Ames, IA, USA, 2005, Doctoral thesis.
 B. F. Zhu, Thermal Stresses and Temperature Control of Mass Concrete, Tsinghua University Press, Beijing, China, 2012.
 Y. Lee, M.S. Choi, S.T. Yi, and J.K. Kim, “Experimental study on the convective heat transfer coefficient of earlyage concrete,” Cement and Concrete Composites, vol. 31, no. 1, pp. 60–71, 2009. View at: Publisher Site  Google Scholar
 J. Zhang, K. Qi, and Y. Huang, “Calculation of moisture distribution in earlyage concrete,” Journal of Engineering Mechanics, vol. 135, no. 8, pp. 871–880, 2009. View at: Publisher Site  Google Scholar
 CEBFIP, CEBFIP Model Code 2010, International Federation for Structural Concrete (fib), Lausanne, Switzerland, 2012.
 M. A. Oliveira, M. Azenha, and P. B. Lourenço, “Simulation of humidity fields in concrete: experimental validation and parameter estimation,” Journal of Advanced Concrete Technology, vol. 13, no. 4, pp. 214–229, 2015. View at: Publisher Site  Google Scholar
 T. Kanstad, T. A. Hammer, Ø. Bjøntegaard, and E. J. Sellevold, “Mechanical properties of young concrete: Part II: determination of model parameters and test program proposals,” Materials and Structures, vol. 36, no. 4, pp. 226–230, 2003. View at: Publisher Site  Google Scholar
 J. Zhang, H. Dongwei, and S. Wei, “Experimental study on the relationship between shrinkage and interior humidity of concrete at early age,” Magazine of Concrete Research, vol. 62, no. 3, pp. 191–219, 2010. View at: Publisher Site  Google Scholar
 Z. P. Bazant and L. Panula, “Practical prediction of time dependent deformation of concrete, Parts I and II,” Mater. Struct., vol. 11, no. 65, pp. 307–328, 1978. View at: Google Scholar
 G. M. Ji, “Cracking risk of concrete structures in the hardening phase: experiments, material modeling and finite element analysis,” The Norwegian University of Science and Technology, Trondheim, Norway, 2008, Doctoral thesis. View at: Google Scholar
 M. Azenha, L. Leitão, J. L. Granja, C. de Sousa, R. Faria, and J. A. O. Barros, “Experimental validation of a framework for hygromechanical simulation of selfinduced stresses in concrete,” Cement and Concrete Composites, vol. 80, pp. 41–54, 2017. View at: Publisher Site  Google Scholar
 J. Lubliner, J. Oliver, S. Oller, and E. Oñate, “A plasticdamage model for concrete,” International Journal of Solids and Structures, vol. 25, no. 3, pp. 299–326, 1989. View at: Publisher Site  Google Scholar
 G. Yao, Z. J. Yu, Y. Yang, and W. Jia, “Study on temperature field and temperature stress of singlecell box girder cured with electric heater,” China Civil Engineering Journal, vol. 51, no. 3, pp. 109–114, 2018. View at: Google Scholar
 J. Zhang, D. W. Hou, and Y. Gao, “Uniform driving force for autogenous and drying shrinkage of concrete,” Journal of Tsinghua University (Science and Technology), vol. 50, no. 9, pp. 1321–1324, 2010. View at: Google Scholar
 I. Khan, A. Castel, and R. I. Gilbert, “Tensile creep and earlyage concrete cracking due to restrained shrinkage,” Construction and Building Materials, vol. 149, pp. 705–715, 2017. View at: Publisher Site  Google Scholar
 J. Q. Wang and Z. Fang, “Study on early crack in the web of longspan PC box girder bridge,” Journal of Hunan University Natural Sciences, vol. 39, no. 5, pp. 6–12, 2012. View at: Google Scholar
Copyright
Copyright © 2020 Shiwei Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.