Skip to main content
Have a personal or library account? Click to login
Beyond Surface Modulus: Shape-Aware Crack-Index Assessment of Mass Concrete Foundations Under Early-Age Thermal Stress Cover

Beyond Surface Modulus: Shape-Aware Crack-Index Assessment of Mass Concrete Foundations Under Early-Age Thermal Stress

Open Access
|Aug 2026

Figures & Tables

Figure 1:

Time-dependent baseline relationships for the hydration-heat and mechanical-property models (K = 51 °C, α = 1.516 day−1, a = 4.0, b = 0.85, fc′(28) = 35 MPa). Values in (b) are normalised to the 28-day reference

Table 1:

Thermal crack-index acceptance bands (Korea Concrete Institute, 2009)

CriterionThermal crack index Icr
To prevent cracksIcr ≥ 1.5
To limit cracks1.2 ≤ Icr < 1.5
To allow cracks0.7 ≤ Icr < 1.2
Table 2:

Baseline material and thermal properties

PropertyConcreteSub-soilUnit
Specific heat c0.250.20kcal·kgf(−1) °C(−1)
Density ρ24001800kgf·m(−3)
Thermal conductivity k2.31.7kcal·m(−1)·hr(−1)·°C(−1)
Convection (atm / steel form)4.5 / 4.5 (vary)12kcal·m(−2)·hr(−1)·°C(−1)
Ambient Tamb27.1 (vary)°C
Casting T015 (vary)°C
28-day compressive strength f′c(28)350kgf·cm(−2)
Strength gain (a, b)4.0, 0.85
28-day modulus Ec2.99×1051.0×10−5kgf·cm(−2)
Thermal expansion αth1.0×10⁻51.0×10⁻5°C(−1)
Poisson’s ratio0.180.20
Cement content350kg·m(−3)
Heat source (K; α)51.0; 1.516°C; day(−1)
Table 3:

Surface-modulus parametric matrix (G-series)

IDGeometryDimensions [m]Acool [m2]V [m3]M [m−1]
G1Small cube2 × 2 × 22082.50
G2Reference cube3 × 3 × 345271.67
G3Large cube4 × 4 × 480641.25
G4Very large cube5 × 5 × 51251251.00
G5Thin slab6 × 6 × 1.572541.33
G6Thick slab6 × 6 × 31081081.00
G7Column2 × 2 × 652242.17
G8Pile cap4 × 4 × 248321.50
Figure 2:

Temperature–time histories at the seven monitoring nodes for the baseline case

Figure 3:

Temperature contour of the quarter model at the instant of peak core temperature

Table 4:

Effect of time step Δt

CaseΔt [h]StepsTmax,core [°C]ΔTcs,max [°C]Icr,minCPU [min]
A10.536063.9325.221.68810
A21.018063.8424.211.73507
A32.0 (ref)9065.1626.171.79303
A44.04567.7327.971.51902
A58.02372.6931.411.28401
Table 5:

Effect of element size

CaseEdge [mm]Grid cellsTmax,core [°C]ΔTcs,max [°C]Icr,minCPU [min]
B1500364.7717.501.46401
B2375464.9821.891.45601
B3300565.1024.491.54302
B4250 (ref)665.1626.171.79303
B5187.5865.2328.231.47008
B61501065.2729.381.39617
Table 6:

Effect of convection coefficient hc.

Casehc [kcal m−2 h−1 °C−1]Tmax,core [°C]ΔTcs,max [°C]Icr,min
C11.565.6416.552.471
C24.5 (ref)65.1626.171.793
C38.064.8729.491.612
C412.064.6831.031.537
Table 7:

Effect of placement temperature T0

CaseT0 [°C]Tmax,core [°C]ΔTcs,max [°C]Icr,min
D115 (ref)65.1626.171.793
D22069.9229.571.647
D32574.6932.971.518
D43079.4836.391.404
D53584.2839.801.302
Table 8:

Effect of ambient temperature Tamb

