Study on Adiabatic Temperature Rise Reflecting Hydration Degree of Concrete

.e thermal model and the relevant parameters of concrete are the most important issues to study the space-time characteristics of temperature field, which are also the theoretical foundation of temperature control and crack prevention for the mass concrete structures. In this research, the improved adiabatic temperature rise test is carried out, and the temperature variation of fly ash concrete is analyzed. Furthermore, a thermal model of concrete considering the hydration degree is established based on the existing achievements. Meanwhile, the thermal conductivity and specific heat of concrete are measured via three approaches: by treating the parameters as constant values, by computing the parameters as variables of the degree of hydration, and by backanalyzing the parameters through BP neural network. Finally, the thermal parameters determined by different methodologies are substituted into the thermal model, respectively, and the finite element analysis of the concrete specimen is performed. By comparing simulated temperatures with various measured results, it can be found that the numerical analysis results of parameters calculated by BP neural network are closest to the measured values in the whole curing ages..erefore, BP neural network method is an effective way to calculate the thermal parameters, and BP inversion algorithm provides a new way for accurately study the temperature profile of mass concrete structures.


Introduction
Cracking is a common phenomenon in the concrete structures.Even some buildings which were carefully designed and constructed also appear with cracks.Among the numerous factors that cause the cracking of concrete, the temperature stress has been generally taken as the uppermost reason [1,2].In the early ages, the temperatures of concrete change quickly with the rapid hydration of cement, and the temperature distributions in concrete structures are uneven because of the varied heat boundary conditions.e changes of both the temperature and the boundary conditions can lead to the obvious temperature gradients of the concrete structures, and the temperature gradients will lead to the temperature stress.e larger the temperature gradients are, the greater the temperature stress will become.Worst of all, the macrocracks will occur, and some of them will continue to expand and eventually form the penetrating cracks when the stress exceeds the tensile strength of concrete.And these penetrating cracks destroy the integrity of concrete structures and even become the potential hazards of the concrete constructions [3].e changes of concrete temperature must be paid great attention in the design and construction process [4][5][6][7][8][9][10][11].It is well known that the thermal model and the relevant parameters of concrete are the basic conditions for studying the space-time character of the temperature field and also are the prerequisite for accurate calculation of the temperature stress.e determinations of the thermal model and the relevant parameters are the key issues of studying the temperature.erefore, it is necessary to carry out a thorough research.
In the simulation of concrete temperature field, the thermal model and its related parameters have big influences on the reliability and accuracy of the calculation results.e important factors affecting the thermal model include thermal conductivity, specific heat, adiabatic temperature rise, and surface heat dissipation coefficient, and so on.e above parameters may be acquired by experimental methods or empirical formulas in most cases [4,5].However, the temperatures of concrete in the hardening process are affected by many complex factors, and the environmental condition of the laboratory is different from that of the construction site, both of which make it difficult to determine parameters precisely.erefore, many researchers had adopted the advanced back analysis method to solve the problem [6][7][8].But there were common shortcomings in the experimental method and in the numerical inversion method; thermal conductivity, specific heat, and other thermal parameters of concrete were regarded as constant values, and the changes of the concrete thermodynamic parameters in the hydration exothermic process were not taken into account.Hence, the temperatures simulated by the traditional thermal model with constant parameters are different from the experimental results because the actual engineering conditions cannot be well reflected.
In order to accurately predict the temperature distribution in the concrete, various factors that affect the temperature must be considered.It is known to all that the process of cement hydration is from the original flow state to plastic state and then to mature solidification.In addition, thermodynamic parameters also change with the transformation of the concrete status.De Schutter and Taerwe [9] studied the evolution of thermal parameters by means of experiments and literature summaries and indicated that thermal conductivity and specific heat of the early age concrete changed with the hydration degree.Kim et al. [10] researched the effects of water-cement ratio, curing age, temperature, and humidity on thermal conductivity through a series of experiments.Cui et al. [11] carried out some experiments to study the temperature change rules of the early-age concrete, introduced the thermal conductivity and specific heat considering the hydration degree, and simulated the temperature field by the finite element method.
ose researchers started to know the relevant rules of changing the concrete thermal parameters and applied them into the thermal performance analysis of concrete structures.However, most existing researches focus on the early stage of concrete hydration, while the research on the real temperature evolution of concrete during the 28-day curing ages is not enough.erefore, the thermal parameters need to be further refined.
Few researchers have addressed the thermal parameters considering the hydration degree during the whole 28 d curing ages.Over the years, the thermal parameters had been regarded as constant values when the temperature field of concrete was simulated, and the influence of the hydration degree on these parameters was ignored in most cases.And with the further study, the thermal parameters considering the hydration degree were only used to accurately simulate the early 7-day concrete hydration [11].
erefore, it is urgent to research the thermal parameters considering the hydration degree during the entire curing ages, in order to provide suitable parameters for the thermal model and simulate the temperature field of concrete structure accurately.
Based on above facts, the adiabatic temperature rise tests of fly ash concrete with different water-cement ratios were firstly performed, the real-time temperatures at the testing points in the specimen were recorded during the whole curing ages, and then the temporal and spatial distribution of the concrete temperature field was analyzed carefully.A concrete thermal model considering the hydration degree was established according to the experiment and the available research achievements.Based on the experimental data, several important thermal parameters considering the hydration degree were retrieved by BP neural network to ensure that the thermal parameters are in line with the actual values of concrete hydration.Finally, the results from three different algorithms were compared and analyzed thoroughly, showing that the temperatures obtained from the approach proposed in this paper are closer to the experimental data than others.

