Research Article  Open Access
Practical SoilShallow Foundation Model for Nonlinear Structural Analysis
Abstract
Soilshallow foundation interaction models that are incorporated into most structural analysis programs generally lack accuracy and efficiency or neglect some aspects of foundation behavior. For instance, soilshallow foundation systems have been observed to show both small and large loops under increasing amplitude load reversals. This paper presents a practical macroelement model for soilshallow foundation system and its stability under simultaneous horizontal and vertical loads. The model comprises three spring elements: nonlinear horizontal, nonlinear rotational, and linear vertical springs. The proposed macroelement model was verified using experimental test results from largescale model foundations subjected to small and large cyclic loading cases.
1. Introduction
Several researchers ([1–5], among others) have investigated extensively the subject of soilstructure interaction (SSI). Rayhani et al. [6] showed that soilstructure interaction might amplify or attenuate the base shear through inertial and kinematic interactions. Stewart et al. [7] concluded that the effects of SSI on rigid structures founded on soil are more significant, when compared to flexible structures. Bobet et al. [8] showed that the soilstructure system is a function of the relative rigidity of the structure compared to that of the ground. SSI is shown to be significant in the presence of soft soils or when structural mass is very large [9]. As a consequence, the SSI problem can be excluded from computations if the soil where the structure is founded is very rigid. When considering the effect of soilstructure interaction in base isolated multistoried structures on elastic layered soil, Spyrakos et al. [10] have found that SSI effects are significant for squat lightweight buildings on low stiffness soilstratum. Although their study dealt with harmonic excitations, it still gives an insight on the danger of neglecting SSI in the design of base isolated buildings.
In the seismic resistant design of structures, we are most interested in the strength reduction factors to account for the nonlinear behavior that might be experienced by a structure subjected to an earthquake ground motion. Few researchers [11, 12] have recently attempted to assess the effect of SSI on the strength reduction factors. Eser et al. [11] have shown that the presence of soft soils reduces the strength reduction factors, which is primarily controlled by the changes in the structural period and displacement ductility.
Incorporation of SSI requires explicit modelling of soilfoundation system adequately. For instance, several models have been proposed depending on the foundation type, its embedment, and its rigidity ([13–15], among others). The Federal Emergency Management Agency [16] required that the foundation stiffness should be determined with one of the following three methods:(i)uncoupled spring model comprising three spring elements, for shallow foundations that are stiffer than the supporting soil;(ii)a finite element formulation of linear (or nonlinear) foundation behavior using Winkler models, for shallow foundations that are less stiff than the supporting soil;(iii)decoupled Winkler model, for shallow foundations that are flexible with respect to the supporting soil.However, El Ganainy and El Naggar [17] have demonstrated that the decoupled technique (Beam on a Nonlinear Winkler Foundation, BNWF) is not capable of predicting accurately the settlement seen in foundations on soft soils. Although it can be used to predict the overall deformation behavior of foundations, the BNWF requires a large number of nonlinear springs, which is considered as a major drawback [17].
To address some of the abovementioned issues, macroelement formulations have been proposed. The first formulation has been developed by Nova and Montrasio [18] and later modified and/or extended by other researchers [19–22] and recently the formulations developed by Gajan et al. [23] and Shirato et al. [3]. The major advantages of such models are their simplicity and their ability to capture the global response of bearing foundations [22, 24, 25]. However, on one hand, calibrating macroelement parameters pose a constraint on their adoption for practical applications. And on the other hand, knowing that most of the available macroelement models are based on specified bounding surfaces poses another problem for their capability to cover a wide range of problems [17].
In a completely different modelling approach, El Shamy and Zamani [26] proposed a new 3D particlebased technique using the discrete element method (DEM) to analyze the seismic performance of soilfoundationstructure systems. In their model, the soil is idealized as a collection of spherical particles using DEM; the footing is considered as a rigid block, whereas the structure is modelled using a number of spherical particles in the form of a column, which can be clamped to simulate a rigid structure or bonded to simulate a flexible structure of predefined rigidity.
To overcome the difficulties in performing complete nonlinear simulations, Seylabi et al. [27] proposed an equivalent linearization of nonlinear soilstructure systems considering both the effect of SSI and the nonlinear behavior of the structure on equivalent linear parameters. In their model the structure is modelled as an elastoplastic singledegreeoffreedom system (SDOF) and the soil beneath the structure is modelled by a discrete model combining different spring and dashpot elements.
In this paper, a new macroelement model is developed for the analysis of the nonlinear response of shallow foundations under cyclic loading. This model may easily be incorporated into available structural analysis programs such as OpenSees [28]. The soilfoundation system is simulated using three spring elements: horizontal and rotational nonlinear springs and a linear vertical spring. The nonlinear springs are assigned appropriate nonlinear model of plasticity with material degrading parameters.
2. Proposed Macroelement Model
The problem being studied here is that of a shallow foundation of any shape embedded in soil and subjected to simultaneous axial and lateral forces, as shown in Figure 1. The foundation is considered to be very stiff. The depth of embedment is . The proposed model incorporates three types of springs:(i)vertical translational elastic spring with stiffness ;(ii)shear inelastic spring with preyield stiffness ;(iii)rotational inelastic spring with preyield stiffness .These equivalent springs represent the foundationsoil system. The macroelement model replaces the system soilshallow foundation, thus decreasing considerably both the overall number of degrees of freedom and the computation effort required to run large models.
2.1. Constitutive Equations
Two material models are considered in this study, namely, the BoucWen model [29, 30] for the shear and rotational springs and the linear model for the vertical spring.
2.1.1. Shear and Rotational Springs
Model Assumptions. The smooth BoucWen model of hysteresis by Bouc [29] and Wen [30] has found many engineering applications. For instance, the use of original BoucWen model and its extensions to soilstructure interaction include [31–33].
In this study, we considered the Baber and Noori [34] extension to the original BoucWen model. This version includes the degrading behavior observed in many engineering materials.
The nonlinear behavior of the soilshallow foundation system is modelled via the nonlinear shear and nonlinear rotational springs. A force is mobilized at the shear spring and a moment is mobilized at the rotational spring. In what follows, we use a general force to denote and , depending on which spring we are modelling.
The constitutive relationship for is expressed in the BoucWen model as a linear part and a hysteretic part:where is the deformation; is the preyield stiffness; is the mobilized force at the beginning of yielding; is the posttopreyield stiffness ratio; and is a hysteretic quantity controlling the nonlinear behavior. The latter is governed by the following differential equation with respect to time :Equation (2) represents the Baber and Noori [34] formulation of the rate of hysteretic deformation, which accommodates degradation. In (2) , , and are parameters that control the shape and transition from elastic to inelastic regions of the hysteretic loop, while , , and are variables that control the stiffness degradation and material deterioration. Note that the yield deformation of the spring does not appear in (2), unlike the original BoucWen model. In fact, the adoption of the BaberNoori [34] formulation in calculating forces is easy by considering that and are expressed asin which and can take values as in the original BoucWen model of hysteresis. In this case, is a dimensional quantity and is bounded between and . The evolution of material degradation is governed by the following equations [34]:where is defined by the following rate equation:And , , and are the parameters that control the degradation.
Incremental Response Equations. The tangent of (i.e., ) is computed using (1):Note that this is the continuum tangent and not the algorithmically consistent tangent. It is clear from (6) that the derivative of with respect to is needed. To get it, (2) can be rewritten asHence,From this point, we will derive incremental response equations to obtain computer implementable equations. In what follows, the same procedure developed by Haukaas and Der Kiureghian [35] is used. The force at time is obtained asThe variable is next discretized by a BackwardEuler solution scheme:It is seen that cancels from the equation, yielding a nonlinear equation in . A Newton scheme of the form to solve the general nonlinear equation of the form is employed to solve for in (10) [35].
The equations describing the degradation behavior are described as follows:where is found by discretization of the rate equation in (5) using the BackwardEuler solution scheme:where again cancels. Note that , , and are history variables that must be stored at each converged step.
For the complete implementable procedure to obtain and the algorithmically consistent tangent (referred to here as the Local Newton Routine), the reader is referred to [35].
2.1.2. Vertical Spring
The vertical force mobilized in the vertical translational spring is given byin which is the elastic vertical stiffness of the soilfoundation system and is the vertical displacement of the spring. Note that an appropriate constitutive model could be assigned to the vertical spring; however, in this paper, the spring is considered linear elastic, but the effect of axial loads is present on the soilshallow foundation global behavior.
2.2. Equilibrium Equations and Model Implementation
The macroelement model described in the previous subsections is used to simulate the soilshallow foundation interaction.
Initially, we assume that the shear and rotational springs are linear; then we replace them by the general forces and , respectively. The total potential energy of the model with reference to Figure 2 is expressed as follows:where and are the foundation’s top horizontal and vertical displacements, respectively, and are dependent on the depth of embedment and the spring’s displacements:In (15b), , positive in case of compression, is the sum of the vertical component of the displacement originating from the axial flexibility of the soil and the additional vertical displacement that happens in the laterally deformed configuration shown in Figure 2.
The above kinematic equations (see (15a), (15b)), which relate the total displacement and the vertical displacement to the internal displacements , , and , assume large displacements and large rotations.
Considering the model subjected to the axial load with the resulting horizontal displacement , the vector of unknowns is then . Using Castigliano’s second theory by imposing the stationary of with respect to , , and yields the following governing equations:The following equation should also be considered to permit the displacementcontrolled cyclic analysis by imposing the top lateral displacement :However, when the lateral forcecontrolled cyclic response is desired, is to be imposed. Hence, the system of governing equations to be solved is reduced to the three equilibrium equations (see (16a), (16b), and (16c)). The lateral and vertical displacements may be calculated at the end of the analysis using the kinematic equations (see (15a), (15b)).
Considering the case of lateral displacementcontrolled analysis, (16a), (16b), (16c), and (16d) are rewritten in the following form:Note that the resultant moment found by the equilibrium of the macroelement in the deformed configuration needs also to be applied at the top of the foundation. Also note that in (17), the terms and were replaced by the general forces and , respectively. These forces are assigned the constitutive model developed in Section 2.1.1.
The system can be solved using Newton’s method following the pseudocode provided below (referred to here as the Global Newton Routine): Previous converged solution (1)While ()(2)Update , , , and using the Local Newton Routine(3)Compute (4)Compute the Jacobian matrix (5)(6) end
where the Jacobian matrix is defined as . The tangents and required in the Jacobian matrix should be algorithmically consistent tangents (Section 2.1.1).
From the current estimation of the top lateral displacement and axial load , the code returns the vector that satisfies equilibrium and kinematics (see (17)). In each global Newton iteration the current values of , , , and based on the current values of and are updated using the Local Newton Routine (Section 2.1.1).
The above pseudocode (Global Newton Routine) along with the Local Newton Routine (that implements the BoucWenBaberNoori model of hysteresis) has been implemented numerically in MATLAB (Mathworks, Inc.).
3. Verification with Experimental Results
To verify the validity of the proposed macroelement model in predicting the cyclic behavior of the soilshallow foundation system, its predictions are compared with experimental results. In the framework of the TRISEE Project (3D Site Effects and SoilFoundation Interaction in Earthquake and Vibration Risk Evaluation) a program of largescale 1 g model has been tested to investigate the response of soilshallow foundation under cyclic and dynamic loads. The experiment was carried out at ELSA (European Laboratory for Structural Assessment) in Ispra, Italy; test results are reported in many references ([36], among others).
The TRISEE experiment consists of three phases; only Phases I and III are considered in the comparison because the proposed macroelement model is restricted to 2D loading conditions only. However, the model may be extended to include the 3D case.
The experimental setup consists of a square steel shallow foundation mounted inside a stiff concrete caisson filled with Ticino sand. The embedment of the foundation was about 1 m, with a steel framework placed around the foundation to retain the sand. Two different soil relative densities have been used, and , representing high (HD) and low density (LD) soil conditions, respectively.
The HD and LD specimens were loaded vertically by and , respectively. The static safety factor was about 5 in both tests. In the HD test, the load was applied at 0.9 m above the foundation and in the LD test at . In Phase I, a series of unidirectional forcecontrolled small amplitude cycles was applied to identify the significance of nonlinear soil behavior. In Phase III, displacementcontrolled, unidirectional, increasing amplitude cycles were imposed to the top of the foundation. For further information on experimental setup and results refer to [36].
The new macroelement model is used to simulate the foundation. The parameters of the numerical model are presented in Table 1. These parameters have been calibrated using experimental momentrotation and horizontal forcehorizontal displacement curves. Parameters for the BoucWenBaberNoori model of hysteresis are given in Tables 2 and 3, for HD and LD tests, respectively. The parameters were again calibrated with experimental results.



