Next Article in Journal
Dynamic Planning of Construction Site for Linear Projects
Next Article in Special Issue
Numerical Analysis of Twin Tunnels Lining under Different Seismic Conditions
Previous Article in Journal
Acknowledgment to Reviewers of Infrastructures in 2020
Previous Article in Special Issue
Seismic Assessment of Six Typologies of Existing RC Bridges
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of the Second Order Effect of the SSI on the Building during a Seismic Load

1
Laboratory of Civil Engineering and Environment, High School of Technologyof Salé Mohammed V University of Rabat, Salé 11060, Morocco
2
Mechanics and Applied Mathematics Department, Industrial and Seismic Engineering Research Team, National School of Applied Sciences of Oujda, Mohammed First University, Oujda 60000, Morocco
3
Doctoral School of Science and Technology, Faculty of Technology & Modeling Center, Lebanese University, Saida 1600, Lebanon
*
Author to whom correspondence should be addressed.
Submission received: 30 December 2020 / Revised: 23 January 2021 / Accepted: 26 January 2021 / Published: 29 January 2021
(This article belongs to the Special Issue Multi-Hazard Approach to Infrastructures Risk Reduction)

Abstract

:
The type and the properties of the soil can potentially intensify the internal forces on buildings during seismic loads. To predict the effects of the soil parameters on the soil–structure interaction of buildings, it is necessary to consider the soil–structure interaction (SSI) in the modeling process. Therefore, this document aims to evaluate the seismic effect on the maximal displacement and inter-story drift, and evaluate the behavior of buildings under the second-order effect known in the literature as the P-delta effect. For this purpose, three cases of buildings with 5, 10 and 15 stories were modelled using a FLAC 2D finite-difference element calculation software with infinite soil conditions, including five types of base with four types of soil (one cohesive soil and three non-cohesive soils) considering the soil–structure interaction and a fixed base (without soil–structure interaction). According to the results for the above-mentioned boundary, as the height of the building increases and due to the weak properties of the soil, we notice that the maximal displacements and inter-story drift increase considerably. To that purpose, we recommend considering the second-order effect in seismic design, especially for non-cohesive soil.

Graphical Abstract

1. Introduction

The latest research has shown that the interaction between soil and structure is a critical factor in the analysis of structures. Today we are moving away from considering a fixed base in earthquake engineering, as this method can cause great damage during an earthquake.
Computer and analytical approaches to solving soil–structure interaction problems have established the importance of SSI (soil–structure interaction) during dynamic loads. Since 1990, much effort has been made to understand the effects of soil on the structure during dynamic loading to improve the seismic design performance.
In 2017, an analysis on the impact of boundary conditions on the maximum displacement of the structure was carried out [1]. Values that more than doubled were found as the difference between the maximum displacements, taking into account the SSI and a fixed base. In 2015, an analysis of the impact of the maximum displacement was carried out [2], which led to more than 100% difference by increasing the height of the floors.
The soil–structure interaction has many effects on the structure, proved by a large amount of research that shows that SSI negatively influences the structure, and in other cases, the SSI can favorably affect the structure [3]. Moreover, an unsafe design can lead to the failure of the whole structure as shown by past earthquakes (e.g., the Kobe earthquake in 1995) [4].
It has been shown that the shear at the base of the flexible base is lower than the fixed base, while the inter-story drift of the structure increases significantly with the SSI effect [5]. The effect of SSI can increase efforts because of the second-order effects which are not taken into account in the structural design. Moreover, in the literature, a major part of the research is focusing on the analysis of the maximum displacement and the inter-story drift, and the second-order effort is practically neglected.
In a numerical analysis developed by Shehu et al. in 2019, with design parameters, the design moments due to P-delta effects vary from about 25 to 100%, which is very significant of the resistance capacity of the elements [6]. P-delta effects have only a small influence on the response. For buildings responding in a highly inelastic manner, P-delta effects can have a significant influence on response [7].
In the structural analysis of structures, displacements and strains are assumed to be small. This means that the geometrical characteristics of the structure do not change during loading, known as the first-order effect; however, it is accepted that large displacements and strains will occur when buildings are under large earthquake loadings, also due to the effect of unfavorable soil–structure interactions. Therefore, the inter-story drift may be large and the second-order effects will produce additional forces and moments in the structure. These second-order effects can reduce structural stiffness [7].
In this paper, an analysis of the impact of SSI was conducted to see the impact of the interaction of the soil structure, by changing the number of floors to 5, 8 and 10 floors and using four soil types, a fixed base and an acceleration spectrum with a maximum value of 0.24 g. For this purpose, several numerical models were used in the finite element calculation program FLAC 2D. In order to consider the dynamic soil–structure interactions in such a precise way, using the infinite elements of the soil to avoid the reflection of seismic waves, the results showed that the effects of the SSI generate values of moments due to P-delta effects that cannot be neglected.

