Pile Lateral Load Analysis Using Finite Difference Method Latpile Function

Pile Lateral Load Analysis Using Finite Difference Method
LatPile Function
Pile bending moments, forces and deflection due to lateral loads using a finite difference analysis
=LatPile(Segments, Loads, Soil table, Sections table, PY tables, Defaults, Output code, Koption)
Optional input in italics
Input:
Segments 2 column range listing pile length and number of segments

Loads 7 row single column range listing applied loads at head of pile:
See below for details of load input
Soil table Range listing soil properties, 7-9 rows x 1 column for each soil type
Effective depths at the top of each layer may optionally be entered in rows 8 and 9
See SoilProps sheet for further details of how Effective Depths are calculated.
Sections table Range listing pile section properties, 3 to 13 rows x 1 column for each section type
Notes: If the depth specified for Soil types and section types does not extend to the bottom
then the last provided value is applied to nodes to the bottom of the pile
Section stiffness (EI) may be specified as constant in row 3, or calculated for reinforc
piles of circular or rectangular cross section, to Eurocode 2 method.
Optional input:
PY tables Contiguous range with 2 columns for each table: pile deflection(Y) and support force per u
PY curves are used for those soil layers with "PY" in row 2 of the soil table, and Table Num
Defaults Change system defaults; available parameters and default values are listed below.
Output Code options 1: Output array includes all 6 output columns (see below for details)
2: Output Pile EI and soil stiffness only
KOption 1: initial stiffness (P/Y) for soil type 2 and 3 = Ki * Depth (as COM625 Manual and LPile)
2: initial stiffness (P/Y) for soil type 2 and 3 = Ki * Diam
Sign Convention Assuming a vertical pile, with nodes listed from top to bottom:
Deflections and loads are positive from left to right
Positive bending moments cause compression in the right side of the pile
Units All units are kN and m except:

Reinforcement bar diameter and cover (mm)
Concrete Shrinkage (microstrain)
Any other consistent units may be used if a fixed value is entered for the pile flexural stiffn
Non-linear pile behaviour The calculated pile flexural stiffness takes account of post-cracking behaviour of reinforced
concrete, including tension stiffening and shrinkage effects if a C (circular) or R (rectangula
is entered in row 3 of the Sections Table. If a numerical value is entered in row 3 this is us
fixed stiffness (EI) and the entries below are ignored.
Note that non-linearity due to steel yielding, or concrete crushing is not included in
A separate check must be made of steel and concrete stresses.
Output:
Output is a 6 column array, enter as an array function:
To enter an array function: Enter the function in the top left cell of the output range
Select the entire output range
Press F2
Press Ctrl-shift-enter
Column 1 Node Position from top, m
Column 2 Deflection, m
Column 3 Bending Moment, kNm
Column 4 Shear force, KN
Column 5 Pile EI, kN/m2
Column 6 Soil spring stiffness, kN/m2
Other Functions
PY624tab Function
Table of soil stiffness values at a specified effective depth
=PY624tab(Soil data, Pile Sections, Loads, Deflections, Koption)
Input:
Soil data As LatPile function, except that rows 8 and 9 are not optional
Pile Sections As LatPile function
Loads As LatPile function
Deflections Single column range listing pile deflection in ascending order
Optional input:
KOption As LatPile function
Output:
Output is a single column array , enter as an array function:
Spring force/m of pile for each input deflection
PfEffDepth Function
Returns the maximum soil resistance per metre due to wedge failure at a specified effective depth
=PfEffDepth(Soil data, Pile EffDepth, PYTabsA)
Input:
Soil data As LatPile function, except that row 1 is pile diameter and rows 8 and 9 are not required
Eff Depth Effective depth at which the soil stiffness is required
PYTabsA As LatPile function
Output:
Maximum Spring force/m due to wedge failure
entered in rows 8 and 9
ths are calculated.
mn for each section type
does not extend to the bottom of the pile
bottom of the pile
row 3, or calculated for reinforced concrete
ode 2 method.
tion(Y) and support force per unit length (P).

of the soil table, and Table Number in row 3.
values are listed below.
COM625 Manual and LPile)
side of the pile
ntered for the pile flexural stiffness.
racking behaviour of reinforced

if a C (circular) or R (rectangular)
ue is entered in row 3 this is used as a
e crushing is not included in the analysis.
t cell of the output range

ows 8 and 9 are not required
Examples
Example 1: Example 2 from COM624 Manual

