Abstract
This paper addresses the issue of crack expansion in adjacent buildings caused by foundation pit construction and develops a predictive model using the response surface method. Nine factors, including the distance between the foundation pit and the building, soil elastic modulus, and density, were selected as independent variables, with the crack propagation area as the dependent variable. An orthogonal test of 32 conditions was conducted, and crack propagation was analyzed using the FEM-XFEM model. Results indicate that soil elastic modulus, Poisson’s ratio, and distance between the pit and building significantly impact crack propagation. A predictive model was developed through ridge regression and validated with additional test conditions. Single-factor analysis showed that elastic modulus and Poisson’s ratio of the silty clay layer, elastic modulus of sandy soil, and pit distance have near-linear effects on crack propagation. In contrast, cohesion, density, and Poisson’s ratio of sandy soil exhibited extremum points, with certain factors showing high sensitivity in specific ranges. This study provides theoretical guidance for mitigating crack propagation in adjacent buildings during excavation.
1 Introduction
The rapid pace of urbanization has significantly increased the demand for underground space utilization in cities, which is essential for accommodating expanding populations and infrastructure needs (; ; ; ; ). However, this rapid development often takes place in densely packed urban environments, where construction of foundation pits is frequently conducted in close proximity to existing structures. This proximity poses a unique challenge, as excavation activities can induce ground movements that impact surrounding buildings, some of which may already exhibit pre-existing damage (; ; Wang et al., 2023a). The consequences of these disturbances are critical, as they can not only exacerbate existing structural damage but also lead to catastrophic failures in the most extreme cases. This issue is particularly critical in urban environments, where the interaction between foundation pit excavation activities and aging infrastructure demands in-depth studies to evaluate its impact on structural safety, crack propagation, and effective mitigation strategies (Xue, 2023).
Extensive research has been conducted on the impact of foundation pit excavation on neighboring structures, particularly focusing on the resulting ground settlement and structural deformations. Early studies, such as those by Peck, using a wealth of field monitoring data, proposed a method to estimate surface settlement that accounts for soil parameters and excavation depth (Wei et al., 2021). Similarly, Hsieh divided the settlement profile into a shoulder type and a concave type and noted that the maximum settlement occurs within a distance of 2H (H being the excavation depth of the pit) (). Research by Wang et al. demonstrated that as excavation depth increases, the lateral displacement and deformation of retaining structures increase, which in turn exacerbates the settlement of adjacent buildings (Wang et al., 2021). While these studies provide valuable insights into ground displacement, they tend to oversimplify the problem by not fully considering the interaction between the soil, the pit, and the surrounding buildings. The presence of nearby building structures can significantly influence soil behavior, including soil loss and surface settlement. Specifically, the structural integrity of neighboring buildings can alter local soil stress distributions and pore pressure gradients, which in turn affect the soil’s compaction and settlement patterns. The impact of these changes is particularly critical in urban environments, where even minor disturbances to soil can lead to substantial damage to adjacent structures.
Ding et al. utilized the finite difference method to simulate the entire excavation process, analysing the effects of isolation pile construction, diaphragm wall construction, and foundation pit dewatering and excavation on adjacent bridge piers (). Lan conducted a numerical analysis of the stress characteristics of buildings under differential settlement and established threshold values for differential settlement (Tirca, 2017). Ou simulated the damage behavior of historical buildings affected by differential settlement via a combined finite-discrete element method (FDEM) (). Zhang, through field monitoring and numerical simulation, analysed the factors influencing the shear deformation of buildings induced by foundation pit construction (Zhang, 2023). Chen adopted numerical simulations to study the distribution patterns of surface settlement caused by excavation and the influence of factors such as diaphragm wall thickness and depth (Zhao et al., 2023). Lu et al., using a deep foundation pit project as a case study, applied the finite element method to identify key factors affecting pit deformation and building settlement (). Zhong, using the Baiyun Station foundation pit project in Guangzhou as a case study, developed a three-dimensional numerical model of the foundation pit, elevated bridge, and pile-beam foundation to analyse the effects of excavation on the lateral displacement of adjacent piles (Zhong et al., 2023). While these studies have taken into account the coupled system of the soil-pit-building structure, they primarily focus on analyzing factors such as stress, strain, and displacement in existing building structures, overlooking issues related to their service condition and crack propagation. In reality, cracks in building structures can significantly affect their load-bearing capacity and safety performance. The presence of cracks may lead to local stress concentration, thereby reducing the overall bearing capacity of the structure. This is because the expansion of cracks will destroy the continuity of the material, resulting in its inability to effectively share the external load under load. Therefore, in order to enhance the resistance of structures, new materials and structures are used in the construction of buildings and infrastructure (; Wang et al., 2023b).
Xu et al. studied the impact of soil parameters on the arching effect in deep foundation pits, emphasizing the key role of cohesion and internal friction angle in shaping the arching characteristics (Xu et al., 2024). Ye et al., through the study of actual engineering cases, confirmed that the impact of foundation pit excavation on adjacent tunnel structures is closely related to the distance from the pit and the geological conditions (Ye et al., 2021). Based on monitoring data from a foundation pit project, Yan et al. investigated the effect of foundation pit excavation on the cracking of nearby bridge pier structures. The results showed that the crack distribution is closely related to the proximity and excavation depth (Yan et al., 2024). Numerous studies have shown that soil parameters such as elasticity, cohesion, and internal friction angle, as well as the distance between the foundation pit and the building, play a crucial role in determining the magnitude of ground settlement and its effects on the surrounding structures (; ; ; ). In particular, factors like the stiffness and cohesion of the soil, and the distance from the excavation, significantly affect the settlement response and the resulting damage to adjacent buildings. While soil types such as silty clay and sandy soils are commonly considered in these analyses, the variability in soil properties can lead to differences in settlement behavior. Therefore, the factors selected for analysis in this study focus on the most influential parameters for ground settlement and building deformation. These parameters include the elastic modulus, Poisson’s ratio, and density of the silty clay and sandy soil layers, as well as the distance between the foundation pit and the building.
A review of the above studies reveals a significant gap in current research regarding the pre-existing damage conditions of adjacent buildings. Structures with pre-existing cracks or other forms of damage are more susceptible to further deterioration under the dynamic loading conditions induced by foundation pit excavation. This makes it essential to integrate the structural health of neighboring buildings into the predictive models. Failure to accurately represent and fully consider the damage status of existing buildings can lead to an underestimation of the impact of foundation excavation, jeopardizing the safety of these structures. This study investigates the effect of foundation pit excavation on the crack propagation of adjacent buildings, with a focus on buildings exhibiting pre-existing damage. A combination of finite element method (FEM), extended finite element method (XFEM), and response surface method (RSM) is employed in this study to model and predict the impact of excavation-induced disturbances on building integrity. An analysis of the disturbances induced by the excavation of a new foundation pit is conducted, and the patterns of damage propagation in nearby existing structures are discussed. This study aims to provide guidance for the rational design and construction of foundation pits.
2 Theoretical basis and methodology
2.1 FEM-XFEM methodology
In damage theory, Kachanov regarded microcracks and microvoids in materials as defects that do not bear load (). Consequently, the effective load-bearing area of the material decreases from to . If we assume that these microdefects are uniformly distributed in all spatial directions, and that is independent of the normal direction, an isotropic damage variable can be defined as in Equation 1.where is the isotropic damage variable, is the nominal load-bearing area, and is the effective load-bearing area.
When studying structural damage in buildings, both the length and depth of cracks influence the load-bearing capacity and operational performance of the structure (; ; ). Existing engineering practices and theoretical research indicate that cracks induced by disturbances from foundation pit excavation typically appear at the junctions of walls and top slabs. Since these components are relatively thin, the resulting macrocracks can be approximated as planes perpendicular to the direction of loading (). The area of these cracks can be approximated as the product of the crack length and depth. In this study, the cracked portion of the building is treated as an area that cannot bear tensile forces. Thus, the larger the crack area, the smaller the effective load-bearing area and the greater the degree of structural damage. Conversely, the smaller the crack area, the larger the effective load-bearing area and the lesser the degree of structural damage. Therefore, the crack propagation area c is used to represent the extent of damage in the building, and it is defined as in Equation 2.where is the crack propagation area, is the crack length, and is the average crack depth.
In the XFEM model, a crack plane with a specified length and width can be set to be inserted at a specified position in the model as pre-damage of the structure. During XFEM analysis, the mesh that intersects the crack plane is considered to be cracked and cannot withstand the tensile force of the crack plane normal. In this way, the Kachanov damage theory can be applied to XFEM analysis.
A 3D finite element model that integrates the soil, foundation pit, and building was constructed via advanced finite element software (Figure 1). The model dimensions are 400 m × 150 m × 50 m, and a total of 60,384 elements are generated. The soil, diaphragm walls, and building are all modelled with solid elements. The contact between the soil and diaphragm walls, as well as between the building foundation and the soil, is set as surface-to-surface contact, with tangential behavior defined as “penalty friction” and normal behavior as “hard contact.” The diaphragm walls and internal supports are modelled with coupled contact conditions. The soil is represented via the Mohr‒Coulomb material model, whereas the internal supports and building structures are modelled via a linear elastic material model (; Zhou et al., 2024; ). Boundary conditions constrain displacements in all three directions at the bottom of the soil and normal displacements along the sides. Dewatering of the foundation pit is simulated by setting pore pressure boundary conditions, and isotropic seepage models are used for groundwater flow (; ). The excavation of the foundation pit is simulated via the *Model Change feature, with excavation steps set according to the conditions listed in Table 1. The specific parameters of the FEM model are shown in Table 2.
FIGURE 1
TABLE 1
| Analysis step | Working conditions |
|---|---|
| Geo-stress balance | Initial stress field and seepage field balance |
| Building | The application of the building |
| Dewatering 1 | First dewatering (−4.5 m) |
| Excavation 1 | Excavating the first layer of soil (−2.5 m) |
| Dewatering 2 | Second dewatering (−10 m) |
| Excavation 2 | Excavating the second layer of soil (−8 m) |
| Dewatering 3 | Third dewatering (−15.5 m) |
| Excavation 3 | Excavating the third layer of soil (−13.7 m) |
| Dewatering 4 | Fourth dewatering (−18.7 m) |
| Excavation 4 | Excavating the fourth layer of soil (−16.7 m) |
Analysis step settings (Yang et al., 2024).
TABLE 2
| Type | Items | Density (kN/m3) | Cohesion (kPa) | Internal friction angle (°) | Elastic modulus (MPa) | Poisson’s ratio | Void ratio |
|---|---|---|---|---|---|---|---|
| Soil | Silty clay | 20 | 34.02 | 18.66 | 9.03 | 0.33 | 0.71 |
| Sand | 19.7 | 0 | 31 | 30 | 0.3 | 0.63 | |
| Siltstone | 23.5 | 400 | 40 | 400 | 0.33 | — | |
| Building | Wall, plates, Beam, Foundation | 24 | — | — | 2.25 × 104 | 0.2 | — |
| Inner support | Diaphragm Wall | 24.5 | — | — | 3.15 × 104 | 0.2 | — |
| Concrete support | 24.5 | — | — | 3.15 × 104 | 0.2 | — | |
| Steel support | 78.5 | — | — | 2.1 × 106 | 0.25 |
FEM model parameters (Yang et al., 2024).
Since conventional finite element methods cannot simulate strong discontinuities such as cracks, the extended finite element method (XFEM) is used to simulate crack propagation in buildings under settlement (Yang et al., 2024; ). Referring to the two-stage analysis method in the literature (; ; ; ; ), an XFEM model containing only the building was established (Figure 2). The settlement calculated by the FEM model at the building’s foundation is applied as displacement vectors to the corresponding positions in the XFEM model (Figure 3). The cohesive crack initiation criterion is used in the model, and on the basis of reference (; ), the maximum principal stress is set to 1.7 MPa and the fracture energy to 150 N/m. The specific parameters are shown in Table 3.
FIGURE 2
FIGURE 3
TABLE 3
| (MPa) | (kN/m3) | (MPa) | (N/m) | |
|---|---|---|---|---|
| 22,500 | 24 | 0.2 | 1.7 | 150 |
XFEM model parameter.
Note: , , , , are the Elastic modulus, density, Poisson’s ratio, Maximum principal stress, Fracture Energy of the building structure.
While XFEM is an effective tool for modeling crack propagation in building structures under foundation pit excavation, it has certain limitations. For example, XFEM assumes that cracks propagate according to predefined failure criteria (such as stress intensity factors or strain criteria). This may not fully capture complex crack behaviors in real-world materials, where cracks can sometimes exhibit more irregular patterns or be influenced by additional factors like soil-structure interaction or dynamic loading. Additionally, XFEM requires a sufficiently refined mesh to accurately capture crack initiation and propagation, which can be computationally expensive. These limitations should be considered when applying XFEM in real-world scenarios, and further research may be needed to improve its applicability to a wider range of structural conditions.
2.2 Response surface method
In actual construction, different working conditions involve varying construction parameters, meaning that there are countless scenarios. If a numerical model were to be built for each scenario, it would undoubtedly increase computational costs. Therefore, this paper establishes a response surface model to predict the computational results under different construction parameters. The response surface method (RSM) is a statistical approach based on scientifically reasonable experimental designs, which quantifies variables and, through multiple trials, produces a simple mathematical model to describe the unknown relationships between these variables.
Given the high stiffness of the argillaceous siltstone layer, its variability has a relatively small effect on ground settlement and is therefore excluded from the parameter analysis. Additionally, since the plasticity parameters of sandy soil are generally fixed and show minimal variability, they are also excluded. Referring to the parameter sensitivity analysis by Lu (), this study selects the following influencing factors: the elastic modulus , density , cohesion , internal friction angle , and Poisson’s ratio of the silty clay layer, as well as the elastic modulus , density , and Poisson’s ratio of the sandy soil layer. Furthermore, the distance between the foundation pit and the building is considered, resulting in a total of nine influencing factors, and the dependent variable is the crack propagation area c in the building. The first-order model formula is as follows Equation 3.where represents the i-th influencing factor (), denotes the coefficient of , and represents the error term.
The third-order model formula includes square and cubic terms, which allows consideration of the strong nonlinear influence of the factors on the predicted variable. The formula is as follows Equation 4.
In this equation, the nine independent variables lead to 28 coefficients, which requires 28 equations to solve. To eliminate the effects of dimensional differences between parameters and ensure the comparability of the data, the orthogonal test results were normalized via the Equation 5 (; ):where represents the normalized data, represents the original data, represents the maximum value, and represents the minimum value.
When the response surface function is fit via the least squares method, high-order functions can easily lead to overfitting. Therefore, ridge regression was used to fit the re-sponse surface function in this study. Ridge regression, a regularized linear regression model, helps mitigate multicollinearity and overfitting issues. The penalty parameter for ridge regression was set to 1.
Since the number of influencing factors in this study is large, analysing the effect of each factor on crack propagation individually would require many experiments. To optimize computational resources, the orthogonal experimental method is used to investigate the impact of various factors on crack propagation. In orthogonal experiments, too few levels for each factor can affect the prediction accuracy, whereas too many levels can result in excessive computational demands. On the basis of the coefficient of variation for each parameter, four levels were selected for each factor. Ultimately, the L32(49) orthogonal experiment table was chosen, and the resulting orthogonal table is shown in Table 4. In this table, Group 1 corresponds to the average values of all the parameters.
TABLE 4
| Test groups | (MPa) | (MPa) | (kN/m3) | (kN/m3) | (kPa) | (°) | (m) | ||
|---|---|---|---|---|---|---|---|---|---|
| 1 | 9.03 | 30 | 20 | 19.70 | 34.02 | 18.66 | 0.33 | 0.3 | 13.8 |
| 2 | 9.03 | 12 | 18.8 | 18.52 | 13.61 | 16.79 | 0.27 | 0.24 | 3 |
| 3 | 9.03 | 21 | 19.4 | 19.11 | 23.81 | 20.53 | 0.3 | 0.27 | 8 |
| 4 | 9.03 | 39 | 20.6 | 20.29 | 44.23 | 22.39 | 0.36 | 0.33 | 20 |
| 5 | 6.24 | 30 | 20 | 18.52 | 13.61 | 20.53 | 0.3 | 0.33 | 20 |
| 6 | 6.24 | 12 | 18.8 | 19.70 | 34.02 | 22.39 | 0.36 | 0.27 | 8 |
| 7 | 6.24 | 21 | 19.4 | 20.29 | 44.23 | 18.66 | 0.33 | 0.24 | 3 |
| 8 | 6.24 | 39 | 20.6 | 19.11 | 23.81 | 16.79 | 0.27 | 0.3 | 13.8 |
| 9 | 11.82 | 30 | 18.8 | 19.11 | 44.23 | 18.66 | 0.27 | 0.27 | 20 |
| 10 | 11.82 | 12 | 20 | 20.29 | 23.81 | 16.79 | 0.33 | 0.33 | 8 |
| 11 | 11.82 | 21 | 20.6 | 19.70 | 13.61 | 20.53 | 0.36 | 0.3 | 3 |
| 12 | 11.82 | 39 | 19.4 | 18.52 | 34.02 | 22.39 | 0.3 | 0.24 | 13.8 |
| 13 | 14.61 | 30 | 18.8 | 20.29 | 23.81 | 20.53 | 0.36 | 0.24 | 13.8 |
| 14 | 14.61 | 12 | 20 | 19.11 | 44.23 | 22.39 | 0.3 | 0.3 | 3 |
| 15 | 14.61 | 21 | 20.6 | 18.52 | 34.02 | 18.66 | 0.27 | 0.33 | 8 |
| 16 | 14.61 | 39 | 19.4 | 19.70 | 13.61 | 16.79 | 0.33 | 0.27 | 20 |
| 17 | 9.03 | 30 | 20.6 | 19.70 | 44.23 | 16.79 | 0.3 | 0.24 | 8 |
| 18 | 9.03 | 12 | 19.4 | 18.52 | 23.81 | 18.66 | 0.36 | 0.3 | 20 |
| 19 | 9.03 | 21 | 18.8 | 19.11 | 13.61 | 22.39 | 0.33 | 0.33 | 13.8 |
| 20 | 9.03 | 39 | 20 | 20.29 | 34.02 | 20.53 | 0.27 | 0.27 | 3 |
| 21 | 6.24 | 30 | 20.6 | 18.52 | 23.81 | 22.39 | 0.33 | 0.27 | 3 |
| 22 | 6.24 | 12 | 19.4 | 19.70 | 44.23 | 20.53 | 0.27 | 0.33 | 13.8 |
| 23 | 6.24 | 21 | 18.8 | 20.29 | 34.02 | 16.79 | 0.3 | 0.3 | 20 |
| 24 | 6.24 | 39 | 20 | 19.11 | 13.61 | 18.66 | 0.36 | 0.24 | 8 |
| 25 | 11.82 | 30 | 19.4 | 19.11 | 34.02 | 16.79 | 0.36 | 0.33 | 3 |
| 26 | 11.82 | 12 | 20.6 | 20.29 | 13.61 | 18.66 | 0.3 | 0.27 | 13.8 |
| 27 | 11.82 | 21 | 20 | 19.70 | 23.81 | 22.39 | 0.27 | 0.24 | 20 |
| 28 | 11.82 | 39 | 18.8 | 18.52 | 44.23 | 20.53 | 0.33 | 0.3 | 8 |
| 29 | 14.61 | 30 | 19.4 | 20.29 | 13.61 | 22.39 | 0.27 | 0.3 | 8 |
| 30 | 14.61 | 12 | 20.6 | 19.11 | 34.02 | 20.53 | 0.33 | 0.24 | 20 |
| 31 | 14.61 | 21 | 20 | 18.52 | 44.23 | 16.79 | 0.36 | 0.27 | 13.8 |
| 32 | 14.61 | 39 | 18.8 | 19.70 | 23.81 | 18.66 | 0.3 | 0.33 | 3 |
Orthogonal test table.
The normalized parameters calculated by Equation 5 are as follows Equation 6.
3 Project overview
Figure 4 shows the plan view of a newly constructed foundation pit. The pit measures 212.4 m in length, 26.4 m in width, and 16.7 m in depth. According to the geotechnical investigation report, the site is composed of silty clay, sandy soil, and argillaceous siltstone from top to bottom (Figure 5). A representative existing building, located approximately 13.8 m from the foundation pit, is a six-story shear wall frame structure. Prior to excavation, construction personnel observed a continuous crack 33 mm deep at the junction between the shear wall and the top slab on the top floor of the building (Figure 6), along with other minor cracks. Clearly, this building is already damaged, and the disturbance caused by new foundation pit excavation is likely to cause further or accelerated damage. To ensure the safety of the existing structure during foundation pit excavation, a thorough assessment and analysis are necessary.
FIGURE 4
FIGURE 5
FIGURE 6
Although the building has multiple cracks, only one crack is relatively deep, classified as a medium-to-deep crack. Such medium-to-deep cracks are typically caused by factors such as significant shear forces, bending moments, or uneven settlement acting on shear walls. If not addressed in a timely manner, these cracks may further develop and affect the structural stability. The other cracks are relatively shallow and are considered surface cracks, which generally have little impact on the structural safety. Therefore, after considering the computational efficiency and cost, this study only takes the most severe crack into account in the calculations, while neglecting the other minor cracks.
4 Analysis and model application
4.1 FEM-XFEM analysis
As shown in
Figure 7, the crack propagation area at the top of the shear wall was extracted, and the results for each group are displayed in
Figure 8. The analysis reveals the following:
Figure 9.
(1) In groups 12 and 16, no crack propagation was observed. In contrast, Group 2 experienced the most severe crack propagation, with a crack length of 9.6 m and an average depth of 0.156 m. The other groups displayed varying degrees of damage propagation. Notably, in several cases, cracks extended to the top slab of the building, with the total crack length in Group 22 exceeding 14 m.
(2) As the building’s differential settlement increases, stress concentration occurs at the crack tip. When the stress intensity factor at the crack tip reaches the fracture toughness of the concrete, cracks begin to form (). Cracks in the shear wall significantly alter the stress distribution in the top slab, eventually leading to crack propagation into the slab. Figure 9 shows the computational results for the top slab under varying degrees of damage. The figure shows that the crack initiates above the shear wall and gradually extends toward the center of the building. The specific distribution and propagation of cracks are detailed in Table 5.
FIGURE 7
FIGURE 8
FIGURE 9
TABLE 5
| Test groups | Initial crack extension area (m2) | Roof crack area (m2) | Total crack extension area (m2) | Test groups | Initial crack extension area (m2) | Roof crack area (m2) | Total crack extension area (m2) |
|---|---|---|---|---|---|---|---|
| 1 | 1.334 | 0.248 | 1.582 | 17 | 0.927 | 1.604 | 2.531 |
| 2 | 1.498 | 0.64 | 2.138 | 18 | 0.576 | 0 | 0.576 |
| 3 | 0.94 | 0.158 | 1.098 | 19 | 0.986 | 0.41 | 1.396 |
| 4 | 0.915 | 0.56 | 1.475 | 20 | 0.79 | 0.208 | 0.998 |
| 5 | 0.664 | 0.23 | 0.894 | 21 | 0.833 | 0.168 | 1.001 |
| 6 | 0.99 | 0.16 | 1.15 | 22 | 1.209 | 2.86 | 4.069 |
| 7 | 1.268 | 0.101 | 1.369 | 23 | 1.111 | 0.78 | 1.891 |
| 8 | 0.463 | 0 | 0.463 | 24 | 0.714 | 0.974 | 1.688 |
| 9 | 0.69 | 0.24 | 0.93 | 25 | 0.789 | 0.574 | 1.363 |
| 10 | 0.953 | 0.194 | 1.147 | 26 | 0.511 | 0.096 | 0.607 |
| 11 | 1.044 | 0.268 | 1.312 | 27 | 0.7 | 0 | 0.7 |
| 12 | 0 | 0 | 0 | 28 | 0.789 | 0.24 | 1.029 |
| 13 | 0.368 | 0 | 0.368 | 29 | 0.899 | 0.161 | 1.06 |
| 14 | 1.328 | 0.328 | 1.656 | 30 | 0.221 | 0 | 0.221 |
| 15 | 1.037 | 0.24 | 1.277 | 31 | 0.154 | 0 | 0.154 |
| 16 | 0 | 0 | 0 | 32 | 0.702 | 0.26 | 0.962 |
Summary of the orthogonal test results.
4.2 RSM prediction model
The response surface model was constructed based on the orthogonal test results, and the regression coefficients were calculated using ridge regression in MATLAB. Table 6 presents the regression coefficients obtained from the model (the intercept is 1.896), which describe the influence of each factor on crack propagation.
TABLE 6
| Linear coefficient () | −0.618 | −0.317 | −0.162 | 1 | −0.692 | 0.039 | −0.704 | −0.498 | −0.166 |
| Quadratic coefficient () | −0.437 | 0.021 | −0.046 | 0.09 | 0.505 | 0.054 | −0.065 | 0.575 | −0.281 |
| Cubic coefficient () | 0.134 | −0.237 | 0.078 | −0.805 | 0.771 | −0.201 | 0.329 | 0.487 | −0.104 |
Regression coefficients for the crack propagation prediction model.
5 Results and discussion
5.1 FEM-XFEM and RSM model results
The comparison between the numerical simulation results of 32 working conditions and the values predicted from the response surface function is shown in Figure 10. In civil engineering and fracture mechanics, data typically exhibit high variability and uncertainty, and if the mean squared error (MSE) is less than 0.4, the fit is considered good (; ; ). In this study, the ridge regression MSE was 0.15, confirming the reliability of the predictions for crack propagation in existing buildings.
FIGURE 10
To verify the predictive effectiveness of the above model, four random working conditions, as shown in Table 7, were selected. Finite element calculations were performed to obtain the damage development for each condition (Figure 11), specifically the crack propagation area, and the results were compared with those predicted by the model, as shown in Figure 12. The analysis reveals that the damage trends predicted by the model closely match the numerical simulation results, indicating that the established response surface model can effectively predict damage propagation in buildings affected by foundation pit excavation.
TABLE 7
| Test groups | (MPa) | (MPa) | (kN/m3) | (kN/m3) | (kPa) | (°) | (m) | ||
|---|---|---|---|---|---|---|---|---|---|
| 33 | 7.08 | 17.40 | 19.34 | 19.26 | 28.92 | 20.15 | 0.333 | 0.312 | 20 |
| 34 | 13.77 | 33.60 | 20.06 | 19.62 | 28.92 | 19.03 | 0.297 | 0.258 | 3 |
| 35 | 10.43 | 22.80 | 19.34 | 18.91 | 16.67 | 17.91 | 0.297 | 0.276 | 8 |
| 36 | 7.08 | 17.40 | 19.34 | 19.26 | 28.92 | 19.03 | 0.297 | 0.258 | 13.8 |
Parameters of the validation conditions.
FIGURE 11
FIGURE 12
The FEM-XFEM model provides detailed, high-accuracy predictions of crack initiation and propagation, but it is computationally expensive, particularly when simulating a large number of scenarios. In contrast, the RSM model offers a computationally efficient alternative, though it sacrifices some precision for speed. The comparison between the two models reveals that both models are effective, but the choice between them depends on the specific application—FEM-XFEM is preferable for detailed analysis, while RSM is more suitable for rapid predictions across multiple scenarios.
5.2 Significance analysis of influencing factors
To analyse the sensitivity of building damage to various influencing factors, further calculations were conducted to determine the Pearson Product-Moment Correlation Coefficient (PPMCC) and the Maximal Information Coefficient (MIC) between the total crack propagation area and each parameter. The PPMCC measures the linear relationship between variables, with values ranging from −1 to 1, whereas the MIC measures the nonlinear relationship between variables, with values ranging from 0 to 1. The results are shown in
Figure 13, and the analysis reveals the following:
(1) In terms of the MIC, the distance between the foundation pit and the building () has the strongest correlation (0.828), followed by the elastic modulus of the sandy soil layer (, 0.596) and the Poisson’s ratio of the silty clay layer (, 0.419). These parameters significantly affect the crack propagation area in a nonlinear manner. In terms of linear correlation (PPMCC), the elastic modulus of the silty clay layer (, −0.45) and the distance between the foundation pit and the building (, 0.27) show strong linear correlations, with both exhibiting a negative correlation.
(2) Considering both MIC and PPMCC, shows high correlation in both metrics and is the most influential factor. This is followed by and , with being stronger in the nonlinear dimension and being stronger in the linear dimension. The overall ranking of influence is as follows: .
(3) Some parameters, such as , show strong negative linear correlation but weak nonlinear correlation (0.423). Conversely, the Poisson’s ratio of the sandy soil layer () shows notable correlation in MIC (0.306), but relatively smaller positive correlation in PPMCC (0.26). These differences suggest that certain parameters exert varying levels of influence in linear and nonlinear dimensions, indicating the need to consider the multidimensional effects of these factors to accurately predict the crack propagation area.
FIGURE 13
5.3 Sensitivity analysis of influencing factors in the prediction model
To analyse the influence of each factor on building crack propagation individually, all other factors were set to their average values while the target factor was varied. For each target factor, standardized values were selected at intervals of 0.01. This study involves nine influencing factors, resulting in a total of 900 scenarios. The prediction model from
Section 4.2was used to calculate the crack propagation area for each scenario.
Figure 14illustrates the impact of each of the nine influencing factors on crack propagation. The analysis reveals:
(1) The trends for , , , and are quite similar, with the crack propagation area decreasing or stabilizing as the independent variable increases. Notably, and show the most significant impact, with the crack propagation area decreasing from approximately 1.6 m2 to approximately 0.4 m2. A higher elastic modulus means less soil deformation under the same construction load, resulting in a smaller impact on existing buildings. Additionally, indicated that attention should be given to crack development in buildings closer to foundation pits. On the other hand, , , and significantly increase in the crack propagation area as the factor increases, with a “turning point” () where the crack area reaches its maximum.
(2) Identifying high-sensitivity regions in crack propagation analysis is crucial. In addition to , , , and , which exhibit high sensitivity across the entire range, and also show high sensitivity in the range [0.6, 1.0], where small increases in the factor cause a sharp increase in the crack propagation area. Furthermore, demonstrates a high sensitivity region between [0.2, 0.4]. These high-sensitivity regions indicate that under specific conditions, small changes in these parameters can have a significant effect on crack propagation. Therefore, in construction, attention should be focused on controlling crack propagation within these high-sensitivity regions.
FIGURE 14
6 Conclusion
This study aimed to address two key research questions:
(1) How does foundation pit excavation impact crack propagation in adjacent buildings?
(2) Can a predictive model effectively quantify the influence of soil properties and distance on crack behavior?
The investigation utilized FEM-XFEM and RSM models to analyze the crack propagation induced by foundation pit excavation. The results demonstrate the significant role of soil properties and spatial parameters in influencing structural damage. The findings highlight that both models effectively predict crack propagation, with the RSM providing a computationally efficient framework for large-scale parameter studies.
(1) This work selects nine parameters, including the distance between the foundation pit and the building, the elastic modulus of the soil, its density, and Poisson’s ratio, as independent variables, and the crack propagation area of existing adjacent buildings as the dependent variable. On the basis of an orthogonal test, 32 working conditions were set up, and the crack propagation areas for different conditions were calculated via the FEM-XFEM model. The results of the orthogonal test indicated that the elastic modulus, Poisson’s ratio of the soil, and distance between the building and the foundation pit significantly influence crack propagation in buildings. Building crack propagation was negatively correlated with the soil elastic modulus, the Poisson’s ratio of silty clay, and the distance from the pit and positively correlated with the Poisson’s ratio of sandy soil and the cohesion of silty clay.
(2) Using the ridge regression model, the results from the 32 conditions were fitted to create a predictive model for building crack propagation. The MSE results confirmed a good fit, and the effectiveness of the predictive model was further validated with four additional working conditions.
(3) The predictive model was used to conduct a single-factor analysis of the nine influencing factors, and the effects of these factors were visualized. The results showed that the behavior of , , , and is close to linear, while the behavior of , , and reveals an extremum in crack propagation area (). Additionally, certain factors showed high-sensitivity intervals for crack propagation area [ and in the range (0.6, 1.0), and in the range (0.2, 0.4)], indicating that particular attention should be given to crack propagation in buildings within these intervals during construction.
While this study provides valuable insights into the relationship between foundation pit excavation and crack propagation, several limitations should be acknowledged. First, the use of simplified assumptions regarding soil behavior and crack propagation limits the ability to capture the full complexity of real-world conditions. Additionally, the absence of comprehensive field data for model validation restricts the study’s generalizability. Future research should focus on incorporating dynamic loading effects, such as seismic activity, and explore interactions between multiple cracks and heterogeneous soil layers. Collaborative efforts to collect real-world excavation data would significantly enhance the model’s applicability and reliability.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
ZW: Conceptualization, Funding acquisition, Resources, Writing – original draft. TW: Methodology, Resources, Writing – original draft. YK: Data curation, Investigation, Validation, Writing – original draft. ZZ: Investigation, Project administration, Validation, Writing – original draft. FX: Formal Analysis, Writing – original draft. PY: Data curation, Formal Analysis, Writing – original draft, Writing – review and editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This research was funded by “Research on Construction Techniques for Subway Deep Foundation Pits Near Structures and Buildings (2021)” and “Science and Technology Plan Project of China Construction Third Engineering Bureau (2022)”.
Conflict of interest
Author ZW was employed by Metro Project Management Branch of Nanchang Rail Transit Group Co., Ltd. Author TW was employed by CCTEB Infrastructure Construction Investment Co., Ltd. Authors YK and ZZ were employed by China Construction Third Engineering Bureau Group Co., Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that Generative AI was used in the creation of this manuscript. Generative AI was used during the preparation of this work the author(s) used ChatGPT in order to retranslation and optimization. After using ChatGPT, the author(s) reviewed and edited the content as needed and take(s) full responsibility for the content of the publication.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Abbreviations
, , Elastic modulus of the silty clay layer and sandy soil layer; , , Density of the silty clay layer and sandy soil layer; , , Poisson’s ratio of the silty clay layer and sandy soil layer; , , Cohesion and internal friction angle of the silty clay layer; , Distance between the foundation pit and the building; Isotropic damage variable; , , Nominal and effective load-bearing area; , Crack propagation area; , , Length and average depth of the crack; , The i-th influencing factor; , Coefficient of ; , Error term; , Normalized data; , Original data; , , Maximum and minimum value of .
References
1
ArapakouA. E.PapadopoulosV. P. (2012). Factors affecting differential settlements of framed structures. Geotech. Geol. Eng.30 (6), 1323–1333. 10.1007/s10706-012-9546-x
2
ArshadA. R. A. (2016). A review on selection of tunneling method and parameters effecting ground settlements. Electron. J. Geotechnical Eng.21, 4459–4475. 10.1016/j.tust.2020.103614
3
ÁvilaF.PuertasE.TorrúsC.GallegoR. (2024). Influence of crack propagation on the seismic behavior of historic rammed earth buildings: the Tower of Muhammad in the Alhambra (Spain). Eng. Struct.10.1016/j.engstruct.2023.117365
4
CaoC.ShiC.LiuL.LiuJ. (2020). Evaluation of the effectiveness of an alternative to control groundwater inflow during a deep excavation into confined aquifers. Environ. EARTH Sci.79 (22), 502. 10.1007/s12665-020-09253-3
5
CaoJ.LiuT.LiuH. M. (2014). Soil nailing reliability optimization design of foundation pit based on response surface method. Appl. Mech. Mater.556-562, 4655–4659. 10.4028/www.scientific.net/amm.556-562.4655
6
ChaiM.LiuP.HeY.HanZ.DuanQ.SongY.et al (2023). Machine learning‐based approach for fatigue crack growth prediction using acoustic emission technique. FATIGUE Fract. Eng. M.46 (8), 2784–2797. 10.1111/ffe.14032
7
ChenL.MaH. (2018). Sensitivity analysis of influencing factors of deep foundation pit engineering support structure design in soft soil area. AIP Conf. Proc.1995 (1), 020027. 10.1063/1.5048758
8
DengY.GuoY.ZouB.WangJ.LiuX.ChenQ.et al (2023). Failure analysis and zoning control of water gushing in foundation pit. Eng. FAIL Anal.145, 107029. 10.1016/j.engfailanal.2022.107029
9
DingL. Y.WuX. G.LiH.LuoH. B.ZhouY. (2011). Study on safety control for Wuhan metro construction in complex environments. Int. J. Proj. Manag.29 (7), 797–807. 10.1016/j.ijproman.2011.04.006
10
DmochowskiG.SzolomickiJ. (2021). Technical and structural problems related to the interaction between a deep excavation and adjacent existing buildings. Appl. Sci.11 (2), 481. 10.3390/app11020481
11
GüçlüerK.ÖzbeyazA.GöymenS.GünaydınO. (2021). A comparative investigation using machine learning methods for concrete compressive strength estimation. MATER TODAY Commun.27, 102278. 10.1016/j.mtcomm.2021.102278
12
HanM.LiZ.MeiG.BaoX.JiaJ.LiuL.et al (2022). Characteristics of subway excavation in soft soil and protective effects of partition wall on the historical building and pile foundation building. B Eng. Geol. Environ.81 (8), 307. 10.1007/s10064-022-02802-z
13
HanW.LiG.SunZ.LuanH.LiuC.WuX. (2020). Numerical investigation of a foundation pit supported by a composite soil nailing structure. SYMMETRY-BASEL12 (2), 252. 10.3390/sym12020252
14
HsiehP.OuC. (1998). Shape of ground surface settlement profiles caused by excavation. Can. Geotech. J.35, 1004–1017. 10.1139/t98-056
15
HuZ.LeiM.ZhangH.JiaC.ShiC. (2024). Ultimate support pressure determination for shield tunnel faces in saturated strata based on seepage flow control. Appl. Sci.14 (15), 6781. 10.3390/app14156781
16
HuZ.WangQ.YangS.ShiZ.LiuB.SongH.et al (2021). Numerical simulation of soil displacement and settlement in deep foundation pit excavations near water. GEOFLUIDS2021, 1–14. 10.1155/2021/5559009
17
JafarifarN.PilakoutasK.BennettT. (2016). The effect of shrinkage cracks on the load bearing capacity of steel-fibre-reinforced roller-compacted-concrete pavements. MATER Struct.49 (6), 2329–2347. 10.1617/s11527-015-0652-0
18
KachanovM. (1992). Effective elastic properties of cracked solids: critical review of some basic concepts. Appl. Mech. Rev.8 (45), 304–335. 10.1115/1.3119761
19
KaklauskasG.GhaboussiJ. (2001). Stress‒strain relations for cracked tensile concrete from RC beam tests. J. Struct. Eng. (New York, N.Y.)127 (1), 64–73. 10.1061/(asce)0733-9445(2001)127:1(64)
20
KimC. Y.BaeG. J.HongS. W.ParkC. H.MoonH. K.ShinH. S. (2001). Neural network based prediction of ground surface settlements due to tunnelling. Comput. Geotech.28 (6-7), 517–547. 10.1016/s0266-352x(01)00011-8
21
KorffM.MairR. J.Van TolF. A. F. (2016). Pile‒soil interaction and settlement effects induced by deep excavations. J. Geotech. GEOENVIRON142 (8). 10.1061/(asce)gt.1943-5606.0001434
22
KumarS.BaraiS. V. (2009). Determining double-K fracture parameters of concrete for compact tension and wedge splitting tests using weight function. Eng. Fract. Mech.76 (7), 935–948. 10.1016/j.engfracmech.2008.12.018
23
LaxmanK. C.TabassumN.AiL.ColeC.ZiehlP. (2023). Automated crack detection and crack depth prediction for reinforced concrete structures using deep learning. Constr. Build. MATER370, 130709. 10.1016/j.conbuildmat.2023.130709
24
LeiM.LiuL.LinY.LiJ. (2019). Study on flexural bearing capacity calculation method of enclosure pile with partial excision in deep foundation pit. Adv. Civ. Eng.2019 (1). 10.1155/2019/4812857
25
LiW.WuM.HuB.WangP. (2024). Overcoming the brittleness of shear failure: a new FRP-RSC strengthening philosophy. COMPOS Struct.330, 117857. 10.1016/j.compstruct.2023.117857
26
LiY.OhkuboT.TeramotoA.SagaK.KawashimaY. (2022). Lab-scale reproduction test method for temperature-driven movement of through-thickness cracks in concrete exterior walls for crack repair evaluation. Constr. Build. MATER331, 127169. 10.1016/j.conbuildmat.2022.127169
27
LiangR.XiaT.HuangM.LinC. (2017). Simplified analytical method for evaluating the effects of adjacent excavation on shield tunnel considering the shearing effect. Comput. Geotech.81, 167–187. 10.1016/j.compgeo.2016.08.017
28
LinF.SongX.GuX.PengB.YangL. (2012). Cracking analysis of massive concrete walls with cracking control techniques. Constr. Build. MATER31, 12–21. 10.1016/j.conbuildmat.2011.12.086
29
LiuB.WuW.LuH.ChenS.ZhangD. (2024). Effect and control of foundation pit excavation on existing tunnels: a state-of-the-art review. Tunn. Undergr. Sp. Tech.147, 105704. 10.1016/j.tust.2024.105704
30
LiuC.PengZ.CuiJ.HuangX.LiY.ChenW. (2023). Development of crack and damage in shield tunnel lining under seismic loading: refined 3D finite element modeling and analyses. THIN WALL Struct.185, 110647. 10.1016/j.tws.2023.110647
31
LiuJ.ShiC.CaoC.LeiM.WangZ. (2020). Improved analytical method for pile response due to foundation pit excavation. Comput. Geotech.123, 103609. 10.1016/j.compgeo.2020.103609
32
LiuL.LeiM.CaoC.ShiC. (2019). Dewatering characteristics and inflow prediction of deep foundation pits with partial penetrating curtains in sand and gravel strata. WATER-SUI11 (10), 2182. 10.3390/w11102182
33
LiuZ.XueJ.YeJ.QianJ. (2021). A simplified two-stage method to estimate the settlement and bending moment of upper tunnel considering the interaction of undercrossing twin tunnels. Transp. Geotech.29, 100558. 10.1016/j.trgeo.2021.100558
34
LuB. (2021). Study on back analysis methodof soil parameter uncertainty in excavation deformation calculation of foundation pit. Nanjing: Southeast University.
35
LuoP.LiC. (2019). Numerical simulation analysis of the influence of foundation pit excavation on the surrounding buildings of an office building. lOP Conf. Ser. Earth Environ. Sci.252, 052077. 10.1088/1755-1315/252/5/052077
36
OuW.ChenX.ChanA.ChengY.WangH. (2022). FDEM simulation on the failure behavior of historic masonry heritages subjected to differential settlement. BUILDINGS-BASEL.12 (10), 1592. 10.3390/buildings12101592
37
SalamiB. A.OlayiwolaT.OyehanT. A.RajiI. A. (2021). Data-driven model for ternary-blend concrete compressive strength prediction using machine learning approach. Constr. Build. MATER301, 124152. 10.1016/j.conbuildmat.2021.124152
38
ShiC.CaoC.LeiM.PengL.AiH. (2016). Effects of lateral unloading on the mechanical and deformation performance of shield tunnel segment joints. Tunn. Undergr. Sp. Tech.51, 175–188. 10.1016/j.tust.2015.10.033
39
ShiC.SunX.LiuS.CaoC.LiuL.LeiM. (2021). Analysis of seepage characteristics of a foundation pit with horizontal waterproof curtain in highly permeable strata. WATER-SUI13 (9), 1303. 10.3390/w13091303
40
Sou-SenL.Hsien-ChuangL. (2004). Neural-network-based regression model of ground surface settlement induced by deep excavation. Autom. Constr.13 (3), 279–289. 10.1016/s0926-5805(03)00018-9
41
TircaLLAH (2017). High-rise building subjected to excessive settlement of its foundation-A case study. lnternational J. Struct. lntegrity2 (8). 10.1108/ijsi-05-2016-0019
42
WangG.ChenW.CaoL.LiY.LiuS.YuJ.et al (2021). Retaining Technology for deep foundation pit excavation adjacent to high-speed railways based on deformation control. Front. EARTH SC-SWITZ9. 10.3389/feart.2021.735315
43
WangK.YangZ.GuoJ.DangY.YanY. (2023a). Numerical analysis of the influence of deep excavation on nearby pile foundation building. BUILDINGS-BASEL13 (11), 2842. 10.3390/buildings13112842
44
WangP.WuH.KeL.LiW.LeungC. K. Y. (2023b). Mechanical properties and microstructure of glass fiber reinforced polymer (GFRP) rebars embedded in carbonated reactive MgO-based concrete (RMC). Cem. Concr. Compos.142, 105207. 10.1016/j.cemconcomp.2023.105207
45
WeiH.Sang-BingT.TsaiS. (2021). Influence of foundation pit excavation and precipitation on settlement of surrounding buildings. Adv. Civ. Eng.2021 (1). 10.1155/2021/6638868
46
XuQ.XieJ.LuL.WangY.WuC.MengQ. (2024). Numerical and theoretical analysis on soil arching effect of prefabricated piles as deep foundation pit supports. Undergr. SPACE16, 314–330. 10.1016/j.undsp.2023.09.011
47
XueH. (2023). Research on the control of excavation deformation of super deep foundation pit adjacent to the existing old masonry structure building. SUSTAINABILITY-BASEL15 (9), 7697. 10.3390/su15097697
48
YanX.TongL.LiH.HuangS.FuY.WangW.et al (2024). Effects of the excavation of deep foundation pits on adjacent approach bridges: a case study of Nanjing Yangtze River Bridge. J. Civ. Struct. HEALTH15, 1185–1205. 10.1007/s13349-024-00875-w
49
YangP.LeiM.HuZ.XiaoF.JiangT.ZhongQ. (2024). Analysis of the impact of new pit construction projects on the crack expansion of existing buildings. Eng. FAIL Anal.162, 108431. 10.1016/j.engfailanal.2024.108431
50
YeS.ZhaoZ.WangD. (2021). Deformation analysis and safety assessment of existing metro tunnels affected by excavation of a foundation pit. Undergr. SPACE6 (4), 421–431. 10.1016/j.undsp.2020.06.002
51
ZhangD. (2023). Influences of deep foundation pit excavation on the stability of adjacent ancient buildings. BUILDINGS-BASEL13 (8), 2004. 10.3390/buildings13082004
52
ZhaoC.FengY.WangW.NiuZ. (2023). Mechanical properties and numerical analysis of underground continuous wall in underground grain silo foundation pit. BUILDINGS-BASEL.13 (2), 293. 10.3390/buildings13020293
53
ZhongW.HuangN.NieN.DingH.GaoF. (2023). Study on the influence of excavation of superlarge and ultra-deep foundation pits on the pile foundation of existing viaducts. Adv. Civ. Eng.2023, 1–13. 10.1155/2023/5834958
54
ZhouS.ChengX.LiQ.ZhengG.GuoJ.MaY.et al (2024). Failure analysis of water gushing in excavation and application of rapid dewatering and recharge emergency measures. Eng. FAIL Anal.159, 108074. 10.1016/j.engfailanal.2024.108074
Summary
Keywords
impact of foundation pit excavation, crack propagation prediction, damage in adjacent buildings, response surface methodology, finite element method
Citation
Wu Z, Wang T, Kuang Y, Zhu Z, Xiao F and Yang P (2025) A combined XFEM and RSM approach for predicting crack propagation in nearby buildings during foundation pit construction. Front. Built Environ. 11:1514217. doi: 10.3389/fbuil.2025.1514217
Received
20 October 2024
Accepted
21 April 2025
Published
01 May 2025
Volume
11 - 2025
Edited by
Peng Wang, Shenzhen University, China
Reviewed by
Masoud Salar, Chabahar Maritime University, Iran
Georgios Eleftherios Stavroulakis, Technical University of Crete, Greece
Tianyu Wang, Hong Kong University of Science and Technology, Hong Kong SAR, China
Shangtao Hu, Guangzhou University, China
Updates
Copyright
© 2025 Wu, Wang, Kuang, Zhu, Xiao and Yang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
* Correspondence: Pengliang Yang, ypengliang@126.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.