Numerical Simulation and Calculation Method Study on Seamless Construction of Super-Length Raft Structures Based on Novel Magnesium Oxide Expansive Strengthening Band Method

: The drive for continuous innovation in large-scale infrastructure necessitates advancements in techniques, addressing the challenges of constructing super-length concrete structures. This study investigated the emerging shift from traditional united expanding agent (UEA) to mag-nesia expansive agent (MEA) in conjunction with expansive strengthening bands (ESBs), marking a pivotal transition in ensuring monolithic integrity. Despite a decade of exploration, MEA–ESB implementation in real-world projects remains underdocumented, with scholarly focus primarily con-ﬁ ned to material characterization. This research integrated empirical on-site tests of MEA–ESB with high-ﬁ delity numerical simulations in ABAQUS. The ﬁ nite element model (FEM) validation against actual test data underscored the precision of our modeling, capturing the complex thermomechanical behavior of the system. We introduced a sophisticated parametric analysis framework, elucidating the in ﬂ uence of critical parameters like the ESB-to-raft-width ratio and MEA concrete expansion rates. This granular understanding facilitated the ﬁ ne-tuning of design parameters, advancing the practical application of MEA methodologies. A groundbreaking contribution entailed the formulation of predictive models for early-stage cracking, anchored in the guidelines of the ACI Commi tt ee 207 and re ﬁ ned through extensive parametric exploration. These formulae empower engineers to anticipate and mitigate cracking risks during the design phase, thereby enhancing project safety and e ﬃ ciency. Notably, this study identi ﬁ ed limitations in current prediction models, highlighting the need for future research to incorporate comprehensive lifecycle considerations, including hydration heat e ﬀ ects and time-dependent mechanical property evolution.


