# A New Safety Factor Prediction Model for Mass Concrete Surface Cracking in Early Age

Authors: Sheng Qiang; Min-zhi Liu
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

The temperature difference between mass concrete inside and outside is an important index to prevent the structure surface cracking, but the definition and the application criterion are not very clear till now. In this paper, the surface concrete temperature spatial gradient in 0.5 m deep is found as a more reasonable and feasible index. Then a new model, which can be used to calculate the surface safety factor simply, is deduced in three steps. First, a basic model is created under the simplest boundary condition. Then the feasibility of the basic model is analyzed under different boundary conditions and the model is improved. At last, according to the comparison of different concrete materials, the improved model is calibrated further. With the help of the proposed model, the surface safety factor can be adjusted by observing and controlling the temperature spatial gradient dynamically, which will bring convenience to designer and constructor.


---

## Body

## 1. Introduction

One of the hard tasks of concrete structure in research and application is cracking prevention or mitigation. In the running period, concrete structure tends to become damaged and cracked because of some extreme conditions, such as fire or freeze-thaw, which may cause obvious strength reduction [1, 2]. On the other hand, many cracks usually tend to form in the early age of concrete [3, 4], even under well-designed construction conditions. The early age cracking reasons include many factors, such as humidity gradient, autogenous shrinkage, temperature gradient, structure restraint, and casting block shape and size. According to the observation on some concrete dams, the surface cracks of mass concrete caused by humidity gradient are not severe as thin concrete structure. The material researchers have created many methods to reduce the drying and autogenous shrinkage [5, 6]. The structure and construction researchers of mass concrete put more emphasis on the later three factors, especially the temperature gradient [7, 8]. Some effective temperature control methods have been implemented in mass concrete structure, such as pipe cooling, surface covering, and low casting temperature. Pipe cooling can decrease the inside temperature of mass concrete. Some researchers proposed pipe simulation models to predict the inside temperature and stress, which increase the reliability of inner cracking prevention of mass concrete [9–11].The surface cracking of concrete in early age is related to the temperature difference between inside and outside. In fact, as the index of the temperature control and crack prevention, the control of the temperature difference between inside and outside is the control of the influence of concrete temperature spatial gradient on concrete thermal stress. As for some mass concrete, especially the concrete dam pouring blocks which are not in the base restraint area, the inside temperature in early age may be high, but the inside stress may not be so high because its temperature spatial gradient is not large. For the surface concrete, the large temperature spatial gradient always leads to large tensile stress which will cause fracture [12–15]. The reason is the inharmonious deformation in different areas.Although the temperature difference between inside and outside is an important index, the definition and the application criterion are not very clear till now, even in the design standard or construction standard. According to the experience of author, the definition and the criterion are correlative. Because of the unclear definition, the different criteria suggested by different researchers for a same object always cause puzzle in engineering application. This paper tries to deduce a new model which is simple to solve such problem.In the published literatures, it is difficult to find a convenient and scientific prediction model for surface cracking control of mass concrete. A common used model in US Bureau of Reclamation to calculate the stress caused by temperature gradient is recommended [16]. The model can predict the stress caused by the temperature gradient in a certain thickness. An example with some data was shown in a table of this book. In this sample, some basic data and an assumed temperature distribution along concrete thickness was shown in a table; then the stresses were computed in the table with the recommended model. We plot the temperature and stress of the sample in Figures 1 and 2. Figure 1 shows the relationship of temperature and distance to surface, and Figure 2 shows the relationship of stress and distance to surface. It can be found that the stress distribution near surface changes from tensile stress to pressure stress with the temperature increasing, which is similar to the finite element solution. However, the stress inside concrete becomes from pressure stress to tensile stress with the temperature keeping increasing, which is sharply different with finite element solution. The common solution of nowadays is that, in the early age, when the surface temperature is lower than inside, tensile stress will occur near surface and pressure stress will occur inside, while in the later age, when the surface temperature is higher than inside, pressure stress will occur near surface and tensile stress will occur inside. On the other hand, the assumed temperature gradient in this example can only describe the mature age. Actually, in the early age, the temperature gradient near surface will be much greater than inside and it should be nonlinear as well.Figure 1
Assumed relation of temperature and distance to surface.Figure 2
Calculated relation of stress and distance to surface.In order to control the tensile stress developing in early age on concrete surface, a cracking safety factorK can be defined as the ratio of tensile strength to tensile stress. This definition means when the K is lower than 1.0, the obvious crack can be found by human eye.As for the common situation, the main influence factors of the surface cracking safety for mass concrete in early age include concrete elastic module, tension strength, and temperature spatial gradient in surface layer.The relation between temperature spatial gradient and surface cracking safety in different situations is researched in this paper. The results show that the safety factor will decrease when the concrete elastic module and the temperature spatial gradient increase, while it will increase when the concrete tension strength increases. The last proposed model is feasible to control the surface safety by observing and controlling the temperature spatial gradient, which will make it easier for the designer or constructor to get the reasonable index of temperature spatial gradient.The deduction of the new model can be separated into three steps which will be shown in the following paragraphs. At first, the basic model is created under the simplest boundary condition. Then the feasibility of the basic model is analyzed under different boundary conditions and the model is improved. At last, according to the comparison of different concrete materials, the improved model is calibrated further and then the practical formula of the concrete surface safety factor model is proposed.


## 2. Deduction of Basic Model

### 2.1. The Research Object

