Next Article in Journal
A Facile and Efficient Bromination of Multi-Walled Carbon Nanotubes
Next Article in Special Issue
Influence of the Precursor, Molarity and Temperature on the Rheology and Structural Buildup of Alkali-Activated Materials
Previous Article in Journal
The Influence of the Environment for Glass-Reinforced Plastic Composite Material Used for Ground Water Transport Pipes
Previous Article in Special Issue
Predictive Hydration Model of Portland Cement and Its Main Minerals Based on Dissolution Theory and Water Diffusion Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Mathematical Model of the Expansion Evolution of Magnesium Oxide in Mass Concrete Based on Hydration Characteristics

1
Guizhou Engineering Technology Research Center for Exploitation and Utilization of Water Resources in Karst Region, Guiyang 550002, China
2
State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, China
3
Physics of AmoRphous and Inorganic Solids Laboratory (PARISlab), Department of Civil and Environmental Engineering, University of California, Los Angeles, CA 90095, USA
*
Author to whom correspondence should be addressed.
Materials 2021, 14(12), 3162; https://doi.org/10.3390/ma14123162
Submission received: 18 April 2021 / Revised: 28 May 2021 / Accepted: 4 June 2021 / Published: 8 June 2021

Abstract

:
The low swelling property of magnesium oxide concrete is a significant feature that can be used to control the cracking of mass concrete. Based on the characteristics of the chemical reaction, this work proposes a coupled hydro-thermo-mechanical model that can be implemented with the finite element method for predicting the autogenous volumetric deformation of magnesium concrete. By introducing the degree of the hydration reaction of magnesia and the degree of the hydration reaction of cementitious materials as intermediate variables of the chemical reaction system, a prediction model of the concrete temperature and chemical fields is established, and using this model, the effect of the temperature on the reaction rate can be considered in real time. In addition, by combining the relationship between the degree of the hydration reaction of magnesium oxide and the comprehensive expansion of concrete, a mathematical model for calculating the expansion stress of magnesia concrete was established. The algorithms were derived by mathematical equations, and the simulation results were compared to the experimental temperature and autogenous volumetric strain curves, which showed that the hydration model provides a relatively high accuracy. The model was also applied to an arch dam, and the coupled thermo-chemical-mechanical responses of mass concrete during construction were investigated. Simulation results show that 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 obviously behind that of the inner regions. Quantitative analysis for differences of internal and external expansion is worthy of further attention and study on a basis of further experimental data as well as monitored data.

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.

2. Coupled Thermo-Chemo-Mechanical Model

2.1. Transient Heat Transfer

The transient heat transfer process of concrete can be described as follows:
ρ C T ( x , t ) t = λ T Δ T ( x , t ) + Q t
Q t = Q ξ t
where ρ , C , λ T , Δ and Q are the density, volumetric heat capacity, thermal conductivity coefficient, Laplacian operator and final volumetric heat of hydration, respectively. The latent heat release due to concrete hydration is a nonlinear and thermally dependent process.

2.2. A Chemical Model for the Hydration Process

According to the Arrhenius law, the time-dependent hydration progress can be described as in Equation (2). In the model, the hydration rate is related to the temperature and the chemical affinity A ξ ( ξ ) :
ξ t = A ξ ( ξ ) exp ( E a R T )
where E a is the activation energy of the reaction and R is the ideal gas constant. Cervera et al. developed an analytical form of the normalized affinity based on thermodynamics [18]. The free energy of the thermo-chemical system can be divided into three parts: the thermal contribution, the thermo-chemical contribution and the 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 [19], a modified form of the chemical affinity can be derived [18,20]:
A ξ ( ξ ) = β 1 ( β 2 + ξ ) ( ξ - ξ ) exp ( η ¯ ξ ξ )
where β 1 , β 2 and β 3 are material coefficients, ξ is the ultimate hydration degree and η ¯ represents the viscosity due to the microdiffusion of free water through hydrates. These parameters can be calibrated using the experimental results.

2.3. Stress–Strain Relationship

Concrete strain mainly includes elastic strain, creep strain, thermal strain and autogenous volumetric strain:
ε ˙ = ε ˙ e l + ε ˙ a s + ε ˙ t h + ε ˙ c r
where ε e l , ε c r , ε t h and ε a u are the elastic strain, creep strain, thermal strain and autogenous volumetric strain, respectively.
Accordingly, the stress evolution can be described as follows:
σ ˜ ˙ = E ( ξ ) ε ˙ e l = E ( ξ ) ( ε ˙ ε ˙ a s ε ˙ t h ε ˙ c r )

