Contrastive Numerical Investigations on Thermo-Structural Behaviors in Mass Concrete with Various Cements

This work is a contrastive investigation of numerical simulations to improve the comprehension of thermo-structural coupled phenomena of mass concrete structures during construction. The finite element (FE) analysis of thermo-structural behaviors is used to investigate the applicability of supersulfated cement (SSC) in mass concrete structures. A multi-scale framework based on a homogenization scheme is adopted in the parameter studies to describe the nonlinear concrete behaviors. Based on the experimental data of hydration heat evolution rate and quantity of SSC and fly ash Portland cement, the hydration properties of various cements are studied. Simulations are run on a concrete dam section with a conventional method and a chemo-thermo-mechanical coupled method. The results show that SSC is more suitable for mass concrete structures from the standpoint of temperature control and crack prevention.


Introduction
Cracking is a major concern that should be considered for mass concrete structures from the standpoints of structures and materials. Thus, for mass concrete structures in which thermal-stress-induced cracks are prone to occur, considerations should include not only rational structure types and reasonable construction sequences but also suitable cement materials. Supersulfated cement (SSC) is particularly suitable for use in the internal part of hydraulic structure mass concrete due to its slow hydration rate and very low hydration heat.
For the cement industry, re-use binders are used for waste materials from other industries; these binders are effective for energy-saving emissions reduction to reduce energy consumption and CO 2 emissions [1,2]. SSC is a type of cement with little or no clinker, which contains nearly 80% ground granulated blast furnace slag, 15% gypsum as the sulfate activator and less than 5% clinker or lime as the alkali activator [3]. The hydration heat of SSC is much lower than that of ordinary Portland cement; therefore, SSC is thought to have promising potential for industrial applications [4].
For mass concrete structures, reasonable joint arrangement and stratification are significant. Concrete is poured in layers, which offers convenience in construction and improves the heat dissipation conditions. Moreover, pre-cooling techniques, such as cooling aggregates and reducing cement content, and the post-cooling of embedded cooling pipes are widely used. The first successful post-cooling technique was employed during the construction of the Hoover Dam [5] and is still used today. Considering the advantage of low hydration heat of SSC concrete, water-cooling is likely to be cancelled, which can bring about construction cost-savings for mass concrete structures.
In most cases, proper materials may be selected to satisfy the strength requirements; reasonable construction and cooling schedules may also be established according to regular calculations. However, because regular calculations typically simplify the hydration process and difficulties arise in predicting concrete properties such as the modulus of elasticity [6], cracks may still occur unpredictably in mass concrete structures and affect their durability and service performance. Therefore, it is necessary to simulate the thermal field in a proper manner that considers the inner temperature, hydration heat generation rate and strength development rate in mass concrete so the effective temperature control measures can be adopted to prevent cracking.
Since the effect of the temperature and the cement hydration can promote or hinder in the concrete formation of micro-cracks which in turn influence the behavior of the material itself, the understanding of this phenomenon is very relevant. Indeed, recent works [7,8] have shown how the presence of micro-cracks can result, sometimes, in the improvement of the dissipative capacity of mass concrete structures. Thus, the accurate prediction of thermal-mechanical behaviors is of great significance for concrete structures.
In general, computer simulations (typically based on the finite element method) allow us to precisely evaluate the thermal evolution of early-age concrete by considering the hydration degree as an internal variable and other related parameters in describing the thermo-chemical-mechanical coupling process [9][10][11][12][13][14][15]. Cervera et al. proposed a coupled thermo-chemo-mechanical model for the behavior of early-age concrete that can be implemented into the finite element method and predict the evolution of the hydration degree and heat over time [16][17][18]. Lackner et al. found a significant linear correlation between hydration degree and the intrinsic material functions for the fracture energy [19]. Gawin et al. proposed a model that not only took the thermal and chemical fields but also the water diffusion process coupled effects into consideration [20]. Taking various chemical reactions explicitly into account, Di Luzio and Cusatis introduced a new hygro-thermo-chemical model for high-performance concrete suitable for the analysis of moisture transport and heat transfer [21,22].
In this paper, one type of environment-friendly cement, SSC, is introduced, and its hydration properties are studied based on hydration heat liberation rate data. A thermo-chemical-mechanical coupled hydration model was adopted to simulate the temperature and stress field in a relatively precise way for mass concrete structures. A multi-scale framework is employed to obtain the nonlinear thermal and mechanical characteristics. The simulation model is applied in a case study on a monolith of an arch dam.