Adiabatic Temperature Rise Test
In this study, the adiabatic temperature rise tests were carried out to accurately evaluate the effect of hydration heat on the temperature of fly ash concrete during the curing ages.

Raw Materials.
An ordinary Portland cement with compressive strength of 42.5 N/mm 2 and the first-grade fly ash were selected as the cementing materials; the two materials complied with the Chinese National Standards GB175-2007 and GB/T 50146-2014, respectively, and the mass ratio of them was 1 : 1. e oxide contents and physical properties of the two materials are tabulated in Table 1.Manufactured sand with a fineness modulus of 3.05 was selected as the fine aggregate, and the coarse aggregate was the crushed limestone of a continuous gradation (5 mm∼80 mm with a maximum nominal size of 80 mm).
Ng et al. [12] have presented that water-cement ratio (w/c) has a great influence on the hydration degree, and the same hydration degree can be achieved when w/c is close to 0.48.It has been proved by Mill too.So two concrete specimens (named as C1 and C2) with w/c of 0.5 and 0.45 were prepared.e specific mix proportion of concrete in the test is indicated in Table 2.

Experimental Schemes.
e adiabatic temperature rise chamber, provided by the Yangtze River Scientific Research Institute of China, was used to carry out the adiabatic tests of fly ash concrete.e specimen was a concrete cylinder with a diameter of 380 mm and a height of 360 mm.In order to describe the temperatures of different positions in the concrete specimen, 13 testing points (T 0 , T 1 , T 2 , . .., T 12 in C1 and T 13 , T 14 , T 15 , . .., T 25 in C2) were evenly arranged into 5 layers along the height direction of the concrete specimen (Figure 1).Micro K thermocouple detectors were adopted to provide quick and accurate temperature measurement, and a multichannel temperature measuring device was used to record the data according to the set intervals automatically.

2
Advances in Materials Science and Engineering e temperature sensors xed on a welded iron frame were wrapped with the epoxy resin asphalt to prevent damage during the concreting process and were put into the cylinder container as soon as the asphalt dried.e concrete was carefully poured through a funnel to ensure that the positions of sensors were not shifted, and at the same time, the concrete specimen was vibrated on the platform vibrator, and then the specimen was put into the outer cylinder after the surface was attened.e whole experimentation was carried out under the room temperature of 25 °C, and the temperatures of the laboratory were adjusted timely according to the concrete specimen by air conditioner in order to ensure that the ambient temperatures are in accordance with the test temperatures during the adiabatic test.Figure 2 shows the concrete pouring, the insulation measures of the adiabatic chamber, and the temperature acquisition system in turn.