2.4. Mathematical Model of Expansion Evolution

2.4.1. Expansion Rate

The expansion of magnesium oxide concrete is essentially due to the hydration of magnesium oxide to form magnesium hydroxide, which follows the general rule of chemical reactions: temperature and concentration are the main factors that influence the reaction process. A higher temperature corresponds to a faster reaction, a higher hydration rate of magnesium oxide and more expansion. Meanwhile, the expansion rate is proportional to the magnesium oxide content involved in the reaction (i.e., the concentration in the liquid chemical reaction), which conforms to the kinetic chemical reaction equation (usually an Arrhenius equation).
At time t, the hydration rate of magnesium oxide can be expressed as a function of chemical reaction affinity and temperature at that time:
ζ t = A ˜ ζ ( ζ ) exp ( E m R T )
where ζ is the degree of hydration of magnesium oxide, Em is the chemical reaction activation energy and R is ideal gas constant.
Similar to the hydration process of cement [19], the chemical reaction affinity of magnesium oxide hydration can be obtained through a thermodynamic derivation:
A ˜ ζ ( ζ ) = B 1 ( B 2 ζ + ζ ) ( ζ - ζ ) exp ( η ¯ m ζ ζ )
where B1 and B2 are material constants, η ¯ m represents the viscosity due to the micro-diffusion of free water through the hydrates and ζ represents the ultimate hydration degree of magnesium oxide. The parameters Em, B1, B2, η ¯ m and ζ can be obtained via experimentation. Equation (8) can reasonably reflect the influences of the residual concentration ( ζ ζ ), temperature (T) and water diffusion rate ( exp ( η ¯ m ζ ζ ) ) on the magnesium oxide hydration.

2.4.2. Volumetric Strain

The autogenous volumetric strain of MgO concrete ε a u can be composed of the strain caused by the hydration of the cementing material and that of the magnesium oxide hydrating material. A linear or piecewise linear function can be used to describe the total volumetric strain:
ε ˙ a s = - ( k a ξ ˙ + k b ζ ˙ ) I ( ξ > ξ 0 )
where ka and kb are material parameters that can be obtained via experimentation, I represents a unit tensor and ξ 0 represents the threshold value of hydration (corresponding to the hydration degree of final condensation).

3. Numerical Implementation of the Thermo-Chemical Model

3.1. Hydration of Cementitious Material and Magnesium Oxide

