Bridge Software Institute 2

Similar documents
LATERALLY LOADED PILE GROUPS

Modeling Foundations in RS

SUBMERGED CONSTRUCTION OF AN EXCAVATION

RSPile. Tutorial 3 Grouped Pile Analysis. Pile Analysis Software. Grouped Pile Analysis

Moment-rotation Behavior of Shallow Foundations with Fixed Vertical Load Using PLAXIS 3D

PLAXIS 2D - SUBMERGED CONSTRUCTION OF AN EXCAVATION

PLAXIS 3D Benchmark for Suction Anchor Bearing Capacity. Richard Witasse Plaxis bv

OP-029. INTERPRETATION OF CIU TEST

Simplified Design Procedure for Piles Affected by Lateral Spreading based on 3D Nonlinear FEA Using OpenSees. University of Washington

Manual. Ansys Exercise. Offshore Wind Farm Design

Settlement Analysis of Piled Raft Foundation System (PRFS) in Soft Clay

A Sliding Block Analysis

MODELLING OF PILES. Modelling of Piles in USFOS

Stability Analysis of Cellular Retaining Structure by PLAXIS Finite Element Code

Lateral Loading of Suction Pile in 3D

OPENSEES Soil-Pile Interaction Study under Lateral Spread Loading

FOUNDATION IN OVERCONSOLIDATED CLAY

New Constitutive Material Models in RS 2

PLAXIS 3D. Tutorial Manual

SETTLEMENT OF A CIRCULAR FOOTING ON SAND

EFFECT OF DREDGING AND TIE-ROD ANCHOR ON THE BEHAVIOR OF BERTHING STRUCTURE

PLAXIS 3D. Tutorial Manual

New developments in numerical modelling of pile installation

17. SEISMIC ANALYSIS MODELING TO SATISFY BUILDING CODES

PLAXIS 2D. Tutorial Manual

FB-MULTIPIER vs ADINA VALIDATION MODELING

CHAPTER 4. Numerical Models. descriptions of the boundary conditions, element types, validation, and the force

PySimple1Gen OpenSees Command

DRY EXCAVATION USING A TIE BACK WALL

Capacity Analysis of Suction Anchors in Clay by Plaxis 3D Foundation

Investigation of the behaviour of single span reinforced concrete historic bridges by using the finite element method

Release notes. D-Geo Pipeline New feature

Analysis and Evaluation of Single Piles in Laterally Spreading Soil

Dam Construction by Stages

Simulations of tunnel excavation in 2D and 3D conditions considering building loads

PHASED EXCAVATION OF A SHIELD TUNNEL

THERMAL EXPANSION OF A NAVIGABLE LOCK

Modelling of Clay Behaviour in Pile Loading Test using One-Gravity Small-Scale Physical Model

Global to Local Model Interface for Deepwater Top Tension Risers

A load-displacement based approach to assess the bearing capacity and deformations of mono-bucket foundations

Terrain settlement analysis

Non-Linear Analysis of Bolted Flush End-Plate Steel Beam-to-Column Connection Nur Ashikin Latip, Redzuan Abdulla

Soil-Structure Interaction for Piled-Raft Foundation

WP1 NUMERICAL BENCHMARK INVESTIGATION

GEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE

Finite Element Analysis of a 10 x 22 FRP Building for TRACOM Corporation

THE COMPUTATIONAL MODEL INFLUENCE ON THE NUMERICAL SIMULATION ACCURACY FOR FORMING ALLOY EN AW 5754

1D Consolidation Analysis & Design

BEARING CAPACITY OF CIRCULAR FOOTING

BEARING CAPACITY OF STRIP FOOTING

ANALYSIS OF BOX CULVERT - COST OPTIMIZATION FOR DIFFERENT ASPECT RATIOS OF CELL

Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction

SAFI Sample Projects. Design of a Steel Structure. SAFI Quality Software Inc. 3393, chemin Sainte-Foy Ste-Foy, Quebec, G1X 1S7 Canada

GEO-SLOPE International Ltd , 700-6th Ave SW, Calgary, AB, Canada T2P 0T8 Main: Fax:

Accuracy of the Rubicon Toolbox Finite Element Model

GEO-SLOPE International Ltd, Calgary, Alberta, Canada Sheet Pile Wall

Supplementary information

FREE VIBRATION AND EARTHQUAKE ANALYSIS OF A BUILDING

DYNAMIC ANALYSIS OF A GENERATOR ON AN ELASTIC FOUNDATION

CONSTRUCTION OF A ROAD EMBANKMENT

GEOSTUDIO Tutorials Results and Procedure Comparison

Simulation of AJWSP10033_FOLDED _ST_FR

