pushover analysis on single piles using p-y,t-z,q-z material

[amin_asareh]

model basic -ndm 2 -ndf 3
file mkdir static
file mkdir dynamic
node 1 0 3.81
node 2 0 3
node 3 0 2
node 4 0 1
node 5 0 0
node 6 0 -0.1
node 7 0 -1.1
node 8 0 -2.1
node 9 0 -3.1
node 10 0 -4.1
node 11 0 -5.1
node 12 0 -6.1
node 13 0 -7.1
node 14 0 -8.1
node 15 0 -9.1
node 16 0 -10.1
node 17 0 -11.1
node 18 0 -12.1
node 19 0 -13.1
node 20 0 -14.1
node 21 0 -15.1
node 22 0 -16.1
node 23 0 -17.1
node 24 0 -17.5
node 105 0 0
node 106 0 -0.1
node 107 0 -1.1
node 108 0 -2.1
node 109 0 -3.1
node 110 0 -4.1
node 111 0 -5.1
node 112 0 -6.1
node 113 0 -7.1
node 114 0 -8.1
node 115 0 -9.1
node 116 0 -10.1
node 117 0 -11.1
node 118 0 -12.1
node 119 0 -13.1
node 120 0 -14.1
node 121 0 -15.1
node 122 0 -16.1
node 123 0 -17.1
node 124 0 -17.5
fix 105 1 1 1
fix 106 1 1 1
fix 107 1 1 1
fix 108 1 1 1
fix 109 1 1 1
fix 110 1 1 1
fix 111 1 1 1
fix 112 1 1 1
fix 113 1 1 1
fix 114 1 1 1
fix 115 1 1 1
fix 116 1 1 1
fix 117 1 1 1
fix 118 1 1 1
fix 119 1 1 1
fix 120 1 1 1
fix 121 1 1 1
fix 122 1 1 1
fix 123 1 1 1
fix 124 1 1 1
#fix 24 1 1 1

geomTransf Corotational 1

# steel material for piles
#uniaxialMaterial Steel02 $matTag $Fy $E $b $R0 $cR1 $cR2 $a1 $a2 $a3 $a4
uniaxialMaterial Steel02 1 290000000 7.00E+010 0.003 20 0.925 0.15

# soil material for lateral bearing P-Y
#uniaxialMaterial PySimple1 $matTag $soilType $pult $Y50 $Cd <$c>.

uniaxialMaterial PySimple1 10 1 1000.5 0.0083375 0.3 282932.5992
uniaxialMaterial PySimple1 11 1 2397.84 0.0083375 0.3 282932.5992
uniaxialMaterial PySimple1 12 1 11468.24 0.0083375 0.3 282932.5992
uniaxialMaterial PySimple1 13 1 25589.64 0.0083375 0.3 282932.5992
uniaxialMaterial PySimple1 14 1 37384.54 0.0083375 0.3 282932.5992
uniaxialMaterial PySimple1 15 1 50679.44 0.0083375 0.3 282932.5992
uniaxialMaterial PySimple1 16 1 76851.84 0.0083375 0.3 282932.5992
uniaxialMaterial PySimple1 17 1 90045.00 0.0083375 0.3 282932.5992
uniaxialMaterial PySimple1 18 2 1694830.00 0.004430 0.1 1099084.328
uniaxialMaterial PySimple1 19 2 2178968.45 0.004992 0.1 1099084.328
uniaxialMaterial PySimple1 20 2 2723706.90 0.005554 0.1 1099084.328
uniaxialMaterial PySimple1 21 2 3329045.35 0.006117 0.1 1099084.328
uniaxialMaterial PySimple1 22 2 3994983.80 0.006679 0.1 1099084.328
uniaxialMaterial PySimple1 23 2 4483273.85 0.006876 0.1 1099084.328
uniaxialMaterial PySimple1 24 2 4853792.35 0.006876 0.1 1099084.328
uniaxialMaterial PySimple1 25 2 5224310.85 0.006876 0.1 1099084.328
uniaxialMaterial PySimple1 26 2 5594829.35 0.006876 0.1 1099084.328
uniaxialMaterial PySimple1 27 2 5965347.85 0.006876 0.1 1099084.328
uniaxialMaterial PySimple1 28 2 6335866.35 0.006876 0.1 1099084.328
uniaxialMaterial PySimple1 29 2 2593629.50 0.006876 0.1 1099084.328


