9/11/2008

Geotechnical Elements and

Models in OpenSees

Pedro Arduino
University of Washington, Seattle

OpenSees Days 2008,

OpenSees User Workshop, Monday Sept 8, 2008

Type of Geotechnical Problems that

can be solved using OpenSees

Static Problems
 Deformation analyses (1D, 2D, or 3D)
 Consolidation problems
 Soil-structure interaction problems

Shallow foundations (e.g. bearing capacity, deformation)

Pile foundations (e.g. vertical and lateral capacity)

Dynamic (earthquake problems)

 Free-field analysis
 Liquefaction induced problems
 Soil structure interaction problems (e.g. response of
pile foundations, bridge bents, or complete structures embedded in
soils to earthquake excitations)

1
 9/11/2008

What do we need??

Solid elements to characterize the soil domain
(continuum).

Appropriate boundary conditions to accurately
represent the soil domain boundaries.

Robust constitutive models to characterize the
soil stress-strain response under monotonic and
cyclic loading conditions

Interface elements to capture the interaction
between the soil and adjacent structures.

Everything else you are learning in this
workshop (i.e., how to create beam elements,
apply loads, record results, perform the analysis,
etc.

Outline

Finite Elements (for solids)
 Single-phase
 Multi-phase (coupled) finite elements
 Zero length element

Material Models
 Elastic
 Elasto-plastic Continuum Models
 Elasto-plastic Uniaxial models

Boundary Conditions
 Equal DOF
 Absorbent boundaries

2
 9/11/2008

Finite Elements (solids)

Single-phase formulation
 To capture the response of dry soils need
one single phase
 Phase 1 – soil skeleton

Multi-phase formulation
 To capture the response of saturated soils 
need two phases
 Phase 1  Soil skeleton
 Phase 2  pore water

Zero-Length element
 To capture interface response between solid
and beam elements

Single Phase Formulations

Small deformation solid elements

 2-D quadrilateral elements (4, 9 nodes)
 3-D solid elements, brick (8, 20 nodes)

Large deformation (total
Lagrangian) solid elements, bricks
(20 nodes)
n4 n3
n8

n5 n7
n6
n4
n3
n1

n2
n1 n2
quad (4 node) stdBrick (8 node)

3
 9/11/2008

quad element definition

n4 n3

n1 n2

quad (4 node)

element quad $eleTag $n1 $n2 $n3 $n4 $thick $type $matTag
<$press $rho $b1 $b2>

Must define first all the required arguments. In particular:

Nodes $n1, $n2, $n3, $n4 and
Material type $matTag
The arguments in <…> are optional

Multi-Phase Formulations

Fully coupled u-p elements (2D & 3D)

Fully coupled u-p-U elements (3D) for
small deformations
n4 n3 n7 n6 n5

n9
n8 n4

n1 n2 n1 n2 n3
quadUP 9_4_quadUP

Degrees of Freedom (DOFs) are:

 u  solid displacement, on
 P  pore fluid pressures, on
 U  pore fluid displacements, on

4
 9/11/2008

quadUP element definition

n4 n3

n1 n2
quadUP

element quadUP $eleTag $n1 $n2 $n3 $n4 $thick $type $matTag
$bulk $fmass $hPerm $vPerm <$b1 $b2 $t>

$bulk  combined undrained bulk modulus Bc=Bf/n

$fmass  fluid mass density
$hperm & $vperm  horiz. And vert. permeability

zerolength element

Connects two points at the same

coordinate
n4 n3
n6

n5
n1 n2
solid element beam element

zero-length element

element zeroLength $eleTag $n1 $n2 –mat $matTag1 $matTag2 …

-dir $dir1 $dir2 … <-orient $x1 $x2 $x3 $yp1 $yp2 $yp3>

5
 9/11/2008

Material Models

Linear Elastic Material model (nDMaterial)

 To characterize the response of the soil (or other
continuum) in its elastic regime

Elasto-Plastic Material models (nDMaterial)
 To characterize the nonlinear stress-strain
response of soils

Elasto-plastic Uniaxial models
 To characterize the interface response between
soil and structural elements (uniaxialMaterial).

nDMaterial
Elastic

Small deformation elasticity

 Linear isotropic
 Nonlinear isotropic
 Cross anisotropic

Elastic Isotropic Material

nDMaterial ElasticIsotropic $matTag $E $v

6
 9/11/2008

nDMaterial
Elasto-Plastic (Small Deformations)

J2-Plasticity Material (von Mises)

Drucker-Prager Material

Template Elasto-Plastic Material (UC
Davis)

Cam-Clay Material (Berkeley)

MutiYield Materials (San Diego)

FluidSolidPorous Material(SanDiego)

nDMaterial
J2Plasticity

von-Mises type
nDMaterial J2Plasticity $matTag $K $G $sig0 $sigInf $delta $H

σd
σd-inf

σd-0

Von-Mises Yield Surface

εd
Stress-strain curve

7
 9/11/2008

nDMaterial
Template Elasto-Plastic Material

Versatile tool to generate multiple types

of elasto-plastic materials by combining
yield surfaces, plastic potentials and
evolution laws

Developed by Boris Jeremic at UC Davis

http://sokocalo.engr.ucdavis.edu/~jeremic

nDMaterial
Template Elasto-Plastic Material
nDMaterial Template3Dep $matTag $ElmatTag
–YS $ys –PS $ps –EPS $eps <-ELS1 $el> <-ELT1 $et>

-YS  Yield surfaces (von Mises, Drucker Prager, Mohr-

Coulomb, Camclay)
-PS  Plastic potentials (von Mises, Drucker-Prager Mohr-
Coulomb, Camclay, Leon)
-EPS  Initial state of stress
-ELS1  Scalar evolution laws for isotropic hardening (linear,
nonlinear Camclay)
-ELT1  Tensorial evolution laws for Kinematic hardening (linear,
nonlinear Armstrong-Frederick)
Von-Mises Drucker Prager Mohr-Coulomb Camclay

8
 9/11/2008

nDMaterial
MultiYield Materials

Material models based on Multiyield

Plasticity (Mroz et al., Prevost et al.)

Two types
 Pressure Independent Multi-yield
 Pressure Dependent Multiy-ield (captures well the
response of liquefiable soils)

 Material to couple solid &fluid phases

Developed by Elgamal et al. at UCSD

http://cyclic.ucsd.edu/opensees/

nDMaterial
PressureDependentMultiYield
nDMaterial PressureDependMultiYield $matTag $nd $rho
$refShearModul $refBulkModul $frictionAng $peakShearStra
$refPress $pressDependCoe $PTAng
$contrac $dilat1 $dilat2, $liquefac1 $liquefac2 $liquefac3
<$noYieldSurf=20 <$r1 $Gs1 …>
$e=0.6 $cs1=0.9 $cs2=0.02 $cs3=0.7 $pa=101>

15 parameters!!??

9
 9/11/2008

nDMaterial
PressureDependentMultiYield
nDMaterial PressureDependMultiYield $matTag $nd $rho
$refShearModul $refBulkModul $frictionAng $peakShearStra
$refPress $pressDependCoe $PTAng
$contrac $dilat1 $dilat2, $liquefac1 $liquefac2 $liquefac3
<$noYieldSurf=20 <$r1 $Gs1 …>
$e=0.6 $cs1=0.9 $cs2=0.02 $cs3=0.7 $pa=101>

nDMaterial
PressureDependentMultiYield02
nDMaterial PressureDependMultiYield02 $matTag $nd $rho
$refBulkModul $frictionAng $peakShearStra $refPress
$pressDepenCoe $PTAng
$contrac1 $contrac3 $dilat1 $dilat3
<$noYieldSurf=20 <$r1 $Gs1 …>
$contrac2=5.0 $dilat2=3.0 $liquefac1=1.0 $liquefac2=0.0
$e=0.6 $cs1=0.9 $cs2=0.02 $cs3=0.7 $pa=101>

10
 9/11/2008

nDMaterial
PressureIndependentMultiYield
nDMaterial PressureIndependMultiYield $matTag $nd $rho
$refShearModul $refBulkModul $cohesi $peakShearStra
$frictionAng $refPress=101 $pressDependCoe=0.
<$noYieldSurf=20 <$r1 $Gs1 …>>

nDMaterial
PressureIndependentMultiYield
nDMaterial PressureIndependMultiYield $matTag $nd $rho
$refShearModul $refBulkModul $cohesi $peakShearStra
$frictionAng $refPress=101 $pressDependCoe=0.
<$noYieldSurf=20 <$r1 $Gs1 …>>

11
 9/11/2008

nDMaterial
FluidSolidPorousMaterial

Couples the response of two

phases:fluid and solid – developed to
simulate the response of saturated porous media

nDMaterial FluidSolidPorousMaterial $matTag $nd

$soilMatTag $combinedBulkModul

$soilMatTag  the tag of previously defined material

$combinedBulkModul  combined undrained bulk modulus,
Bc=Bf/n

Additional commands for multiyield

materials

Help perform stage analysis

updateMaterialStage –material $matTag –stage $sNum

$MatTag  the tag of previously defined material

$sNum  (0 - elastic, 1-plastic, 2 – linear elastic constant f(σ
σ3) )

updateParameter –material $matTag –refG $newVal

$MatTag  the tag of previously defined material

$sNewVal  new parameter value

12
 9/11/2008

Elasto-plastic Uniaxial models

To capture interface response between
solid (soil) and beam elements (pile)

PyTzQ Uniaxial Materials

•PySimple1
•TzSimple1
•QzSimple1

•PyLiq1
•TzLiq1

uniaxialMaterial
PySimple1
uniaxialMaterial PySimple1 matTag $soilType $pult $Y50 $Cd
<$c>

$soilType  =1 Matlock (clay), =2 API (sand)

Cd=1.0
$pult  ultimate capacity of p-y material
$Y50  displ. @ 50% of pult
Cd  drag resistance (=1 no gap, <1 gap)
$c  viscous damping

p
pult (Reese 1974)
Cd=0.3

y50 (API 1993) y

13
 9/11/2008

uniaxialMaterial
TzSimple1 & QzSimple1
uniaxialMaterial TzSimple1 matTag $tzType $tult $z50 <$c>

$tzType =1 Reese & O’Neill (clay), =2 Mosher (sand)

$tult  ultimate capacity of t-z material
$z50  displ. @ 50% of tult
$c  viscous damping

uniaxialMaterial QzSimple1 matTag $qzType $qult $z50

<$suction $c>

$qzType =1 Reese & O’Neill (clay), =2 Vijayvergiya (sand)

$qult = ultimate capacity of q-z material
$z50 = displ. @ 50% of qult
$suction  uplift resistance = suction*qult
$c viscous damping

uniaxialMaterial
PyLiq1
uniaxialMaterial PyLiq1 $matTag $soilType $pult $Y50 $Cd $c
$pRes $solidElem1 $solidElem2

$soilType  =1 Matlock (clay), =2 API (sand)

$pult  ultimate capacity of p-y material
$Y50  displ. @ 50% of pult
Cd  drag resistance (=1 no gap, <1 gap)
$c  viscous damping
$pRes  residual (minimum) p-y resistance as ru=1.0
$solidElem1 & $solidElem2  solid elements from which PyLiq1
will obtain effective stresses and pore pressures

14
 9/11/2008

uniaxialMaterial
PyLiq1

Boundary Conditions

EqualDof
equalDOF $rNodeTag $cNodeTaga $dof1 $dof2 …
$rNodeTag  master node
$cNodeTag  slave node
$dof1 $dof2 …  constrained dof’s Same lateral
deformation

15
 9/11/2008

Absorbent Boundaries
Lysmer (1969)

Quad Element

1. set DampP 755

2. set DampN 1216
3. uniaxialMaterial Elastic 1 0 $DampP
4. uniaxialMaterial Elastic 2 0 $DampN
5. node 1 16.0 0.0
6. node 2 16.0 0.0
7. element zeroLength 1 1 2 -mat 1 2
-dir 1 2 –orient 1 –2 0 2 1 0

C N = a ρ VS
C P = b ρV P
zeroLength Element &
uniaxial material

Many more capabilities currently

under development!!

16