After the calculation of a concrete pouring block in 60.0 days with finite element method (FEM), the temperature field in the 1.0 m deep area near surface is gotten. Figure3 shows the temperature development curves of the survey points with different distances from surface during 60.0 days. The distance of the 16 points to surface varies from 0.00 m to 1.00 m with an increment of 0.05 m. It can be seen form Figure 3 that the deeper points have higher temperature and lower temperature spatial gradient. Figure 4 shows the temperature contours near surface in different ages. It can be seen when the age is over 10.0 d the temperature spatial gradient in the surface of 1.0 m deep area is nearly same; else the temperature field in the area of only 0.5 m deep changes evenly. Therefore, the research object of this paper will be selected as the area of 0.5 m deep.Figure 3
The temperature development curves of the survey points with different distances from surface in 60.0 days.Temperature contours in 1.0 m deep area at different ages (°C).
(a)
(b)
(c)
(d)
(e)
(f)


### 2.2. The Finite Element Test Model and Calculation Case Series 1

The finite element test model is a large-scale concrete block with the dimension of 40 m × 20 m × 20 m on a base, which is shown in Figure5. In order to improve the calculation precision, the surface concrete mesh is refined.Figure 5
The finite element test model.The main material parameters of the concrete can be found in Table1. During the deduction course of the basic formula, some factors are not considered to simplify the boundary conditions, including the daily variation of air temperature, pipe cooling, and surface insulation measures. The influence of these factors will be analyzed further in the later paragraphs.Table 1
Main concrete material parameters in basic model deduction.
Material parameters
Adiabatic temperature rising (°C)θ(τ)=T0×(1-e-aτb)
Elastic module (GPa)E(τ)=E0×(1-e-aτb)
Tensile strength (MPa)ft(τ)=f0×(1-e-aτb)
θ
0
a
b
E
0
a
b
f
t
a
b
20.0
0.50
1.00
35.0
0.75
0.50
3.25
0.75
0.50Table2 shows 6 cases of series 1 with the same initial temperature 20°C, in which the air temperature is changed to form the different temperature spatial gradient in the surface concrete.Table 2
The detailed condition of case series 1.
Case
Air temperatureTa (°C)
Pouring temperatureT0 (°C)
Note
1-1
30.0
20.0
No pipe cooling, no surface insulation, no daily variation of air temperature
1-2
25.0
1-3
20.0
1-4
15.0
1-5
10.0
1-6
5.0The three-dimensional finite element method for unsteady temperature field and stress field [9] are applied in the construction simulation for the test block.


### 2.3. Deduction of the Basic Formula

On the basis of calculation result of case series 1, the relation between the surface safety factorK and concrete temperature spatial gradient ∂T/∂l at different ages is described in Figure 6. The safety factor K equals σt/σ1, where σ1 is the first principle stress of surface point and σt is the tensile strength at the same age. ∂T/∂l is the average temperature spatial gradient in the 0.5 m deep area, which can be measured easily by the temperature sensors in concrete. The following conclusion can be gotten according to Figure 6.(1)
With the growing of concrete age, the temperature spatial gradient increases first and then decreases. The concrete inside temperature keeps rising after concrete pouring because of hydration reaction, while the air temperature is assumed to be steady. So the temperature spatial gradient keeps rising. When the inside temperature arrives at the peak, it begins to drop. So temperature spatial gradient begins to decrease, too.(2)
With the growing of concrete age, the concrete surface safety factor decreased first and then increases a little. At the beginning of the age, the concrete elastic module and the stress are very small, so the surface safety factor is of large value. The growth of both the elastic module and the temperature spatial gradient induces the thermal stress increase rapidly, while the tensile strength increases relatively slowly. So the safety factor goes down fast. With the developing of concrete strength and the dropping of temperature spatial gradient, the safety factor grows gently in the later age.(3)
In every case, the relation curves between the temperature spatial gradient and the safety factor have the regular rule. The curves’ shape is very similar and shows significant grading rule with the change of the air temperature.Relation between safety factor and temperature spatial gradient in case series 1.
(a)
(b)Table3 shows the safety factor and temperature spatial gradient at age 10 d in every case. The data in Table 3 can be fitted by a curve which can be described by different functions well. To simplify the new model, the power function is adopted. The fit line can be seen in Figure 7.Table 3
Safety factor and temperature spatial gradient at age 10 d in case series 1.
Case number
Case  1-1
Case  1-2
Case  1-3
Case  1-4
Case  1-5
Case  1-6
∂
T
/
∂
l (°C/m)
6.562
9.387
12.209
15.033
23.964
26.787
K
1.906
1.526
1.258
1.070
0.656
0.601Figure 7
Relation between safety factor and temperature spatial gradient at age 10 d in case series 1.From Figure7 the temperature spatial gradient function f(∂T/∂l) can be expressed by
(1)f(∂T∂l)=9.9166×(∂T∂l)-0.8438.Equation (1) is based on the relation curve of age 10 d only. Then the age factor will be introduced into the basic formula. If the relation curves as Figure 7 at every calculation step are drawn, the safety factor K corresponding to the temperature spatial gradient of 10°C/m can be interpolated. Then these K values can be divided by the K value of age 10 d; a standardized safety factor K′ will be obtained. The values of the interpolated K and the standardized K′ at different ages are shown in Table 4. Based on Table 4, the relation between age and K′ can be described by another power function, which is shown in Figure 8.Table 4
The interpolatedK and standardized K′ at different ages.
Age (d)
0.25
0.5
0.75
1
1.25
1.5
1.75
2
2.25
2.5
2.75
3
K
/
2.80
2.87
2.68
2.65
2.52
2.39
2.30
2.25
2.18
2.11
2.05
K
′
1.91
1.96
1.83
1.81
1.72
1.63
1.57
1.53
1.49
1.44
1.40
Age (d)
3.5
4
4.5
5
5.5
6
6.5
7
8
9
10
12
K
1.94
1.86
1.81
1.77
1.72
1.68
1.65
1.62
1.56
1.51
1.47
1.41
K
′
1.32
1.27
1.24
1.20
1.17
1.15
1.12
1.10
1.06
1.03
1.00
0.96
Age (d)
14
16
20
25
30
35
40
45
50
55
60
K
1.35
1.30
1.24
1.18
1.14
1.10
1.06
1.02
0.98
0.94
0.90
K
′
0.92
0.88
0.85
0.80
0.78
0.75
0.72
0.70
0.67
0.64
0.62Figure 8
Relation between age and standardize safety factorK′ when temperature spatial gradient is 10°C/m in case series 1.From Figure8, the expression g(τ) to describe the age factor can be expressed by
(2)g(τ)=1.8245×τ-0.2557.Multiplying (1) and (2), the basic formula (3) can be obtained, which is based on assuming the same developing law of concrete elastic module and tensile strength:
(3)K(τ)=f(∂T∂l)g(τ).Comparing the calculated safety factorK with (3) and with finite element method in Figure 9, it can be seen that they have a good agreement in most areas. To show the curves of the 6 cases more clearly, they are separately shown in two plots in Figure 9.The safety factor comparison between two methods. (a) Results of Cases 1-1, 1-3, and 1-5. (b) Results of Cases 1-2, 1-4, and 1-6.
(a)
(b)