Conventional Method for Concrete Hydration Heat and Elastic Modulus Estimation
For mass concrete structures, the hydration heat is the main heat source, which is significant in simulating the temperature and stress field. The conventional method to calculate the hydration heat of concrete is to employ the adiabatic temperature rise model, which considers the hydration heat as a function (usually composite exponential or hyperbolic type) of hydration time. For mass concrete structures, one of the most widely used models is given by [23]: where Qptq is the hydration heat at time t, and Q 8 is the ultimate hydration heat; a and b are material parameters that can be calibrated with adiabatic test data. Adiabatic tests are conducted with an initial temperature of 20-25˝C. From the viewpoint of cement hydration kinetics, concrete hydration occurs quickly under adiabatic conditions. Thus, variations in the adiabatic temperature rise are generally too low to measure after 28 days [24]. Therefore, using the above model, hydration heat evolution presents a fat trend after of 28 days, which shows a different method than standard engineering practices. Owing to lower placing temperatures (approximately 10˝C or lower) in mass concrete structures, hydration processes are retarded, then hydration heat is still generated at a considerably slower rate after 28 days.
Similar to the hydration heat calculation, the elastic modulus of mass concrete is described as follows: Eptq where Eptq is the elastic modulus at time t, and E 8 is the ultimate elastic modulus; α and β are material parameters.

Long-Term Hydration Model
Hydration models based on the Arrhenius concept are widely accepted by researchers because they use rigorous thermodynamic theory and are consistent with experimental data [25]. According to Arrhenius's law, the time-dependent hydration progress could be described as in Equation (1). In the modeling, the hydration rate is related to the temperature and the chemical affinity A ξ pξq: where E a is the activation energy of the reaction, and R is a constant for ideal gases. Cervera et al.
developed an analytical form of the normalized affinity based on thermodynamics [17]. The free energy for the thermo-chemical system can be divided into three parts: thermal contribution, thermo-chemical coupling, and chemical contribution. The chemical contribution in this work will be considered as a quartic function instead of a cubic function. Similar to the strategy presented in [17], a fixed form of the chemical affinity can be derived as follows [26]: where β 1 , β 2 and β 3 are material coefficients, ξ 8 is the ultimate hydration degree, and η represents the viscosity due to micro-diffusion of the free water through the already-formed hydrates. The advantage of the fixed hydration model is discussed in [26]. Notably, ξ 8 is set as one in this paper because the construction periods of mass concrete structures are quite long (typically several months or even years). In this way, the number of arguments the function takes can be reduced. Parameters β 1 , β 2 , β 3 and η can be calibrated using the experimental results.

Transient Heat Transfer Process
The transient heat transfer process of concrete can be described as follows: where ρ, C, λ T , ∆ and Q 8 are the density, volumetric heat capacity, thermal conductivity coefficient, Laplacian operator and the final volumetric heat of hydration, respectively. The latent heat release due to concrete hydration is actually a nonlinear and thermally dependent process.

Thermal-Creep Stress Calculation
According to the studies of Di Luzio and Cusatis [27], based on the assumption that the strain additivity holds, the composition of concrete strain mainly includes elastic strain, cracking strain, creep (or viscoelastic) strain, thermal strain, and shrinkage strain [28]: where . ε e , . ε c , . ε cr , . ε th and . ε sh are the elastic strain rate, cracking strain rate, creep strain rate, thermal strain rate and shrinkage strain rate, respectively.
Based on engineering facts, the thermal stress of mass concrete structures is usually calculated by a viscoelastic model. This paper aims to conduct contrastive analysis of conventional computations with the coupled hydration model; thus, thermal-chemical coupled effects which used to be ignored in engineering practice are our major considerations. The constitutive law in this paper for the mechanical behavior is the conventional one for engineering calculations, which is simplified.
Accordingly, the stress evolution can be described as follows:

