Thermo-Hydro-Mechanical Combined Effect Analysis Model for Early-Age Concrete Bridges and Its Application

(e early cracking of concrete beam bridges remains a concern in civil engineering. An analytical model considering the combined effect of thermo-hydro-mechanical 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 thermo-hydro-mechanical 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 long-span 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. (e results show that the rate of moisture conduction is significantly lower than the rate of temperature conduction; even for thin-walled components, there exists a significant humidity gradient on the surface layer. (e humidity-induced shrinkage and restraint of the precast members are the main causes of web cracking.


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 thermo-hydro-mechanical processes [1 ,2]. is 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 thermohydro-mechanical processes accumulates to a certain extent, early-age cracks are induced in the concrete. e early-age 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. e 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. e cracking phenomenon remains to be quite frequent during construction [7]. e 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]. erefore, it is necessary to establish a suitable method to analyze the combined effect of the thermo-hydro-mechanical 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 thermo-chemomechanical 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 thermo-chemo-mechanical 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 chemo-mechanical model for early-age cracking of concrete in which the effect of relative humidity is also neglected. Gawin et al. [16,17] developed a hydro-thermalchemo-mechanical model of concrete at early ages which comprehensively accounts for the aging and hydrothermal coupling but lacks a comprehensive constitutive law for early-age concrete. Di Luzio and Cusatis [18] proposed the solidification-microprestress-microplane (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 thermo-hydromechanical 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.
e corresponding numerical analysis method was proposed and implemented into ABAQUS; the results were validated by conducting some typical tests. Finally, a long-span 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.

Coupling Analysis Theory of Temperature and Moisture Fields
2.1. eory 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.

Hydration
Heat. e active chemical components in the cement clinker react with the mixing water, releasing a lot of heat. e total heat release (H total ) considering the influence of clinker composition can be calculated as [20] H total � 500C 3 S + 260C 2 S + 866C 3 A + 420C 4 AF 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] H total � Q 0 (W + kF), (2) 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#. e hydration heat release (H(t)) at different times can be expressed as a function of the equivalent age (t e ), as shown in where α(t e ) is a function of the hydration degree and is related to t e . It can be expressed as where η 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 by where Blaine is the fineness coefficient of cement (m 2 ·kg −1 ); FA_CaO is the mass ratio of CaO to FA; w/c is the watercement ratio; HI is the hydraulic index and can be expressed as where 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 by where 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.

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] ρc . (11) 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 W w are the amounts of cement, aggregate, and water, respectively, by weight; c c , c a , and c w 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 . e 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 . e 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). us, the key performance parameters of the concrete temperature field evolve with the hydration process.

Initial and Boundary Conditions.
e 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]. e 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. e heat transfer between the successive segments is considered based on the fourth type of the boundary conditions.

2.2.
eory 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].

Humidity Reduction due to Cement Hydration.
Zhang et al. [23] proposed a two-stage calculation model for humidity reduction due to cement hydration based on a humidity conduction test. e model is expressed in where H su is the relative humidity (RH) considering self-desiccation 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.

Moisture Conduction.
According to the second Fick's law, the differential equation for concrete moisture conduction can be expressed as [23] zh zt where H s is the humidity reduction due to cement hydration and D(h) is the moisture transfer coefficient. 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); ; 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).

Advances in Civil Engineering
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.

Initial and Boundary Conditions.
e initial condition for moisture conduction in poured concrete is generally 100%. e 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.

Analysis of Stress Field
In addition to the corresponding self-stresses 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 thermo-hydro-mechanical processes, its mechanical properties evolve with time.

Mechanical Parameters of Early-Age
Concrete. e elastic modulus E, tensile strengths f t , and compressive strength f c of concrete varying with the equivalent age can be expressed as [26] E t e � E 28 · e 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.

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 in where α 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).

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. e 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]. e relationship between moisture-induced 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.

Tensile Creep of Early-Age
Concrete. e double power law (DPL) creep model proposed by Bazant and Panula [28] is selected because it has shown good accuracy in many early-age tensile creep tests [29,30]. e model is expressed in where 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.

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 in where < > are Macaulay brackets and < x > is taken as ((|x| + x)/2); I1 is the first stress invariant; J 2 represents the second deviatoric stress invariant; σ max 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; c is equal to 3 in the low confining pressure stress state.

