Eleload command

[Paolo]

Hello!
I’m working on a 3d structure similar to example 5.
I’m using eleload type beamuniform for gravity load, but I get an error message (expected integer but got “-type”).
Anyway gravity load anlysis run, axial force on the pillars is right, but I obtain a big axial force on the beams and this is wrong.
I tryed this command on a simple portal in 3D space and I got wrong results too.
Even if I insert two element for the beam and just a nodal load on the middle node I get an axial force.

The script is

model BasicBuilder -ndm 3 -ndf 6
set h 300.0
set by 500.0
set bx 500.0
node 1 [expr -$bx/2] [expr $by/2] 0
node 2 [expr $bx/2] [expr $by/2] 0
node 5 [expr -$bx/2] [expr $by/2] $h
node 56 0 [expr $by/2] $h
node 6 [expr $bx/2] [expr $by/2] $h
fix 1 1 1 1 1 1 1
fix 2 1 1 1 1 1 1
geomTransf PDelta 1 0 1 0
element elasticBeamColumn 1 1 5 900 290000 126087 114210 67500 67500 1
element elasticBeamColumn 2 2 6 900 290000 126087 114210 67500 67500 1
# Define beam elements
geomTransf PDelta 3 0 0 1
element elasticBeamColumn 171 5 56 900 290000 126087 114210 67500 67500 3
element elasticBeamColumn 172 56 6 900 290000 126087 114210 67500 67500 3
system BandGeneral
constraints Transformation
numberer RCM
test NormDispIncr 1.0e-12 10 3
algorithm Newton
integrator LoadControl 0.01
analysis Static
initialize
pattern Plain 1 "Linear" {
eleLoad -ele 171 -type -beamUniform 0 -24
eleLoad -ele 172 -type -beamUniform 0 -24
}
# I also tryed with:
# load 56 0 0 -12000 0 0 0
integrator LoadControl 0.01
analyze 100
print ele 1 2 171 172


Do you know what am I mistakin?
Thank you.
Best greatings.

[silvia]

i didn't get the error you got, please test it with the latest version of opensees.
I do get some results, you might want to check your numbers with the node/reaction recorder.
also, you shouldn't have PDelta in your beam transformation!