Shrinkage
Total shrinkage of concrete may be composed of three parts: the autogeneous shrinkage, carbonation shrinkage and drying shrinkage. In good-quality concrete, carbon dioxide penetrates only a very thin surface layer, and the carbonation shrinkage is usually small and can be neglected [29].
As concrete hardens, the pore structure would change, which leads to the development of autogenous shrinkage. As a consequence, autogenous shrinkage is closely related to the evolution of hydration in the material. Experimental results indicate that autogenous shrinkage evolution is in a linear relationship with the hydration degree [10], so autogenous shrinkage ε as can be modeled by linear function or piecewise linear function [30]: where k a and k b are material parameters, 1 is the unit tensor and ξ 0 is a threshold. Considering that the drying shrinkage behaviors usually take place on the surface of structures, and humidity diffusion coefficient is quite small (under the environment with relative humidity of 50%, it may take a month to reach a drying depth of 7 cm, which is far smaller than the size of mass concrete structures); as a consequence, for the inner part of mass concrete structures, the humidity can be treated as constant [23]; thus, the drying creep effect is usually neglected for mass concrete structures.

Creep
The mass concrete creep can be modeled by means of the solidification theory. Bofang Zhu compared the generalized Kelvin rheological model with the Dirichlet series of creep compliance and summarized a creep model with eight parameters, which has been widely adopted in China [23]. The formulations can be deduced from Kelvin chain model (with two units in series): where t is the time, τ is the age of concrete, Epτq is the elastic modulus of the concrete, Jpt, τq is the creep compliance function, Cpt, τq is the creep function, and A i -D i are material parameters.
Since the creep experiments usually show the total creep strain in practice, and as has been discussed in Section 2.4.1, the humidity is usually treated as constant within mass concrete structures. In this work, additional creep due to drying (drying creep), also called the stress-induced shrinkage (or Pickett effect) [31], would be simplified, as well.

Thermal Strain
The thermal strain ε th is related to the temperature variation: where α t is the thermal dilatation coefficient (kept constant), and 1 is the unit tensor.

Thermal and Mechanical Characteristics
In a fully coupled thermo-chemical-mechanical model, the hydration degree can be used as an intermediate variable to determine the thermal or mechanical characteristics. In this work, the heat capacity and thermal conductivity of concrete are estimated using a homogenization scheme that is well-suited for heterogeneous material.

Heat Capacity
Given the heat capacity and volume fraction values of cement, water and aggregates, the heat capacity of the matrix can be estimated using the law of mixtures [32]: where c f reshpaste , c paste pξq, c water , c cement , c f iller , and c agg are the heat capacity (per unit volume) of fresh paste, paste in the hydration process, water, cement, the filler, and the aggregates, respectively. f water , f cement , f paste , f f iller , and f agg are the volume fraction of water, cement, the paste, the filler and the aggregates, respectively. A and B are material constants; in Bentz's study [32], these values were measured with w/c varying between 0.3 and 0.5.

Thermal Conductivity
The paste matrix is made of water and cementitious powder, so a reasonable approximation of thermal conductivity can be calculated from the well-known Hashin-Shtrikman (H-S) bounds. An estimate can be determined from the measured thermal conductivities (λ 1 ,λ 2 ) and phase volume fractions ( f 1 , f 2 ). For λ 2 ě λ 1 , the Hashin-Shtrikman lower (λ l ) and upper (λ u ) bounds for the thermal conductivity of a two-phase composite are given by [32]: An estimate can then be determined by the average of the H-S bounds: where λ hom is the effective thermal conductivity of the composite. According to [32,33], the thermal conductivity of cement paste is a function of volume fraction rather than hydration degree. Thus, this work ignores the nonlinear behaviors of thermal conductivity.

