# Define SECTIONS -------------------------------------------------------------
set SectionType FiberSection ; # options: Elastic FiberSection
#set SectionType Elastic ; # options: Elastic FiberSection
# define section tags:
set ColSecTag 1
set BeamSecTag 2
set BeamSecTagCen 3
set GirdSecTag 4
set ColSecTagFiber 5
set BeamSecTagFiber 6
set GirdSecTagFiber 7
set SecTagTorsion 70
# Section Properties:
set HCol 0.5; # square-Column width
set BCol $HCol
set HBeam 0.3; # Beam depth -- perpendicular to bending axis @sides
set BBeam $HBeam; # Beam width -- parallel to bending axis @sides
set HBeam2 0.4; # Beam Depth ---- Frame 2/Axis 2-2
set BBeam2 0.3; # Beam width ---- Frame 2/Axis 2-2
set HGird 0.6; # Girder depth -- perpendicular to bending axis
set BGird 0.3; # Girder width -- parallel to bending axis
if {$SectionType == "Elastic"} {
# material properties:
set fc 4000*$psi; # concrete nominal compressive strength
set Ec [expr 57*$ksi*pow($fc/$psi,0.5)]; # concrete Young's Modulus
set nu 0.2; # Poisson's ratio
set Gc [expr $Ec/2./[expr 1+$nu]]; # Torsional stiffness Modulus
set J $Ubig; # set large torsional stiffness
# column section properties:
set AgCol [expr $HCol*$BCol]; # rectuangular-Column cross-sectional
area
set IzCol [expr 0.5*1./12*$BCol*pow($HCol,3)]; # about-local-z Rect-Column
gross moment of inertial
set IyCol [expr 0.5*1./12*$HCol*pow($BCol,3)]; # about-local-z Rect-Column
gross moment of inertial
# beam sections:
set AgBeam [expr $HBeam*$BBeam]; # rectuangular-Beam cross-sectional
area
set IzBeam [expr 0.5*1./12*$BBeam*pow($HBeam,3)]; # about-local-z Rect-
Beam cracked moment of inertial
set IyBeam [expr 0.5*1./12*$HBeam*pow($BBeam,3)]; # about-local-y Rect-
Beam cracked moment of inertial
# girder sections:
set AgGird [expr $HGird*$BGird]; # rectuangular-Girder cross-
sectional area
set IzGird [expr 0.5*1./12*$BGird*pow($HGird,3)]; # about-local-z Rect-
Girder cracked moment of inertial
set IyGird [expr 0.5*1./12*$HGird*pow($BGird,3)]; # about-local-y Rect-
Girder cracked moment of inertial
section Elastic $ColSecTag $Ec $AgCol $IzCol $IyCol $Gc $J
section Elastic $BeamSecTag $Ec $AgBeam $IzBeam $IyBeam $Gc $J
section Elastic $GirdSecTag $Ec $AgGird $IzGird $IyGird $Gc $J
set IDconcCore 1; # material numbers for recorder (this
stressstrain recorder will be blank, as this is an elastic section)
set IDSteel 2; # material numbers for recorder (this
stressstrain recorder will be blank, as this is an elastic section)
} elseif {$SectionType == "FiberSection"} {
# MATERIAL parameters
source LibMaterialsRC.tcl; # define library of Reinforced-concrete
Materials
# FIBER SECTION properties
# Column section geometry:
set cover [expr 2.5*$in]; # rectangular-RC-Column cover
set numBarsTopCol 8; # number of longitudinal-reinforcement bars on
top layer
set numBarsBotCol 8; # number of longitudinal-reinforcement bars on
bottom layer
set numBarsIntCol 6; # TOTAL number of reinforcing bars on the
intermediate layers
set barAreaTopCol [expr 1.*$in*$in]; # longitudinal-reinforcement bar
area
set barAreaBotCol [expr 1.*$in*$in]; # longitudinal-reinforcement bar
area
set barAreaIntCol [expr 1.*$in*$in]; # longitudinal-reinforcement bar
area
set numBarsTopBeam 6; # number of longitudinal-reinforcement bars on
top layer
set numBarsBotBeam 6; # number of longitudinal-reinforcement bars on
bottom layer
set numBarsIntBeam 2; # TOTAL number of reinforcing bars on the
intermediate layers
set barAreaTopBeam [expr 1.*$in*$in]; # longitudinal-reinforcement bar
area
set barAreaBotBeam [expr 1.*$in*$in]; # longitudinal-reinforcement bar
area
set barAreaIntBeam [expr 1.*$in*$in]; # longitudinal-reinforcement bar
area
set numBarsTopGird 6; # number of longitudinal-reinforcement bars on
top layer
set numBarsBotGird 6; # number of longitudinal-reinforcement bars on
bottom layer
set numBarsIntGird 2; # TOTAL number of reinforcing bars on the
intermediate layers
set barAreaTopGird [expr 1.*$in*$in]; # longitudinal-reinforcement bar
area
set barAreaBotGird [expr 1.*$in*$in]; # longitudinal-reinforcement bar
area
set barAreaIntGird [expr 1.*$in*$in]; # longitudinal-reinforcement bar
area
set nfCoreY 20; # number of fibers in the core patch in the y
direction
set nfCoreZ 20; # number of fibers in the core patch in the z
direction
set nfCoverY 20; # number of fibers in the cover patches with long
sides in the y direction
set nfCoverZ 20; # number of fibers in the cover patches with long
sides in the z direction
# rectangular section with one layer of steel evenly distributed around the
perimeter and a confined core.
BuildRCrectSection $ColSecTagFiber $HCol $BCol $cover $cover $IDconcCore
$IDconcCover $IDSteel $numBarsTopCol $barAreaTopCol $numBarsBotCol $barAreaBotCol
$numBarsIntCol $barAreaIntCol $nfCoreY $nfCoreZ $nfCoverY $nfCoverZ
BuildRCrectSection $BeamSecTagFiber $HBeam $BBeam $cover $cover $IDconcCore
$IDconcCover $IDSteel $numBarsTopBeam $barAreaTopBeam $numBarsBotBeam
$barAreaBotBeam $numBarsIntBeam $barAreaIntBeam $nfCoreY $nfCoreZ $nfCoverY
$nfCoverZ
BuildRCrectSection $GirdSecTagFiber $HGird $BGird $cover $cover $IDconcCore
$IDconcCover $IDSteel $numBarsTopGird $barAreaTopGird $numBarsBotGird
$barAreaBotGird $numBarsIntGird $barAreaIntGird $nfCoreY $nfCoreZ $nfCoverY
$nfCoverZ
# assign torsional Stiffness for 3D Model
uniaxialMaterial Elastic $SecTagTorsion $Ubig
section Aggregator $ColSecTag $SecTagTorsion T -section $ColSecTagFiber
section Aggregator $BeamSecTag $SecTagTorsion T -section $BeamSecTagFiber
section Aggregator $GirdSecTag $SecTagTorsion T -section $GirdSecTagFiber
} else {
puts "No section has been defined"
return -1
}
set GammaConcrete [expr 150*$pcf]; # Reinforced-Concrete weight
density (weight per volume)
set QdlCol [expr $GammaConcrete*$HCol*$BCol]; # self weight of Column, weight per
length
set QBeam [expr $GammaConcrete*$HBeam*$BBeam]; # self weight of Beam, weight per
length
set QGird [expr $GammaConcrete*$HGird*$BGird]; # self weight of Gird, weight per
length