Numerical Analysis Method of Thermo-Hydro-Mechanical Combined Effect
To realize the analysis of the thermo-hydro-mechanical 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. erefore, the whole iteration method is used.
is method does not significantly increase the calculation cost and shows sufficient calculation accuracy. e 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. e 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. e 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. e 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. e 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. e finite element software ABAQUS and its secondary development platform were selected for this complex process. e 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. e 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. e 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.

Temperature Field Simulation of a Prefabricated Box
Girder by Electric Heating. A prefabricated concrete box girder bridge of the Chongqing-Guizhou 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.
e 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 beam-end 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. e 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. e subroutine UMATHT was used for the temperature field analysis. e 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: e concrete molding temperature was 20°C. e 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). e 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). e equivalent surface conductance between concrete and air is χ s � 15 W/(m 2 ·°C). e heat conductivity k u is calculated using equation (9) and taken as 2.618 W/(m 2 ·°C). e specific heat capacity c is calculated using equation (10) and taken as 913.22 J/(kg·°C). e density of concrete is 2497 kg·m −3 . e 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. is demonstrates that the proposed method can effectively simulate the temperature field of concrete bridges.

Simulation of Moisture Conduction Process and Shrinkage
Deformation. Zhang et al. [33] conducted a moisture conduction test on early-age 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. e specimens were cuboid-shaped 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. e digital temperature and humidity sensors, which are located at the geometric center of the specimen, were used. e free deformation of the concrete was measured using LVDT.

Advances in Civil Engineering
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. erefore, the moisture field analysis method is consistent with the temperature field analysis method. e concrete was also modeled using DC3D8 with 2300 elements. e ambient temperature and humidity vary only slightly and are taken as 23°C and 16%, respectively. Table 3 lists the other relevant calculation parameters. e 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. e 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. e results show that the Heat radiation zone A A Figure 2: Simulation of temperature load [32]. humidity conduction rate is much lower than the temperature conduction rate; even for the thin-walled components, there remains a significant humidity gradient beneath the surface.

Simulation of Tensile Creep in Early-Age Concrete.
To understand the role of tensile creep in early-age concrete, Khan et al. [34] evaluated the creep deformation of concrete dog-bone 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. e loading ages for the dog-bone creep specimens were 1 and 2 d. e 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. e DPL creep

Early-Age Cracking Analysis of a PC Box
Girder Bridge with Cantilever Casting  [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. e 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.

Calculation Model.
e shrinkage deformation of 2# block is constrained by 1# block. erefore, 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. e 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. e foregoing analysis results show that a significant humidity gradient is formed beneath the surface. erefore, 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 thermohydro-mechanical analysis. ere are 193596 concrete elements and 19800 rebar elements, as shown in Figure 14.
e active and inactive element methods are used to simulate the construction process. e changes in the ambient temperature and humidity are simulated using the subroutine FILM. e hydration heat function is calculated using (27). e 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. e density of concrete is 2459.7 kg·m −3 .
e humidity reduction due to cement hydration reaction is calculated using equation (12). e humidity conductivity coefficient is calculated using the CEB-FIP (2010) formula considering temperature correction, as Simulation results of humidity-I Ambient humidity Simulation results of temperature-I Ambient temperature  Advances in Civil Engineering shown in equation (15). For the humidity boundary conditions, the main girder concrete should be kept wet during hardening. erefore, 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). e stress relaxation due to the creep is calculated using (24). q, d, and p are taken as 0.6, 0.2,    Figure 15 shows a comparison between the measured data and the simulation results at measuring points 1 and 3. e simulated and measured results are consistent, indicating that the temperature field of the bridge is reflected in the simulation. Figure 16 shows the humidity simulation results at points 1 and 3 under various calculation conditions. e 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. e lower the ambient humidity is, the higher the humidity gradient along the thickness of the web increases with time.   Figure 17 shows the web cracking simulation results under three calculation conditions; the web cracks are colored in red. e results show the following:

Analysis of Stress Field Calculation Results.
(1) e cracks first appear at the interface, mainly because the free deformation of 2# block is restrained by 1# block. (2) e 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.
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). e 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),    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. erefore, 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. e 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.

Conclusion
A new method is proposed to analyze the thermo-hydromechanical 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. e humidity conduction rate is much lower than humidity conduction rate; even for thin-walled 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 thin-walled components, the combined effect of humidity-induced 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
e data used to support the findings of this study are available from the corresponding author upon request. 14 Advances in Civil Engineering