Introduction
To meet the demands of societal production and daily life, modern architectural structures continue to evolve toward greater spans and larger volumes, resulting in the proliferation of super-length concrete structures.Traditional construction methods for such structures involve the incorporation of expansion joints or post-cast strips [1].However, expansion joints disrupt the continuity of the building structure and introduce water seepage issues that are difficult to resolve, leading to their gradual phase-out [2].The implementation of post-cast strips necessitates separate formwork support and chipping operations, adding complexity to the construction process.Moreover, post-cast strips must be cast after the main structure has been cured, thereby extending the construction duration.Additionally, engineering experience has shown that the interface between new and old concrete at post-cast strips is prone to leakage and other quality issues [1].To address these concerns, You et al. proposed an innovative seamless construction method for superlength concrete structures using the united expanding agent (UEA)-expansive strengthening band (ESB) [1,3].Following its introduction, this method has garnered extensive attention from scholars and practitioners worldwide, prompting a plethora of studies on its efficacy.
Yang et al. [4] employed the UEA-ESB method in the design and construction of a 95.8 m long beam-slab foundation for the Xinfei Sales Building, where the structure remained in excellent condition years after commissioning.Zhang et al. [5] employed a construction method in which post-cast strips were replaced with UEA-ESB during the Bank of China Jining Branch Head Office project.This substitution resulted in an expedited construction schedule and cost savings, and it upheld high-quality engineering standards throughout the execution process.Jin [6] elucidated the design philosophy and specific practices of the UEA-ESB method and pointed out its sensitivity to climate conditions, highlighting that while expansion could be maintained in humid areas, issues might arise in dry regions.Therefore, he recommended moisture retention curing and suggested cautious adoption in dry areas.Li et al. [7] stated that the UEA-ESB method can achieve the goal of self-sealing structures and control the occurrence of cracks in mass concrete structures, and they proposed relevant design and construction measures.Dong and Lin [8] designed and implemented post-cast UEA-ESB in a super-length and ultra-wide foundation structure with plan dimensions of 97.8 × 76.1 m.This resolved the problems of dry and cold shrinkage effects in super-length and ultra-wide slabs and achieved self-sealing for both the large foundation and roof structures.
Despite the widespread use of UEA in engineering practice, existing research indicates that it possesses certain limitations.For instance, the incorporation of UEA may actually increase the risk of concrete contraction cracking in projects where wet curing is challenging [9], and the products of UEA are unstable, decomposing at around 80 °C.To address these issues, Liu et al. [10] proposed the use of magnesia expansive agent (MEA) as a substitute for UEA.Studies have demonstrated that under different curing conditions, MEA effectively suppressed the drying shrinkage of concrete.Studies by Mo and coauthors investigated MgO and MgO-based expansive agents in cementitious materials for shrinkage mitigation and performance enhancement.Key findings include the following: MEA expansion properties were strongly influenced by calcination conditions, with microstructure playing a critical role [11]; quaternary blended cements with reactive MgO, slag, and fly ash effectively reduced autogenous shrinkage and exhibited improved lateage strengths [12]; homogeneous distribution of MgO additives with varying reactivities led to faster expansion and shrinkage compensation in expansive cements [13]; curing temperature significantly impacted MgO expansion behavior [14]; MgO-based additives compensated for autogenous and thermal shrinkage under non-wet conditions [15]; biochar combined with MEA synergistically mitigated shrinkage and induced expansion, while enhancing late-age compressive strength [16].Jiao et al. [17] conducted systematic research on the impact of MEA on the microstructure, mechanical properties, contraction performance, and durability of concrete.Their findings revealed that while MEA reduced the early strength and chloride ion permeability resistance of concrete, it improved its contraction characteristics and enhanced later-stage strength.Li et al. [18,19] reviewed the influence of various expansive agents, including MEA, on the performance of compensatory contraction concrete.The MEA-ESB methodology innovates by employing MEAcompensated shrinkage concrete and introducing ESBs to supplant traditional post-tensioning or construction joints, thereby enabling near-seamless continuity.The method offers several advancements over conventional practices: it simplifies construction through reduced jointing requirements, accelerating formwork cycling and overall project timelines.Moreover, it markedly improves waterproofing efficacy by limiting potential weak points, ensuring enhanced durability and resilience against water infiltration, especially vital in subterranean infrastructures.The reduction in joint treatments significantly decreases construction time and costs, while the promotion of a monolithic structure bolsters structural integrity and mitigates crack formation attributed to environmental factors.Furthermore, this technique's adaptability to intricate designs fosters greater architectural creativity without compromising structural dependability, positioning the MEA-ESB method as a highly advantageous strategy within contemporary construction practices.
In the field of thermal crack control in the seamless construction of super-long concrete structures, Klemczak, Zmij [20], and Jędrzejewskaa et al. [21] delved into the reliability and application of existing guidelines for designing and evaluating early-age cracking in concrete structures, focusing particularly on tank walls, bridge abutments, and wallon-slab structures.These studies compared international guidelines (American [22], European [23,24], and Japanese [25]) and practical standards (like EN 1991-4 [26] and CIRIA C766 [27]), identifying discrepancies between predicted and actual crack widths.The need for refined predictive models was underscored, with the literature [21] proposing adjustments to the calculation methods to better anticipate cracking behavior under restrained hardening conditions.Smolana et al. [28] and Klemczak et al. [29] emphasized the importance of realistic estimations of early-age thermal loads and induced stresses in mass foundation slabs and reinforced concrete walls.They reviewed analytical and numerical methods for assessing cracking risks, highlighting the significance of considering the exothermic reactions of cement hydration and subsequent temperature variations.Notably, the literature [28] presented a case study involving real-world application of these methods alongside on-site measurements, underlining the practical implications and limitations of each approach.Klemczak, Knoppik-Wróbel [30], and Li et al. [31] presented case studies, the former detailing the construction of a mass concrete tunnel in Shantou City, China, and the latter discussing the behavior of tank walls and bridge abutments.Both highlighted the importance of understanding thermal-shrinkage mechanisms and showcased successful mitigation techniques, such as gradient concrete and plastic anti-crack grids.The superiority of gradient concrete in controlling cracking while reducing temperatures was emphasized, suggesting its broader implementation.Kanavaris et al. [32] and Azenha et al. [33] contributed to the understanding of thermal and shrinkage stress development by offering comprehensive analytical models and expert recommendations.The literature [32], as part of RILEM TC 254-CMS's work, improved the concept of 'massivity' to predict thermal cracking potential, incorporating binder type, content, and casting conditions.Meanwhile, the literature [33] provided recommendations from RILEM TC 287-CCS focusing on the simulation of thermo-chemo-mechanical behavior in massive concrete structures, covering aspects from material properties to numerical modeling peculiarities.
In summary, while the ESB technology has existed for over three decades, the adoption of MEA as a substitute for UEA is a relatively recent development, with merely a decade of exploration.Past research predominantly focused on material-level investigations, with limited documentation of real-world implementations.This study broke new ground by integrating on-site MEA-ESB test data with numerical simulations, affirming the precision of our Finite Element Model constructed in ABAQUS.By introducing advanced parametric analysis and devising calculation formulas for predicting cracking in super-long raft structures with ESBs, our work transcended conventional material-centric studies.The significance of our approach lies in three advancements.Firstly, our simulation-driven exploration of MEA-ESB dynamics, validated against empirical evidence, illuminated the underlying mechanics governing these systems.It bridged the gap between theoretical constructs and practical engineering manifestations, reinforcing the foundation for future design and construction practices.Secondly, the creation of a parametric analysis framework offered a quantitative lens to decipher the impacts of critical variables, including the ESB-to-raft-width ratio and expansion rates of MEA concrete, empowering engineers with a toolset for fine-tuning and enhancing the MEA methodology.Lastly, the introduction of formulaic predictions for early-stage cracking, informed by the framework of the ACI Committee 207 and findings from parametric analysis, streamlined the comprehension of intricate parameter interactions.It empowered proactive decision-making during the design stage by enabling the anticipation of structural crack, thereby mitigating construction hazards and safeguarding overall project integrity.Our findings thus chart a progressive path for the broader adoption and optimization of the MEA-ESB strategy in large-scale engineering projects.

