1. Introduction
For mass concrete structures, thermal stresses that may induce cracking are a major concern, and temperature control measures must be taken during the construction period. For hydraulic mass concrete structures, cracks are mainly caused by temperature deformation shortly after construction [
1]. This early cracking behavior of concrete is essentially the result of the trade-off between strength and stress during the hardening process. The properties of concrete change with age, and the strength gradually increases, which causes the concrete to release a large amount of heat [
2]. Under the combined action of the heat of hydration and the cooling of the external environment, a temperature gradient is generated inside the concrete structure.
In the past, to prevent cracks induced by thermal stresses, mass concrete structures were usually poured in layers to artificially increase the construction interval. Moreover, other temperature measures, such as the use of ice to precool aggregates, embedded cooling water pipes, and surface heat preservation, have been used to prevent temperature-induced cracks [
3]. However, these measures are not economical, and some measures, such as pipe cooling, are very complex and may seriously affect the construction schedule.
In research and practice, during the dam concrete preparation process, an appropriate amount of magnesium oxide that has been burned under high temperatures can be added to make the concrete expand, which can compensate for shrinkage deformation [
4,
5] during a temperature drop. In combination with regular auxiliary measures (e.g., surface heat preservation and maintenance), MgO can be added as an expansion agent for mass concrete to control self-generated volume deformation and thereby reduce and eliminate dam concrete cracks. This method can prevent dam cracks and reduce the need for most traditional concrete temperature control measures (aggregate precooling, ice mixing, embedded cooling water pipe, etc.) [
6].
The addition of MgO in mass concrete is a superior technique that can reduce project construction and management costs by simplifying the temperature control measures. However, the expansion mechanism of MgO concrete is very complex [
7], not only can the hydration of the hydration of cement and the MgO itself have an impact on the autogenous deformation, but the fly ash in the cementing material can have an important influence on the MgO concrete, as well [
8]. Thus, further research is required. Usually, the autogenous volume deformation (expansion) and expansion stress during construction need to be simulated numerically, and appropriate models and reasonable parameter studies are crucial for such simulations. Therefore, a model that properly describes the expansion characteristics of concrete based on material properties is necessary.
Since the 1970s, MgO concrete has been used in a large number of hydraulic engineering projects in China, and a large amount of deformation data has been collected [
8]. B. Zhu (2002, 2003) [
9,
10] proposed a method for calculating the difference between the indoor test results and the measured values for the autogenous volumetric deformation of micro-expansion concrete and established a calculation model of the autogenous volume deformation by considering this difference. Based on engineering test data and relevant laboratory tests, G. Yang et al. (2004) [
11] suggested that the autogenous volumetric deformation characteristics of MgO concrete under constant temperature conform to a hyperbolic model. Using simulation techniques, G. Zhang et al. (2002, 2005) [
12,
13,
14] comprehensively analysed the influences of concrete elastic modulus growth, creep degree with age and temperature on the amount of magnesium oxide expansion and studied the influences of curing temperature and dam concrete temperature on the effective compensation. S. Liu et al. (2006) [
15] established a mathematical model of the relationships between MgO concrete autogenous volume deformation and MgO content, fly ash content and temperature. P. Xu et al. proposed the equivalent age method for calculating the spontaneous volume deformation of MgO concrete. (2008) [
16]. Other models, such as the arctangent curve model of C Chen et al. (2008) [
17], have also been proposed.
The establishment of an expansion and deformation model of magnesium oxide concrete is very important, and many studies that have been carried out can serve as references for practical use and research and are conducive to the establishment of a simulation model and the technical development and popularisation of MgO expansion concrete. However, a unified standard model that considers the chemical reactions of MgO concrete has not yet been developed. Although some models can consider the effect of temperature on the expansion characteristics of magnesium oxide, the calculation methods are controversial because they do not consider the nature of the chemical reaction (i.e., the material properties are not directly related to the chemical reaction degree). Therefore, additional research and relevant projects are needed.
In this paper, based on hydration kinetics, a coupled thermo-chemical-mechanical model is proposed for describing the hydration properties of magnesium oxide concrete. In this model, two variables are introduced as the intermediate variables of the chemical reaction system within concrete: the hydration degree of cementitious materials, which determines the hydration heat release, and the hydration degree of magnesium oxide, which determines the autogenous volumetric deformation. In addition, the model is numerically implemented with the finite element method, and the temperature, hydration and stress field within concrete can be simulated. Based on a parameter study, the model is applied in a case study that considers an arch dam to evaluate the coupled thermo-chemical-mechanical responses of mass concrete during construction.
4. Validation
4.1. Experimental Study of Magnesium Oxide Concrete
To scrutinize the deformation characteristics, which are closely related to the hydration properties, necessary experiments should be conducted to obtain the corresponding thermal and mechanical parameters. In this work, the experimental study comes from the engineering field (the Sanjiang Archdam Projiect in Guizhou, China). The usual test indexes of concrete include compressive strength, split tensile strength, bulk density, elastic modulus, ultimate tensile value, linear expansion coefficient, adiabatic temperature rise, thermal conductivity, specific heat and autogenous volume deformation. For MgO concrete, the low expansive characteristic caused by the hydration of magnesium oxide is significant and complicated.
The autogenous volume deformation can be monitored by using a digital electric bridge. A strainmeter is placed in the center of a concrete specimen that has been poured in a bottle. The height of the concrete cylinder is 10 cm and the radius is 60 cm. The resistance can be measured by the digital electric bridge, which is connected to the wire of the strainmeter, and then, the autogenous volume deformation can be calculated. Because the ambient temperature has a considerable impact on the test results, a thermostatic water tank is used to adjust the curing temperature of the specimens.
4.1.1. Materials and Methods
Special ordinary Portland cement is used for tests which is produced by Guizhou cement plant in China. The purity of magnesium oxide is greater than 90% and it is produced by Haicheng, Liaoning, China. The fly ash is secondary ash produced by Guizhou Thermal Power Plant.
Table 1 summarizes the mass fractions of the components of two types of MgO concrete with a gradation degree of three.
4.1.2. Test Results
The corresponding adiabatic temperature rise of the two types of MgO concrete at different ages are listed in
Table 2, and the initial temperature for the adiabatic test is the indoor temperature of 20 °C.
Table 3 gives the measured autogenous volume deformation results of the concrete specimen for curing temperatures of 20 °C and 30 °C.
4.2. Activation Energy of the Magnesium Oxide Hydration
While the activation energies of the cement hydration of concrete for calculations are plentiful in the literature, they are not easily found for magnesium oxide applications. Moreover, the activation energy of magnesium oxide hydration is closely related to the magnesium oxide production technique; thus, the activation energy values in different studies are quite different. Therefore, in this work, a method that can be used to obtain the activation energy of magnesium oxide hydration within concrete is proposed.
According to Equation (7), for an arbitrary degree of hydration of magnesium oxide (
), the expansion rate is related to only the temperature of concrete. Thus, for given curing temperatures, the ratio of the expansion rate is a function of temperature:
Therefore, to reach the same level of hydration, the required time is related to the curing temperature as well. Equation (7) can be transformed to
By integrating both sides of this equation,
Then, the time required to reach a certain degree of hydration under different temperatures can be described as:
Ultimately, the chemical reaction activation energy of MgO hydration
Em can be determined:
According to Equation (38) and combined with the experimental data of the autogenous volume deformation of oxygen magnesia concrete samples at different ages with varying curing temperatures, the activation energy can be obtained by choosing two different curing ages at which the hydration degrees are equivalent. In practice, the measurement of autogenous volume deformation cannot be conducted at every age of the magnesium oxide concrete; therefore, close values may be chosen, and the interpolation method should be used to obtain appropriate estimates of Em.
4.3. Numerical Verification
As discussed in
Section 2, the exothermic and autogenous volumetric 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. Based on the results of the adiabatic experiments and by adopting the strategy in reference [
21], the results of the parameter calibration of the thermo-chemical model that describes the properties of cement hydration are listed in
Table 4. Due to a lack of experimental data of the adiabatic temperature rise with different initial temperatures, the reaction activation energy can be set to 41570 J mol
−1 based on references [
22,
23], and
is equal to 5000 K.
The temperature evolution of the simulations and the adiabatic tests on the concrete samples with fly ash content 30% and 40% are plotted in
Figure 1. The results show that the predictions of the hydration model agree with the experimental values.
Based on the data in
Section 4.1 and by adopting the strategy from
Section 4.2, the activation energy of the magnesium oxide hydration in concrete can be obtained since the experimental data of autogenous volume deformation of oxygen magnesia concrete samples at different ages with curing temperatures are given.
Once the activation energy is obtained, the other parameters of the model can be calibrated on the basis of the experimental data of autogenous volume deformation. The results of the parameter calibration of the thermo-chemical model that describes the properties of MgO hydration are listed in
Table 5. As shown in the table, given the same MgO content and cement type, the reaction activation energy may vary.
The autogenous volumetric strain evolution of the simulation values and the experimental data for the concrete samples with fly ash content 30% and 40% are plotted in
Figure 2. Comparing the simulation values and experimental data, the hydration model proposed in this work shows a good accuracy under different curing temperatures.
4.4. Engineering Application
In
Section 4.3, the thermo-chemical model parameters are calibrated, and the simulation values show good agreement with the experimental data. In this section, an arch dam in Guizhou, China, was selected for thermo-chemical-mechanical analysis. Thus, the hydration properties and mechanical responses of magnesium oxide in mass concrete can be explored.
Mechanical tests on concrete are typically conducted under isothermal conditions; therefore, the hydration process is simulated under isothermal conditions in this research, and the mechanical characteristics can be linked to the hydration degree rather than to time. The evolution of the Young’s modulus for the isothermal test of the fly ash concrete specimens with fly ash content 30% and 40% are plotted in
Figure 3. For simplification, Poisson’s ratio is assumed to be constant (0.167).
The Sanjiang Arch Dam, which was constructed with magnesium oxide concrete during 2002 and 2004, was selected as a typical structure for thermo-chemical-mechanical analysis. 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 is poured. Thus, the incremental pouring process of concrete can be simulated in the FE model.
For mass concrete structures, thermal stresses are a major concern, and temperature control measures must be taken during construction. There are multiple material and structural factors that impact thermal stress, such as concrete hydration properties, thermal and mechanical boundary conditions, the construction schedule and temperature control measures; for simplification, these conditions are not described in detail. In subsequent studies, these details will be discussed, and the simulation results will be compared with the monitored values.
4.5. Results and Discussion
Based on the parameter studies and numerical validation in
Section 4.3 and
Section 4.4, temperature and stress fields can be simulated using the hydration model.
Figure 4 and
Figure 5 show the coupled thermo-chemical-mechanical field simulations on different dates. As shown in
Figure 4 and
Figure 5, the temperature is uneven within the mass concrete structure, and the hydration degree is also present in an inhomogeneous distribution. For the inner region with poor heat dissipation conditions, the temperature appears to be higher than that of exterior region on the same height. The hydration degree of the cement and that of MgO present the same pattern as the temperature distribution, as well.
The temperature and stress evolution processes are discussed for further comparative analysis. Thus, twelve typical points on the crown cantilever are selected, and points ①, ②, ③, ④, ⑤ and ⑥ are on the upstream profile (exterior areas of the dam), whereas points (1), (2), (3), (4), (5) and (6) are distributed over the midplane (interior areas of the dam). Accordingly, points ① and (1), points ② and (2), points ③ and (3), points ④ and (4), points ⑤ and (5), and points ⑥ and (6) are at six different elevations. The position distribution of the typical points is shown in
Figure 6. Thus, the pouring layers of the typical points are different, as is the age of the concrete.
Figure 7 shows the temperature history of each typical point, and
Figure 8 shows the maximum principal stress histories of the typical points. The temperature procedure differences give rise to different effective stress development histories. Under normal conditions, the tensile stress increases as the temperature decreases. However, due to the expansion characteristics of MgO concrete, the autogenous volumetric strain can compensate for the shrinkage strain caused by the decrease in temperature, and the tensile stress can be reduced effectively. As shown in the stress evolution figure (
Figure 8), for the typical points in both interior and exterior areas of the dam, the tensile stress does not show an upward trend for quite a few months, while the temperature decreases during that period, as shown in
Figure 7.
For the nodes on the downstream profile of the dam section, the temperature rises to a certain extent and then changes with air temperature, while the temperature of the nodes on the midplane rises higher and drops more slowly due to the poor heat dissipation conditions. Therefore, for points distributed in the interior areas of the dam, the shrinkage compensating effect is superior to that of the points distributed on the exterior areas of the dam. As shown in
Figure 8, for the exterior points, the period during which the tensile stress does not increase is obviously shorter than that in the interior areas of the dam.
Figure 9 and
Figure 10 show the cementitious material and magnesium oxide hydration degree histories of the selected typical points, respectively. As shown in the figures, the hydration degrees of both the cementitious material and the magnesium oxide are inhomogeneously distributed within the mass concrete structure and are affected by the distribution of temperature. The poor heat emission condition of the inner part of the dam accelerates the hydration process of the concrete in these regions. Thus, the increase in temperature (hydration of cementitious material) and expansion volumetric deformation (hydration of MgO) of the concrete on the upstream and downstream surfaces lags behind that of the inner regions.
Moreover,
Figure 9 and
Figure 10 show that the differences of cement hydration degree for the inner and exterior regions appear to be smaller than that of magnesium oxide hydration degree, especially at the late stage. This phenomenon may result from the activation energy of MgO hydration being larger than that of cement hydration. However, quantitative analysis for differences of internal and external expansion needs further experimental data as well as monitored data which would be scrutinized in future work.
5. Conclusions
In this work, a coupled thermo-chemical-mechanical model is proposed to study the expansion characteristics and mechanical responses of mass concrete structures. Focusing on the chemical-physical responses of mass concrete structures, parameter studies and numerical simulations were conducted. The major achievements are as follows:
The coupled thermo-chemical-mechanical hydration model separates the hydration degree of the cementitious material and the magnesium oxide. In addition, the heat release is determined by the hydration degree of the cementitious material, while the autogenous volumetric deformation is determined by that of magnesium oxide.
The thermo-chemical model is numerically implemented in the finite element method. In each load step of the finite element analysis, the hydration of both the cementitious material and the magnesium oxide can be updated synchronously.
The activation energy of magnesium oxide hydration, which is an important parameter in hydration kinetics, has rarely been reported in the literature. To overcome the problem, a mathematical derivation method for obtaining the activation energy is introduced.
Numerical simulations are conducted, and the results of the comparison of the experimental temperature and autogenous volumetric strain curves show that the hydration model in this work provides a relatively high accuracy.
The model is also applied in a case study on an arch dam, and the hydration degree of magnesium oxide within mass concrete, which can lead to strain compensation for shrinkage caused by temperature drop, is highlighted. The increase in temperature (hydration of cementitious material) and expansion volumetric deformation (hydration of MgO) of the concrete on the upstream and downstream surfaces lags behind that of the inner regions.
The differences of cement hydration degree for the inner and exterior regions and that of magnesium oxide hydration degree show not to be synchronized; this may be caused by the difference of cement hydration activation energy and MgO hydration activation energy. Meticulous numerical simulations and quantitatively contrastive analyses with current traditional model as well as monitoring data should be scrutinized in future work.