Heat of Hydration Analysis

February 7, 2022


In a certain concrete structure with considerable mass or where a construction progresses rapidly with a number of construction joints, the rate and amount of heat generation due to hydration are important. Non-uniform thermal expansion and contraction due to heat of hydration and cooling of concrete accompanied by changing constraints create undesirable stresses. The stresses may cause detrimental cracking in the concrete, thereby reducing its strength and durability.


Heat of hydration analysis thus becomes important when casting mass concrete structures. It enables us to predict and control temperature and stress distribution within a structure to avoid potential problems. Mass concrete structures requiring heat of hydration analysis depend on their dimensions, shapes, cement types and construction conditions.


In practice, hydration analyses are normally carried out for slabs or mats in excess of 800~1000mm in thickness and walls confined at bottom in excess of about 500mm. Surface cracking may develop initially due to the temperature difference between the surface and center. Through-cracks can also develop as a result of contraction restrained by external boundary conditions in the cooling process of high heat of hydration.


The heat of hydration analysis is largely classified into several sub-analyses. It entails temperature distribution analysis for conduction, convection, heat source, etc.; change in modulus of elasticity due to curing and maturity; and stress analysis for creep and shrinkage. The following outlines the various components affecting the analysis.



1. Heat Transfer Analysis

2. Thermal Stress Analysis

3. Procedure for Heat of Hydration Analysis


Heat Transfer Analysis

MIDAS/Gen calculates changes in nodal temperatures with time due to conduction, convection and heat source in the process of cement hydration. The following outlines pertinent items considered in MIDAS/Gen and some of the main concepts in heat transfer analysis:



Conduction is a type of heat transfer accompanied by energy exchange. In the case of a fluid, molecular movements or collisions and in the case of solid, movements of electrons cause the energy exchange from a high-temperature zone to a low-temperature zone. The rate of heat transfer through conduction is proportional to the area perpendicular to heat flux multiplied by the temperature gradient in that direction (Fourier’s law).

Heat of Hydration Analysis EQ1

where,Heat of Hydration Analysis eq2

In general, the thermal conductivity of saturated concrete ranges between 1.21~ 3.11, and its unit is kcal/h⋅m⋅℃. The thermal conductivity of concrete tends to decrease with increasing temperature, but the effect is rather insignificant in the ambient temperature range.



Convection is another form of heat transfer whereby heat is transmitted between a fluid and the surface of a solid through a fluid’s relative molecular motion. Heat transfer by forced convection occurs in the case where a fluid is forced to flow on a surface such that an artificial fluid current is created.

If the fluid current is naturally created by a difference in density due to a temperature difference within the fluid thereby inducing a buoyancy effect, the form of heat transfer is referred to as free convection. Because the fluid’s current affects the temperature field in this type of heat transfer, it is not a simple task to determine the temperature distribution and convection heat transfer in practice.

From an engineering perspective, the heat transfer coefficient, hc is defined to represent the heat transfer between a solid and a fluid, where T  represents the surface temperature of the solid, and the fluid flowing on the surface retains an average temperature, T.

Heat of Hydration Analysis EQ3

The heat transfer coefficient (hc) widely varies with the current type, geometric configuration and area in contact with the current, physical properties of the fluid, average temperature on the surface in contact with convection, location and many others, and as such it is extremely difficult to formulate the coefficient. In general, convection problems associated with temperature analyses of mass concrete structures relate to the type of heat transfer occurring between the concrete surface and atmosphere. Accordingly, the following empirical formula is often used, which is a function of an atmospheric wind speed.

 Heat of Hydration Analysis EQ4

The unit for heat transfer coefficients (Convection coefficients) is kcal/m2⋅h⋅℃.



Heat Source

Heat source represents the amount of heat generated by a hydration process in mass concrete. Differentiating the equation for adiabatic temperature rise and multiplying the specific heat and density of concrete obtain the internal heat generation expressed in terms of unit time and volume. Adiabatic conditions are defined as occurring without loss or gain of heat; i.e., as isothermal.

Internal heat generation per unit time & volume (kcal/m3⋅h)

Heat of Hydration Analysis EQ5

Equation for adiabatic temperature rise (℃)

Heat of Hydration Analysis EQ6where,

T∶ Adiabatic Temperature (℃)
K∶ Maximum Adiabatic Temperature Rise (℃)
α∶ Response Speed
t∶ Time (days)




Pipe Cooling

Pipe cooling is accomplished by embedding pipes into a concrete structure through which a low-temperature fluid flows. The heat exchange process between the pipes and concrete reduces the temperature rise due to the heat of hydration in the concrete, but increases the fluid temperature. The type of heat exchange is convection between the fluid and pipe surfaces. The amount of the heat exchange is expressed as follows:

Heat of Hydration Analysis EQ7





  1. Ambient Temperature

Initial temperature is an average temperature of water, cement and aggregates at the time of concrete casting, which becomes an initial condition for analysis.

  • Ambient Temperature