Layers Loads
Length Subdivisions Type
500 50 Shear
Moment
Axial at head
Axial at base
Depth of zero axial
Number of cycles
Water table depth
Shear and moment are specified at the head of the pile. There is currently no provision for
applying loads at other locations.
Axial load is assumed to vary linearly from the head to toe of the pile, or bi-linearly from the
head to the depth of zero axial, and depth of zero axial to the base.
If Depth of zero axial is <=0, >= pile length, or blank, it is ignored.
Soil Tables
Depth of layer base 20
Soil cohesion, Cu; or "PY" to use PY table 0.00
Soil friction angle, Phi, or PY Table number 0
Soil Density 0.000
Strain at 50% ultimate stress, Eps50 0.0050
Soil Type 0
Initial stiffness, ki 0
Layer top Effective depth for wedge failure (optional)
Layer top Effective depth for initial stiffness (optional) 0
Soil Type, 0 = "void", 1 = "Soft to firm clay", 2 = "Stiff Clay", 3 = "Sand"
PY Tables
Y P
0.000E+00 0.000E+00
4.800E-01 9.568E+02
9.600E-01 1.138E+03
1.440E+00 1.259E+03
1.920E+00 1.353E+03
2.400E+00 1.431E+03
2.880E+00 1.497E+03
3.360E+00 1.556E+03
3.840E+00 1.609E+03
4.320E+00 1.657E+03
4.800E+00 1.701E+03
5.280E+00 1.742E+03
5.760E+00 1.781E+03
6.240E+00 1.817E+03
6.720E+00 1.851E+03
7.200E+00 1.883E+03
9.000E+00 1.883E+03
Page 5
Examples
Sections Table
Depth of base of section ii
Pile diam or width in
EI or Section Type (C or R) for RC lb-in2
Bar diameter Diat }
Number of bars* Numt }
Bar cover Covt }
Concrete compressive strength f'c }
Steel yield stress fys }N/A with non-metric
E steel Es } units
E concrete Ec }
Concrete flex. tens. strength f'ct }
Concrete shrinkage strain Epscs }
Section depth (rectangular sect.) }
* Total number for circular section, or bars on tensile face for rectangular section
Defaults Table (default values in brackets)

Tension stiffening Beta value (0.5)
Maximum iterations (20)
Tolerance on deflection norm (0.001)
Tolerance on soil stiffness norm (0.001)
Tolerance on section stiffness norm (0.001)
Output
X Deflection Moment
#NAME? #NAME? #NAME?
Page 6
Examples

Soil Resistance
-600.000 -400.000 -200.000 0.000 200.000 400.000 600.000 800.000
0.000
-50.000
-100.000
Depth
-150.000
-200.000
-250.000
-300.000
LatPile COM624
Bending Moment
0.000
0.00E+00 5.00E+05 1.00E+06 1.50E+06 2.00E+06 2.50E+06 3.00E+06 3.50E+06
-20.000
-40.000 Page 7
-60.000
Bending Moment
Examples
0.000
0.00E+00 5.00E+05 1.00E+06 1.50E+06 2.00E+06 2.50E+06 3.00E+06 3.50E+06
-20.000
-40.000
-60.000
-80.000
-100.000
-120.000
-140.000
-160.000
-180.000
LatPile COM624
Page 8
Examples
Loads
lb, in-lb
50750
110833.333333333
100000
100000
1
600
ently no provision for
or bi-linearly from the
80 200 260 404 600

13.90 0.00 5.56 20.83 100.00
0 30 0 0 0
0.072 0.035 0.035 0.035 0.035
0.0050 0.0200 0.0050 0.0010
2 3 1 2 2
1000 60 200 1200 2000
0 160 278 342 404
0 60 180 240 384
Y P Y P Y P
0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00
5.000E-02 3.710E+02 1.400E-01 1.800E+02 4.800E-01 2.790E+03
1.000E-01 7.420E+02 4.500E-01 5.670E+02 9.600E-01 3.317E+03
1.500E-01 1.094E+03 9.000E-01 7.150E+02 1.440E+00 3.671E+03
2.000E-01 1.216E+03 1.350E+00 8.180E+02 1.920E+00 3.944E+03
2.500E-01 1.321E+03 1.800E+00 9.010E+02 2.400E+00 4.171E+03
3.000E-01 1.413E+03 2.250E+00 9.700E+02 2.880E+00 4.365E+03
3.500E-01 1.496E+03 2.700E+00 1.031E+03 3.360E+00 4.537E+03
4.000E-01 1.572E+03 3.150E+00 1.085E+03 3.840E+00 4.691E+03
4.500E-01 1.642E+03 3.600E+00 1.134E+03 4.320E+00 4.831E+03
5.000E-01 1.707E+03 4.050E+00 1.180E+03 4.800E+00 4.960E+03
5.500E-01 1.768E+03 4.500E+00 1.222E+03 5.280E+00 5.080E+03
6.000E-01 1.826E+03 4.950E+00 1.261E+03 5.760E+00 5.191E+03
1.350E+00 2.669E+03 5.400E+00 1.299E+03 6.240E+00 5.296E+03
1.500E+00 2.669E+03 1.440E+01 1.801E+03 6.720E+00 5.395E+03
7.200E+00 5.489E+03
9.000E+00 5.489E+03
Page 9
Examples
300.0 600.0
36.000 36.0
2.56E+11 1.91E+11
0.5
100
0.001
0.001
0.001
Spring
Shear Force Pile EI stiffness layer EffDepth1 EffDepth2
#NAME? #NAME? #NAME? #NAME? #NAME? #NAME?
Page 10
Examples

