1. Introduction
In recent years, the reduction in suitable urban spaces and the need for large-scale housing for the growing population have led to an increase in the construction of buildings with minimal distances between them. This has resulted in the accumulation of ground stress and an increase in deformation within the foundations of adjacent buildings. Existing regulations, such as FEMA 273 [
1], allow for the consideration of this effect in foundation design by adjusting the subgrade reaction modulus in 2D analysis, but it needs to be expanded and modified in terms of 3D analysis [
2]. Additionally, the interaction between buildings has the potential to impact foundation performance. Moreover, soil–structure interaction (SSI) is influenced by the flexibility of the soil below the foundation and the relative vibrations between the foundation and the free surface. By accounting for these effects, it is possible to determine the inertial forces and actual displacements of a structure under seismic ground motion or quasi-dynamic loads [
3].
Traditionally, the effects of soil deformability are neglected even if the motion of the foundation differs from the free-field ground motion [
4]. Wong et al. [
4] and Luco et al. [
5] investigated dynamic structure–soil–structure interaction (SSSI) between two shear walls placed on rigid foundations subjected to SH waves. They demonstrated that structures with high frequencies positioned near each other have minimal mutual influence, whereas structures with lower frequencies can exert significant effects on one another [
5,
6]. As the number of structures increases, the maximum displacement experiences an increase due to SSSI effects. This depends on the frequencies of the structures that are relatively close, and their masses are larger compared to the foundation area, thereby magnifying interaction effects [
6]. However, these effects were previously neglected in a study [
7] when Young’s modulus exceeded 6.9 × 106 kPa. Various types of soil were analysed to consider SSSI effects, including sandy and soft soils [
8]. These types of soil amplify seismic wave interactions with soil–structure systems. In most studies, maximum principal stresses were observed on the contact surface between the soil and the foundation beneath the columns, whereas minimum stress values appeared at the centre of the foundations [
9]. Çelebi et al. [
10] observed a noticeable change in the seismic response of structures with only a few stories when SSI effects were disregarded in loose soil. Furthermore, there is a significant discrepancy between analyses using linear soil behaviour and a nonlinear elasto-plastic Mohr–Coulomb soil behaviour.
Cayci et al. [
11] indicated that linear models may provide inaccurate estimates of soil behaviour when subjected to seismic loads. The influence of neighbouring buildings can be particularly strong for embedded foundations, leading to reduced responses at the top of the buildings [
12]. Alexander et al. [
13] demonstrated that the effect of an earthquake is often transferred from a taller new building to a shorter, older one. Buildings with higher aspect ratios (tall and narrow) are more susceptible to failure on medium and dense soils, while buildings with very small aspect ratios (short and wide) are at higher risk on loose soils in terms of SSSI effects. When there is a significant difference in height (ε > 1.5, where ε is the relation between the height of a taller building to a shorter one) between buildings, maximum building displacement occurs, and it may increase due to earthquake effects transmitted from taller buildings to shorter ones [
14].
Some researchers have explored these effects in the context of three or more buildings [
4,
6,
15,
16,
17], demonstrating that these effects are more pronounced compared to situations with only two adjacent buildings. Aldaikh et al. [
16] revealed that SSSI effects are more severe when buildings are in close proximity, and these effects can be considered negligible when the spacing between buildings is more than 2.5 times the foundation width [
15]. Li and Liu [
17] noted that SSSI effects become significant when the spacing between adjacent buildings is less than 9 m (or 0.3 × B). SSSI effects should also be taken into account in seismic designs [
18]. Aldaikh et al. [
15] emphasized that the height of structures plays a significant role in terms of SSSI effects under seismic loads. Ignoring SSI and SSSI effects may lead to erroneous estimations of seismic capacity and drift parameters, as SSI effects increase drift parameters while reducing seismic capacity [
19]. Pile foundations could be a reliable strategy for mitigating these effects [
18]. Furthermore, SSSI effects increase the maximum drift in adjacent buildings, but when the spacing between buildings is larger and the columns are made of stone, the maximum drift is reduced [
20].
The effects of soil–structure interaction can be considered using both direct and substructure methods. In the direct method, the structure and the soil are modelled simultaneously and analysed together, allowing for the consideration of nonlinear soil behaviour [
21,
22,
23]. The substructure method consists of dividing the two sub-systems that are analysed separately, and the results are combined in the final analysis stage by assuming linear behaviours.
There is currently a gap in using both SSI and SSSI effects on foundation design in engineering projects [
24]. The integration and application of both soil–structure interaction (SSI) and structure–soil–structure interaction (SSSI) effects within the realm of structural and geotechnical engineering have been considered in this study. This paper introduces a methodology that utilizes diverse subgrade reaction modulus values, as outlined in the provided formulation and table, to account for and apply these effects in the analysis and design process. This approach represents a step forward in bridging the gap between structural and geotechnical considerations, ensuring a more holistic and accurate assessment of the interaction between buildings and their underlying soft soil.
2. Materials and Methods
Three different plans were assessed to achieve the required results, each proposing varying locations for the shear walls. All structures were analysed to obtain the optimum locations of the shear walls via SAP 2000 and ETABS 9.2 software (Computers and Structures, INC., New York, NY, USA) (see
Figure 1).
The structural models of the buildings presented in this article are described below.
Two buildings A1 and A2 were considered, each with the same plan dimensions of 18 m × 12 m (
Figure 1a). A1 had a height of 36 m (12-story), while A2 had a height of 18 m (6 storeys), creating a significant height difference (ε > 1.5, ε is representative of the proportion of building heights) [
13]. The irregular distribution of mass and stiffness components in building plans introduces complexities in how these structures interact with the underlying soil. This asymmetry amplifies the sensitive behaviour of the buildings under varying loads and environmental conditions. The significance lies in its potential to significantly alter deformation patterns and load transfer mechanisms, consequently impacting the accuracy of SSI and SSSI analyses. A concrete elasticity modulus E = 25.2 GPa and Poisson’s ratio υ = 0.2 were assumed for both structures. To model the beams and columns of the buildings, a three-dimensional element (BEAM 4) with six degrees of freedom at each node (three for translational movements and three for rotational movements) was employed. The slabs were modelled using a four-node element (SHELL 181), where each node also possessed six degrees of freedom. The SHELL 181 element was also used to model the shear walls [
2].
The foundation model for both buildings had dimensions of 13 m × 19 m × 1 m. The modulus of elasticity (E) and Poisson’s ratio (υ) values are presented in
Table 1 for both foundations. To represent the foundation elements, a three-dimensional eight-node SOLID 45 element was utilized, with each node having three degrees of freedom. Furthermore, a four-node contact element (CONTA 173) was employed to model the interface between the foundation and the soil. The following properties in
Table 1 were attributed to the contacting elements. These properties correspond to the coefficient of friction (MU: representing the frictional ratio between contact surfaces), normal stiffness (FKN: indicating the normal stiffness of the contact interface), cohesive strength (COHE: defining adhesive forces between the contact surfaces), and maximum shear stress (TAUMAX: representing the maximum allowable shear stress at the contact interface), respectively [
25].
A 2D soil domain (plan dimensions: 28 m × 34 m and height: 20 m) was considered for the analysis of buildings. The chosen dimensions aim to capture the influence of nearby soil conditions, neighbouring structures, or geological features that could impact the building’s behaviour [
26]. SOLID 45 element was employed to model the soil. The modulus of elasticity for the soils (E) and the Poisson ratio (υ) values are shown in
Table 1. The Drucker–Prager model with an internal friction angle (φ), cohesion (c), and dilation angle (Ψ) was used for soil plasticity. The use of negative dilation angles in soil material modelling allows for the accurate representation of specific soil behaviour, especially those associated with contraction tendencies or minimal lateral strain during shearing, providing a more realistic depiction of soil behaviour in geotechnical analyses [
27].
The selection of specific building dimensions and heights (12-storey building (A1) and 6-storey building (A2)) allows for the examination of a substantial height difference (ε > 1.5). This deliberate difference aids in analysing the impact of varied building heights on structural behaviour and soil–structure interaction. The assumptions of concrete properties (elasticity modulus E = 25.2 GPa and Poisson’s ratio υ = 0.2) for both structures ensure consistency in material behaviour across the models. These values are often standard and representative of typical concrete properties used in construction [
28,
29]. The choice of three-dimensional elements (BEAM 4, SHELL 181, and SOLID 45) with specified degrees of freedom provides an appropriate level of detail and accuracy in modelling the structural components (beams, columns, slabs, and shear walls) and the foundation. These elements offer the necessary complexity to capture realistic behaviour and interactions within the structural system. The foundation’s dimensions and material properties (elasticity modulus E = 25.2 GPa and Poisson’s ratio υ = 0.2) mirror the properties assumed for the building materials [
30]. This consistency ensures compatibility between the structural elements and their supports. The soil area dimensions and properties are chosen to encapsulate a substantial area around the buildings, allowing for an accurate representation of the soil–structure interaction. The selection of the SOLID 45 element for soil modelling, along with the Drucker–Prager model for soil plasticity, enables the simulation of realistic soil behaviour under varying loads and conditions [
27,
31]. The structural discretization was performed using the automatic meshing procedure in ANSYS to create a comprehensive model. To ensure accuracy, a grid sensitivity analysis was conducted, generating five different meshes ranging from 175,352 to 1,649,803 cells for the initial design within ANSYS Meshing software (version 11, Ansys, Inc., Canonsburg, PA, USA). It was observed that the standard deviation of the base support reactions stabilized at approximately 0.82% for the configuration containing 1,123,492 cells. Consequently, this mesh density was selected for subsequent finite element method (FEM) analyses.
In foundation design incorporating structure–soil–structure interaction (SSSI), key boundary conditions encompass structural elements’ loads, fixity, and rigidity, while soil aspects include properties like elasticity, shear strength, and interface characterization. These conditions focus on ensuring compatibility between structural and soil deformations, accounting for the effects in seismic areas, and employing accurate finite element analysis (FEA) parameters [
32,
33]. These considerations underpin a holistic approach to accurately model SSSI, ensuring realistic simulations and enhancing the reliability of foundation design. Flexibility impacts how the structure responds to external loads and how this response interacts with the soil. It influences the extent of settlements, rotations, and differential movements in the foundation, crucial aspects of SSSI. Utilizing advanced structural analysis software (ANSYS) that is capable of accurately modelling various structural elements and soil–structure interaction phenomena would be a method for considering flexibility [
34,
35].
In the first step, the shorter structure (A2) (comprising beams, columns, shear walls, and slabs) with plan dimensions of 18 m × 12 m and a height of 18 m (6 storeys) was independently modelled. It was analysed in ANSYS, and subsequently, its foundation was designed using SAFE geotechnical software (version 8.0.6, Computers and Structures, INC., New York, NY, USA), with the weight of the reinforcement used calculated (further details will be provided later). In the second step, the structure’s foundation was incorporated into the previous model within ANSYS. Following that, the soil and the contact surface between the soil and the foundation were modelled. In the final step, various load combinations were applied to the structure, and the model was analysed.
To model the taller structure (A1), with plan dimensions of 18 m × 12 m and a height of 36 m (12 storey), all steps pertaining to the analysis of the first structure were repeated. After independently modelling the first and second structures and assessing the total deformation associated with the foundation of each structure, both buildings were modelled side by side, with the distance between the two foundations assumed to be zero (i.e., the two foundations are in contact) to account for the most severe effects [
17].
To validate a conceptual design in ANSYS software, the buildings were initially modelled and designed using structural software (ETABS). Subsequently, the sections of beams, columns, and shear walls were modelled in ANSYS software. A comparison was made between the results of the building analysis in ANSYS and ETABS software, with the same loads applied to validate the modelling processes in ANSYS [
2]. The reaction results for A2, represented in terms of moment Mx along the foundation’s width under earthquake load, are presented in
Figure 2 for both ETABS and ANSYS software.
To account for the interaction between adjacent buildings with SSSI effects, both buildings were concurrently modelled in ANSYS. The model was then subjected to various loads, such as gravity loads (comprising dead load and live load) and earthquake loads in the X direction (parallel to the foundation width) and Y direction (parallel to the foundation length). Subsequently, both building foundations were modelled in SAFE for foundation design.
In the first step of foundation design, the foundations were analysed with the same subgrade reaction modulus (Ks). In the second step, they were analysed with different values of subgrade reaction modulus for various sections of the foundations, and the quantities of reinforcement used were calculated for both steps. The variation in Ks values for foundation design was necessary to account for SSI and SSSI effects [
2]. The results of the analysis of foundations with SSSI (employing different Ks values) and without SSSI effects (using a single Ks value) were then compared [
36].
3. Results and Discussion
By combining the gravity loads (dead load and live load) with seismic loads as the inertia loads in both the X and Y directions (seismic loads applied with a coefficient of 0.1 g (0.981 m/s
2) as quasi-dynamic loads in one step time), deformations were obtained for both A2 and A1, both without SSSI effects (analysed separately) and with SSSI effects (when adjacent to each other).
Figure 3 illustrates the total deformation in the Z direction for both A2 and A1 (with the Z-axis representing the vertical direction). These results include the effects of the soil–structure interaction (SSI) and structure–soil–structure interaction (SSSI) resulting from the combination of gravitational loads and seismic loads in both the X and Y directions. SMN represents the maximum stress, SMX indicates the minimum stress values (negative in the direction of applied forces), and DMX denotes the maximum displacement. In both the gravity load combination and earthquake load combination along the X direction, both buildings displayed a noticeable inclination towards each other. This discernible deviation was more accentuated and visible in
Figure 3b. The structural response during seismic loading resulted in a more pronounced convergence of the buildings towards each other along the X-axis. This convergence indicates a significant interaction between the buildings, suggesting increased vulnerability under seismic forces compared to gravitational loads. The observed inclination during the seismic scenario along the X direction underscores the importance of thoroughly assessing structural behaviour, especially under dynamic forces, highlighting the necessity of designing structures for extreme loading conditions. However, there was no leaning between the buildings in terms of earthquake forces along the Y direction, as depicted in
Figure 3c.
The positions of the maximum total deformation contours vary significantly between the two buildings across all load combinations. The maximum deformation occurred within the taller building (A1) in the gravity load combination. In the case of the shorter building (A2), this occurred during the earthquake load combination in the direction of the width of the foundations.
Figure 4 displays the deformation of the foundations for both A2 and A1 in the Z direction, considering interactions resulting from the three load combinations. It is evident that the maximum deformations are in the zone between the foundations and in one of the building corners, with the specific location depending on the load combination directions. Therefore, the adjacent taller building exerts a greater influence on the shorter building when an earthquake load combination in the X-axis is considered (as depicted in
Figure 4b).
Figure 5 illustrates the deformation of the foundations for A2 and A1 in the Z direction, analysed separately without structure–soil–structure interaction (SSSI) effects, due to the combination of gravitational loads. The maximum total deformation is lower in A2 when analysed without SSSI effects (
Figure 5b almost 0.029 m) compared to when SSSI effects are considered (
Figure 4a: almost 0.038 m). Conversely, the difference in the total deformation for A1 with and without SSSI effects is not significantly different (
Figure 4a and
Figure 5a: almost 0.037 m). In fact, the settlement of the taller building (A1) had a greater effect on the settlement of the shorter building (A2) when SSSI effects were considered. These impacts have also changed the area of overall deformations occurring on the foundations. Soil–structure interactions relate to the influence of soil characteristics beneath a structure on its seismic performance. The flexibility of the building’s foundation and its interaction with the neighbouring soil are key factors in shaping the building’s reaction during earthquakes. To adequately address soil–structure interaction effects, it is imperative to consider the flexibility of both the structure and the soil in the analysis [
37].
It is worth noting that the concentration of total deformations is primarily in the middle of the foundations when SSSI effects are disregarded.
Figure 6 illustrates the deformation of the foundation of A1 without SSSI effects, using SAFE software, and with SSSI effects, using both ANSYS and SAFE software, due to the combination of gravitational loads. The difference between foundational deformations in terms of gravitational loads for A1 in
Figure 6a–c is evident. By adjusting the values of the subgrade reaction modulus in foundation design to account for SSSI effects, the foundation’s deformation in
Figure 6a converts into
Figure 6c, resembling
Figure 6b, which represents foundation deformations with SSSI effects. This shows that by changing the values of the subgrade reaction modulus, it is possible to apply SSSI effects in foundation design.
The maximum total deformation of the A1 foundation without SSSI effects, induced by the gravitational load, is situated on the right side of the building in the X-axis. However, it shifts to the left side of the foundation when influenced by the presence of A2. In fact, the maximum total deformation in the A1 foundation occurs in the corner of the foundation near the length of the other foundation without the SSSI effect (top right). Conversely, it is on the opposite side when considering the SSSI effect (top left). The total deformation with SSSI effects is 12 cm greater than that occurring in the case without SSSI effects in A1.
Figure 6a shows the foundation deformation of A1 by considering the same Ks for all elements of the foundation, while
Figure 6c shows it by considering the different values of Ks for each element of the foundation, paying attention to the deformation and pattern of deformation shown in
Figure 6b, which is related to foundation deformations that consider SSSI effects.
Figure 7 illustrates the total deformation of the A2 foundation under the influence of gravitational loads, showing three scenarios: (a) without SSSI effects, with structure–soil–structure interactions analysed (b) in ANSYS software, and (c) in SAFE software by varying values for Ks to account for SSSI effects. By adjusting the subgrade reaction modulus values beneath the foundation, one can manipulate the foundation’s deformation to achieve the desired results, taking into account SSSI effects in direct finite element modelling. With SSSI effects considered, the maximum deformation occurs on the left side (in the X-axis) of the A2 foundation, whereas without SSSI effects, it is located on the top left side. The maximum total deformation value decreases by 11 cm when SSSI effects are not applied. The differential settlements of A1 cause A2 to incline towards A1.
Figure 8 illustrates the deformation of the foundation of A1 under the influence of SSSI effects, considering the combination of gravitational and earthquake loads in the X direction using ANSYS and SAFE software (ANSYS software for analysing the structure and foundation and SAFE software for designing the foundation). It is evident that maximum deformation occurs on the right side of the foundation, adjacent to the other building’s foundation. In the context of this load combination, both buildings incline towards each other.
Finally,
Figure 9 shows the total deformation of the A2 foundation when subjected to a combination of gravitational and earthquake loads in the X direction, accounting for SSSI effects. The maximum deformation occurs in the areas of the foundations adjacent to each other, specifically on the right side of the A1 foundation and the left side of the A2 foundation, demonstrating that A2 tends to incline towards A1.
The magnitude and location of the maximum deformation change due to SSSI effects are shown above. These effects have led to a reduction between 10% and 15% in the weight of the reinforcement used (Wr) in the foundation’s design.
There have been several studies focusing on the utilization of varying values of the subgrade reaction modulus to account for soil–structure interaction effects in foundation design [
1,
2,
23,
24]. Employing different values for the subgrade reaction modulus in foundation design allows for the incorporation of SSI and SSSI effects, providing convenient tools for civil engineers in their simulations. For example, Equations (1) and (2),
Table 2, and
Figure 10 demonstrate a strong agreement with our results [
1,
2]:
where a denotes the coefficient of the storey, and b denotes the coefficient of different zones for a foundation, and they are calculated using Equation (2),
Figure 10, and
Table 2, respectively, with ‘n’ representing the number of storeys [
2].