e Temperature of Specimen with Di erent Water-Cement Ratios.
In each of the two specimens, three testing points were selected as representatives, and the historic curves of adiabatic temperature rise were plotted in order to give out the changes of temperature over time.Figure 3 presents the temperatures of the testing points in the two specimens which have taken account into the heat loss of the adiabatic chamber.It can be seen from Figure 3 that the temperatures of the both specimens are increasing with the increase of curing age, and the concrete is rapidly hydrated in the rst few hours to reach a higher temperature, after which the rate of hydration becomes slow and the rate of temperature rise diminishes signi cantly.Meanwhile, we can also acquire from Figure 3 that the temperatures of C2 are lower than those of C1 in the whole process on account of the lower w/c.Moreover, the rising rates of C1 and C2 are slightly di erent; the temperature rising of C1 is faster than that of C2 at the beginning, and the increasing rates of the both are opposite at the later stage.e w/c of C1 is relatively large which can provide adequate conditions for the hydration of cement leading to a faster hydration rate at the beginning while the w/c of C2 is smaller, and the cement cannot contact with water fully and the hydration rate is slow accordingly.With the curing ages increasing, the hydration of C1 is completed basically, while, by contrast, the internal temperature of C2 increases with the cement hydration, which is an important factor that promotes the hydration rate of cement, so the rate of C2 is slightly larger than that of C1 at a later period.3 shows the temperatures at the testing points in the specimens C1 and C2 at di erent ages, and the spatial distribution of temperature is nonuniform obviously.e data of the testing points T 0 , T 6 , and T 7 and T 25 , T 18 , and T 19 indicate that the temperature near concrete core is higher while it slightly decreases along the radius at the same curing ages, and this is mainly because that the heat loss compensation of the equipment is not su cient.

e Temperature at Di erent Positions of Specimen. Table
Figure 4 shows the temperatures along the height section of 1 d, 3 d, 7 d, 14 d, and 28 d after pouring, taking C1 as an example.Both of them manifest that the temperatures of the testing points vary a lot in the rst curing days.e temperature of testing points near concrete core is highest among the ve points at the same moment because more heat is dissipated through the top and bottom of the adiabatic chamber.What is more, the temperatures of the bottom are lower than those of the top in the early hydration process, but it' is the opposite after 3 days.It is just due to the fact that the top of the concrete in the pouring process is exposed to the air and the ambient temperature is higher,  Advances in Materials Science and Engineering which improves the temperature of the top because of the heat transfer and leads to a faster hydration.As a result, the temperature of the top is higher at the early stage than later.After pouring, the external conditions of the top and the bottom are the same while the heat insulation effect of the bottom is better than that of the top during the hydration process.From Figure 4, it also can be shown that all the temperature changes have a similar trend, and the difference of the temperatures of the testing points becomes small from the third curing day.

Thermal Governing Equations considering the Hydration Degree
e appropriate thermal governing equations considering the hydration degree must be required in the numerical analysis.is section is intended to investigate a theoretical model which will be implemented in the numerical method.

ermal Governing Equations considering the Hydration
Degree.A cylinder specimen is frequently used in the adiabatic temperature rise test, so a calculation model of transient temperature field for concrete in cylindrical coordinates (r, φ, z) needs to be established.According to the energy conservation law and the second law of thermodynamics, a domain undergoing thermal changes considering the hydration degree is subjected to the following field equation: where T is the temperature ( °C), τ is the curing age (d), ρ is the density of concrete (kg/m 3 ), λ(α) is the thermal conductivity at the different hydration degree (kJ/(m•d• °C)), c(α) is the specific heat at the different hydration degree (kJ/(kg• °C)), and Q(α) is the released heat associated with an internal source and is provided by the hydration reaction itself (kJ/(m 3 •d)).e critical issues to simulate the temperature field of concrete is how to determine the hydration degree α, thermal conductivity λ(α), specific heat c(α), and the released heat Q(α).Among the above parameters, the calculation of the hydration degree is the first step in determination of other three parameters.