Ambient temperature represents a curing temperature, which may be a constant, sine function or time-variant function.

  • Prescribed Temperature

A prescribed temperature represents a boundary condition for a heat transfer analysis and always maintains a constant temperature. The nodes that are not specified with convection conditions or constant temperatures are analyzed under the adiabatic condition without any heat transfer. In a symmetrical model, the plane of symmetry is typically selected as an adiabatic boundary condition. The basic equilibrium equations shown below are used for heat transfer analysis. Analysis results are expressed in terms of nodal temperatures varying with time.

Heat of Hydration Analysis EQ8-1

Heat of Hydration Analysis EQ8-2where,
T∶ Nodal Temperatue
ρ∶ Density
c∶ Specific Heat
kxx kyy kzz ∶ Heat Conductivity
h∶ Convection Coefficient
Q∶ Rate of Heat Flow-Quantity of Heat Penetrating Per Unit Time
q∶ Heat Flux-Quantity of Heat Penetrating a Unit Surface Area Per Unit Time




Thermal Stress Analysis

Stresses in a mass concrete at each stage of construction are calculated by considering heat transfer analysis results such as nodal temperature distribution, change in material properties due to changing time and temperature, time-dependent shrinkage, time and stress-dependent creep, etc. The following outlines some important concepts associated with thermal stress analysis and pertinent items considered in MIDAS/Gen.


Equivalent Concrete Age Based on Temperature and Time & Accumulated Temperature

Change in material properties occurring from the process of maturing concrete can be expressed in terms of temperature and time. In order to reflect this type of phenomenon, equivalent concrete age and accumulated temperature concepts have been incorporated. Equivalent concrete age is calculated on the basis of CEB-FIP MODEL CODE, and the Ohzagi equation is adopted for calculating accumulated temperature, which bases on a maturity theory.

- Equivalent Concrete Age as per CEB-FIP MODEL CODE

Heat of Hydration Analysis EQ9.1-1

- Ohzagi’s Equation for Accumulated Temperature

Heat of Hydration Analysis EQ10

Concrete Compressive Strength Calculation Using Equivalent Concrete Age and Accumulated Temperature



Heat of Hydration Analysis EQ11-1


Heat of Hydration Analysis EQ12

- Ohzagi’s Equation

Heat of Hydration Analysis EQ13

- KS Concrete Code (1996)

Heat of Hydration Analysis EQ14




  • Deformations resulting from temperature changes

Thermal deformations and stresses are calculated by using the nodal temperature changes at each stage obtained through a heat transfer analysis.

  • Deformations due to shrinkage

Additional deformations and stresses develop due to shrinkage after initial curing. MIDAS/Gen adopts ACI CODE and CEB-FIP MODEL CODE to include the shrinkage effects in thermal stress analyses, which reflect the cement type, structural configuration and time.

  • Deformations due to creep

Additional deformations and stresses develop as a result of sustained stresses in concrete structures. MIDAS/Gen adopts ACI CODE and CEB-FIP MODEL CODE to consider the effects of creep.




Procedure for Heat of Hydration Analysis

  1. Select Model > Properties > Time-Dependent Material (Creep/Shrinkage) and Time Dependent Material (Comp. Strength), and specify the time-dependent material properties. Link the general material properties and time-dependent material properties in Model > Properties > Time-Dependent Material Link.
  2. Enter the relevant data required for Heat of Hydration Analysis in the submenus of Load > Hydration Heat Analysis Data.
  3. Enter the Integration factor, Initial temperature, stress output points and whether or not to consider the effects of creep and shrinkage in Analysis > Hydration Heat Analysis Control.
  4. Select the Analysis > Perform Analysis menu or click perform analysis icon Perform analysis.
  5. Once the analysis is completed, check the results in contours, graphs, animations, etc.
Model of a Pier Cap of an Extradosed Prestressed Concrete Box for Heat of Hydration Analysis
Figure 1. Model of a Pier Cap of an Extradosed Prestressed Concrete Box for Heat of Hydration Analysis Reflecting the Concrete Pour Sequence
Heat properties and Time-Dependent Material Properties
Figure 2. Heat properties and Time-Dependent Material Properties Dialog Box
Construction Stage Dialog Box to Reflect the Concrete Pour Sequence
Figure 3. Construction Stage Dialog Box to Reflect the Concrete Pour Sequence (Element, Boundary, and Load Groups are Defined)
Graphs of Analysis Results for Each Construction Stage
Figure 4. Graphs of Analysis Results for Each Construction Stage


Subscribe S.O.S Newsletter

Author Information
Young-Wook Hur | Technical Manager, MIDAS IT

Young-Wook has over 8+ years of experience in building design, especially apartment, complex and palnt structures. He is currently working at MIDAS IT as a technical manager, mainly assisting with various projects related to structural engineering.

E-BOOK Heat of Hydration

Heat of hydration analysis is important when casting mass concrete structures,
which is classified into several sub-analyses.

Read all