Draft on Laterally Loaded Piles M.

Budhu April, 2002

Laterally loaded piles

Structures founded on piles are often subjected to lateral loads and moments in addition to
vertical loads. Lateral loads may come from wind, traffic, seismic events, waves and earth
pressures. Moments may come from the eccentricity of the vertical force, fixity of the
superstructure to the pile or piles and the location of the lateral forces on the pile with reference
to the ground surface.

When a pile is subjected to lateral forces and moments, the pile tends to bend or deflect as
illustrated in Figure 1. The deflection of the pile causes strains in the soil mass. To satisfy
equilibrium, the soil must provide reactions along the length of the pile to balance the applied
loads and moments.

Figure 1 Pile – soil response to lateral loads and moments

Because soil is a non-linear material, the soil reaction is not linearly related to the pile deflection.
Consequently at every point along the length of the pile, a non-linear relationship between soil
resistance (p) and pile deflection (y) exists as illustrated in Figure 1.

In designing laterally loaded piles, we need to know the pile deflection, particularly the
pile head deflection, to satisfy serviceability requirement and the bending moments for sizing the
pile. The pile head deflection depends on soil type, pile installation, pile flexibility (or pile
stiffness), loading condition and on how the pile is attached to the superstructure and pile cap. A
Draft on Laterally Loaded Piles M. Budhu April, 2002

pile that is attached to the pile cap such that no rotation occurs is called a fixed head pile (Figure
2a). A pile that is attached to the pile cap such that rotation is unrestricted is called a free head
pile (Figure 2b).

The mechanism of failure depends on the length to diameter or width ratio, soil type and the
fixity of the pile head. Free head piles tend to fail by rotation. Lateral loads and moments
applied to a free head pile is initial resisted by the soil near ground level. For very small pile
deflections, the soil behaves elastically and as the deflection increases the soil yields and then
permanent soil displacement occurs. The soil resistance is shifted to the lower part of the pile as
yielding progressively occurs from the top to the bottom of the pile. Fixed head piles tend to fail
by translation. Piles in general are neither fixed head nor free head. They have undermined
fixity somewhere between free head and fixed head conditions. You can view fixed head and
free head as two limiting conditions in which piles in practice will respond somewhere within
these limits.

Passive soil
pressure Translation of a
short, fixed head
Passive soil pile
pressure

(a) (b)

Passive soil Passive soil

pressure pressure

Yield Yield
and and
failure failure

(c) (d)

(a) Short, free head (b) Short, fixed head (c) Long, free head (d) Long, fixed head
Figure 2 Possible failure modes in short and long piles.

Two types of piles are normally defined to distinguish failure mechanisms. One is called a
short pile that is characterized by length to diameter (or width) ratio of less than 10 (sometimes a
ratio of up to 15 is quoted in the literature). Short, free head piles tend to fail by rotation about a
Draft on Laterally Loaded Piles M. Budhu April, 2002

point near the pile tip (or base). Rotation occurs when the sum of the soil resistance (called
passive pressures) at the top part of the pile is lower than the sum of the soil resistance at the
lower part of the pile. Short, fixed head piles tend to fail by translation. The failure of short
piles is due to soil failure. Long piles tend to fail by structural failure near ground level rather
than by soil failure.

Figure 3 Group piles subjected to lateral loads

Laterally loaded pile, particularly group piles are particularly difficult to analyze mainly because
of the complexity of the soil–structure (pile) interaction. The displacement and rotation are in the
directions of the resultant lateral load and resultant moment (Figure 3). Outer piles in a group
are subjected uplift (pull) and compressive (push) forces while the piles in the center translate at
the level of the superstructure connection. The response of a pile group to lateral loads and
moments is influenced by

1. Geometry of the group

2. Pile-soil interaction
3. Stiffness or flexibility of the piles
4. Load conditions
5. Individual pile response
6. Pile group response resulting from individual pile responses.

Analyses of lateral loaded piles

Basic structural mechanics

Consider an element of thickness, dz, of the pile shown in Figure 1 at a depth z. The free body
diagram for this element is shown in Figure 4. We will assume that the pile is symmetrical about
the z-axis, the loads lie in the yz plane, deflection of the pile occurs only along the y-axis, i.e., no
out of plane deflection, and shear deflection is negligible.
Draft on Laterally Loaded Piles M. Budhu April, 2002

Figure 4 Free body diagram of a pile section under lateral loads and moments.

Taking moments about O, we get

dz
( M + dM ) − M + Pz dy − Vdz + p dz =0 Eq.1
2
Neglecting higher terms (e.g. dz2) and rearranging Eq.1 gives
dM dy
+ Pz −V =0 Eq. 2
dz dz
Let us differentiate equation with respect to z, then
d2M d 2 y dV
+ Pz 2 − =0 Eq. 3
dz 2 dz dz
From basic mechanics, we have the following identities
d2M d4 y
= E I
p p Eq. 4
dz 2 dz 4
dV
=p Eq. 5
dz
p = − ky Eq. 6
where Pz is the vertical load, p is the soil resistance, y is the soil compression (or pile deflection),
Ep is the elastic modulus of the pile, k is a soil (stiffness) parameter, and Ip is the second moment
of area of the pile. Substituting the above Eqs. 4-6 into Eq. 3 gives the governing equation for a
laterally loaded elastic pile as
d4 y d2 y
E p I p 4 + Pz 2 + ky = 0 Eq. 7
dz dz
If the pile above the ground level is subjected to a distributed load, for example from water, then
the governing equation becomes
d4 y d2 y
E p I p 4 + Pz 2 + ky − W = 0 Eq. 8
dz dz
Draft on Laterally Loaded Piles M. Budhu April, 2002

where W is the resultant of the distributed load. Eq. 8 can be solved using numerical methods
such as the finite difference method or finite element method. The finite difference equations are
shown in Appendix D.

The term, ky, in the governing equation changes depending on the assumption made regarding
the soil response. One can assume an elastic soil or an elasto-plastic soil or some other type of
load-deformation response. The solution of Eq. 8 therefore varies on the assumption of how the
soil will respond.

Three methods of obtaining solutions to the governing equation based on three different types of
soil responses are described below.
.
1. Elastic beam on elastic foundation method

In this method, the pile is assumed to be an elastic beam that is attached to discrete springs
representing the soil (Figure 5). The soil is assumed to be elastic but discontinuous. The soil
stiffness parameter, k, is taken as the lateral subgrade modulus of the soil. It is assumed to be a
constant value for fine-grained soils and linearly increasing with depth for coarse-grained soils.
The lateral subgrade modulus, Kh, is the ratio of the horizontal soil reaction per unit area to the
lateral soil displacement. The lateral subgrade modulus can be obtained from field tests using
instrumented test piles. The lateral soil reactions and the bending moments are inferred from
grain gauges attached to the test pile. The determination of the lateral subgrade modulus is often
difficult and expensive. Empirical relationships are often used to relate results from simpler and
routine tests such as SPT, cone tests and or laboratory shear tests to Kh. A relationship between
Kh and elastic properties was proposed by Valsangkar et. al. (1973) as shown in Table 1. The
main problem with the elastic beam on elastic foundation approach is that it neglects soil
continuity and soil shearing resistance. This method is satisfactory for very small strain levels
(<0.001%)

Discrete springs
representing adjacent
soil mass

Figure 4 Simulation of adjacent soil mass as a set of discrete springs

2. p-y analysis
Draft on Laterally Loaded Piles M. Budhu April, 2002