e Hydration Degree.
e solidifying process of concrete is accompanied by significant temperature changes, and the temperature rise is caused by the heat release of cementing materials reacting with water.rough the amount of the heat release, it is possible to calculate how much the hydration reaction has been completed at a certain curing age [13], which can be de ned as the hydration degree.e ambient temperature has a signi cant e ect on the hydration rate of concrete during the curing process, but it does not a ect the maximum adiabatic temperature rise of concrete.Many thermal models adopt the concept of "equivalent age" to simplify the in uence of temperature on the hydration rate.Equivalent age of the concrete is an arti cial age which is concluded according to the actual temperature history, and it is determined through judging the result that the temperature has accelerated or delayed the aging process.For a certain type of concrete, there is a unique relationship between the adiabatic temperature rise and the equivalent age.
Saul [14] puts forward that regardless of the combination of time and temperature, the thermal performance of different positions in concrete specimen should be exactly identical as long as the maturity of concrete is the same.Based on the Arrhenius law which describes the e ect of the temperature on the rate of chemical reaction, Professor Bazant [15] proposes the equation of the equivalent age, which can be used to calculate the hydration quantity in di erent hydration temperatures.
e expression of the equivalent age is shown in the following: where t e designates the equivalent age (d), E A is the reaction activation energy (J/mol), R is the ideal gas constant (J/(K•mol)), T r is the reference temperature (K), and T is the actual temperature at any point of concrete (K).Refer to the literature [11], E A /R 2700 K, T r 20 °C, that is, 293 K. e thermal properties of concrete such as elastic modulus, hydration reaction, speci c heat, thermal conductivity, and hydration heat are related to the hydration degree by which all of these indicators can be expressed.e models used frequently to calculate the hydration degree include complex exponential, hyperbolic, and exponential model.
e hyperbolic model selected in this paper is expressed as follows: where α(t e ) is the hydration degree of concrete in the equivalent age (d) and both m and n are coe cients.From the experimental results, it can be seen that the temperatures of concrete at di erent positions are similar, so the temperature of the concrete core is selected as the representative to study the relationship between the hydration degree and the equivalent age, shown in Figure 5.According to the measured temperatures expressed by black dots in Figure 5, the parameters m and n involved in (3) can be determined.For the specimen C1, m 0.97 and n 2.78; for the specimen C2, m 0.98 and n 3.82.

Heat Generation Rate.
From (1), it can be seen that the hydration heat Q(α) is the sole heat-generated source leading to the temperature rise of concrete under adiabatic condition.erefore, the heat generation rate of concrete considering the hydration degree under di erent temperatures can be formulated as where Q ∞ is the ultimate hydration heat quantity of concrete (kJ/m 3 ).

Determination of the Thermal Parameters
In order to obtain the accurate temperatures of concrete during the curing ages, not only the governing equation of heat transfer needs to be established in the numerical simulation but also the thermal parameters based on the hydration degree of concrete need to be carefully taken into account.In this section, the relevant parameters of concrete were calculated by three methods listed below, respectively.Advances in Materials Science and Engineering taken as constant values in some existing researches, and these values could be obtained by the test or by the empirical formula.e thermal conductivity and speci c heat of y ash concrete were determined by the method recommended in the literature [2] with the consideration of the ambient temperature.Moreover, the adiabatic temperature rise of y ash concrete was tted by the hyperbolic model, and the speci c parameters are given in Table 4.

ermal Parameters considering the Hydration Degree:
Method 2. e thermal parameters in Table 4 characterize the properties of the solidi ed concrete, and the changes of thermal parameters during the concrete hardening cannot be re ected.In order to predict the concrete temperatures of the early curing ages more accurately, the variation of thermal properties during the hardening process should be taken into account.