Model Establishment
A finite element model (FEM) of the raft structure was constructed using the generalpurpose finite element software ABAQUS 2022.The implicit analysis was run in the simulations.User subroutines were developed in Fortran language for secondary development.The heat generation in concrete was realized through the HETVAL user subroutine interface.Concrete shrinkage and expansion of the strengthening bands were achieved using the UEXPAN user subroutine interface.Lastly, the change in concrete's elastic modulus over time was implemented via the UMAT user subroutine.

Material Constitutive Relations
The standard region of the raft was constructed with C35 concrete, and the ESB employed C40 concrete; detailed mix proportions are outlined in Reference [34].The testing procedures for assessing the physical and mechanical properties of the concrete, not detailed in [34], adhered to the requirements outlined in the 'Standard for test methods of concrete physical and mechanical properties' (GB/T50081-2019) [35], ensuring a rigorous evaluation of the materials' performance.According to the 'Test code for hydraulic concrete' SL/T 352-2020 [36], the thermodynamic properties of both C35 and C40 concrete were calculated separately, as presented in Table 1.Additionally, the 28-day cube compressive strengths of the concrete are also listed in Table 1.The calculation of contraction for concrete without MEA adopted the contraction model stipulated by CEB-FIP (1990) [37].The variation in the elastic modulus of C35 and C40 concrete with age was computed according to the specifications of 'Standard for Construction of Mass Concrete' [38].The heat source function for concrete was determined using the function model specified in the 'Test code for hydraulic concrete' SL/T 352-2020 [36], with the relevant thermodynamic parameters calculated based on the concrete mix proportion.
Experiments were conducted on the concrete with MEA to obtain their strain-time curves, as shown in Figure 1.As depicted in the figure, in the absence of expansion agents, noticeable shrinkage of the concrete commenced on the second day post-pouring, with shrinkage strain escalating progressively over time, culminating in a maximum value of 411 µε at 56 days.Conversely, concrete admixed with CaO as an expansion agent initially exhibited expansion, peaking at a maximum expansion strain of 83 µε on the eighth day, after which shrinkage set in as the age increased, leading to a final maximum shrinkage strain of 120 µε at 56 days.Compared with CaO-injected concrete, concrete incorporating MgO as the expansion agent sustained expansion until the 17th day, achieving a maximum expansion strain of 374 µε, which was 4.5 times that of the CaO counterpart.Thereafter, the concrete underwent marginal contraction; nonetheless, at 56 days, it still registered an expansion strain of 298 µε, indicative of a sustained expansion effect.