## 2.1. The Research Object

After the calculation of a concrete pouring block in 60.0 days with finite element method (FEM), the temperature field in the 1.0 m deep area near surface is gotten. Figure3 shows the temperature development curves of the survey points with different distances from surface during 60.0 days. The distance of the 16 points to surface varies from 0.00 m to 1.00 m with an increment of 0.05 m. It can be seen form Figure 3 that the deeper points have higher temperature and lower temperature spatial gradient. Figure 4 shows the temperature contours near surface in different ages. It can be seen when the age is over 10.0 d the temperature spatial gradient in the surface of 1.0 m deep area is nearly same; else the temperature field in the area of only 0.5 m deep changes evenly. Therefore, the research object of this paper will be selected as the area of 0.5 m deep.Figure 3
The temperature development curves of the survey points with different distances from surface in 60.0 days.Temperature contours in 1.0 m deep area at different ages (°C).
(a)
(b)
(c)
(d)
(e)
(f)


## 2.2. The Finite Element Test Model and Calculation Case Series 1

The finite element test model is a large-scale concrete block with the dimension of 40 m × 20 m × 20 m on a base, which is shown in Figure5. In order to improve the calculation precision, the surface concrete mesh is refined.Figure 5
The finite element test model.The main material parameters of the concrete can be found in Table1. During the deduction course of the basic formula, some factors are not considered to simplify the boundary conditions, including the daily variation of air temperature, pipe cooling, and surface insulation measures. The influence of these factors will be analyzed further in the later paragraphs.Table 1
Main concrete material parameters in basic model deduction.
Material parameters
Adiabatic temperature rising (°C)θ(τ)=T0×(1-e-aτb)
Elastic module (GPa)E(τ)=E0×(1-e-aτb)
Tensile strength (MPa)ft(τ)=f0×(1-e-aτb)
θ
0
a
b
E
0
a
b
f
t
a
b
20.0
0.50
1.00
35.0
0.75
0.50
3.25
0.75
0.50Table2 shows 6 cases of series 1 with the same initial temperature 20°C, in which the air temperature is changed to form the different temperature spatial gradient in the surface concrete.Table 2
The detailed condition of case series 1.
Case
Air temperatureTa (°C)
Pouring temperatureT0 (°C)
Note
1-1
30.0
20.0
No pipe cooling, no surface insulation, no daily variation of air temperature
1-2
25.0
1-3
20.0
1-4
15.0
1-5
10.0
1-6
5.0The three-dimensional finite element method for unsteady temperature field and stress field [9] are applied in the construction simulation for the test block.


## 2.3. Deduction of the Basic Formula

