|
|
|
# load pattern test
# -----------------
# Jian Zhao, Iowa State University
# June 20, 2004
# -----------------
#----------Create Model (2-D 3dof problem-plane truss)
model basic -ndm 2 -ndf 3
# Create nodes
# tag X Y
node 1 0.0 0.0
node 2 1.0 0.0
node 3 1.0 1.0
# Fix supports at node 1
# Tag DX DY RZ
fix 1 1 1 1
fix 2 0 1 1
fix 3 0 0 1
# Define faked material
# unconfined CONCRETE tag f'c ec0 f'cu ecu
uniaxialMaterial Concrete01 1 -4.5 -0.002 0.0 -0.006
# confined CONCRETE tag f'c ec0 f'cu ecu
uniaxialMaterial Concrete01 2 -6.0 -0.004 -5.0 -0.014
# unconfined CONCRETE tag f'c ec0 f'cu ecu ratio ft et0 ft0 beta etu
uniaxialMaterial Concrete03 3 -6.0 -0.004 -5.0 -0.014 0.5 0.45 0.0015 0.1 2.5 0.0025
# STEEL02 tag $Fy $E $b $R0 $cR1 $cR2 $a1 $a2 $a3 $a4
uniaxialMaterial Steel02 4 68.1 29000 0.010 22 0.925 0.15 0.04 1.0 0.04 1.0
#uniaxialMaterial StrPen01 tag $sy $fy $su $fu $Kz $R $Cd $db $fc $la
uniaxialMaterial StrPen01 5 0.022 60 0.77 90 0.5 0.65 0.0 1.0 4.35 25.0
set mat 5
# Define truss element
# tag inode jnode A Mattag
element truss 1 1 2 1 $mat
element truss 2 1 3 1 $mat
element truss 3 2 3 1 $mat
#----------DONE Create Model
#----------Create Analysis
# Create displacement command recorders
recorder Node -file commanddisp.out -time -node 2 -dof 1 disp
# Create reaction force recorders (for checking material models)
recorder Element -file reaction.out -time -ele 1 axialForce
pattern Plain 1 "Series -dt 0.01 -filePath cycload.txt" {
sp 2 1 1
}
integrator LoadControl .01
system SparseGeneral -piv
test NormUnbalance 1.0e-6 100
numberer Plain
constraints Transformation
algorithm Newton
analysis Static
set numIncr 310
analyze $numIncr
wipe