## 1 Introduction

Observations based on the effects of previous strong earthquakes on structural damage have concluded that the seismic response of many engineering structures could be affected by deformability of the underlying foundation. Currently, structures become higher and larger and their seismic response may be affected by the medium on which they are founded. In addition, the response of a structure under seismic excitation is affected by the interaction between three connected structural parts: the structure, the foundation and the soil underlying and surrounding the foundation. Furthermore, the interaction between soil, foundation and structure sub-systems cannot be disregarded because the problem becomes increasingly important especially for structures placed on relatively soft soil in seismically active areas. The effects of soil-structure interaction (SSI) on the seismic response of structural systems are significant when the shear wave velocity of the supporting soil is less than 600 m/s.^{[1,2,3]} Hence, the effect of SSI should be taken into consideration in order to ensure SSI safety and economic efficiency of the structure during an earthquake.^{[3]}

On the other hand, the capacity of the soil to support various kinds of loadings depends on the strength and stability of the supporting soil materials. Consequently, not only the behaviour and responses of soils and foundations are important, but how soils and foundations affect seismic wave propagations and the subsequent loading transferred to main structures are also of great interest.^{[4]}

Various methods have been developed to implement SSI during the last fifty years and the problems can be solved by two main methods: the direct method and the substructure method. Current design codes provide guidelines to include its effects in seismic design of structure^{[5,6,7]} analysing it by the substructure method.

Furthermore, the soil exhibits spatial heterogeneities, which are among the main sources of structural damage and disfunctioning of the built systems.^{[8]} This heterogeneity or variability resulting from its deposition history is in many cases imperfectly known and large errors might appear during the identification of the parameters of the soil model. So, it is necessary to use probabilistic approaches to provide reliable design of the structure. In fact, the spatial variability of soils is a major source of SSI problems and the dynamic SSI problem should be modelled into a stochastic from word.^{[9,10,11]} However, to what extent can uncertainty in soil parameters affect the structure response to address the SSI problem?

The response of a structure to a seismic excitation might increase or decrease depending on the characteristics of the ground motion and dynamic properties of the structure and the underlying soil.^{[3,12]} Under earthquake excitation, the strain would invariably be larger than that measured during tests where the soil non-linearity plays an important role in the dynamic response of structures and soils, since each seismic record has different levels of shear stress in soil deposits.^{[13]} Even for a moderate earthquake, the strain range increases; degradation in soil stiffness becomes significant and has a major contribution to the overall response. It is obvious that strain induced in soil will depend upon dynamic loading, geological conditions of a site, stress history of soil and a number of other factors.^{[4]}

Widely, in order to evaluate the seismic demand of structures, it is assumed that the structure is founded on a rigid soil. However, real soil profiles are featured by complex shapes and nonlinear behaviour. Adequate information on the dynamic soil properties, in particular dynamic shear modulus and damping coefficient, is essential for accurate calculation of soil response and better analysis of SSI problems. Both shear modulus and damping coefficient are the most important material properties to characterize the dynamic behaviour of soils^{[14]} and both of these properties are affected by shear strain amplitude, effective stress level, void ratio and number of loading cycles in the case of clean sands. However, in the case of clay, they are affected by the number of loading cycles, over-consolidation ratio and plasticity index.^{[15]} Consequently, the shear modulus and damping coefficient of a soil medium exhibiting nonlinear behaviour should be sufficiently reflected in SSI analysis.^{[2,3]}

Over and above, in the dynamic analysis of SSI, soils are represented by linear, equivalent linear or nonlinear models. In linear analysis, the shear modulus and damping ratio are assumed constant. However, the relationship between cyclic secant shear modulus and cyclic shear strain amplitude is widely used in design practice to evaluate the seismic response and site effects through subsurface soils.^{[16]} According to EC8-2004,_{[5]} if the ground acceleration is equal to/or greater than 0.1 g, shear modulus must be multiplied by an average reduction factor. So, to what extent can correction in soil parameters (*G* and *ξ*) due to seismic event affect the structure response to address the SSI problem?

