Early Cracking Risk Prediction Model of Concrete under the Action of Multifield Coupling

(rough the adiabatic temperature rise experiment, the adiabatic temperature rise of concrete with hydration time was recorded. Based on the maturity degree theory, the relationship between the hydration degree of the concrete and the equivalent age was determined. (en, the hydration degree prediction model of the concrete’s early elastic modulus and tensile strength was established. (e local temperature and humidity of the concrete were measured by the shrinkage experiment, and based on the capillary water tension theory, a temperature-humidity prediction model for the early shrinkage of the concrete was designed. According to the ratio of the creep deformation and elastic deformation of concrete which were obtained through the restraint ring experiment, a model for predicting the early creep coefficient of concrete was proposed. Based on the coupling effect of “hydration-temperature-humidity,” a predictionmodel of early cracking risk coefficient of concrete under multifield coupling was proposed. Finally, several groups of slab cracking frame experiments were carried out, and the cracking risk prediction results of concrete were consistent with the actual situation, which indicated the correctness of the early cracking risk prediction model of concrete.


Introduction
e large use of superplasticizers and mineral admixtures, the higher fineness of cement, and the lower water-binder ratio are the typical characteristics of modern concrete. erefore, the strength of concrete can be significantly improved. However, there are some problems of modern concrete, such as the rapid hydration, the large shrinkage, the rapid internal temperature rises, and the rapid internal humidity decrease, which led to early cracking of concrete. If the early cracks of concrete cannot be strictly controlled, especially the cracks caused by nonload factors in the early stage of concrete structure, the durability and safety of concrete would be reduced [1][2][3][4][5]. Meanwhile, fiber [6] would be used to improve the anticracking property of concrete.
A large number of studies have shown that [7][8][9][10][11][12] the anticracking property of concrete would be directly or indirectly affected by the change of internal factors such as hydration, temperature, humidity, shrinkage, and creep. Shen et al.'s [10] research found that the curing condition of concrete at 20°C had better anticracking properties than that of 45°C. Zhang et al. [8] established a relationship between the internal humidity of concrete and the free shrinkage strain of concrete at the macro level. e results show that the faster the internal humidity of concrete decreased, the greater the shrinkage at the macro level. Bentz [9] defined the ratio of creep strain to effective shrinkage strain as the creep coefficient of concrete, and the shrinkage of concrete was one of the main factors affecting the cracking, so the early anticracking property of concrete would also be affected by the early creep of concrete. e above research showed that the early anticracking property of concrete had a complex relationship with cement hydration, temperature, humidity, shrinkage, and creep. Many scholars have studied the early anticracking property of concrete with mineral admixtures such as fly ash and slag but have not yet reached a consistent conclusion [2,[13][14][15]. erefore, it is necessary to propose an accurate method to predict the early cracking risk of modern concrete.
Liu et al. [16] used the dumbbell method to evaluate the early anticracking performance of concrete materials under constraint conditions. Although this method can quickly evaluate the anticracking property of different concrete materials, it was only aimed at the rapid qualitative detection and analysis of the construction site and cannot quantitatively evaluate the cracking of concrete. Meagher [6] used HIPERPAV software to simulate the influence of admixtures on the early cracking risk of concrete repair plate; Zhang [6] used the ring restraint experiment to evaluate the anticracking property of low shrinkage cement-based materials. Although the above method can evaluate the early anticracking property of concrete, the root cause of the early cracking of concrete cannot be obtained.
Based on (1), Bendimerad et al., [17] neglected the early creep of concrete and considered that the deformation tensile stress σ t (t) was only a function of the elastic modulus E (t) and the plastic shrinkage strain ε sh (t), that is, σ t (t) � E (t) × ε sh (t), and the plastic shrinkage and cracking risk of recycled aggregate concrete within 24 hours was explored. Which method can predict the risk factor η of concrete cracking within 24 hours and the mathematical model of early cracking prediction of concrete has also been defined. However, as the age of concrete increased, the creep of concrete would become larger [18,19], and it also led to the limitation of Ahmed's method to predict the early cracking of concrete due to neglecting creep, which was only applicable to the evaluation of concrete cracking resistance within 24 hours.
erefore, based on the adiabatic temperature rise experiment, shrinkage experiment, and creep experiment of concrete, a prediction model that described the cracking of concrete due to early multifield coupling effects such as "hydration-temperature-humidity" was proposed in this study. e basic cause of concrete early cracking from the mechanism was explained, and the early cracking prediction model of concrete proposed by Bendimerad et al., [17] has been further improved by this prediction model. At the same time, the slab cracking frame experiment method was used to test the cracking resistance of several groups of concrete. Finally, the test results of Khan et al. [20] and Zhang et al. [6] were selected, and the cracking risk prediction model proposed in this study was used to verify the test results of the literature [6,20].

Materials and Experiment Design
e raw materials used for concrete were cement, fly ash, slag, sand, coarse aggregate, superplasticizer, and tap water. e physical and chemical properties of these raw materials were as follows: P·O42.5 ordinary Portland cement conformed with the current Chinese national standard GB 175-2007 "Common Portland Cement." e 3 d and 28 d compressive strength of P·O42.5 ordinary Portland cement was 27 MPa and 49 MPa, respectively, and the specific gravity and specific surface area were 3.16 g/cm 3 and 354 m 2 /kg respectively. Grade I fly ash conformed with the current Chinese national standard GB/T 1596-2017 "Fly ash used for cement and concrete." e ignition loss, water demand ratio, and moisture content of fly ash were 3.7%, 88%, and 0.4%, respectively; e specific surface area, 7-day activity index, and ignition loss of slag were 425 m 2 /kg, 78%, and 0.26%, respectively. And the main chemical composition is shown in Table 1.
e fineness modulus of the river sand was 2.56. e particle size of calcareous coarse aggregate was a continuous gradation of 5-20 mm. e polycarboxylate superplasticizer had a water reduction rate of 31%.

Mix Proportion Design of Concrete.
According to the concrete mix proportion shown in Table 2, the elastic modulus experiment, axial tensile experiment, adiabatic temperature rise experiment, shrinkage experiment, and creep experiment of concrete were carried out in this study.
e above experiments were mainly used to determine the parameters in the prediction model of concrete cracking risk coefficient. Finally, in order to verify the applicability of the early cracking prediction model of concrete, 3 groups of slab cracking frame experiments were carried out according to the concrete mix proportion shown in Table 2

Adiabatic Temperature Rise Experiment of Concrete.
According to the current Chinese national standard for "Test code for hydraulic concrete" (SL 352-2006), based on the 3 types of concrete mix proportions shown in Table 2, 50 L concrete samples were selected for the adiabatic temperature rise experiment. At the early stage of the experiment, the adiabatic temperature rise was recorded every 0.5 h, and every 5 h after 24 h, until the end of the adiabatic temperature rise experiment, which lasted a total of 7 days. e 2 Advances in Materials Science and Engineering adiabatic temperature rise experiment device is shown in Figure 1(a). e temperature sensor in the figure was used to record the adiabatic temperature rise value of the concrete. e temperature compensation circuit and the heat insulation layer were mainly used to avoid the loss of hydration heat of concrete.

Shrinkage Experiment of Concrete.
e coupling effect of the early "temperature-humidity" field of concrete was considered in this study, according to the current Chinese national standard for "Test code for hydraulic concrete" (SL 352-2006), based on the 3 types of concrete mix proportions shown in Table 2, the concrete shrinkage deformer as shown in Figure 1(b) [8] was used, the size of the concrete specimen was 100 mm × 100 mm × 515 mm, and a temperature and humidity senor was installed in the center of the specimen, which can automatically collect the deformation, temperature, and humidity at the inner center of the concrete. e room temperature was about 22°C, and the relative humidity was about 60%. e shrinkage of concrete in 1 d, 3 d, 5 d, 7 d, 14 d, 21 d, and 28 d was recorded.

Creep Experiment of Concrete.
Gawin et al. [21] proposed a mathematical calculation method for creep of concrete, which decomposed the total creep strain rate tensor into two parts: viscoelastic strain rate and viscous flow strain rate. Li et al. [11] analyzed the early creep performance of concrete by using the ring restraint experiment, which showed that this method can not only ensure the accuracy of the early creep of concrete measured by the experiment but also be simpler than the calculation method of concrete early creep proposed by Gawin et al. [21]. erefore, the ring restraint experiment was also adopted in this study, and the early creep experiment of concrete was carried out according to the concrete mix proportion shown in Table 2.

Early Cracking Experiment of Concrete.
According to the current Chinese national standard for "Standard for test method for long-term performance and durability of ordinary concrete" (GB/T 50082-2009) in this study, a slab cracking frame as shown in Figure 1(c) was selected. e size of the cracking frame was 900 mm × 700 mm × 100 mm. e side plates, bolts, and stiffening rib were mainly used to fix the steel formwork of the cracking frame, and the crack inducer was mainly used to induce the early cracking of concrete due to shrinkage. is experiment was mainly to observe the cracking conditions after the initial setting of concrete and was used to verify the correctness of the cracking risk prediction model proposed in this paper. e room temperature was about 22°C, and the relative humidity was about 60%. When the concrete cracked, the test observation stopped.

e Establishment of Early Cracking Risk Prediction Model of Concrete.
According to the criterion of concrete cracking [17], the cracking risk coefficient η was used to evaluate the early cracking resistance of modern concrete, as shown in the following equation: where σ t (t) is the deformation tensile stress (cracking driving force) of concrete at the early age t, MPa; f t (t) is the tensile strength of concrete at the early age t, MPa; E (t) is the elastic modulus of the concrete at the early age t, MPa; ε sh−e (t) is the effective shrinkage strain of the concrete at the early age t. Among them, when η < 0.7, the concrete does not crack; when 0.7 ≤ η ≤ 1.0, the concrete may crack; when η > 1.0, the concrete cracks. Gao et al.'s [22,23] research showed that the effective shrinkage strain ε sh−e (t) was the result of the interaction of the shrinkage strain and the creep strain (strain of concrete due to creep) of the concrete, as shown in where t and t 0 are the age of concrete and the age of concrete under loading, respectively, d; ε sh (t) is the shrinkage strain of the concrete at the age t; ε creep (t, t 0 ) is the shrinkage strain of the concrete due to creep in the early age (t, t 0 ). Note. TS-1, TS-2, and TS-3 represented 3 experiment mix proportions of concrete, respectively.
where ϕ (t,t 0 ) is the creep coefficient of concrete at an early age (t, t 0 ). It can be known from equations (3) and (4) that the expression of the relationship between the effective shrinkage strain ε sh−e of concrete at an early age (t, t 0 ) is shown as follows:   Advances in Materials Science and Engineering Equation (5) was substituted into equation (2), and the risk prediction model of early cracking of modern concrete was obtained as follows: erefore, based on the prediction model of concrete cracking proposed by Ahmed [17], the effect of creep on the early cracking of concrete was further considered, and the prediction model of modern concrete cracking risk with creep coefficient was put forward innovatively. In Section 3.2, based on the coupling effect of multifields such as "hydration-temperature-humidity-force" in the early stage of concrete, the elastic modulus E (t) , tensile strength f t (t), shrinkage strain ε sh (t), and creep coefficient ϕ (t,t 0 ) of concrete at age t would be determined.

Experiment Results of Adiabatic Temperature Rise and Determination of Cement Hydration
Degree. e adiabatic temperature rise experiment results of TS-1, TS-2, and TS-3 concrete are shown in Figure 2. Based on the adiabatic temperature rise experiment of concrete, the hydration degree α of cement can be expressed as [7] α � where ΔT (t) is the adiabatic temperature rise value of the concrete temperature rise sample at the age t,°C; ΔT ∞ is the adiabatic temperature rise value reached when the concrete is fully hydrated,°C.
However, the cement in the concrete cannot be completely hydrated [7,24]. It was assumed that when the adiabatic temperature rise of the concrete reached ΔT max , the cement hydration was completed, and the final hydration degree α u of the cement can be expressed as An important typical characteristic of modern concrete was used of mineral admixtures. Schindler and Folliard [25] had studied a lot of research on the final hydration degree of cement in modern concrete mixed with slag and fly ash, and the final hydration degree model of cement containing slag and fly ash was proposed as follows: where W/C represents the mass ratio of water and cement in concrete; P F represents the mass percentage of fly ash in cement in total cementitious materials; P S represents the percentage of slag in cement in total cementitious materials. Luzio and Cusatis [26] concluded that the final hydration degree of cement in concrete mixed with silica fume was expressed as follows: where α ∞ s is the final reaction degree of silica fume, which was expressed as follows [27]: where SF eff is the "effective coefficient" of silica fume, that is, the mass fraction of SiO2 in silica fume; (S/C) req � min(0.16, 0.4(W/C)). According to equations (7) and (8), the relationship between the hydration degree α(t) of cement at age t and the final hydration degree α u of the cement can be expressed as where ΔT (t) and ΔT max can be obtained from the adiabatic temperature rise curve. eoretically, the longer the adiabatic temperature rise experiment duration, the more accurate the measured value of ΔT max was. However, under different temperature histories of the same mix ratio, the internal cement hydration degree was different at the same age. erefore, based on the maturity theory, the concept of equivalent age was introduced to eliminate the influence of the temperature history of concrete on its concrete internal hydration. e age t required for cement hydration in concrete under different temperature histories was equivalent to the age t eq required for concrete under reference temperature (usually 20°C), and its calculation expression was as follows [7]: Advances in Materials Science and Engineering where t eq is the equivalent age, d; R is the ideal gas constant, 8.314 J/(mol · K); U ar is the activation energy of the internal cement hydration reaction at the reference temperature of the concrete; T is the temperature of the concrete,°C; U aT is the activation energy of the internal cement hydration reaction at temperature T of the concrete. Kim [28] found that the apparent activation energy can be regarded as a constant in the early stage of cement hydration. However, with the hydration reaction proceeds of the cement, the effects of temperature T and curing time on the apparent activation energy of the cement need to be considered. At this time, U aT can be expressed as where t is the curing age, d. e continuous curve of (13) was divided into n segments, the interval of each segment is (t i − t i−1 ), and the expression as shown in (15) was obtained: Pane and Hansen [29] explored the relationship between the hydration degree α of cement and the equivalent age t eq , a model of the equivalent age and adiabatic temperature rise experiment of the cement hydration degree was established, and the relationship expression as shown in equation (16) was proposed: where the constants of A and B were obtained by fitting the model according to the above equation. erefore, the research ideas were as follows in this section, the fresh concrete was placed in the adiabatic temperature rise experiment chamber, and then the temperature inside the adiabatic temperature rise experiment chamber was recorded by the temperature sensor, and the relationship curve between adiabatic temperature rise values and hydration time was drawn. e final hydration degree α u of concrete was calculated according to equation (10), and then the hydration degree α(t) corresponding to each hydration time was determined according to equation (12). At the same time, based on the maturity theory, the concept of equivalent age was introduced to eliminate the influence of the temperature history of concrete on its internal cement hydration degree. According to equation (15), the equivalent age t eq corresponding to the actual hydration time was determined. Finally, based on (16), the parameters A and B were determined; that is, the mathematical model of concrete hydration degree and equivalent age was determined.
To sum up, according to the adiabatic temperature rise curve shown in Figure 2 and equations (7)- (16) in Section 3.1, the test value and fitting curve of cement hydration degree and equivalent age are obtained as shown in Figure 3, and the relevant parameters and model fitting parameters are listed in Table 3.
e R 2 values of the model result and the test result in Figure 3 were not less than 0.98, indicating that the correlation was very high. It can be seen that the model of cement hydration degree and equivalent age proposed by Pane [29] can be well described the development law of cement hydration degree and equivalent age in different concrete mix proportions, and the effect of temperature history on cement internal hydration degree of the concrete can be eliminated.

Determination of Early Elastic Modulus Prediction
Model of Concrete. Zhang and Schutter [7,30] have explored the relationship between the mechanical properties (elastic modulus, tensile strength, etc.) of concrete and the hydration degree of cement. erefore, according to the test results of the elastic modulus of concrete, the hydration degree method was used to describe the early elastic modulus of concrete. e prediction model of the hydration degree of the elastic modulus was as follows: where E(α) is the elastic modulus of concrete when the degree of cement hydration is α, MPa; β E is the subsequent increase coefficient of the elastic modulus of concrete after 28 days, 1.05; E 28 is the elastic modulus of concrete after 28 days of curing, GPa; α 0 is the hydration degree at the initial setting of the concrete; α u is the final hydration degree of the concrete; b is a constant obtained by fitting the hydration degree model of the concrete's elastic modulus to the test values. e measured value and fitting curve of elastic modulus of concrete are shown in Figure 4, and the relevant parameters and model fitting parameters are listed in Table 3.
In Figure 4, R2 of TS-1, TS-2, and TS-3 model results and test values were 0.94, 0.95, and 0.98, respectively, with high correlation. It can be seen that the early elastic modulus prediction model of concrete proposed by Zhang and Schutter [7,30] can be well described the relationship between the elastic modulus and hydration degree of concrete.

Determination of Early Tensile Strength Prediction
Model of Concrete. According to the test results of the tensile strength of concrete, the hydration degree of concrete was used to describe the early tensile strength.
e prediction model of the tensile strength was [17,31] where f(α) is the tensile strength of concrete when the degree of cement hydration is α, MPa; β f is the subsequent increase coefficient of the tensile strength of concrete after 28 days, 1.05; f 28 is the tensile strength of concrete after 28 days of curing, GPa; the meaning of α 0 and α u was shown in equation (17); d is a constant obtained by fitting the hydration degree model of the concrete's tensile strength to the test values. e measured tensile strength and fitting curve of concrete are shown in Figure 5, and the relevant parameters and model fitting parameters are listed in Table 3.
e R2 values of model results and test values in Figure 5 were not less than 0.95, and the correlation between the predicted results and test results was very high. It can be seen that the early tensile strength prediction model of concrete proposed by Zhang and Schutter [7,30] can be well described the relationship between concrete tensile strength and hydration degree.

Early Shrinkage Prediction Model of Concrete.
Zhang [8] considered that the change of the internal humidity of concrete was the root cause of the early autogenous shrinkage and dry shrinkage of concrete, and the relationship between early shrinkage of concrete and internal humidity was explored by using a normal experiment environment without temperature control. e results showed that although the internal humidity of the concrete changed gradually along the thickness direction of the specimen, the local humidity can be used to estimate the shrinkage strain of the concrete specimen at the macro level, and then the shrinkage stress of the structure can be calculated. Bentz et al., [32] used capillary water tension theory to predict the shrinkage deformation of concrete. e results showed that the concrete shrinkage value can be calculated by knowing the internal relative humidity of the concrete. erefore, based on the coupling effect of the "temperaturehumidity" field in the early stage of the concrete, the capillary water tension theory was used to establish the early shrinkage strain prediction model of the concrete according to the local temperature and humidity of the concrete.

(1) Shrinkage Prediction Model of Concrete Based on Humidity
Field. e change of internal relative humidity of concrete would lead to the change of capillary water tension. e relationship between the internal humidity of concrete and the age measured by the shrinkage experiment in this study is   shown in Figure 6. e mechanical model of the concrete element under the tension of capillary water is shown in Figure 7. e calculation process of concrete element shrinkage caused by capillary water tension was as follows: where d x , d y , and d z are the side lengths of the concrete element.
When the volume of a concrete element changes due to internal humidity, the shrinkage strains of the concrete element along the X-axis, Y-axis, Z-axis directions were assumed as ε x , ε y , and ε z , respectively, and the volume V of the concrete element at this time was erefore, the volume strain ε V (change of unit volume) of the concrete was After the concrete element was loaded, the shrinkage strain of the concrete element along the X-axis, Y-axis, and Z-axis directions was where E is the static elastic modulus of the concrete, GPa; μ is Poisson's ratio of the concrete, taking 0.2 [33]; σ x , σ y , and σ z represent the capillary water tension of the concrete element in the X-axis, Y-axis, and Z-axis directions. e higher-order terms of equation (21) were ignored, equations (22)-(24) were substituted into (21), and the volume strain of concrete was obtained as follows: where the bulk modulus K � E/3(1 − 2μ), GPa.
Assuming that the strains in each direction of concrete shrinkage were equal. Under the action of capillary water tension σ, the shrinkage strain ε w in the X-axis, Y-axis, and Z-axis directions was Lura et al. [33] considered that the capillary pores caused by chemical shrinkage of the cementitious material internal the concrete were completely wetted, and the contact angle between water and solids was 0. At this time, the capillary water tension of the concrete was σ � − ρR(273 + T)ln RH M , (27) where σ is the capillary water tension of concrete; ρ is the density of water, kg/m 3 ; R is the ideal gas constant, 8.314 J/ (molK); T is the internal temperature of the concrete as shown in Figure 6,°C; RH is the internal relative humidity of the concrete as shown in Figure 6; M is the molar mass of water, g/mol. However, concrete was a heterogeneous material and contains pores in it, so the bulk modulus K of the concrete in equations (25)- (26) needed to be modified [34]. Meanwhile, because the water in the pores was consumed at the hydration of cement in the early stage, making the pores tend to be unsaturated [31], it was necessary to introduce saturation S to further modify the shrinkage strain of concrete in (26). erefore, the equation after modify of concrete was (shrinkage prediction model of concrete based on humidity field) where σ is the capillary water tension of the concrete, which was calculated according to equation (27); the specific calculation process of the modified bulk modulus K and the saturation S was shown in equations (29)- (30). e reciprocal of the modified bulk modulus K [4, 34] was where K is the bulk modulus used for calculating the shrinkage strain of concrete in equation (28); the bulk modulus , where E C is the elastic modulus of concrete (calculated according to (17)); µ is Poisson's ratio of concrete, taking 0.2 [33], the units of K C and E C are GPa; K S is the bulk modulus of the solid material, taking 44 GPa [33]. e saturation S [35] was where p � (W/C)/(W/C) + (ρ W /ρ C ), W/C is the mass ratio of water and cement, ρ W /ρ C is the density ratio of water and cement, and α is the hydration degree of concrete.
(2) Shrinkage Prediction Model of Concrete Based on Temperature Field. Based on the shrinkage experiment of concrete, the relationship between the internal temperature of concrete and its age was measured in this study as shown in Figure 6. e finite element method was used by Liu et al. [36,37] to predict the anticracking property of concrete. e thermal expansion coefficient of concrete was regarded as a constant, and the values were 9.94 με/°C and 10 με/°C. Li [11] determined that the thermal expansion coefficient of C30 and C80 8 Advances in Materials Science and Engineering concrete was between 10 με/°C and 12 με/°C through the variable temperature experiment of concrete. erefore, the thermal expansion coefficient was taken as 10 με/°C in this study, and it was applied to the prediction model of concrete shrinkage. At this time, the temperature deformation ε T of concrete with temperature history between (T 0 , T) was (shrinkage prediction model of concrete based on temperature field) where β T is the thermal expansion coefficient of concrete, με/°C; T 0 is the internal temperature of the concrete at the initial setting,°C; T is the internal temperature of the Advances in Materials Science and Engineering 9 concrete at time t,°C; and ΔT is the difference between T and T 0 ,°C.

(3) Determination and Verification of Early Shrinkage Prediction Model of Concrete.
In this study, the shrinkage prediction model of concrete was determined based on the humidity field (based on the capillary water tension theory) and the temperature field in Section 3.2.4 (2), and the shrinkage prediction model of concrete in the early stage under the coupling effect of "temperature-humidity" field was obtained as follows: erefore, based on equation (32), according to the relationship curve of temperature-humidity and age shown in Figure 6, combined with the hydration degree prediction model of concrete elastic modulus, the prediction value of concrete shrinkage strain shown in Figure 8 was obtained, which showed that the test results almost coincided with the prediction results, indicating that the method of using the local temperature and humidity inside the concrete to predict the macro shrinkage of the component was feasible and also verified the accuracy of the early shrinkage prediction model of concrete, which was consistent with the conclusion of Zhang [8].

Determination of Early Creep Prediction Model of
Concrete. In this study, based on the ring restraint test results in Section 2.2.4, the creep coefficient of concrete was defined as the ratio of creep deformation to elastic deformation under circumferential constraint [9]. e early creep coefficient of concrete obtained by the test is shown in Figure 9. According to the early creep results of concrete explored by Jiang et al. [38], the concrete creep coefficient was fitted with an exponential model, and the specific fitting equation was as follows: where φ t represents the creep coefficient of concrete at time t; D 1 , D 2 , and φ 0 are parameters obtained by fitting the test values based on the exponential model (equation (33)), and the fitting results are shown in Table 4. e R 2 values of the predicted results and test values in Figure 9 were not less than 95%, and the correlation between the predicted results and test results was very high. It can be seen that the early creep prediction model of concrete proposed by Schutter [7,30] can well describe the relationship between concrete early creep and age.

Determination of Cracking Risk Prediction Model of Concrete under the Action of Multifield Coupling.
According to the coupling effect of multiple-fields such as "hydration-temperature-humidity-force" in the early stage of concrete, and based on equation (6), from the perspective of the change mechanism of the internal hydration degree, temperature, and humidity of the concrete, the cracking risk prediction model of modern concrete was proposed as follows:

Verification of the Cracking Risk Prediction Model of Modern Concrete
In this study, a slab cracking frame was used for the early cracking experiment of concrete, and the cracking of TS-1, TS-2, and TS-3 concrete after the initial setting was recorded. Based on equation (6), according to the elastic modulus, tensile strength, macro shrinkage strain, and creep coefficient of concrete, the experimental value of the concrete cracking risk coefficient was calculated. Meanwhile, based on equation (34), according to the hydration degree, temperature, humidity, and creep of concrete in the early stage, the predicted value of the concrete cracking risk coefficient was calculated. Finally, according to the cracking risk coefficients calculated by equations (6) and (34), they were compared with the cracking state measured by the slab cracking frame experiment to verify the accuracy of the early cracking risk prediction model of concrete. In addition, the experiment results of Khan [20] and Zhang [6] were selected in this paper, and the cracking risk prediction model proposed in this study was used to verify the test results of the literature [6,20].

Cracking the State of Concrete in the Slab Cracking Frame
Experiment. e cracking state of the concrete in the slab cracking frame experiment is shown in Figure 10. e TS-1 mix proportion was used for the slab cracking frame experiment, and the concrete began to crack on the third day. e TS-2 mix proportion was used for the slab cracking frame experiment, and the concrete cracks appeared on the first day. e TS-3 mix proportion was used for the slab cracking frame experiment, and the concrete had cracks on the first day.

Verification of the Cracking Risk Prediction Model of
Concrete. Based on the test results of the elastic modulus, tensile strength, shrinkage strain, creep coefficient, and adiabatic temperature rise values of the concrete, the cracking risk coefficients calculated according to equations (6) and (34) are shown in Table 5. Khan et al. [20] believed that the early cracking of concrete was caused by creep and shrinkage. Finally, the correctness of his opinion was verified by the finite element method. Based on the restraint ring experiment, the concept of equivalent shrinkage strain (equivalent shrinkage strain was K × ε, where K is the creep coefficient of concrete) was proposed by Zhang et al. [6], and the anticracking property of ECC (traditional engineered cementitious composite) and LSECC (low shrinkage engineered cementitious composite) concrete was evaluated by mathematical calculation model. erefore, the test results of Khan et al. [20] and Zhang et al. [6] were selected in this study, the cracking risk prediction model proposed in this paper was adopted, and the test results of literature [6,20] were further verified. e calculation data and results of cracking risk are listed in Table 5.
It can be seen from Table 5 that the early cracking risk coefficient of concrete calculated based on equations (6)   and cracking state of concrete, which showed that the early cracking risk prediction model of the modern concrete can be well predicted the cracking resistance of concrete. Khan et al. [20] used the restraint ring experiment to test the creep of concrete and found that the concrete just cracked on the 15.5 days. According to the relevant data of literature [20] (as shown in Table 5), and based on the prediction model of cracking risk coefficient proposed in this paper, the cracking risk coefficient of concrete on the 15.5 days was calculated to be 1.0033. According to the cracking criteria of concrete, it can be seen that the concrete just cracks at this time, which was consistent with the results of the literature [20].
Based on the restraint ring experiment, Zhang et al. [6] proposed a mathematical model for predicting the cracking of ECC and LSECC concrete, which predicted that the ECC concrete would crack on the 5th day, while the cracking time of the experiment was on the 7th day; LSECC concrete mixed with fiber had not been cracked. According to the relevant data of literature [6] (as shown in Table 5), and based on the prediction model of cracking risk coefficient proposed in this paper, the cracking risk coefficient of ECC concrete on the 5th day was 0.798, and the cracking risk coefficient on the 7th day was 2.135; the cracking risk coefficient of LSECC concrete on the 28th day was 0.697. erefore, the prediction results showed that ECC concrete could crack on the 5th day, and it was certain that the cracking would occur on the 7th day, while the cracking of LSECC concrete had not been cracked within 28 d, and the prediction results were consistent with the literature [6], which further showed that the concrete cracking risk prediction model proposed in this paper was correct and reasonable.

Conclusion
According to the adiabatic temperature rise experiment, shrinkage experiment, creep experiment, and theoretical analysis of concrete, the early cracking risk prediction model of modern concrete was established. e early cracking risk coefficient of concrete calculated by the model was consistent with the actual cracking state based on the slab cracking frame experiment of concrete. e following conclusions were drawn: (1) Based on the capillary water tension theory, the shrinkage strain of the specimen at the macro level can be estimated by using the local temperature and humidity inside the concrete, and then the shrinkage stress of the structure can be calculated (2) e concrete cracking risk prediction model with creep coefficient can well predict the early cracking of the concrete (3) When the concrete mix proportion, adiabatic temperature rise curve, 28 d elastic modulus, 28 d tensile strength, internal temperature and humidity curve, and early creep coefficient are known, the cracking risk coefficient of concrete can be well predicted, and then the early cracking resistance of concrete can be evaluated

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest. Note. "TS-1 (equation (6))" denotes the cracking risk coefficient of concrete calculated according to equation (6); "TS-1 (equation (34))" denotes the cracking risk coefficient of concrete calculated according to equation (34); "C" denotes cracking of concrete; "PC" denotes possible cracking of concrete; "NC" denotes not cracking of concrete.