Scripts Problem:earthquake simulation of 3D steel structure

[bluestar]

I am struggling to model a 1 story steel structure under earthquake loading, but I could not run the analysis successfully! Really need your suggestions!

Below are my scripts, I also provided the OpenSees system Prompts.

wipe

#units kip sec inch
file mkdir EQ


model basic -ndm 3 -ndf 6

node 1 144 0 0
node 2 144 144 0
node 3 0 144 0
node 4 0 0 0
node 5 144 0 144
node 6 144 144 144
node 7 0 144 144
node 8 0 0 144
node 9 72 72 144


mass 5 [expr 10/386] [expr 10/386] [expr 10/386] 0 0 0
mass 6 [expr 10/386] [expr 10/386] [expr 10/386] 0 0 0
mass 7 [expr 10/386] [expr 10/386] [expr 10/386] 0 0 0
mass 8 [expr 10/386] [expr 10/386] [expr 10/386] 0 0 0

fix 1 1 1 1 1 1 1
fix 2 1 1 1 1 1 1
fix 3 1 1 1 1 1 1
fix 4 1 1 1 1 1 1

geomTransf Linear 1 -1 0 0
geomTransf Linear 2 0 0 1



set Ubig 1.0E+10

#column
element elasticBeamColumn 1 1 5 28.2 29500 11346 $Ubig 833 270 1
element elasticBeamColumn 2 2 6 28.2 29500 11346 $Ubig 833 270 1
element elasticBeamColumn 3 3 7 28.2 29500 11346 $Ubig 833 270 1
element elasticBeamColumn 4 4 8 28.2 29500 11346 $Ubig 833 270 1

#Beam
element elasticBeamColumn 5 5 6 28.2 29500 11346 $Ubig 833 270 2
element elasticBeamColumn 6 6 7 28.2 29500 11346 $Ubig 833 270 2
element elasticBeamColumn 7 7 8 28.2 29500 11346 $Ubig 833 270 2
element elasticBeamColumn 8 8 5 28.2 29500 11346 $Ubig 833 270 2



#recorder

recorder Node -file EQ/node0.out -time -nodeRange 5 8 -dof 1 2 3 disp
recorder Element -file EQ/ele0.out -time -eleRange 1 8 localForce
recorder Drift -file EQ/drift0.out -time -iNode 2 -jNode 6 -dof 2 -perpDirn 3


#I thought belowing scripts has problem.

set accelSeries "Series -dt 0.005 -filePath A10000.tcl -factor 1"
pattern UniformExcitation 1 2 -accel $accelSeries

#Expected integer but got "dt"

rayleigh 0. 0. 0. [expr 2*0.02/pow([eigen 1],0.5)]

#ProfileSPDLinDirectSolver::solve<> - aii<0 <i, aii>:<0,0>
#ArpackSolver::Error with _saupd info = -9
#Starting vector is zero
#Warning DirectIntegrationAnalysis:;eigen<> - EigenSOE failed in solve<> in expression "2*0.02/pow<_@_,0.5>"

constraints Plain
numberer Plain
system BandGeneral
test NormDispIncr 1.0e-8 10
algorithm Newton

integrator Newmark 0.5 0.25
analysis Transient

analyze 3995 0.01
#WARNING BandGenLinLapackSolver::solve<> -LAPACK routine returned 25
#WARNING NewtonRaphson::solveCurrentStep<> -the LinearSysOfEqn failed in solve<> #DirecIntegrationAnalysis::analyze<> the Algorithm failed at time 0.02
#Opensees > analyze failed, returned: -3 error flag -3



Looking forward to your reply, Thanks!!!

[fmk]

node 9 is not connected to anything or fixed

[bluestar]

Thanks FMK,

Sorry for my inferior mistake, but after I fix this problem, I still could not run the analysis

Below is my script:

wipe

#units kip sec inch
file mkdir EQ


model basic -ndm 3 -ndf 6

node 1 144 0 0
node 2 144 144 0
node 3 0 144 0
node 4 0 0 0
node 5 144 0 144
node 6 144 144 144
node 7 0 144 144
node 8 0 0 144
node 9 72 72 144



mass 5 [expr 10/386] [expr 10/386] [expr 10/386] 0 0 0
mass 6 [expr 10/386] [expr 10/386] [expr 10/386] 0 0 0
mass 7 [expr 10/386] [expr 10/386] [expr 10/386] 0 0 0
mass 8 [expr 10/386] [expr 10/386] [expr 10/386] 0 0 0

fix 1 1 1 1 1 1 1
fix 2 1 1 1 1 1 1
fix 3 1 1 1 1 1 1
fix 4 1 1 1 1 1 1
rigidDiaphragm 3 9 5 6 7 8

geomTransf Linear 1 -1 0 0
geomTransf Linear 2 0 0 1



set Ubig 1.0E+10

#column
element elasticBeamColumn 1 1 5 28.2 29500 11346 $Ubig 833 270 1
element elasticBeamColumn 2 2 6 28.2 29500 11346 $Ubig 833 270 1
element elasticBeamColumn 3 3 7 28.2 29500 11346 $Ubig 833 270 1
element elasticBeamColumn 4 4 8 28.2 29500 11346 $Ubig 833 270 1

#Beam
element elasticBeamColumn 5 5 6 28.2 29500 11346 $Ubig 833 270 2
element elasticBeamColumn 6 6 7 28.2 29500 11346 $Ubig 833 270 2
element elasticBeamColumn 7 7 8 28.2 29500 11346 $Ubig 833 270 2
element elasticBeamColumn 8 8 5 28.2 29500 11346 $Ubig 833 270 2



#recorder

recorder Node -file EQ/node0.out -time -nodeRange 5 8 -dof 1 2 3 disp
recorder Element -file EQ/ele0.out -time -eleRange 1 8 localForce
recorder Drift -file EQ/drift0.out -time -iNode 2 -jNode 6 -dof 2 -perpDirn 3


#I thought belowing scripts has problem.

set accelSeries "Series -dt 0.005 -filePath A10000.tcl -factor 1"

#A10000.tcl file link:
#http://opensees.berkeley.edu/wiki/image ... A10000.tcl

pattern UniformExcitation 1 2 -accel $accelSeries
rayleigh 0. 0. 0. [expr 2*0.02/pow([eigen 1],0.5)]


constraints Plain
numberer Plain
system BandGeneral
test NormDispIncr 1.0e-8 10
algorithm Newton

integrator Newmark 0.5 0.25
analysis Transient

analyze 3995 0.01

[fmk]

expr 10/386 is equal to ZERO (integer math) .. you want 10/386.0

and you will have to constrain node 9 some as the rigid diaphragm only will limit some displacements and rotations, not all.