To answer the questions raised above, a nonlinear and stochastic analysis of SSI effects for shallow foundations is carried out in the present study. The main emphasis is on the effects of the correction of shear modulus and damping ratio values due to the nonlinear soil behaviour under seismic excitation on the response of the SSI system in terms of impedance functions and structure displacements. The soil under the structure is considered homogenous and is modelled as a discrete model based on the concept of cone model. The model cannot only give the impedance of the soil but can also model the soil responses.^{[10,11]} Furthermore, the response spectrum being a very important tool for seismic analysis or structural design and providing a very handy tool for engineers to quantify the demands of earthquake ground motion on the capacity of buildings to resist earthquakes, a comparative analysis is carried out on design response spectra and base shear forces for four seismic codes (Algerian Seismic Rules (RPA99 -2003),^{[17]} Eurocode 8 (EC8-2004),[5] International Building Code (IBC-2015)^{[18]} and the Indian Seismic Code (IS-1893-2002)^{[19]} considering three cases: the initial case of damping ratio or fixed base, the case of uncertain shear modulus of soil and the case of total foundation damping for SSI (before and after the correction of shear modulus and damping ratio).

## 2 Nonlinear and stochastic formulation of SSI problem

### 2.1 Simulation of the nonlinear behaviour

The nonlinear stress-strain behaviour of a soil can be expressed in two ways, with an increase in the shear strain amplitude: (i) degradation of shear modulus *G* and (ii) an increase in the damping coefficient *ξ*.

The soil nonlinearity can be taken into account by many models.^{[20,21,22,23,24,25]} In the equivalent linear method, the nonlinear stress-strain behaviour of a soil is approximated by a secant stiffness *G*_{sec} and an equivalent damping *ξ*_{eq} that are compatible with the strain in the soil induced by the ground shaking. Okada et al.^{[26]} studied the seismic pile response of a structure-pile-soil system to evaluate the strain-dependent nonlinearity of the ground using an equivalent linearization method.

In the present study, the shear modulus and damping coefficient are modified using a linearized iterative method based on design response spectra (Sa) for four seismic codes (RPA99-2003,^{[17]} EC8-2004,[5] IBC-2015[18] and IS-1893-2002^{[19]}). All the steps of the method are summarized in **Appendix 1** according to Chowdhury and Dasgupta.^{[27]}

### 2.2 Stochastic based formulation

The soil properties can be delineated using deterministic or probabilistic model. They greatly depend on the soil deposition conditions and loading history, and may be derived from a common set of in situ or laboratory test data.^{[28,29]} In other words, soil is nonhomogeneous and its formation is an outcome of a variety of random geological processes,^{[30]} and consequently, the stochasticity of soil properties remains largely unknown. Several methods are used to treat stochastic problems where Monte-Carlo simulation (MCS)^{[31,10]} provide a robust conceptually simple way to account rationally for various uncertainties.^{[32,33,34]} The most used probability distribution function to estimate uncertainty in geotechnical properties for performance based earthquake engineering are uniform, normal, lognormal, gamma and exponential.^{[35]} The lognormal distribution is commonly used to model certain types of data that appear in several fields of engineering. Since this includes most, if not all, engineering (mechanical) systems, particular properties of the lognormal random variable (such as non-negativeness and skewness) and of the lognormal hazard function (which increases initially and then decreases) make lognormal distribution a suitable fit for some engineering data sets. However, the lognormal distribution can have widespread application; a generalized lognormal distribution can be used to provide better fits for many types of experimental or observational data.^{[36]}

On other hand, structural engineering design is replete with uncertainties, some of which are obvious and some of which many engineers may never have considered as stated by Bulleit.^{[37]} Otherwise, most shallow foundations supporting structures are generally not circular in shape. So equivalent circular foundations, where foundation radii are computed separately for translational (r_{a}) and rotational (r_{m}) deformation modes,^{[38]} are used.

In addition, the effective height of the structure (h) shall be taken as 0.7 times the structure height for multi-story structures, although in reality, there are irregular structures in plan and/or elevation (Fig. 1a). Thus, the equivalent foundation radius and the effective height of the structure may be assumed subject to uncertainties and may be considered as random variables.^{[11]}

In this study, the effects of soil nonlinearity by correction of the shear modulus and the damping ratio together with soil stochasticity are considered. Several thousands of random draws are generated for soil according to the lognormal distribution with a coefficient of variation (COV) for G taken as 20% according to the literature.^{[39,40,41,42]}

*μ*

_{ln,x}and

*σ*

_{ln,x}are the lognormal mean value and standard deviation, respectively.

## 3 Modelling of the SSI problem

The schematic view of a SSI system and the equivalent system of single degree of freedom (SDOF) are shown in Figs. 1a and 1b, respectively. The superstructure is attached to the ground by springs-dampers elements, which are^{[43]} able to move in two directions to give translational and rotational motions. The structure is described by its stiffness *k*, mass *m*, equivalent height *h* and damping *c* (Fig. 1b). The SSI problem is solved in the frequency domain using the Cone model, which is a simple model for foundation vibration analysis based on the wave propagation.^{[11,44]} Also, cone model is well suited for evaluating the impedance functions.^{[42,45,46,47]}

SSI has both kinematic and inertia interaction effects on a structure, but for shallow foundation, only the inertia effect is considered^{[10]} to calculate the dynamic response of a structure supported on a shallow foundation. The soil-foundation interface is modelled by discrete-elements including horizontal and rocking equivalent linear springs and viscous dashpots with frequency-independent coefficients. The parameters of this interface element (impedance functions) are obtained by the simple physical cone model. The cone model is simple one-dimensional model for foundation vibration analysis, which consists of replacing the soil deposit, for each degree of freedom of the foundation by a truncated semi-infinite elastic cone. The cone model was originally developed by Ehlers (1942) to represent a surface disk under translational motions, and since that time, the model shows multi uses and development in dynamic soil structure interaction problems.^{[41,42,43,45,46,47]} The concept of cone model has been rearranged and extended by Wolf^{[48,49,50,51,52]} to include a complete range of dynamic excitations and physical situations.^{[53]}

Most of the published results using cone model are confined to the determination of the dynamic response of the foundation in the form of complex impedance functions, which represent the dynamic stiffness of the soil medium surrounding the foundation. The real part of the impedance function, denoted dynamic stiffness, reflects the stiffness and inertia of the supporting soil, and its dependency on frequency is solely attributed to the influence of frequency on inertia. While, the imaginary part represents the energy dissipation in the system generated as a result of the wave propagation away from the foundation (radiation damping).

To study the impedance function of a shallow foundation resting on the surface of a soil layer underlain by a half-space, the foundation was replaced by a rigid massless foundation of radius *r*_{0} as illustrate in Fig. 2. The layer with depth *d* has a shear modulus *G*_{1}, Poisson's ratio *ν*_{1} , mass density *ρ*_{1} and hysteretic damping ratio *ξ*_{1}. The corresponding parameters of the half-space are *G*_{2}, *ν*_{2}, *ρ*_{2}, and *ξ*_{2}.

*a*

_{0}= ω

*r*

_{0}/ν

_{s}, where: ν

_{s}is the shear wave velocity of soil) as:

*K*(

*K*

_{h}for horizontal motion and

*K*

_{r}for rocking motion), are as follows:

*a*

_{0}) and c(

*a*

_{0}) are the dimensionless dynamic stiffness and damping factors.

*u*

_{0}

*) and the total displacement of the SSI system (*

^{t}*u*

^{t}) (Fig. 1(b)) are expressed, respectively, by

*u*

_{g}is the ground displacement.

^{[10,11,44]}

*S*

_{ξx}(

*ω*) and

*S*

_{θx}(

*ω*) are horizontal and rocking components of the impedance functions, respectively.

*ω/ω*

_{s}is the frequency ratio and

**is the complex number.**

*i*A FORTRAN computer program is realized in order to assess the effects of the stochasticity (shear modulus G with COV =20%) and the nonlinear behaviour of soil (corrected shear modulus and damping ratio) on the foundation and superstructure responses.

The steps followed to realize the program are shown in **Appendix II**.

## 4 Results and Discussions

The correction of dynamic soil characteristics (*G* and *ξ*) is performed based on an iteration method,^{[27]} which uses response spectra as explained in **Appendix I**.

The mean values of soil parameters before the correction are selected such that *G*_{1}/*G*_{2} = 0.544 and the damping ratio *ξ*_{1} = *ξ*_{2} = 5%. The other parameters are mass density ratio *ρ*_{1}/*ρ*_{2} = 0.85, Poisson's ratio *ν*_{1} = *ν*_{2} = 0.25 and depth to radius ratio d/*r*_{0} = 1. Two values of the coefficient of variation (COV) for shear modulus parameters are assumed here (10% and 20%). Then, by using the linear iterative method^{[27]} and RPA99-2003^{[17]} code design spectrum, the soil characteristics (*G*_{1}, *ξ*_{1}) are corrected according to **Appendix I** and subsequently: *G*_{1}/*G*_{2} = 0.309 and *ξ*_{1}/*ξ*_{2} = 1.32.

### 4.1 Effects of the correction of shear modulus and damping coefficient on SSI response

The impedance functions for horizontal and rotational motions are determined for a circular foundation of radius *r*_{0} on the surface of a soil layer resting on a half-space.

The results in terms of impedance functions (written in terms of real and imaginary parts versus the dimensionless frequency) for horizontal and rotational vibration modes are shown in Fig. 3. This figure shows the influence of the correction of shear modulus and damping coefficient on the real and imaginary part of impedance functions (spring and damping coefficients). It appears from this figure that the correction of the shear modulus and the damping coefficient due to the imposed seismic motion is very remarkable. Both real and imaginary parts of impedance functions changing as *G* and *ξ* values are corrected.

The change is apparent for all the frequencies from 0 to 8. This result is consistent with the results of other researchers like those of Lutes et al.,^{[39]} when they used a non-classical modal-analysis-based formulation to quantify the variation of stochastic response in an SSI system with both uncertain soil properties in the soil–foundation system and uncertain structural properties in the superstructure on structural response for seismically excited soil-structure interacting (SSI) systems, those of Mittal et al.,^{[54]} when carrying out a parametric study on the effect of seismic zone on structural response of tall chimney by incorporating strain dependent shear modulus and determining the effect of variability in soil type and those of Zafarkhan and Dehkordi^{[55]} while presenting the effect of variability in soil type (different shear modulus) and structure height on soil-structure system responses.

The displacements of the structure *u* (*ω*) and the mass with respect to the free field (*u*_{0} (*ω*) + *hθ*_{0} (*ω*) + *u* (*ω*)) and the total displacement of the base (*u*_{g} (*ω*) + *u*_{0} (*ω*)) are plotted versus the excitation frequency normalized by that of the structure *ω*/*ω*_{s} in Fig 4. The displacements before and after the correction of G and *ξ* are compared.

The displacement relative to the free field of the structure shows a decrease of about 33%, while the mass displacement and the total displacement of the base increase in absolute values by about 4% and 6%, respectively. But on average, all the displacements show the same dimensionless frequency shift of 18%. Fattah et al.^{[56]} studied the effect of sand density and whether the sand is dry or saturated on its dynamic response for the foundation-soil system. They found that for dry and saturated conditions, the maximum amplitude of displacement decreases with increasing the relative density of sand and contact area of footing and increase with increasing the amplitude of loading. The maximum displacement amplitude response of the foundation resting on dry sand models is more than that on the saturated sand by about 5.0–10%.