ermal Conductivity Coe cient.
e thermal conductivity coe cient is an indicator which describes the heat capacity of concrete in the temperature eld calculation, and it is a ected by water-cement ratio, cement content, density,   Advances in Materials Science and Engineering temperature and humidity, and so on [9,10,16,17].e formula of thermal conductivity coefficient proposed by Schindler is widely used, in which the effect of the hydration degree on thermal conductivity coefficient is taken into account [18], as shown in the following equation: where λ u denotes the thermal conductivity coefficient of hardened concrete (kJ/(m•d• °C)).e parameter is further found to be linear with the temperature [2,10].Kim et al. [10] put forward a formula for correcting the thermal conductivity coefficient at different pouring temperatures with water-cement ratio above 0.4: where λ T denotes the correction of the thermal conductivity coefficient at different pouring temperatures (kJ/(m•d• °C)) and T 1 denotes the pouring temperatures ( °C).
According to the abovementioned equation, an integrated relationship about the thermal conductivity coefficient of concrete, pouring temperatures, and the hydration degree can be obtained: where λ(α, T 1 ) is the thermal conductivity coefficient at different pouring temperatures and the hydration degree (kJ/(m•d• °C)).

Specific Heat.
It has been pointed out that the temperatures and content have an important effect on specific heat [9,19], and Van Breugel [20] further proposes a calculation model of specific heat considering temperature, mixing ratio, and the hydration degree.From the formula, it can be seen that the specific heat decreases linearly with the hydration degree: where ω c , ω f , ω a , and ω w are the mass of cement, fly ash, aggregate, and water in per cubic meter of concrete, respectively (kg/m 3 ); c c , c f , c a , and c w are the specific heat of cement, fly ash, aggregate, and water, respectively (kJ/(kg• °C)); c cef is the artificial specific heat of cement (kJ/(kg• °C)); and T j is the actual temperature of cement at time t j ( °C). e equations of thermal conductivity coefficient and specific heat mentioned above are suitable to the early age concrete.In method 2, thermal parameters were calculated by ( 7) and ( 8) in the early curing ages; after 7 days, the corresponding parameters of concrete were taken as constant values.e parameters involved are listed in Table 5.

4.3.
ermal Parameters Inverted by BP: Method 3. e thermal parameters of concrete are affected by many factors such as ambient temperature, water-cement ratio, cement content, and aggregate content, and the influence of these factors may be dependent on each other.As a result, the thermal parameters are difficult to be accurately determined by the experiment.Furthermore, the incorporation of the fly ash can reduce the adiabatic temperature rise of concrete, which is due to the pozzolanic reaction of the fly ash, so the influence of the fly ash should be considered when the temperature of the fly ash concrete is calculated.However, the effect of the fly ash dosage on the calculation results cannot be quantitatively given.Nevertheless, BP neural network provides an effective approach for such complex situations.In this section, the thermal conductivity coefficient, specific heat, and the released heat considering the hydration degree will be predicted by BP neural network.

Target Function for BP Neural Network.
e BP inversion method is to establish the relationship between the thermal parameters and the measured temperatures.e measured temperatures are the actual values obtained in the adiabatic test, which are the results of combined action of various factors including the pozzolanic reaction of the pulverised fuel ash.It can be expressed as where T denotes the temperature ( °C) and g(x) is the function of the various factors.And the difference between the measured temperatures and the theoretical ones can be calculated by the following equation: where x denote the factors that a ect the temperature, such as c(α), λ(α), and Q(α); α and u i are the measured temperatures; g i (x) is the theoretical temperature; and ϕ(x) is the di erence.
If the minimum value of (10) can be acquired through adjusting the range of a ecting factors, then the theoretical values will match well with the measured ones.erefore, the parameter inversion becomes an optimization problem in essence [21,22].

Back Analysis Scheme.
According to the mixing ratio of concrete and practical experience in engineering, the ranges of thermal parameters were determined.Considering the temperature changes of the di erent positions in the adiabatic temperature rise test, the temperatures of seven testing points from the top, bottom, and center of the specimens were selected as the representatives; the speci c points were T 0 , T 3 , T 6 , T 7 , T 12 , T 1 , and T 10 of the specimen C1 and T 13 , T 15 , T 18 , T 19 , T 22 , T 24 , and T 25 of the specimen C2.
e corresponding points of the above sensors in the temperature eld calculation were chosen, and the temperatures at 21 di erent curing ages were selected at the same time; then 147 sets of training samples which were composed of thermal parameters and temperatures within the parameter range were generated for training networks (shown in Table 6).
e samples were rst normalized when BP neural network was trained, and the temperatures and four thermal parameters were taken as the input values and the outputs, respectively, after which the network mapping relationship of the temperatures to thermal parameters could be achieved.In addition, (11) was used to estimate the hidden layer nodes, and the repeated trail calculations proved that the training error was minimal when the number of nodes was 10.
where n h denotes the hidden layer nodes, n i denotes the input layer nodes, n o denotes the output layer nodes, m is the constant, and m ∈ [1,10].