Aging Degree
Experimental evidences show that the evolution of the concrete strength depends not only on the degree of hydration, but also on the kinetics of the hydration reaction. On the basis of this evidence, Cervera et al. proposed an aging degree model to link strength evolution with hydration degree and current temperature, and the aging degreeκ, can be defined as [17,34]: where T max represents the maximum temperature during hydration of concrete which is possibly under standard conditions (100˝C); T ref is the reference temperature for the experimental calibration of the aging model; ξ set defines the value of the hydration degree at the end of the setting phase; and n κ , B κ and A κ are material parameters.

Compressive and Tensile Strength
The Young's modulus values vary due to the aging degree's evolution as follows: where E c pξq and E t pξq are compressive and tensile strength at the degrees of hydration ξ and ξ 8 , respectively, E c 8 and E t 8 are compressive and tensile strength, respectively [17,35] .

Young's Modulus
The Young's modulus varies due to the aging degree's evolution as follows: where E pξq and E 8 are Young's modulus at degree of hydration ξ and ξ 8 , respectively.

Poisson's Ratio
Based on the work of De Schutter et al. [36,37], the Poisson's ratio of concrete can be estimated as follows: where µ pξq and µ 8 are Poisson's ratio at the degrees of hydration ξ and ξ 8 , respectively.

Calibration Method of the Hydration Model
As discussed in Section 2, the value of thermo-chemical-mechanical characteristics are associated with the temperature or hydration degree values. In this section, based on the experimental values, the hydration model in this paper is numerically verified, and the parameters are calibrated.
At present, the kinetic parameters of the hydration reaction of concrete were achieved by testing temperature or the rate of hydration heat evolution with the adiabatic tests with the following equations: ξptq " where Q is the hydration heat evolution quantity; Q 8 is the ultimate hydration heat evolution quantity; ξptq is the hydration degree; t 0 is the end time of the induction period, t 0.5 is the time moment when the hydration degree is 0.5. Once the hydration heat evolution is obtained, Equations (18) and (19) are combined, and an inverse analysis is adopted to determine the parameters of the hydration model.

Fly Ash Concrete
Using the approach in Section 3.1, the temperature evolution of the simulations and adiabatic tests on concrete samples are plotted in Figures 1 and 2. The results show that the hydration model simulation conforms to the experimental values (the error is no more than 2%). Tables 1 and 2 summarize the mass fractions of the components and the thermal parameters of two types of fly ash concrete with degree four of gradation. Table 3 lists the results of the parameter calibration for the hydro-thermal model.  Table 2. Thermal parameters of the components.
Parameters     Mechanical tests on concrete are typically conducted under isothermal conditions; therefore, the hydration process is simulated under isothermal condition in this research, so the mechanical characteristics can be linked to the hydration degree rather than to time. Young's modulus and the Poisson's ratio evolution for isothermal test of fly ash concrete specimen are plotted in Figures 3 and 4. Mechanical tests on concrete are typically conducted under isothermal conditions; therefore, the hydration process is simulated under isothermal condition in this research, so the mechanical characteristics can be linked to the hydration degree rather than to time. Young's modulus and the Poisson's ratio evolution for isothermal test of fly ash concrete specimen are plotted in Figures 3 and 4.
The autologous volumetric deformation of fly ash concrete specimen are plotted in Figures 5 and 6 as has discussed in Section 2.4.1; piecewise linear functions are adopted to simulate the autologous volumetric deformation.

SSC Concrete
The hydration properties of SSC are different from slag Portland cement and ordinary Portland cement [38]. SSC achieves a lower heat release compared to slag Portland cement and ordinary Portland cement. Conversely, slag cannot be simulated without a suitable alkaline solution, so the second exothermal peaks of the SSC present later than that of slag Portland cement and ordinary Portland cement. On the basis of the experimental studies in [3,38], the hydration properties are simulated in a numerical method.