2. Numerical Modeling

2.1. Model Description

In order to show the influence of the soil–structure interaction, a coupled model was created between structure, soil and infinite elements as shown in Figure 1. The soil was surrounded by the free field boundary on FLAD 2D 2005. The frame structural elements were idealized as 2D elastic beam elements for BEAM and COLUMN. The behavior of the reinforced concrete RC frame structure was assumed to be elastic and was modeled by using the parameters in the next section. For the meshing, a coarse mesh was used [8].
In numerical studies performed by different authors (e.g., [9,10,11,12]), concerning the 2D versus 3D modeling of the SSI, the difference between the results for the 2D and 3D models were negligible where the artificial model was adopted for rigid soil [12]. A large amount of other research [13,14,15,16,17] has applied a similar analysis.
In this research, we used three structural models, composed of 5-, 8- and 10-story models, and four selected soil types including one non-cohesive soil and three cohesive soils, represented as types 1, 2, 3 and 4.
In order to avoid the reflection of the seismic waves that generate an amplification of the stresses, free-field boundary conditions were used to take into account the continuity of the ground model as shown in Figure 1.
Soil characteristics of the soil during earthquakes are shown in Table 1.

2.2. Material Properties

Concerning the types of soil used, Table 1 shows the four types used in the different models with the following characteristics.
All types of soil can be described as follows:
Type 1: Dense permeable soils with a sandy-gravelly texture.
Type 2: High compactness soils with moderate cohesion of medium permeability and sand-loamy texture.
Type 3: Marls and marly limestone, massive brown limestone.
Type 4: Marly limestone and grey-brown marls.
These soil types represent the properties of the soils in critical regions in Morocco.
The properties of the structure are the elasticity modulus E = 30 GPa, the Poisson’s ratio n = 0.2, mass density ρ = 2400 kg/m3, yield stress б = 25 MPa. The bay length of the frame along the axis (ox) is taken to be 4.0 m, the story height is considered to be 3.0 m. The cross section of the beams is 30 × 30 cm, the column cross section is 30 × 30 cm.

2.3. Earthquake Description

In order to realize an appropriate study on the seismic response of the structural and artificial earthquake record, a maximum acceleration of 0.24 g was selected, which represents properties close to the riskiest zone in Morocco (Zone 3) [18] and which, for this reason, is used for the design of civil engineering structures.
The characteristics of the earthquake chosen for this study are presented in the Figure 2 and Figure 3 below.
The dynamic input is applied at internal gridpoints of the bottom of the model. Wave reflections at model boundaries are minimized by specifying free-field boundary conditions.

2.4. Damping

In this study, critical damping of 5% was taken into account according to the seismic code [18] which suggests 5% critical damping values for concrete structures and a damping interval for the construction that varies between 2 and 8% [19]. In addition, it was defined in the literature as an estimated value of the critical damping for the soil [20,21].
Roy and Dutta in 2010 studied the damping due to SSI on the response of the structure over a critical damping range of 2% to 30% [22], and it was noted that high values of damping do not give significant differences with critical damping of 5% [23].

2.5. Concept of the Second-Order Effect (P-delta Effect)

Basic parameters associated with the second-order effect are shown in Figure 4 and include a mass, m and a weight, P, associated with a rigid column with a flexible connection at the base. In order to take into account the viscous damping, a dumper is connected to the structure. Under the action of the lateral loads V, the overturning moment applied, without taking into account the effect P-delta equal to V.h. δ marks the horizontal displacement caused by the force V. The additional overturning moment induced by the effect of the second order will be generated by the normal force P and the distance δ, and the resulting moment M’ is given by the following expression:
M   =   V h   +   P δ
Consequently, a reduction in effective stiffness affects the structure from the inclusion of the P-Δ effect [24].

3. Results and Discussion

3.1. Lateral Deflection Results

The lateral deflection shown below represents the maximum at the last story level. The limit of the displacement is given by the code, which is about 0.4% of the overall height of the structure [24].
As shown in Figure 5, Figure 6 and Figure 7, the displacement increases with the decrease of the ground stiffness in the different cases and also with the increase of the height of the structure.
It can also be noticed that the non-cohesive soil (type 1) has practically the same value as the fixed base case.
Comparing the results of the maximum displacements of the fixed base and the other soil type models, we can observe that the effects of the SSI increase the displacements considerably between the 5-story and 8-story cases. It increases by 96% at the 5-story level for the soil type 4, which is less than between the 8-story and 10-story cases which increase by 37% at the 5-story level of soil type 4.
It is observable that the case of soil type 1 does not present any difference in displacement values with the fixed base for the three cases, which proves that the fixed base is only valid with a non-cohesive soil.
In addition, we can observe that in all cases, we exceeded the limit value of the code but it is amplified with the increasing story structure and the decreasing soil rigidity.