Deflection
000 600.000 800.000 0.000 5.000 10.000 15.000 20.000 25.000 30.000 35
0.000
-20.000
-40.000
Depth
-60.000
-80.000
-100.000
-120.000
-140.000
LatPile COM624
Shear Force
0.000
50E+06 3.00E+06 3.50E+06 -30000.0 -20000.0 -10000.0 0.0 10000.0 20000.0 30000.0 40000.0 5000
-50.000
Page 11
-100.000
Shear Force
Examples
0.000
50E+06 3.00E+06 3.50E+06 -30000.0 -20000.0 -10000.0 0.0 10000.0 20000.0 30000.0 40000.0 5000
-50.000
-100.000
-150.000
-200.000
-250.000
-300.000
LatPile COM624
Page 12
Examples
pf Spring force Spring force Depth

#NAME? #NAME? #NAME? 0
#NAME? #NAME? #NAME? -10
Page 13
Examples

25.000 30.000 35.000
M624
0.0 30000.0 40000.0 50000.0 60000.0
Page 14
Examples
0.0 30000.0 40000.0 50000.0 60000.0
M624
Page 15
Examples
m 36.000 36.000 36.000 36.000 36.000

kPa 0.00 13.90 0.00 5.56 20.83
degrees 0.00 0.00 30.00 0.00 0.00
kN/m3 0.00 0.07 0.03 0.03 0.03
0.01 0.01 0.00 0.02 0.01
1 to 3 0.00 2.00 3.00 1.00 2.00
KN/m3 0.00 1000.00 60.00 200.00 1200.00
Page 16
Examples
Deflection, in Moment, in-lb Shear, lb Soil Resist, lb/in

0.251 1.00E+05 50200 0.00E+00
0.23 6.05E+05 50700 0.00E+00
0.209 1.12E+06 50700 6.02E+02
0.189 1.60E+06 44900 6.43E+02
0.17 2.02E+06 38400 6.63E+02
0.151 2.37E+06 31600 6.80E+02
0.133 2.66E+06 24800 6.95E+02
0.117 2.87E+06 17800 7.06E+02
0.101 3.02E+06 11200 5.07E+02
0.0865 3.10E+06 6210 4.86E+02
0.0733 3.14E+06 1490 4.56E+02
0.0613 3.14E+06 -2880 4.18E+02
0.0505 3.09E+06 -6840 3.75E+02
0.0409 3.00E+06 -10400 3.28E+02
Page 17
Examples
2.89E+06 -13400 2.81E+02

2.74E+06 -16000 2.33E+02
2.57E+06 -18100 1.87E+02
-19700 1.45E+02
-21000 1.06E+02
-21900 7.52E+01
-22600 1.13E+02
-23600 8.72E+01
-24300 4.75E+01
-24100 -6.29E+01
-23400 -7.62E+01
-22600 -7.99E+01
-4.82E+02
-4.76E+02
-4.57E+02
Page 18
Examples
36.000
100.00
0.00
0.03
0.00
2.00
2000.00
Page 19
Recommended Typical Soil Properties
Source: COM624 Manual
Strain at 50% ultimate stress, Eps50
Average
Clay Consistency Undrained Shear Eps50
0 Strength, kPa
Soft 10 to 25 0.02
Medium 25 to 50 0.01
Stiff 50 to 100 0.007
100 to 200 0.005
200 to 400 0.004
Initial soil stiffness, Ki

0 100 200
Average Undrained Shear Strength, kPa
50 to 100 100 to 200 200 to 400
Stiff Clay
Static load, kN/m3 140,000 270,000 540,000 140000
Cyclic load, kN/m3 50,000 110,000 220,000 50000
Sand (static or cyclic) Loose Medium Dense