SSC Concrete
The hydration properties of SSC are different from slag Portland cement and ordinary Portland cement [38]. SSC achieves a lower heat release compared to slag Portland cement and ordinary Portland cement. Conversely, slag cannot be simulated without a suitable alkaline solution, so the second exothermal peaks of the SSC present later than that of slag Portland cement and ordinary Portland cement. On the basis of the experimental studies in [3,38], the hydration properties are

SSC Concrete
The hydration properties of SSC are different from slag Portland cement and ordinary Portland cement [38]. SSC achieves a lower heat release compared to slag Portland cement and ordinary Portland cement. Conversely, slag cannot be simulated without a suitable alkaline solution, so the second exothermal peaks of the SSC present later than that of slag Portland cement and ordinary Portland cement. On the basis of the experimental studies in [3,38], the hydration properties are simulated in a numerical method. Figure 7 shows the cumulative heat release quantity of the model simulations and experimental values for SSC; the simulations are in accordance with the experimental data (the error is within 4%). The parameters of the hydrothermal model of SSC are listed in Table 4.   The hydration properties can be analyzed, given the kinetic parameters of the hydration reaction of the cement paste and the composition of concrete, which is listed in Table 5. Assuming the aggregates have no impact on the hydration process, the ultimate hydration heat evolution quantity can be obtained as below: where e Q  is the ultimate hydration heat evolution quantity of the cement paste, and cement f is the fraction of cement paste in the concrete. The compressive strength evolution for isothermal tests of the SSC specimen is plotted in Figure  8. The hydration properties can be analyzed, given the kinetic parameters of the hydration reaction of the cement paste and the composition of concrete, which is listed in Table 5. Assuming the aggregates have no impact on the hydration process, the ultimate hydration heat evolution quantity can be obtained as below: where Q e 8 is the ultimate hydration heat evolution quantity of the cement paste, and f cement is the fraction of cement paste in the concrete. The compressive strength evolution for isothermal tests of the SSC specimen is plotted in Figure 8.

Project Background and Assumptions
The Dagangshan Hydropower Station is one of the largest hydropower projects built on the main stream of Dadu River in Southwest China. The water retaining structure of the Dagangshan water complex is a double-curvature arch dam with a maximum height of 210 m, 29 dam sections, and total concrete dosages of approximate 3,220,000 m 3 [39]. In this study, one of the dam sections, #18, was selected for the thermo-chemical-mechanical analysis.
For mass concrete structures, thermal stresses are a major concern, and temperature control measures are Gordian techniques in the construction period. There are multiple factors that impact thermal stress, such as concrete material properties, boundary conditions, construction schedule and temperature control measures.
The cooling effects of the cooling pipes can be estimated according to [24,40,41]; however, in this work for simplification, the cooling effects are not included in the calculations.

FE Model Description and Initial Conditions
The #18 dam monolith located on the riverbed was selected as a typical structure to conduct the thermo-chemical-mechanical analysis. The 3D FE mesh model for the simulation is presented in Figure 9. A birth and death technique is adopted to address the existence and nonexistence of the concrete. Corresponding elements in the finite element model are made active when a fresh concrete slab was poured. Thus, the incremental pouring process of concrete can be simulated in the FE model.

Project Background and Assumptions
The Dagangshan Hydropower Station is one of the largest hydropower projects built on the main stream of Dadu River in Southwest China. The water retaining structure of the Dagangshan water complex is a double-curvature arch dam with a maximum height of 210 m, 29 dam sections, and total concrete dosages of approximate 3,220,000 m 3 [39]. In this study, one of the dam sections, #18, was selected for the thermo-chemical-mechanical analysis.
For mass concrete structures, thermal stresses are a major concern, and temperature control measures are Gordian techniques in the construction period. There are multiple factors that impact thermal stress, such as concrete material properties, boundary conditions, construction schedule and temperature control measures.
The cooling effects of the cooling pipes can be estimated according to [24,40,41]; however, in this work for simplification, the cooling effects are not included in the calculations.

