Modeling and Analysis of Soil-Pile Interaction
Modeling and Analysis of Soil-Pile Interaction
Modeling and Analysis of Soil-Pile Interaction
1 Introduction
The behavior of any structural system depends on the behavior of its key components.
Generally, for civil engineering structures, the complete structural system comprises of
3D frames in superstructure, its foundation and the surrounding soil. In-depth under-
standing of the interaction among the aforementioned components is of utmost
importance for an efficient, safe and economic design of overall system. Due to rapid
increase in the urbanization, the construction trend has been shifted towards high-rise
buildings, widespread road networks, and long span bridges etc. In most of these
structures, generally, piles are selected as an appropriate foundation type to resist, the
static as well as the anticipated dynamic loading. Therefore, for the accurate predictions
of design demands, a thorough understanding of soil-pile interaction is of vital
importance. The soil-pile interaction comprises of two phenomena i.e., kinematic
interaction and inertial interaction. The effect of ground motions on the foundation is
referred to as the kinematic interactions, whereas the effect of superstructure loads on
foundation is termed as inertial interactions (Fig. 1). Under dynamic loading, both
interaction phenomena are needed to be considered in the analysis for an appropriate
response prediction of soil-pile interaction (Peiris 2014).
Kinematic Interaction
Inertial Interaction
For the dynamic analysis and design of structures, a common practice is to model
the structure in such a way that its movement at the base is restricted. However, in
actual scenario, there is always some flexibility at the supports which decreases the
stiffness of the system and enhances the natural period. Depending upon the type of
soil/site and the characteristics of anticipated ground motions, the design response
spectrum may show a significant variation in the spectral acceleration with slight
change in natural period of the overall system. In this way, an increase in natural period
can significantly alter the overall dynamic response of soil-pile interaction. Also, apart
from affecting the shift in natural period, different soil media can exhibit different
130 S.A. Suhail et al.
The most important part in modeling of soil-pile interaction is the modeling of soil
media. Two traditional modeling approaches are available in the literature for this
purpose, (a) modeling using the Winkler model (1867) and, (b) the continuum-based
approach. The main objective of foundation is to uniformly distribute the loads
transferred from the superstructure to the soil. Considering this viewpoint, both
approaches have their own advantages and limitations. The basic principle and the
merits and de-merits of both approaches are discussed below.
The basic Winkler (1867) model idealizes the soil as identical and linearly elastic
springs. Each spring behaves independently and deforms only under the loaded area
(Fig. 2). The elastic load-deformation relationship is given by the following form of
Hook’s law.
p ¼ kw ð1Þ
external load becomes confined. Furthermore, springs used to model soil are linear
elastic whereas under dynamic loading, soil behaves nonlinearly which is another
serious limitation of Winkler model. To address these limitations, many researchers
proposed refined/improved Winkler models which are illustrated in the later sections of
this article.
Load
Foundation
Rigid Layer
Springs
In order to address the shortcomings of the classical modeling, various modified and
improved foundation models were proposed for both the Winkler approach and the
continuum-based approach. In general, two improvement techniques are employed by
most of the improved models. In the first technique, structural elements are introduced
to develop continuity in soil springs, whereas in the second approach, continuum-based
models are simplified to develop more realistic soil representation.
132 S.A. Suhail et al.
T T
Rigid Layer
Springs
p ¼ kw Tr2 w ð2Þ
d2w
p ¼ kw T ð3Þ
dx2
p ¼ kw Dr4 w ð4Þ
Here D is the flexural rigidity of the plate and can be evaluated using Eq. (5).
2
D ¼ Ep h3p = 12 1 lp ð5Þ
Modeling and Analysis of Soil-Pile Interaction 133
Rigid Layer
Springs
Where; Ep is the Young’s modulus of the plate material, hp is the plate thickness
and lp is the Poisson’s ratio of the plate material.
(c) Pasternak Foundation Model: Pasternak’s (1954) foundation model can be con-
sidered as modified version of Hetenyi (1950) foundation model as it also
incorporates the shear interaction mechanism among the spring elements by
introducing beam/plate element that undergoes shear deformations only (Fig. 5).
Vertical equilibrium is considered and pressure-deflection (p w) relationship is
defined by Eq. (6) (Pasternak 1954).
p ¼ kw Gr4 w ð6Þ
Rigid Layer
Springs
p ¼ kw ð7Þ
dw
mn ¼ k1 ð8Þ
dn
k2 G
1þ p ¼ r2 p þ k2 w Gr2 w ð9Þ
k1 k1
Foundation
First Layer k1
Shear Layer with
Shear Modulus G
k2
Second Layer
Rigid Layer
(f) Beam Column Analogy Model: A new subgrade model is introduced in the clas-
sical “beam on elastic foundation” model (Fig. 7) and the prevailing differential
equations are articulated as shown in Eq. 10. Modeling with this approach shows
significant improvement over Winkler hypothesis, however, it may also lead to
inaccurate results due to introduction of a fabricated shear force in the model.
Applied Force
Springs
Rigid Layer
d 4 wð xÞ d 2 w ð xÞ
Eb Ib 4
Cp2 þ Cp1 wð xÞ ¼ qð xÞ ð10Þ
dx dx2
Where Eb Ib is the flexural stiffness of the beam, wð xÞ is the settlement of beam,
qð xÞ is the load applied, Cp1 and Cp2 are constants (Horvath 1993).
(g) New Continuous Winkler Model: An idea was proposed to incorporate interac-
tions among springs by intermeshing of springs instead of introducing new
structural element in the model. The meshing is carried in such a way that the
interconnection among the springs is achieved (Fig. 8). This powerful idea leads
to the development of new continuous Winkler model. Physically, interaction
among springs without additional structural element is achieved by adding
another spring having specific axial stiffness which is not in contact with the
foundation. More details and the case studies for this model are available in the
literature (Kurian 2001).
(a) Vlasov Foundation: The variational principle had been adopted in continuum
approach to develop a foundation model that represents soil-structure interaction
in a more realistic way. This model limits the deformations of an elastic layer by
introducing a certain function in the formulation, as illustrated below:
(i) Vertical displacement is computed as function of hðzÞ which governs the
variation of deformation in vertical directions. Vertical displacement is
computed using Eq. (11).
d2w
p ¼ kw 2t ð12Þ
dx2
Where;
2
Eo H
Z dh
k¼ dz ð13Þ
ð1 vo Þ2 0 dz
Eo H
Z
t¼ h2 dz ð14Þ
4ð 1 v o Þ 0
E
Eo ¼ ð15Þ
ð1 v2 Þ
vo ¼ vð1 vÞ ð16Þ
(b) Reissner Model: This model is based on two assumptions, (a) the plane stresses in
the foundation are very small, hence can be ignored, and (b) the horizontal
stresses are zero (same as in Vlasov foundation model). The pressure-deformation
relationship at the interface between foundation and subgrade is determined as
disturbance of foundation below the slab and is shown in Eq. (17).
C2 2
C1 w C2 r2 w ¼ p r p ð17Þ
4C1
using the rigid stump test (Siva et al. 1969; Pranesh 1970; Siva et al. 1971; Siva et al.
1973; Vallabhan et al. 1991). The major shortcoming of Winkler approach is the lack
of interaction among the soil springs, as it models the soil media using independent
springs. In the case of continuum-based modeling approach, the young’s modulus of
soil is the necessary parameter. Since an accurate determination of the young’s mod-
ulus of soil is relatively complex, the Winkler modeling approach gained more
appreciation among the practicing geotechnical engineers (Bowles 1997).
The response of soil to dynamic loading is complex enough due to a number of reasons.
As discussed earlier, soil is composed of granular particles that are joined together by
intergranular forces. The stress-strain characteristics of soil govern the overall dynamic
soil-pile interactions. Under a general loading, the soil mass tends to attain explicit
structural configurations so that the potential energy becomes minimum. Under certain
level of loading, the soil stress-strain behavior remains elastic, after then it enters into the
plastic range, with an unproportioned increase in deformation with increasing loads.
This soil deformation is followed by viscoplastic deformation which is dominant in fine
grained soil. This initial deformation in the soil is mostly occurred due to the expulsion
of pore water under the loading while simultaneously transferring the pore pressure to
the solid soil mass. As a result of this behavior, the rate of strain decreases with time.
This whole process of expulsion of water from the soil under the action of loading and
squeezing of soil mass is termed as the primary consolidation. After attaining primary
consolidation, the soil still resists load for certain interval of time and secondary
compression of soil occurs following the logarithmic function of time (Bowles 1997). It
must be noted that the behavior of unsaturated soil having sufficient elastic range will be
different, as the strain will no longer be the function of time. On the other hand, for a
fully saturated soil sample, the strain will be a function of time under dynamic loading.
This hardening behavior of fine grained saturated soil must be incorporated in the model
for the accurate prediction of response of soil-pile structure system, since this
time-dependent property causes a reasonable increase in the strength characteristics of
the soil (cohesion and angle of internal friction). Furthermore, under dynamic loading,
the stress-strain behavior of soil no longer remains linear, as the soil exhibits inelastic
behavior. Hence, the nonlinearity of soil must also be considered in the modeling and
analysis of the dynamic soil-pile interaction problems.
Elastic Element
Plastic Element
elasto-plastic behavior of soil. This conjunction forms the schematic system referred to
as the St. Venant’s unit (Fig. 9). It is worth mentioning that St. Venant’s unit should be
used in large number and connected in parallel for smooth transition of soil from elastic
to plastic state. In case if only one element is used, the transition is abrupt and the
accuracy of results can be compromised (Dutta et al. 2002). The mathematical
expression for strain moduli for elastic and plastic range is shown in Eq. (18):
ru
eep ¼ Me r þ Mp log ð18Þ
ru r
Where; Mp and Me are the plastic and elastic strain moduli, respectively, and ru is
the ultimate load. This model appears to be useful, however, practical problems can
occur with regard to the selection of input parameters and arrangement of springs.
Numerous yield criterion such as Von Mises, Mohr-Coulomb, Drucker-Prager and
Tresca yield criterion (Harr 1966; Das 1983) had been available in the literature, which
can be adopted for modeling elasto-plastic idealization of soil. The post-yielding
behavior can be incorporated in the model using either deformation theory, flow theory
or incremental theory (Desai et al. 1987). In deformation theory, plastic strains are
uniquely defined by the state of stress, whereas in incremental theory, plastic strains are
dependent on multiple factors including the state of stress, increment in stress and strain
etc. Generally, the incremental theory is applied in modeling the post-yielding
behavior. The constitutive modeling of soil can be developed using one of the fol-
lowing methods; iterative method, initial strain method, initial stress method and
incremental method. The detailed formulations are available in the literature (Desai
et al. 1987; Ziendkiewicz et al. 1989). The key advantage of these formulations
includes relatively simple coding of yield function and flow rule, as only simple soil
parameters (which can be easily obtained) are required to implement the aforemen-
tioned formulations in computer programs. The 3D elasto-plastic model for the
unsaturated soil is also proposed recently (Sun et al. 2000).
Load P
Spring K
Dashpot
Coefficient θ
Various models attempting to idealize the rheological characteristics of clayey soil are
well documented in the literature. These include mechanical models, theory of plastic
flow and the molecular theory of flow (Vyalov 1978). The elastic element in con-
junction with viscous and plastic elements is a common way to model rheological soil
properties. Springs and dashpots are used in combinations for this purpose. These
springs and dashpots can be connected either in series i.e. Maxwell model (Fig. 10) or
in parallel i.e., Kelvin model (Fig. 11). Various conceptual details and explanation of
these models are available in the literature (Vyalov 1978).
Load P
Dashpot Spring K
Coefficient θ
The complete analysis of soil-pile interaction system involves the modeling and
analysis of piles and pile cap, the selection of boundary conditions, the modeling of soil
media, soil-pile interaction, non-linear soil response and various other important
components and phenomenon. Various analytical methods, based on idealization of
140 S.A. Suhail et al.
soil behavior, were developed over the past few decades to predict the response of
soil-pile system subjected to dynamic loading. The analytical methods can be cate-
gorized into two types i.e. the “beam on elastic foundation” approach and the
“continuum-based” approach. Their brief introduction and underlying concepts are
reviewed as follows.
Where; Ep and Ip are the elastic modulus and moment of inertia of pile foundation,
respectively, Q is the axial load, x is the vertical depth and y is the lateral deflection of
the pile at a depth x. It is worth mentioning here, that in the case of foundation springs,
the spring deflection multiplied by the spring constant results into resistive force of the
ground, unlike the conventional spring which provides the spring force. Using this
approach, various important characteristics such as the deflection of beam, bending
moment and shear force can be evaluated. The basic parameters required to calculate
the aforementioned characters are elastic modulus, beam dimensions, loading condi-
tions and the spring constant. The spring constant can be obtained conveniently using
empirical methods. Due to its simplicity, the beam on foundation approach has been
widely used for the dynamic analysis of soil-pile interaction. However, the prediction
of dynamic behavior of soil-pile interaction is relatively complex due to soil nonlin-
earity (Matlock et al. 1960). Also, several important factors are ignored in this
approach. These include soil-pile gapping, soil slippage, damping and analysis over the
wide range of frequencies, which are not possible to consider using elastic springs on
one side of the pile foundation.
Various refinements in the beam on foundation approach were proposed by the
introduction of dashpot in addition to springs (Fig. 13). Due to the addition of the
dashpot, the dynamic response of the soil-pile interaction can be obtained over a wide
range of frequencies (Nogami et al. 1989). Furthermore, to account for the nonlinearity
of soil media, soil is divided into two regions i.e., near field element and far field
element. Near field element accounts for the soil nonlinearity, whereas the far field
element governs the elastic behavior (Fig. 14) (Nogami et al. 1992). To account for the
soil-pile gapping and soil slippage phenomena, the near and far field element concept
was extended and used at both sides of pile foundation. The near field element com-
prises of nonlinear springs and stiffness is calculated assuming that the plane stress
conditions govern. The near field element is modeled by a homogenous isotropic
viscoelastic material. The second part is far field model, which accounts for the wave
propagation away from the pile.
Spring
k1 k2 k3
Mass
c1 c2 c3
Dashpot
Fig. 13. (a) Winkler model (b) Near and far field element & interface model
Interface Model
Pile
Expandable Frame
Pile
(a) (b)
General
Modelling
Details
Loading Boundary
Conditions Conditions
Finite
Element
Modelling
Soil
Soil-Pile
Behaviour
Interface
Modelling
Damping
finite element. The unknown displacements are evaluated using global loading con-
ditions and boundary conditions. For the purpose of dynamic analysis of a soil-pile
interacting system, the components of a typical finite element model are shown in
Fig. 15.
The general finite element modeling starts with the selection of appropriate finite
elements to represent the soil and pile. For 3D soil-pile interaction models, generally
8-noded brick elements are used to model the soil media. The boundary conditions are
used in dynamic analysis to reduce the analysis domain. Unlike static analysis, the
fixed boundary conditions cannot be applied at some distance from the interested
region, as it can reflect outward propagating waves back to the model and hence can
affect the accuracy of results. Therefore, the selection of appropriate boundary con-
ditions is quite important. Various options for the boundary conditions are available
such as the transmitting boundary conditions, absorbing boundary conditions, contin-
ued fraction absorbing conditions etc. The detailed discussion on the boundary con-
ditions is beyond the scope of this review. The adequacy of interface modeling between
the soil and pile can also significantly affect the dynamic response. It can be done in
two ways i.e., the perfectly bonded interface or frictional interface. The two important
parameters which must be considered while modeling the interface between soil and
pile are the soil slippage and soil gapping. The frictional modeling of interface is a
good option since it directly incorporates the effects of soil slippage and gapping in the
analysis. Various studies are available in the literature to develop interface elements
which can account for discontinuity at the interface (Goodman et al. 1968; Zienkiewicz
et al. 1970; Ghaboussi et al. 1973; Pande et al. 1979; Beer 1985).
Modeling and Analysis of Soil-Pile Interaction 145
With the propagation of waves in the soil media, their energy dissipates due to the
inherent damping offered by soil. For dynamic analysis using the FEA, two types of
damping are considered i.e., the material damping and geometrical/radiation damping.
The magnitude of damping should be appropriately selected based on the empirical
investigations and existing data. Similarly, the selection of appropriate constitutive soil
model is also a pre-requisite for the accurate dynamic analysis of soil-pile interaction
problems. Two types of material models (linear/nonlinear elastic and elasto-plastic) are
generally used in the FEA. It is imperative to apply the initial loading conditions, e.g.
the dead loads and other static loads of similar nature, prior to the dynamic loading on
the model. These initial loading conditions should represent the actual loading that
foundation would resist in the absence of any dynamic load, and thus, should provide a
right starting point for the lateral load analysis.
The earliest applications of the finite element method involved the determination of
behavior of elastic pile considering the 2D soil models. Although the response of piles
was evaluated considering soil nonlinearity and soil-pile slippage, the early studies
were limited to kinematic interaction (Trochanis et al. 1991). Wu et al. (1997a; 1997b)
performed the dynamic analysis of soil-pile interaction by considering both the kine-
matic and inertial interaction phenomenon. A quasi 3D finite element model was
proposed which evaluates the time-dependent stiffness and damping factors in addition
to the dynamic response of soil-pile system (Wu et al. 1997a, 1997b). The application
of the FEM for dynamic analysis was extended by introducing the plasticity-based
hierarchical single surface (HiSS) to model the soil media. Various important char-
acteristics such as the deformation modes of bonding, soil slippage, pile-soil separation
and rebounding of piles can be efficiently evaluated by adopting the combination of the
FEM and the HiSS (Cai et al. 2000). The FEM was also employed to predict the
dynamic response of socketed and floating piles as well as micro piles (Naggar et al.
2000; Sadek et al. 2003). The 3D finite element analysis was also used to investigate
the effects of plasticity of soil (Maheshwari et al. 2004). The major disadvantage of
performing dynamic analysis using the FEM is high computation cost. This short-
coming makes it less common in practical design applications. However, the
reduced-dimension finite element formulations are also proposed recently to carry out
the nonlinear analysis of vibrating discs in 2D unbounded domain, considering the soil
nonlinearity. More recently, a simplified finite element formulation is presented as a
fair compromise between the realistic representation of soil-pile interaction and the
accuracy in numerical results (Thammarak et al. 2011).
7 Concluding Remarks
The review suggests that the beam on foundation approach remained the most widely
used analysis approach for soil-pile interaction problems in last few decades. However,
one of the primary shortcomings associated with this approach is that it does not
considers the soil media as 3D continuum, rather it treats it as the one-dimensional
discretized springs. Hence, the radial and 3D interaction components are not consid-
ered. Many improvements had been carried out to the traditional beam on foundation
approach to evaluate the dynamic response of soil-pile interaction system, however,
146 S.A. Suhail et al.
such improvements are not commonly used in practical design due to their complexity.
The beam on foundation approach is also criticized for not accounting the interaction
among the discretized springs, since the displacements at one point does not affect any
other point in the analysis domain. On the other hand, continuum-based approach
addresses most of the short comings associated with beam on foundation approach.
This approach treats the soil as a continuum material and also considers the 3D
interaction among the soil-pile-structure. The most common continuum-based method
(i.e., the finite element method, FEM) is discussed in detail. The FEM have recently
gained more popularity due to its generality and diverse options for modeling the
complex conditions. The complex soil behaviors such as the material nonlinearity,
non-homogenous material conditions, geometrical changes etc., can be modeled with
reasonable accuracy. However, the use of FEM demands significant amount of com-
putational cost to generate results with high precision. Based on the review of modeling
and analysis techniques for dynamic analysis of soil-pile interaction, following are
some observations.
(a) To anticipate the dynamic response of soil-pile interaction, the combined effect of
kinematic and inertial interaction must be considered in the modeling and
analysis.
(b) Although Winkler hypothesis have obvious shortcomings and limitations, it can
still yield reasonably accurate results. For the practical point of view, Winkler
idealization can be used in combination of modeling structures with fixed base.
(c) In the case of clayey soil, the consolidation behavior is of utmost importance.
Under such conditions, the nonlinear stress-strain relationship must be considered
for the dynamic analysis. The clayey soil has low permeability and
time-dependent behavior under loading, which can be more accurately modeled
using the viscoelastic idealization.
(d) The complete soil-pile interaction system can be modeled efficiently using the
discretization techniques. The finite element modeling can provide reasonably
accurate response predictions and is also quite useful to incorporate the nonlin-
earity, non-homogeneity and anisotropy of soil.
(e) The beam on foundation approach have obvious limitations for the dynamic
analysis, but are still used widely for the design of soil-pile interaction systems.
The main reason of its widespread usage is its simplicity and less computational
cost. However, the accuracy of results can be lower compared to the detailed finite
element analysis, which can now be regarded as the state-of-the-art practice.
This article is an attempt to briefly summarize the widely-used modeling and
analysis methods for the dynamic analysis of soil-pile systems. Further details of any
particular modeling technique or analysis method can be found in the bibliography list
provided as the references.
Modeling and Analysis of Soil-Pile Interaction 147
References
Allotey, N., El Naggar, M.H.: A numerical study into lateral cyclic nonlinear soil–pile response.
Can. Geotech. J. 1281, 1268–1281 (2008a). doi:10.1139/T08-050
Allotey, N., El Naggar, M.H.: Generalized dynamic Winkler model for nonlinear soil–structure
interaction analysis. Can. Geotech. J. 573, 560–573 (2008b). doi:10.1139/T07-106
Ashour, M., Norris, G.: Modelling lateral soil-pile response based on soil-pile interaction.
J. Geotech. Geo-environ. Eng. 126(5), 420–428 (2000)
Basu, D., Salgado, R., Prezzi, M.: Analysis of laterally loaded piles in multilayered soil deposits.
Joint Transportation Research Program, 151 (2008). doi:10.5703/1288284313454
Bowles, J.E.: (1997). Foundation Analysis and Design, 5th (edn.) Engineering Geology, vol.
20 (1997). doi:10.1016/0013-7952(84)90010-3
Cai, Y.X., Gould, P.L., Desai, C.S.: Nonlinear analysis of 3D seismic interaction of soil–pile–
structure systems and application. Eng. Struct. 22(2), 191–199 (2000). doi:10.1016/S0141-
0296(98)00108-4
Chen, W.-F., Duan, L.: Bridge Engineering Handbook, 2nd edn. CRC Press, Boca Raton (2014).
doi:10.1201/b15616
Cox, W.R., Reese, L.C., Grubbs, B.R.: Field testing of laterally loaded piles in sand. In: Offshore
Technology Conference, pp. 459–464 (1974). doi:10.4043/2079-MS
Dutta, S.C., Roy, R.: A critical review on idealization and modeling for interaction among
soil-foundation-structure system. Comput. Struct. 80(20–21), 1579–1594 (2002). doi:10.
1016/S0045-7949(02)00115-3
El Naggar, M.H., Novak, M.: Nonlinear analysis for dynamic lateral pile response. Soil Dyn.
Earthq. Eng. 15(4), 233–244 (1996). doi:10.1016/0267-7261(95)00049-6
Goerguelue, U.: Beam Theories: The difference between Euler-Bernoulli and Timoshenko.
Lecture Handouts (2009). Accessed
Heidari, M., El, H., Jahanandish, M., Ghahramani, A.: Generalized cyclic p – y curve modeling
for analysis of laterally loaded piles. Soil Dyn. Earthq. Eng. 63, 138–149 (2014). doi:10.1016/
j.soildyn.2014.04.001
Hetenyi, M.: A general solution for the bending of beams on an elastic foundation of arbitrary
continuity. J. Appl. Phys. 21(1), 55–58 (1950). doi:10.1063/1.1699420
Horvath, J.S.: Beam-column-analogy model for soil-structure interaction analysis. J. Geotech.
Eng. 119(2), 358–364 (1993)
Jonathan, P., Angeles, L., Raymond, B., Greggory, L., Uildings, B., Stewart, B.J.P., Fenves, G.L.:
Seismic soil-structure interaction in buildings II: empirical findings. J. Geotech. Geo-environ.
Eng. 125, 38–48 (1999)
Kerr, A.D.: A study of a new foundation model. Acta Mech. 1(1), 135–147 (1965)
Kerr, A.D.: Elastic and viscoelastic foundation models. J. Appl. Mech. 31, 491 (1964). doi:10.
1115/1.3629667
Maheshwari, B.K., Truman, K.Z., El Naggar, M.H., Gould, P.L.: Three-dimensional finite
element nonlinear dynamic analysis of pile groups for lateral transient and seismic excitations.
Can. Geotech. J. 41(1), 118–133 (2004). doi:10.1139/t03-073
Mcgann, C.R., Arduino, P., Mackenzie-helnwein, P.: Simplified procedure to account for a
weaker soil layer in lateral load analysis of single piles. J. Geotech. Geoenvironmental Eng.
138, 1129–1137 (2012). doi:10.1061/(ASCE)GT.1943-5606.0000684
Mechanik, M.U.N.D.: Zeitschrift for Angewandte, 1 (1925)
Mindlin, R.D.: Force at a point in the interior of a semi-infinite solid. J. Appl. Phys. 7(5), 195
(1936)
148 S.A. Suhail et al.
Naggar, M.H.E., Bentley, K.J.: Dynamic analysis for laterally loaded piles and dynamic p–y
curves. Can. Geotech. J. 37(6), 1166–1183 (2000)
Nogami, T., Konagai, K.: Time domain flexural response of dynamically loaded single pile.
J. Eng. Mech. 114(9), 1512–1525 (1989)
Nogami, T., Lam, Y.C.: Two-parameter layer model for analysis of slab on elastic foundation.
J. Eng. Mech. 113(9), 1279–1291 (1987)
Nogami, T., Otani, J., Konagai, K., Chen, H.: Nonlinear soil-pile interaction model for dynamic
lateral motion. J. Geotech. Eng. 118(1), 89–106 (1992)
Pande, G.N., Sharma, K.G.: On joint/interface elements and associated problems of numerical
ill-conditioning. Int. J. Numer. Anal. Methods Geomech. 3(3), 293–300 (1979)
Pasternak, P.L.: On a new method of analysis of an elastic foundation by means of two
foundation constants. Gas. Izd. Lit. po Strait I Arkh, Moscow, Russia (in Russian) (1954)
Peiris, T.P.: Soil-Pile Interaction of Pile embedded in Deep Layered Marine Sediment under
Seismic Excitation. Ph.D. Thesis, Queensland University of Technology (2014)
Prezzi, M., Salgado, R., Basu, D.: A continuum-based model for analysis of laterally loaded piles
in layered soils. Géotechnique 59(2), 127–140 (2009). doi:10.1680/geot.2007.00011
Sadek, M., Shahrour, I.: Influence of piles inclination on the seismic behavior of groups of
flexible piles. Transp. Res. Rec. 1–14 (2003)
Stewart, J.P., Fenves, G.L., Seed, R.B.: Seismic soil-structure interaction in buildings I:
analytical methods. J. Geotech. Geoenvironmental Eng. 125, 26–37 (1999)
Terzaghi, K.: Evalution of conefficients of subgrade reaction. Geotechnique 5(4), 297–326
(1955). doi:10.1680/geot.1955.5.4.297
Thammarak, P., Tassoulas, J.L.: Reduced-dimension nonlinear finite-element analysis of a
vibrating disk in an unbounded domain. Earthq. Eng. Struct. Dynam. 40(January), 1353–1366
(2011). doi:10.1002/eqe
Thomas, H.R., Bendani, K.: Primary/secondary compression solution algorithm. J. Comput. Civ.
Eng. 114(5), 578–593 (1989)
Trochanis, B.A.M., Member, A., Bielak, J., Christiano, P.: Three-dimensional nonlinear study of
piles. J. Geotech. Eng. 117(3), 429–447 (1991)
Vallabhan, C.V.G., Das, Y.C.: Modified vlasov model for beam on elastic foundations.
J. Geotech. Eng. 117(6), 956–966 (1991)
Winkler, E.: Die Lehre Von Elasticitaet Und Festigkeit, 1st edn. H. Dominicus, Prague (1867)
Wu, G., Finn, W.D.L.: Dynamic elastic analysis of pile foundations using finite element method
in the frequency domain. Can. Geotech. J. 34, 34–43 (1997a)
Wu, G., Finn, W.D.L.: Dynamic nonlinear analysis of pile foundations using finite element
method in the time domain. Can. Geotech. J. 34, 44–52 (1997b)
Zeevart, L.: Foundation Engineering for Difficult Sub-Soil Conditions, 2nd edn. Van Nostrand
Reinhold Company, New York (1972)
Zienkiewicz, O.C., Taylor, R.L.: The Finite Element Method. Butterworth-Heinemann, Oxford
(2005). ISBN 978-0-7506-6431-8
Zienkiewicz, O.C., Best, B., Dullage, C., Stagg, K.G.: Analysis of nonlinear problems in rock
mechanics with particular reference to jointed rock systems. Supplemental Notes, vol. 3(8–14),
pp. 501–509 Report/Paper Numbers: No. 364 (1970)