Submerged kN/m3 5,000 15,000 35,000
Above water table kN/m3 7,000 25,000 60,000
Unit conversion: 1 ton/ft2 = 95.8 kPa

1 lb/in3 = 272 kN/m3
Effective Depth
When a pile with a layered soil is analysed, the upper soil layers affect the lower layers in three respects:
If the ultimate resistance is controlled by wedge sliding, the resistance at any level is the sum of the
resistances from each layer above.
If the ultimate resistance is controlled by flow with a granular soil the reisitance is related to the vertical
pressure at the level of analysis.
The initial stiffness for type 2 and 3 soils is proportional to the vertical pressure at the level of analysis.
To allow for these effects the LatPile function calculates two effective depths:
1) The effective depth for wedge failure for a point within any layer, L, has two components:
i) The depth of soil with the same properties as Layer L that would have the same wedge sliding
resistance as the actual soil layers above the top of Layer L.
ii) The depth of the point below the top of Layer L
2) The effective depth for vertical pressure for a point within a layer, L, is the depth of material that
would have the same vertical pressure as that due to the actual soils layers above the point.
Depths 1) and 2) are calculated by the LatPile function, or may optionally be specified for the top
of each layer as rows 8 and 9 of the Soil Properties Table.
Example PY Curves; Data from Examples sheet
Soil Tables
Depth of point 80 100 220 340
Soil cohesion, Cu; or "PY" to use PY table 13.90 0.00 5.56 20.83
Soil friction angle, Phi, or PY Table number 0 30 0 0
Soil Density 0.072 0.035 0.035 0.035
Strain at 50% ultimate stress, Eps50 0.0050 0.0200 0.0050
Soil Type 2 3 1 2
Initial stiffness, ki 1000 60 200 1200
Effective depth for wedge failure 40 190 298 278
Effective depth for initial stiffness 60 120 240 360
Depth = 40 in below surface COM 624 Output PY624tab

Deflection Force Force 2000.0
in lb/in lb/in
1800.0
0 0.0 #NAME?
0.48 956.8 #NAME? 1600.0
0.96 1137.8 #NAME?
1.44 1259.2 #NAME? 1400.0
1.92 1353.0 #NAME? 1200.0
2.4 1431.0 #NAME?
2.88 1497.0 #NAME? 1000.0
3.36 1556.0 #NAME? 800.0
3.84 1609.0 #NAME?
4.32 1657.0 #NAME? 600.0
4.8 1701.0 #NAME? 400.0
5.28 1742.0 #NAME?
5.76 1781.0 #NAME? 200.0
6.24 1817.0 #NAME?
0.0
6.72 1851.0 #NAME?
0 1 2 3
7.2 1883.0 #NAME?
9 1883.0 #NAME?

Deflection Force Force
in lb/in lb/in
0 0.0 #NAME? 3000.0
0.05 371.0 #NAME?
0.1 742.0 #NAME?
0.15 1094.0 #NAME? 2500.0
0.2 1216.0 #NAME?
0.25 1321.0 #NAME? 2000.0
0.3 1413.0 #NAME?
0.35 1496.0 #NAME?
1500.0
0.4 1572.0 #NAME?
0.45 1642.0 #NAME?
0.5 1707.0 #NAME? 1000.0
0.55 1768.0 #NAME?
0.6 1826.0 #NAME? 500.0
1.35 2669.0 #NAME?
1.5 2669.0 #NAME?
0.0
0 0.2 0.4

2000.0
in lb/in lb/in 1800.0
0 0.0 #NAME? 1600.0
0.14 180.0 #NAME?
0.45 567.0 #NAME? 1400.0
1200.0
1000.0
800.0
1800.0
1600.0
1400.0
0.9 715.0 #NAME? 1200.0
1.35 818.0 #NAME?
1000.0
1.8 901.0 #NAME?
2.25 970.0 #NAME? 800.0
2.7 1031.0 #NAME?
600.0
3.15 1085.0 #NAME?
3.6 1134.0 #NAME? 400.0
4.05 1180.0 #NAME? 200.0
4.5 1222.0 #NAME?
4.95 1261.0 #NAME? 0.0
5.4 1299.0 #NAME? 0 2 4
14.4 1801.0 #NAME?