3.2. Inter-Story Drift Results

The inter-story drift mentioned below is calculated by the following:
D r i f t   =   d   ( i + 1 ) d   ( i )   h
The limitation of the code in our case is given by 1% [25].
As the figures show, the inter-story drift increases considerably between levels 1 and 2 (Figure 8, Figure 9 and Figure 10), with a small difference between the case of rigid soil (non-cohesive soil) and the embedded base case.
This indicates that the non-cohesive soil (type 1) has practically a similar behavior compared to the fixed base.
Comparing the results of the inter-story drifts between the fixed base and the other soil cases, the following results of the difference between the maximum and minimum values can be observed:
For the case of 5 stories, a difference of 0.4%; for the case of embedding and soil type 1, as well as the case of soil type 4, a difference of 1.3%.
For the case of 8 stories, a difference of 0.4%; for the embedded case and floor type 1, as well as the case of floor type 4, a difference of 1.9%.
For the case of 10 stories, a difference of 0.8%; for the embedded case and floor type 1, as well as the case of floor type 4, a difference of 2.3%.
The above results show that the inter-story drifts are amplified between the intermediate levels, which will generate second-order solicitations that are not inconsiderable.

3.3. Maximum Second-Order Moment Results

As shown in the figures, the values of the second-order moments increase considerably in the intermediate levels (Figure 11, Figure 12 and Figure 13), with identical values between the case of a non-cohesive floor (type 1) and the embedded base case, except for a small increase on the third floor for the 10 stories case. This indicates that the non-cohesive soil has behavior almost similar to a fixed base.
By comparing the results of the second-order moments between the fixed base and the other soil cases, the following results can be observed:
For the case of 5 stories, an increase of 233% between the fixed base case and the soil type 4 case, with a maximum value at the second-floor level.
For the case of 8 stories, an increase of 228% between the fixed base case and the soil type 4 case, with a maximum value at the third-floor level.
For the case of 10 stories, an increase of 250% between the fixed base case and the soil type 4 case, with a maximum value on the third-floor.
The above results show that the increases in second-order loads between the fixed base case and the soil type 4 case are proportional but present significant values when the height increases and the soil stiffness decreases.

4. Conclusions

In this paper, a soil–structure interaction model has been developed, using finite difference software (FLAC 2D), which is a powerful tool for this type of analysis.
Our purpose was to understand the second-order effects through a step-by-step analysis of displacements, according to the results presented previously. We find that the rigid non-cohesive soil type has almost identical properties to an embedding base during a seismic load. This result may somewhat relax structural engineers during their engineering studies, but the dynamic soil–structure interaction has dangerous effects on the increase in the number of floors and the decrease in soil stiffness.
We also find that the effect of the SSI amplified the second-order loads that proportionally change with the case the embedded base. This result can be taken into account during the engineering analysis in order to predict a structure more resistant to these effects.
It should be noted that, in this study, the effects of SSI in the case of rigid and non-cohesive soils is negligible, and on the other hand, for cohesive soil types, these effects cannot in any case be underestimated.

Author Contributions