On the basis of calculation result of case series 1, the relation between the surface safety factorK and concrete temperature spatial gradient ∂T/∂l at different ages is described in Figure 6. The safety factor K equals σt/σ1, where σ1 is the first principle stress of surface point and σt is the tensile strength at the same age. ∂T/∂l is the average temperature spatial gradient in the 0.5 m deep area, which can be measured easily by the temperature sensors in concrete. The following conclusion can be gotten according to Figure 6.(1)
With the growing of concrete age, the temperature spatial gradient increases first and then decreases. The concrete inside temperature keeps rising after concrete pouring because of hydration reaction, while the air temperature is assumed to be steady. So the temperature spatial gradient keeps rising. When the inside temperature arrives at the peak, it begins to drop. So temperature spatial gradient begins to decrease, too.(2)
With the growing of concrete age, the concrete surface safety factor decreased first and then increases a little. At the beginning of the age, the concrete elastic module and the stress are very small, so the surface safety factor is of large value. The growth of both the elastic module and the temperature spatial gradient induces the thermal stress increase rapidly, while the tensile strength increases relatively slowly. So the safety factor goes down fast. With the developing of concrete strength and the dropping of temperature spatial gradient, the safety factor grows gently in the later age.(3)
In every case, the relation curves between the temperature spatial gradient and the safety factor have the regular rule. The curves’ shape is very similar and shows significant grading rule with the change of the air temperature.Relation between safety factor and temperature spatial gradient in case series 1.
(a)
(b)Table3 shows the safety factor and temperature spatial gradient at age 10 d in every case. The data in Table 3 can be fitted by a curve which can be described by different functions well. To simplify the new model, the power function is adopted. The fit line can be seen in Figure 7.Table 3
Safety factor and temperature spatial gradient at age 10 d in case series 1.
Case number
Case  1-1
Case  1-2
Case  1-3
Case  1-4
Case  1-5
Case  1-6
∂
T
/
∂
l (°C/m)
6.562
9.387
12.209
15.033
23.964
26.787
K
1.906
1.526
1.258
1.070
0.656
0.601Figure 7
Relation between safety factor and temperature spatial gradient at age 10 d in case series 1.From Figure7 the temperature spatial gradient function f(∂T/∂l) can be expressed by
(1)f(∂T∂l)=9.9166×(∂T∂l)-0.8438.Equation (1) is based on the relation curve of age 10 d only. Then the age factor will be introduced into the basic formula. If the relation curves as Figure 7 at every calculation step are drawn, the safety factor K corresponding to the temperature spatial gradient of 10°C/m can be interpolated. Then these K values can be divided by the K value of age 10 d; a standardized safety factor K′ will be obtained. The values of the interpolated K and the standardized K′ at different ages are shown in Table 4. Based on Table 4, the relation between age and K′ can be described by another power function, which is shown in Figure 8.Table 4
The interpolatedK and standardized K′ at different ages.
Age (d)
0.25
0.5
0.75
1
1.25
1.5
1.75
2
2.25
2.5
2.75
3
K
/
2.80
2.87
2.68
2.65
2.52
2.39
2.30
2.25
2.18
2.11
2.05
K
′
1.91
1.96
1.83
1.81
1.72
1.63
1.57
1.53
1.49
1.44
1.40
Age (d)
3.5
4
4.5
5
5.5
6
6.5
7
8
9
10
12
K
1.94
1.86
1.81
1.77
1.72
1.68
1.65
1.62
1.56
1.51
1.47
1.41
K
′
1.32
1.27
1.24
1.20
1.17
1.15
1.12
1.10
1.06
1.03
1.00
0.96
Age (d)
14
16
20
25
30
35
40
45
50
55
60
K
1.35
1.30
1.24
1.18
1.14
1.10
1.06
1.02
0.98
0.94
0.90
K
′
0.92
0.88
0.85
0.80
0.78
0.75
0.72
0.70
0.67
0.64
0.62Figure 8
Relation between age and standardize safety factorK′ when temperature spatial gradient is 10°C/m in case series 1.From Figure8, the expression g(τ) to describe the age factor can be expressed by
(2)g(τ)=1.8245×τ-0.2557.Multiplying (1) and (2), the basic formula (3) can be obtained, which is based on assuming the same developing law of concrete elastic module and tensile strength:
(3)K(τ)=f(∂T∂l)g(τ).Comparing the calculated safety factorK with (3) and with finite element method in Figure 9, it can be seen that they have a good agreement in most areas. To show the curves of the 6 cases more clearly, they are separately shown in two plots in Figure 9.The safety factor comparison between two methods. (a) Results of Cases 1-1, 1-3, and 1-5. (b) Results of Cases 1-2, 1-4, and 1-6.
(a)
(b)


## 3. The Feasibility of the Basic Formula under Different Boundary Conditions

### 3.1. Considering Daily Variation of Air Temperature

In the practical site situation, the air temperature always keeps changing every day and night. A daily variation of air temperature of about 10°C is considered in case series 2 which is based on case series 1. The detailed conditions can be seen in Table5.Table 5
Calculation condition in case series 2.
Case
Air temperatureTa (°C)
Note
2-1
25
±
5
The other condition is same as case 1
2-2
20
±
5
2-3
15
±
5
2-4
10
±
5Figure10 shows the influence of daily variation of air temperature on the surface layer concrete. It can be seen that although the temperature of the point in 0.5 m deep has nearly no daily wave, the temperature of the point in 0.05 m deep keeps changing with the air temperature. Consequently, the temperature spatial gradient defined in the basic model will be influenced by the daily variation of air temperature. The surface safety factor will be influenced accordingly.Figure 10
Surface layer concrete temperature change with the daily variation of air temperature.Figure11 shows the relation between concrete surface safety factor and temperature spatial gradient in case series 2. In this case series, the time step is 0.25 d. When the calculation result in one case is observed separately by 0.0 d, 0.25 d, 0.50 d, and 0.75 d in every day, it can be found in Figure 12 that the maximum temperature spatial gradient or the minimum safety factor always appears at the time of 0.25 d every day, which means this time is the most dangerous moment. For this reason, if the concrete at this time does not crack, then at other times the concrete will be safe enough. So the calculation result of 0.25 d in every day from every case of series 2 is extracted to be studied further in Figure 13.Figure 11
Relation between concrete surface safety factor and temperature spatial gradient in case series 2.Figure 12
Relation between concrete surface safety factor and temperature spatial gradient at different times of every day in case 2-1.Figure 13
Relation between minimum surface safety factor and temperature spatial gradient in case series 2.The minimum safety factor of case series 2 calculated by FEM and (3) is compared in Figure 14. It can be seen that they meet well in most areas. In order to increase the precision further, a parameter α is introduced into the model. The new formula is expressed by
(4)K(τ)=αf(∂T∂l)g(τ).Figure 14
Comparison between calculated values by FEM and (3).Figure15 shows that when α equals 1.05, the calculated values by FEM and (4) meet better. In fact, the value by (3) is a little less than FEM which will cause a little more safety measure in engineering. Therefore, it will be safer to set the value of α between 1.0 and 1.05.Figure 15
Comparison between calculated values by FEM and (4).


### 3.2. Considering Pipe Cooling