in lb/in lb/in
6000.0
0 0.0 #NAME?
0.48 2790.0 #NAME?
0.96 3317.0 #NAME? 5000.0
1.44 3671.0 #NAME?
1.92 3944.0 #NAME?
4000.0
2.4 4171.0 #NAME?
2.88 4365.0 #NAME?
3.36 4537.0 #NAME? 3000.0
3.84 4691.0 #NAME?
4.32 4831.0 #NAME? 2000.0
4.8 4960.0 #NAME?
5.28 5080.0 #NAME?
5.76 5191.0 #NAME? 1000.0
6.24 5296.0 #NAME?
6.72 5395.0 #NAME? 0.0
7.2 5489.0 #NAME? 0 1 2 3
9 5489.0 #NAME?
the vertical
0.0
0.0
0.0
0.0
0.0
0.0 COM624
PY624Tab
0.0
0.0
0.0
0.0
0.0
0 1 2 3 4 5 6 7 8 9 10
000.0
500.0
000.0
500.0 COM624
PY624Tab
000.0
500.0
0.0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6
000.0
800.0
600.0
400.0
200.0
000.0 COM624
PY624Tab
800.0
800.0
600.0
400.0
200.0
000.0 COM624
PY624Tab
800.0
600.0
400.0
200.0
0.0
0 2 4 6 8 10 12 14 16
000.0
000.0
000.0
000.0 COM624
PY624Tab
000.0
000.0
0.0
0 1 2 3 4 5 6 7 8 9 10
12
10
0
0 2
12
10
6 Column AL
Column AM
4
0
0 2 4 6 8 10 12
Disclaimer
This program is free software; you can redistribute it and/or
modify it under the terms of the GNU General Public License
as published by the Free Software Foundation; either version 2
of the License, or (at your option) any later version.
This program is distributed in the hope that it will be useful,

but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
A copy of the GNU General Public License may be obtained from:

The Free Software Foundation, Inc.
59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
[email protected]
www.interactiveds.com.au
http://newtonexcelbach.wordpress.com/
ã Copyright 2012 Interactive Design Services Pty Ltd. all rights reserved
Revision History
Current Version 1.04 4-Jul-12
Rev Date By Description
1.00 17-Jul-10 DAJ First release

1.01 11-Jan-11 DAJ Code to read first row in input PY tables corrected (row and column indeces reversed)
1.02 21-Jan-12 DAJ KOption added; PY624tab function added; PY_624 function modified, effective depth inp
1.03 26-Jan-12 DAJ PY curve corrected for Soil type 1, 3*Eps50 < Y < 15*Eps50 and Type 2 below water table
1.04 4-Jul-12 DAJ Effective depth calculation code modified so that blank rows 8 and 9 in the soil table are ig
mn indeces reversed)
dified, effective depth input added
Type 2 below water table, As*Y50 < Y < 6*As*Y50
nd 9 in the soil table are ignored.

Pile Lateral Load Analysis Using Finite Difference Method Latpile Function

Uploaded by

Copyright:

Available Formats

Pile Lateral Load Analysis Using Finite Difference Method Latpile Function

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Pile Lateral Load Analysis Using Finite Difference Method Latpile Function

Uploaded by

Copyright:

Available Formats

Pile Lateral Load Analysis Using Finite Difference Method

Segments 2 column range listing pile length and number of segments

Units All units are kN and m except:

KOption As LatPile function

tion(Y) and support force per unit length (P).

COM625 Manual and LPile)

side of the pile

ntered for the pile flexural stiffness.

racking behaviour of reinforced

e crushing is not included in the analysis.

t cell of the output range

Example 1: Example 2 from COM624 Manual

Defaults Table (default values in brackets)

#NAME? #NAME? #NAME?

ently no provision for

or bi-linearly from the

80 200 260 404 600

#NAME? #NAME? #NAME? #NAME? #NAME? #NAME?

pf Spring force Spring force Depth

#NAME? #NAME? #NAME? -140

25.000 30.000 35.000

0.0 30000.0 40000.0 50000.0 60000.0

0.0 30000.0 40000.0 50000.0 60000.0

m 36.000 36.000 36.000 36.000 36.000

Deflection, in Moment, in-lb Shear, lb Soil Resist, lb/in

2.89E+06 -13400 2.81E+02

Strain at 50% ultimate stress, Eps50

Initial soil stiffness, Ki

Sand (static or cyclic) Loose Medium Dense

Unit conversion: 1 ton/ft2 = 95.8 kPa

Example PY Curves; Data from Examples sheet

Depth = 40 in below surface COM 624 Output PY624tab

Depth = 100 in below surface COM 624 Output PY624tab

Depth = 220 in below surface COM 624 Output PY624tab

Depth = 340 in below surface COM 624 Output PY624tab

This program is distributed in the hope that it will be useful,

A copy of the GNU General Public License may be obtained from:

Current Version 1.04 4-Jul-12

Rev Date By Description

1.00 17-Jul-10 DAJ First release

You might also like