opensees has stopped working

Forum for OpenSees users to post questions, comments, etc. on the use of the OpenSees interpreter, OpenSees.exe

Moderators: silvia, selimgunay, Moderators

Post Reply
ali2012
Posts: 8
Joined: Sat Feb 04, 2012 5:16 am
Location: Sharif University of Technology

opensees has stopped working

Post by ali2012 »

hello to all
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
Top
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: opensees has stopped working

Post by vesna »

You are not doing your analysis correctly.

First you have to do the static analysis with a gravity load and then you can do the transient analysis.
Top
ali2012
Posts: 8
Joined: Sat Feb 04, 2012 5:16 am
Location: Sharif University of Technology

Re: opensees has stopped working

Post by ali2012 »

thanks a lot vesna
Top
samcivil
Posts: 30
Joined: Tue Apr 17, 2012 10:57 pm
Location: Mohaghegh Ardabili University
Contact:
Contact samcivil

Re: opensees has stopped working

Post by samcivil »

ali2012 wrote:
> hello to all
> 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






hi
please help me about modeling explosive load in opensees!

could i have your email?
Mohaghegh University/Iran
Top
Post Reply

Return to “OpenSees.exe Users”