Cooling pipe is used widely in the construction of mass concrete [17, 18]; will it influence obviously the calculation result of the proposed model? It is necessary to research this feasibility further.The unsteady temperature field including cooling pipe is always a hard problem. An explicit cooling pipe iterative algorithm [19] is applied in the simulation of this paper.Insert four pipes into the finite element test model. The inlet and outlet of pipes can be seen in the local mesh on the left of Figure16 and the three-dimensional pipelines on the right. The distance between concrete surface and pipe is 1.0 m. Both of the horizontal and vertical distances of pipelines are 2.0 m. The pipe is adopted as one of the most common types in concrete engineering with an inner diameter of 0.028 m and the material of high-density polyethylene.The layout of cooling pipe in the test model. (a) Local mesh section with pipe inlets and outlets. (b) Three-dimensional pipeline layout.
(a)
(b)Based on case series 1, pipe cooling measure is added in case series 3 with the inlet water temperature of 15.0°C. In order to analyze the influence of pipe cooling separately, the daily variation of air temperature will not be considered in case series 3. The detailed condition is shown in Table6.Table 6
Calculation condition in case series 3.
Case
Air temperatureTa (°C)
Pipe inlet water temperatureTw (°C)
Note
3-1
30.0
15.0
The other condition is same as case 1
3-2
25.0
3-3
20.0
3-4
15.0
3-5
10.0
3-6
5.0Figure17 shows the comparison between the calculation results of case series 3 by FEM and (3). Figure 18 shows the comparison between the calculation results of case series 3 by FEM and (4) in which α equals 1.10.Comparison of case series 3 results between FEM and (3).
(a)
(b)Comparison of case series 3 results between FEM and (4).
(a)
(b)


### 3.3. Considering Both Daily Variation of Air Temperature and Pipe Cooling

In order to meet the requirement of practical situation, daily variation of air temperature and pipe cooling should be considered in the test block at the same time. In case series 4, pipe cooling measure is added on the basis of case series 2 calculation condition. The detailed condition can be seen in Table7.Table 7
Calculation condition of case series 4.
Case
Air temperatureTa (°C)
Pipe inlet water temperatureTw (°C)
Note
4-1
25
.
0
±
5
.
0
15.0
The other calculation condition is same as case series 1.
4-2
20
.
0
±
5
.
0
4-3
15
.
0
±
5
.
0
4-4
10
.
0
±
5
.
0Figure19 shows the relation between the temperature spatial gradient and the surface safety factor at different ages in case series 4. To observe it clearly, the results of the four times in every day of case 4-1 are drawn separately in Figure 20. It can be found that the temperature spatial gradient of 0.25 d is maximum and the safety factor of the time is minimum. Therefore, the results of 0.25 d in every case are drawn in Figure 21; the calculated value of (3) is shown as well. From this figure, the results by FEM and (3) do not meet very well, especially in the later age. While in Figure 22, the results become better when (4) is applied and α equals 1.10.Figure 19
Relation between concrete surface safety factor and temperature spatial gradient in case series 4.Figure 20
Relation between concrete surface safety factor and temperature spatial gradient at different times of every day in case 4-1.Figure 21
Comparison of case series 4 results at 0.25 d every day between FEM and (3).Figure 22
Comparison of case series 4 results at 0.25 d every day between FEM and (4).


### 3.4. Considering Both Daily Variation of Air Temperature, Pipe Cooling, and Surface Insulation

In winter or cold wave time, surface insulation measure is always taken to protect the concrete in early age. In case series 5, a layer of 2.0 cm thick expandable polyethylene (EPE), which is one of the most common application concrete insulation materials in construction, is added on the basis of case series 4 calculation condition. The detailed calculation condition is shown in Table8.Table 8
Calculation condition of case series 5.
Case
Air temperatureTa (°C)
Pipe inlet water temperatureTw (°C)
Surface insulation
Note
5-1
25
.
0
±
5
.
0
15.0
2.0 cm EPE material on all free face
The other condition is same as case series 1
5-2
20
.
0
±
5
.
0
5-3
15
.
0
±
5
.
0
5-4
10
.
0
±
5
.
0Figure23 shows the relation between the temperature spatial gradient and the surface safety factor at different ages in case series 5. To observe it clearly, the results of the four times in every day of case 5-1 are drawn separately in Figure 24. Being similar to case 4-1, the most dangerous situation appears at 0.25 d of every day. Therefore, the results of 0.25 d in every case are drawn in Figure 25; the calculated value of (3) is shown as well. From the figure, the results by FEM and (3) have obvious difference. When (4) is applied and α equals 1.15, the results become much closer, as shown in Figure 26.Figure 23
Relation between concrete surface safety factor and temperature spatial gradient in case series 5.Figure 24
Relation between concrete surface safety factor and temperature spatial gradient at different times of every day in case 5-1.Figure 25
Comparison of case series 4 results at 0.25 d every day between FEM and (3).Figure 26
Comparison of case series 5 results at 0.25 d every day between FEM and (4).


## 3.1. Considering Daily Variation of Air Temperature

In the practical site situation, the air temperature always keeps changing every day and night. A daily variation of air temperature of about 10°C is considered in case series 2 which is based on case series 1. The detailed conditions can be seen in Table5.Table 5
Calculation condition in case series 2.
Case
Air temperatureTa (°C)
Note
2-1
25
±
5
The other condition is same as case 1
2-2
20
±
5
2-3
15
±
5
2-4
10
±
5Figure10 shows the influence of daily variation of air temperature on the surface layer concrete. It can be seen that although the temperature of the point in 0.5 m deep has nearly no daily wave, the temperature of the point in 0.05 m deep keeps changing with the air temperature. Consequently, the temperature spatial gradient defined in the basic model will be influenced by the daily variation of air temperature. The surface safety factor will be influenced accordingly.Figure 10
Surface layer concrete temperature change with the daily variation of air temperature.Figure11 shows the relation between concrete surface safety factor and temperature spatial gradient in case series 2. In this case series, the time step is 0.25 d. When the calculation result in one case is observed separately by 0.0 d, 0.25 d, 0.50 d, and 0.75 d in every day, it can be found in Figure 12 that the maximum temperature spatial gradient or the minimum safety factor always appears at the time of 0.25 d every day, which means this time is the most dangerous moment. For this reason, if the concrete at this time does not crack, then at other times the concrete will be safe enough. So the calculation result of 0.25 d in every day from every case of series 2 is extracted to be studied further in Figure 13.Figure 11
Relation between concrete surface safety factor and temperature spatial gradient in case series 2.Figure 12
Relation between concrete surface safety factor and temperature spatial gradient at different times of every day in case 2-1.Figure 13
Relation between minimum surface safety factor and temperature spatial gradient in case series 2.The minimum safety factor of case series 2 calculated by FEM and (3) is compared in Figure 14. It can be seen that they meet well in most areas. In order to increase the precision further, a parameter α is introduced into the model. The new formula is expressed by
(4)K(τ)=αf(∂T∂l)g(τ).Figure 14
Comparison between calculated values by FEM and (3).Figure15 shows that when α equals 1.05, the calculated values by FEM and (4) meet better. In fact, the value by (3) is a little less than FEM which will cause a little more safety measure in engineering. Therefore, it will be safer to set the value of α between 1.0 and 1.05.Figure 15
Comparison between calculated values by FEM and (4).