This variation explains that the structure becomes more rigid and the soil more flexible and shows the interest of taking into account the nonlinear soil behaviour in SSI analysis. The figure clearly shows that SSI alters the frequency content of the structure response. Also, the figure reported herein shows that SSI is more significant when soil nonlinearity case will take place during earthquake compared to the linear case (G and ξ are constants), where the nonlinearity reduces the stiffness of the soil (reduce shear modulus) and increases the hysteretic damping. Similar observation has been made by Li et al.^{[57]} when studying the effect of different soil properties on structure response, and Farghaly and Ahmed^{[2]} in the case of seismic loading. Fattah et al.[58] found that the displacement amplitude of vibration response of foundation system placed on the surface of dry dense sand models is less than that of dry loose sand models; according to the authors, these results are attributed to the increase in the stiffness and the modulus of elasticity of dense sandy soil that makes the soil stiffer and resist vibrations. In 2017, Fattah et al.^{[59]} found that the amplitude of displacement of a single pile embedded in saturated sandy soil due to vertical vibration is more than that of the dry soil, and Fattah et al.^{[60]} explained that this trend is due to the increase in the pore water pressure during dynamic load that causes reduction in the inter-particle forces between the solid particles of the soil skeleton, hence causing an increase in displacement response.

### 4.2 Effects of stochastic shear modulus and correction of soil characteristics on the SSI response

In this subsection, the stochastic and the nonlinear effects are considered by means of uncertainty of G and correction of *G* and *ξ*, respectively, on the response of the SSI system. A lognormal distribution is used to generate random numbers of shear modulus using Monte Carlo Simulations for coefficient of variation of 20%.

Fig. 5 displays the stochastic as well as the nonlinear displacements of the structure, the mass and the base due to the random variations of shear modulus and the correction of shear modulus and damping, respectively. Fig. 5 shows that the random variations of *G* greatly affect the response of the structure. It is observed that 20% COV of *G* reduces the peak structure displacement (16.6%) by half that is reduced by the correction of *G* and *ξ* (33.1%).

Additionally, the frequency is shifted about 18% due to the correction of *G* and *ξ*, which means that the soil becomes more rigid. However, Fig. 5b indicates that the peak value of mass displacement is reduced by 17.4% due to 20% COV of G compared to the nonlinear as well as to the deterministic ones (initial G). For the displacement of the base (Fig. 5c), for random variations of the shear modulus (COV = 20%), the displacement remains more or less unchanged compared to the initial case but it reduces by 10.4% due to nonlinear effects. Table 1 summarizes the peak SSI responses and the error margin between stochastic and deterministic response and between nonlinear and deterministic (initial G) response.

Error margin (%) in displacements and frequency SSI system due to random G and corrected (G and ξ).

SSI response | Deterministic | Stochastic (20% COV of G) | Nonlinear | Error margin (%) | |
---|---|---|---|---|---|

20% COV of G vs. initial G | Corrected (G & ξ) vs initial G | ||||

|u|/|ug| | 1.81 | 1.51 | 1.21 | − 16.6% | − 33.1% |

|u_{0}+hθ_{0}+u|/|ug | 7.35 | 6.07 | 7.37 | − 17.4% | + 0.3% |

|u_{0}+ug|/|ug| | 0.86 | 0.86 | 0.95 | ±00.0% | + 10.4% |

Frequency (ω/ω_{s}) | 0.5 | 0.5 | 0.41 | ±00.0% | −18.0% |