Element Type
The raft and the foundation were simulated using the eight-node thermally coupled hexahedral element C3D8T.A sensitivity analysis was performed on the mesh to establish an appropriate grid size of 1000 mm.The raft was divided into 10 layers in the thickness direction, resulting in a total of 44,427 elements.Meanwhile, the foundation contained 85,791 elements.The established FEM is illustrated in Figure 2.

Boundary Conditions
The bottom of the foundation was subjected to consolidation boundary conditions.Equally important are the thermal boundary conditions for the structure.An equivalent coefficient of heat dissipation was introduced, treating the insulation layer in contact with air as a third-type boundary condition.The heat dissipation coefficient was calculated in accordance with the formula specified in the 'Test code for hydraulic concrete' SL/T 352-2020 [36], based on the thickness of the insulation layer.The ambient air temperature was taken as the average value of 25 °C from on-site measurements, and the mold temperature of the concrete was also determined based on on-site measurements.In the initial step of the analysis, the raft structure and foundation were subjected to an initial temperature by defining the temperature field variables.Subsequently, the temperatures of the raft and foundation were computed using the HETVAL user subroutine.To facilitate comparison with experimental results, a thermal analysis step was used with a fixed time step of 1 h.

Numerical Simulation Results and Validation
This section aims to rigorously validate the predictive capability of the established FEM for the mechanical and thermodynamic performance of super-length raft structures with ESBs through systematic comparative analyses.Specifically, this section presents in detail the outcomes from finite element analysis (FEA), including the principle of temperature distribution, principle of strain distribution, temperature-time characteristics, and strain-time characteristics, juxtaposing them directly with data obtained from on-site tests [34].The validation process encompassed not only quantitative comparisons of results, such as the congruence between modeled temperatures and strains with actual measurements, but also qualitative assessments of characteristic features to ensure that the simulation accurately replicates physical phenomena observed in the on-site tests.Furthermore, any notable discrepancies were meticulously examined to identify their potential causes and evaluate their implications on the overall efficacy of the model.Through this comprehensive validation exercise, the objective was to ascertain the accuracy and validity of the constructed FEM, thereby laying a robust foundation for subsequent parametric studies and practical engineering applications.
While the fundamental parameters of the model were rigorously calculated based on material properties, we acknowledged that any theoretical model inevitably involves a degree of approximation.Consequently, a parameter sensitivity analysis was conducted to investigate the impacts of minor variations in key parameters, such as the maximum adiabatic temperature rise of concrete and the reaction rate coefficient, on the model outputs.This process aided in understanding the model's responsiveness to input variations, not by directly 'tuning' parameters to fit on-site test data, but rather through scientific analysis to delineate the optimal range of parameter estimates.In light of these analyses, when systematic discrepancies emerged between model predictions and test findings, rigorous scientific methodologies were employed to scrutinize the validity of model assumptions, boundary conditions, and computational methodologies.For instance, enhancements to numerical simulation parameters, like optimizing grid density and time steps, were pursued to minimize numerical errors rather than blindly manipulating parameters to achieve data concordance.

Principle of Temperature Distribution
According to the literature [34], during the thermal accumulation phase, the temperature of the raft continuously rose, reaching its peak at 48 h post-pouring.The principle of temperature distribution across the surface, center, and bottom of the raft, as compared between FEA and experimental data, are depicted in Figure 3.The figure also presents the temperature cloud maps obtained through numerical simulation.It is apparent from the figure that the principle of surface temperature distribution simulated by the finite element method deviated most significantly from the experimental results, with a maximum discrepancy of 11.2 °C.This is attributable to the fact that surface temperatures are heavily influenced by environmental factors such as ambient temperature, wind speed, sunlight, and curing conditions-influences that finite element methods cannot accurately simulate.Nevertheless, the finite element method still managed to reflect the trend of temperature distribution across the surface of the raft.The simulated temperature distributions at the center and bottom of the raft aligned well with the experimental findings.The finite element approach accurately captured the localized temperature peaks of the central temperature within the ESB and accurately reflected the high-temperature plateaus observed at the bottom of the raft in areas with dense pits or wells.During the thermal balance phase, the temperature of the raft tended to stabilize.As illustrated in Figure 4, the FEA closely mirrored the observed trends in temperature distribution across the raft's surface, center, and bottom after 28 days of curing.The FEA accurately captured the characteristic pattern of lower temperatures at both ends and a peak in the center, which aligns well with the empirical data obtained from the experiments.