## 3.2. Considering Pipe Cooling

Cooling pipe is used widely in the construction of mass concrete [17, 18]; will it influence obviously the calculation result of the proposed model? It is necessary to research this feasibility further.The unsteady temperature field including cooling pipe is always a hard problem. An explicit cooling pipe iterative algorithm [19] is applied in the simulation of this paper.Insert four pipes into the finite element test model. The inlet and outlet of pipes can be seen in the local mesh on the left of Figure16 and the three-dimensional pipelines on the right. The distance between concrete surface and pipe is 1.0 m. Both of the horizontal and vertical distances of pipelines are 2.0 m. The pipe is adopted as one of the most common types in concrete engineering with an inner diameter of 0.028 m and the material of high-density polyethylene.The layout of cooling pipe in the test model. (a) Local mesh section with pipe inlets and outlets. (b) Three-dimensional pipeline layout.
(a)
(b)Based on case series 1, pipe cooling measure is added in case series 3 with the inlet water temperature of 15.0°C. In order to analyze the influence of pipe cooling separately, the daily variation of air temperature will not be considered in case series 3. The detailed condition is shown in Table6.Table 6
Calculation condition in case series 3.
Case
Air temperatureTa (°C)
Pipe inlet water temperatureTw (°C)
Note
3-1
30.0
15.0
The other condition is same as case 1
3-2
25.0
3-3
20.0
3-4
15.0
3-5
10.0
3-6
5.0Figure17 shows the comparison between the calculation results of case series 3 by FEM and (3). Figure 18 shows the comparison between the calculation results of case series 3 by FEM and (4) in which α equals 1.10.Comparison of case series 3 results between FEM and (3).
(a)
(b)Comparison of case series 3 results between FEM and (4).
(a)
(b)


## 3.3. Considering Both Daily Variation of Air Temperature and Pipe Cooling

In order to meet the requirement of practical situation, daily variation of air temperature and pipe cooling should be considered in the test block at the same time. In case series 4, pipe cooling measure is added on the basis of case series 2 calculation condition. The detailed condition can be seen in Table7.Table 7
Calculation condition of case series 4.
Case
Air temperatureTa (°C)
Pipe inlet water temperatureTw (°C)
Note
4-1
25
.
0
±
5
.
0
15.0
The other calculation condition is same as case series 1.
4-2
20
.
0
±
5
.
0
4-3
15
.
0
±
5
.
0
4-4
10
.
0
±
5
.
0Figure19 shows the relation between the temperature spatial gradient and the surface safety factor at different ages in case series 4. To observe it clearly, the results of the four times in every day of case 4-1 are drawn separately in Figure 20. It can be found that the temperature spatial gradient of 0.25 d is maximum and the safety factor of the time is minimum. Therefore, the results of 0.25 d in every case are drawn in Figure 21; the calculated value of (3) is shown as well. From this figure, the results by FEM and (3) do not meet very well, especially in the later age. While in Figure 22, the results become better when (4) is applied and α equals 1.10.Figure 19
Relation between concrete surface safety factor and temperature spatial gradient in case series 4.Figure 20
Relation between concrete surface safety factor and temperature spatial gradient at different times of every day in case 4-1.Figure 21
Comparison of case series 4 results at 0.25 d every day between FEM and (3).Figure 22
Comparison of case series 4 results at 0.25 d every day between FEM and (4).


## 3.4. Considering Both Daily Variation of Air Temperature, Pipe Cooling, and Surface Insulation

In winter or cold wave time, surface insulation measure is always taken to protect the concrete in early age. In case series 5, a layer of 2.0 cm thick expandable polyethylene (EPE), which is one of the most common application concrete insulation materials in construction, is added on the basis of case series 4 calculation condition. The detailed calculation condition is shown in Table8.Table 8
Calculation condition of case series 5.
Case
Air temperatureTa (°C)
Pipe inlet water temperatureTw (°C)
Surface insulation
Note
5-1
25
.
0
±
5
.
0
15.0
2.0 cm EPE material on all free face
The other condition is same as case series 1
5-2
20
.
0
±
5
.
0
5-3
15
.
0
±
5
.
0
5-4
10
.
0
±
5
.
0Figure23 shows the relation between the temperature spatial gradient and the surface safety factor at different ages in case series 5. To observe it clearly, the results of the four times in every day of case 5-1 are drawn separately in Figure 24. Being similar to case 4-1, the most dangerous situation appears at 0.25 d of every day. Therefore, the results of 0.25 d in every case are drawn in Figure 25; the calculated value of (3) is shown as well. From the figure, the results by FEM and (3) have obvious difference. When (4) is applied and α equals 1.15, the results become much closer, as shown in Figure 26.Figure 23
Relation between concrete surface safety factor and temperature spatial gradient in case series 5.Figure 24
Relation between concrete surface safety factor and temperature spatial gradient at different times of every day in case 5-1.Figure 25
Comparison of case series 4 results at 0.25 d every day between FEM and (3).Figure 26
Comparison of case series 5 results at 0.25 d every day between FEM and (4).