FE Model Description and Initial Conditions
The #18 dam monolith located on the riverbed was selected as a typical structure to conduct the thermo-chemical-mechanical analysis. The 3D FE mesh model for the simulation is presented in Figure 9. A birth and death technique is adopted to address the existence and nonexistence of the concrete. Corresponding elements in the finite element model are made active when a fresh concrete slab was poured. Thus, the incremental pouring process of concrete can be simulated in the FE model.
The construction monolith #18 of the Dagangshan Dam began on 11 March 2012, and it was completed on 25 September 2014. Next, 72 layers of concrete gradually rose with the progress of construction. The initial temperature of the concrete was assumed the placing temperature [23], according to the monitoring records. The construction monolith #18 of the Dagangshan Dam began on 11 March 2012, and it was completed on 25 September 2014. Next, 72 layers of concrete gradually rose with the progress of construction. The initial temperature of the concrete was assumed the placing temperature [23], according to the monitoring records.

Boundary Conditions
In normal conditions, the construction of the selected monolith was always behind the adjacent monoliths; the boundary condition for both sides was assumed adiabatic. The upstream and downstream surfaces of the placed concrete as well as the top surface of the concrete block and foundation were assumed directly exposed to the outside environment, where heat convection and radiation occurred. Now that the project has been completed, the distribution characters of the air temperature can be simulated based on the measured data.
According to the monitoring temperature of the meteorological station at the Dagangshan Dam site, the outside air temperature (Ta) of the measured temperature values (Figure 10

Boundary Conditions
In normal conditions, the construction of the selected monolith was always behind the adjacent monoliths; the boundary condition for both sides was assumed adiabatic. The upstream and downstream surfaces of the placed concrete as well as the top surface of the concrete block and foundation were assumed directly exposed to the outside environment, where heat convection and radiation occurred. Now that the project has been completed, the distribution characters of the air temperature can be simulated based on the measured data.
According to the monitoring temperature of the meteorological station at the Dagangshan Dam site, the outside air temperature (T a ) of the measured temperature values (Figure 10 The construction monolith #18 of the Dagangshan Dam began on 11 March 2012, and it was completed on 25 September 2014. Next, 72 layers of concrete gradually rose with the progress of construction. The initial temperature of the concrete was assumed the placing temperature [23], according to the monitoring records.

Boundary Conditions
In normal conditions, the construction of the selected monolith was always behind the adjacent monoliths; the boundary condition for both sides was assumed adiabatic. The upstream and downstream surfaces of the placed concrete as well as the top surface of the concrete block and foundation were assumed directly exposed to the outside environment, where heat convection and radiation occurred. Now that the project has been completed, the distribution characters of the air temperature can be simulated based on the measured data.
According to the monitoring temperature of the meteorological station at the Dagangshan Dam site, the outside air temperature (Ta) of the measured temperature values (Figure 10

Results and Discussion
Based on the parameter studies in Section 2 and numerical validation in Section 3, temperature and stress fields are simulated using both the conventional method and the hydration model. Figures 11 and 12 show the coupled thermo-chemical-mechanical fields simulations on different dates, accordingly; Figures 13 and 14 show the solutions of temperature and stress field using the conventional model. The results indicate that temperature distribution only differed within the relatively fresh concrete layers; temperature computations of the coupled model are smaller than that of the conventional model. However, different distribution patterns are present in the stress fields. Moreover, the maximum temperature resulting from the coupled model is lower than that of the conventional model, mainly because the hydration heat release process in the coupled model takes a long period, which corresponds to the engineering facts.

Results and Discussion
Based on the parameter studies in Section 2 and numerical validation in Section 3, temperature and stress fields are simulated using both the conventional method and the hydration model. Figures 11 and 12 show the coupled thermo-chemical-mechanical fields simulations on different dates, accordingly; Figures 13 and 14 show the solutions of temperature and stress field using the conventional model. The results indicate that temperature distribution only differed within the relatively fresh concrete layers; temperature computations of the coupled model are smaller than that of the conventional model. However, different distribution patterns are present in the stress fields. Moreover, the maximum temperature resulting from the coupled model is lower than that of the conventional model, mainly because the hydration heat release process in the coupled model takes a long period, which corresponds to the engineering facts.

Results and Discussion
Based on the parameter studies in Section 2 and numerical validation in Section 3, temperature and stress fields are simulated using both the conventional method and the hydration model. Figures 11 and 12 show the coupled thermo-chemical-mechanical fields simulations on different dates, accordingly; Figures 13 and 14 show the solutions of temperature and stress field using the conventional model. The results indicate that temperature distribution only differed within the relatively fresh concrete layers; temperature computations of the coupled model are smaller than that of the conventional model. However, different distribution patterns are present in the stress fields. Moreover, the maximum temperature resulting from the coupled model is lower than that of the conventional model, mainly because the hydration heat release process in the coupled model takes a long period, which corresponds to the engineering facts.     The temperature and stress evolution processes are discussed for further comparative analysis. Six typical point are selected, and points ①, ③ and ⑤ are on the downstream profile of the dam section, whereas, points ②, ④ and ⑥ are distributed over the midplane. Points ① and②, points ③ and ④, and points ⑤ and ⑥ are on the three different elevations. The position distribution of the typical points is shown in Figure 15. Thus, the pouring layers of the typical points are different, as is the age of the concrete. Figure 16 shows the temperature history of each typical point, and Figure 17 shows the maximum principle stress histories of the typical points.
For the nodes on the downstream profile of the dam section, the temperature rises to a certain extent and then changes with air temperature. Compared with the conventional model, the temperature prediction values of the coupled model present a slower growth trend inside the mass concrete. The temperature procedure differences give rise to differentiate the effective stress development history. For points on the surface, the calculated stress values of the coupled model are relatively larger than that of conventional model. However, the opposite occurs for the inner part of dam section.  The temperature and stress evolution processes are discussed for further comparative analysis. Six typical point are selected, and points ①, ③ and ⑤ are on the downstream profile of the dam section, whereas, points ②, ④ and ⑥ are distributed over the midplane. Points ① and②, points ③ and ④, and points ⑤ and ⑥ are on the three different elevations. The position distribution of the typical points is shown in Figure 15. Thus, the pouring layers of the typical points are different, as is the age of the concrete. Figure 16 shows the temperature history of each typical point, and Figure 17 shows the maximum principle stress histories of the typical points.
For the nodes on the downstream profile of the dam section, the temperature rises to a certain extent and then changes with air temperature. Compared with the conventional model, the temperature prediction values of the coupled model present a slower growth trend inside the mass concrete. The temperature procedure differences give rise to differentiate the effective stress development history. For points on the surface, the calculated stress values of the coupled model are relatively larger than that of conventional model. However, the opposite occurs for the inner part of dam section. The temperature and stress evolution processes are discussed for further comparative analysis. Six typical point are selected, and points 1 , 3 and 5 are on the downstream profile of the dam section, whereas, points 2 , 4 and 6 are distributed over the midplane. Points 1 and 2 , points 3 and 4 , and points 5 and 6 are on the three different elevations. The position distribution of the typical points is shown in Figure 15. Thus, the pouring layers of the typical points are different, as is the age of the concrete. Figure 16 shows the temperature history of each typical point, and Figure 17 shows the maximum principle stress histories of the typical points.
For the nodes on the downstream profile of the dam section, the temperature rises to a certain extent and then changes with air temperature. Compared with the conventional model, the temperature prediction values of the coupled model present a slower growth trend inside the mass concrete. The temperature procedure differences give rise to differentiate the effective stress development history. For points on the surface, the calculated stress values of the coupled model are relatively larger than that of conventional model. However, the opposite occurs for the inner part of dam section.           Figure 18 shows the hydration degree histories of the selected typical points. Owing to the existence of temperature heterogeneity, the hydration degree is present in an inhomogeneous distribution within the mass concrete structure. The poor heat emission condition of the inner part of dam accelerates the hydration process of the concrete in these regions. Thus, the growth of cement strength of concrete on the upstream and downstream surfaces lags behind that of the inner regions. Figure 18 shows the hydration degree histories of the selected typical points. Owing to the existence of temperature heterogeneity, the hydration degree is present in an inhomogeneous distribution within the mass concrete structure. The poor heat emission condition of the inner part of dam accelerates the hydration process of the concrete in these regions. Thus, the growth of cement strength of concrete on the upstream and downstream surfaces lags behind that of the inner regions.  Figure 19 shows the temperature evolution of the SSC concrete under adiabatic conditions. The results show that the adiabatic temperature of the SSC concrete increase less than 15 °C. For mass concrete structures such as the Dagangshan Dam, though pre-cooling techniques are adopted, water cooling techniques are still required to control the maximum temperature. For example, for Dagangshan Dam, the design maximum temperatures are 27 °C (for the foundation constrained region) and 30 °C (for the non-constraint region). Therefore, when using SSC concrete, the maximum temperature of mass concrete can satisfy the design requirements without the need of water cooling.   Figure 19 shows the temperature evolution of the SSC concrete under adiabatic conditions. The results show that the adiabatic temperature of the SSC concrete increase less than 15˝C. For mass concrete structures such as the Dagangshan Dam, though pre-cooling techniques are adopted, water cooling techniques are still required to control the maximum temperature. For example, for Dagangshan Dam, the design maximum temperatures are 27˝C (for the foundation constrained region) and 30˝C (for the non-constraint region). Therefore, when using SSC concrete, the maximum temperature of mass concrete can satisfy the design requirements without the need of water cooling. Figure 18 shows the hydration degree histories of the selected typical points. Owing to the existence of temperature heterogeneity, the hydration degree is present in an inhomogeneous distribution within the mass concrete structure. The poor heat emission condition of the inner part of dam accelerates the hydration process of the concrete in these regions. Thus, the growth of cement strength of concrete on the upstream and downstream surfaces lags behind that of the inner regions.  Figure 19 shows the temperature evolution of the SSC concrete under adiabatic conditions. The results show that the adiabatic temperature of the SSC concrete increase less than 15 °C. For mass concrete structures such as the Dagangshan Dam, though pre-cooling techniques are adopted, water cooling techniques are still required to control the maximum temperature. For example, for Dagangshan Dam, the design maximum temperatures are 27 °C (for the foundation constrained region) and 30 °C (for the non-constraint region). Therefore, when using SSC concrete, the maximum temperature of mass concrete can satisfy the design requirements without the need of water cooling.

Conclusions
In view of the disadvantages of the conventional temperature prediction model, this work introduces a fixed hydration model to study the non-linear thermo-chemical-mechanical responses in mass concrete structures. Focusing on the chemical-physical responses of mass concrete structures, the authors conducted numerical simulations and parameter studies. The major achievements are as follows:

‚
The conventional heat release model of mass concrete only relates the heat liberation rate to time, which ignores the effect of real-time temperature on the hydration process. To overcome the problem, a chemical-physical coupled model is introduced.
‚ Based on the chemical-physical coupled model, the hydration properties of fly ash concrete and SSC concrete are studied.
‚ A multi-scale framework is employed to obtain the nonlinear thermal characteristics. Using a homogenization scheme that is well suited for heterogeneous materials, the hydration process of mass concrete structures can be estimated.

‚
Contrastive analysis of the coupled model and the conventional model are conducted, based on an engineering project.

‚
Temperature prediction values of the conventional model and the coupled hydration model show a difference, especially at an early age when cracking is prone to occur within mass concrete structures.

‚
The cooling effect of the cooling pipes used inside the mass concrete was neglected. Thus, the chemical-physical responses under actual working conditions should be scrutinized in future research.
‚ When using proper materials such as super sulfated slag concrete, water cooling is likely to be abrogated, which can induce cost-savings.