The transient heat transfer equation and hydration equations are listed in (1), (3) and (7); to implement the coupled equations in the finite element method, an implicit algorithm of a backward-Euler temporal integration and a Gaussian quadrature spatial integration scheme is adopted. At time t + δ t , where δ t is the time increment, these equations can be discretized as follows:
ξ t + δ t t = ( ξ t + δ t ξ t ) δ t
T t + δ t t = ( T t + δ t T t ) δ t
ζ t + δ t t = ( ζ t + δ t ζ t ) δ t
Therefore, the hydration equation of the cementitious material can be written as
( ξ t + δ t ξ t ) δ t = A ξ ( ξ t + δ t ) exp ( E a R T t + δ t )
ξ t + δ t ξ t A ξ ( ξ t + δ t ) exp ( E a R T t + δ t ) δ t = 0
Similarly, the hydration equation of magnesium oxide can be written as
( ζ t + δ t ζ t ) δ t = A ζ ( ζ t + δ t ) exp ( E a R T t + δ t )
ζ t + δ t ζ t A ζ ( ζ t + δ t ) exp ( E a R T t + δ t ) δ t = 0
Based on the principle of virtual work, the basic energy balance equation for heat transfer can be written as
V ρ C T ( x , y , z , t ) t d V = S λ T Δ T ( x , y , z , t ) d S + V Q ξ t d V
Combining this equation with (10) and (11) gives
V ρ C ( T t + δ t T t ) δ t d V = S λ T Δ T ( x , y , z , t ) d S + V Q ( ξ t + δ t ξ t ) δ t d V
T t + δ t T t δ t V N T ρ C N d V = V λ T N T T ( x , y , z , t ) d V + ( V N T Q ( ξ t + δ t ξ t ) δ t d V )
T t + δ t T t δ t C t + δ t = K t + δ t T t + δ t + Q ( 1 δ t C + K ) T t + δ t = Q + 1 δ t C T t
where N is the shape–function matrix, the superscript T indicates the transpose of a vector or matrix, and K and C are the conduction and heat capacity matrix, respectively. At each time increment, the coupled Equations (11) and (20) are used to obtain the hydration and temperature fields. To form the Q matrix, ξ t + δ t and ζ t must be determined using Equations (14) and (16); thus, the Newton–Raphson iterative method is adopted. In an implicit time-integration algorithm, the nodal temperatures and hydration degrees of cementitious material and magnesium oxide can be solved as follows:
(1)
Input variables: δ t , T t , ξ t , ζ t , K t and C t ;
(2)
Set a proper temperature increment value and calculate T t + δ t ;
(3)
Use a Newton–Raphson iterative method to obtain the values of ξ t + δ t and ζ t + δ t :
ξ t + δ t k + 1 = ξ t + δ t k f ( ξ t + δ t k ) f ( ξ t + δ t k )
f ( ξ t + δ t k ) = ξ t + δ t k ξ t A ξ ( ξ t + δ t k ) exp ( E a R T t + δ t ) δ t
ζ t + δ t k + 1 = ζ t + δ t k f ( ζ t + δ t k ) f ( ζ t + δ t k )
f ( ζ t + δ t k ) = ζ t + δ t k ζ t A ζ ( ζ t + δ t k ) exp ( E a R T t + δ t ) δ t
where k is the iteration counter, f ( ξ t + δ t k ) and f ( ζ t + δ t k ) are the derivatives of f ( ξ t + δ t k ) and f ( ζ t + δ t k ) , respectively, ξ t + δ t 0 = ξ t , and ζ t + δ t 0 = ζ t ; this process continues while | ξ t + δ t k + 1 ξ t + δ t k ξ t + δ t k + 1 | > T o l and | ζ t + δ t k + 1 ζ t + δ t k ζ t + δ t k + 1 | > T o l ( T o l is the tolerance of the iteration).
(4)
Calculate the Q , K t + δ t and C t + δ t matrices and verify the balance of Equation (20). If a state of equilibrium can be reached for Equation (20), then proceed to step (5); otherwise, return to step (2) and correct the temperature increment.
(5)
Output variables T t + δ t , ξ t + δ t , ζ t + δ t , K t + δ t and C t + δ t .

3.2. Expansion and Stress Evolution of Magnesium Oxide Concrete