## 4. The Applicability of the Model for Different Concrete Materials

### 4.1. The Limitation of Basic Model and Calibration

The developing laws of concrete elastic module and tensile strength are assumed same as in Section3. In most concrete engineering, their developing laws are different, which will influence the concrete surface safety factor obviously. In this section, (4) is calibrated further to adapt the authentic different concrete.A calibration factorηγ is introduced into (4) to consider the influence of tensile strength and elastic module.η is the ratio of axial tensile strength and elastic module developing law, as shown in
(5)η=(1-e-a1τb1)(1-e-a2τb2),
in which (1-e-a1τb1) is the developing law of axial tensile strength, and (1-e-a2τb2) is the developing law of elastic module. a1,b1,a2,b2 are the parameters of the composite index formula. As for different age strength, the rebound number method [20] can be easily applied to obtain the values which will decide the parameters. For some prediction cases, concrete creep effect must be included. As in the finite element method simulation, the creep effect can always be considered in the expression of elastic module with an equivalent formula whose form is same but parameters’ values are different. In calculation, the essential function of the creep effect is to lower the elastic module to a certain extent.γ is a constant for a certain concrete material. Its value usually is between 0.9 and 1.2, which will vary with concrete material. The detailed information will be shown in later content of this paper.After the above modification, the proposed model can be expressed by(6)K(τ)=αf(∂T∂l)ηγg(τ),
in which f(∂T/∂l), g(τ), and η can be obtained separately by (1), (2), and (5). The definition of the parameters is same as the above equations, which will not be repeated again.


### 4.2. Verification of the Calibrated Model

Based on two mass concrete engineering in West China [21, 22], 7 types of concrete material are selected to verify the applicability of (6). Table 9 shows the detailed parameters’ values of different concrete materials, in which materials 1~3 are the normal concrete; materials 4~7 are the rolled compacted concrete (RCC). Case series 6 is planned in Table 10.Table 9
Main materials’ parameters of different concrete materials.
Material number
Adiabatic temperature (°C)θ(τ)=T0×(1-e-aτb)
Elastic module (GPa)E(τ)=E0×(1-e-a1τb1)
Tensile strength (MPa)ft(τ)=f0×(1-e-a2τb2)
θ
0
a
b
E
0
a
1
b
1
f
0
a
2
b
2
1
23.0
0.202
0.900
33.72
0.652
0.219
2.96
0.078
0.749
2
28.0
0.209
0.936
47.64
0.545
0.270
4.20
0.100
0.700
3
27.5
0.385
0.613
44.64
0.287
0.408
2.73
0.252
0.500
4
17.5
0.198
0.804
38.88
0.431
0.307
2.39
0.040
0.880
5
21.4
0.302
0.689
43.80
0.428
0.315
2.57
0.060
0.825
6
22.9
0.221
0.665
47.04
0.269
0.436
2.76
0.087
0.737
7
20.4
0.282
0.621
45.96
0.245
0.450
2.62
0.064
0.805Table 10
The calculation condition of case series 6.
Case
Air temperatureTa (°C)
Initial temperatureT0 (°C)
Note
6-1
20.0
25.0
No pipe cooling and  no surface insulation measure
6-2
20.0
6-3
15.0


### 4.3. Test Results and Analysis

The 7 kinds of material in Table9 are separately applied in the 3 different calculation conditions in Table 10. Therefore there are indeed 21 calculation cases in total. Figures 27, 28, 29, 30, 31, 32, and 33 are the comparison of calculation result by FEM and (6) with different materials and different calculation conditions.Figure 27
Comparison of case series 6 results using material 1 between FEM and (6).Figure 28
Comparison of case series 6 results using material 2 between FEM and (6).Figure 29
Comparison of case series 6 results using material 3 between FEM and (6).Figure 30
Comparison of case series 6 results using material 4 between FEM and (6).Figure 31
Comparison of case series 6 results using material 5 between FEM and (6).Figure 32
Comparison of case series 6 results using material 6 between FEM and (6).Figure 33
Comparison of case series 6 results using material 7 between FEM and (6).According to the comparison, it is found that the calculation result by FEM and (6) meet well, especially in early age.Table11 lists the values of parameter α and γ in different conditions. It can be seen that the value of α ranges from 0.7 to 1.1 and the value γ ranges from 0.9 to 1.2. For normal concrete materials 1~3, the parameter value varies little, while for the RCC materials 4~7, the parameter value changes more.Table 11
The relation between model parameter value and concrete material.
Material number
1
2
3
4
5
6
7
Parameterα
0.8
0.85
0.85
0.75
0.7
0.95
1.1
Parameterγ
0.9
0.9
1.0
1.0
1.0
1.0
1.2In the engineering application, the two parameters can be decided from Table11 by comparing the concrete material properties with the value in Table 9. And the concrete temperature spatial gradient ∂T/∂l can be obtained by two sensors in 0.05 m deep and 0.50 m deep. Then the surface safety factor to prevent crack at any age can be calculated by (6) easily. The precision of the equation calculation result is high as FEM, especially in early age when the surface cracking appears most often. During the dynamic observation and calculation in construction, if the safety factor tends to be lower than the expected value, necessary temperature control method should be taken to reduce the temperature difference between concrete inside and outside, such as thickening cover insulation material or enhancing the pipe cooling.


## 4.1. The Limitation of Basic Model and Calibration