Inversion Result.
e suitable weights between the input layer and the hidden layer and the hidden layer and the output layer were attained through training BP neural network, and the temperatures of T 0 , T 1 , and T 10 and T 25 , T 13 , and T 22 measured by the adiabatic test were taken as the input values; then the thermal conductivity coe cient, speci c heat, the released heat, and the hydration degree      Advances in Materials Science and Engineering were predicted.At last, the thermal parameters required were obtained after all the data were processed through antinormalization, listed in Table 7.

Finite Element Analysis.
e specimen used in Section 2 was simulated by finite element method, and the numerical model was meshed with 24000 hexahedral solid elements and 26061 nodes.e initial pouring temperature was 26.9 °C, and the ambient temperature was 25 °C.Moreover, the time step in the computation should be small enough in order to capture the released heat of the very early age in the sample.Previous studies suggest that the increase of time step within 0.2 day to 1 day is fine enough to achieve the required accuracy [23].
erefore, the time step of 1 day was selected over the whole 28 days in the temperature field analysis.
All the surfaces of the sample should be insulated from the air, and there is no heat exchange between the specimen and the outside environment in the adiabatic temperature rise test.However, the simulation results showed that the temperature field of concrete was not completely uniform.It has also been proposed in the literature [12] that the adiabatic temperature rise chamber cannot sensitively and accurately track and adjust the temperatures during the entire experiment process, and the temperatures throughout the specimen are inconsistent.erefore, the surface heat transfer coefficient β was chosen to simulate the uneven temperature field of concrete in this study.Considering the insulation material used and the hysteresis of tracking temperature of the adiabatic test chamber, the coefficient β was calculated by the formula mentioned in the literature [2], and the value of β was taken as 24.1 kJ/(m 2 •d• °C).
Figure 6 shows the temperature distribution on the 28th day in a 3D cutoff view.From it, we can see that the temperatures are different within the specimen, and the temperatures near the center are higher than the ones near the concrete surface.Also the maximum temperature difference between the core and the edge is approximately 6.8 °C.Furthermore, the calculation results are consistent with the measured ones.

Result Analysis.
In order to express the temperature variation of the specimen along the radial and the height sections, the temperatures of six positions corresponding to the sensors in the adiabatic temperature rise test were selected.Figure 7 indicates the comparison between the measured temperatures and the numerical calculation results of the specimens with the change of the curing ages at different positions.
(1) In different positions of the concrete, the temperature changing rules found in Method 1, Method 2, and Method 3 mentioned in Section 4 are consistent with the changing rules of the measured temperatures.e increase of temperature gradually slows down after a sharp growth in the initial period.Since a large amount of heat is generated by the initial cement hydration, the temperature increases rapidly in a short time.As the hydration reaction continues, the content of anhydrate cement decreases and the amount of hydration products increases, both of which will hinder the contact between the cement and water, and the increasing rate of temperature slows down.(2) e numerical simulation results of the specimens C1 and C2 are compared, and it can be seen that Method 1 is fitting the specimen C1 well while the fitting deviation for the specimen C2 is larger.e curvature characteristic of the hyperbolic curve itself shows that the rate changes rapidly in the early stage and becomes gentler later.Obviously, the fitting form coincides with the hydration process of the specimen C1 better, but it is not applicable to the specimen C2 which has a significant temperature rise at the later stage.(3) rough the analysis of calculation results obtained from the above three methods, a significant difference among the temperatures of six measuring points can be found; in the initial 7-day stage, the temperature historic curve obtained by Method 2 is the highest while that of Method 1 is the least.Afterward, the temperature historic curves obtained by those three methods are similar, and the temperatures calculated by Method 3 are always closest to the measured values during the whole curing ages.It can be explained that when the hydration degree is considered, the thermal conductivity coefficient is larger, specific heat is lower in the early age, and the temperatures of the same period rise rapidly.With the increase of the curing age, the thermal conductivity coefficient decreases and specific heat increases at the same time; thus, the temperature rise rate slows down.BP neural network in Method 3 is used to invert the thermal parameters which based on the hydration degree and modifies the thermal conductivity coefficient and specific heat of the whole process.erefore, the simulation temperatures are in good agreement with the measured ones.(4) From comprehensively studying the numerical simulation results of the di erent methods, the results of Method 2 are close to the measured values only in the early age, whereas the deviation between the simulation results and the measured ones at the later stage is slightly larger.Moreover, the temperatures calculated by Method 3 are similar to the measured values in the entire curing process.erefore, it is necessary to consider the thermal parameters which change in accord with the hydration degree in the analysis of concrete temperature eld in order to accurately reect the temperature change of concrete hydration process and provide scienti c basis for taking appropriate measures in design and construction of the practical projects.
(5) e results of the adiabatic temperature rise test show that the temperature near concrete core is the highest, and it decreases gradually along the radius, whereas the temperature near the bottom is slightly higher than that near the top.Consequently, the surface heat dissipation coe cient is necessary to be considered in the nite element numerical calculation.