Principle of Strain Distribution
Upon completion of 48 h post-casting, the distribution principle of strain across the surface, center, and bottom of the raft, as depicted in Figure 5, were analyzed through FEA.The FEA results accurately captured the characteristic of compressive strains being minimal at the ends and maximal at the center.Moreover, the FEA was able to reflect the unique distribution principle of increased strains on the surface within the lowering plate region.The results obtained from the FEA demonstrated a commendable congruence with the data collected from the experiments.After casting for a duration of 28 days, the FEA results concerning the strain distribution principle across the surface, center, and bottom of the raft were compared with experimental data, as illustrated in Figure 6.The FEA demonstrated a high level of accuracy in simulating the strain distribution characteristics, including the special trend where the surface strain increases in the central region of the ESB.The congruence between the finite element predictions and the experimental findings was satisfactory.

Temperature-Time Characteristic
The temperature-time curves of the ESB 1 are presented in Figure 7, comparing the FEA results with those from the experimental tests.The correlation between the FEA outcomes and the data derived from experimental trials exhibits a noteworthy degree of alignment.The maximum temperatures predicted by the FEA at the surface, center, and bottom of the ESB were 60.8 °C, 76.6 °C, and 55.8 °C, respectively, showing deviations of 0.7%, 0.8%, and 0% when compared with the experimental values of 61.2 °C, 77.2 °C, and 55.8 °C.Furthermore, the FEA accurately simulated the temporal evolution of temperatures within the ESB.As shown in Figure 8, the FEA and experimental data exhibited a commendable degree of concurrence regarding the temperature-time curves for standard region 2. The FEA predicted maximum surface, center, and bottom temperatures of 58.8 °C, 74.4 °C, and 74.2 °C, respectively.These values reflect deviations of 5.2%, 0.5%, and 1.7% from the corresponding experimental measurements of 62 °C, 74.8 °C, and 75.5 °C.Consequently, the FEA effectively captured the temporal progression of temperatures within the standard region, underscoring its reliability in such simulations.

Strain-Time Characteristic
The strain-time curves for the ESB 1 are illustrated in Figure 9, where the FEA demonstrates a satisfactory alignment with the experimental data.The FEA simulations determined the maximum compressive strains at the surface, center, and bottom of the band to be −171.7 µε, −373.2 µε, and −191.9 µε, respectively.These values exhibit discrepancies of 23.5%, 6.1%, and 0% when compared with the experimental findings of −131.4 µε, −351.6 µε, and −191.9 µε.This indicates that the FEA was capable of accurately replicating the temporal evolution of strains within the ESB.The strain-time curves for standard region 2, as depicted in Figure 10, exhibit a commendable correlation between the FEA and the experimental data.The FEA simulations yielded maximum compressive strains of −72.6 µε at the surface, −334.6 µε at the center, and −168.4 µε at the bottom.These values reflect discrepancies of 10.3%, 9.4%, and 7.5%, respectively, when compared with the experimental results of −80.9 µε, −306.0 µε, and −182.2 µε.Consequently, the FEA was deemed to accurately simulate the temporal development of strains within the standard region.