#uniaxialMaterial TzSimple1 $matTag $tzType $tult $z50 <$c>

uniaxialMaterial TzSimple1 110 2 533.9642244 0.0020677 282932.5992
uniaxialMaterial TzSimple1 111 2 2094.38 0.0020677 282932.5992
uniaxialMaterial TzSimple1 112 2 7737.876745 0.0020677 282932.5992
uniaxialMaterial TzSimple1 113 2 11123.89869 0.0020677 282932.5992
uniaxialMaterial TzSimple1 114 2 14505.75511 0.0020677 282932.5992
uniaxialMaterial TzSimple1 115 2 19443.92433 0.0020677 282932.5992
uniaxialMaterial TzSimple1 116 2 22842.18819 0.0020677 282932.5992
uniaxialMaterial TzSimple1 117 2 34663.1718 0.00127 1099084.328
uniaxialMaterial TzSimple1 118 2 39545.30867 0.00127 1099084.328
uniaxialMaterial TzSimple1 119 2 44427.44554 0.00127 1099084.328
uniaxialMaterial TzSimple1 120 2 49309.58242 0.00127 1099084.328
uniaxialMaterial TzSimple1 121 2 54191.71929 0.00127 1099084.328
uniaxialMaterial TzSimple1 122 2 59073.85616 0.00127 1099084.328
uniaxialMaterial TzSimple1 123 2 63955.99304 0.00127 1099084.328
uniaxialMaterial TzSimple1 124 2 68838.12991 0.00127 1099084.328
uniaxialMaterial TzSimple1 125 2 73720.26678 0.00127 1099084.328
uniaxialMaterial TzSimple1 126 2 78602.40366 0.00127 1099084.328
uniaxialMaterial TzSimple1 127 2 83484.54053 0.00127 1099084.328
uniaxialMaterial TzSimple1 128 2 34174.95811 0.00127 1099084.328
uniaxialMaterial QzSimple1 129 2 341654.922 0.008671 0 1099084.328
#--------------------- ---------------------------
#creat section Fiber
set indiam 0.3145
set outdiam 0.3335

section Fiber 1 {

# Create the pile fibers
patch circ 1 20 4 0 0 $indiam $outdiam 0 360
}


# pile elements
#element dispBeamColumn $eleTag $iNode $jNode $numIntgrPts $secTag $transfTag <-mass $massDens>
element dispBeamColumn 1 1 2 10 1 1 -mass 3700
element dispBeamColumn 2 2 3 10 1 1 -mass 3700
element dispBeamColumn 3 3 4 10 1 1 -mass 3700
element dispBeamColumn 4 4 5 10 1 1 -mass 3700
element dispBeamColumn 5 5 6 10 1 1 -mass 3700
element dispBeamColumn 6 6 7 10 1 1 -mass 3700
element dispBeamColumn 7 7 8 10 1 1 -mass 3700
element dispBeamColumn 8 8 9 10 1 1 -mass 3700
element dispBeamColumn 9 9 10 10 1 1 -mass 3700
element dispBeamColumn 10 10 11 10 1 1 -mass 3700
element dispBeamColumn 11 11 12 10 1 1 -mass 3700
element dispBeamColumn 12 12 13 10 1 1 -mass 3700
element dispBeamColumn 13 13 14 10 1 1 -mass 3700
element dispBeamColumn 14 14 15 10 1 1 -mass 3700
element dispBeamColumn 15 15 16 10 1 1 -mass 3700
element dispBeamColumn 16 16 17 10 1 1 -mass 3700
element dispBeamColumn 17 17 18 10 1 1 -mass 3700
element dispBeamColumn 18 18 19 10 1 1 -mass 3700
element dispBeamColumn 19 19 20 10 1 1 -mass 3700
element dispBeamColumn 20 20 21 10 1 1 -mass 3700
element dispBeamColumn 21 21 22 10 1 1 -mass 3700
element dispBeamColumn 22 22 23 10 1 1 -mass 3700
element dispBeamColumn 23 23 24 10 1 1 -mass 3700

