This warning ocurrs every time when I go through my code below. As a beginner, I wonder what's wrong with my codes? please help me point them out, thank you!puts "set parameters"
set PI 3.14159
wipe;
model basic -ndm 3 -ndf 6;
## Define the model builder, ndm = #dimension, ndf = #dofs
puts "start modeling"
node 111 0 0 0
node 211 5000 0 0
node 311 10000 0 0
node 121 0 6000 0
node 221 5000 6000 0
node 321 10000 6000 0
node 112 0 0 3700
node 212 5000 0 3700
node 312 10000 0 3700
node 122 0 6000 3700
node 222 5000 6000 3700
node 322 10000 6000 3700
node 113 0 0 7200
node 213 5000 0 7200
node 313 10000 0 7200
node 123 0 6000 7200
node 223 5000 6000 7200
node 323 10000 6000 7200
node 114 0 0 10700
node 214 5000 0 10700
node 314 10000 0 10700
node 124 0 6000 10700
node 224 5000 6000 10700
node 324 10000 6000 10700
node 115 0 0 14200
node 215 5000 0 14200
node 315 10000 0 14200
node 125 0 6000 14200
node 225 5000 6000 14200
node 325 10000 6000 14200
fix 111 1 1 1 1 1 1
fix 211 1 1 1 1 1 1
fix 311 1 1 1 1 1 1
fix 121 1 1 1 1 1 1
fix 221 1 1 1 1 1 1
fix 321 1 1 1 1 1 1
mass 112 6.05 6.05 6.05 0 0 0
mass 212 12.10 12.10 12.10 0 0 0
mass 312 6.05 6.05 6.05 0 0 0
mass 122 6.05 6.05 6.05 0 0 0
mass 222 12.10 12.10 12.10 0 0 0
mass 322 6.05 6.05 6.05 0 0 0
mass 113 6.03 6.03 6.03 0 0 0
mass 213 12.07 12.07 12.07 0 0 0
mass 313 6.03 6.03 6.03 0 0 0
mass 123 6.03 6.03 6.03 0 0 0
mass 223 12.07 12.07 12.07 0 0 0
mass 323 6.03 6.03 6.03 0 0 0
mass 114 6.08 6.08 6.08 0 0 0
mass 214 12.16 12.16 12.16 0 0 0
mass 314 6.08 6.08 6.08 0 0 0
mass 124 6.08 6.08 6.08 0 0 0
mass 224 12.16 12.16 12.16 0 0 0
mass 324 6.08 6.08 6.08 0 0 0
mass 115 8.05 8.05 8.05 0 0 0
mass 215 16.11 16.11 16.11 0 0 0
mass 315 8.05 8.05 8.05 0 0 0
mass 125 8.05 8.05 8.05 0 0 0
mass 225 16.11 16.11 16.11 0 0 0
mass 325 8.05 8.05 8.05 0 0 0
node 2 5000 3000 3700
node 3 5000 3000 7200
node 4 5000 3000 10700
node 5 5000 3000 14200
fix 2 0 0 1 1 1 0
fix 3 0 0 1 1 1 0
fix 4 0 0 1 1 1 0
fix 5 0 0 1 1 1 0
rigidDiaphragm 3 2 112 122 212 222 312 322
rigidDiaphragm 3 3 113 123 213 223 313 323
rigidDiaphragm 3 4 114 124 214 224 314 324
rigidDiaphragm 3 5 115 125 215 225 315 325
#BCR295
uniaxialMaterial Steel01 1 295.0 205000 0.05
#uniaxialMaterial Hysteretic 1 295.0 1.44e-3 400.0 0.0117 -295.0 -1.44e-3 -400.0 -0.0117 1 1 0 0 0
#SN400B SS400
uniaxialMaterial Steel01 2 235.0 205000 0.05
#add shear and torsion properties to cols and beams
#col V T
uniaxialMaterial Elastic 101 4.258e8
uniaxialMaterial Elastic 102 1.749e13
#beam1 VY VZ T
uniaxialMaterial Elastic 201 1.932e8
uniaxialMaterial Elastic 202 2.300e8
uniaxialMaterial Elastic 203 1.180e10
#beam2 VY VZ T
uniaxialMaterial Elastic 301 3.154e8
uniaxialMaterial Elastic 302 3.416e8
uniaxialMaterial Elastic 303 3.181e10
#col
section Fiber 1 {
# #material numY numZ I J#
patch rect 1 5 1 -300 291 300 300
patch rect 1 5 1 -300 -300 300 -291
patch rect 1 1 5 -300 -291 -291 291
patch rect 1 1 5 291 -291 300 291
}
#beam
#RG1,11,12
section Fiber 2 {
# material numY numZ I J
patch rect 2 1 5 164 -87 173 87
patch rect 2 1 5 -173 -87 -164 87
patch rect 2 5 1 -164 -3 164 3
}
#4G1,11
section Fiber 3 {
# material numY numZ I J
patch rect 2 1 5 164 -87.5 175 87.5
patch rect 2 1 5 -175 -87.5 -164 87.5
patch rect 2 5 1 -164 -3.5 164 3.5
}
#4G12
section Fiber 4 {
# material numY numZ I J
patch rect 2 1 5 156 -87.5 170 87.5
patch rect 2 1 5 -170 -87.5 -156 87.5
patch rect 2 5 1 -156 -4.5 156 4.5
}
#3G1
section Fiber 5 {
# material numY numZ I J
patch rect 2 1 5 187 -99.5 198 99.5
patch rect 2 1 5 -198 -99.5 -187 99.5
patch rect 2 5 1 -187 -3.5 187 3.5
}
#3G11,12
#2G1,11
section Fiber 6 {
# material numY numZ I J
patch rect 2 1 5 187 -100 200 100
patch rect 2 1 5 -200 -100 -187 100
patch rect 2 5 1 -187 -4 187 4
}
#2G12
section Fiber 7 {
# material numY numZ I J
patch rect 2 1 5 179 -100 195 100
patch rect 2 1 5 -195 -100 -179 100
patch rect 2 5 1 -179 -5 179 5
}
section Aggregator 1001 101 Vy 101 Vz 102 T -section 1
section Aggregator 1002 201 Vy 202 Vz 203 T -section 2
section Aggregator 1003 201 Vy 202 Vz 203 T -section 3
section Aggregator 1004 201 Vy 202 Vz 203 T -section 4
section Aggregator 1005 301 Vy 302 Vz 303 T -section 5
section Aggregator 1006 301 Vy 302 Vz 303 T -section 6
section Aggregator 1007 301 Vy 302 Vz 303 T -section 7
geomTransf PDelta 1 0 1 0
element dispBeamColumn 111112 111 112 5 1001 1
element dispBeamColumn 211212 211 212 5 1001 1
element dispBeamColumn 311312 311 312 5 1001 1
element dispBeamColumn 121122 121 122 5 1001 1
element dispBeamColumn 221222 221 222 5 1001 1
element dispBeamColumn 321322 321 322 5 1001 1
element dispBeamColumn 112113 112 113 5 1001 1
element dispBeamColumn 212213 212 213 5 1001 1
element dispBeamColumn 312313 312 313 5 1001 1
element dispBeamColumn 122123 122 123 5 1001 1
element dispBeamColumn 222223 222 223 5 1001 1
element dispBeamColumn 322323 322 323 5 1001 1
element dispBeamColumn 113114 113 114 5 1001 1
element dispBeamColumn 213214 213 214 5 1001 1
element dispBeamColumn 313314 313 314 5 1001 1
element dispBeamColumn 123124 123 124 5 1001 1
element dispBeamColumn 223224 223 224 5 1001 1
element dispBeamColumn 323324 323 324 5 1001 1
element dispBeamColumn 114115 114 115 5 1001 1
element dispBeamColumn 214215 214 215 5 1001 1
element dispBeamColumn 314315 314 315 5 1001 1
element dispBeamColumn 124125 124 125 5 1001 1
element dispBeamColumn 224225 224 225 5 1001 1
element dispBeamColumn 324325 324 325 5 1001 1
geomTransf Linear 2 0 -1 0
geomTransf Linear 3 1 0 0
element dispBeamColumn 112212 112 212 5 1006 2
element dispBeamColumn 212312 212 312 5 1006 2
element dispBeamColumn 122222 122 222 5 1006 2
element dispBeamColumn 222322 222 322 5 1006 2
element dispBeamColumn 113213 113 213 5 1005 2
element dispBeamColumn 213313 213 313 5 1005 2
element dispBeamColumn 123223 123 223 5 1005 2
element dispBeamColumn 223323 223 323 5 1005 2
element dispBeamColumn 114214 114 214 5 1003 2
element dispBeamColumn 214314 214 314 5 1003 2
element dispBeamColumn 124224 124 224 5 1003 2
element dispBeamColumn 224324 224 324 5 1003 2
element dispBeamColumn 115215 115 215 5 1002 2
element dispBeamColumn 215315 215 315 5 1002 2
element dispBeamColumn 125225 125 225 5 1002 2
element dispBeamColumn 225325 225 325 5 1002 2
element dispBeamColumn 112122 112 122 5 1006 3
element dispBeamColumn 312322 312 322 5 1006 3
element dispBeamColumn 113123 113 123 5 1006 3
element dispBeamColumn 313323 313 323 5 1006 3
element dispBeamColumn 114124 114 124 5 1003 3
element dispBeamColumn 314324 314 324 5 1003 3
element dispBeamColumn 115125 115 125 5 1002 3
element dispBeamColumn 315325 315 325 5 1002 3
element dispBeamColumn 212222 212 222 5 1007 3
element dispBeamColumn 213223 213 223 5 1006 3
element dispBeamColumn 214224 214 224 5 1004 3
element dispBeamColumn 215225 215 225 5 1002 3
puts "static analysis"
pattern Plain 1 Linear {
load 112 0 0 -59313.07 0 0 0 ;
load 122 0 0 -59313.07 0 0 0 ;
load 312 0 0 -59313.07 0 0 0 ;
load 322 0 0 -59313.07 0 0 0 ;
load 113 0 0 -59125.52 0 0 0 ;
load 123 0 0 -59125.52 0 0 0 ;
load 313 0 0 -59125.52 0 0 0 ;
load 323 0 0 -59125.52 0 0 0 ;
load 114 0 0 -59561.91 0 0 0 ;
load 124 0 0 -59561.91 0 0 0 ;
load 314 0 0 -59561.91 0 0 0 ;
load 324 0 0 -59561.91 0 0 0 ;
load 115 0 0 -78937.40 0 0 0 ;
load 125 0 0 -78937.40 0 0 0 ;
load 315 0 0 -78937.40 0 0 0 ;
load 325 0 0 -78937.40 0 0 0 ;
load 212 0 0 -118626.14 0 0 0 ;
load 222 0 0 -118626.14 0 0 0 ;
load 213 0 0 -118251.04 0 0 0 ;
load 223 0 0 -118251.04 0 0 0 ;
load 214 0 0 -119123.83 0 0 0 ;
load 224 0 0 -119123.83 0 0 0 ;
load 215 0 0 -157874.81 0 0 0 ;
load 225 0 0 -157874.81 0 0 0 ;
}
set Tol 1.0e-6;
constraints Transformation;
numberer RCM;
system BandGeneral
test NormDispIncr $Tol 200;
algorithm Newton;
integrator LoadControl 0.1;
analysis Static;
analyze 10
puts "transient analysis"
# maintain constant gravity loads and reset time to zero
loadConst -time 0.0
set pDamp 0.02
set lambda [eigen 10]
set omegaI [expr pow([lindex $lambda 0],0.5)];
set omegaJ [expr pow([lindex $lambda 1],0.5)];
set alphaM [expr $pDamp*(2*$omegaI*$omegaJ)/($omegaI+$omegaJ)];
set betaKcomm [expr 2.*$pDamp/($omegaI+$omegaJ)];
#rayleigh $alphaM 0. 0. $betaKcomm
rayleigh $alphaM $betaKcomm 0. 0.
recorder Node -file output/eigen1.out -time -nodeRange 112 115 -dof 1 2 "eigen 1"
recorder Node -file output/eigen2.out -time -nodeRange 112 115 -dof 1 2 "eigen 2"
set period "period.txt"
set periods [open $period "w"]
puts $periods "$lambda"
close $periods
#record
set rey "reyleigh.txt"
set reys [open $rey "w"]
puts $reys "$alphaM"
puts $reys "$betaKcomm"
close $reys
#record
puts " Period T1: [expr 2*$PI/sqrt([lindex $lambda 0])] sec"
puts " Period T1: [expr 2*$PI/sqrt([lindex $lambda 1])] sec"
set fac {0.5 1.0 1.25 2.0 4.0 6.0 10.0}
timeSeries Path 1 -filePath TAKNS.txt -dt 0.01 -factor [lindex $fac 5]
timeSeries Path 2 -filePath TAKEW.txt -dt 0.01 -factor [lindex $fac 5]
timeSeries Path 3 -filePath TAKUD.txt -dt 0.01 -factor [lindex $fac 5]
pattern UniformExcitation 2 1 -accel 1
pattern UniformExcitation 3 2 -accel 2
pattern UniformExcitation 4 3 -accel 3
set colnodeList []
for {set floor 1} {$floor <= 5} {incr floor 1} {
lappend colnodeList 11$floor 21$floor 31$floor
}
set col1nodeList []
set col2nodeList []
set nodeList1 []
set nodeList2 []
set coleleList []
for {set floor1 1; set floor2 2;} {$floor1 <= 4 && $floor2 <= 5} {incr floor1 1; incr floor2 1;} {
lappend col1nodeList 11$floor1 21$floor1 31$floor1
lappend col2nodeList 11$floor2 21$floor2 31$floor2
lappend nodeList1 $floor1
lappend nodeList2 $floor2
lappend coleleList 11[expr $floor1]11$floor2 21[expr $floor1]21$floor2 31[expr $floor1]31$floor2
}
set cmd "recorder EnvelopeNode -file output/floorAccEnvx.out -timeSeries 1 -node $colnodeList -dof 1 accel"; eval $cmd;
set cmd "recorder EnvelopeNode -file output/floorAccEnvy.out -timeSeries 2 -node $colnodeList -dof 2 accel"; eval $cmd;
set cmd "recorder EnvelopeNode -file output/floorDispEnv.out -node $colnodeList -dof 1 2 disp"; eval $cmd;
set cmd "recorder EnvelopeNode -file output/floorReactionEnv.out -node $colnodeList -dof 1 2 reaction"; eval $cmd;
#set cmd "recorder Node -file output/floorDisp.out -time -node $nodeList -dof 1 2 disp"; eval $cmd;
#set cmd "recorder Node -file output/floorReaction.out -time -node $nodeList -dof 1 2 reaction"; eval $cmd;
#set cmd "recorder Drift -file output/RIGIdriftx.out -time -iNode 1 2 3 4 -jNode 2 3 4 5 -dof 1 -perpDirn 3"; eval $cmd;
#set cmd "recorder Drift -file output/RIGIdrifty.out -time -iNode 1 2 3 4 -jNode 2 3 4 5 -dof 2 -perpDirn 3"; eval $cmd;
#
set cmd "recorder Node -file output/floorDisp.out -time -node $colnodeList -dof 1 2 disp"; eval $cmd;
set cmd "recorder Node -file output/floorReaction.out -time -node $colnodeList -dof 1 2 3 4 5 6 reaction"; eval $cmd;
set cmd "recorder Drift -file output/COLdriftx.out -time -iNode $col1nodeList -jNode $col2nodeList -dof 1 -perpDirn 3"; eval $cmd;
set cmd "recorder Drift -file output/COLdriftxbot.out -time -iNode 111 -jNode 112 -dof 1 -perpDirn 3"; eval $cmd;
set cmd "recorder Drift -file output/COLdrifty.out -time -iNode $col1nodeList -jNode $col2nodeList -dof 2 -perpDirn 3"; eval $cmd;
set cmd "recorder Element -file output/COLELEforce.out -time -ele $coleleList section 1 force"; eval $cmd;
set cmd "recorder Element -file output/COLELEdeformation.out -time -ele $coleleList section 1 deformation"; eval $cmd;
#bottom rotation
set cmd "recorder Drift -file output/COLrotation.out -time -iNode 111 211 311 -jNode 112 212 312 -dof 5 -perpDirn 3"; eval $cmd;
constraints Transformation
numberer RCM
system UmfPack
test NormDispIncr 1.0e-1 200
algorithm NewtonLineSearch 0.7
integrator Newmark 0.5 0.25
#integrator HHT 0.9
analysis Transient
analyze 4000 0.01