Parameter Analysis Model
According to ACI Committee 207 [22], cracking in raft structures initiates when the tensile stress σt(t) within the raft equals the tensile strength ft(t) of the concrete.The criterion they propose for predicting early-age cracking in the raft structures is as follows: ( The ACI Committee 207 [22] recommends applying the aforementioned criterion at 7 days after pouring and suggests that the maximum tensile stress σt in the raft concrete can be calculated using the following formula: ( The meanings and calculation methods of the parameters in this equation will be elaborated upon in Section 4. This formula pertains to ordinary super-length concrete structures without considering the influence of ESBs.The ESBs impose a pre-compressive strain on the raft concrete, which in turn reduces the tensile stress experienced by the concrete, as illustrated in Figure 11.Consequently, taking the impact of ESBs into account, the adjusted formula for calculating the maximum tensile stress σ(t) in the raft concrete becomes: ( Figure 11.The mechanical analysis model for the raft structure.
The paramount design parameters for ESBs are the ESB width and the MEA concrete expansion rate [1,3].To facilitate dimensionless analysis, the ESB width was normalized, and the ESB-to-raft-width ratio and the MEA concrete expansion rate were adopted as the analytical parameters.An FEM of the basic load-bearing elements was established using ABAQUS, with relevant parameter values consistent with those detailed in Section 2.1, as shown in Figure 12.The dimensions of the standard region of the raft structure were specified as 20 × 20 × 2 m, with the ESB width set at 2 m and the MEA concrete expansion rate at 0.04%.

The ESB-to-Raft-Width Ratio r
The range for the ESB-to-raft-width ratio was selected to be between 0 and 0.2.When the ESB-to-raft-width ratio was 0, it indicated that no ESB had been implemented.The strain-time curves for the surface of the raft at varying ESB-to-raft-width ratios are depicted in Figure 13a.It is observed that as the ESB-to-raft-width ratio increased, the surface strain of the raft became less pronounced at the same time point.Figure 13b illustrates that the 7-day strain decreased linearly from −307 µε to −411 µε as the ESB-to-raft-width ratio rose from 0 to 0.2, representing a 34.1% decrease.A broader ESB width resulted in greater pre-compressive stress exerted upon the raft, thereby effectively reducing the surface strain [39].

The MEA Concrete Expansion Rate e
Pursuant to the provisions stipulated in the current national standard 'Code for Concrete Admixture Application' [40], the value for limiting the expansion rate of ESB was generally categorized at a level of 0.005%, with a minimum limit of 0.015%.Concrete with an expansion rate exceeding 0.06% is classified as self-stress concrete [40].In the present study, the parameter range for the MEA concrete expansion rate was set between 0% and 0.1%, with each level being 0.02%.The conditions where the MEA concrete expansion rate was 0 represented scenarios without the installation of ESB.The strain-time curves on the surface of the raft at varying MEA concrete expansion rates are presented in Figure 14a, indicating that as the MEA concrete expansion rate increased, the surface strain of the raft decreased at the same time point.Figure 14b reveals that when the MEA concrete expansion rate increased from 0% to 0.1%, the 7-day strain decreased linearly from −307 µε to −404 µε, a reduction of 31.7%.

Computational Formula
Section 3.2 presented selected outcomes from the parametric analysis, which encompassed numerical simulations of a total of 31 raft structures.Numerical simulation results under varying parameter conditions are presented in Table 2. ) The MEA concrete expansion rate (%) Compared with the raft structures without ESBs, those with ESBs exhibited reduced 7-day strain, as the ESB imposed a pre-compressive strain on the raft structure.The differential 7-day strain between these two conditions is denoted by εESB.The response surface methodology (RSM) is a statistical technique capable of establishing an implicit functional relationship between research factors and response values through a judiciously selected finite number of simulations [41].In this study, RSM was employed to determine the correlation between the 7-day strain of the raft and the parametric variables.
The calculation formula for the εESB, as obtained using the RSM, is presented as follows: where εESB is the pre-compressive strain imposed on the raft concrete by the ESB; r is the ESB-to-raft-width ratio; and e is the MEA concrete expansion rate.The εESB of the raft calculated using Equation (4) was compared with the results from the FEA, as depicted in Figure 15.Inspection of the data revealed a high level of concordance between the strains computed using the specified formula and those determined through FEA, as evidenced by an adjusted R 2 value of 0.946, indicating a strong predictive model fit.By substituting Equation (4) into Equation (3), the resultant formula for calculating the maximum tensile stress in raft structures with ESBs was obtained: ( where σ(t) is the maximum tensile stress in the raft concrete with ESBs 7 days after pouring.; KR is the coefficient pertaining to the vertical distribution of restraint intensity along the element's height; Kf is the factor quantifying the extent of restraint implementation; αT is the thermal expansion factor, µε/°C; ΔT is the inner-surface temperature difference of the raft, °C; and Ecm,eff(t) is the enduring Young's modulus associated with concrete 7 days after casting.For raft structures, the ratio of length L to height H is generally greater than 2.5.According to ACI Committee 207 [22], KR = [(L/H − 2)/(L/H + 1)] h/H , where h is the vertical distance from the plane of the restraining element.
The factor Kf is calculated using the following formula: where Ac is the sectional area under consideration for the concrete component, m 2 ; Ec is the Young's modulus of the concrete, MPa; Ar is the restraining element's cross-sectional area that imposes confinement on the subject component, m 2 ; and Er is the Young's modulus of the restraining element, MPa.The inner-surface temperature difference ΔT is computed using the one-dimensional finite difference method, as prescribed in reference [38], and is therefore not elaborated upon further in this paper.
The enduring Young's modulus Ecm,eff(t) associated with concrete 7 days after casting is calculated using the following formula: where β is the correction coefficient for fly ash and slag content, according to the table in reference [38] and E0 is the elastic modulus of concrete at 28 days under standard curing conditions.Following these calculations, the maximum tensile stress in the raft structures with ESBs was obtained, which was then compared against the tensile strength of the concrete.Employing the cracking criterion outlined in Equation (1), a determination was made regarding the potential for cracking in the raft structures.This process provided a valuable reference for the design and construction of super-length raft structures with ESBs.
This study, grounded in the specifications of ACI Committee 207 [22], established crack control criteria for super-length raft structures with ESBs and proposed a formula to calculate the maximum tensile stress 7 days after concrete pouring.Nonetheless, the methodology put forth herein is limited to utilizing the maximum tensile stress at 7 days post-pouring as the sole criterion for cracking assessment, and does not extend to determinations at other time intervals.Additionally, the study assumes the age-dependent mechanical behavior of concrete follows the general trends prescribed by standards, overlooking the potential acceleration in the concrete's mechanical property development due to the influence of hydration heat.

Conclusions
This study utilized on-site test data based on the MEA-ESB to construct an FEM of a raft structure using the finite element software ABAQUS.Comparison with test data revealed that the FEM accurately simulated the spatiotemporal distribution and variation principles of temperature and strain within the raft.A parametric analysis model was established, and the cracking calculation formula for the super-long raft structure with ESBs was studied.The key conclusions drawn from this research are as follows: (1) An FEM of the super-length raft structure with ESBs was developed using the ABAQUS software.The HEATVAL, UEXPAN, and UMAT user subroutines were utilized to simulate the thermal hydration, autogenous shrinkage and expansion behavior, and development of the elastic modulus of concrete, respectively.This approach effectively replicated the spatial and temporal distribution principles of temperature and strain within the raft, thereby demonstrating the viability of this modeling methodology.
(2) Based on the framework of the ACI Committee 207 and findings from parametric analysis, a crack determination formula for super-long raft structures with ESBs was proposed.The employment of formula-based methods circumvents the substantial workload associated with finite element approaches, enabling designers and engineers to make relatively accurate assessments of the cracking status in super-length raft structures with ESBs.
(3) The proposed formula-based method was unable to predict cracking for all times following the pouring of super-long raft structures with ESBs and did not account for the effects of hydration heat on the development of the concrete's mechanical properties.Future research necessitates the development of a more comprehensive mathematical model, one which simultaneously incorporates the evolution of hydration heat in concrete, temperature field distribution, variations in stress and strain, as well as shrinkage effects, and integrates a time variable into the model.This enhancement will enable the model to dynamically reflect the performance changes in concrete throughout its entire lifecycle.

Figure 1 .
Figure 1.The impact of different expansive agents on unrestrained volume deformation of concrete.

Figure 2 .
Figure 2. The finite element model of raft structure and foundation.

Figure 12 .
Figure 12.The parametric analysis model for the raft structure.

Figure 13 .
Figure 13.The impact of the ESB-to-raft-width ratio on the strain-time curves and the 7-day strain.(a) Strain-time curves; (b) 7-day strain.

Figure 14 .
Figure 14.The impact of the MEA concrete expansion rate on the strain-time curves and the 7-day strain.(a) Strain-time curves; (b) 7-day strain.

Figure 15 .
Figure 15.Comparison of calculation results from the formula and FEA.

Table 1 .
The mechanical and thermodynamic properties of concrete.

Table 2 .
Parameter analysis computational results.