CONSTITUTIVE MODELING OF GEOMATERIALS

Oasys Frew. Copyright Oasys 2013

Liquefaction Analysis in 3D based on Neural Network Algorithm

Eurocode Training EN 1997 in SCIA Engineer

Numerical Modelling of Engineering Structures subjected to Dynamic Seismic Loading

Distance Between Two Snaked-lay of Subsea Pipeline. Yuxiao Liu1, a *

Stability with FE Pore-Water Pressure

Application Note - Gravity Wall Analysis

RAPID DRAWDOWN ANALYSIS

FREE VIBRATION AND EARTHQUAKE ANALYSIS OF A BUILDING

9 FREE VIBRATION AND EARTHQUAKE ANALYSIS OF A BUILDING

3 SETTLEMENT OF A CIRCULAR FOOTING ON SAND (LESSON 1) Figure 3.1 Geometry of a circular footing on a sand layer

Modeling passive earth pressures on bridge abutments for nonlinear Seismic Soil - Structure interaction using Plaxis

LARGE DEFORMATION FINITE ELEMENT ANALYSIS OF OFFSHORE PENETRATION TESTS

Design of diaphragm and sheet pile walls. D-Sheet Piling. User Manual

Prediction of One-Dimensional Compression Test using Finite Elements Model

Scientific Manual FEM-Design 17.0

Engineering Effects of Boundary Conditions (Fixtures and Temperatures) J.E. Akin, Rice University, Mechanical Engineering

Investigation of Traffic Load on the Buried Pipeline by Using of Real Scale Experiment and Plaxis-3D Software

Application Note - Reinforced Soil Analysis

A05 Steel Catenary Riser Systems

What makes Bolt Self-loosening Predictable?

UTEXAS4 A COMPUTER PROGRAM FOR SLOPE STABILITY CALCULATIONS. Stephen G. Wright. May Revised September Shinoak Software Austin, Texas

Anchored Retaining Walls in Granite Residual Soils II. A Method for Preliminary Design

A comparison between large-size shaking table test and numerical simulation results of subway station structure

( ) ( ) tutorial question 1. Given. γ soil := 21 kn m 3. x 2 = 3.3 m. at rest. k 0 = 0.5 lateral earth pressure co-eff. at 0.5 m below grade σ t u

Stress Analysis of Cross Groove Type Constant Velocity Joint

Internal Forces and Moments

OpenSees2DPS (Beta 0.1): 2D Plane Strain Analysis

FE-Analysis of piled and piled raft foundations. Jean-Sébastien LEBEAU

Global Journal of Engineering Science and Research Management

This is NOT a truss, this is a frame, consisting of beam elements. This changes several things

Engineeringmanuals. Part3

Chapter 20. Finite Element Method (FEM) Introduction

Analysis of Pile Behaviour due to Damped Vibration by Finite Element Method (FEM)

A Multiple Constraint Approach for Finite Element Analysis of Moment Frames with Radius-cut RBS Connections

Optimal Support Solution for a Meniscus Mirror Blank

General Applications

Transcription:

EXECUTIVE SUMMARY This report summarizes six API soil models that have been implemented into the FB- MultiPier computer program according to the American Petroleum Institute Recommended Practice 2A LRFD (API RP 2A LRFD). This report provides review of the theory, validation of the numerical models, and examples of modeling soil-structure interaction using the FB-MultiPier API soil models. As a further aid for the usage of FB- MultiPier, this report serves as a detailed user s guide to current users of the program. Bridge Software Institute 2

TABLE OF CONTENTS 1. API Soil Models...4 1.1 Axial load transfer (T-z) curves...4 1.2 Tip load - displacement (Q-z) curves...8 1.3 Axial soil resistance - deflection (P-y) curves...11 2. Validation of FB-MultiPier API Soil Models...16 2.1 Lateral soil response curves (P-y) of API sand soil model...16 2.2 Axial soil response curves (T-z, Q-z) of API sand soil model...17 2.3 Lateral soil response curves (P-y) of API clay soil model...19 2.4 Axial soil response curves (T-z, Q-z) of API clay soil model...20 3. Case Study and Comparison...22 3.1 Case 1: Lateral soil-structure interaction analysis using an API sand soil model.22 3.2 Case 2: Axial soil-structure interaction analysis using an API sand soil model...23 3.3 Case 3: Lateral soil-structure interaction analysis using an API clay soil model..25 3.4 Case 4: Axial soil-structure interaction analysis using an API clay soil model...27 3.5 Case 5: Lateral soil-structure interaction analysis using a multi-layered API soil model...29 3.6 Case 6: Combined lateral/axial soil-structure interaction analysis using a multilayered API sand soil model...30 4. Summary and Conclusion...31 5. References...32 Appendix A. Showcase FB-MultiPier Analysis...33 Appendix B. Application of the API Soil Models in FB-MultiPier...34 Bridge Software Institute 3