It can be seen from Table 1 that, for corrected *G* and *ξ*, the frequency ratio (*ω*/*ω*_{s}) goes on decreasing about 18% or time period goes on increasing for the three displacements (Figs. 5a to 5c). This trend means that soil becomes more flexible when shear degradation and damping variation are taken into account. Also, when the shear modulus goes on decreasing and damping coefficient increasing the displacements of structure, masse and the base reduce by 33%, 3% and 10%, respectively, compared to SSI responses for deterministic (initial *G* and *ξ*).

### 4.3 Effects of the correction of the foundation damping on the SSI response

Usually, the seismic calculation of the structure is performed using the response spectrum with a well-fixed damping ratio. However, in reality, there is not only the soil damping but there is also the foundation damping that includes both effects of energy loss from waves propagating away from a vibrating foundation (radiation damping) and hysteretic action in supporting soil (material damping).

Accordingly, the effects of the foundation damping of flexible base on the design spectrum (Fig. 6) of four design codes (RPA99-2003,^{[17]} EC8-2004,[5] IBC-2015[18] and IS-1893-2002^{[19]}) are studied. Furthermore, to further enrich the information, the base shear force of structure is also studied (Fig. 7). The design spectrum and the base shear values are evaluated according to the four cases for each code: (i) fixed base, (ii) SSI effects (initial G and x), (iii) nonlinear SSI (corrected G and x), and (iv) stochastic SSI (random G with COV = 20%). Table 2 in **Appendix III** summaries the ordinates of design spectra and shear base force for the four codes.

^{[61]}

_{0}is contributed from both the structural viscous damping ξ and the foundation damping ξ

_{f}consisting generally of radiation and material damping components.

*T*is the fixed base period and

_{rd}of the soil-structure system is calculated based on a combination of translational damping ξ

_{y}and rotational damping ξ

_{xx}:

Figs. 6 and 7 show the reduction in design spectrum and base shear, respectively, due to SSI using Eq. (4) for the three cases (ii, iii and iv) mentioned above, where they are compared to the fixed base case for the four codes.

It appears from Figs. 6a–6d (red line) that when including SSI effects with the foundation damping (flexible base case) of the system, a significant decrease appears in design spectral values. This decrease is of 9%, 10%, 6% and 2% for RPA99-2003,^{[17]} EC8-2004,[5] IBC-2015[18] and IS-2002^{[19]} seismic codes, respectively, compared to the fixed base (black line).

Besides, Figs. 6a–6d (green line) indicate that the same structure responds more close to the previous case when the shear modulus is taken as random variable and only the IBC-2015^{[18]} code shows a remarkable decrease. The design spectral values decrease, compared to the fixed base (black line) for the four codes, by 7%, 9%, 13% and 5%, respectively.

In the last instance, Figs. 6a–6d (blue line) show the effects of the correction of both shear modulus and damping coefficient of the soil on the design response spectra where the correction effect clearly appears by a huge decrease observed in the responses compared to the fixed base and especially for the three codes EC8-2004,^{[5]} IBC-2015^{[18]} and IS-1893.[19] The observed reductions in the response spectra for the four codes are 18%, 36%, 37% and 43%, respectively.

On the other hand, Fig. 7 shows the effect of SSI on the base shear for the four codes for the three cases included: SSI effects (initial G and x) in red line, nonlinear SSI (corrected G and x) in blue line, and stochastic SSI (random G with COV = 20%) in green line compared to the fixed base in black line. It is observed that for SSI analysis including flexible base damping, the base shear forces of the structure are noticeably decreased by 9%, 3%, 5% and 4.5%, respectively. In the stochastic case, the results are very close to the previous case results, but only the IBC-2015^{[18]} code shows remarkable decrease, as 7%, 3%, 11% and 7.5%, respectively, compared to the fixed base for the four codes.

Lastly, the effect of the correction of G and ξ of soil on the base shear force for all codes’ cases, where remarkable decreases of 18%, 9%, 47% and 46%, respectively, are displayed. For the two last codes, results are similar to those obtained by Mittal and Gajinkar^{[62]} when studying the response of chimney by incorporating strain dependent shear modulus, where they find that the design shear force reduces between 31% to 47% when the soil become more flexible. Whereas Worku,^{[63]} when studying the importance of inertial SSI on the design spectral and on the shear base for EC8^{[5]} and NEHRP-2003,^{[61]} found a reduction between 7% to 39% for two different flexible soil types and four different varying structure heights.

On other hand, Jayalekshmi and Chinmayi^{[64]} showed that the design base shear forces obtained as per conventional design practice are higher compared to SSI values and decrease with an increase in foundation stiffness, based on a comparative study on seismic provisions of IS-1893^{[19]} and IBC-2006^{[65]} codes for a multi-storey reinforced concrete framed buildings. Nevertheless, it is important to note that some of the SSI codes cap the maximum allowed base shear reduction to 30%.

## 5 Conclusions

The present study demonstrated the importance of taking into account of inertial SSI. A comparative analysis on design response spectra and base shear forces was carried out for four seismic codes (Algerian Seismic Rules RPA99 – 2003,^{[17]} EC8-2004,[5] International Building Code IBC-2015^{[18]} and the Indian code IS-1893-2002[19]) considering the three cases of SSI: (*i*) SSI effects (initial G and x), (*ii*) nonlinear SSI (corrected G and x), and (*iii*) stochastic SSI (random G with COV = 20%).

Firstly, the great influence of the corrected shear modulus and damping coefficient of soil on the impedance functions and the response of structure in the case of SSI analysis was shown. Results showed a shift in the frequency of flexible base by 18% and a decrease in the response of the structure.

Secondly, the influence of the modification of damping of the SSI system for the three cases (SSI effects – initial G and x), nonlinear SSI (corrected G and x) and stochastic SSI (random G with COV = 20%) compared to the fixed base on the design response spectra and shear base forces for the four codes was shown schematically in Figs. 6 and 7. From the results and figures, it was seen that design response spectra and base shear forces for the four codes (RPA99-2003,^{[17]} EC8-2004,[5] IBC-2015[18] and IS-1893-2002^{[19]}) are very sensitive to foundation damping of the flexible base.

A comparative study between seismic codes has been conducted to demonstrate the similarities and dissimilarities between them. The IBC-2015^{[18]} and IS-1893^{[19]} codes show considerable decreased values in design spectra and shear base force compared to the other codes. These results agree with the results in literature, where the shearing strains due to earthquakes ranging between 0.01 and 0.5% reduce shear modulus by 0.9 and 0.2 G_{max}, respectively.

Hence, ignoring the SSI effects in seismic design could lead to large overestimation of design response spectra for buildings. This could lead to uneconomic designs, which can be avoided if SSI is accounted for. Furthermore, the effects of correction of shear modulus and damping ratio on the response spectra of SSI system is very important, and in this case study, the total system damping exceeds 20%. Moreover, results of such a study may encourage to include SSI effects in some seismic code like Algerian Seismic Rules.

## References

- [1]↑
Çelebi, E., Goktepe, F., Karahan, N. (2012). Non-linear finite element analysis for prediction of seismic response of buildings considering soil-structure interaction.

*Natural hazards and earth system science, 12*: 3495–3505. - [2]↑
Farghaly, A.A., Ahmed, H.H. (2013). Contribution of soil-structure interaction to seismic response of buildings.

*Geotechnical engineering, 17(5)*: 959–971. - [3]↑
Park, J.H., Choo, J.F., Cho, J.R. (2013). Dynamic soil-structure interaction analysis for complex soil profiles using unaligned mesh generation and nonlinear modeling approach.

*Structural engineering, 17(4)*: 753–762. - [4]↑
Jia, J. (2018). Soil dynamics and foundation modeling, offshore and earthquake engineering Springer nature, Bergen, Norway.

- [5]↑
EC8-2004 (English): Eurocode 8: Design of Structures for Earthquake Resistance. Part 5 : Foundation, retaining structures and geotechnical aspects.

- [6]↑
FEMA 440 (2005), Recommended Improvements of Nonlinear Static Seismic Analysis Procedures, Applied Technology Council: California.

- [7]↑
ASCE (2010): Applied Technology Council. ‘Tentative Provisions for the Development of Seismic Regulations for Buildings’ ATC-3-06: California.

- [8]↑
Breysse, D., La Borderie, C., Elachachi, S.M., Niandou, H. (2007). Spatial variation in soil properties and their influence on structural reliability.

*International Journal of Civil Engineering and Environmental Systems, 2(24)*: 73–86. - [9]↑
Cottereau, R., Clouteau, D., Soize, C. (2007). Probabilistic impedance of foundation: impact on the seismic design on uncertain soils.

*Earthquake engineering and structural dynamics, 1*: 1–17. - [10]↑
Guellil, M.E., Harichane, Z., Djilali, B., Sadouki, A. (2017). Soil and structure uncetrtainty effects on the soil foundation structure dynamic response.

*Earthquake and structure, 2*: 153–163. - [11]↑
Harichane, Z., Guellil, M.E., Gadouri, H. (2018). Benefits of probabilistic soil-foundation-structure interaction analysis.

*Inter Jour of Geotechnical earthquake engineering, 9(1)*: 43–64. - [12]↑
Pitilakisa, D., Dietz, M., Wood, D.M., Clouteau, C., Modaressi, A. (2008). Numerical simulation of dynamic soil–structure interaction in shaking table testing.

*Soil Dynamics and Earthquake Engineering, 28*: 453–467. - [13]↑
Van Nguyen, Q., Ftahi, B., Hokmabadi, A.S. (2017). Influence of Size and Load-Bearing Mechanism of Piles on Seismic Performance of Buildings Considering Soil–Pile Structure Interaction.

*International Journal of Geomechanics (ASCE*), DOI: . - [14]↑
Chehata, A., Harichane, Z., Karray, M. (2017). Non-linear soil modelling by correction of the hysteretic damping using a modified Iwan model together with Masing rules.

*InternatIonal Journal of Geotechnical Engineering*, 1–13. - [15]↑
Brennan, A.J. Thusyanthan, N.I. Madabhushi, S.P.J. (2005). Evaluation of shear modulus and damping in dynamic centrifuge tests. Journal of geotechnical and geoenvironmental engineering. DOI: .

- [16]↑
Kishida, T. (2016). Comparaison and correction of modulus reduction models for clays and silts.

*J. geotechnical and geoenvironemental Engineering*. - [17]↑
RPA99, version (2003) Règlement parasismique Algérienne. Central national de recherche appliquée en génie parasismique.

- [18]↑
International Building Code (IBC) (2015) Chapter 16 – Section 1613 ICC Central Regional Office, Country Club Hills (IL, USA).

- [19]↑
IS 1893-2002 Part-1 Criteria For Earthquake Resistant Deign of Structures Part-1 General Provisions And Buildings Fifth Revision.

- [20]↑
Zhang, L., Ahmari, S. (2011). Nonlinear analysis of laterally loaded rigid piles in cohesive soil. Int. J. Numer.

*Anal. Meth. Geomech*, DOI: - [21]↑
Ravikumar, C.R., Gunneswara, T.D.R. (2012). Study of soil interaction in a model building frame with plinth beam supported by pile group.

*International Journal of Advanced Structural Engineering, 4(11)*: 1–15. - [22]↑
Truty, A. (2018). On consistent nonlinear analysis of soil–structure interaction problems.

*Studia Geotechnica et Mechanica, 40(2)*: 86–95. - [23]↑
Haldar, S., Basu, D. (2016). Analysis of Beams on Heterogeneous and Nonlinear Soil. International Journal of Geomechanics.

*International Journal of Geomechanics*, DOI: . - [24]↑
Tsai, C.C., Liu in, H.W. (2017). Site response analysis of vertical ground motion in consideration of soil nonlinearity.

*Soil Dynamics and Earthquake Engineering, 102*: 124–136. - [25]↑
Jastrzębska, M., Lupieżowiec, M., Uliniarz, R., Jaroń, A. (2014). Analysis of the vibration propagation in the subsoil.

*Studia Geotechnica et Mechanica, 26(3)*: 9–19. DOI: . - [26]↑
Okada, T., Fujita, K., Takewaki, I. (2016). Robustness evaluation of seismic pile response considering uncertainty mechanism of soil properties.

*Innovative infrastructure solutions*, DOI . - [27]↑
Chowdhury, I., Dasgupta, S.P. (2009). Dynamics of structure and foundation – A Unified Approach, 2. Applications. Taylor & Francis Group, London, UK.

- [28]↑
Raychowdhury, P., Singh, P. (2012). Effect of nonlinear soil-structure interaction on seismic response of low-rise SMRF buildings.

*Earthquake Engineering And Engineering Vibration, 11(4)*: 541–551. DOI: . - [29]↑
Hu, Q., Li, H., Yang, G., Cai, Y. (2019). Effects of uncertainty of dynamic shear modulus ratio on design ground motion.

*Soil Mechanics and Foundation Engineering, 56(2)*: 82–90. DOI . - [30]↑
Uzielli, M., Lacasse, S., Nadim, F., Phoon, K.K. (2016). Soil variability analysis for geotechnical practice. Conference: Proceedings of the 2nd International Workshop on Characterisation and Engineering Properties of Natural Soils, At Singapore. DOI: .

- [31]↑
Manolis, G.D. (2002). Stochastic soil dynamic.

*International Journal of Soil Dynamics and Earthquake Engineering, 22(1)*: 3–15. - [32]↑
Cao, Z., Wang, Y., Li, D. (2016). Efficient Monte Carlo Simulation of Parameter Sensitivity in Probabilistic Slope Stability Analysis. Probabilistic Approaches for Geotechnical Site Characterization and Slope Stability Analysis, 169–184.

- [33]↑
Djilali Berkane, H., Harichane, Z., Guellil, M.E., Sadouki, A. (2018). Investigation of Soil Layers Stochasticity Effects on the SpatiallyVarying Seismic Response Spectra.

*Indian geotechnical journal.*https://doi.org/10.1007/s40098-018-0301-y. - [34]↑
Górska, K., Muszyński, Z., Rybak, J. (2012). Displacement monitoring and sensitivity analysis in the observational mathod.

*Studia Geotechnica et Mechanica, 35(3)*, 14(4): 25–43. DOI: . - [35]↑
Baecher, G.R., Christian, J.T. (2003). Reliability and statistics in geotechnical engineering. John Wiley & Sons Ltd, England.

- [36]↑
Martín, J. Pérez, C.J. (2008). Application of a generalized lognormal distribution to engineering data fitting. Proceedings of the European safety and reliability conference, Esrel 2008, and 17th SRA-Europe, Valencia, Spain, September, 22–25, 2008.

- [37]↑
Bulleit, W. (2008). Uncertainty in Structural Engineering.

*Practice Periodical on Structural Design and Construction, ASCE, 13(1)*: 24–30. DOI: . - [38]↑
Stewart, J.P., Kim, S., Bielak, J., Dobry, R., Power, M.S. (2003). Revisions to soil-structure interaction procedures in NEHRP design provisions.

*Earthquake Spectra, 19(3)*: 677–96. - [39]↑
Lutes, L.D., Sarkani, S., Jin, S. (2000). Response variability of an SSI system with uncertain structural and soil properties.

*Engineering Structures, 22(6)*: 605–620. - [40]↑
Badaoui, M., Berrah, M.K., Mébarki, A. (2009). Soil height randomness influence on seismic response: Case of an Algiers site.

*Computers and Geotechnics, 36*: 102–12. - [41]↑
Moghaddasi, M., Cubrinovski, M., Chase, J.G., Pampanin, S., Carr, A. (2011). Effects of soil-foundation-structure interaction on seismic structural response via robust Monte Carlo simulation.

*Engineering Structure, 33*: 1338–1347. - [42]↑
Mirzai, F., Mahsuli, M., Ghannad, M.A. (2016). Probabilistic analysis of soil-structure interaction effects on the seismic performance of structures.

*Earthquake Engineering & Structural Dynamics*, DOI: . - [43]↑
Kumar, M., Dass, Goel, M., Matsagar, V.A., Rao, K.S. (2014). Response of semi-buried structures subjected to multiple blast loading considering soil–structure interaction.

*Indian geotechnical journal, DOI*. - [44]↑
Wolf, J. P., Deeks, A. J. (2004). Foundation Vibration Analysis : a Strength of Materials Approach. Amsterdam: Elsevier.

- [45]↑
Pradhan, P.K., Baidya, D.K., Ghosh, D.P. (2004). Dynamic response of foundation resting on layered soil by cone model.

*Soil Dyn. Earthq. Eng, 24(6)*: 425–34. - [46]↑
Pradhan, P.K., Mandal, A., Baidya, D.K., Ghosh, D.P. (2008). Dynamic response of machine foundation on layered soil.

*Geotech. Geol. Eng, 26(4)*: 453–68. - [47]↑
Azarhoosh, Z., Ghodrati Amiri, G.R. (2010). Elastic Response of Soil-Structure Systems Subjected to Near-Fault Rupture Directivity Pulses. Soil Dynamics and Earthquake Engineering - Proceedings of Sessions of Geoshanghai 2010. Geotechnical special publication; 201: American Society of Civil Engineers 1801 Alexander Bell Drive Reston, Virginia, 20191-4400.

- [48]↑
Meek JW, Wolf JP. Cone models for homogeneous soil I. Geotechnical Engng Div, ASCE 1992; 118(5): 667–85.

- [49]↑
Meek JW, Wolf JP. Cone models for soil layer on rigid rock II. Geotechnical Engng Div, ASCE 1992; 118(5): 686–703.

- [50]↑
Wolf JP, Meek JW. Cone models for a soil layer on a flexible rock half-space. Earthquake Engng Struct Dyn 1993; 22: 185–93.

- [51]↑
Meek JW, Wolf JP. Cone models for a embedded foundation. Geotechnical Engng Div, ASCE 1994; 120(1): 60–80.

- [52]↑
Wolf JP, Meek JW. Dynamic stiffness of foundation on layered soil half-space using cone frustums. Earthquake Engng Struct Dyn 1994; 23: 1079–95.

- [53]↑
Mohasseb, S. Ghazanfari, N. Rostami, M. Rostami, S. (2020). Effect of Soil–Pile–Structure Interaction on Seismic Design of Tall and Massive Buildings Through Case Studies. Transportation Infrastructure Geotechnology 1: 13–45. https://doi.org/10.1007/s40515-019-00086-7

- [54]↑
Mittal, R.K., Rawat, A., Rawat, S. (2016). Soil Structure Interaction in Indian Seismic code: Recommendations for Inclusion of Potential Factors.

*International journal of earth sciences and engineering, 09(03)*: 124–130. - [55]↑
Zafarkhah, E., Dehkordi, M.R. (2017). Evaluation and numerical simulation of soil type effects on seismic soil-structure interaction response of RC structures.

*Journal of Vibroengineering, 19(7)*: 5208–5230. https://doi.org/10.21595/jve.2017.18286. - [56]↑
Fattah, M. Y., Al-Mosawi, M. J., Al-Ameri, A. F. I., (2017), ‘Dynamic Response of Saturated Soil - Foundation System Acted upon by Vibration’, Journal of Earthquake Engineering, Vol. 21, No. 7, pp. 1158–1188, Taylor & Francis Group, LLC, DOI: .

- [57]↑
Li, P., Lu, X. Chen, Y. (2004). Computer simulation on dynamic soil-structure interaction system. 13th World conference on Earthquake Engineering. Vancouver, Canada.

- [58]↑
Fattah, M.Y. Al-Mosawi, M.J. Al-Ameri, F.I. (2016). Vibration response of saturated sand - foundation system. Int J of Earthquakes and Structures; 11(1): 83–107.

- [59]↑
Fattah, M.Y. Zabar, B.Z. Mustafa, F.M. (2017). Effect of saturation on response of a single pile embedded in saturated sandy soil to vertical vibration. European Journal of Environmental and Civil Engineering; 24(3): 381–400. https://doi.org/10.1080/19648189.2017.1391126.

- [60]↑
Fattah, M.Y. Al-Neami, M.A. Jajjawi, N.H. (2014). Prediction of liquefaction potential and pore water pressure beneath machine foundations. Central European Journal of Engineering. DOI: .

- [62]↑
Mittal, R.K., Gajinkar, V. (2014). Evaluation of soil-structure interaction guidelines in indian seismic codes.

*International journal of civil engineering and technology, 5(5)*: 97–104. - [63]↑
Worku, A. (2014). Soil-structure- interaction provisions: A potential tool to consider for economical seismic design of buildings.

*Journal of the South African institution of civil engineering, 1(56):*54–62. - [64]↑
Jayalekshmi, B.R., Chinmayi, H.K. (2014). Soil–Structure Interaction effect on seismic force evaluation of rc framed buildings with various shapes of shear wall: as Per IS 1893 and IBC.

*Indian geotechnical journal. DOI*.

# Appendix I

To arrive at the corrected ** G** and damping

**value, these steps are followed:**

*ξ*- Identify the bedrock level (
*d*) of the site. - Find out the shear wave velocity from the expression:
*G*=*ρ ν*_{s}^{2}. - Find out the free field time period of the site for the n
^{th}mode from the expression: .${T}_{n}=\frac{4.d}{\left(2n-1\right).{\nu}_{s}}$ - Based on the site response spectra/spectra given in code and damping value as obtained in soil report obtain the acceleration
*Sa (RPA 99 – version 2003 [16]).* - Obtain shear strain for the soil profile based on the expression:
.${\gamma}_{z}=-\frac{16.\beta .{S}_{a}.d.\rho}{\left(2n-1\right).\pi .\left(\pi +2\right).G}sin\frac{\left(2n-1\right).\pi .z}{2.d}$ - Check if this strain is near or equal to the initial strain (10−3 to 10−4)%.
- If there exists a significant variation correct
*G*based on the equation*G*=*G*_{max}*(*1+*ψ/ψr).*where:- –
*G is the new value of shear modulus;* - –
*G*_{max}*is the 1st shear modulus calculate in eq. 2;* - –
*ψ is strain range; and* - –
*ψr is reference strain range*.

- –
- Find out the ratio
*G*/*G*_{max}. - Obtain new damping ratio based on Zhang's expression:
where:$\xi =0.333\frac{1+{e}^{\left(-0.0145P{I}^{1.3}\right)}}{2}\left[0.586.{\left(\frac{G}{{G}_{max}}\right)}^{2}-1.547\frac{G}{{G}_{max}}+1\right]$ - –
*PI is plastic Index.*

- –
- Repeat the steps as mentioned from
to*2*till the strain is same as the previous cycle.*7*- –When the value for which the strain becomes constant is the corrected dynamic shear modulus of the soil.

# Appendix II

# Appendix III

Ordinate of elastic and inelastic design spectra for RPA99-2003, EC8-2004, IBC-2015 and IS-1893-2002.

T ≤ T_{B} | T_{B} ≤ T ≤ TC | T ≥ T_{C} | Seismic base shear | |
---|---|---|---|---|

RPA99v2003 | ||||

EC8-2004 | S=2.5._{e}a._{g}S.η | |||

S = 2.5._{d}a._{g}S.η | ||||

IBC-2015 | S = _{a}S (_{DS}T_{0}≤T≤T)_{S} | |||

IS-1893-2002 |