Axial Pile Response
Axial Pile Response
Outline
The previous lectures on piled foundations covered different forms of offshore pile
construction (driven, drilled and grouted, grouted driven), and the main methods of
estimating axial capacity. Offshore design of piles is dominated by design criteria
in terms of capacity, with settlement rarely being critical (apart from estimation of
foundation stiffness, to enable assessment of the dynamic response of the structure).
However, it is still necessary to consider the load-displacement response of pile
foundations, in order to assess the potential for degradation of shaft friction due to
progressive failure or cyclic loading.
Element of pile
Force per unit length, P
Ultimate
force/unit
length y P / d p net
Elastic range: ≈ =
d 4G 4G
Displacement, y
The schematic shows the main features of load transfer springs, and also indicates
the relationships between the initial spring stiffness and the shear modulus, G, of the
soil continuum. It is important that the gradient of the load transfer curve is
expressed in units of modulus (e.g. stress divided by displacement/diameter; or
force per unit length of pile divided by displacement), in order to avoid problems of
different scales (or size of pile).
In particular, for laterally loaded piles, it is common to refer to the horizontal load
transfer as P (load per unit length), rather than an average (net) pressure, p, given by
p = P/d where d is the pile diameter.
Axial loading
M
H dH dM d2y M
y ≈ −P ; = H; =
dz dz dz 2 (EI )p
dz P
d4y P
Hence =−
dz 4 (EI )p
Lateral loading
Actual modelling of the pile response is normally achieved through a computer code
that implements the beam-column equations for axial or lateral loading. The
governing differential equations are standard for structural engineering, and a
variety of codes are available in the industry.
l l/2 db Ep
η = λ=
d Gl
db l For solid cylindrical pile :
Ep =
(EA )steel
Depth πd 2 / 4
2η 2π tanh µl l
++ ρ
P (1 − ν )ξ ζ µl d
=
G l dw 1 8 η tanh µl l
1+
πλ (1 − ν ) ξ µl d
Offshore Foundation Systems 406
The final solution for the elastic response of a compressible pile is given above.
This solution may be incorporated into spreadsheet calculations of pile response
under axial load.
The main load transfer parameter, ζ, is written as:
ζ = ln(rm/ro)
where rm (the maximum radius of influence of the pile) is expressed as
rm = {0.25 + [2.5(1 - ν) ρ - 0.25](Gl/Gb)}l
or, for floating piles (Gl/Gb = 1)
rm = 2.5(1 - ν) ρ l
The numerator of the expression for the pile stiffness is made up of the base
stiffness (first term) and shaft stiffness (second term). The second term in the
denominator is generally small (except for very compressible piles) and may often
be ignored.
The compressibility of the pile is controlled by the tanµl/µl term, where the product
µl is given by
µl = (8/ζλ)0.5 l/d
When µl is greater than about 2, the tanhµl terms becomes close to unity and the
whole expression may be simplified (see next page). At the other extreme, for stiff
piles tanµl/µl becomes unity (as µl approaches zero) and the expression reverts to
that given previously for a rigid pile.
40
P 3000
G d
Glwd 30 ξ= l =1 η= b =1
Gb d
1000
G avg
20 ρ= = 0.75
Gl
300
Ep
10 100 λ=
30 Gl
10
0
1 10 100
Pile slenderness ratio, l/d
The design chart provides a simple way of estimating the non-dimensional pile
stiffness for a range of stiffness ratios, λ = Ep/Gl, and slenderness ratios, l/d. The
typical range of stiffness ratio, λ, is between 100 and 1000, although lower values
are appropriate for rock sockets.
The division of piles into rigid, intermediate, and compressible, is achieved through
the parameter, µl, which is proportional to the slenderness ratio, l/d, divided by the
square root of the stiffness ratio, λ. Thus, for piles where
l/d < 0.25λ0.5
the pile may be considered as essentially rigid.
On the other hand, for piles where
l/d > 1.5λ0.5
the pile is very compressible and the response becomes independent of the overall
pile length. For this case, the pile head stiffness may be estimated from:
P 1
≈ πρd E p G l ≈ ρd E p G l
w 2ς
For consistency, the shear modulus, Gl, should be replaced by the shear modulus at
a depth corresponding to z = 1.5dλ0.5, which is essentially the depth limit to which
any significant load transfer occurs.
τd
Shear stress, τ w=ς
2G
Single pile: ς ≈ 4
n
Group of n piles: ς* = nς − ∑ ln( 2si / d ) )
Displacement, w 2
1 − ν Pb
wb = ξ
Base load, Pb 2 db Gb
Single pile: ξ = 1
1 n d
Group of n piles: ξ* = ξ 1 + ∑ b
π 2 s i
Displacement, wb
Offshore Foundation Systems 406
The same solution for a single pile can be used to estimate the stiffness of a group of
piles, making use of modified load transfer parameters that make allowance for
interaction between neighbouring piles. Interactive effects can lead to a reduction in
effective stiffness of each pile by a factor of 2 or 3.
w τ/τs
0.025 m 1 Sand
The American Petroleum Institute (API) guidelines suggest a default load transfer
curve for axial analysis. This curve gives an approximately parabolic shaped curve,
reaching peak shear stress at a displacement of 1 % of the pile diameter. Beyond
the peak, a degree of strain softening is assumed, with the shear transfer reducing by
up to 30 % by a displacement of 2 % of the pile diameter.
In practice, strain softening for large diameter piles could continue over
significantly larger displacements, with much lower residual values of load transfer
(particularly under cyclic loading).
τ
B
C
w
Strain softening effects can lead to progressive failure of long, slender piles, such
that the actual capacity is less than the ideal (rigid pile) capacity. A simple design
chart may be constructed, showing the potential reduction in capacity as a function
of the relative compressibility of the pile. The non-dimensional factor, K, is
essentially the ratio of elastic shortening of the pile under the maximum shaft
capacity, to the displacement needed to strain soften between peak and residual.
Cyclic loading can lead to significant degradation of shaft friction. At the element
level (modelling one particular point down the pile), 1-way cyclic loading will tend
to lead to increasing average displacements, followed by failure at a value of shaft
friction that is lower than obtained under monotonic loading.
Under 2-way cyclic loading, the average displacement may not change much, but
the cyclic component will increase, followed by failure at a reduced level of shear
stress.
200
Shear stress Monotonic test
(kPa) 150
100
0
0 5 10 15 20
Displacement (mm)
Offshore Foundation Systems 406
A 1-way cyclic loading test on a model pile is shown here, and compared with the
response under monotonic loading. Failure occurs once the cumulative
displacement under cyclic loading reaches the monotonic loading envelope.
100
Monotonic response
75
Relative shear
stress (%) 50
25 2-way cyclic
1-way cyclic loading
loading 0
0 0.01 0.02 0.03 0.04 0.05 0.06
-25
-50
Displacement/Diameter
The load transfer code, RATZ, has in-built algorithms to model degradation under
cyclic loading. The plastic displacement undergone within each cycle is treated as
similar to additional monotonic displacement, with resulting degradation of the
available shaft friction.
Interaction Diagram
Normalised 1
cyclic shear
stress 0.8
Gerber Monotonic failure
0.6 (parabolic)
0.4
0.2 Goodman
(linear)
0
0 0.2 0.4 0.6 0.8 1
Normalised mean shear stress
500
200
100
0
0 5 10 15 20
-100
Displacement (mm)
Offshore Foundation Systems 406
This shows a typical 'element' response, obtained from a test on a model grouted
pile installed in a sample of calcarenite (so-called 'rod shear test'). The large
amount of strain softening is noteworthy, as is the very low cyclic resistance.
Two key sets of field tests were undertaken on grouted piles, in order to establish
design rules for axial response, and to provide data for calibrating load transfer
codes such as RATZ.
35 0.44 x 2.4 m
40
Depth 2.08 x 5.1 m 0.95 x 10.0 m
(m) 0.41 x 10.0 m
45
0.44 x 15.6 m
Grouted Driven
Grouted Section Tests Pile Tests
50
A total of 6 drilled and grouted pile tests were undertaken, 5 on short sections of
different diameter, and 1 on a scaled 'long' pile representative of a true prototype. In
the separate study by Esso Australia, 3 grouted driven piles were tested.
The zone of testing covered a weak limestone with strength and deformation
characteristics similar to those on the North-West Shelf.
An extensive site investigation was carried out in order to characterise the site, with
similar types of testing to an offshore site.
100
0
0 0.2 0.4 0.6 0.8 1
Displacement/diameter (%)
Offshore Foundation Systems 406
A consistent trend of decreasing peak shaft friction with increasing diameter of pile
was noted. This arises due to the decreasing effect of dilation at larger diameters, as
the normal stress change around the pile is proportion to the amount of radial
movement as a ratio of the pile diameter. For a given degree of roughness, the
radial movement would be similar for each pile size, and hence the ratio would
decrease as the diameter increased.
Load
(MN) 1
RATZ simulation
0.5
0
0 50 100 150 200
-0.5
Displacement (mm)
Offshore Foundation Systems 406
The monotonic tests on the short elements were used to calibrate the RATZ model
for application to the long (model) pile, which was subjected to cyclic loading.
Note the simulation of the very low cyclic shear stress.
1
0
0 10 20 30 40 50 60 70
Displacement (mm)
Offshore Foundation Systems 406
35
Depth
(m)
40
RATZ
simulation
45 Field
data
50
The profiles of load down the pile shows the gradual transfer of load to the lower
part of the pile, and also a reduction in shear transfer in the upper elements.