The p-y (soil resistance – pile deflection) method (Matlock and Ripperger, 1958; Reese, Cox
and Koop, 1975 and others) is conceptually similar to the beam on elastic foundation method
except that the soil stiffness parameter k is not constant i.e. there is a non-linear relationship
between the soil resistance and its displacement (equal to the pile deflection. The non-linear
relationships are obtained from results of well-instrumented test piles. General procedures have
been developed to construct p-y curves for many soil types lessening the need to conduct lateral
load tests, at least, in the preliminary stages of design. However, it is recommended that pile
load tests be conducted whenever such tests are feasible and economical. Computer programs
(e.g. Com 624) are available for routine use. Like the elastic beam on elastic foundation
approach, the p-y method neglects soil continuity and soil shearing resistance.

3. Continuum Analysis

In this method of analysis, soil is treated as a continuous media with assumptions made on its
stress-strain behavior or constitutive relationships. The simplest stress-strain behavior is elastic
described by Hooke’s law. The solution gives the load-deformation response of the pile for
small strain levels (Poulos,1971; Randolph, 1981). Soil yielding or pile yielding cannot be
obtained from this analysis. To account for soil yielding, the soil can be assumed to be an elasto-
plastic material (Davies and Budhu, 1986; Budhu and Davies,1986, 1987).

H Shear stress
e on side faces

Soil resistance
Soil active stress
on bearing face
at back face

Figure 6 Soil stresses on a pile segment resulting from lateral loads that are taken into account in
Davies and Budhu method

We will adopt the elasto-(rigid) plastic solution using boundary element method proposed by
Davies and Budhu (1986) and Budhu and Davies (1986, 1987) for the general analysis of lateral
loaded single piles. The soil is assumed to behave linearly elastic at small strain levels. Yielding
occurs when the normal and shear stresses on the pile reach values determined from a limiting
equilibrium analysis. Bearing failure on the front face, shear stresses along the sides and tension
Draft on Laterally Loaded Piles M. Budhu April, 2002

failure along the back face (Figure 6) are incorporated into the analysis. Failure on the bearing
face occurs when
σ bf = N c su
where σbf is the limiting compressive stress, Nc is the bearing capacity factor that is assumed to
vary from 2 at the surface to a constant value of 9 at a depth of three pile diameters and below.
The limiting shear stress at the sides of the pile is
τ ss = α u su
where the adhesion factor, αu, is taken as 1 for soft clays. On the back face or tension face, the
normal stress cannot exceed the in situ lateral effective stress, i.e.,
σ bf ≤ Kisσ zo′
where Kis is the lateral earth pressure coefficient that accounts for pile installation. For soft
clays, Kis, is usually 1.0 or greater.
The important findings from the continuum analysis are as follows.
1. The analysis is valid only for a pile with a length greater than its effective length. The
effective length is the length of pile that is effective in resisting the loads and moments. Piles
longer than their effective lengths will behave identically.
2. The maximum bending moment in free head piles would normally occur at depth less than 5
diameters below ground level.
3. Eccentric lateral loads increase the maximum bending moment (regardless of the type of
analysis).
4. Fixed head piles reduce the pile head deflection by about 50%.
5. Laterally loaded piles, in general, fail by yielding of the pile section well before the soil fails.
6. The lower the flexibility of the pile the greater the pile deflection and the greater the
maximum bending moment. Flexibility, K, is the ratio of the elastic modulus of the pile to
the elastic modulus of the soil.
The predictions from boundary element analyses are curve-fitted to give equations that can be
easily used in practice. A summary of the relevant equations for single, solid piles embedded in
homogenous stiff clays (undrained shear strength constant with depth), homogenous soft clays
(undrained shear strength varies linearly with depth) and sand are summarized in Table 1.

Table 1 Design equations for laterally loaded piles

su
Soft Es
Clay

m
c z

Free head pile Fixed head pile

Effective pile length 2

Ep 2

(Le) Le = 1.3 K 9 ; K = Le = 1.3 K 9

mD
Draft on Laterally Loaded Piles M. Budhu April, 2002

Collapse load 1.2 ⎡ Hc 1.2 ⎡ W L⎤

(1 + 0.25α u ) + 0.008 r ⎤⎥
Hc WL = ( 4.5 + α u ) + 0.025 r ⎥
= ⎢ cDL2 g + 0.88 ⎢⎣ D ⎦
cDL 2
g + 0.88 ⎣ D ⎦
2
2 ∞>g≥ − ; 0 ≤ αu ≤ 1
∞>g≥ − ; 0 ≤ αu ≤ 1 3
3
Yield load/Working Hy 2
Hy 2
= 0.5n 3 βe = 0.8n 3
load cD 3
cD 3

3 1
−
βe = exp(−2f n ) 4 4

Replace Hy by Hw to get working load

Pile deflection uE = IuH
H
+ Ium
M
uE = IuH
H
2
(elastic) at ground mD mD 3 mD 2
1 1
surface
IuH = 3.2K − 3 IuH = 1.4K − 3
5
−
IuM = 5.0 K 9

Pile deflection u = Iuy uE u = Iuy uE

including yield at h − 1.4k 0.32 h n − 32k 0.43
Iuy = 1 + n n
Iuy = 1 + n

ground surface 40k 0.53

n 105k 0.54
n

H K H K
hn = ; kn = hn = ; kn =
cD 3 1000 cD 3 1000
Pile rotation (None)
(elastic) at ground θE = IθH H
+ Iθm M 4
surface mD3 mD
5
IθH = 5.0 K − 9
7
Iθm = 13.6 K − 9
h n − 14 k 0.32
Iθy = 1 + n

54 k 0.53
n

Pile rotation θ = Iθy θE (None)

including yield at h n − 14 k 0.32
Iθy = 1 + n
ground surface 54 k 0.53
n

Maximum bending ME = IMH HD ME = -IMH HD

moment (elastic) 2 2
IMH = 0.3 K for e = 0 9 IMH = 0.4 K 9

For eccentric loading

ImH = aKb
a = 0.6f
b = 0.17f –0.3
8 > f> 0.5

Maximum bending Mm = IMy ME Mm = IMy ME

moment including IMy = 1 + hn − 8k 0.32
n
IMy = 1 + hn − 30kn0.53
yield 96k 0.48
n 312kn0.56
Draft on Laterally Loaded Piles M. Budhu April, 2002

Maximum bending 2

moment occurs at Lm = 0.53 K 9 D; L > Le

depth, Lm, from
ground surface

su
Stiff Es
Clay

Free head pile Fixed head pile

Effective pile length L 4

(Le) e =0.5K 11
E
D ; K= p
Es
Collapse load 4.1 ⎡
(1 + 0.25αu ) + 0.008 r ⎤⎥
Hc WL Hc 1.2 ⎡ W L⎤
= = ( 4.5 + α u ) + 0.025 r ⎥
cDL2 g + 0.88 ⎢⎣
s u DL 1.2g + 1 ⎢⎣ D ⎦ D ⎦
2
2 ∞>g≥ − ; 0 ≤ αu ≤ 1
∞>g≥ − ; 0 ≤ αu ≤ 1 3
3
Yield load/Working Hy 2
Hy 2
= 0.5n 3 βe = 0.8n 3
load cD3 cD3
3 1
−
βe = exp(−2f n ) 4 4

Replace Hy by Hw to get working load

Pile deflection uE = IuH
H
+ Ium
M
uE = IuH
H
(elastic) at ground mD 2 mD 3 mD 2
1 1
surface
IuH = 3.2K − 3 IuH = 1.4K − 3
5
−
IuM = 5.0 K 9

Pile deflection u = Iuy uE u = Iuy uE

including yield at h n − 1.4k 0.32
n h n − 32k 0.43
n
Iuy = 1 + ; Iuy = 1 + ;
ground surface 40k 0.53
n 105k 0.54
n

Pile rotation (None)

(elastic) at ground θE = IθH H
+ Iθm M 4
surface mD3 mD
Draft on Laterally Loaded Piles M. Budhu April, 2002

5
IθH = 5.0 K − 9
7
Iθm = 13.6 K − 9
h n − 14 k 0.32
Iθy = 1 + n

54 k 0.53
n

Pile rotation θ = Iθy θE (None)

including yield at h n − 14 k 0.32
Iθy = 1 + n
ground surface 54 k 0.53
n

Maximum bending ME = IMH HD ME = -IMH HD

moment (elastic) 2 2
IMH = 0.3 K for e = 0
9 IMH = 0.4 K 9
For eccentric loading
ImH = aKb
a = 0.6f
b = 0.17f –0.3
8 > f> 0.5

Maximum bending Mm = IMy ME Mm = IMy ME

moment including IMy = 1 + hn − 8k n0.32
IMy = 1 + hn − 30kn0.53
yield 96k n0.48 312kn0.56

Maximum bending 2

moment occurs at Lm = 0.53 K 9 D; L > Le

depth Lm, from
ground surface

Ep
K= , Ep = elastic modulus of pile, m = slope of line depicting the linear variation of the
mD
elastic modulus of the soil with depth, c = slope of line depicting the linear variation of the
undrained shear strength ( su ) with depth and D is the diameter of the pile. For non-cylindrical
E pR I pR
and annual piles, the equivalent pile stiffness is, Ep = where Ip is the second moment of
⎛ π D4 ⎞
⎜ ⎟
⎝ 64 ⎠
area and the superscript R refers to the real pile.
K γK o
kn = , hn = H 3 , where c slope of the undrained shear strength with depth, Wr = , γ = soil
1000 cD c
unit weight (use effective weight for submerged section), Ko is the lateral earth pressure
coefficient, which should account for the method of pile installation ( a value of 1 is reasonable
e e
for soft clays), g = and f = . The above equations are for solid, cylindrical piles of diameter,
L D
Draft on Laterally Loaded Piles M. Budhu April, 2002

σ
D, with n = where σ is either the yield stress or working stress. For other types of geometry,
cD
⎛I⎞
⎜ ⎟
σ ⎝y⎠
an equivalent n value can be calculated from n = where I is second moment of area, D is
cD πD3
32
the projected width and y is the distance from the neutral axis to the extreme fibers.

Procedures
1. Obtain soils information for the site, particularly the variation of undrained shear strength
with depth.
2. Select pile (geometry and material, e.g. concrete).
3. Obtain loads and moments
4. Determine fixity condition – free head or fixed head.
5. Determine effective length and check that embedded length is greater than the effective
length.
6. Determine K, f, g, Wr, n, kn, hn
7. Calculate collapse load, yield load and working load.
8. Calculate elastic influence factors, IuH etc.
9. Calculate elastic ground surface deflection and rotation. The latter for only free head
pile.
10. Calculate elastic maximum bending moment.
11. Calculate yield influence factors, Iuy etc
12. Calculate total ground surface deflection and rotation.
13. Calculate total maximum bending moment.
14. Scale ground surface deflection as follows.
Neglecting the bending of the free-standing portion of the pile, the pile head deflection = u +
eθ

ILLUSTRATIVE EXAMPLE

A 460-mm diameter timber pile of total length 9.75 m is embedded to a depth, L, of 9.14 m in a
soft, normally consolidated, clay deposit of medium plasticity. The working stress, σw, for
timber of this quality is 11 MN/m2, and its Young’s modulus of elasticity is 10.3 GN/m2. A
lateral load is to be applied at a height, e, of 690 mm above ground level. Required are: (1) The
working load; (2) the lateral deflection; and (3) the maximum bending moment. Assume: αu = 1;
Ko = 1; γ = 9.4 kN/m3; c = 2.7 kN/m3; and m = 2.7 MN/m3, i.e., m/c =1,000.

SOLUTION
Step 1: Calculate normalized values
γ ' K o 9.4 × 1 Ep
Wr = = = 3.5 , L/D = 20, Wr L/D = 70, K = = 8300 , kn = 8300/1000 = 8.3, g = 0.075, f
c 2.7 mD
3 1
−
= 1.5, n = σy/cD = 8900, βe = exp(−2f 4 n 4 ) = 0.76

Step 2: Calculate collapse load and working load

Draft on Laterally Loaded Piles M. Budhu April, 2002

1.2 ⎡
(1 + 0.25α u ) + 0.008 r ⎤⎥ ;
Hc WL
= ⎢ Hc = 236 kN
cDL 2
g + 0.88 ⎣ D ⎦
2
Hw
= 0.5n 3
βe ; H w = 42kN
cD3
Step 3: Calculate working moment
M = Hw e = 29 kNm.

Step 4: Calculate elastic influence factors

IUH = 0.16, IθH = 0.033, and IθM = 0.012

Step 5: Calculate elastic deflection and rotation

Ground level:
Hw
uE = IuH 2
+ Ium M 3 ; uE = 16 mm
mD mD
Hw
θE = IθH 3
+ Iθm M 4 ; θE = 0.0083 rad.
mD mD

Step 6: Calculate yield influence factors

Iuy = 2.1, Iθy = 1.8

Step 7: Calculate deflection and rotation at pile head

Ground level:
u = Iuy uE ; u = 33 mm,
θ = Iθy θE; θ = .015 rad.
Pile head:
Neglecting the bending of the free-standing portion of the pile, upile head = 33 + 690 x 0.015 = 43
mm

Step 8: Calculate bending moment on pile

ImH = aKb; a = 0.6f, b = 0.17f –0.3 ; ImH = 3.5
ME = IMH HwD, ME = 68 kN.m.

hn − 8k n0.32
IMy = 1 + ; IMY = 1.54
96k n0.48
Thus, MMY = 104 kN m.
Draft on Laterally Loaded Piles M. Budhu April, 2002

Appendix D
dM dy
+ Px − Vv = 0 2.2
dx dx

d 2M d 2 y dVv
+ Px − =0 2.3
dx 2 dx 2 dx

d 2M d4y
= EI 2.4
dx 2 dx 4

dVv
=p 2.5
dx

p = −Es y 2.6

d4y d2y
EI + Px + Es y = 0 2.7
dx 4 dx 2

Vn = Vv cos S − Px sin S 2.8

dy
Vn = Vv − Px 2.9
dx

d 2M d2y
+ Px 2 + ky − W = 0 2.13
dx 2 dx

⎛ d 2M ⎞ ⎡ y R + y m −1 (− 2 Rm − 2 Rm −1 ) + y m (4 Rm + Rm −1 + Rm +1 )⎤ 1
⎜⎜ ⎟⎟ = ⎢ m − 2 m −1 ⎥ 4 2.14
⎠ m ⎣+ Ym +1 (− 2 Rm − 2 Rm +1 ) + Ym + 2 Rm +1
2
⎝ dx ⎦h

⎛d2y⎞ Px( y m −1 − 2 y m + y m +1 )
Px⎜⎜ 2 ⎟⎟ = 2.15
⎝ dx ⎠ m h2

Rm = E m I m 2.16

( ) (
Ym − 2 Rm −1 + Ym −1 − 2 Rm −1 − 2 Rm + Pxh 2 + Ym Rm −1 + 4 Rm + Rm +1 − 2 Pxh 2 + K m h 4 ) 2.17
( )
+ y m +1 − 2 Rm − 2 Rm +1 + Pxh 2 + y m + 2 Rm +1 − Wm h 4 = 0

⎛d2y⎞
⎜⎜ 2 ⎟⎟ = 0 2.18
⎝ dx ⎠ 0

⎛d3y⎞ ⎛ dy ⎞
R0 ⎜⎜ 2 ⎟⎟ + Px ⎜ ⎟ = 0 2.19
⎝ dx ⎠ 0 ⎝ dx ⎠ 0
Draft on Laterally Loaded Piles M. Budhu April, 2002

R−1 = R0 = R1 2.20

y −1 − 2 y 0 + y1 = 0 2.21

⎛ Pxh 2 ⎞ ⎛ Pxh 2 ⎞
y −2 ⎜
= y −1 ⎜ 2 − ⎟ ⎜
⎟ − y1 ⎜ 2 − R ⎟⎟ + y 2 2.22
⎝ R0 ⎠ ⎝ 0 ⎠

y 0 = a0 y1 − b0 y 2 + d 0 2.23

2 Ro + 2 R1 − 2 Px h 2
a0 = 2.24
R0 + R1 + k 0 h 4 − 2 Pxh 2

Ro + R1
b0 = 2.25
R0 + R1 + k 0 h 4 − 2Pxh 2

W0 h 4
d0 = 2.26
R0 + R1 + k 0 h 4 − 2Pxh 2

y m = a m y m +1 − bm y m + 2 + d m 2.27

− 2bm −1 Rm −1 + a m − 2 bm −1 Rm −1 + 2 Rm − 2bm −1 Rm + 2 Rm +1 + Pxh 2 (1 − bm −1 )

am = 2.28
cm

Rm +1
bm = 2.29
cm

dm =
( )
Wm h 4 − d m −1 a m − 2 Rm −1 − 2 Rm −1 − 2 Rm + Pxh 2 − d m − 2 Rm −1
2.30
cm

c m = Rm −1 − 2a m −1 Rm −1 − bm − 2 Rm −1 + a m − 2 a m −1 Rm −1 + 4 Rm − 2a m −1 Rm
2.31
+ k m h 4 − Pxh 2 (2 − a m −1 )

J 1 = 2hS t 2.32

M t h2
J2 = 2.33
Rt
2 Pt h 3
J3 = 2.34
Rt
Draft on Laterally Loaded Piles M. Budhu April, 2002

h Mt
J4 = 2.35
2 Rt S t

− Pxh 2
E= 2.36
Rt

Rt
3
( yt −2 − 2 yt −1 + 2 yt +1 − yt + 2 ) + Px ( yt −1 − yt +1 ) = Pt 2.37
2h 2h

Rt
( yt −1 − 2 yt + yt +1 ) = M t 2.38
h2

Q2
yt = 2.39
Q1

J 2 + G1 yt − d t −1
y t +1 = 2.40
G2

at y t +1 − y t + d t
yt + 2 = 2.41
bt

G1 H 2 ⎛ G ⎞1
Q1 = H 1 + + ⎜⎜1 − at 1 ⎟⎟ 2.42
G2 ⎝ G2 ⎠ bt

at (J 2 − d t −1 ) H 2 (d t −1 − J 2 ) d t
Q2 = J 3 + + + + d t −1 (2 + E + at − 2 ) − d t − 2 2.43
bt G2 G2 bt

G1 = 2 − at −1 2.44

G2 = 1 − bt −1 2.45

H 1 = −2a t −1 − Eat −1 − bt − 2 + at −1 at − 2 2.46

H 2 = −a t − 2 bt −1 + 2bt −1 + 2 + E (1 + bt −1 ) 2.47

y t −1 − y t +1 = J 1 2.48

Q4
yt = 2.49
Q3
Draft on Laterally Loaded Piles M. Budhu April, 2002

at −1 y t − J 1 + d t −1
yt +1 = 2.50
G4

at y t +1 − y t + d t
yt + 2 = 2.51
bt

H 2 at −1 at at −1 1
Q3 = H 1 + − + 2.52
G4 bt G4 bt

J 1 H 2 at ( J 1 − d t − 1) − G4 d t + bt d t −1 H 2
Q4 = J 3 + − + d t −1 (2 + E − at − 2 ) − d t − 2 2.53
G4 bt G4

G4 = 1 + bt −1 2.54

yt −1 − 2 y t + y t +1
= J4 2.55
y t −1 − yt +1

at d t −1 (1 − J 4 ) dt d H (1 − J 4 )
J3 − + + d t −1 (2 + E − at − 2 ) − d t − 2 + t −1 2
bt (G2 + J 4 G4 ) bt G2 + J 4 G4
yt = 2.56
a 1
H1 + H 2 H 3 − t H 3 +
bt bt

y t (G1 + J 4 at −1 ) − d t −1 (1 − J 4 ) d (1 − J 4 )
y t +1 = = H 3 y t − t −1 2.57
G2 + J 4 G4 G2 + J 4 G4

1
yt +2 = (at y t +1 − y t + d t ) 2.58
bt

G1 + J 4 at −1
H3 = 2.59
G2 + J 4 G4

y t = y −t 2.60