Through the above steps, the temperatures and hydration degrees of the cementitious material and magnesium oxide can be calculated at each moment. Combined with the relationship between the mechanical properties and the hydration degree of concrete, the stress and strain of concrete at each moment can be obtained. In fact, this method can calculate the temperature, stresses and hydration degree synchronously. However, considering that the stress results have little effect on temperature and hydration (i.e., the stress field is weakly coupled with the temperature and chemical fields), a sequential coupling method is essentially the same as a synchronous computation.
The relationship between the apparent stress and strain of concrete can be described as follows:
σ ˜ ˙ = E ε ˙ e l = E ( ε ˙ ε ˙ a s ε ˙ t h ε ˙ c r )
E is the elasticity matrix, which can be written as
E = E ( ξ ) [ 1 μ 1 μ μ 1 μ 0 0 0 μ 1 μ 1 μ 1 μ 0 0 0 μ 1 μ μ 1 μ 1 0 0 0 0 0 0 1 2 μ 2 ( 1 + μ ) 0 0 0 0 0 0 1 2 μ 2 ( 1 + μ ) 0 0 0 0 0 0 1 2 μ 2 ( 1 + μ ) ]
where E(ξ) is Young’s modulus and μ is Poisson’s ratio.
The evolution of Young’s modulus can be described as follows:
E ( ξ ) = ( ξ ξ 0 1 ξ 0 ) r e E
where E ( ξ ) and E are the Young’s moduli at degrees of hydration ξ and ξ , respectively, ξ 0 represents the hydration degree of the final set of concrete (the concrete completely loses its plasticity and the strength begins to increase) and ξ 0 and γ e are parameters that depend on the concrete composition.
The thermal strain of concrete can be written as
ε ˙ t h = - α T ˙ I
where α is the coefficient of linear expansion and I is the unit tensor {1,1,1,0,0,0}T.
Considering the influence of magnesium oxide hydration on concrete expansion, the strain produced by the autogenous deformation of concrete can be written as
{ ε ˙ a s = - k b ζ ˙ I ( ζ > ζ 0 ) ε ˙ a s = 0 ( ζ ζ 0 )
The creep strain can be calculated by using solidification theory. Therefore, the calculation of volume deformation, strain and the ultimate apparent stress of magnesium oxide concrete during the hydration process can be completed by using the above steps.

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:
ζ t ( T = T 1 ) ζ t ( T = T 2 ) = A ˜ ζ ( ζ ) exp ( E m R T 1 ) A ˜ ζ ( ζ ) exp ( E m R T 2 ) = exp [ E m R ( 1 T 1 1 T 2 ) ]
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
d ζ = A ˜ ζ ( ζ ) exp ( E m R T ) d t
d ζ A ˜ ζ ( ζ ) exp ( E m R T ) = d t
By integrating both sides of this equation,
d ζ A ˜ ζ ( ζ ) exp ( E m R T ) = d t
t = 0 ζ d ζ A ˜ ζ ( ζ ) exp ( E m R T )
Then, the time required to reach a certain degree of hydration under different temperatures can be described as:
t 1 = 0 ζ d ζ A ˜ ζ ( ζ ) exp ( E m R T 1 )
t 2 = 0 ζ d ζ A ˜ ζ ( ζ ) exp ( E m R T 2 )
t 1 t 2 = exp ( E m R T 2 ) exp ( E m R T 1 ) = exp [ E m R ( 1 T 2 1 T 1 ) ]
Ultimately, the chemical reaction activation energy of MgO hydration Em can be determined:
E m R = ln t 1 t 2 1 T 1 1 T 2
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 E a / R 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.

Author Contributions

Conceptualization, C.F., and C.Z.; Methodology, C.F.; Software, C.F., and C.Z.; Validation, C.F. and X.Y.; Formal analysis, C.F. and C.Z.; Investigation, C.F., C.Z. and J.X.; Resources, C.F., X.Y. and L.T.; Data curation, C.F. and C.Z.; Writing—original draft preparation, C.F.; Writing—review and editing, C.Z., X.Y., J.X. and L.T.; Visualization, X.Y. and J.X.; Supervision, X.Y.; Project administration, C.Z.; Funding acquisition, C.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China under Grant No. 51809200 and China postdoctoral science foundation under Grant No. 2018M633636XB.

Data Availability Statement

Some or all data, models or code that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Zhu, B. Thermal Stresses and Temperature Control of Mass Concrete; Butterworth-Heinemann: London, UK, 2013. [Google Scholar]
  2. Ulm, F.-J.; Coussy, O. Modeling of Thermochemomechanical Couplings of Concrete at Early Ages. J. Eng. Mech. 1995, 121, 785–794. [Google Scholar] [CrossRef]
  3. Liu, X.; Zhang, C.; Chang, X.; Zhou, W.; Cheng, Y.; Duan, Y. Precise simulation analysis of the thermal field in mass concrete with a pipe water cooling system. Appl. Therm. Eng. 2015, 78, 449–459. [Google Scholar] [CrossRef]
  4. Mo, L.; Deng, M.; Wang, A. Effects of MgO-based expansive additive on compensating the shrinkage of cement paste under non-wet curing conditions. Cem. Concr. Compos. 2012, 34, 377–383. [Google Scholar] [CrossRef]
  5. Deng, M.; Tang, M. Formation and expansion of ettringite crystals. Cem. Concr. Res. 1994, 24, 119–126. [Google Scholar]
  6. Mo, L.; Deng, M.; Tang, M.; Al-Tabbaa, A. MgO expansive cement and concrete in China: Past, present and future. Cem. Concr. Res. 2014, 57, 1–12. [Google Scholar] [CrossRef]
  7. Gao, P.-W.; Xu, S.-Y.; Chen, X.; Li, J.; Lu, X.-L. Research on autogenous volume deformation of concrete with MgO. Constr. Build. Mater. 2013, 40, 998–1001. [Google Scholar] [CrossRef]
  8. José, N.; Ahmed, H.; Miguel, B.; Luís, E.; Jorge, D.B. Magnesia (MgO) Production and Characterization, and Its Influence on the Performance of Cementitious Materials: A Review. Materials 2020, 13, 4752. [Google Scholar] [CrossRef] [PubMed]
  9. Yang, Y.; Li, C. Study on long term autogenous volume deformation of MgO concrete. J. Hydraul. Eng. 1999, 65, 54–58. [Google Scholar]
  10. Zhu, B. Incremental Type of Computing Model for the Volume Expansion of Concrete with Gentle Volume Expansion. Hydroelec. Eng. 2003, 29, 20–23. [Google Scholar]
  11. Zhu, B. Computational model and experimental method for bulk expansion of gentle expansive concrete. J. Hydraul. Eng. 2002, 33, 18–21. [Google Scholar]
  12. Yang, G.; Yuan, M. The hyperbola model for autogenous expansion volume deformation of MgO concrete. J. Hydroel. Eng. 2004, 4, 38–44. [Google Scholar]
  13. Zhang, G.; Chen, X.; Du, L. Analysis of expansion effect of magnesium oxide (MgO) in mass concrete. J. Hydraul. Eng. 2005, 2, 185–189. [Google Scholar]
  14. Zhang, G. MgO micro-expansion thermal integral model considering temperature process effect. Water Power 2002, 11, 28–32. [Google Scholar]
  15. Zhang, G.; Feng, J.; Luo, X.; Yang, B. Simulation of expansion process of MgO concrete considering the temperature history effect. J. Hydraul. Eng. 2002, 8, 29–34. [Google Scholar]
  16. Liu, S.H.; Fang, K.H. Study on Autogenous Deformation of Concrete Incorporating MgO as Expansive Agent. Damage Assess. Struct. VII 2006, 302, 155–161. [Google Scholar] [CrossRef]
  17. Xu, P.; Zhu, Y.; Ben, N. Computing model for autogenous volume deformation of MgO concrete based on degree of hy-dration. Water Resour. Hydrop. Eng. 2008, 2, 22–25. [Google Scholar]
  18. Cervera, M.; Oliver, J.; Prato, T. Thermo-Chemo-Mechanical Model for Concrete. I: Hydration and Aging. J. Eng. Mech. 1999, 125, 1018–1027. [Google Scholar] [CrossRef] [Green Version]
  19. Chen, C.; Feng, L.; Fang, K. Arc tangent curve model of the self-volume deformation of MgO concrete. J. Hydroel. Eng. 2008, 27, 106–110. [Google Scholar]
  20. Zhou, W.; Feng, C.; Liu, X.; Liu, S.; Zhang, C.; Yuan, W. Contrastive Numerical Investigations on Thermo-Structural Behaviors in Mass Concrete with Various Cements. Materials 2016, 9, 378. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Zhou, W.; Feng, C.; Liu, X.; Liu, S.; Zhang, C. A macro–meso chemo-physical analysis of early-age concrete based on a fixed hydration model. Mag. Concr. Res. 2016, 68, 981–994. [Google Scholar] [CrossRef]
  22. Lackner, R.; Mang, H.A. Chemoplastic material model for the simulation of early-age cracking: From the constitutive law to numerical analyses of massive concrete structures. Cem. Concr. Compos. 2004, 26, 551–562. [Google Scholar] [CrossRef]
  23. Cervera, M.; Faria, R.; Oliver, J.; Prato, T. Numerical modelling of concrete curing, regarding hydration and temperature phenomena. Comput. Struct. 2002, 80, 1511–1521. [Google Scholar] [CrossRef]
Figure 1. Temperature evolution of magnesium oxide concrete (fly ash content 30% and 40%) under adiabatic conditions.
Figure 1. Temperature evolution of magnesium oxide concrete (fly ash content 30% and 40%) under adiabatic conditions.
Materials 14 03162 g001
Figure 2. Comparison of the simulation and experimental values of the autogenous volumetric strain of magnesium oxide concrete with various curing temperatures (fly ash content 30% and 40%).
Figure 2. Comparison of the simulation and experimental values of the autogenous volumetric strain of magnesium oxide concrete with various curing temperatures (fly ash content 30% and 40%).
Materials 14 03162 g002
Figure 3. Young’s modulus evolution of magnesium oxide concrete (fly ash content 30% and 40%) under adiabatic conditions.
Figure 3. Young’s modulus evolution of magnesium oxide concrete (fly ash content 30% and 40%) under adiabatic conditions.
Materials 14 03162 g003
Figure 4. Coupled thermo-chemical-mechanical field computations on 1 March 2003. (a) Temperature (°C); (b) hydration degree of cementitious material (-); (c) hydration degree of magnesium oxide (-); (d) first principal stress (MPa).
Figure 4. Coupled thermo-chemical-mechanical field computations on 1 March 2003. (a) Temperature (°C); (b) hydration degree of cementitious material (-); (c) hydration degree of magnesium oxide (-); (d) first principal stress (MPa).
Materials 14 03162 g004
Figure 5. Coupled thermo-chemical-mechanical field computations on 1 June 2003. (a) Temperature (°C); (b) hydration degree of cementitious material (-); (c) hydration degree of magnesium oxide (-); (d) first principal stress (MPa).
Figure 5. Coupled thermo-chemical-mechanical field computations on 1 June 2003. (a) Temperature (°C); (b) hydration degree of cementitious material (-); (c) hydration degree of magnesium oxide (-); (d) first principal stress (MPa).
Materials 14 03162 g005
Figure 6. Position distribution of typical points.
Figure 6. Position distribution of typical points.
Materials 14 03162 g006
Figure 7. The temperature histories of different typical points.
Figure 7. The temperature histories of different typical points.
Materials 14 03162 g007
Figure 8. Max principal stress histories of different typical points.
Figure 8. Max principal stress histories of different typical points.
Materials 14 03162 g008
Figure 9. Cement hydration degree histories of different typical points.
Figure 9. Cement hydration degree histories of different typical points.
Materials 14 03162 g009
Figure 10. Magnesium oxide hydration degree histories of different typical points.
Figure 10. Magnesium oxide hydration degree histories of different typical points.
Materials 14 03162 g010
Table 1. Component mass of oxygen magnesia concrete samples with varying fly ash contents.
Table 1. Component mass of oxygen magnesia concrete samples with varying fly ash contents.
Fly Ash ContentComposition of Concrete (kg/m3)Set-Retarding
Superplasticizer
w/bWaterCementFly AshMgOSandStone (mm)
5–2020–4040–8080–150
0.5730%100123528.85413433435145141.22
0.5740%100105708.85413433435145141.22
Table 2. Adiabatic temperature increases of oxygen magnesia concrete samples with varying water–binder ratios.
Table 2. Adiabatic temperature increases of oxygen magnesia concrete samples with varying water–binder ratios.
Fly Ash Content1 d3 d7 d14 d21 d28 d45 d60 d90 d
30%3.638.6414.2418.8121.0622.4024.1524.9525.80
40%4.9310.3815.1818.3619.7320.5021.4521.8722.30
Table 3. Autogenous volume deformation of oxygen magnesia concrete samples at different ages with varying water–binder ratios and curing temperatures.
Table 3. Autogenous volume deformation of oxygen magnesia concrete samples at different ages with varying water–binder ratios and curing temperatures.
Fly Ash ContentCuring TemperatureAutogenous Volumetric Strain (10−6)
1 d3 d7 d14 d28 d60 d90 d120 d150 d180 d
30%20 °C2.97.316.127.148.166.77783.790.694.9
40%20 °C3.47.516.123.63858.366.673.582.888.8
30%30 °C7.324.144.36483.2102116.4125.2129.1132
40%30 °C17.430.338.951.971.693.5108120125130.3
Table 4. Parameters describing the cement hydration properties.
Table 4. Parameters describing the cement hydration properties.
Parameters β 1 (106 h−1) β 2 (10−2) ξ (-) η ¯ (-) Q (J kg−1) E a / R (K)
Fly ash content 30%1.885.616.626,3135000
Fly ash content 40%4.089.016.822,0625000
Table 5. Parameters for describing the MgO hydration properties.
Table 5. Parameters for describing the MgO hydration properties.
Parameters B 1 (109 h−1) B 2 (104) ζ (-) η ¯ m (-) E m / R (K)kakb
Fly ash content 30%15.21014.012,0000204
Fly ash content 40%4.52.814.511,5000204
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Feng, C.; Zhao, C.; Yu, X.; Xiong, J.; Tang, L. A Mathematical Model of the Expansion Evolution of Magnesium Oxide in Mass Concrete Based on Hydration Characteristics. Materials 2021, 14, 3162. https://doi.org/10.3390/ma14123162

AMA Style

Feng C, Zhao C, Yu X, Xiong J, Tang L. A Mathematical Model of the Expansion Evolution of Magnesium Oxide in Mass Concrete Based on Hydration Characteristics. Materials. 2021; 14(12):3162. https://doi.org/10.3390/ma14123162

Chicago/Turabian Style

Feng, Chuqiao, Cheng Zhao, Xiaomin Yu, Jie Xiong, and Longwen Tang. 2021. "A Mathematical Model of the Expansion Evolution of Magnesium Oxide in Mass Concrete Based on Hydration Characteristics" Materials 14, no. 12: 3162. https://doi.org/10.3390/ma14123162

APA Style

Feng, C., Zhao, C., Yu, X., Xiong, J., & Tang, L. (2021). A Mathematical Model of the Expansion Evolution of Magnesium Oxide in Mass Concrete Based on Hydration Characteristics. Materials, 14(12), 3162. https://doi.org/10.3390/ma14123162

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop