3.1. Compressive strength and modulus of elasticity
For Compressive strength tests, cylindrical concrete specimens of 100-mm diameter and 200-mm height were cast according to JIS A1108: 2006—method of test for compressive strength of concrete [
13], and ASTM C39/C39M-18: 2018—standard test method for compressive strength of cylindrical concrete specimens [
14].
After casting, the specimens for steam curing were stored in a temperature-controlled room at 20°C and 80% R.H. After 2 hours of casting time, the specimens for 1, 14, and 91 days of compressive strength test were subjected to steam curing. The temperature increased gradually during 2 hours and 15 minutes, up to 65°C as a stable maximum temperature for 3 hours, the target temperature used in the experiments was similar to steam curing in general precast concrete factories. The temperature gradually reduced during 10 hours to stable at 20°C for 6 hours and 45 minutes, the total time for the steam curing procedure was 24 hours.
Figure 1 shows the general temperature regime for a typical steam curing profile.
After steam curing, the specimens were demolded after 24 hours and compressive strength was measured. Then the steam-cured specimens were subjected to drying conditions in a temperature-controlled room at 20°C and 60% R.H. to observe compressive strength at 14 and 91 days.
The other cylinder specimens were stored in water at 20°C to observe 7, 28, and 91 days of compressive strength under standard conditions.
The modulus of elasticity of concrete was also measured on the cylinder specimens used for compressive strength tests.
3.2. Autogenous and Drying shrinkage
Autogenous shrinkage identified as microscopic shrinkage happens in concrete after the initial setting because of cement hydration [
15]. Concrete has drying shrinkage as a result of the different levels of moisture inside and outside the concrete [
16]. Embedded strain gauges connected to a data logger were positioned in the center of beam concrete specimens sized 100mmx100mmx400mm, and they were used to measure concrete strain and temperature during the experiments.
Figure 2 shows the cross and side sections in the mold. Immediately after casting, the specimens were stored in a temperature-controlled room at 20°C and 80% R.H. Then after 2 hours, they were subjected to steam curing protocol as shown in
Figure 1. All specimens were demolded after 24 hours. Then the accelerated cured specimens were subjected to drying conditions in a temperature-controlled room at 20°C and 60% R.H. The strain and temperature were measured in all specimens up to 230 days.
The initial setting time, which is defined as the starting point of autogenous shrinkage, was 3.5 hours after casting, the concrete’s thermal expansion coefficient was assumed to be 10x10
-6/°C in accordance with JCI (2016) [
17], and the total strain from which the thermal strain was subtracted to find autogenous shrinkage was observed.
3.4. Thermal stress analysis using three-dimensional finite element method (3-D FEM)
During the cement hydration process at an early concrete age, thermal cracking occurs due to restrained temperature deformations caused by excessive temperature differences within a massive concrete member of the structure or outer restraint from other attached structural members [
18]. This thermal difference generates tensile stresses in concrete. As a rule, when external restraint is predominant, cracks penetrating through a concrete section (through cracks) are formed [
17]. Guidelines for control of cracking of mass concrete JCI (2016) highlights some other factors that lead to thermal cracking in concrete structures including volume change due to heat of cement hydration, autogenous shrinkage, combined erects of the type of structure, boundary conditions, materials, mixture proportions, construction method, weather conditions, etc. The above factors are prominent in mass concrete structures. However, in thin concrete members, the effects of heat of cement hydration are low due to the insignificant temperature gradients that enable easy and uniform heat dissipation into the surrounding.
JCI guidelines provide an indication of concrete behavior, used as an index of cracking probability, the thermal cracking index, which is the percentage of tensile strength from the maximum principal stresses. The chance of cracking is higher when the thermal cracking index is low. when the cracking index value is 1.0 the probability of cracking is 50%, and when it decrees to be less than 0.6 the cracking probability will increase to be 100%. Equation (1) was used to calculate the thermal cracking index [
17]:
where
Icr: thermal cracking index,
ft(te): splitting tensile strength (N/mm
2),
σmax(te): maximum principal stress (N/mm
2), and
te : temperature adjusted age (day).
In this study, thermal cracking index Icr includes the influence of thermal, autogenous, and drying shrinkage in different ratios, regarding to that, from here on, it will be referred as cracking index.
In this study, thermal analysis and thermal stress analysis were conducted by 3-D FEM for a precast concrete product for mix proportions (N 45%) and (A+FA 45%) using the modified fly ash cement proposed in this study which is subjected to steam curing (which is shown in
Figure 1), then air curing conditions for 6 months (at 20°C and 60% R.H.).
For FEM analysis, a precast concrete model similar to a standard precast box culvert was selected from JIS A 5372- 2016 - precast reinforced concrete products [
19].
Figure 3 shows the actual design of a culvert box and quarter portion for numerical analysis, computer program of thermal stress analysis for mass concrete structures (JCMAC-3) proposed by JCI, which is a 3-D finite element method (FEM) simulation tool, was used. In order to simulate the time-dependent distributions of temperature and relative humidity in a concrete member, the transient heat transfer and moisture transfer analyses based on diffusion equations were carried out by 3D FEM analysis. In thermal stress analyses, autogenous shrinkage strain was added to thermal strain and linear elastic analyses were carried out, in which stress relaxation due to creep of concrete was taken into account by using the reduction coefficient of elastic modulus as shown below.
Dr. Ishikawa proposed a drying shrinkage model using the capillary tension theory for the pore size distribution of the hardened cement paste to calculate the drying shrinkage strain of unrestrained concrete element from its relative humidity [
20]. In this study, Ishikawa’s model was adopted to calculate the restraint stress due to drying shrinkage. In FEM analysis for mix proportion (N 45%) the values of the coefficients in the model to consider the influence of materials used and concrete mix proportion on drying shrinkage were determined so that calculated restraint stress coincided with the observed values by the stress release method [
21]. Drying shrinkage strain of concrete specimens was also measured and the coefficient for mix proportion (A+FA 45%) was decided as the coefficient for mix proportion (N 45%) multiplied by the ratio of drying strain of (A+FA 45%) to that of (N 45%).
In the FEM analysis, the real properties of the concrete which were obtained as experimental results such as adiabatic temperature rise, compressive strength, modulus of elasticity, autogenous and drying shrinkage, and casting temperature were used. Splitting tensile strength was calculated utilizing compressive strength data and constants as specified in JCI guidelines for Equation (2) [
17]:
where
ft(te): splitting tensile strength of concrete at
(te) (N/mm
2 ),
fc(te): compressive strength of concrete at
(te) (N/mm
2 ),
te: temperature adjusted age (days), C
1 = 0.13, C
2 = 0.85,
Some other properties were also obtained from JCI Guidelines [
17] such as specific heat (1.15 J/g °C) and coefficient of thermal conductivity (2.7 w/m °C). Poisson’s ratio was at 0.23, creep of concrete influence was evaluated by using the effective modulus of elasticity, which was obtained by multiplying modulus of elasticity by a reduction coefficient using Equation (3).
where
Ee(te): effective modulus of elasticity of concrete at
(te),
ϕ(te): reduction coefficient for modulus of elasticity due to creep,
EC (te): modulus of elasticity of concrete at
(te).
At the early age during hardening process until reaching the maximum temperature the reduction constant ϕ(
te) was taken to be 0.42, and then it increses during one day of reaching the maximum temperature to be 0.65 as a stable value to the later ages as recommended by JCI guidelines [
17].
For equations (1), (2), and (3) the temperature-adjusted age
(te) can be calculated by using Equation (4) provided by JCI [
17].
where ∆ti: Period of constant temperature continuing in concrete (day), T(∆ti): concrete temperature for ∆ti (°C), and T0 is 1 °C.