Conceptualization, E.M.E., H.B. and A.E.O.; methodology, E.M.E. and A.E.O.; software, E.M.E. and A.E.O.; validation, M.R., M.C. and F.H.C.; formal analysis, E.M.E. and A.E.O.; investigation, E.M.E. and A.E.O.; resources, E.M.E. and A.E.O.; data curation, M.R. and M.C.; writing—original draft preparation, E.M.E.; writing—review and editing, E.M.E. and A.E.O.; visualization, E.M.E. and A.E.O.; supervision, M.R., M.C. and F.H.C.; project administration, M.R.; funding acquisition, M.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This paper and the research behind it would not have been possible without the exceptional support of my supervisor, Rougui Mohamed, Chourak Mimoun and Fadi Hage Chehade. there enthusiasm, knowledge and exacting attention to detail have been an inspiration and kept my work on track from my first to the final draft of this paper. Abdelday El Omari, for his help for the numerical modal and answered with unfailing patience numerous questions. Hasnae Boubel, my colleague at Mohamed V University, have also looked over my transcriptions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Edip, K.; Garevski, M.; Sheshov, V.; Bojadjieva, J. Boundary Effects in Simulation of Soil-Structure Interaction Problems. Soil Mech. Found. Eng. 2017, 54, 239–243. [Google Scholar] [CrossRef]
  2. Shanmugam, J.; Dode, P.A.; Chore, H.S. Analysis of Soil Structure Interaction in Framed Structure. Int. J. Comput. Appl. 2015, 975, 8887. [Google Scholar]
  3. Kwag, S.; Ju, B.; Jung, W. Beneficial and Detrimental Effects of Soil-Structure Interaction on Probabilistic Seismic Hazard and Risk of Nuclear Power Plant. Adv. Civ. Eng. 2018, 2018, 1–18. [Google Scholar] [CrossRef]
  4. Mylonakis, G.; Gazatas, G. Seismic soil-structure interaction: Beneficial or detrimental? J. Earthq. Eng. 2000, 4, 277–301. [Google Scholar] [CrossRef]
  5. Far, H. Dynamic behaviour of unbraced steel frames resting on soft ground. Steel Constr. 2019, 12, 135–140. [Google Scholar] [CrossRef]
  6. Saha, R.; Dutta, S.C.; Haldar, S.; Kumar, S. Effect of soil-pile raft-structure interaction on elastic and inelastic seismic behaviour. Structures 2020, 26, 378–395. [Google Scholar] [CrossRef]
  7. Shehu, R.; Angjeliu, G.; Bilgin, H. A Simple Approach for the Design of Ductile Earthquake-Resisting Frame Structures Counting for P-Delta Effect. Buildings 2019, 9, 216. [Google Scholar] [CrossRef] [Green Version]
  8. Montgomery, C.J. Influence of P-Delta effects on seismic design. Can. J. Civ. Eng. 1981, 8, 31–43. [Google Scholar] [CrossRef]
  9. Lysmer, J.; Kuhlemeyer, R.L. Finite dynamic model for infinite media. J. Eng. Mech. Div. ASCE 1969, 95, 859–877. [Google Scholar]
  10. Dutta, C.H.; Roy, R. A critical review on idealization and modelling for interaction among soil–foundation–structure system. Comput. Struct. 2002, 80, 1579–1594. [Google Scholar] [CrossRef]
  11. Zhang, X.; Wegner, J.L.; Haddow, J.B. Three-dimension dynamic soil–structure interaction analysis in the time domain. Earthq. Eng. Struct. Dyn. 1999, 28, 1501–1524. [Google Scholar] [CrossRef]
  12. Seo, C.G.; Yuna, C.B.; Kimb, J.M. Three-dimensional frequency-dependent infinite elements for soil–structure interaction. Eng. Struct. 2007, 29, 3106–3120. [Google Scholar] [CrossRef]
  13. Zheng, J.; Takeda, T. Effects of soil-structure interaction on seismic response of PC cable-stayed bridge. Soil Dyn. Earthq. Eng. 1995, 14, 427–437. [Google Scholar] [CrossRef]
  14. Galal, K.; Naimi, M. Effect of conditions on the response of reinforced concrete tall structures to near fault earth-quakes. Struct. Des. Tall Spec. Build. 2008, 17, 541–562. [Google Scholar] [CrossRef]
  15. Tabatabaiefar, S.H.R.; Fatahi, B.; Samali, B. An empirical relationship to determine lateral seismic response of mid-rise building frames under influence of soil-structure interaction. Struct. Des. Tall Spec. Build. 2014, 23, 526–548. [Google Scholar] [CrossRef]
  16. Far, H.; Fatahi, B.; Samali, B. Numerical and Experimental Investigations on Seismic Response of Building Frames under Influence of Soil-Structure Interaction. Adv. Struct. Eng. 2014, 17, 109–130. [Google Scholar]
  17. Far, H.; Fatahi, B. Idealisation of soil–structure system to determine inelastic seismic response of mid-rise building frames. Soil Dyn. Earthq. Eng. 2014, 66, 339–351. [Google Scholar] [CrossRef]
  18. Moroccan Earthquake Regulations; RPS 2000; Ministry of Housing and the Development of Space: Rabat, Morocco, 2000.
  19. Veletsos, A.S.; Meek, J.W. Dynamic behavior of building—Foundation systems. Earthq. Eng. Struct. Dyn. 1974, 3, 121–138. [Google Scholar] [CrossRef]
  20. Velez, A.; Gazetas, G.; Krishnan, R. Lateral Dynamic Response of Constrained-Head Piles. J. Geotech. Eng. 1983, 109, 1063–1081. [Google Scholar] [CrossRef]
  21. Gazetas, G. Formulas and Charts for Impedances of Surface and Embedded Foundations. J. Geotech. Eng. 1991, 117, 1363–1381. [Google Scholar] [CrossRef]
  22. Roy, R.; Dutta, S.C. Inelastic seismic demand of low-rise buildings with soil-flexibility. Int. J. Non Linear Mech. 2010, 45, 419–432. [Google Scholar] [CrossRef]
  23. Fatahi, B.; Far, H.; Samali, B. Soil-structure interaction vs. Site effect for seismic design of tall buildings on soft soil. Geomech. Eng. 2014, 6, 293–320. [Google Scholar] [CrossRef]
  24. Shehu, R. The p-δ-ductility effect: Overview the effectof the second order in the ductil structures. In Proceedings of the International Scientific Forum, Tirana, Albania, 12–14 December 2013; Volume 3, pp. 143–155. [Google Scholar]
  25. Moroccan Earthquake Code; RPS 2000; Ministry of Housing and the Development of Space: Rabat, Morocco, 2011.