CaseTamb [°C]Tmax,core [°C]ΔTcs,max [°C]Icr,min
E12064.8231.021.614
E22565.0627.601.737
E327.1 (ref)65.1626.171.793
E43065.3124.201.876
E53565.5920.842.033
Table 9:

Effect of cement content (F-series)

CaseC [kg m−3]K [°C]Tmax,core [°C]ΔTcs,max [°C]Icr,min
F130043.758.0721.262.16
F2350 (ref)51.065.1626.171.79
F340058.372.2531.101.53
F445065.679.3536.021.33
Table 10:

Effect of surface modulus M on the thermal response (G-series)

CaseGeometryM [m−1]Tmax,core [°C]ΔTcs,max [°C]Icr,mint(Tmax,core) [h]Pcr,JCI [%]
G12 × 2 × 2 cube2.5061.5421.152.103442.8
G23 × 3 × 3 cube (ref)1.6765.1626.171.79625.6
G34 × 4 × 4 cube1.2566.5230.921.2957820.6
G45 × 5 × 5 cube1.0067.0133.671.0959437.7
G56 × 6 × 1.5 slab1.3358.8213.564.523450.1
G66 × 6 × 3 slab1.0064.4023.002.553711.2
G72 × 2 × 6 column2.1761.8224.691.743506.2
G84 × 4 × 2 pile cap1.5063.2323.411.778545.7
Figure 4:

Diagnostic M-only plots for (a) Tmax,core, (b) ΔTcs,max, and (c) Icr,min

Table 11:

Regression-model performance for geometry descriptors and crack-index response

ResponseCandidate modelPredictorsR2RMSEMAELOOCV RMSEInterpretation
Tmax,coreM-onlyM0.2272.276 °C1.623 °C2.742 °CWeak mixed-geometry descriptor
Tmax,coreQuadratic heat-flow modeldmin, dmin20.9400.633 °C0.449 °C1.108 °CStrong screening relation
ΔTcs,maxM-onlyM0.0805.494 °C4.195 °C6.824 °CNot adequate
ΔTcs,maxTwo-variable modelM, dmin0.8692.073 °C1.710 °C3.788 °CAcceptable for screening
Icr,minM-onlyM0.0031.0020.7171.220Invalid for crack-index prediction
Icr,minInverse gradient modelΔTcs,max0.9510.2230.1770.302Best performing and physically interpretable
Figure 5:

Relationship between Icr,min and ΔTcs,max for the updated G-series

Table 12:

Comparative sensitivity of peak core temperature Tmax,core and minimum crack index Icr,min

ParameterRange studiedΔTmax,core [°C]ΔIcr,minDominant response
Time step Δt0.5–8 h8.790.51Numerical stability
Element size150–500 mm0.500.39Numerical resolution
Convection coefficient hc1.5–12 kcal m−2 h−1 °C−10.960.93Icr,min
Placement temperature T015–35 °C19.140.49Tmax,core
Ambient temperature Tamb20–35 °C0.770.42Icr,min
Cement content C300–450 kg m−321.280.83Tmax,core
Surface modulus M (cubic subset G1–G4)1.00–2.50 m−15.471.008Icr,min
Geometry/surface-modulus portfolio (full G-series)1.00–2.50 m−18.193.428Shape-dependent Icr,min response
DOI: https://doi.org/10.2478/cee-2027-0013 | Journal eISSN: 2199-6512 (formerly 1336-5835) | Journal ISSN: 1336-5835
Language: English
Submitted on: May 25, 2026
Accepted on: Jun 19, 2026
Published on: Aug 18, 2026
Published by: University of Žilina
In partnership with: Paradigm Publishing Services
Publication frequency: 4 issues per year

© 2026 Tien-Toi Pham, Ngoc-Tuyen Tran, Chi-Cong Vu, Hong-Hai Tran, Hong-Ha Le, published by University of Žilina
This work is licensed under the Creative Commons Attribution 4.0 License.