The developing laws of concrete elastic module and tensile strength are assumed same as in Section3. In most concrete engineering, their developing laws are different, which will influence the concrete surface safety factor obviously. In this section, (4) is calibrated further to adapt the authentic different concrete.A calibration factorηγ is introduced into (4) to consider the influence of tensile strength and elastic module.η is the ratio of axial tensile strength and elastic module developing law, as shown in
(5)η=(1-e-a1τb1)(1-e-a2τb2),
in which (1-e-a1τb1) is the developing law of axial tensile strength, and (1-e-a2τb2) is the developing law of elastic module. a1,b1,a2,b2 are the parameters of the composite index formula. As for different age strength, the rebound number method [20] can be easily applied to obtain the values which will decide the parameters. For some prediction cases, concrete creep effect must be included. As in the finite element method simulation, the creep effect can always be considered in the expression of elastic module with an equivalent formula whose form is same but parameters’ values are different. In calculation, the essential function of the creep effect is to lower the elastic module to a certain extent.γ is a constant for a certain concrete material. Its value usually is between 0.9 and 1.2, which will vary with concrete material. The detailed information will be shown in later content of this paper.After the above modification, the proposed model can be expressed by(6)K(τ)=αf(∂T∂l)ηγg(τ),
in which f(∂T/∂l), g(τ), and η can be obtained separately by (1), (2), and (5). The definition of the parameters is same as the above equations, which will not be repeated again.


## 4.2. Verification of the Calibrated Model

Based on two mass concrete engineering in West China [21, 22], 7 types of concrete material are selected to verify the applicability of (6). Table 9 shows the detailed parameters’ values of different concrete materials, in which materials 1~3 are the normal concrete; materials 4~7 are the rolled compacted concrete (RCC). Case series 6 is planned in Table 10.Table 9
Main materials’ parameters of different concrete materials.
Material number
Adiabatic temperature (°C)θ(τ)=T0×(1-e-aτb)
Elastic module (GPa)E(τ)=E0×(1-e-a1τb1)
Tensile strength (MPa)ft(τ)=f0×(1-e-a2τb2)
θ
0
a
b
E
0
a
1
b
1
f
0
a
2
b
2
1
23.0
0.202
0.900
33.72
0.652
0.219
2.96
0.078
0.749
2
28.0
0.209
0.936
47.64
0.545
0.270
4.20
0.100
0.700
3
27.5
0.385
0.613
44.64
0.287
0.408
2.73
0.252
0.500
4
17.5
0.198
0.804
38.88
0.431
0.307
2.39
0.040
0.880
5
21.4
0.302
0.689
43.80
0.428
0.315
2.57
0.060
0.825
6
22.9
0.221
0.665
47.04
0.269
0.436
2.76
0.087
0.737
7
20.4
0.282
0.621
45.96
0.245
0.450
2.62
0.064
0.805Table 10
The calculation condition of case series 6.
Case
Air temperatureTa (°C)
Initial temperatureT0 (°C)
Note
6-1
20.0
25.0
No pipe cooling and  no surface insulation measure
6-2
20.0
6-3
15.0


## 4.3. Test Results and Analysis

The 7 kinds of material in Table9 are separately applied in the 3 different calculation conditions in Table 10. Therefore there are indeed 21 calculation cases in total. Figures 27, 28, 29, 30, 31, 32, and 33 are the comparison of calculation result by FEM and (6) with different materials and different calculation conditions.Figure 27
Comparison of case series 6 results using material 1 between FEM and (6).Figure 28
Comparison of case series 6 results using material 2 between FEM and (6).Figure 29
Comparison of case series 6 results using material 3 between FEM and (6).Figure 30
Comparison of case series 6 results using material 4 between FEM and (6).Figure 31
Comparison of case series 6 results using material 5 between FEM and (6).Figure 32
Comparison of case series 6 results using material 6 between FEM and (6).Figure 33
Comparison of case series 6 results using material 7 between FEM and (6).According to the comparison, it is found that the calculation result by FEM and (6) meet well, especially in early age.Table11 lists the values of parameter α and γ in different conditions. It can be seen that the value of α ranges from 0.7 to 1.1 and the value γ ranges from 0.9 to 1.2. For normal concrete materials 1~3, the parameter value varies little, while for the RCC materials 4~7, the parameter value changes more.Table 11
The relation between model parameter value and concrete material.
Material number
1
2
3
4
5
6
7
Parameterα
0.8
0.85
0.85
0.75
0.7
0.95
1.1
Parameterγ
0.9
0.9
1.0
1.0
1.0
1.0
1.2In the engineering application, the two parameters can be decided from Table11 by comparing the concrete material properties with the value in Table 9. And the concrete temperature spatial gradient ∂T/∂l can be obtained by two sensors in 0.05 m deep and 0.50 m deep. Then the surface safety factor to prevent crack at any age can be calculated by (6) easily. The precision of the equation calculation result is high as FEM, especially in early age when the surface cracking appears most often. During the dynamic observation and calculation in construction, if the safety factor tends to be lower than the expected value, necessary temperature control method should be taken to reduce the temperature difference between concrete inside and outside, such as thickening cover insulation material or enhancing the pipe cooling.


## 5. Conclusion

The reasonable definition and criteria of the temperature difference between concrete inside and outside are researched in this paper in order to prevent the surface cracking in early age. The surface concrete temperature spatial gradient in 0.5 m deep is found as a more reasonable and feasible index. Then a new model, which can be used to calculate the surface safety factor simply, is deduced. If the necessary parameters are known, the surface safety factor of a certain age can be obtained according to the concrete temperature spatial gradient of the age. During concrete construction the concrete temperature spatial gradient can be recorded easily by two or more temperature sensors. With the help of the proposed model, the surface safety factor can be adjusted by observing and controlling the temperature spatial gradient dynamically. The precision of the model under different conditions is tested and improved by comparison with FEM. The last proposed model (6) can meet the practical engineering requirement, which will bring convenience to designer and constructor. In the future, it will be further verified in more engineering applications.



---

*Source: 183209-2014-03-18.xml*