Figure 1. The modeling concept in the FLAC 2D program.
Figure 1. The modeling concept in the FLAC 2D program.
Infrastructures 06 00020 g001
Figure 2. Time history of the artificial earthquake.
Figure 2. Time history of the artificial earthquake.
Infrastructures 06 00020 g002
Figure 3. The response spectrum of acceleration.
Figure 3. The response spectrum of acceleration.
Infrastructures 06 00020 g003
Figure 4. Concept of P-delta effect: example of single degree of freedom.
Figure 4. Concept of P-delta effect: example of single degree of freedom.
Infrastructures 06 00020 g004
Figure 5. Maximum lateral deflection case of a 5-story building.
Figure 5. Maximum lateral deflection case of a 5-story building.
Infrastructures 06 00020 g005
Figure 6. Maximum lateral deflection case of an 8-story building.
Figure 6. Maximum lateral deflection case of an 8-story building.
Infrastructures 06 00020 g006
Figure 7. Maximum lateral deflection case of a 10-story building.
Figure 7. Maximum lateral deflection case of a 10-story building.
Infrastructures 06 00020 g007
Figure 8. Inter-story drift case of a 5-story building.
Figure 8. Inter-story drift case of a 5-story building.
Infrastructures 06 00020 g008
Figure 9. Inter-story drift case of an 8-story building.
Figure 9. Inter-story drift case of an 8-story building.
Infrastructures 06 00020 g009
Figure 10. Inter-story drift case of a 10-story building.
Figure 10. Inter-story drift case of a 10-story building.
Infrastructures 06 00020 g010
Figure 11. Maximum second-order moment case of a 5-story building.
Figure 11. Maximum second-order moment case of a 5-story building.
Infrastructures 06 00020 g011
Figure 12. Maximum second-order moment case of an 8-story building.
Figure 12. Maximum second-order moment case of an 8-story building.
Infrastructures 06 00020 g012
Figure 13. Maximum second-order moment case of a 10-story building.
Figure 13. Maximum second-order moment case of a 10-story building.
Infrastructures 06 00020 g013
Table 1. Geotechnical properties.
Table 1. Geotechnical properties.
SoilRegionYoung’s Modulus E (GPa)Poisson’s Ratio (ν)Cohesion c (kPa)Friction Angle φ(◦)Dilation Angle ψ(◦)Unit Weight γ (KN/m3)
Type 1Casablanca 168.250.25035521
Type 2El hoceima28.350.351020019
Type 3Rabat240.415025018
Type 4Agadir100.2418020018
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Echebba, E.M.; Boubel, H.; El Omari, A.; Rougui, M.; Chourak, M.; Chehade, F.H. Analysis of the Second Order Effect of the SSI on the Building during a Seismic Load. Infrastructures 2021, 6, 20. https://0-doi-org.brum.beds.ac.uk/10.3390/infrastructures6020020

AMA Style

Echebba EM, Boubel H, El Omari A, Rougui M, Chourak M, Chehade FH. Analysis of the Second Order Effect of the SSI on the Building during a Seismic Load. Infrastructures. 2021; 6(2):20. https://0-doi-org.brum.beds.ac.uk/10.3390/infrastructures6020020

Chicago/Turabian Style

Echebba, El Mehdi, Hasnae Boubel, Abdelhay El Omari, Mohamed Rougui, Mimoun Chourak, and Fadi Hage Chehade. 2021. "Analysis of the Second Order Effect of the SSI on the Building during a Seismic Load" Infrastructures 6, no. 2: 20. https://0-doi-org.brum.beds.ac.uk/10.3390/infrastructures6020020

Article Metrics

Back to TopTop