Figures 3–6 compare the experimental and the numerical hysteretic horizontal forcehorizontal displacement curves and momentrotation curves for Phases I and III. The numerical model reproduces correctly the overall behavior of the foundation, verifying the ability of the proposed macroelement model to simulate the cyclic behavior of shallow foundations.
The uplift is important for the HD Phase III test, and this can be observed from the S shaped momentrotation curve (see Figure 5). For the LD sand, only plasticity is developed and the uplift is not present.
4. Foundation Stiffness Matrix
In the previous section, large rotations have been considered in formulating the macroelement model. However, assuming that small rotations permits the construction of the foundation stiffness matrix. In this case, (17) can be rewritten asRyan et al. [37, 38] studied the stability of bearing isolators and constructed the leadrubber bearing stiffness matrix relating a change in the isolator forces to the change in displacements. Following the same idea, the foundation matrix , relating a change in foundation forces to the change in displacements , is derived in three steps:(1)Differentiate the equations of equilibrium (in (18): , and ) resulting in in which and the matrices and are given by(2)Differentiate the equations of kinematics (see (15a) and (15b)):(3)Substitute from (19) into (21): where . The resultant flexibility matrix, , relates the displacement increment to the force increment and is given byThen foundation stiffness matrix is the inverse of the flexibility matrix, . The stiffness matrix is computed after calculation of displacements and computation of forces to satisfy the governing equations (see (17)). In a structural analysis program, the routine to compute the foundation stiffness matrix may be incorporated into existing structural analysis software and used for analysis of structures on shallow foundations.
5. Conclusions
A practical macroelement model is presented in this paper to simulate the response of soilshallow foundation systems. The proposed model was verified against experimental test results of largescale model foundations subjected to small and large loading cycles. A summary of the main points presented in this paper is given below:(1)The proposed macroelement model can simulate with a good accuracy the lateral response and rocking of shallow foundations under quasistatic cyclic loadings.(2)The uplift could be simulated adequately using wellchosen parameters of the BoucWenBaberNoori model of hysteresis.(3)The soil squeezeout phenomenon observed by El Ganainy and El Naggar [17] can be included by proper choice of degradation parameters of the BaberNoori model.(4)The proposed model does not take into consideration full coupling between the different springs. Nevertheless, the model represents a first step for future improvements.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This work was supported by the Sustainable Construction Material and Structural Systems (SCMASS) Research Group of the University of Sharjah, UAE.
References
 A. S. Veletsos and A. M. Prasad, “Seismic interaction of structures and soils: stochastic approach,” Journal of Structural Engineering, vol. 115, no. 4, pp. 935–956, 1989. View at: Publisher Site  Google Scholar
 T. Balendra and A. C. Heidebrecht, “Influence of different sites on seismic base shear of buildings,” Earthquake Engineering & Structural Dynamics, vol. 14, no. 4, pp. 623–642, 1986. View at: Publisher Site  Google Scholar
 M. Shirato, T. Kouno, S. Nakatani, and R. Paolucci, “Most impacting conventional design parameters on the seismic permanent displacement of shallow foundations,” Soils and Foundations, vol. 52, no. 2, pp. 346–355, 2012. View at: Publisher Site  Google Scholar
 M. Naseri and E. S. Hosseininia, “Elastic settlement of ring foundations,” Soils and Foundations, vol. 55, no. 2, pp. 284–295, 2015. View at: Publisher Site  Google Scholar
 N. Zamani and U. El Shamy, “A microscale approach for the seismic response of mdof structures including soilfoundationstructure interaction,” Journal of Earthquake Engineering, vol. 18, no. 5, pp. 785–815, 2014. View at: Publisher Site  Google Scholar
 M. H. T. Rayhani, M. H. El Naggar, and S. H. Tabatabaei, “Nonlinear analysis of local site effects on seismic ground response in the bam earthquake,” Geotechnical and Geological Engineering, vol. 26, no. 1, pp. 91–100, 2008. View at: Publisher Site  Google Scholar
 J. P. Stewart, S. Kim, J. Bielak, R. Dobry, and M. S. Power, “Revisions to soilstructure interaction procedures in NEHRP design provisions,” Earthquake Spectra, vol. 19, no. 3, pp. 677–696, 2003. View at: Publisher Site  Google Scholar
 A. Bobet, R. Salgado, and D. Loukidis, “Seismic design of deep foundations,” Publication FHWA/IN/JTRP2000/22, Joint Transportation Research Program, Indiana Department of Transportation and Purdue University, West Lafayette, Ind, USA, 2001. View at: Google Scholar
 Q. Gao, J. H. Lin, W. X. Zhong, W. P. Howson, and F. W. Williams, “Isotropic layered soilstructure interaction caused by stationary random excitations,” International Journal of Solids and Structures, vol. 46, no. 34, pp. 455–463, 2009. View at: Publisher Site  Google Scholar
 C. C. Spyrakos, Ch. A. Maniatakis, and I. A. Koutromanos, “Soilstructure interaction effects on baseisolated buildings founded on soil stratum,” Engineering Structures, vol. 31, no. 3, pp. 729–737, 2009. View at: Publisher Site  Google Scholar
 M. Eser, C. Aydemir, and I. Ekiz, “Soil structure interaction effects on strength reduction factors,” Structural Engineering and Mechanics, vol. 41, no. 3, pp. 365–378, 2012. View at: Publisher Site  Google Scholar
 B. Ganjavi and H. Hao, “Strength reduction factor for MDOF soilstructure systems,” Structural Design of Tall and Special Buildings, vol. 23, no. 3, pp. 161–180, 2014. View at: Publisher Site  Google Scholar
 M. Ornek, M. Laman, A. Demir, and A. Yildiz, “Prediction of bearing capacity of circular footings on soft clay stabilized with granular soil,” Soils and Foundations, vol. 52, no. 1, pp. 69–80, 2012. View at: Publisher Site  Google Scholar
 Z. Y. Ai, Z. X. Li, and Y. C. Cheng, “BEM analysis of elastic foundation beams on multilayered isotropic soils,” Soils and Foundations, vol. 54, no. 4, pp. 667–674, 2014. View at: Publisher Site  Google Scholar
 S. C. Dutta, K. Bhattacharya, and R. Roy, “Effect of flexibility of foundations on its seismic stress distribution,” Journal of Earthquake Engineering, vol. 13, no. 1, pp. 22–49, 2009. View at: Publisher Site  Google Scholar
 American Society of Civil Engineers, FEMA 356Prestandard and Commentary for the Seismic Rehabilitation of Buildings, Federal Emergency Management Agency, Washington, DC, USA, 2000.
 H. El Ganainy and M. H. El Naggar, “Efficient 3D nonlinear Winkler model for shallow foundations,” Soil Dynamics and Earthquake Engineering, vol. 29, no. 8, pp. 1236–1248, 2009. View at: Publisher Site  Google Scholar
 R. Nova and L. Montrasio, “Settlements of shallow foundations on sand,” Geotechnique, vol. 41, no. 2, pp. 243–256, 1991. View at: Publisher Site  Google Scholar
 R. Paolucci, “Simplified evaluation of earthquakeinduced permanent displacements of shallow foundations,” Journal of Earthquake Engineering, vol. 1, no. 3, pp. 563–579, 1997. View at: Google Scholar
 S. Pedretti, Nonlinear seismic soilstructure interaction: analysis and modeling method [Ph.D. thesis], Politecnico di Milano, Milano, Italy, 1998.
 C. Cremer, A. Pecker, and L. Davenne, “Cyclic macroelement for soilstructure interaction: material and geometrical nonlinearities,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 25, no. 13, pp. 1257–1284, 2001. View at: Publisher Site  Google Scholar
 C. Cremer, A. Pecker, and L. Davenne, “Modelling of nonlinear dynamic behaviour of a shallow strip foundation with macroelement,” Journal of Earthquake Engineering, vol. 6, no. 2, pp. 175–211, 2002. View at: Publisher Site  Google Scholar
 S. Gajan, J. M. Thomas, and B. L. Kutter, “Physical and analytical modeling of cyclic load deformation behavior of shallow foundations,” in Proceedings of the of the 57th Annual Meeting of Earthquake Engineering Research Institute, Ixtapa, Mexico, 2005. View at: Google Scholar
 C. T. Chatzigogos, A. Pecker, and J. Salençon, “Macroelement modeling of shallow foundations,” Soil Dynamics and Earthquake Engineering, vol. 29, no. 5, pp. 765–781, 2009. View at: Publisher Site  Google Scholar
 C. T. Chatzigogos, R. Figini, A. Pecker, and J. Salençon, “A macroelement formulation for shallow foundations on cohesive and frictional soils,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 35, no. 8, pp. 902–931, 2011. View at: Publisher Site  Google Scholar
 U. El Shamy and N. Zamani, “Discrete element method simulations of the seismic response of shallow foundations including soilfoundationstructure interaction,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 36, no. 10, pp. 1303–1329, 2012. View at: Publisher Site  Google Scholar
 E. Esmaeilzadeh Seylabi, H. Jahankhah, and M. Ali Ghannad, “Equivalent linearization of nonlinear soilstructure systems,” Earthquake Engineering and Structural Dynamics, vol. 41, no. 13, pp. 1775–1792, 2012. View at: Publisher Site  Google Scholar
 F. McKenna, G. L. Fenves, M. H. Scott, and B. Jeremic, Open System for Earthquake Engineering Simulation (OpenSees), Pacific Earthquake Engineering Research Center, University of California, Berkeley, Calif, USA, 2000.
 R. Bouc, “Modèle mathématique d'hystérésis: application aux systèmes à un degré de liberté,” Acustica, vol. 21, pp. 16–25, 1971 (French). View at: Google Scholar
 Y. K. Wen, “Method for random vibration of hysteretic systems,” Journal of Engineering MechanicsASCE, vol. 102, no. 2, pp. 249–263, 1976. View at: Google Scholar
 A. M. Trochanis, J. Bielak, and P. Christiano, “Simplified model for analysis of one or two piles,” Journal of Geotechnical EngineeringASCE, vol. 117, no. 3, pp. 448–466, 1991. View at: Publisher Site  Google Scholar
 D. Badoni and N. Makris, “Nonlinear response of single piles under lateral inertial and seismic loads,” Soil Dynamics and Earthquake Engineering, vol. 15, no. 1, pp. 29–43, 1996. View at: Publisher Site  Google Scholar
 N. Gerolymos and G. Gazetas, “Development of Winkler model for static and dynamic response of caisson foundations with soil and interface nonlinearities,” Soil Dynamics and Earthquake Engineering, vol. 26, no. 5, pp. 363–376, 2006. View at: Publisher Site  Google Scholar
 T. T. Baber and M. N. Noori, “Random vibration of degrading, pinching systems,” Journal of Engineering Mechanics, vol. 111, no. 8, pp. 1010–1026, 1985. View at: Publisher Site  Google Scholar
 T. Haukaas and A. Der Kiureghian, “Finite element reliability and sensitivity methods for performancebased earthquake engineering,” PEER 2003/14, Pacific Earthquake Engineering Research Center, College of Engineering, University of California, Berkeley, Calif, USA, 2004. View at: Google Scholar
 P. Negro, R. Paolucci, S. Pedretti, and E. Faccioli, “Large scale soilstructure interaction experiments on sand under cyclic load,” in Proceedings of the 12th World Conference on Earthquake Engineering, Auckland, New Zealand, 2000. View at: Google Scholar
 K. L. Ryan, J. M. Kelly, and A. K. Chopra, “Experimental observation of axial load effects in isolation bearings,” in Proceedings of the of the 13th World Conference on Earthquake Engineering, Vancouver, Canada, August 2004. View at: Google Scholar
 K. L. Ryan, J. M. Kelly, and A. K. Chopra, “Nonlinear model for leadrubber bearings including axialload effects,” Journal of Engineering Mechanics, vol. 131, no. 12, pp. 1270–1278, 2005. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Moussa Leblouba et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.