i wrote a code. in this code i wanted get the maximum displacement of a beam on the two simple supported in the explosion load. also the natural frequencies of the beam. i devided the beam into the 4 elements by elasticbeamcolumn element.analysis is elastic performed. but when i run the program i face to this error:
openSees has stopped working. the length of the beam is four meter. what is the problem. please help me!!!!
code:
# vahedha=meter, kg, sec, N
wipe
model basic -ndm 2 -ndf 3
file mkdir data;
set A 0.05
set IZ 0.0004
set E 2e11
#dead load= N/m
set deadload [expr $A*7800*9.81]
geomTransf Linear 1
#define bare enfejar
set enfejar "Series -dt 0.001 -filePath load.txt -factor 1"
source elements.tcl

fix 1 1 1 0
fix 5 1 1 0
#define Recorders
recorder Node -file data/disp3TH.out -time -node 3 -dof 2 disp
recorder Node -file data/R1TH.out -time -node 1 -dof 2 reaction
recorder Node -file data/R2TH.out -time -node 5 -dof 2 reaction
recorder Node -file data/eigenvector1.out -time -dof 2 "eigen 1"
recorder Node -file data/eigenvector2.out -time -dof 2 "eigen 2"
recorder Node -file data/eigenvector3.out -time -dof 2 "eigen 3"
eigen 5

constraints Plain
numberer Plain
system BandSPD
test EnergyIncr 1e-8 20
algorithm ModifiedNewton
integrator Newmark 0.5 0.25
analysis Transient
analyze 40 0.001
puts "Done!"
element.tcl code:
node 1 0 0
node 2 1 0
node 3 2 0
node 4 3 0
node 5 4 0
element elasticBeamColumn 1 1 2 $A $E $IZ 1
element elasticBeamColumn 2 2 3 $A $E $IZ 1
element elasticBeamColumn 3 3 4 $A $E $IZ 1
element elasticBeamColumn 4 4 5 $A $E $IZ 1
pattern Plain 1 Constant {
eleLoad -ele 1 -type -beamUniform $deadload
}
pattern Plain 2 Constant {
eleLoad -ele 2 -type -beamUniform $deadload
}
pattern Plain 3 Constant {
eleLoad -ele 3 -type -beamUniform $deadload
}
pattern Plain 4 Constant {
eleLoad -ele 4 -type -beamUniform $deadload
}
pattern Plain 5 $enfejar {
eleLoad -ele 1 -type -beamUniform 1
}
pattern Plain 6 $enfejar {
eleLoad -ele 2 -type -beamUniform 1
}
pattern Plain 7 $enfejar {
eleLoad -ele 3 -type -beamUniform 1
}
pattern Plain 8 $enfejar {
eleLoad -ele 4 -type -beamUniform 1
}
mass 2 0 [expr $deadload*1] 1e-9
mass 3 0 [expr $deadload*1] 1e-9
mass 4 0 [expr $deadload*1] 1e-9
mass 1 0 [expr $deadload*0.5] 1e-9
mass 5 0 [expr $deadload*0.5] 1e-9