Openseespl: 3D Lateral Pile-Ground Interaction User Manual (Beta 1.0)
Openseespl: 3D Lateral Pile-Ground Interaction User Manual (Beta 1.0)
Openseespl: 3D Lateral Pile-Ground Interaction User Manual (Beta 1.0)
http://cyclic.ucsd.edu/openseespl/
December 2011
For conversion between SI and English Units, please check:
http://www.unit-conversion.info/
(1 kPa = 0.1450378911491 psi)
(1 psi = 6.89475 kPa)
i
Table of Contents
ii
7.1.2 Input Motion 62
7.1.3 Model Inclination 64
7.2 TIME HISTORY OUTPUT 65
7.2.1 Soil Response Time Histories 65
8. PILE GROUP 67
8.1 PILE GROUP PARAMETERS 67
8.2 PILE GROUP MESHING 67
8.3 OUTPUT FOR A PILE GROUP MODEL 69
APPENDIX A HOW TO DEFINE THE SOIL FINITE ELEMENT MESH 73
APPENDIX B OWN WEIGHT APPLICATION WITH DRY AND SATURATED SOIL
CASES 77
APPENDIX C BENCHMARK LINEAR FINITE ELEMENT ANALYSIS OF LATERALLY
LOADED SINGLE PILE USING OPENSEES & COMPARISON WITH ANALYTICAL
SOLUTION 83
APPENDIX D FINITE ELEMENT ANALYSIS OF ARKANSAS TEST SERIES PILE #2
USING OPENSEES (WITH LPILE COMPARISON) 96
APPENDIX E FINITE ELEMENT ANALYSIS OF STANDARD CALTRANS 16" CIDH
PILE USING OPENSEES FOR GENERAL COMPARISON WITH LPILE (WITH DEFAULT
P-Y MULTIPLIER = 1.0) 107
APPENDIX F FINITE ELEMENT ANALYSIS OF CALTRANS 42" CIDH PILE USING
OPENSEES FOR GENERAL COMPARISON WITH LPILE (WITH DEFAULT P-Y
MULTIPLIER = 1.0) 115
APPENDIX G FINITE ELEMENT ANALYSIS OF STANDARD CALTRANS 16" CIDH
PILE SUBJECTED TO AXIAL LOAD 128
APPENDIX H MOMENT-CURVATURE ANALYSIS OF CIRCULAR NONLINEAR RC
BEAM (FIBER SECTION) 133
REFERENCES 145
OPENSEESPL-RELATED REFERENCES 146
iii
Table of Figures
Figure 2.1: OpenSeesPL main window. ......................................................................................... 4
Figure 2.2: OpenSeesPL’s menu bar and submenu bars: a) menu bar; b) menu File; c) menu
Execute; d) menu Display; and e) menu Help.......................................................................... 5
Figure 2.3: OpenSeesPL copyright message. ................................................................................. 6
Figure 2.4: Buttons available in the Finite Element Mesh window. ............................................... 7
Figure 3.1: Definition of pile model. .............................................................................................. 8
Figure 3.2: Definition of linear pile properties. .............................................................................. 9
Figure 3.3: Definition of nonlinear pile properties (Aggregator Section). ................................... 10
Figure 3.4: Definition of nonlinear pile properties (Fiber Section). ............................................. 12
Figure 3.5: Moment-curvature response for the pile (with default steel and concrete parameters)
................................................................................................................................................ 14
Figure 3.6: Material Parameters of the Concrete01 material (Mazzoni et al. 2006) . .................. 15
Figure 3.7: Typical hysteretic stress-strain relation of the Concrete01 material (Mazzoni et al.
2006) ...................................................................................................................................... 15
Figure 3.8: Material Parameters of the Steel01 material (Mazzoni et al. 2006) . ......................... 16
Figure 3.9: Typical hysteretic behavior of model with Isotropic hardening of the Steel01 material
(Mazzoni et al. 2006) . ........................................................................................................... 16
Figure 3.10: Schematic of fiber section definition for a circular cross section (Mazzoni et al.
2006) ...................................................................................................................................... 17
Figure 3.11: Schematic of patch definition for a circular cross section (Mazzoni et al. 2006) .... 17
Figure 3.12: Schematic of layer definition for a circular cross section (Mazzoni et al. 2006) ..... 18
Figure 4.1: Soil strata definition. .................................................................................................. 19
Figure 4.2: Analysis options. ........................................................................................................ 21
Figure 4.3: Rayleigh damping coefficients. .................................................................................. 22
Figure 4.4: OpenSees parameters. ................................................................................................ 23
Figure 4.5: Multi-yield surfaces in principal stress space and deviatoric plane (Prevost 1985;
Parra 1996; Yang 2000) ......................................................................................................... 24
Figure 4.6: Shear-effective confinement and shear stress-strain response (Yang and Elgamal
2002; Yang et al. 2003). ......................................................................................................... 25
Figure 4.7: Von Mises multi-surface kinematic plasticity model (Yang 2000; Yang et al. 2003).
................................................................................................................................................ 26
Figure 4.8: Soil materials in OpenSeesPL. ................................................................................... 27
Figure 4.9: Soil backbone curve and yield surfaces. .................................................................... 30
Figure 4.10: Backbone curves for Medium Sand. ........................................................................ 31
Figure 4.11: U-Sand1A. ................................................................................................................ 32
Figure 4.12: Initial yield domain at low levels of effective confinement (Yang et al. 2003). ...... 36
Figure 4.13: Schematic of constitutive model response showing (a) octahedral stress - effective
confinement p response, (b) octahedral stress - octahedral strain response, and (c)
configuration of yield domain (Yang et al. 2003). ................................................................. 36
Figure 4.14: U-Sand1B. ................................................................................................................ 37
Figure 4.15: U-Clay1. ................................................................................................................... 40
Figure 4.16: U-Clay2. ................................................................................................................... 41
Figure 4.17: U-Sand2A. ................................................................................................................ 42
iv
Figure 4.18: U-Sand2B. ................................................................................................................ 43
Figure 4.19: Pile zone material. .................................................................................................... 44
Figure 4.20: Pile-soil interfacing layer material. .......................................................................... 45
Figure 4.21: Outermost zone material. ......................................................................................... 45
Figure 5.1: Definition of meshing parameters. ............................................................................. 48
Figure 6.1: Pushover analysis. ...................................................................................................... 51
Figure 6.2: Pushover load pattern. ................................................................................................ 52
Figure 6.3: User-defined pushover load pattern (U-Push) ............................................................ 53
Figure 6.4: Analysis running progress window. ........................................................................... 54
Figure 6.5: Zoom in. ..................................................................................................................... 54
Figure 6.6: Zoom out. ................................................................................................................... 55
Figure 6.7: Response time histories and profiles for pile. ............................................................ 56
Figure 6.8: Response relationships for pile................................................................................... 57
Figure 6.9: Deformed mesh and contour fill ................................................................................. 58
Figure 6.10: 2D plane (Y = 0) view of the longitudinal displacement contour in the deformed
mesh window.......................................................................................................................... 59
Figure 6.11: Output for an Eigenvalue analysis............................................................................ 60
Figure 7.1: Newmark Time Integration ........................................................................................ 61
Figure 7.2: Definition of 3D base excitation and boundary conditions. ....................................... 63
Figure 7.3: User-defined input motion (U-Shake). ....................................................................... 64
Figure 7.4: Response time histories window. ............................................................................... 66
Figure 8.1: Pile group definition. .................................................................................................. 67
Figure 8.2: Pile group horizontal meshing. ................................................................................... 68
Figure 8.3: Sample mesh of a 3 by 3 pile group model (half mesh configuration). ..................... 69
Figure 8.4: Pile response profiles for a pile group model. ............................................................ 70
Figure 8.5: Pile response time histories for a pile group model. .................................................. 70
Figure 8.6: Pile response relationships for a pile group model..................................................... 71
Figure 8.7: Pile response relationships at the pile cap for a pile group model. ............................ 71
Figure 8.8: Deformed mesh of a pile group model. ...................................................................... 72
Figure A.1: Finite element mesh created with default values. ...................................................... 73
Figure A.2: Mesh refinement example 1: a) Change “Num of Slices” to 32; b) the resulting mesh
................................................................................................................................................ 74
Figure A.3: Mesh refinement example 2: a) Change “Number of Mesh Layers” in the vertical
direction; b) the resulting mesh .............................................................................................. 75
Figure A.4: Mesh refinement example 3: a) Change meshing controlling parameters in the
horizontal direction; b) the resulting mesh ............................................................................. 76
Figure C.1: Finite element mesh employed in this study.............................................................. 85
Figure C.2: Comparison of pile deflection profiles (νs=.25, l/a=50). ........................................... 86
Figure C.3: Comparison of pile bending moment profiles (νs=.25, l/a=50). ................................ 87
Figure C.4: Sample pile deflection (h/a=.1, l/a=50) under an applied pure pile-head horizontal
load (Abedzadeh and Pak, 2004)............................................................................................ 89
Figure C.5: Sample pile bending moment (h/a=.1, l/a=50) under an applied pure pile-head
horizontal load (Abedzadeh and Pak, 2004). ......................................................................... 90
Figure C.6: Comparison of the load-deflection curves for the linear and nonlinear runs. ........... 91
Figure C.7: Comparison of the pile deflection profiles for the linear and nonlinear runs. ........... 92
v
Figure C.8: Comparison of the pile bending moment profiles for the linear and nonlinear runs. 93
Figure C.9: Stress ratio contour fill of the nonlinear run at different load levels (red color shows
yielded soil elements). ............................................................................................................ 95
Figure D.1: Finite element mesh employed in this study. ............................................................ 98
Figure D.2: Comparison of the load-deflection curves for the linear and nonlinear runs. ........... 99
Figure D.3: Comparison of the pile deflection profiles for the linear and nonlinear runs. ......... 100
Figure D.4: Comparison of the pile bending moment profiles for the linear and nonlinear runs.
.............................................................................................................................................. 101
Figure D.5: Stress ratio contour fill of the nonlinear run at different load levels (red color shows
yielded soil elements). .......................................................................................................... 103
Figure D.6: Comparison of the pile deflection profiles for the linear and nonlinear runs. ......... 104
Figure D.7: Comparison of the pile bending moment profiles for the linear and nonlinear runs.
.............................................................................................................................................. 105
Figure E.1: Finite element mesh employed in this study. ........................................................... 109
Figure E.2: Comparison of pile deflection profiles for load case 1. ........................................... 110
Figure E.3: Comparison of pile rotation profiles for load case 1................................................ 111
Figure E.4: Comparison of bending moment profiles for load case 1. ....................................... 111
Figure E.5: Comparison of shear force profiles for load case 1. ................................................ 112
Figure E.6: Comparison of pile deflection profiles for load case 2. ........................................... 112
Figure E.7: Comparison of pile rotation profiles for load case 2................................................ 113
Figure E.8: Comparison of bending moment profiles for load case 2. ....................................... 113
Figure E.9: Comparison of shear force profiles for load case 2. ................................................ 114
Figure E.10: Stress ratio contour fill for load case 1 (red color shows yielded soil elements)... 114
Figure E.11: Stress ratio contour fill for load case 2 (red color shows yielded soil elements)... 114
Figure F.1: Finite element mesh employed in this study. ........................................................... 117
Figure F.2: Comparison of pile deflection profiles for the fixed-head condition. ...................... 118
Figure F.3: Comparison of pile rotation profiles for the fixed-head condition. ......................... 118
Figure F.4: Comparison of bending moment profiles for the fixed-head condition. .................. 119
Figure F.5: Comparison of shear force profiles for the fixed-head condition. ........................... 119
Figure F.6: Comparison of pile deflection profiles for the free-head condition. ........................ 120
Figure F.7: Comparison of pile rotation profiles for the free-head condition. ........................... 120
Figure F.8: Comparison of bending moment profiles for the free-head condition. .................... 121
Figure F.9: Comparison of shear force profiles for the free-head condition. ............................. 121
Figure F.10: Stress ratio contour fill of the nonlinear run for the fixed-head condition (red color
shows yielded soil elements). ............................................................................................... 122
Figure F.11: Stress ratio contour fill of the nonlinear run for the free-head condition (red color
shows yielded soil elements). ............................................................................................... 123
Figure F.12: Comparison of pile deflection profiles for the fixed-head condition. .................... 124
Figure F.13: Comparison of pile rotation profiles for the fixed-head condition. ....................... 124
Figure F.14: Comparison of bending moment profiles for the fixed-head condition. ................ 125
Figure F.15: Comparison of shear force profiles for the fixed-head condition. ......................... 125
Figure F.16: Comparison of pile deflection profiles for the free-head condition. ...................... 126
Figure F.17: Comparison of pile rotation profiles for the free-head condition. ......................... 126
Figure F.18: Comparison of bending moment profiles for the free-head condition. .................. 127
vi
Figure F.19: Comparison of shear force profiles for the free-head condition. ........................... 127
Figure G.1: Finite element mesh employed in this study. .......................................................... 130
Figure G.2: Pile profile response at the axial load of 243 kips. .................................................. 131
Figure G.3: Close-up of final deformed mesh (factor of 120). ................................................... 131
Figure G.4: Stress ratio contour fill for the nonlinear analysis (red color shows yielded soil
elements). ............................................................................................................................. 132
Figure H.1: Schematic of the fiber section definition for the circular pile cross section. .......... 135
Figure H.2: Material properties for the Fiber section. ................................................................ 135
Figure H.3: Finite element mesh employed in this study. .......................................................... 136
Figure H.4: Comparison of the moment-curvature curves calculated by using OpenSeesPL and
OpenSees Example. .............................................................................................................. 136
Figure H.5: Displacement response profiles histories of the pile. .............................................. 138
Figure H.6: Lateral (longitudinal) shear versus displacement at the pile head. .......................... 138
Figure H.7: Moment-curvature relation at the maximum moment location (ground surface) in
OpenSeesPL. ........................................................................................................................ 139
vii
1. Introduction
OpenSeesPL is a graphical user interface (GUI) for three dimensional (3D) ground and ground-
structure response. The OpenSees Finite Element (FE) Computational Analysis framework
(http://opensees.berkeley.edu) is employed to conduct all analyses. The OpenSeesPL graphical
interface (pre- and post-processor) is focused on facilitating a wide class of 3D studies (with
additional capabilities yet under development). In the current version, OpenSeesPL may be
employed to study a number of geometries and configurations of interest including:
- Linear and nonlinear (incremental plasticity based) 3D ground seismic response with
capabilities for 3D excitation, and multi-layered soil strata. Multi-yield surface cohesionless
(Drucker-Prager cone model), and (Mises or J2) soil models are available. The coupled solid-
fluid analysis option allows for conducting liquefaction studies.
- Inclusion of a pile or shaft in the above 3D ground mesh (circular or square pile in a soil
island). The pile can extend above ground and can support a bridge deck, or a point mass at the
pile top. The bridge deck can be specified to only translate laterally, or to undergo both lateral
translation and rotation. In addition to the seismic excitation option, the pile system may be
subjected to monotonic or cyclic lateral push-over loading (in prescribed displacement or
prescribed force modes). Soil within the zone occupied by the pile (as specified by pile diameter
for instance) can be specified independently, allowing for a variety of useful modeling scenarios.
- Slopes and pile systems embedded in sloping ground are also currently being simulated.
1.1 Overview
OpenSeesPL is a FE user-interface for 3D lateral pile-ground interaction response. This interface
allows conducting pushover pile analyses as well as seismic (earthquake) simulations. The FE
analysis engine for this interface is the Pacific Earthquake Engineering Research (PEER) Center
OpenSees Framework, developed under the leadership of Professor Gregory Fenves of UC
Berkeley. For more information, please visit http://opensees.berkeley.edu/.
OpenSeesPL allows simulations for any size of pile and pile diameter. The pile cross section can
be circular or square. Linear and nonlinear material properties options are available for pile
definition.
1
OpenSeesPL allows for definition of multiple soil strata. Nonlinearity of soil materials is
simulated by incremental plasticity models to allow for modeling permanent deformation and for
generation of hysteretic damping. In addition, OpenSeesPL allows including user-defined soil
materials.
OpenSeesPL allows for convenient pre-processing and graphical visualization of the analysis
results including the deformed mesh, ground response time histories and pile responses.
OpenSeesPL makes it possible for geotechnical and structural engineers/researchers to quickly
build a model, run FE analysis and evaluate the performance of the pile-ground system.
NOTE: Seismically-induced deformations are complex mechanisms. Much expertise and sound
engineering judgment are necessary in interpreting the OpenSeesPL computational results.
Internet Explorer 3.0 or above (or compatible Browser) with Java Applet enabled is needed to
view the graphic results. For best results, your system’s video should be set to 1024 by 768 or
higher.
1.3 Installation
After downloading the OpenSeesPL installation file (OpenSeesPL_Setup.exe), double-click on
the icon and the installation procedure will start. Once installed, the default case in OpenSeesPL
is a good way to go through the steps involved in conducting an OpenSeesPL analysis. The
interface will allow the user to prepare and save an input file, to run the analysis, and to display
the response.
Note: Tcl/tk 8.5 must be installed in order to run OpenSeesPL. Please restart the computer after
the installation of Tcl/tk 8.5 for the change to take effect.
2
1.4 Acknowledgments
OpenSeesPL is based on research underway since the early 1990s, and a partial list of related
publications is included in the Appendix section. The OpenSeesPL graphical interface is written
in Microsoft Visual C++ 2005 with the Microsoft Foundation Class (MFC) libraries. The Java
Applet package used to display graphical results in OpenSeesPL is obtained from the website
http://ptolemy.eecs.berkeley.edu/. GIF images are generated with GNUPLOT for MS-Windows
32 bit Version 3.7, available at http://www.gnuplot.org/.
3
2. Getting Started
2.1 Start-Up
On Windows start OpenSeesPL from the Start button, or from an icon on your desktop. To Start
OpenSeesPL from the Start button:
2.2 Interface
There are 3 main regions in the OpenSeesPL window – menu bar, the model input window, and
the finite element mesh window.
4
2.2.1 Menu Bar
The menu bar, shown in Figure 2.2, offers rapid access to most of OpenSeesPL’s main features.
a)
c)
b)
e)
d)
Figure 2.2: OpenSeesPL’s menu bar and submenu bars: a) menu bar; b) menu File; c) menu
Execute; d) menu Display; and e) menu Help.
File: Controls reading, writing and printing of model definition parameters, and exiting
OpenSeesPL.
Execute: Controls running analyses.
Display: Controls displaying of the analysis results.
Help: Visits OpenSeesPL website and display the copyright info (Figure 2.3).
5
Figure 2.3: OpenSeesPL copyright message.
The model input window controls definitions of the model and analysis options, which are
organized into four regions (Figure 2.1):
Model Definition: Controls definitions of pile and soil strata including material properties.
Meshing parameters are also defined.
Analysis Types: Controls analysis options: pushover analysis, Eigenvalue analysis or base
shaking simulation.
Boundary Conditions: Controls boundary conditions.
Model Inclination: Controls the inclination angles for the ground surface and the whole
model.
The finite element mesh window (Figure 2.1) displays the mesh generated. Once the mesh
window is focused, the mesh can be rotated by dragging the mouse, moved in 4 directions by
pressing keys of LEFT ARROW, RIGHT ARROW, UP ARROW or DOWN ARROW
6
respectively. The view can be zoomed in (by pressing key ‘F9’), out (by pressing key ‘F10’) or
frame (by pressing key ‘F11’).
To display a 2D view, press key ‘F2’ (for Plane XY, where X is the longitudinal directon, Y the
transverse direction), ‘F3’ (for Plane YZ, where Z is the vertical direction) or ‘F4’ (for Plane
XZ). An isometric view of the mesh can be achieved by pressing key ‘F5’.
Alternatively, users can press the corresponding button shown in Figure 2.4.
7
3. Pile Model
To define pile geometry, click Pile Parameters in the Model Input window. The pile geometry
is defined by the following parameters (Figure 3.1):
Pile Diameter/Side Length (D) The diameter (if a circular pile is chosen), or the side length (if
a square pile is chosen) of the pile cross section. The value entered must be greater than zero.
Total Pile Length The total length of the pile. The value entered must be greater than zero.
Pile Length above Surface The height of the pile above the ground surface. The value entered
must be greater than zero.
8
Axial Load The axial load applied at the pile head (positive as compression).
If checkbox Pile Group is enabled (note that the pile group option might not be available in the
version you have), users can activate pile group by checking Pile Group. Please see Chapter 8
for the detailed information.
The material properties of the pile for the linear beam element (elasticBeamColumn) are
defined by the following parameters (Figure 3.2):
Moment of Inertia (I) The Moment of Inertia of the pile. This can be specified directly or
calculated based on the pile diameter.
OpenSees uses the Section command to define the nonlinear beam-column element (a section
defines the stress resultant force-deformation response at a cross section of a beam-column
element). Two types of sections are available in OpenSeesPL for the nonlinear beam element
(nonlinearBeamColumn): Aggregator Section or Fiber Section. Detail information can be
found in the OpenSees User Manual (Mazzoni et al. 2006).
9
The Aggregator Section is defined by the following parameters in OpenSeesPL (Figure 3.3):
Flexural Rigidity My & Mz The Flexural Rigidity of the pile which is equal to the product of
Young’s Modulus (E) and the Moment of Inertia (I). My corresponds the moment-curvature
about section local y-axis and Mz corresponds the moment-curvature about section local z-axis.
Shear Rigidity Vy & Vz The Shear Rigidity of the pile which is equal to the product of the
Shear Modulus (G) and the area of the pile cross section (A). Vy corresponds the shear force-
deformation along section local y-axis and Vz corresponds the shear force-deformation along
section local z-axis.
Torsional Rigidity T The Torsional Rigidity of the pile which is equal to the product of the
Shear Modulus (G) and J.
Axial Rigidity P The Axial Rigidity of the pile which is equal to the product of Young’s
Modulus (E) and the area of the pile cross section (A).
10
The dialog of defining Fiber Section is shown in Figure 3.4 (the Fiber Section is only available to
circular pile in this version of OpenSeesPL). Two materials are available: Concrete01 and
Steel01 in this version of OpenSeesPL. Concrete01 (Figure 3.6) is defined by the following
parameters (for Core and Cover, see Figure 3.10):
Concrete Compressive Strength The concrete compressive strength at 28 days ($fpc in Figure
3.6).
Concrete Strain at Maximum Strength The concrete strain at maximum strength ($epsc0 in
Figure 3.6).
Concrete Crushing Strength The concrete crushing strength ($fpcu in Figure 3.6).
Concrete Strain at Crushing Strength The concrete strain at crushing strength ($epsU in
Figure 3.6).
Note that the compressive concrete parameters should be input as negative values. Typical
hysteretic stress-strain relation of the Concrete01 material is shown in Figure 3.7).
Steel01 is defined by the following parameters (Figure 3.8 and Figure 3.9):
Strain-hardening Ratio The strain-hardening ratio (ratio between post-yield tangent and initial
elastic tangent)
Patch (Figure 3.10) is defined by the following parameters (for both Core and Cover):
Number of Subdivisions (fibers) in the Radial Direction The number of subdivisions (fibers)
in the radial direction of the pile circular cross section ($numSubdivRad in Figure 3.11).
Internal Radius The internal radius of the patch ($intRad in Figure 3.11).
External Radius The external radius of the patch ($extRad in Figure 3.11).
The values of $yCenter and $zCenter (y & z-coordinates of the center of the circle) as shown in
Figure 3.11 are zeros. And the $startAng (starting angle) and $endAng (ending angle) are set to 0
and 360 degrees respectively in OpenSeesPL since only a full mesh is available for fiber section
nonlinear beam element).
11
Layer is defined by the following parameters (Figure 3.12):
Number of Reinforcing Bars along Layer The number of reinforcing bars along layer
($numBars in Figure 3.12).
Radius of Reinforcing Layer The radius of reinforcing layer ($radius in Figure 3.12) .
The values of $yCenter and $zCenter (y & z-coordinates of the center of the circle) as shown in
Figure 3.12 are zeros. And the $startAng (starting angle) and $endAng (ending angle) are set to 0
and 360 degrees respectively in OpenSeesPL since only a full mesh is available for fiber section
nonlinear beam element).
12
The moment-curvature response for the pile is shown in Figure 3.5 (for default steel and concrete
parameters).
a)
13
b)
Figure 3.5: Moment-curvature response for the pile (with default steel and concrete parameters)
14
Figure 3.6: Material Parameters of the Concrete01 material (Mazzoni et al. 2006) .
Figure 3.7: Typical hysteretic stress-strain relation of the Concrete01 material (Mazzoni et al.
2006) .
15
Figure 3.8: Material Parameters of the Steel01 material (Mazzoni et al. 2006) .
Figure 3.9: Typical hysteretic behavior of model with Isotropic hardening of the Steel01 material
(Mazzoni et al. 2006) .
16
Figure 3.10: Schematic of fiber section definition for a circular cross section (Mazzoni et al.
2006) .
Figure 3.11: Schematic of patch definition for a circular cross section (Mazzoni et al. 2006)
17
Figure 3.12: Schematic of layer definition for a circular cross section (Mazzoni et al. 2006)
18
4. Soil Parameters
To define soil strata, click Soil Parameters in the Model Input window (Figure 4.1).
Thickness The thickness for a soil layer. Definitions following a zero height will be ignored. In
other words, the total number of soil layers in use will be equal to the number of the last soil
layer that contain no zero values, e.g., if you need 5 strata, enter nonzero heights for Stratum #1
through Stratum #5.
To perform a liquefaction analysis, check the checkbox Saturated Soil Analysis (Figure 4.1)
and specify the water table depth:
19
Water Table Depth The Water Table Depth refers to the depth below ground surface.(e .g., 0.0
corresponds to a fully saturated soil profile, 1.0 is 1m below ground surface). Dry sites should
specify water table depth to be equal to the entire model depth.
First, some important master control options are defined by clicking Analysis Options as shown
in Figure 2.1. This will display the interface shown in Figure 4.2. Here you can:
1. Select to keep the soil properties as defined by their linear properties, or opt to conduct
nonlinear soil computations (note that the default is Linear),
2. Select among a number of available Brick elements in OpenSees,
3. Apply own weight of the soil using a global lateral stress coefficient, and a single value of
Young’s modulus that is user defined (this will reduce initial shear stresses in the mesh due to
own weight application, but generally will have minimal impact on the subsequent earthquake
computations anyway),
4. Apply own weight of the soil using a global permeability (horizontal & vertical), e.g., one can
specify a large permeability value for the application of own weight in a saturated soil analysis,
5. by clicking Rayleigh Damping (Figure 4.3) you can change the viscous damping
characteristics of the model, and
6. by clicking OpenSees Parameters (Figure 4.4) you can OpenSees analysis parameters
(advanced feature, please exercise with care).
C = Am M + Ak K
where M is the mass matrix, C is the viscous damping matrix, K is the initial stiffness matrix.
Am and Ak are two user-specified constants.
Am
Ak f
4 f
where f is frequency.
20
Figure 4.2: Analysis options.
The user can define damping coefficients (Figure 4.3) by specifying two frequencies, f1 and f2
(must be between 0.1 and 50 Hz), and two damping ratios, 1 and 2 (suggested values are
between 0.2% and 20%).
The Rayleigh damping parameters Am and Ak are obtained by solving the follow equations
simultaneously:
A
1 m Ak f1
4 f1
A
2 m Ak f 2
4 f 2
The user can also directly define Rayleigh damping coefficients Am and Ak (Figure 4.3).
21
Figure 4.3: Rayleigh damping coefficients.
22
Figure 4.4: OpenSees parameters.
In OpenSees, the soil model (Figure 4.5) for cohesionless soils is developed within the
framework of multi-yield-surface plasticity (e.g., Prevost 1985). In this model, emphasis is
placed on controlling the magnitude of cycle-by-cycle permanent shear strain accumulation
(Figure 4.6) in clean medium to dense sands (Parra 1996; Yang 2000; Yang et al. 2003).
Furthermore, appropriate loading-unloading flow rules were devised to reproduce the observed
strong dilation tendency, and resulting increase in cyclic shear stiffness and strength (the “Cyclic
Mobility” mechanism). The material types for the cohesionless soils in OpenSees are called
PressureDependMultiYield and PressureDependMultiYield02.
23
Clay material is modeled as a nonlinear hysteretic material (Parra 1996; Yang 2000; Yang et al.
2003) with a Von Mises multi-surface (Iwan 1967; Mroz 1967) kinematic plasticity model
(Figure 4.7). In this regard, focus is on reproduction of the soil hysteretic elasto-plastic shear
response (including permanent deformation). In this material, plasticity is exhibited only in the
deviatoric stress-strain response. The volumetric stress-strain response is linear-elastic and is
independent of the deviatoric response. This constitutive model simulates monotonic or cyclic
response of materials whose shear behavior is insensitive to the confinement change. Plasticity is
formulated based on the multi-surface (nested surfaces) concept, with an associative flow rule
(according to the well-known Provost approach). In the clay model, the nonlinear shear stress-
strain back-bone curve is represented by the hyperbolic relation (Kondner 1963), defined by the
two material constants, low-strain shear modulus and ultimate shear strength. The material type
for the cohesive soils in OpenSees is called PressureIndependMultiYield.
3
2
1
1 p
p0 3
2 3
3
2 2
3 Deviatoric plane
Figure 4.5: Multi-yield surfaces in principal stress space and deviatoric plane (Prevost 1985;
Parra 1996; Yang 2000)
24
Figure 4.6: Shear-effective confinement and shear stress-strain response (Yang and Elgamal
2002; Yang et al. 2003).
1
1 2 3
3
2
(a) Von Mises multi-surface.
25
Figure 4.7: Von Mises multi-surface kinematic plasticity model (Yang 2000; Yang et al. 2003).
As shown in Figure 4.1, the soil materials can be selected from an available menu of
cohesionless and cohesive soil materials (Figure 4.8). There are 18 predefined materials using
the PressureDependMultiYield soil model. Basic model parameter values for these materials
are listed in Table 4.1.
If ‘Cohesionless very loose’ is chosen, the user is allowed to define the residual shear strength
(0.2 kPa is specified by default). The cohesionless very loose soil is same as the cohesionless
loose soil except the user is allowed to specify the residual shear strength for the very loose one.
As shown in Figure 4.1, parabolic variation of soil modulus with depth is used if P is selected.
Linear variation of soil modulus with depth is used if L is selected. And the constant soil
modulus with depth is used if C is selected;
26
Figure 4.8: Soil materials in OpenSeesPL.
27
Table 4.1: Predefined soil materials in OpenSeesPL
Reference
Reference
bulk Friction Mass
shear modulus Permeability
Cohesionless Soil modulus Br angle density4
Gr (kPa, at coeff.3 (m/s)
p r =80kPa) 1
(kPa, at (degrees) 2 (ton/m3)
p r =80kPa)
Very loose, silt
5.5E+04 1.5E+05 29 1.0E-07 1.7
permeability
Very loose, sand
5.5E+04 1.5E+05 29 6.6E-05 1.7
permeability
Very loose, gravel
5.5E+04 1.5E+05 29 1.0E-02 1.7
permeability
Loose, silt permeability 5.5E+04 1.5E+05 29 1.0E-07 1.7
Loose, sand permeability 5.5E+04 1.5E+05 29 6.6E-05 1.7
Loose, gravel
5.5E+04 1.5E+05 29 1.0E-02 1.7
permeability
Medium, silt permeability 7.5E+04 2.0E+05 33 1.0E-07 1.9
Medium, sand
7.5E+04 2.0E+05 33 6.6E-05 1.9
permeability
Medium, gravel
7.5E+04 2.0E+05 33 1.0E-02 1.9
permeability
Medium-dense, silt
1.0E+05 3.0E+05 37 1.0E-07 2.0
permeability
Medium-dense, sand
1.0E+05 3.0E+05 37 6.6E-05 2.0
permeability
Medium-dense, gravel
1.0E+05 3.0E+05 37 1.0E-02 2.0
permeability
Dense, silt permeability 1.3E+05 3.9E+05 40 1.0E-07 2.1
Dense, sand permeability 1.3E+05 3.9E+05 40 6.6E-05 2.1
Dense, gravel
1.3E+05 3.9E+05 40 1.0E-02 2.1
permeability
Bulk Mass
Shear modulus Cohesion Permeability
Cohesive Soil modulus B density4
G (kPa) c (kPa) 5 coeff.3 (m/s)
(kPa) (ton/m3)
Soft 1.3E+04 6.5E+04 18.0 1.0E-09 1.3
Medium 6.0E+04 3.0E+05 37.0 1.0E-09 1.5
Stiff 1.5E+05 7.5E+05 75.0 1.0E-09 1.8
1. Where p r is the reference mean effective confining pressure at which soil appropriate soil properties are defined.
2. Friction angles for cohesionless soils are based on Table 7.4 (p.425) of Das, B.M. (1983).
3. Permeability values are based on Fig. 7.6 (p.210) of Holtz and Kovacs (1981).
4. Mass density is based on Table 1.4 (p.10) of Das (1995).
5. Cohesion for cohesive soils are based on Table 7.5 (p.442) of Das (1983).
28
Backbone Curve
Gr
(4.1)
1
r
where Gr is the low-strain shear modulus (see 4.2.3.1), and r satisfies the following equation at
p r :
2 2 sin G r max
f p r (for sands) (4.2a)
3 sin 1 max / r
and,
2 2 sin 2 2 G r max
f p r c (for clays) (4.2b)
3 sin 3 1 max / r
where f is the peak (octahedral) shear strength, is the friction angle, c is the cohesion, and
max is the maximum shear strain (10% is employed in OpenSeesPL).
1
3
xx yy 2 yy zz 2 xx zz 2 6 xy2 6 yz2 6 xz2 1/ 2
2
3
xx yy 2 yy zz 2 xx zz 2 6 xy2 6 yz2 6 xz2
1/ 2
The number of yield surfaces used for the predefined sands and clays is 20.
29
Shear
stress
Shear
strength
Number of yield
surfaces = 5
Shear
strength
Number of yield
surfaces = 0 Number of yield
surfaces = 1
f max
r (4.3a)
G r max f
Or
f r
max (4.3b)
Gr r f
Gr
Gr 1 (4.4)
1 ( )
f max
Take Medium Sand (Table 4.1) as an example, Gr = 75,000 kPa, p r = 80 kPa, = 33°.
30
2 2 sin(33)
f (80 ) = 50.2 kPa (4.5)
3 sin(33)
(75000 )
1 (4.6)
1 (1494 )
max
Figure 4.10 shows the backbone curves at max = 2%, 5% and 10%, based on Eq. (4.6).
60
50
Shear stress (kPa)
40
30
20
γmax = 2%
10 γmax = 5%
γmax = 10%
0
0 1 2 3 4 5
Shear strain (%)
User-defined materials include user-defined sand (U-Sand1 and U-Sand2) with confinement-
dependent material properties, and user-defined clay (U-Clay1 and U-Clay2) with properties
independent of confinement variation. To define the parameters of a user-defined material, click
the list of soil materials and select U-Sand1A, U-Sand1B, U-Clay1, U-Clay2, U-Sand2A, or U-
Sand2B accordingly (Figure 4.8).
31
4.2.3.1 User Defined Sand1A (U-Sand1A)
Saturated Mass Density The saturated mass density of the cohesionless soil.
Reference Pressure The reference mean effective confining pressure ( p r ) at which soil
appropriate soil properties below are defined.
Gmax The reference low-strain shear modulus Gr, specified at a reference mean effective
confining pressure p r .
32
Bmax The reference bulk modulus Br, specified at a reference mean effective confining pressure
p r .
p d p d
G Gr ( ) B Br ( ) (4.7)
p r p r
Peak Shear Strain An octahedral shear strain at which the maximum shear strength is reached,
specified at a reference mean effective confining pressure p r . The suggested values are between
0.001% and 20%.
Friction Angle The friction angle ( ) at peak shear strength in degrees. The suggested values
are between 5 and 65 degrees.
Fluid Mass Density The mass density of the fluid, which is usually 1.0 ton/m3..
Combined Bulk Modulus The combined undrained bulk modulus Bc relating changes in pore
pressure and volumetric strain, may be approximated by:
Bc ≈ Bf / n (4.8)
where Bf is the bulk modulus of fluid phase (2.2×106 kPa for water typically), and n the initial
porosity.
The nonlinear shear stress-strain backbone curve can be defined by specifying a G/Gmax curve
(Figure 4.11). To specify the G/Gmax curve, first enter “number of points defining G/Gmax curve”
and then enter pairs of shear strain and G/Gmax values. The maximum number of points that can
be entered is 13 (the backbone curve becomes horizontal after point 13). In addition:
- If the number of points is zero, then the built-in hyperbolic curve will be used instead.
- If the number of points is 1, the material is elastic-perfectly-plastic.
The user-defined backbone curve is activated if the number of points is greater than zero. In this
case, the user specified friction angle is ignored. Instead, is defined as follows:
3 3 m / p r
sin (3.9)
6 3 m / p r
33
where m is the product of the last modulus and strain pair in the modulus reduction curve.
Therefore, it is important to adjust the backbone curve so as to render an appropriate . If the
resulting is smaller than the phase transformation angle PT, PT is set equal to .
Also remember that improper modulus reduction curves can result in strain softening response
(negative tangent shear modulus), which is not allowed in the current model formulation. Finally,
note that the backbone curve varies with confinement, although the variations are small within
commonly interested confinement ranges. Backbone curves at different confinements can be
obtained using the OpenSees element recorder facility (Mazzoni et al. 2006).
Phase Transformation (PT) Angle The transformation angle (degrees) of the cohesionless soil.
1 ( / PT ) 2
P" c1 (4.10)
1 ( / PT ) 2
where is the stress ratio and PT is the stress ratio along the PT surface (Yang et al. 2003).
1 ( / PT ) 2
P" d1 exp(d 2 d ) (4.11)
1 ( / PT ) 2
where d is the octahedral shear strain accumulated during a dilation phase (Yang et al. 2003).
34
p
y l 2 cos 3 ( ) (3.12)
2 l1
l1 defines the effective confining pressure (e.g., 10 kPa) below which the mechanism is in effect
( l1 is actually p y in Figure 4.12). Smaller values should be assigned to denser sands. l 2 defines
the maximum amount of perfectly plastic shear strain developed at zero effective confinement
during each loading phase ( l 2 is actually s max in Figure 4.12). Smaller values should be assigned
to denser sands.
ry l 3 y (4.13)
l 3 defines the maximum amount of biased perfectly plastic shear strain b accumulated at each
loading phase under biased shear loading conditions, as b= l 2 x l 3 ( ry is actually r , and l 3 is R
in Figure 4.13). Typically, l 3 takes a value between 0.0 and 3.0. Smaller values should be assigned
to denser sands.
d2 0 2 3 5
l1 (kPa) 10 10 5 0
l3 1 1 1 0
35
py p
r ic
v iato ace
De in sp
s tr
a s
smax
Figure 4.12: Initial yield domain at low levels of effective confinement (Yang et al. 2003).
7
e
lure
su rfac 3
Fai rface
2 PT su
(a) 6 1
5 0
p
4
PT su
rface
7
3
(b)
1 2
5 6
0
4
s d r Rs
(c)
3
2 1
3 3
2 3
2 2
Translated
Initial
Enlarged
Figure 4.13: Schematic of constitutive model response showing (a) octahedral stress - effective
confinement p response, (b) octahedral stress - octahedral strain response, and (c)
configuration of yield domain (Yang et al. 2003).
36
4.2.3.2 User Defined Sand1B (U-Sand1B)
Note: All parameters shown in Figure 4.14 are defined at the reference mean confinement p r .
Mass Density The mass density of the cohesionless soil (ρ). The suggested range of values are
between 1 and 3 ton/m3.
Reference Shear Wave Velocity The reference shear wave velocity (Vsr). The suggested range
is between 10 and 6000 m/s. The reference shear modulus G = ρ Vsr2.
Reference Mean Confinement The reference mean confinement. This is the confinement level
at which shear wave velocity and peak shear strain are defined. The suggested range is between
10 kPa or larger.
37
Confinement Dependence Coeff. The confinement dependence coefficient. The suggested
range is between 0.1 and 10.
Initial Lateral/Vertical Confinement Ratio The initial lateral/vertical stress ratio (also known
as coefficient of lateral earth pressure at rest K0). K0 is related to Poisson’s ratio by the following
relation K0 = ν / (1 - ν). The suggested range for K0 is between 0.1 and 0.9.
Cohesion The suggested range is between 0 and 5000000 kPa. See Section 4.2.3.1 for more
information.
Friction Angle The suggested range is between 5 and 65 degrees. See Section 4.2.3.1 for more
information.
Peak Shear Strain The suggested range is between 0.001% and 20%. See Section 4.2.3.1 for
more information.
Number of Yield Surfaces NYS The suggested range is between 0 and 30. In particular, NYS
= 0 dictates an elastic shear response (Cohesion, Friction Angle and Peak Shear Strain are
ignored, see Figure 4.9), NYS = 1 indicates an elastic-perfectly plastic shear response (Peak
Shear Strain is ignored, see Figure 4.9).
Advanced Options:
Use K0 for Elastic Own Weight If checked, users can specify the initial lateral/vertical
confinement ratio K0 which will be used for the application of own weight at the elastic stage
(first run).
Young’s Modulus for Elastic Own Weight The elastic modulus used for the application of
own weight at the elastic stage.
Non-liquefiable clay with shear response properties independent of confinement variation can be
defined as shown in Figure 4.7 and Figure 4.15.
The nonlinear shear stress-strain backbone curve can be defined by specifying a G/Gmax curve
(Figure 4.15). The user-defined backbone curve is activated if the number of points is greater
than zero. In this case, if the user specifies =0, cohesion c will be ignored. Instead, c is defined
by c=sqrt(3)*m/2, where m is the product of the last modulus and strain pair in the modulus
reduction curve. Therefore, it is important to adjust the backbone curve so as to render an
appropriate c.
If the user specifies , this will be ignored. Instead, is defined as follows:
38
3 ( 3 m 2c) / p r
sin (3.14)
6 ( 3 m 2c) / p r
Also remember that improper modulus reduction curves can result in strain softening response
(negative tangent shear modulus), which is not allowed in the current model formulation. Finally,
note that the backbone curve varies with confinement, although the variation is small within
commonly interested confinement ranges. Backbone curves at different confinements can be
obtained using the OpenSees element recorder facility (Mazzoni et al. 2006).
The second type of user-defined clay (U-Clay2) can be defined as shown in Figure 4.16. See
Section 4.2.3.2 for information about parameters defining U-Clay2.
39
Figure 4.15: U-Clay1.
40
Figure 4.16: U-Clay2.
41
Figure 4.17: U-Sand2A.
42
Figure 4.18: U-Sand2B.
The pile zone refers to the pile domain under the ground surface. The material for the pile zone
(Figure 4.19) can be selected from an available menu of cohesionless and cohesive soil materials
including the elastic isotropic material. In addition, user-defined cohesionless and cohesive soil
materials (U-Sand1A, U-Sand1B, U-Clay1, U-Clay2, U-Sand2A, and U-Sand2B) are also
available to choose.
If an elastic isotropic material is selected, the user is requested to specify, Young’s Modulus,
Poisson’s Ratio, Mass Density Permeability of the material used for the pile zone.
43
Figure 4.19: Pile zone material.
The material for the pile-soil interfacing layer (Figure 4.20) can be selected from an available
menu of cohesionless and cohesive soil materials including the elastic isotropic material. In
addition, user-defined cohesionless and cohesive soil materials (U-Sand1A, U-Sand1B, U-
Clay1, U-Clay2, U-Sand2A, and U-Sand2B) are also available to choose.
If an elastic isotropic material is selected, the user is requested to specify, Young’s Modulus,
Poisson’s Ratio, Mass Density Permeability of the material used for the pile-soil interfacing
layer.
The material for the outermost zone (Figure 4.21) can be selected from an available menu of
cohesionless and cohesive soil materials including the elastic isotropic material. In addition,
user-defined cohesionless and cohesive soil materials (U-Sand1A, U-Sand1B, U-Clay1, U-
Clay2, U-Sand2A, and U-Sand2B) are also available to choose.
If an elastic isotropic material is selected, the user is requested to specify, Young’s Modulus,
Poisson’s Ratio, Mass Density Permeability of the material used for the pile-soil interfacing
layer.
44
Figure 4.20: Pile-soil interfacing layer material.
45
5. Mesh Generation
To define the finite element mesh, click Mesh Parameters button in the Model Input window
(Figure 5.1).
Number of Beam Elements above Ground Surface The number of beam elements used for
the pile section above the ground surface.
This section controls mesh refinement along the horizontal direction. Length of each soil
horizontal layer is defined in the left column. Number of mesh elements in each defined is
specified in the column “Number of Mesh Layers”. Note that the first mesh layer is starting
from the center of the mesh when the pile is located and the length of the first mesh layer is equal
to the pile radius. Ratio of Element Length over Next is used to obtain a gradually changing
element size within a layer if Uniform Meshing is unchecked (obviously this option is only
valid if the # of mesh layers is 2 or larger).
This section defines the soil profile (layering) along the vertical direction starting from the
ground surface downwards (looking at the side view from the top downwards. Height
(thickness) of each soil layer is defined in the left column. Number of mesh elements in each
defined is specified in the column “Number of Mesh Layers” (at least equal to 1 to define a soil
profile consisting of a single type of soil). Height (thickness) of this layer must be equal to the
entire soil stratum height. Note that the number of mesh layers in the upper zone (where the pile
foundation is embedded) will automatically define the number of beam column elements of this
pile (below ground surface). As such, it is generally advisable to select an adequate number of
46
mesh layers in this zone. Note: If there is any error during mesh generation, please follow the
error message instructions to adjust the controlling parameters and then try again.
Note: Element size is a parameter that affects frequency content of the ground response. Smaller
size elements (particularly along the soil domain height), will permit higher frequencies (if
present in the input motion) to propagate to the ground surface with more fidelity.
Model Length The length of the soil domain (along the longitudinal direction) to be scaled.
Model Width The width of the soil domain (along the transverse direction) to be scaled.
a) General Definition
47
b) Horizontal Meshing for Single Pile Models
48
c) Vertical Meshing
d) Mesh Scaling
49
6. Pushover & Eigenvalue Analyses
In a pushover or base-shaking analysis, four runs are conducted in sequence in order to achieve
convergence and simulate the actual loading situation.
1) 1st run: Gravity of soil domain is applied in this run; all soil materials are prescribed as
linear during this run.
2) 2nd run: Soil elements are changed to nonlinear if “Nonlinear” is chosen in Analysis
Options (see Section 4.1.1Error! Reference source not found.).
3) 3rd run: Pile elements are added and gravity of the pile structure is applied in this run.
4) 4th run: Pushover or base-shaking analysis is started.
To conduction a pushover analysis, click Pushover and then click Define Pattern in Figure 2.1.
Two types of pushover analyses are available (Figure 6.1): Static and Dynamic Pushover.
The force-based method is used if the Force-Based Method radio button is chosen.
Moment of X The applied bending moment about the longitudinal direction (Mx).
Moment of Y The applied bending moment about the longitudinal direction (My).
Moment of Z The applied bending moment about the longitudinal direction (Mz).
50
Vertical Displacement The displacement applied in the vertical direction.
Rotation around X The applied rotation around the longitudinal axis (X).
Rotation around Y The applied rotation around the transverse axis (Y).
Rotation around Z The applied rotation around the vertical axis (Z).
51
6.1.2 Load Pattern
To conduct a pushover analysis, a load pattern must be defined. The load pattern is shown in
Figure 6.2.
If Static Pushover is chosen, the pushover options include monotonic pushover as well as
pushover by a user-defined loading pattern (U-Push). Please see Section 6.1.2.3 for how to
define a U-Push. In a monotonic pushover, the pushover load/displacement is linearly increasing
with steps. In a dynamic monotonic pushover, users are allowed to define the loading duration.
If Dynamic Pushover is chosen, a Sine Wave loading pattern is also available (Figure 6.2).
To define your own load pattern (U-Push), click U-Push in Figure 6.2. The U-Push window is
shown in Figure 6.3. Click Select/Change Pushover File to change file. The user-defined
pushover file should contain single-column data.
52
Figure 6.3: User-defined pushover load pattern (U-Push)
To run the analysis, click “Save Model & Run Analysis” in Menu “Analyze”.
Upon the user requests to run the analysis, OpenSeesPL will check all the entries defined by the
user to make sure the model is valid. Thereafter, a small window (Figure 6.4) will show the
progress of the analysis.
53
By default, graphical output windows will be opened upon completion of the analysis.
Response time histories and profiles are displayed by X-Y plot using Java Applet. Therefore,
make sure to enable Java Applet in your web browser (Internet Explorer). You may also view the
digital data by clicking on the link under the X-Y plot. If occasionally the graph becomes
crooked, you can click on the Fill button to refresh it.
To zoom in on any region of the plot, select a box with the mouse pointer (Figure 6.5). Start at
the upper left corner of the region you wish to view in more detail and drag downwards and to
the right. To bring the graph to the original scale, click on the "fill" button at the upper right
corner.
(a) Select a box using the mouse pointer… (b) Then release the mouse.
54
(a) Select a box using the mouse pointer… (b) Then release the mouse.
To zoom out, drag the mouse pointer upwards (Figure 6.6). When zooming out, a reference box
is drawn that will represent the current view, and dragging will cause a box to be displayed that
represents the new view. Again, click on the "fill" button at the upper right corner to bring the
graph to the original scale.
To view the pile response, click Pile Response Profiles in Menu Display. The figures show the
response time histories and response profiles of the pile. Seven types of response are available
(Figure 6.7):
Displacement
Acceleration
Rotation
Moment
Shear
Pressure
55
Figure 6.7: Response time histories and profiles for pile.
To view the pile response relationships, click Pile Response Relationships in Menu Display.
The figures show the response relationships of the pile. Two types of response are available
(Figure 6.8):
Load-displacement
Moment-curvature
To zoom-in or zoom-out, use mouse to select a window. Click "fill" to get back to the original
figure.
56
Figure 6.8: Response relationships for pile.
By default, the deformed mesh is for the dynamic analysis (if ‘Due to Seismic Excitation’ is
chosen) or the pushover analysis (if “Due to Pushover’ is chosen). However, the deformed mesh
due to gravity is also available (‘Due to Gravity’ is chosen)
The deformed mesh can be viewed in 3D or 2D (can be selected from a list of 2D cut planes, see
Figure 6.10).
57
To view the animation of any given type, click the “Play Animation” button. The text of the
button will change to “Stop Animation” when the animation is being played. To stop the
animation, click the “Stop Animation” button.
The Scale Factor can be changed to improve the viewing effects. The time between playing two
frames can be defined by specifying the Animation Playing Delay (in millisecond).
Note that the animation will not be played if the current time step is in the last step and “Endless
Playing” is unchecked.
At any time, the deformed mesh can be rotated by dragging the mouse, moved in 4 directions by
pressing keys of LEFT ARROW, RIGHT ARROW, UP ARROW or DOWN ARROW
respectively. The view can be zoomed in (by pressing key ‘F9’), out (by pressing key ‘F10’) or
frame (by pressing key ‘F11’).
58
Figure 6.10: 2D plane (Y = 0) view of the longitudinal displacement contour in the deformed
mesh window.
To conduct an Eigenvalue analysis, click Eigenvalue and then specify Number of Frequencies
in Figure 2.1. And then click Save Model & Run Analysis. Figure 6.11 shows the output
window for an Eigenvalue analysis, which can be accessed by clicking menu Display and then
choosing Deformed Mesh.
59
Figure 6.11: Output for an Eigenvalue analysis
60
7. Base Shaking Analysis
OpenSeesPL employs the Newmark time integration procedure with two user defined
coefficients and (Newmark 1959, Chopra 2004). Standard approaches may be adopted by
appropriate specification of these constants (Figure 7.1). Default values in OpenSeesPL are
0.55, and ( ( + ½)2 ) / 4).
Computations at any time step are executed to a convergence tolerance of 10-4 (Euclidean Norm
of acceleration vector), normalized by the first iteration Error Norm (predictor multi-corrector
approach).
Note: An additional fluid-phase (Chan 1988) time integration parameter is set to 0.6 in the data
file.
61
7.1.2 Input Motion
One, two and three directions of excitations are available: longitudinal, transverse and vertical
directions (Figure 2.1 and Figure 7.2).
The bedrock is assumed to be rigid, the input motion is total motion; Base seismic excitation can
be defined by either of the following two methods:
i) Via a built-in input motion library. This library includes near-fault soil surface motions as well
as long-duration rock outcrop motions recorded during past strong earthquakes worldwide.
ii) ‘U-Shake’, a user-defined input motion (Figure 7.3). The input motion file to be defined
should consist of two columns, Time (seconds) and Acceleration (g), delimited by SPACE(S).
0.00 0.000
0.02 0.005
... .….....
19.98 0.004
20.00 0.000
Note that the user-defined input motion file must be placed in the subfolder “motions/”.
(This subfolder also contains all provided built-in input motion files).
The amplitude of the input motion can be scaled by a factor ranging from 0.01 to 1.0. In addition,
if ‘0.2g sinusoidal motion’ is chosen, the user must specify excitation frequency and number of
cycles (Figure 2.1).
Scale Factor The amplitude of the input motion is multiplied by the Scale Factor. The Scale
Factor may be positive or negative.
Frequency The Frequency (in Hz) has to be specified if harmonic “sinusoidal motion” is chosen
Number of Cycles The Number of Cycles has to be specified if “sinusoidal motion” is chosen.
62
Figure 7.2: Definition of 3D base excitation and boundary conditions.
63
Figure 7.3: User-defined input motion (U-Shake).
Ground Surface Inclination Angle along Longitudinal Direction The inclination angle of the
ground surface along the longitudinal direction (in degrees) (zero degree represents level ground
surface).
64
Whole Model Inclination Angle along Longitudinal Direction The inclination angle in
degrees of the whole model (zero degree represents level ground). For mildly-inclined infinite-
slopes, suggested values are from 0 to 10 degrees.
To view the soil response time histories, click Soil Response Histories in Menu Display.
The figures show the response time histories of the soil domain from the ground surface till the
bottom, at a number of locations which are along the longitudinal direction crossing the pile
center.
Pile response and deformed mesh output are also available in a base shaking analysis.
Please refer to Section 6.1.4.
65
Figure 7.4: Response time histories window.
66
8. Pile Group
Number of Piles The number of piles along X-direction (longitudinal) and Y-direction
(transverse). Note that both numbers do not have to be the same. Therefore, one can easily build
a m by n pile group model in OpenSeesPL. If “1” is entered for both, single pile will be
considered.
Spacing The spacing (specified as a factor of the pile diameter) between pile centers along X-
direction (longitudinal) and Y-direction (transverse). Obviously, the spacing must be greater than
1.
If Fixed is chosen for the pile head, a rigid pile cap will be employed. If Free/Pinned is chosen,
the pinned connection is considered for the pile heads of the pile group.
67
To define the finite element mesh for a pile group model, click Mesh Parameters button in the
Model Input window (Figure 5.1). And click Pile Group in the Horizontal Meshing tab to
define the controlling parameters in the horizontal directions (Figure 8.2).
For General Definition and Vertical Meshing Tabs, please refer to Chapter 5. Figure 8.3 shows
a sample mesh of a 3 x 3 pile group model (half mesh configuration).
68
Figure 8.3: Sample mesh of a 3 by 3 pile group model (half mesh configuration).
69
Figure 8.4: Pile response profiles for a pile group model.
Figure 8.5: Pile response time histories for a pile group model.
70
Figure 8.6: Pile response relationships for a pile group model.
Figure 8.7: Pile response relationships at the pile cap for a pile group model.
71
Figure 8.8: Deformed mesh of a pile group model.
72
Appendix A How to Define the Soil Finite Element
Mesh
Step 1
In the user interface, click Pile Parameters. With reference to Figure 3.1, define the following
parameters according to your preference:
Total Pile Length: Starting from the pile head all the way to the pile tip.
Pile Length above Surface: from pile head to mud-line (ground surface).
Soil Parameters: make sure at least the total “Thickness” of soil layers is defined: This is the
total thickness of the ground stratum from the ground surface all the way down to the base of the
soil mesh. Make sure that the pile tip is within the defined soil domain depth.
Note: Earthquake input motion is imparted along the base of the soil mesh. This base is assumed
to represent rigid bedrock. As such, this input earthquake excitation constitutes total motion
imparted at this Bedrock level.
Step 2
Click Mesh Parameters to define additional meshing parameters (please refer to Chapter 5 and
Figure 5.1).
The finite element mesh created with the above default values is shown in Figure A.1. Examples
of mesh generation are shown in Figures A.2-A.4.
73
a)
b)
Figure A.2: Mesh refinement example 1: a) Change “Num of Slices” to 32; b) the resulting mesh
74
a)
b)
Figure A.3: Mesh refinement example 2: a) Change “Number of Mesh Layers” in the vertical
direction; b) the resulting mesh
75
a)
b)
Figure A.4: Mesh refinement example 3: a) Change meshing controlling parameters in the
horizontal direction; b) the resulting mesh
76
Appendix B Own Weight Application with Dry and
Saturated Soil Cases
Boundary Conditions:
The boundary conditions available in OpenSeesPL include Shear Beam, Rigid Box, and Periodic
Boundary.
1) Shear Beam
In this case, the front and back nodes at any depth move together (horizontal and vertical
directions). The Shear Beam boundary condition, if it’s chosen, is enforced for all runs.
Rollers are used for lateral and base boundaries for all gravity runs. The base nodes are fixed
after the first run.
If Fixed Vert is checked, all nodes at lateral boundaries will be fixed in vertical direction before
the dynamic run.
2) Rigid Box
In gravity runs, lateral boundaries are fixed in both horizontal directions and free in vertical
direction. Rollers are used for base nodes, which will be fixed after the first run.
If Fixed Vert is checked, all nodes at lateral boundaries will be fixed in vertical direction before
the dynamic run.
3) Periodic Boundary
In this case, each node on the front boundary moves the same as the analogous node on the back
boundary (horizontal and vertical directions). The Periodic boundary condition, if it’s chosen, is
enforced for all runs.
For gravity runs, rollers are used for lateral and base boundaries. The base nodes are fixed after
the first run.
If Fixed Vert is checked, all nodes at lateral boundaries will be fixed in vertical direction before
the dynamic run.
At first the defined soil properties are used to set up the soil constitutive model. A static solver is
used and own weight is applied in one step with elastic soil properties. Default is global elastic
modulus (600,000 kPa by default) and global initial lateral/vertical confinement ratio (Ko = 0.9
77
by default) for the entire soil domain. These elastic soil properties are used to define an elastic
stiffness matrix (Kmatrix1). A default convergence tolerance of 0.0001 is used (displacement
norm), which the user can specify in the OpeSees Parameters section (from Analysis Options).
Base: Rollers are used to prevent vertical displacement, but lateral deformation is allowed.
The actual defined soil properties in every part of the mesh are activated, and nonlinear (if
specified) properties are activated as well.
The static solver is used, and Kmatrix1 is used for convergence. A convergence tolerance of
0.0001 is used (displacement norm). The boundary conditions for this step remain those of BC1.
A new mass and stiffness matrix is built based on the latest tangent soil stress-strain state, and
the linear properties of the beam column elements. A convergence tolerance of 0.0001 is used
(displacement norm). The load is applied in 20 steps by default (the user can modify this value in
the OpeSees Parameters section (from Analysis Options). The stiffness matrix is not updated.
The boundary conditions for this step remain those of BC1.
4) Solution phase
Solution is started with a stiffness matrix based on the latest soil and beam column stress-strain
state. Four different analysis scenarios are possible:
The static solver is used with a convergence tolerance of 0.0001 that the user can modify in the
OpeSees Parameters section (from Analysis Options) (displacement norm).
Boundary conditions for this case are: Default is fixed boundaries everywhere, but the user can
change that to Shear Beam or Periodic Boundary.
In this case, a dynamic solver is used (modified Newton-Raphson) with the time integration
parameters γ = 0.6 and β = 0.3025, and the actual user specified time step. Note that the user can
also modify the Rayleigh mass and stiffness proportional viscous damping parameters (which are
set by default to 2% at the frequencies of 1 Hz and 6 Hz).
78
After the dynamic load has been applied, analysis can proceed for a user specified number of
seconds so that the “free vibration response” can be assessed if so desired.
Boundary conditions for this case are: Default is fixed boundaries everywhere, but the user can
change that to Shear Beam or Periodic Boundary.
In this case, a dynamic solver is used (modified Newton-Raphson) with the time integration
parameters γ = 0.6 and β = 0.3025, and the actual user specified time step. The convergence
tolerance of 0.0001 is the default but the user can modify this value in the OpeSees Parameters
section (from Analysis Options) (displacement norm). Note that the user can also modify the
Rayleigh mass and stiffness proportional viscous damping parameters which are set by default to
2% at the frequencies of 1 Hz and 6 Hz).
After the dynamic load has been applied, analysis can proceed for a user specified number of
seconds so that the “free vibration response” can be assessed if so desired.
Boundary conditions for this case are: Default is fixed boundaries everywhere, with the base
moving according to the applied base excitation. The user might wish to activate alternate
boundary conditions along the lateral boundaries in the form of Shear beam boundary conditions
where the front and back nodes at any depth move together, or a periodic boundary condition
where each node on the front boundary moves the same as the analogous node on the back
boundary (and the vertical is free, but can be fixed by the user).
Eigenvalue analysis:
In this step the mass and stiffness matrices corresponding to the latest stress-strain state (after
application of own weight of the beam-column elements) are used to compute natural
frequencies and mode shapes (using the static solver).
Boundary conditions for this case are: Default is fixed boundaries everywhere, with the base
moving according to the applied base excitation. The user might wish to activate alternate
boundary conditions along the lateral boundaries in the form of Shear beam boundary conditions
where the front and back nodes at any depth move together, or a periodic boundary condition
where each node on the front boundary moves the same as the analogous node on the back
boundary (and the vertical is free, but can be fixed by the user).
Dry soil case with mildly inclined ground and soil with water table specified
1) Application of soil own weight with elastic soil properties
A dynamic solver is used and own weight is applied in 5 steps (time step is set to 50,000 secs
and gamma γ and beta β parameters are set to 1.5 and 1 in order to obtain a static solution) with
79
elastic soil properties (elastic modulus = 600,000 kPa and initial lateral/vertical confinement
ratio Ko = 0.9 by default) and a default global very large permeability coefficient (100 m/s by
default; the permeability will be changed to the user-specified value before the dynamic run).
Default is global elastic modulus (600,000 kPa by default) and global initial lateral/vertical
confinement ratio (Ko = 0.9 by default) for the entire soil domain. These specified global values
will be used for the top soil layer. For all other soil layers including the pile zone and the
interfacing zone, the elastic modulus employed is equal to the above global value (600,000 kPa
by default) times the ratio of the mass density of the current soil layer over the top soil layer.
These elastic soil properties are used to define an elastic stiffness matrix (Kmatrix1). A default
convergence tolerance of 0.0001 is used (displacement norm), which the user can specify in the
OpeSees Parameters section (from Analysis Options).
Base: Rollers are used to prevent vertical displacement, but lateral deformation is allowed.
2) Model inclination
If the model is inclined, an extra run for applying the horizontal gravity factor is added. The
horizontal gravity factor is applied at the based nodes as acceleration input (the base nodes have
to be fixed before this run). The vertical gravity factor is applied at the first run (through the soil
element body force factor).
The actual defined soil properties in every part of the mesh are activated, and nonlinear (if
specified) properties are activated as well.
The dynamic solver is used (similar to item 1 above), and Kmatrix1 is used for convergence.
Own weight is applied in 5 steps (time step is set to 50,000 secs). A convergence tolerance of
0.0001 is used (displacement norm). The boundary conditions for this step remain those of BC1.
A new mass and stiffness matrix is built based on the latest tangent soil stress-strain state, and
the linear properties of the beam column elements. A convergence tolerance of 0.0001 is used
(displacement norm). The load is applied in 20 steps by default (the user can modify this value in
the the OpeSees Parameters section (from Analysis Options). The stiff matrix is not updated.
The dynamic solver is used (similar to Section 2, and 5 time steps are allowed with no additional
input excitation to ensure convergence to a stable static solution). The boundary conditions for
this step remain those of BC1.
80
5) Solution Phase
Solution is started with a stiffness matrix based on the latest soil and beam column stress-strain
state. Four different analysis scenarios are possible:
The dynamic solver is used (similar to item 1 above) with a convergence tolerance of 0.0001
(displacement norm) that the user can modify in the OpeSees Parameters section (from
Analysis Options).
Boundary conditions for this case are: Default is fixed boundaries everywhere, but the user can
change that to Shear Beam or Periodic Boundary.
In this case, a dynamic solver is used (modified Newton-Raphson) with the time integration
parameters γ = 0.6 and β = 0.3025, and the actual user specified time step. Note that the user can
also modify the Rayleigh mass and stiffness proportional viscous damping parameters (which are
set by default to 2% at the frequencies of 1 Hz and 6 Hz).
After the dynamic load has been applied, analysis can proceed for a user specified number of
seconds so that the “free vibration response” can be assessed if so desired.
Boundary conditions for this case are: Default is fixed boundaries everywhere, but the user can
change that to Shear Beam or Periodic Boundary.
In this case, a dynamic solver is used (modified Newton-Raphson) with the time integration
parameters γ = 0.6 and β = 0.3025, and the actual user specified time step. The convergence
tolerance of 0.0001 is the default but the user can modify this value (displacement norm) in the
OpeSees Parameters section (from Analysis Options). Note that the user can also modify the
Rayleigh mass and stiffness proportional viscous damping parameters which are set by default to
2% at the frequencies of 1 Hz and 6 Hz).
After the dynamic load has been applied, analysis can proceed for a user specified number of
seconds so that the “free vibration response” can be assessed if so desired.
Boundary conditions for this case are: Default is fixed boundaries everywhere, with the base
moving according to the applied base excitation. The user might wish to activate alternate
boundary conditions along the lateral boundaries in the form of Shear beam boundary conditions
where the front and back nodes at any depth move together, or a periodic boundary condition
where each node on the front boundary moves the same as the analogous node on the back
81
boundary (and the vertical is free, but can be fixed by the user).
Eigenvalue analysis:
In this step the mass and stiffness matrices corresponding to the latest stress-strain state (after
application of own weight of the beam-column elements) are used to compute natural
frequencies and mode shapes (using a dynamic solver).
Boundary conditions for this case are: Default is fixed boundaries everywhere, with the base
moving according to the applied base excitation. The user might wish to activate alternate
boundary conditions along the lateral boundaries in the form of Shear beam boundary conditions
where the front and back nodes at any depth move together, or a periodic boundary condition
where each node on the front boundary moves the same as the analogous node on the back
boundary (and the vertical is free, but can be fixed by the user).
82
Appendix C Benchmark Linear Finite Element
Analysis of Laterally Loaded Single Pile Using
OpenSees & Comparison with Analytical Solution
Introduction
In this study:
I) The response of a laterally loaded pile obtained using the OpenSeesPL
interface is compared with the analytical elastic solution proposed by
Abedzadeh and Pak (2004). Detailed information about the analytical elastic
solution is provided in Appendix C-I (please see this end of Appendix C).
II) Based on the linear analysis presented below, nonlinear soil response is
addressed in Appendix C-II (please see this end of Appendix C).
The pile employed in the OpenSees simulation is circular with a diameter of 16" (radius a
= 8") while the one for the analytical elastic solution is a cylindrical pipe pile of the same
radius and a wall thickness h = 0.1a. Both cases have the same pile length l = 33.3 ft (l/a
= 50). The cross-sectional moment of inertia of the pipe pile I = πa3h = 1286.8 in4, which
will be used for the circular pile in the OpenSees simulation.
In summary, the geometric and elastic material properties of the pile are listed below:
Radius a = 8"
Pile length l = 33.3 ft
Young’s Modulus of Pile Ep = 29000 ksi
Moment of Inertia of Pile I = 1286.8 in4
Soil Domain
The ratio of Young’s Modulus of Pile (Ep) to the Shear Modulus of Soil (Gs):
Ep/Gs = 3634 (which will be used later to obtain the analytical elastic solution by
83
interpolation).
Lateral Load
The pile head (free head condition), which is located at the ground surface, is subjected to
a horizontal load (H) of 31.5 kips.
The floating pile is modeled by beam-column elements, and rigid beam-column elements
are used to model the pile size (diameter).
The lateral load is applied at the pile head (ground level) in x (longitudinal) direction.
The above simulations were performed using OpenSeesPL (Lu et al., 2006).
The pile head deflection and the maximum bending moment from OpenSees and the
elastic solution are also listed in Table C.1. In general, the numerical results match well
with the analytical elastic solution. The pile head deflection from the 20-node element
mesh (0.043") is almost identical to the elastic solution (0.042").
84
For nonlinear run, please see Appendix C-II.
a) Isometric view
85
Table C.1: Comparison of OpenSees results and the analytical elastic solution.
OpenSees Results
Elastic solution by
20-node Abedzadeh and Pak (2004)
8-node element
element
Pile head deflection (in) 0.039 0.043 0.042
Maximum moment Mmax (kip-ft) 30 31 27
Depth where Mmax occurs (ft) 2.87 2.87 2.7
10
Pile depth (ft)
15
20
25
86
0
10
Pile depth (ft)
15
20
25
87
Appendix C-I: Elastic Solution of the Response of a Laterally Loaded
Pile in a Semi-Infinite Soil Medium with Constant Modulus along Depth
(For details, please see: Farzad Abedzadeh and Y. S. Pak. (2004). “Continuum
Mechanics of Lateral Soil–Pile Interaction”, Journal of Engineering Mechanics, Vol.
130, No. 11, November, pp. 1309-1318).
Consider a flexible cylindrical pipe pile of radius a, length l, a wall thickness h « a (note
that the moment of inertia I = πa3h). The pile is assumed to be fully embedded in a
homogenous, isotropic, linearly elastic half-space with a shear modulus Gs and a
Poisson’s ratio νs = 0.25.
Using Eqs. (78)-(83) and Figure 9 of the above reference, the pile response (h/a=0.1,
l/a=50) under an applied pure pile-head horizontal load is shown in Figure C.4 and
Figure C.5, where,
Ep – Young’s Modulus of Pile
Gs – Shear Modulus of Soil
w – Pile deflection (in)
H – Horizontal load (kip)
z – Pile depth (ft)
88
Normalized Deflection w /(HaG s )
0 0.02 0.04 0.06 0.08 0.1
0
10
15
Normalized Depth z/a
20
25
30
35
40
Ep/Gs = 25000
Ep/Gs = 2500
45
50
Figure C.4: Sample pile deflection (h/a=.1, l/a=50) under an applied pure pile-head
horizontal load (Abedzadeh and Pak, 2004).
89
Normalized Moment M /(Ha )
0 0.5 1 1.5 2 2.5 3
0
10
15
Normalized Depth z/a
20
25
30
35
40
Ep/Gs = 25000
Ep/Gs = 2500
45
50
Figure C.5: Sample pile bending moment (h/a=.1, l/a=50) under an applied pure pile-
head horizontal load (Abedzadeh and Pak, 2004).
90
Appendix C-II: Nonlinear Response of the Single Pile Model
In the nonlinear run, the same material properties of the linear run are employed except
the soil now assumed to be a clay material with a maximum shear strength or cohesion =
5.1 psi, in the range of a Medium Clay. This maximum shear strength is achieved at a
specified strain γmax = 10%.
The lateral load (H) is applied at an increment of 0.7875 kips and the final load is 94.5
kips (= 3 x 31.5 kips). The 8-node brick element mesh is employed in this nonlinear
analysis (Figure C.1).
Simulation Results
Figure C.6 shows the load-deflection curve for the nonlinear run, along with the linear
result (for the 8-node brick element mesh; the final lateral load is also extended to 94.5
kips) as described in the previous sections for comparison. It is seen from Figure C.6 that
nearly linear behavior is exhibited in the nonlinear run for only low levels of applied
lateral load (less than 10 kips).
100
90
80
70
60
Load (kips)
50
40
30
20
10
Linear
Nonlinear
0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
Pile head deflection (in)
Figure C.6: Comparison of the load-deflection curves for the linear and nonlinear runs.
91
The pile deflection profiles for both linear and nonlinear cases are displayed in Figure 7.
For comparison, the linear and nonlinear responses at the lateral load of 31.5 kips, 63 kips
(= 2 x 31.5), and 94.5 kips (= 3 x 31.5) are shown (Figure C.7). The bending moment
profiles for the 3 load levels are shown in Figure C.8a-c.
10
Depth (ft)
15
20
25
30 Linear
Nonlinear
35
0 0.1 0.2 0.3 0.4 0.5
Pile deflection (in)
a) H = 31.5 kips
10
Depth (ft)
15
20
25
30 Linear
Nonlinear
35
0 0.1 0.2 0.3 0.4 0.5
Pile deflection (in)
b) H = 63 kips
Figure C.7: Comparison of the pile deflection profiles for the linear and nonlinear runs.
92
0
10
Depth (ft) 15
20
25
30 Linear
Nonlinear
35
0 0.1 0.2 0.3 0.4 0.5
Pile deflection (in)
c) H = 94.5 kips
The pile head deflections and the maximum bending moments for both linear and
nonlinear analyses are listed in Table C.2. The stress ratio contour fill of the nonlinear
run is displayed in Figure C.9.
0
10
Depth (ft)
15
20
25
30 Linear
Nonlinear
35
0 50 100 150 200
Bending moment (kip−ft)
a) H = 31.5 kips
Figure C.8: Comparison of the pile bending moment profiles for the linear and nonlinear
runs.
93
0
10
Depth (ft)
15
20
25
30 Linear
Nonlinear
35
0 50 100 150 200
Bending moment (kip−ft)
b) H = 63 kips
10
Depth (ft)
15
20
25
30 Linear
Nonlinear
35
0 50 100 150 200
Bending moment (kip−ft)
c) H = 94.5 kips
94
Table C.2: OpenSees simulation results for the linear and nonlinear runs.
Figure C.9: Stress ratio contour fill of the nonlinear run at different load levels (red color
shows yielded soil elements).
95
Appendix D Finite Element Analysis of
Arkansas Test Series Pile #2 Using Opensees (with
LPILE Comparison)
Introduction
In this study, we conduct a finite element simulation of Pile No. 2 of the Arkansas test
series (Alizadeh and Davisson 1970) using the OpenSeesPL interface. This pipe pile is
subjected to lateral loads. Comparison with LPILE is also included in Appendix D-I
(please see the end of Appendix D).
Pile Data
The pile employed in the OpenSees simulation is circular with a diameter of 16" (radius a
= 8") while the one for the experimental test is a cylindrical pipe pile of the same radius
and a wall thickness h = 0.312". The cross-sectional moment of inertia of the pipe pile I =
838.2 in4 (Bowles 1988, pages 777-778), which will be used for the circular pile in the
OpenSees simulation.
The geometric and elastic material properties of the pile are listed below (Bowles 1988):
Soil Domain
In this section, the pile is embedded in a uniform soil layer (pile top is 0.1' above the
ground line). Linear and nonlinear soil responses are investigated. The Medium density
(relative) granular soil type (Lu et al. 2006) is selected in this initial attempt. The material
properties of the soil are listed below:
At the reference confinement of 80 kPa (or 11.6 psi), the Shear Modulus of Soil Gs =
10.88 ksi and the Bulk Modulus of Soil B = 29 ksi (i.e., Poisson’s ratio νs = 0.33), see Lu
et al. 2006.
Submerged Unit Weight γ' = 62.8 pcf (Bowles 1988)
For nonlinear analysis, the Friction Angle φ = 32° (Bowles 1988) and the peak shear
stress occurs at a shear strain γmax = 10% (at the 11.6 psi confinement)
96
Lateral Load
The pile head (with a free head condition), which is 0.1' above the ground surface, is
subjected to horizontal loads (H) of 21 kips, 31.5 kips and 43 kips (Bowles 1988).
The floating pile is modeled by beam-column elements, and rigid beam-column elements
are used to model the pile size (diameter).
The lateral load is applied at the pile head (ground level) in x (longitudinal) direction.
The above simulations were performed using OpenSeesPL (Lu et al. 2006).
Simulation Results
The pile deflections at the ground line and the maximum bending moments for the linear
and nonlinear analyses are listed in Table D.1, along with the experimental measurements
for comparison (Alizadeh and Davisson 1970; Bowles 1988).
Figure D.2 shows the load-deflection curve for the linear and nonlinear runs. Comparison
of the pile deflection profiles for the linear and nonlinear analyses are displayed in Figure
D.3a-c. The bending moment profiles for the 3 load levels are shown in Figure D.4a-c,
along with the observed for comparison (Alizadeh and Davisson 1970). The stress ratio
contour fill of the nonlinear run is displayed in Figure D.5.
97
(a) Isometric view
98
Table D.1: OpenSees Simulation Results and Experimental Measurements.
Pile
Max. bending
deflection at Mmax Profile
Analysis type moment Mmax
ground line depth (ft) displays
(kip-ft)
(in)
H = 21 kips
Experimental 0.17 62 4
Case 1 Linear soil 0.085 35.1 3.1 Figures
Case 2 Nonlinear soil 0.31 70.5 6.8 3a & 4a
H = 31.5 kips
Experimental 0.26 85 5
Case 3 Linear soil 0.13 52.6 3.1 Figures
Case 4 Nonlinear soil 0.56 115.5 6.8 3b & 4b
H = 43 kips
Experimental 0.4 120 5
Case 5 Linear soil 0.17 70.1 3.1 Figures
Case 6 Nonlinear soil 0.89 164.7 6.8 3c & 4c
60
50
40
Load (kips)
30
20
10
Linear
Nonlinear
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
Pile deflection at ground line (in)
Figure D.2: Comparison of the load-deflection curves for the linear and nonlinear runs.
99
0
Depth (ft) 10
15
20
25
Linear
Nonlinear
30
0 0.2 0.4 0.6 0.8
Pile deflection (in)
d) H = 21 kips
10
Depth (ft)
15
20
25
Linear
Nonlinear
30
0 0.2 0.4 0.6 0.8
Pile deflection (in)
e) H = 31.5 kips
Figure D.3: Comparison of the pile deflection profiles for the linear and nonlinear runs.
100
0
10
Depth (ft)
15
20
25
Linear
Nonlinear
30
0 0.2 0.4 0.6 0.8
Pile deflection (in)
f) H = 43 kips
10
Depth (ft)
15
20
25
Linear
30 Nonlinear
Experimental
35
0 50 100 150
Bending moment (kip−ft)
a) H = 21 kips
Figure D.4: Comparison of the pile bending moment profiles for the linear and nonlinear
runs.
101
0
10
Depth (ft)
15
20
25
Linear
30 Nonlinear
Experimental
35
0 50 100 150
Bending moment (kip−ft)
b) H = 31.5 kips
10
Depth (ft)
15
20
25
Linear
30 Nonlinear
Experimental
35
0 50 100 150
Bending moment (kip−ft)
c) H = 43 kips
102
a) First step b) H = 21 kips
Figure D.5: Stress ratio contour fill of the nonlinear run at different load levels (red color
shows yielded soil elements).
103
Appendix D-I: Comparison with LPILE
In the LPILE run, a p-y modulus of 90 psi is employed (p-y multiplier = 1.0). All other
properties are the same as described earlier.
0
10
Depth (ft)
15
20
a) H = 21 kips
0
10
Depth (ft)
15
20
b) H = 31.5 kips
Figure D.6: Comparison of the pile deflection profiles for the linear and nonlinear runs.
104
0
10
Depth (ft)
15
20
c) H = 43 kips
10
20
Depth (ft)
30
a) H = 21 kips
Figure D.7: Comparison of the pile bending moment profiles for the linear and nonlinear
runs.
105
0
10
Depth (ft) 20
30
b) H = 31.5 kips
10
20
Depth (ft)
30
c) H = 43 kips
106
Appendix E Finite Element Analysis of
Standard CalTrans 16" CIDH Pile Using Opensees for
General Comparison with LPILE (with Default P-Y
Multiplier = 1.0)
Introduction
In this study, we conduct a finite element simulation of the standard Caltran 16" CIDH
pile using the 3D OpenSeesPL interface. The simulated pile responses are compared with
LPILE results.
Pile Data
The geometric and elastic material properties of the pile are listed below:
Diameter D = 16"
Pile length l = 35 ft
Moment of Inertia of Pile I = 850 in4
Young’s Modulus of Pile Ec = 4030 ksi
In this initial study, the pile was modeled to remain linear (also in view of the applied
load levels).
Soil Domain
Linear and nonlinear soil responses are investigated. The Medium relative-density
granular soil type (Lu et al. 2006) is selected in the analyses. The material properties of
the soil are listed below:
At the reference confinement of 80 kPa (or 11.6 psi), the Shear Modulus of Soil Gs =
10.88 ksi and the Bulk Modulus of Soil B = 29 ksi (i.e., Poisson’s ratio νs = 0.33), see Lu
et al. 2006.
For nonlinear analysis, the Friction Angle φ = 33° (given by CalTrans) and the peak shear
stress occurs at a shear strain γmax = 10% (at the 11.6 psi confinement). The parameter
γmax along with the shear modulus define the nonlinear soil stress-strain curve. Other
values of γmax should be explored in the future.
107
Lateral Load
Two load cases (Table 1) are studied. The loads are applied at the pile head.
The floating pile is modeled by beam-column elements (Mazzoni et al. 2006), and rigid
beam-column elements are used to model the pile size (diameter).
The lateral load is applied at the pile head (ground level) in x (longitudinal) direction.
The above simulations were performed using OpenSeesPL (Lu et al. 2006).
Simulation Results
The pile head deflections and the maximum bending moments for the linear and
nonlinear analyses are listed in Table 2, along with LPILE results for comparison (see
Appendix for partial output of LPILE results).
Figures C.2-C.5 show comparisons of the pile deflection, rotation, bending moment and
shear force profiles, respectively, for load case 1. Figures C.6-C.9 show comparisons of
108
the pile deflection, rotation, bending moment and shear force profiles, respectively, for
load case 2. The stress ratio contour fill of the nonlinear runs for load cases 1 & 2 are
displayed in Figures C.10 & C.11.
109
Table E.2: CalTrans CIDH Pile OpenSees Simulation and LPILE Results.
Max. bending
Analysis
Pile head Mmax Profile
moment Mmax
typedeflection (in) depth (ft) displays
(kip-ft)
LPILE 0.24 -48.2 0
Load Case 1
Linear soil 0.038 -20.8 0 Figures 2
Fixed Head
Nonlinear &4
H = 16 kips 0.092 -32.3 0
soil
Load Case 2 LPILE -0.094 -100 0
Free Head Linear soil -0.06 -96.7 0
Figures 6
M = -100 kip-ft
Nonlinear &8
applied opposite to -0.094 -96.9 0
soil
shear
10
Depth (ft)
15
20
25
110
0
10
Depth (ft)
15
20
25
10
Depth (ft)
15
20
25
111
0
10
Depth (ft)
15
20
25
0
OpenSees Linear Soil
5 OpenSees Nonlinear Soil
LPILE
10
Depth (ft)
15
20
25
30
35
−0.1 −0.05 0 0.05 0.1
Pile deflection (in)
112
0
10
Depth (ft)
15
20
25
10
Depth (ft)
15
20
25
113
0
Depth (ft) 10
15
20
25
Figure E.10: Stress ratio contour fill for load case 1 (red color shows yielded soil
elements).
Figure E.11: Stress ratio contour fill for load case 2 (red color shows yielded soil
elements).
114
Appendix F Finite Element Analysis of Caltrans
42" CIDH Pile Using OpenSees for General Comparison
with LPILE (with Default P-Y Multiplier = 1.0)
Introduction
In this study, we conduct a finite element simulation of a CalTrans 42" CIDH pile using
the 3D OpenSeesPL interface. The simulated pile responses are compared with LPILE
results.
Pile Data
The geometric and elastic material properties of the pipe pile are listed below:
In this initial study, the pile was modeled to remain linear (also in view of the applied
load levels).
Soil Domain
Linear and nonlinear soil responses are investigated. The Medium relative-density
granular soil type (Lu et al. 2006) is selected in the analyses. The material properties of
the soil are listed below:
At the reference confinement of 80 kPa (or 11.6 psi), the Shear Modulus of Soil Gs =
10.88 ksi and the Bulk Modulus of Soil B = 29 ksi (i.e., Poisson’s ratio νs = 0.33), see Lu
et al. 2006.
For nonlinear analysis, the Friction Angle φ = 33° and the peak shear stress occurs at a
shear strain γmax = 10% (at the 11.6 psi confinement). The parameter γmax along with the
shear modulus define the nonlinear soil stress-strain curve. Other values of γmax should be
explored in the future.
115
Lateral Load
A total of six load cases (Table 1) are studied. The loads are applied at the pile head.
The floating pile is modeled by beam-column elements (Mazzoni et al. 2006), and rigid
beam-column elements are used to model the pile size (diameter).
The lateral load is applied at the pile head (ground level) in x (longitudinal) direction.
The above simulations were performed using OpenSeesPL (Lu et al. 2006).
Simulation Results
Figures D.2-D.5 show comparisons of the pile deflection, rotation, bending moment and
shear force profiles, respectively, for the fixed-head condition (load cases 1, 2 & 3), along
with LPILE results for comparison. Figures D.6-D.9 show comparisons of the pile
deflection, rotation, bending moment and shear force profiles, respectively, for the free-
head condition (load cases 4, 5 & 6), also along with LPILE results for comparison. The
116
stress ratio contour fill of the nonlinear runs for the fixed and free head conditions are
displayed in Figures D.10 & D.11.
Comparisons of the linear and nonlinear responses using OpenSees are shown in
Appendix (Figures D.12-D.19).
117
0
10
Depth (ft)
15
20
OpenSees Nonlinear Soil, 64 kips
25 OpenSees Nonlinear Soil, 128 kips
OpenSees Nonlinear Soil, 256 kips
30 LPILE, 64 kips
LPILE, 128 kips
LPILE, 256 kips
35
−0.1 0 0.1 0.2 0.3 0.4 0.5 0.6
Pile deflection (in)
Figure F.2: Comparison of pile deflection profiles for the fixed-head condition.
10
Depth (ft)
15
20
OpenSees Nonlinear Soil, 64 kips
25 OpenSees Nonlinear Soil, 128 kips
OpenSees Nonlinear Soil, 256 kips
30 LPILE, 64 kips
LPILE, 128 kips
LPILE, 256 kips
35
−2.5 −2 −1.5 −1 −0.5 0 0.5
Rotation (rad) −3
x 10
Figure F.3: Comparison of pile rotation profiles for the fixed-head condition.
118
0
10
Depth (ft)
15
20
OpenSees Nonlinear Soil, 64 kips
25 OpenSees Nonlinear Soil, 128 kips
OpenSees Nonlinear Soil, 256 kips
30 LPILE, 64 kips
LPILE, 128 kips
LPILE, 256 kips
35
−2500 −2000 −1500 −1000 −500 0 500 1000
Bending moment (kip−ft)
Figure F.4: Comparison of bending moment profiles for the fixed-head condition.
10
Depth (ft)
15
20
OpenSees Nonlinear Soil, 64 kips
25 OpenSees Nonlinear Soil, 128 kips
OpenSees Nonlinear Soil, 256 kips
30 LPILE, 64 kips
LPILE, 128 kips
LPILE, 256 kips
35
−100 −50 0 50 100 150 200 250 300
Shear force (kips)
Figure F.5: Comparison of shear force profiles for the fixed-head condition.
119
0
10
Depth (ft)
15
20
OpenSees Nonlinear Soil, 64 kips
25 OpenSees Nonlinear Soil, 128 kips
OpenSees Nonlinear Soil, 256 kips
30 LPILE, 64 kips
LPILE, 128 kips
LPILE, 256 kips
35
−0.5 0 0.5 1 1.5 2
Pile deflection (in)
Figure F.6: Comparison of pile deflection profiles for the free-head condition.
10
Depth (ft)
15
20
OpenSees Nonlinear Soil, 64 kips
25 OpenSees Nonlinear Soil, 128 kips
OpenSees Nonlinear Soil, 256 kips
30 LPILE, 64 kips
LPILE, 128 kips
LPILE, 256 kips
35
−10 −8 −6 −4 −2 0 2
Rotation (in) −3
x 10
Figure F.7: Comparison of pile rotation profiles for the free-head condition.
120
0
10
Depth (ft)
15
20
Figure F.8: Comparison of bending moment profiles for the free-head condition.
10
Depth (ft)
15
20
Figure F.9: Comparison of shear force profiles for the free-head condition.
121
a) lateral load = 64 kips (left: plan view; right: side view)
b) lateral load = 128 kips (left: plan view; right: side view)
c) lateral load = 256 kips (left: plan view; right: side view)
Figure F.10: Stress ratio contour fill of the nonlinear run for the fixed-head condition (red
color shows yielded soil elements).
122
a) lateral load = 64 kips (left: plan view; right: side view)
b) lateral load = 128 kips (left: plan view; right: side view)
c) lateral load = 256 kips (left: plan view; right: side view)
Figure F.11: Stress ratio contour fill of the nonlinear run for the free-head condition (red
color shows yielded soil elements).
123
Appendix F-I: OpenSees Simulation Results
10
Depth (ft)
15
20
OpenSees Linear Soil, 64 kips
25 OpenSees Linear Soil, 128 kips
OpenSees Linear Soil, 256 kips
30 OpenSees Nonlinear Soil, 64 kips
OpenSees Nonlinear Soil, 128 kips
OpenSees Nonlinear Soil, 256 kips
35
0 0.05
0.1 0.15 0.2 0.25
Pile deflection (in)
Figure F.12: Comparison of pile deflection profiles for the fixed-head condition.
10
Depth (ft)
15
20
OpenSees Linear Soil, 64 kips
25 OpenSees Linear Soil, 128 kips
OpenSees Linear Soil, 256 kips
30 OpenSees Nonlinear Soil, 64 kips
OpenSees Nonlinear Soil, 128 kips
OpenSees Nonlinear Soil, 256 kips
35
−14 −12 −10 −8 −6 −4 −2 0 2
Rotation (rad) −4
x 10
Figure F.13: Comparison of pile rotation profiles for the fixed-head condition.
124
0
10
Depth (ft)
15
20
OpenSees Linear Soil, 64 kips
25 OpenSees Linear Soil, 128 kips
OpenSees Linear Soil, 256 kips
30 OpenSees Nonlinear Soil, 64 kips
OpenSees Nonlinear Soil, 128 kips
OpenSees Nonlinear Soil, 256 kips
35
−2000 −1500 −1000 −500 0 500
Bending moment (kip−ft)
Figure F.14: Comparison of bending moment profiles for the fixed-head condition.
10
Depth (ft)
15
20
OpenSees Linear Soil, 64 kips
25 OpenSees Linear Soil, 128 kips
OpenSees Linear Soil, 256 kips
30 OpenSees Nonlinear Soil, 64 kips
OpenSees Nonlinear Soil, 128 kips
OpenSees Nonlinear Soil, 256 kips
35
−50 0 50 100 150 200 250 300
Shear force (kips)
Figure F.15: Comparison of shear force profiles for the fixed-head condition.
125
0
10
Depth (ft)
15
20
OpenSees Linear Soil, 64 kips
25 OpenSees Linear Soil, 128 kips
OpenSees Linear Soil, 256 kips
30 OpenSees Nonlinear Soil, 64 kips
OpenSees Nonlinear Soil, 128 kips
OpenSees Nonlinear Soil, 256 kips
35
−0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
Pile deflection (in)
Figure F.16: Comparison of pile deflection profiles for the free-head condition.
10
Depth (ft)
15
20
OpenSees Linear Soil, 64 kips
25 OpenSees Linear Soil, 128 kips
OpenSees Linear Soil, 256 kips
30 OpenSees Nonlinear Soil, 64 kips
OpenSees Nonlinear Soil, 128 kips
OpenSees Nonlinear Soil, 256 kips
35
−5 −4 −3 −2 −1 0 1
Rotation (in) −3
x 10
Figure F.17: Comparison of pile rotation profiles for the free-head condition.
126
0
10
Depth (ft)
15
20
OpenSees Linear Soil, 64 kips
25 OpenSees Linear Soil, 128 kips
OpenSees Linear Soil, 256 kips
30 OpenSees Nonlinear Soil, 64 kips
OpenSees Nonlinear Soil, 128 kips
OpenSees Nonlinear Soil, 256 kips
35
−200 0 200 400 600 800 1000 1200 1400
Bending moment (kip−ft)
Figure F.18: Comparison of bending moment profiles for the free-head condition.
10
Depth (ft)
15
20
Figure F.19: Comparison of shear force profiles for the free-head condition.
127
Appendix G Finite Element Analysis of
Standard CalTrans 16" CIDH Pile Subjected to Axial
Load
Introduction
In this study, we conduct a finite element simulation of the standard Caltran 16" CIDH
pile using the 3D OpenSeesPL interface. The simulated pile is subjected to axial load.
Pile Data
The geometric and elastic material properties of the pile are listed below:
Diameter D = 16"
Pile length l = 35 ft
Moment of Inertia of Pile I = 850 in4
Young’s Modulus of Pile Ec = 4030 ksi
In this initial study, the pile was modeled to remain linear (also in view of the applied
load levels).
Soil Domain
Nonlinear soil response is investigated. The Medium relative-density granular soil type
(Lu et al. 2006) is selected in the analyses. The material properties of the soil are listed
below:
At the reference confinement of 80 kPa (or 11.6 psi), the Shear Modulus of Soil Gs =
10.88 ksi and the Bulk Modulus of Soil B = 29 ksi (i.e., Poisson’s ratio νs = 0.33), see Lu
et al. 2006.
For nonlinear analysis, the Friction Angle φ = 33° (given by CalTrans) and the peak shear
stress occurs at a shear strain γmax = 10% (at the 11.6 psi confinement). The parameter
γmax along with the shear modulus define the nonlinear soil stress-strain curve. Other
values of γmax should be explored in the future.
Axial Load
128
An axial load of 243 kips is applied at the pile head (free head connection).
The floating pile is modeled by beam-column elements (Mazzoni et al. 2006), and rigid
beam-column elements are used to model the pile size (diameter).
The axial load is applied at the pile head (ground level) in z (vertical) direction.
The above simulations were performed using OpenSeesPL (Lu et al. 2006).
Simulation Results
The pile vertical displacement and axial force profiles at the axial load of 243 kips are
shown in Figure G.2. The final deformed mesh is shown in Figure G.3. Figure G.4
displays the stress ratio contour fill.
129
(a) Isometric view
130
0
10
Depth (ft)
15
20
25
30
35
−0.18 −0.16 −0.14 −0.12 −0.1 −0.08 0 100 200
Vertical displacement (in) Axial force (kips)
Figure G.2: Pile profile response at the axial load of 243 kips.
131
a) plan view
b) Side view
Figure G.4: Stress ratio contour fill for the nonlinear analysis (red color shows yielded
soil elements).
132
Appendix H Moment-Curvature Analysis of
Circular Nonlinear RC Beam (Fiber Section)
Introduction
In this study, we compare with an OpenSees moment-curvature pushover analysis input
file (see appendix). A single circular reinforced concrete column is rigidly attached to the
soil mesh for that purpose. The soil domain is assumed rigid so as to simulate a cantilever
beam scenario.
The OpenSees input file is Example 9 listed in the OpenSees Example Manual
(http://opensees.berkeley.edu/OpenSees/manuals/ExamplesManual/HTML/). This
OpenSees example introduces the moment-curvature procedures for sections in 3D space.
The moment-curvature analysis of the section in this OpenSees example is by creating a
zero-length rotational-spring element. This section is subjected to a user-defined constant
axial load and to a linearly-increasing moment to a user-defined maximum curvature
(Mazzoni et al. 2006).
Fiber section is used to model the nonlinear behavior of the pile. The fiber section
properties are listed in Tables F.1-4. The schematic of the fiber section definition is also
shown in Figure F.1 (also see Figure F.2 for the input interface for fiber section in
OpenSeesPL):
Core Cover
Concrete Compressive Strength (ksi) -5.2 -4
Concrete Strain at Maximum Strength -0.002885 -0.003
Concrete Crushing Strength (ksi) -1.04 -0.8
Concrete Strain at Crushing Strength -0.0144 -0.01
Steel
Yield Strength (ksi) 66.8
Initial Elastic Tangent (ksi) 29000
Strain-hardening Ratio 0.01
133
Table H.3: Patch information for the pile circular cross section.
Core Cover
Number of Subdivisions (fibers) in the Circumferential Direction 8 8
Number of Subdivisions (fibers) in the Radial Direction 8 4
Internal Radius (in) 0 25
External Radius (in) 25 30
Table H.4: Layer information for the pile circular cross section.
Pile head (lateral) displacement of 0.69 in is applied in 25 equal steps. An axial load of
1800 kips is applied at the pile head (free head connection) during loading.
Simulation Results
The finite element mesh employed is shown in Figure F.3. As mentioned before, the soil
domain is rigid therefore the meshing of the soil domain is insignificant. 10
nonlinearBeamColumn elements are used for the pile.
The comparison of the moment-curvature curves is shown in Figure F.4. Both curves
match quite well.
Response profiles of the single pile are shown in Figure F.5. A shear load of 662 kips is
reached at the pile head longitudinal displacement of 0.69 in (Figure F.5 & F.6). The
maximum moment reached 6609 kip-ft, occurring at the ground surface (Figure F.6).
The moment-curvature curve of the single pile at the ground surface location is shown in
Figure F.7.
134
Cover patch
Core patch
External radius
for cover
Figure H.1: Schematic of the fiber section definition for the circular pile cross section.
135
Figure H.3: Finite element mesh employed in this study.
8000
7000
6000
Moment (kip-ft)
5000
4000
3000
OpenSeesPL
2000 OpenSees Example
1000
0
0.0000 0.0002 0.0004 0.0006 0.0008 0.0010 0.0012
Curvature (rad/in)
136
a) Longitudinal displacement
137
c) Longitudinal shear force
Figure H.6: Lateral (longitudinal) shear versus displacement at the pile head.
138
Figure H.7: Moment-curvature relation at the maximum moment location (ground
surface) in OpenSeesPL.
# --------------------------------------------------------------------------------------------------
# build a section
# Silvia Mazzoni & Frank McKenna, 2006
#
# SET UP ----------------------------------------------------------------------------
wipe; # clear memory of all past model definitions
model BasicBuilder -ndm 3 -ndf 6; # Define the model builder, ndm=#dimension, ndf=#dofs
set dataDir Data; # set up name of data directory -- simple
file mkdir $dataDir; # create data directory
source LibUnits.tcl; # define units
139
set IDconcCover 2; # material ID tag -- unconfined cover concrete
set IDreinf 3; # material ID tag -- reinforcement
# nominal concrete compressive strength
set fc [expr -4.0*$ksi]; # CONCRETE Compressive Strength, ksi (+Tension, -
Compression)
set Ec [expr 57*$ksi*sqrt(-$fc/$psi)]; # Concrete Elastic Modulus
# confined concrete
set Kfc 1.3; # ratio of confined to unconfined concrete strength
set fc1C [expr $Kfc*$fc]; # CONFINED concrete (mander model), maximum stress
set eps1C [expr 2.*$fc1C/$Ec]; # strain at maximum stress
set fc2C [expr 0.2*$fc1C]; # ultimate stress
set eps2C [expr 5*$eps1C]; # strain at ultimate stress
# unconfined concrete
set fc1U $fc; # UNCONFINED concrete (todeschini parabolic model), maximum
stress
set eps1U -0.003; # strain at maximum strength of unconfined concrete
set fc2U [expr 0.2*$fc1U]; # ultimate stress
set eps2U -0.01; # strain at ultimate stress
set lambda 0.1; # ratio between unloading slope at $eps2 and initial slope $Ec
# tensile-strength properties
set ftC [expr -0.14*$fc1C]; # tensile strength +tension
set ftU [expr -0.14*$fc1U]; # tensile strength +tension
set Ets [expr $ftU/0.002]; # tension softening stiffness
# -----------
set Fy [expr 66.8*$ksi]; # STEEL yield stress
set Es [expr 29000.*$ksi]; # modulus of steel
set Bs 0.01; # strain-hardening ratio
set R0 18; # control the transition from elastic to plastic branches
set cR1 0.925; # control the transition from elastic to plastic branches
set cR2 0.15; # control the transition from elastic to plastic branches
uniaxialMaterial Concrete01 $IDconcCore $fc1C $eps1C $fc2C $eps2C ; # build core concrete
(confined)
uniaxialMaterial Concrete01 $IDconcCover $fc1U $eps1U $fc2U $eps2U ; # build cover concrete
(unconfined)
uniaxialMaterial Steel01 $IDreinf $Fy $Es $Bs ; # build reinforcement material
puts "Ec = $Ec"
puts "uniaxialMaterial Concrete01 $IDconcCore $fc1C $eps1C $fc2C $eps2C ; # build core concrete
(confined)"
puts "uniaxialMaterial Concrete01 $IDconcCover $fc1U $eps1U $fc2U $eps2U ; # build cover concrete
(unconfined)"
puts "uniaxialMaterial Steel01 $IDreinf $Fy $Es $Bs ; # build reinforcement material"
#uniaxialMaterial Concrete02 $IDconcCore $fc1C $eps1C $fc2C $eps2C $lambda $ftC $Ets; # build
core concrete (confined)
#uniaxialMaterial Concrete02 $IDconcCover $fc1U $eps1U $fc2U $eps2U $lambda $ftU $Ets; # build
cover concrete (unconfined)
#uniaxialMaterial Steel02 $IDreinf $Fy $Es $Bs $R0 $cR1 $cR2; # build reinforcement
material
140
# with one layer of steel evenly distributed around the perimeter and a confined core.
# confined core.
# by: Michael H. Scott, 2003
#
#
# Notes
# The center of the reinforcing bars are placed at the inner radius
# The core concrete ends at the inner radius (same as reinforcing bars)
# The reinforcing bars are all the same size
# The center of the section is at (0,0) in the local axis system
# Zero degrees is along section y-axis
#
set ri 0.0; # inner radius of the section, only for hollow sections
set ro [expr $DSec/2]; # overall (outer) radius of the section
set nfCoreR 8; # number of radial divisions in the core (number of "rings")
set nfCoreT 8; # number of theta divisions in the core (number of "wedges")
set nfCoverR 4; # number of radial divisions in the cover
set nfCoverT 8; # number of theta divisions in the cover
source ex9.tcl
# define procedure
source MomentCurvature3D.tcl
141
Source code of file MomentCurvature3D. tcl:
proc MomentCurvature3D { secTag axialLoad maxK {numIncr 100} } {
##################################################
# A procedure for performing section analysis (only does
# moment-curvature, but can be easily modified to do any mode
# of section reponse.)
#
# MHS
# October 2000
# modified to improve convergence by Silvia Mazzoni, 2006
#
# Arguments
# secTag -- tag identifying section to be analyzed
# axialLoad -- axial load applied to section (negative is compression)
# maxK -- maximum curvature reached during analysis
# numIncr -- number of increments used to reach maxK (default 100)
#
# Sets up a recorder which writes moment-curvature results to file
# section$secTag.out ... the moment is in column 1, and curvature in column 2
# Define element
# tag ndI ndJ secTag
element zeroLengthSection 2001 1001 1002 $secTag
# Create recorder
recorder Node -file data/Mphi.out -time -node 1002 -dof 6 disp; # output moment (col 1) &
curvature (col 2)
142
pattern Plain 3002 "Linear" {
load 1002 0.0 0.0 0.0 0.0 0.0 1.0
}
143
algorithm $algorithmTypeStatic
}
if {$ok != 0} {
puts "Trying NewtonWithLineSearch .."
algorithm NewtonLineSearch 0.8
set ok [analyze 1]
algorithm $algorithmTypeStatic
}
if {$ok != 0} {; # stop if still fails
to converge
puts [format $fmt1 "PROBLEM" $IDctrlNode
$IDctrlDOF [nodeDisp $IDctrlNode $IDctrlDOF] $LunitTXT]
return -1
}; # end if
}; # end for
integrator DisplacementControl $IDctrlNode $IDctrlDOF $Dincr;
# bring back to original increment
}; # end if
}; # end while loop
}; # end if ok !0
# -----------------------------------------------------------------------------------------------------
global LunitTXT; # load time-unit text
if { [info exists LunitTXT] != 1} {set LunitTXT "Length"}; # set blank if it has not
been defined previously.
set fmt1 "%s Pushover analysis: CtrlNode %.3i, dof %.1i, Curv=%.4f /%s"; # format for
screen/file output of DONE/PROBLEM analysis
if {$ok != 0 } {
puts [format $fmt1 "PROBLEM" $IDctrlNode $IDctrlDOF [nodeDisp $IDctrlNode
$IDctrlDOF] $LunitTXT]
} else {
puts [format $fmt1 "DONE" $IDctrlNode $IDctrlDOF [nodeDisp $IDctrlNode
$IDctrlDOF] $LunitTXT]
}
144
References
M. Alizadeh and M. T. Davisson (1970). “Lateral Load Tests on Piles – Arkansas River
Project”, JSMFD, ASCE, Vol. 96, SM5, September, pp. 31-40
J. E. Bowles (1988). Foundation Analysis and Design, 4th Edition, McGraw-Hill Book
Co., New York, NY 10020.
Iwan, W. D. (1967). "On a class of models for the yielding behavior of continuous and
composite systems." J. Appl. Mech., ASME 34, 612-617.
Mazzoni, S., McKenna, F., and Fenves, G. L. (2006). Open system for earthquake
engineering simulation user manual, Pacific Earthquake Engineering Research Center,
University of California, Berkeley
(http://opensees.berkeley.edu/OpenSees/manuals/usermanual/).
Prevost, J. H. (1985). "A simple plasticity theory for frictional cohesionless soils." Soil
Dynamics and Earthquake Engineering, 4(1), 9-17.
Yang, Z. (2000). "Numerical modeling of earthquake site response including dilation and
liquefaction," PhD Thesis, Department of Civil Engineering and Engineering Mechanics,
Columbia University, New York, NY.
Yang, Z., Elgamal, A., and Parra, E. (2003). "A computational model for cyclic mobility
and associated shear deformation." Journal of Geotechnical and Geoenvironmental
Engineering, 129(12), 1119-1127.
145
OpenSeesPL-Related References
"Analyses and Modeling of Site Liquefaction Using Centrifuge Tests," (1996), E. Parra,
K. Adalier, A. -W. Elgamal, M. Zeghal, and A. Ragheb, Eleventh World Conference on
Earthquake Engineering, Acapulco, Mexico, June 23-28.
"Soil Dilation and Shear Deformations During Liquefaction," (1998a), A.-W. Elgamal, R.
Dobry, E. Parra, and Z. Yang, , Proc. 4th Intl. Conf. on Case Histories in Geotechnical
Engineering, S. Prakash, Ed., St. Louis, MO, March 8-15, pp1238-1259.
146
Mechanics, Columbia University, NY, New York.
"Computational Model for Cyclic Mobility and Associated Shear Deformation," (2003),
Z. Yang, A. Elgamal and E. Parra, J. Geotechnical and Geoenvironmental Engineering,
ASCE, 129(12), 1119-1127.
"A Web-based Platform for Live Internet Computation of Seismic Ground Response,"
(2004), Z. Yang, J. Lu, and A. Elgamal, Advances in Engineering Software, 35, 249-259.
147