#zero length Elements for soil
#element zeroLength $eleTag $iNode $jNode -mat $matTag1 $matTag2 -dir $dir1 $dir2
element zeroLength 101 5 105 -mat 10 -dir 1
element zeroLength 102 6 106 -mat 11 110 -dir 1 2
element zeroLength 103 7 107 -mat 12 111 -dir 1 2
element zeroLength 104 8 108 -mat 13 112 -dir 1 2
element zeroLength 105 9 109 -mat 14 113 -dir 1 2
element zeroLength 106 10 110 -mat 15 114 -dir 1 2
element zeroLength 107 11 111 -mat 16 115 -dir 1 2
element zeroLength 108 12 112 -mat 17 116 -dir 1 2
element zeroLength 109 13 113 -mat 18 117 -dir 1 2
element zeroLength 110 14 114 -mat 19 118 -dir 1 2
element zeroLength 111 15 115 -mat 20 119 -dir 1 2
element zeroLength 112 16 116 -mat 21 120 -dir 1 2
element zeroLength 113 17 117 -mat 22 121 -dir 1 2
element zeroLength 114 18 118 -mat 23 122 -dir 1 2
element zeroLength 115 19 119 -mat 24 123 -dir 1 2
element zeroLength 116 20 120 -mat 25 124 -dir 1 2
element zeroLength 117 21 121 -mat 26 125 -dir 1 2
element zeroLength 118 22 122 -mat 27 126 -dir 1 2
element zeroLength 119 23 123 -mat 28 127 -dir 1 2
element zeroLength 120 24 124 -mat 29 129 -dir 1 2
recorder Node -file static/staticReaction1.out -time -nodeRange 1 89 -dof 1 2 3 reaction
recorder Node -file static/staticanalysisdisp.out -time -nodeRange 1 89 -dof 1 2 3 disp

print node 1
puts "Gravity Load Analysis Complete"
eigen 3
print NodeProp.txt -node
puts "T1 = [expr 2*3.1416/(pow([eigen 1] ,0.5))] sec"

recorder Node -file dynamic/Node1.out -time -node 1 -dof 1 2 disp
recorder Node -file dynamic/dispNode.out -time -nodeRange 1 24 -dof 1 2 disp
recorder Node -file dynamic/accelNodepile.out -time -nodeRange 1 24 -dof 1 2 accel
recorder Node -file dynamic/accelNodesoil.out -time -nodeRange 105 124 -dof 1 2 accel
recorder Node -file dynamic/reactionNodesoil.out -time -nodeRange 105 124 -dof 1 2 reaction
recorder Drift -file dynamic/driftrecorder(1-7).out -time -iNode 1 2 3 4 6 -jNode 2 3 4 5 7 -dof 1 -perpDirn 2
recorder Drift -file dynamic/driftrecorder(8-13).out -time -iNode 8 9 10 11 12 -jNode 9 10 11 12 13 -dof 1 -perpDirn 2
recorder Drift -file dynamic/driftrecorder(13-17).out -time -iNode 13 14 15 16 17 -jNode 14 15 16 17 18 -dof 1 -perpDirn 2
recorder Drift -file dynamic/driftrecorder(17-24).out -time -iNode 18 19 20 21 22 -jNode 19 20 21 22 23 -dof 1 -perpDirn 2
recorder Element -file dynamic/globalforce.out -time -eleRange 1 23 globalForce
recorder Element -file dynamic/sectionstressstrain.out -time -eleRange 1 23 section 1 fiber 0 350 stressStrain

recorder plot dynamic/Node1.out dispoftopofthepile 0 0 600 600 -columns 1 2

recorder display pile 600 0 400 400 -wipe
prp -100 20 30

vup 0 1 0

display 1 5 100

pattern Plain 1 Linear { load 1 1.0 0.0 0.0 0 0 0}

test EnergyIncr 1.0e-8 30 0
algorithm Newton
system BandGeneral
numberer RCM
constraints Plain

integrator DisplacementControl 1 1 0.001

analysis Static
analyze 2000

[/code]