A5 StructModeling 2010
A5 StructModeling 2010
With Examples
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
ModelBuilder Objects
model Command
node Command
mass Command
constraints objects
uniaxialMaterial Command
nDMaterial Command
section Command
Geometric Transformation Command
element Command
block Command
region Command
Time Series
pattern Command
Model Command
This command is used to construct the BasicBuilder object.
$ndf
ModelBuilder Objects
model Command
node Command
mass Command
constraints objects
uniaxialMaterial Command
nDMaterial Command
section Command
Geometric Transformation Command
element Command
block Command
region Command
Time Series
pattern Command
Nodal Coordinates
sample command
$coords
$MassValues
unconstrained
constrained
# Define nodes;
2.
node 1
0.0
0.0
3.
node 2
100.
0.0
0.0
4.
node 3
0.0
$Lcol
0.0
5.
node 4
$Lbeam
$Lcol
0.0
6.
# Boundary conditions;
7.
fix
1 1 1 1 1 1;
8.
fix
1 1 1 1 1 1
9.
fix
0 0 1 1 1 0
10.
fix
0 0 1 1 1 0
11.
12.
13.
14.
--1--
--2--
boundary conditions
3-----------------------4
ModelBuilder Objects
model Command
node Command
mass Command
constraints objects
uniaxialMaterial Command
nDMaterial Command
section Command
Geometric Transformation Command
element Command
block Command
region Command
Time Series
pattern Command
uniaxialMaterials:
Contributed Materials
BARSLIP Material
Bond_SP01
Concrete04 Material Popovics
Concrete07 Chang & Manders
Fatigue Material
Hyperbolic Gap Material
Limit State Material
PINCHING4 Material
PyTzQz Uniaxial Materials
Reinforcing Steel Material
SelfCentering Material
materials
$E
tangent
$eta
$fpc
compressive strength*
$epsc0
$fpcu
crushing strength*
$epsU
strain at crushing strength*
*NOTE: Compressive concrete parameters should be input as negative values.
Silvia Mazzoni, OpenSees Days 2010
tcl if statement
if {logical statement} {
.series of commands..
}
set a 7
if {$a==7} {
puts $a
}
materials
concrete
2.
if {$ConcreteMaterialType == elastic} {
3.
# options: "elastic","inelastic"
4.
5.
} else {;
# $ConcreteMaterialType == inelastic
6.
7.
8.
9.
# Core
# Cover
$s1p
$e1p
stress and strain (or force & deformation) at first point of the envelope
in the positive direction
$s2p
$e2p
stress and strain (or force & deformation) at second point of the
envelope in the positive direction
$s3p
$e3p
stress and strain (or force & deformation) at third point of the envelope
in the positive direction (optional)
$s1n
$e1n
stress and strain (or force & deformation) at first point of the envelope
in the negative direction*
$s2n
$e2n
stress and strain (or force & deformation) at second point of the
envelope in the negative direction*
$s3n
$e3n
stress and strain (or force & deformation) at third point of the envelope
in the negative direction (optional)*
$pinchX
$pinchY
$damage1
$damage2
$beta
materials
reinforcing steel
ModelBuilder Objects
model Command
node Command
mass Command
constraints objects
uniaxialMaterial Command
nDMaterial Command Presented in Geotech
section Command
Geometric Transformation Command
element Command
block Command
region Command
Time Series
pattern Command
OpenSees Sections:
beam-column element:
Elastic Section
Uniaxial Section
Fiber Section
Section Aggregator
Uniaxial Section
user-defined moment-curvature relationship (use
unaxialMaterial)
uncoupled P-M and anything else
Fiber Section
user-defined section geometry/materials via fibers
coupled P-M interaction
coupled bi-directional response
Section Aggregator
combine all uncoupled responses (e.g., Uniaxial flexure +
Uniaxial Axial, Fiber flexure/axial + shear)
Silvia Mazzoni, OpenSees Days 2010
core patch
external
radius
internal
radius
steel layer
$yLoc
$zLoc
$A
area of fiber
$matTag
$numSubdivCirc
$numSubdivRad
irc=4
$numSubdivC
number of subdivisions (fibers) in the radial
direction.
$intRad
internal radius
$extRad
external radius
$startAng
$endAng
$e
xt
Ra
d
$n
um
Su
bd
ivR
ad
=3
$endAng
$i
nt
Ra
d
$startAng
($yCenter, $zCenter)
patch quad $matTag $numSubdivIJ $numSubdivJK $yI $zI $yJ $zJ $yK $zK $yL $zL
Silvia Mazzoni, OpenSees Days 2010
$numBar
$areaBar
$yCenter
$zCenter
$radius
$startAng
$endAng
rs
Ba
um
$n
$endAng
$startAng
$rad
ius
($yCenter, $zCenter)
tcl procedure
used for repeated series of commands on a set of input variables
proc multiply {a b} {
set c [expr $a*$b]
return $c
}
set a 3; set b 5
set result [multiply $a $b]
Silvia Mazzoni, OpenSees Days 2010
CONCRETE
# Core radius
STEEL
if {$Nbars<= 0} { return }
}
}
core patch
cover
CONCRETE
cover patch
Ri
call:
Ro
360o/Nbars
Ab
STEEL
RCcircSection $id $Ri $Ro $cover $coreID $coverID $steelID $Nbars $Ab $nfCoreR $nfCoreT
$nfCoverR $nfCoverT
UniaxialSection
Use uniaxialMaterial to define
section moment-curvature
relationship
Hysteretic Material
section aggregator
section Aggregator $secTag $matTag1 $string1 $matTag2 $string2 ....... <-section $sectionTag>
$secTag
the force-deformation quantities corresponding to each section object. One of the following strings
is used:
P
Axial force-deformation
Mz
Vy
My
Vz
T
<-section $sectionTag>
Vy
Vy
M
Mz, Vy
Mz
P, Mz, Vy
P, Mz
ModelBuilder Objects
model Command
node Command
mass Command
constraints objects
uniaxialMaterial Command
nDMaterial Command
section Command
Geometric Transformation Command
element Command
block Command
region Command
Time Series
pattern Command
geometric transformation
performs a linear geometric transformation of beam
stiffness and resisting force from the basic system to the
global-coordinate system
geomTransf $Type $transfTag $vecxzX $vecxzY $vecxzZ
<-jntOffset $dXi $dYi $dZi $dXj $dYj $dZj>
$Type
$transfTag
$vecxzX
$vecxzY
$vecxzZ
ONLY IN 3D. X, Y, and Z components of vecxz, the vector used to define the local
x-z plane of the local-coordinate system. The local y-axis is defined by taking the
cross product of the x-axis and the vecxz vector.
These components are specified in the global-coordinate system X,Y,Z and define a
vector that is in a plane parallel to the x-z plane of the local-coordinate system.
These items need to be specified for the three-dimensional problem.
$dXi $dYi
$dZi
joint offset values -- absolute offsets specified with respect to the global
coordinate system for element-end node i (the number of arguments depends on
the dimensions of the current model) (optional)
$dXj $dYj
$dZj
joint offset values -- absolute offsets specified with respect to the global
coordinate system for element-end node j (the number of arguments depends on
the dimensions of the current model) (optional)
e
lan
zp
x
al
loc
node i
element orientation:
element 2
Y
X
y
x
z
y
1
t
n
e
m
le
e
ModelBuilder Objects
model Command
node Command
mass Command
constraints objects
uniaxialMaterial Command
nDMaterial Command
section Command
Geometric Transformation Command
element Command
block Command
region Command
Time Series
pattern Command
elements
Truss Element
Corotational Truss Element
Elastic Beam Column Element
NonLinear Beam-Column Elements
Zero-Length Elements
Quadrilateral Elements
Brick Elements
FourNodeQuadUP Element
BeamColumnJoint Element
3D:
element elasticBeamColumn $eleTag $iNode $jNode
$A $E $G $J $Iy $Iz $transfTag
end nodes
$numIntgrPts
$secTag
$transfTag
$massDens
element mass density (per unit length), from which a lumped-mass matrix is formed
(optional, default=0.0)
$maxIters
$tol
element node
element node
element
force-deformation
(stress-strain)
section
elements
set ColumnType inelastic;
Tcl procedure
source RCcircSection.tcl;
RCcircSection
$IDcolFlex $riCol $roCol $cover $IDcore $IDcover $IDsteel $NbCol $AbCol $nfCoreR $nfCoreT $nfCoverR $nfCoverT
# Define torsion
$IDcolTors
T -section $IDcolFlex;
$Acol
$Ec
$G
element elasticBeamColumn 2
2 4
$Acol
$Ec
$G
} else {;
# $ColumnType == inelastic
element nonlinearBeamColumn 1
element nonlinearBeamColumn 2
geomTransf Linear
$IDbeamTrans
element elasticBeamColumn
3 3 4
BEAM
element elasticBeamColumn 1
COLUMN
if {$ColumnType == elastic} {
beamWithHinges Element
2D:
3D:
inelastic-section
Silvia Mazzoni, OpenSees Days 2010
beamWithHinges Element
$eleTag
$iNode
$jNode
end nodes
$secTagI
$HingeLengthI
$secTagJ
$HingeLengthJ
$E
Young's Modulus
$A
$Iz
$Iy
$G
Shear Modulus
$J
$transfTag
$massDens
element mass density (per unit length), from which a lumped-mass matrix is
formed (optional, default=0.0)
$maxIters
ZeroLength Element
DOF spring
Use uniaxialMaterial to define
force-deformation response.
Translational dof: force-deformation
Rotational dof: moment-rotation
e.g. Hysteretic Material
1M applications
ModelBuilder Objects
model Command
node Command
mass Command
constraints objects
uniaxialMaterial Command
nDMaterial Command
section Command
Geometric Transformation Command
element Command
block Command
region Command
Time Series
pattern Command
TimeSeriesType
arguments
load ...
sp ...
eleLoad ...
$nodeTag
$LoadValues
$nodeTag
$DOFtag
$DOFvalue
Structural Example
Reinforced-Concrete Frame:
Building the Model
Analysis
Model Building
problem statement
Reinforced-Concrete Portal Frame
use kip, inch and sec as basic units
GEOMETRY
4000kip
Weight
B
A
Y
Z
A
36
LBeam
42
LCol
5 BBeam
BCol
5
HBeam
HCol
section A-A
section B-B
Ex4.Portal2D.build.InelasticFiberSection.tcl
LibUnits.tcl
element
node
node
element
element
node
Y
X
node
Materials
Section Command
concrete
steel
l
stee
material stress-strain
Silvia Mazzoni, OpenSees Days 2010
HSec
concrete
BSec
HCol
CoverCol
BCol
element
node
Silvia Mazzoni, OpenSees Days 2010
node
element
element
node
Y
X
node
Recorders
Gravity Load
node
element
element
node
3
2
element
node
Y
X
node
Gravity Analysis
Run OpenSees
source Ex4.Portal2D.build.InelasticSection.tcl
Structural Example
Reinforced-Concrete Frame:
Lateral-Load Analysis
Analysis
Model Building
LibAnalysisStaticParameters.tcl
LibAnalysisDynamicParameters.tcl
Questions,
or statements:
The OpenSees Community Forum:
http://opensees.berkeley.edu/community/in
dex.php
which can be accessed from:
http://opensees.berkeley.edu