1. API Soil Models 1-1. Axial Load Transfer (T-z) Curves API Clay and API Sand Models (Refer to G.4.2, G.4.3, and G.7.2 API RP2A LRFD) The unit skin friction and the T-z curves for clay and sand soils can be determined as per Section G.4.2 and Section G.4.3 API RP2A LRFD respectively. Figure 1.1 shows the T-z curves for non-carbonate soils, recommended by API RP2A LRFD. For pipe piles in cohesive soils, the skin friction can be calculated by the equation. f c where c α = undrained shear strength of the soil in stress units = a dimensionless factor, which is defined as 0.5 0.5 1.0 for 1.0 0.25 0.5 1.0 for 1.0 c/ p' o where p'o = effective overburden pressure in stress units For clay soils the ratio of residual stress to ultimate stress (tres/tmax) ranges from 0.70 to 0.90 as shown in the Figure 1.1. In FB-MultiPier analysis, the tres/tmax ratio is set equal to 0.90. In FB-MultiPier analysis, the T-z curves for API clay are generated using piecewise linear function based on the data given in the table and is shown in Figure 1.2 z/d 0.0016 0.0031 0.0057 0.0080 0.0100 0.0200 t/tmax 0.30 0.50 0.75 0.90 1.00 0.90 0.90 Bridge Software Institute 4

Figure 1.1 Axial Pile Load Transfer- Displacement (T-z) Curves (source: API RP2A LRFD 1993) where Z = local pile deflection D = pile diameter t = mobilized soil pile adhesion in stress units tmax = maximum soil pile adhesion or unit skin friction capacity computed For pipe piles in cohesionless soils, the unit skin friction is calculated as f K p' o tan( ) where K = dimensionless coefficient of lateral earth pressure (ratio of horizontal to vertical normal effective stress, for unplugged K=0.8 and for plugged K=1.0) p'o = effective overburden pressure in stress units δ = friction angle between the soil and pile wall, which is defined as Bridge Software Institute 5

0 5 where φ = internal friction angle It is recommended that the ultimate (limiting) values of unit skin friction, fult, be considered, which are given in Table G.4.3-1 API RP2A LRFD. In FB-MultiPier analysis, the T-z curves for API sand are generated using piecewise linear function based on the data given in the table and shown in Figure 1.3 z (in) 0.00 0.10 t/tmax 0.00 1.00 1.00 1.2 1 0.8 t/tmax 0.6 0.4 0.2 0 0 0.01 0.02 0.03 0.04 z/d Figure 1.2 Normalized T-z curve for API Clay Bridge Software Institute 6

1.2 1 0.8 t/tmax 0.6 0.4 0.2 0 0 0.1 0.2 0.3 Notation used for axis labels is defined as t = mobilized skin friction z (in) Figure 1.3 Normalized T-z curve for API Sand tmax = maximum unit skin friction capacity (f) z = local pile deflection D = diameter of pile Bridge Software Institute 7

