I am trying to run this example script, but I can't seem to get it to run.
Here's the code:
####################################################################
# DISPL.-BASED FIBER ELEMENTS, REINFORCED-CONCRETE, REF: 'RCFRAME' #
####################################################################
#
#
# _________________________ _
# | | |
# | | | 12'
# | | |
# |____________|____________| _
# | | |
# | | |
# | | | 15'
# | | |
# | | | _
# === === ===
# | 24' | 24' |
#
#
set load 225.0
#
#
#
# comment out one of lines
#set displayMode "displayON"
set displayMode "displayOFF"
# CREATE THE MODEL BUILDERS
model BasicBuilder -ndm 2 -ndf 3
reliability
# DEFINE NODES
# tag X Y
node 1 0 0
node 2 0 180
node 3 0 324
node 4 288 0
node 5 288 180
node 6 288 324
node 7 576 0
node 8 576 180
node 9 576 324
# SINGLE POINT CONSTRAINTS
# node DX DY RZ
fix 1 1 1 1
fix 4 1 1 1
fix 7 1 1 1
# DEFINE MATERIALS
# Cover concrete
# tag -f'c -epsco -f'cu -epscu
uniaxialMaterial Concrete01 1 -4 -0.002 0.0 -0.006
# Core concrete
uniaxialMaterial Concrete01 2 -5.2 -0.005 -4.70 -0.02
# Steel model
# tag fy E b
uniaxialMaterial Steel01 3 60 30000 0.02
# DEFINE SECTIONS
# Interior column section
section fiberSec 1 {
# mat nfIJ nfJK yI zI yJ zJ yK zK yL zL
patch quadr 2 1 12 -11.5 10 -11.5 -10 11.5 -10 11.5 10
patch quadr 1 1 14 -13.5 -10 -13.5 -12 13.5 -12 13.5 -10
patch quadr 1 1 14 -13.5 12 -13.5 10 13.5 10 13.5 12
patch quadr 1 1 2 -13.5 10 -13.5 -10 -11.5 -10 -11.5 10
patch quadr 1 1 2 11.5 10 11.5 -10 13.5 -10 13.5 10
# mat nBars area yI zI yF zF
layer straight 3 6 1.56 -10.5 9 -10.5 -9
layer straight 3 6 1.56 10.5 9 10.5 -9
}
# Exterior column section
section fiberSec 2 {
patch quadr 2 1 10 -10 10 -10 -10 10 -10 10 10
patch quadr 1 1 12 -12 -10 -12 -12 12 -12 12 -10
patch quadr 1 1 12 -12 12 -12 10 12 10 12 12
patch quadr 1 1 2 -12 10 -12 -10 -10 -10 -10 10
patch quadr 1 1 2 10 10 10 -10 12 -10 12 10
layer straight 3 6 0.79 -9 9 -9 -9
layer straight 3 6 0.79 9 9 9 -9
}
# Girder section
section fiberSec 3 {
patch quadr 1 1 12 -12 9 -12 -9 12 -9 12 9
layer straight 3 4 1.00 -9 9 -9 -9
layer straight 3 4 1.00 9 9 9 -9
}
# NUMBER OF INTEGRATION POINTS
set nP 4
# GEOMETRIC TRANSFORMATION
geomTransf Linear 1
# DEFINE DISPLACEMENT BEAM-COLUMN ELEMENTS
# Columns
# tag ndI ndJ nPts secID transf
element dispBeamColumn 1 1 2 $nP 2 1
element dispBeamColumn 2 2 3 $nP 2 1
element dispBeamColumn 3 4 5 $nP 1 1
element dispBeamColumn 4 5 6 $nP 1 1
element dispBeamColumn 5 7 8 $nP 2 1
element dispBeamColumn 6 8 9 $nP 2 1
# Beams
element dispBeamColumn 7 2 5 $nP 3 1
element dispBeamColumn 8 5 8 $nP 3 1
element dispBeamColumn 9 3 6 $nP 3 1
element dispBeamColumn 10 6 9 $nP 3 1
# CONSTANT GRAVITY LOADS
pattern Plain 1 Linear {
# node FX FY MZ
load 2 0.0 -192.0 0.0 -const
load 5 0.0 -384.0 0.0 -const
load 8 0.0 -192.0 0.0 -const
load 3 0.0 -96.0 0.0 -const
load 6 0.0 -192.0 0.0 -const
load 9 0.0 -96.0 0.0 -const
}
pattern Plain 2 Linear {
# node FX FY MZ
load 2 $load 0.0 0.0
load 3 $load 0.0 0.0
}
# FINITE ELEMENT ANALYSIS MODEL
integrator LoadControl 1 3 1 1
test NormDispIncr 1.0e-06 10 0
algorithm Newton
numberer RCM
constraints Transformation
#system SparseGeneral -piv
system ProfileSPD
sensitivityIntegrator -static
sensitivityAlgorithm -computeAtEachStep
analysis Static
# RECORDERS
recorder Node -file Drift.txt -time -node 3 -dof 1 disp
######################################################################
# THIS IS AN AUXILIARY FILE FOR THE STATIC RC-FRAME EXAMPLES, #
# IT SIMPLY CONTAINS THE COMMANDS THAT SHOULD BE EXECUTED FOR EACH #
# LIMIT-STATE FUNCTION EVALUATION. #
######################################################################
integrator LoadControl 0.05 3 0.05 0.05
analyze 16
integrator LoadControl 0.02 3 0.02 0.02
analyze 10
-------------------------------------------------
what's wrong with this file? Any help would be appreciated!! Thanks!!