Practical Application of the Proposed Method.
Xiaowan dam is located in Yunnan Province.Block 22# with temperature monitor is selected, and its nite element model is presented in Figure 8.And the related parameters are listed in Table 8. Figure 9 shows the temperatures of tted and measured ones.It is easy to notice that the temperature variation between them is within 3 °C.e good agreement between them further con rms the practice values of the proposed method.

Conclusions
e thermal model and the relevant parameters are the key issues in the study of temperature control and anticracking of the concrete structures.A further research about them has been accomplished through theoretical analysis, experiments, and numerical simulation, and the following conclusions can be drawn: (1) e space-time evolution characteristics of y ash concrete are revealed in this paper.e space distribution of temperatures is uneven, and the higher temperature near concrete core and lower temperature near concrete surface are found in the adiabatic temperature rise test.e temperature keeps increasing in the curing process, and the temperature rise rate in the rst 3 days is faster while the rate at the later stage is slower.
(2) e thermal governing equation of the temperature rise of concrete as a function of the concrete age, heat generation of concrete, thermal conductivity, and speci c heat which are dependent on the degree of hydration is presented, and the dynamic changes of concrete thermal properties during the hydration process are considered thoroughly.(3) e BP neural network method is used to inverse the thermal parameters; the inversion results are substituted into the thermal model, and a practical BP algorithm which considers the degree of pozzolanic reaction of the y ash is established.(4) On the basis of the thermal parameters determined via the three methods mentioned in this paper, the nite element analyses of the temperature pro le of y ash concrete specimens are performed.rough comparing the temperature evolutions which are dependent on the thermal parameters calculated by three di erent methods, it can be found that the results calculated by the BP method are closer to the experimental results than others.

Figure 4 :
Figure 4: e temperature of the specimen C1 at di erent ages along the height section ( °C).(a) e temperature near concrete core.(b) e temperature near concrete surface.

Figure 5 :
Figure 5: e relationship between the hydration degree and the equivalent age.

YFigure 6 :
Figure 6: Temperature distribution of concrete specimen at the 28th day.

Figure 7 :
Figure 7: Historic curve of temperature at speci ed points.

Figure 8 :
Figure 8: e finite element model of the block 22#.

Table 1 :
Chemical composition and physical properties of the cementing materials.

Table 2 :
Mixture proportion of concrete.

Table 3 :
e temperatures of specimens at di erent sections ( °C).

Table 5 :
e thermal parameters involved in Method 2.

Table 4 :
e thermal parameters involved in Method 1.
n i 1

Table 6 :
Sample data of back analysis with BP network.

Table 7 :
e thermal parameters calculated by Method 3.

Table 8 :
e related parameters of the dam. Figure 9: Historic curve of temperatures of block 22#.Advances in Materials Science and Engineering