1-2. Tip Load-Displacement (Q-z) Curves API Clay and API Sand Models (Refer to G.4.2, G.4.3, and G.7.3 API RP2A LRFD) The unit end bearing capacity and Q-z curves for clay and sand soils can be determined as per Section G.4.2 and Section G.4.3 API RP2A LRFD respectively. A relatively large pile tip movement which is required to fully mobilize the end bearing resistance may be achieved by a pile tip displacement up to 10% of the pile diameter. Figure 1.4 shows the Q-z curves of both sand and clay soils as recommended in API RP2A LRFD. The unit end bearing of pipe pile founded in cohesive soil is given by q 9c where c = undrained shear strength of the soil, in stress units, at the pile tip The unit end bearing of pipe pile founded in cohesionless soil is given by q p N ' o q where p'o = effective overburden pressure, in stress units, at the pile tip Nq = dimensionless bearing capacity factor, which is defined as tan( ') 2 0 ' Nq e tan 45 2 where φ = effective internal friction angle at the pile tip The ultimate end bearing capacity is then calculated Qp q A where A = sectional area at the tip of pile, which is based on the pile end condition; plugged a gross sectional area is used to compute ultimate end bearing capacity unplugged a cross sectional area is used to compute ultimate end bearing capacity Bridge Software Institute 8

Figure 1.4 Pile Tip Load Displacement (Q-z) Curve (source: API RP2A LRFD 1993) where z = axial pile deflection D = pile diameter Q = mobilized end bearing capacity in force units Qp = total end bearing computed It is recommended that the ultimate (limiting) values of unit end bearing for cohesionless soils, qult, be considered, which are given in Table G.4.3-1 API RP2A LRFD. In FB-MultiPier analysis, Q-z curves of API clay and API sand are defined as a piecewise linear function based on the data given in the table and is shown in Figure 1.5. z/d 0.002 0.013 0.042 0.073 0.100 Q/Qp 0.25 0.50 0.75 0.90 1.00 1.00 Bridge Software Institute 9

1.2 1 0.8 Q/Qp 0.6 0.4 0.2 0 0 0.05 0.1 Figure 1.5 Normalized Q-z curve for API Clay and API Sand Notation used for axis labels is defined as Q = mobilized end bearing capacity Qp = total end bearing capacity z = axial deflection of pile tip D = diameter of pile z/d Bridge Software Institute 10

1-3. Lateral Soil Resistance - Deflection (P-y) Curves API Clay and API Sand Models (Refer to G.8.2 to G.8.7 API RP2A LRFD) The ultimate lateral soil resistance and P-y curves may be calculated using section G.8.2 to G.8.7 API RP2A LRFD. The soil resistance equations are not applicable if the strength variation along the depth of the soil is inconsistent. The ultimate unit lateral resistance, pu, of soft clay under static loading conditions can vary between 8c to 12c except at the shallow depths. In the absence of more definitive criteria, use the empirical equation given by API RP2A LFRD: where c p 3 c p ' J cz for z X D u o R p 9 c for z X u R = undrained shear strength of undisturbed clay soil samples in stress units p'o = effective overburden pressure in stress units z D J = depth below ground surface = diameter of the pile = dimensionless empirical constant XR = depth from the ground surface to the bottom of reduced resistance zone, which is defined as; X R 6 D ' D J c 2.5D where γ' = effective unit weight of soil in weight density units In FB-MultiPier analysis, the value of J is set equal to 0.5, which is recommended for Gulf of Mexico clays. The P-y curves of API clay under cyclic and static loading conditions are defined as a piecewise linear function based on the data given in the table and is shown in Figure 1.6 and Figure 1.7. Bridge Software Institute 11

Static loading z > XR Cyclic loading z < XR p/pu y/yc p/pu y/yc p/pu y/yc 0.00 0.50 0.72 1.00 1.00 0.0 1.0 3.0 8.0 0.00 0.50 0.72 0.72 0.0 1.0 3.0 0.00 0.50 0.72 0.72 z/xr 0.72 z/xr 0.0 1.0 3.0 15.0 1 P/Pu 0.5 0 0 5 10 15 y/y c Static loading Cyclic loading Figure 1.6 Normalized P-y curve for API Clay (for z < XR) Bridge Software Institute 12

1 P/Pu 0.5 Figure 1.7 Normalized P-y curve for API Clay (for z XR) Notation used for axis labels is defined as P = actual lateral resistance in stress units Pu = ultimate lateral bearing capacity in stress units y 0 0 5 10 15 = actual lateral deflection yc = defined as; y 2.5 D c c y/y c Static loading Cyclic loading where εc = strain occurring at one-half the maximum stress on laboratory undrained compression tests of undisturbed soil samples Bridge Software Institute 13

The ultimate lateral bearing capacity for sand at a given depth is used as the smaller value between pus (ultimate lateral resistance at shallow depths) to pud (ultimate lateral resistance at greater depths), which are determined by the following equations us p C z C D p 1 2 ' o p ud C D p 3 ' o where z p'o D C1,C2, & C3 = depth below ground surface = effective overburden pressure in stress units = diameter of the pile = coefficients determined from Figure 1.8, which is a function of φ φ = effective internal friction angle The variation of coefficients of C1, C2, and C3 with φ is given in the graph below. Figure 1.8 Variation of coefficients C1, C2 and C3 with φ (source: API RP2A LRFD 1993) Bridge Software Institute 14

Using the ultimate lateral resistance, the lateral soil resistance deflection (P-y) relationship for sand is approximated as kz P A pu tanh y Apu where pu = ultimate bearing capacity, which is defined as smaller value of pus and pud k = subgrade modulus, force per volume units, is determined from Figure G.8-2, which is a function of φ φ = effective internal friction angle z = depth below ground surface Y = lateral deflection A = factor to account for cyclic or static loading conditions, which is defined as For cyclic loading, A 0.9 z For static loading, A 3.0 0.8 0.9 D Bridge Software Institute 15

2. Validation of FB-MultiPier API Soil Models 2-1. Lateral soil response curves (P-y) of API sand soil model Lateral soil response curves (P-y) obtained from FB-MultiPier V4.16 and GROUP V7.0.24 are plotted for the API sand soil type and presented at three depths shown in Figure 2.1. Pile dimensions and soil properties are given in Section 3. P (KN/m) API Sand P-y Curves Depth = 2.205 m GROUP 600.0 500.0 400.0 300.0 200.0 100.0 0.0 0.00 0.02 0.04 0.06 0.08 0.10 0.12 y (m) FB-Multipier (a) P-y plots at a depth of 2.2 m API Sand P-y Curves Depth = 25.35 m GROUP FB-Multipier P (KN/m) 12000.0 10000.0 8000.0 6000.0 4000.0 2000.0 0.0 0.00 0.05 0.10 0.15 y (m) (b) P-y plots at a depth of 25.4 m API Sand P-y Curves Depth = 37.48 m GROUP FB-Multipier 25000.0 P (KN/m) 20000.0 15000.0 10000.0 5000.0 0.0 0.00 0.05 0.10 0.15 y (m) (c) P-y plots at a depth of 37.5 m Figure 2.1 Soil reaction for a laterally loaded pile Bridge Software Institute 16

2-2. Axial soil response curves (T-z, Q-z) of API sand soil model Axial soil response curves (T-z, Q-z) obtained from FB-MultiPier V4.16 and APILE Plus V4.0 are plotted for API sand soils and presented at three depths shown in Figure 2.2. Pile dimensions and soil properties are given in Section 3. It is assumed that the friction angle between the pipe wall and soil,, is determined in FB-MultiPier by subtracting 5 from the angle of internal friction in sand, = - 5, which is based on feedback from engineers at Exxon. For comparison purposes of axial soil response, of 30.75 is used for FB-MultiPier T-z model whereas of 30 is used for the APILE Plus V4.0 T-z model. It is evident that the APILE Plus V4.0 T- z model does not use the relationship: = - 5. Although the two sets of Q-z curves are not identical, their general form is in good agreement. Because the Q-z curve of FB-MultiPier has been developed using a piecewise linear curve obtained from Section G.7.3 of API RP 2A LRFD with fourteen significant figures in a double precision, the interpretation of the piecewise linear relationship between the normalized displacement (abscissa) by pile diameter and the corresponding end bearing load may cause this minor discrepancy at the early stages of loading. However, the agreement shown in Figure 2.2 suggests a degree of validity in the relationships presented by both FB-MultiPier and APILE. (a) T-z plots at a depth of 2.2 m Figure 2.2 Soil reaction for an axially loaded pile Figure 2.2 continued on the next page. Bridge Software Institute 17

Figure 2.2 cont. (b) T-z plots at a depth of 25.4 m (c) T-z plots at a depth of 37.5 m (d) Q-z plots at the tip of pile Figure 2.2 Soil reaction for an axially loaded pile Bridge Software Institute 18

2-3. Lateral soil response curves (P-y) of API clay soil model Lateral soil response curves (P-y) obtained from FB-MultiPier V4.16 and GROUP V7.0.24 are plotted for the API clay soil type and presented at three depths shown in Figure 2.3. Pile dimensions and soil properties are given in Section 3. API Clay P-y Curves Depth = 2.205 m P (KN/m) GROUP FB-Multipier 250 200 150 100 50 0 0.00 0.50 1.00 1.50 2.00 y (m) (a) P-y plots at a depth of 2.2 m API Clay P-y Curves Depth = 25.35 m GROUP FB-Multipier 500.0 P (KN/m) 400.0 300.0 200.0 100.0 0.0 0.00 0.50 1.00 1.50 2.00 y (m) (b) P-y plots at a depth of 25.4 m API Clay P-y Curves Depth = 37.48 m GROUP FB-Multipier 500 P (KN/m) 400 300 200 100 0 0.00 0.50 1.00 1.50 2.00 y (m) (c) P-y plots at a depth of 37.5 m Figure 2.3 Soil reaction for a laterally loaded pile Bridge Software Institute 19

Since the API clay soil model described as a piecewise linear curve in API RP 2A LRFD is not available in GROUP V7.0.24, a submerged soft clay model of GROUP V7.0.24 is used for comparison purpose. The submerged soft clay model describes soil reaction for a laterally loaded pile in a hyperbolic relationship. 2-4. Axial soil response curves (T-z, Q-z) of API clay soil model Axial soil response curves (T-z, Q-z) obtained from FB-MultiPier V4.16 and APILE Plus V4.0 are plotted for the API clay soil type and presented at three depths shown in Figure 2.4. Pile dimensions and soil properties are given in Section 3. It is noted that a value of 0.9 is used for the normalized residual shear stress (t/tres) in the FB- MultiPier for validation purpose. The end condition of a pile is assumed as unplugged for Q-z validation. This takes into account the annulus of the pile tip for the determination of the ultimate end bearing capacity. T (KPa) T (KPa) 10.0 8.0 6.0 4.0 2.0 API Clay T-z Curves Depth = 2.20 m Apile FB-Multipier 0.0 0.00 0.02 0.04 0.06 0.08 25.0 20.0 15.0 10.0 5.0 Displacement, z (m) (a) T-z plots at a depth of 2.2 m API Clay T-z Curves Depth = 25.35 m APILE FB-Multipier 0.0 0.00 0.02 0.04 0.06 0.08 Displacement, z (m) (b) T-z plots at a depth of 25.4 m Figure 2.4 Soil reaction for an axially loaded pile Figure 2.4 continued on the next page. Bridge Software Institute 20

Figure 2.4 cont. API Clay T-z Curves Depth = 37.48 m T (KPa) APILE FB-Multipier 25.0 20.0 15.0 10.0 5.0 0.0 0.00 0.02 0.04 0.06 0.08 Displacement, z (m) (c) T-z plots at a depth of 37.5 m API Clay Q-z Curves Depth = 48.5 m Q (KN) 80.0 70.0 60.0 50.0 40.0 30.0 APILE FB-Multipier 20.0 10.0 0.0 0.00 0.10 0.20 0.30 0.40 Displacement, z (m) (d) Q-z plots at the tip of pile Figure 2.4 Soil reaction for an axially loaded pile Bridge Software Institute 21

3. Case Study and Comparison 3-1. Case 1: Lateral soil-pile interaction analysis using API sand soil Problem description: A single pipe pile embedded in a single layer of API sand soil is subjected to incremental lateral loads at the pile head shown in Figure 3.1. Pile dimension and properties: Outer Diameter 2.4384 m Thickness 0.0508 m Unit Weight 77 KN/m3 Length 51.55 m Yield Stress 415000 KPa Modulus 200000000 KPa P Figure 3.1 FB-MultiPier lateral soil-structure interaction analysis model Lateral loads are incrementally applied to the head, ranging from 1,000 KN to 10,000 KN. Lateral displacements at the pile head obtained from FB-MultiPier V4.16 and GROUP V7.0.24 are plotted in Figure 3.2. API Lateral Sand 12000 10000 8000 Load (kn) 6000 4000 Group 7.0 FB-Multipier 2000 0 0 0.1 0.2 0.3 0.4 Pile Head Displacement (m) Figure 3.2 Lateral load-displacement results Bridge Software Institute 22

3-2. Case 2: Axial soil-pile interaction analysis using API sand soil Problem description: A single pipe pile embedded in a single layer of API sand soil is subjected to incremental axial loads at the pile head shown in Figure 3.3. Pile dimensions and properties are the same as used in Case 1. P Figure 3.3 FB-MultiPier axial soil-structure interaction analysis model Axial loads are incrementally applied to the pile head, ranging from 2,500 KN to 42,500 KN. Vertical displacements at the pile head obtained from FB-MultiPier V4.16 and APILE Plus V4.0 are plotted in Figure 3.4. The agreement in results is reasonable, but differences are observed and shown in Figure 3.4a. Particularly, predicted axial behavior of API sand soil can vary noticeably with respect to the values of the dimensionless coefficient of lateral earth pressure (i.e., ratio of horizontal to vertical normal effective stress). In order to demonstrate, two different values of the coefficient, i.e., 0.8 and 1.0, are used in both FB-MultiPier and APILE simulations. Results obtained from FB- MultiPier simulations show that full displacement piles (plugged) or open-ended piles (unplugged) can be modeled using a value of either 1.0 or 0.8, which is according to Equation G.4-5 of API RP 2A-LRFD 93, whereas no significant differences are observed in simulations using APILE Plus V4.0 (Figure 3.4b). Bridge Software Institute 23

API Axial Sand 40000 35000 30000 Load (kn) 25000 20000 15000 APILE 1.0 FB-Multipier K = 1.0 FB-Multipier K = 0.8 10000 5000 0 0 0.02 0.04 0.06 0.08 0.1 Pile Head Displacement (m) (a) Load - displacement curve comparison 40000 API Axial Sand APILE Prediction 35000 30000 Load (kn) 25000 20000 15000 K= 1.0 K = 0.8 10000 5000 0 0 0.02 0.04 0.06 0.08 0.1 Pile Head Displacement (m) (b) Results from APILE simulations Figure 3.4 Axial load-displacement results Bridge Software Institute 24

3-3. Case 3: Lateral soil-pile interaction analysis using API clay soil Problem description: A single pipe pile embedded in a single layer of API clay soil is subjected to incremental lateral loads at the pile head shown in Figure 3.5. Pile dimension and properties are same as used in Case 1. P Figure 3.5 FB-MultiPier lateral soil-structure interaction analysis model Lateral loads are incrementally applied to the head, ranging from 250 KN to 2,500 KN. Lateral displacements at the pile head obtained from FB-MultiPier V4.16 and GROUP V7.0.24 are plotted in Figure 3.6. API Lateral Clay 3500 3000 Load (kn) 2500 2000 1500 1000 API Clay FB-Multipier Soft Clay Submerged Group 7.0 Soft Clay < WT FB-MultiPier 500 0 0 0.1 0.2 0.3 0.4 0.5 Pile Head Displacement (m) Figure 3.6 Lateral load-displacement results It is noted that the lateral behavior of the API clay soil is modeled as piecewise linear P-y curves in FB-MultiPier whereas the same type of clay model is not available in GROUP V7.0.24. Instead, a hyperbolic P-y curve model (submerged soft clay) of GROUP V7.0.24 is used (Also see Section 2-3 for the difference in soil reaction Bridge Software Institute 25

predicted by FB-MultiPier V4.16 and GROUP V7.0.24). For comparison purpose, additional simulations using an existing hyperbolic P-y curve (submerged soft clay model of FB-MultiPier) are performed. It is evident that the lateral behavior of the API clay soil is softer than that of submerged soft clay models of both FB-MultiPier and GROUP. Bridge Software Institute 26

3-4. Case 4: Axial soil-pile interaction analysis using API clay soil Problem description: A single pipe pile embedded in a single layer of API clay soil is subjected to incremental axial loads at the pile head shown in Figure 3.7. Pile dimension and properties are same as used in Case 1. P Figure 3.7 FB-MultiPier axial soil-structure interaction analysis model Axial loads are incrementally applied to the head, ranging from 600 KN to 6,000 KN. Axial displacements at the pile head obtained from FB-MultiPier V4.16 and APILE Plus V4.0 are plotted in Figure 3.8. API Axial Clay 7000 6000 Axial Load (kn) 5000 4000 3000 2000 APILE FB-Multipier 1000 0 0 0.02 0.04 0.06 0.08 0.1 Displacement, z (m) Figure 3.8 Axial load-displacement results To ensure the accuracy of the numerical solution, a parametric sensitivity study is performed using several different tolerances by trial and error. Since APILE Plus V4.0 does not provide an option for user s specified tolerance in a control parameter, the length of element in APILE models is used as a means to mimic a tolerance used for a Bridge Software Institute 27

force equilibrium check in FB-MultiPier, i.e., tolerance = 0.01 KN. However, minor differences are observed where the applied load approaches a failure load. As a result, even minor differences in the length of element will yield significantly different predicted displacements at a near-failure load (Figure 3.8). Bridge Software Institute 28

3-5. Case 5: Lateral soil-pile interaction analysis using a multi-layered API clay/sand soil model Problem description: A single pipe pile embedded in a double-layer of API clay and sand soils is subjected to a lateral load at the pile head shown in Figure 3.9. Pile dimensions and properties are the same as used in Case 1. Figure 3.9 FB-MultiPier multi-layered soil-structure interaction analysis model Case 5 represents an extreme lateral load for the soft upper clay layer. The 7500 KN load is possible because of the strong sand layer located beneath the clay. This loading develops the plastic bending moment in the pipe pile equal to 130,591 KN-m at 26.5 meters below the soil surface. The near equal lateral displacements of the cantilever pile, 1.33m for FB-MultiPier V4.16 and 1.29m for GROUP V7.0.24, demonstrates that the soil-pile lateral stiffness modeled by the two computer programs compare very well. The difference in displacement prediction is believed to be due to a different numerical solution procedure and a tolerance used in convergence checks in the programs. Additional simulations are conducted using different tolerances and results are presented in the table given below. Load (KN) Displacement (m) Model Tolerance 0.92 GROUP 0.1 m 5500 0.87 FB-MultiPier 0.01 KN 0.80 GROUP 0.0001 m 1.50 GROUP 0.1 m 7500 1.33 FB-MultiPier 0.01 KN 1.28 GROUP 0.0001 m Bridge Software Institute 29

3-6. Case 6: Combined lateral/axial soil-structure interaction analysis using a multilayered API sand soil model Problem description: A single pipe pile embedded in a double-layer of API sand soil is subjected to both lateral and axial loads at the pile head shown in Figure 3.10. Pile dimensions and properties are the same as used in Case 1. Figure 3.10 FB-MultiPier soil-structure interaction analysis model for a combined loading condition Case 6 represents a two-layer sand model and a combination of axial load (25,000KN) and lateral load (9000 KN) is applied to validate the soil pile stiffness represented in FB-MultiPier V4.16 against GROUP V7.0.24. The lateral load is selected so as to develop the first yield moment in the pipe-pile. Note that the lateral load is considerably larger than the lateral load (7500KN) which would develop the plastic moment in case 5, where the weak clay layer is present in the surface layer. Loads and corresponding displacements of the pile head are investigated and presented in the table below. Although minor discrepancies are observed, it is considered mainly due to numerical issues related to convergence (tolerance). Additional simulations using different tolerances and results are compared in three different load cases. Displacements predicted by FB-MultiPier V4.16 and by GROUP V7.0.24 demonstrate that both the soilpile lateral/axial stiffness modeled by the two computer programs is in good agreement. Lateral Load (KN) 4500 7500 9000 Vertical Load (KN) 25000 25000 25000 FB-MultiPier Lateral Displ. (m) 0.08 0.20 0.28 FB-MultiPier Vertical Displ. (m) 0.013 0.013 0.013 GROUP Lateral Displ. (m) 0.08 0.21 0.30 GROUP Vertical Displ. (m) 0.011 0.011 0.011 Bridge Software Institute 30

4. Summary and Conclusion Six soil curves that represent the behaviors of the API sand and clay soils have been implemented and validated. Using FB-MultiPier, soil-pile interaction has been investigated in connection with the development of the API sand and clay soil models. The numerical model presented here also accounts for three important factors involved in load transfer of piles: (1) pile end condition, (2) the concept of equivalent depth, and (3) effects of effective vertical stress on soil behavior of sub-layer soils in determination of ultimate resistance of sub-layer soil in vicinity of the interface. The new soil models presented here offer an improved tool for studying multidimensional soil-pile interaction subjected to various loading conditions. In addition, it serves as a basis for developing improved FB-MultiPier soil models. For example, one factor influencing soil behavior under dynamic loading, i.e., soil gap forming, is currently being investigated for the application of the FB-MultiPier soil models to seismic analysis. In the future, an enhanced soil model incorporating such effects may be implemented once additional experimental data become available. Bridge Software Institute 31

5. References Hoit, M., McVay, M., Hays, C., (2007) FB-MultiPier Version 4.11 User s manual, Bridge Software Institute, University of Florida, Gainesville Florida. American Petroleum Institute, (1993) Recommended Practice for Planning, Designing and Constructing Fixed Offshore Plattforms Load and Resistance Factor Design, API Recommended practice 2A-LRFD, 1 st edition. Bowles, J.E., (1996) Foundation Analysis and Design, 5 th edition, McGraw-Hill Company, New York, NY. Reese, L.C., Wang, S.T., Vasquez, L., (2006) Analysis of a Group of piles Subjected to Axial and Lateral Loading, GROUP Version 7.0 Technical manual, Ensoft, Inc., Austin Texas. Reese, L.C., Wang, S.T., Vasquez, L., (2004) A Program for the Analysis of the Axial Capacity of Driven Piles, APILE Plus Version 4.0 Technical manual, Ensoft, Inc., Austin Texas. Bridge Software Institute 32

Appendix A Showcase model: soil-pile interaction analysis using FB-MultiPier Problem description: A single pile embedded in soil that consists of 16 layers of both sand and clay soil is subjected to multiple load at the pile head shown in the Figure below. Pile dimensions and properties are the same as used in the case study in Section 3. The soil properties provided by a consultant are used to develop the model. Comparison is made between results obtained from the new API soil models and existing soil models Sand (Reese) and Soft Clay (below water table) of FB- MultiPier. In general, the lateral soil resistance of the API sand soil model seems greater than that of an existing sand model Sand (Reese) of FB-MultiPier. The pile head displacements in lateral and vertical directions are given in the table. Analysis results from using the new API soil models Analysis results from using the existing soil models Lateral Displacement (m) Y-rotation (radian) Vertical Displacement (m) 0.283-0.0231 0.012 0.329-0.0250 0.012 Bridge Software Institute 33

Appendix B Application of API Soil Models in FB-MultiPier Select the Soil Page. Choose a soil Type. Bridge Software Institute 34

Select each of the four soil models: Lateral, Axial, Torsional, and Tip. Enter properties for each soil Model. To do this, click on the soil model, then click the edit button. (The Soil Table can also be used to enter all soil properties) Bridge Software Institute 35

To display a soil curve, select the desired soil model. Then click a node on pile in the Soil Edit Window. Bridge Software Institute 36

This will display a curve at the elevation of the selected node. Bridge Software Institute 37