eleLoad beamuniform - Diagonals

[sobeli]

Hi everyone!,
I have been doing some calibrations with a 3D frame with one and two diagonals where all the elements share the same characteristics, I have gotten the reactions in 3 models, the first and the second have only one diagonal and the reactions are pretty simmilar to etabs, but when I put both diagonals in a third model my force in Z an moments in X are totally different, I've checked my geomtransf, I 've changed my conectivities in node i to node j, and I have rechecked my eleload beamuniform in several combinatios $wy,$wz,$wx, (these examples are only for dead load of the elements) but my results are the same following the local axis system. I hope someone can guide me, I think this is not a real problem in modelling but I can't see my mistake yet.

I share my 3 basic models..

regards,

********************************************************************************************************
MODEL 1

puts " DIMENSIONS...."
wipe
model BasicBuilder -ndm 3 -ndf 6
set NumModos 1;
set carpeta Resultados;
puts " NODES...."
node 1 0 0 0
node 2 3 0 0
node 3 3 0 3
node 4 0 0 3
node 5 0 3 0
node 6 3 3 0
node 7 3 3 3
node 8 0 3 3
puts " DOF...."
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
puts " GEOTRANSF..."
geomTransf PDelta 1 -3 0 0
geomTransf PDelta 2 -3 0 0
geomTransf PDelta 3 -3 0 0
geomTransf PDelta 4 -3 0 0
geomTransf PDelta 5 0 3 0
geomTransf PDelta 6 3 0 0
geomTransf PDelta 7 0 3 0
geomTransf PDelta 8 3 0 0
geomTransf PDelta 9 -3 0 0
puts " MATERIAL..."
set E1 20389019158.0
set Fy1 35153481.31
uniaxialMaterial Steel02 1000 $Fy1 $E1 1e-9 20 0.925 0.15 0.0001 0.0001 0.0001 0.0001
#[expr {Dens_mat} * {Area_sec}]
set g 9.81
set Dens1 [expr 7849.05 * 0.005969026041820605]
set Dens2 [expr 7849.05 * 0.005969026041820605]
set Dens3 [expr 7849.05 * 0.005969026041820605]
set Dens4 [expr 7849.05 * 0.005969026041820605]
set Dens5 [expr 7849.05 * 0.005969026041820605]
set Dens6 [expr 7849.05 * 0.005969026041820605]
set Dens7 [expr 7849.05 * 0.005969026041820605]
set Dens8 [expr 7849.05 * 0.005969026041820605]
set Dens9 [expr 7849.05 * 0.005969026041820605]
puts " SECTION TAGS..."
set DL1 2000
puts " ELEMENTS..."
section Fiber $DL1 -GJ 423464.244 {
patch circ 1000 5 5 0.0 0.0 0.09 0.10 0.0 360.0
}

puts " ELEMENTS..."
element dispBeamColumn 3000 1 5 4 $DL1 1 -mass [expr $Dens1/$g] -integration Lobatto
element dispBeamColumn 3001 2 6 4 $DL1 2 -mass [expr $Dens2/$g] -integration Lobatto
element dispBeamColumn 3002 3 7 4 $DL1 3 -mass [expr $Dens3/$g] -integration Lobatto
element dispBeamColumn 3003 4 8 4 $DL1 4 -mass [expr $Dens4/$g] -integration Lobatto
element dispBeamColumn 3004 5 6 4 $DL1 5 -mass [expr $Dens5/$g] -integration Lobatto
element dispBeamColumn 3005 6 7 4 $DL1 6 -mass [expr $Dens6/$g] -integration Lobatto
element dispBeamColumn 3006 8 7 4 $DL1 7 -mass [expr $Dens7/$g] -integration Lobatto
element dispBeamColumn 3007 5 8 4 $DL1 8 -mass [expr $Dens8/$g] -integration Lobatto
element dispBeamColumn 3008 1 8 4 $DL1 9 -mass [expr $Dens9/$g] -integration Lobatto
puts " TIMESERIES... "
set timeLinear 1
timeSeries Linear $timeLinear
puts " LOAD PATTERN... "
pattern Plain 1 $timeLinear {
#(- - +) original
eleLoad -ele 3000 -type -beamUniform 0.0 0.0 -46.85118385355202
eleLoad -ele 3001 -type -beamUniform 0.0 0.0 -46.85118385355202
eleLoad -ele 3002 -type -beamUniform 0.0 0.0 -46.85118385355202
eleLoad -ele 3003 -type -beamUniform 0.0 0.0 -46.85118385355202
eleLoad -ele 3004 -type -beamUniform 0.0 -46.85118385355202 0.0
eleLoad -ele 3005 -type -beamUniform 46.85118385355202 0.0 0.0
eleLoad -ele 3006 -type -beamUniform 0.0 -46.85118385355202 0.0
eleLoad -ele 3007 -type -beamUniform 46.85118385355202 0.0 0.0
eleLoad -ele 3008 -type -beamUniform -33.12878980946432 0.0 -33.12878980946432
}
puts " STATIC OUTPUT...";
#records
recorder Node -file ReaccionBase1.out -time -node 1 -dof 1 2 3 4 5 6 reaction -timeSeries 1;
recorder Node -file ReaccionBase2.out -time -node 2 -dof 1 2 3 4 5 6 reaction -timeSeries 1;
recorder Node -file ReaccionBase3.out -time -node 3 -dof 1 2 3 4 5 6 reaction -timeSeries 1;
recorder Node -file ReaccionBase4.out -time -node 4 -dof 1 2 3 4 5 6 reaction -timeSeries 1;
puts "STATIC ANALYSIS...";
system UmfPack
numberer RCM
constraints Plain
test NormDispIncr 1.0e-5 10 0;
integrator LoadControl 0.1
algorithm Newton
analysis Static
analyze 10
loadConst -time 0.0

********************************************************************************************************
#MODEL 2

puts " DIMENSIONS...."
wipe
model BasicBuilder -ndm 3 -ndf 6
set NumModos 1; # Se define el número de modos de vibrar a obtener.
set pi 3.1416;
set carpeta Resultados; # Se define el nombre de la carpeta donde se guardarán los resultados.
puts " NODES...."
node 1 0 0 0
node 2 3 0 0
node 3 3 0 3
node 4 0 0 3
node 5 0 3 0
node 6 3 3 0
node 7 3 3 3
node 8 0 3 3
puts " DOF...."
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
puts " GEOTRANSF..."
geomTransf PDelta 1 -3 0 0
geomTransf PDelta 2 -3 0 0
geomTransf PDelta 3 -3 0 0
geomTransf PDelta 4 -3 0 0
geomTransf PDelta 5 0 3 0
geomTransf PDelta 6 3 0 0
geomTransf PDelta 7 0 3 0
geomTransf PDelta 8 3 0 0
geomTransf PDelta 9 -3 0 0
puts " MATERIAL..."
set E1 20389019158.0
set Fy1 35153481.31
uniaxialMaterial Steel02 1000 $Fy1 $E1 1e-9 20 0.925 0.15 0.0001 0.0001 0.0001 0.0001
#[expr {Dens_mat} * {Area_sec}]
set g 9.81
set Dens1 [expr 7849.05 * 0.005969026041820605]
set Dens2 [expr 7849.05 * 0.005969026041820605]
set Dens3 [expr 7849.05 * 0.005969026041820605]
set Dens4 [expr 7849.05 * 0.005969026041820605]
set Dens5 [expr 7849.05 * 0.005969026041820605]
set Dens6 [expr 7849.05 * 0.005969026041820605]
set Dens7 [expr 7849.05 * 0.005969026041820605]
set Dens8 [expr 7849.05 * 0.005969026041820605]
set Dens9 [expr 7849.05 * 0.005969026041820605]
puts " SECTION TAGS..."
set DL1 2000
puts " ELEMENTS..."
section Fiber $DL1 -GJ 423464.244 {
patch circ 1000 5 5 0.0 0.0 0.09 0.10 0.0 360.0
}
puts " ELEMENTS..."
element dispBeamColumn 3000 1 5 4 $DL1 1 -mass [expr $Dens1/$g] -integration Lobatto
element dispBeamColumn 3001 2 6 4 $DL1 2 -mass [expr $Dens2/$g] -integration Lobatto
element dispBeamColumn 3002 3 7 4 $DL1 3 -mass [expr $Dens3/$g] -integration Lobatto
element dispBeamColumn 3003 4 8 4 $DL1 4 -mass [expr $Dens4/$g] -integration Lobatto
element dispBeamColumn 3004 5 6 4 $DL1 5 -mass [expr $Dens5/$g] -integration Lobatto
element dispBeamColumn 3005 6 7 4 $DL1 6 -mass [expr $Dens6/$g] -integration Lobatto
element dispBeamColumn 3006 8 7 4 $DL1 7 -mass [expr $Dens7/$g] -integration Lobatto
element dispBeamColumn 3007 5 8 4 $DL1 8 -mass [expr $Dens8/$g] -integration Lobatto
element dispBeamColumn 3008 4 5 4 $DL1 9 -mass [expr $Dens9/$g] -integration Lobatto
puts " TIMESERIES... "
set timeLinear 1
timeSeries Linear $timeLinear
puts " LOAD PATTERN... "
pattern Plain 1 $timeLinear {
#(- - +) original
eleLoad -ele 3000 -type -beamUniform 0.0 0.0 -46.85118385355202
eleLoad -ele 3001 -type -beamUniform 0.0 0.0 -46.85118385355202
eleLoad -ele 3002 -type -beamUniform 0.0 0.0 -46.85118385355202
eleLoad -ele 3003 -type -beamUniform 0.0 0.0 -46.85118385355202
eleLoad -ele 3004 -type -beamUniform 0.0 -46.85118385355202 0.0
eleLoad -ele 3005 -type -beamUniform 46.85118385355202 0.0 0.0
eleLoad -ele 3006 -type -beamUniform 0.0 -46.85118385355202 0.0
eleLoad -ele 3007 -type -beamUniform 46.85118385355202 0.0 0.0
eleLoad -ele 3008 -type -beamUniform 33.12878980946432 0.0 33.12878980946432
}
#
puts " STATIC OUTPUT...";
recorder Node -file ReaccionBase1.out -time -node 1 -dof 1 2 3 4 5 6 reaction -timeSeries 1;
recorder Node -file ReaccionBase2.out -time -node 2 -dof 1 2 3 4 5 6 reaction -timeSeries 1;
recorder Node -file ReaccionBase3.out -time -node 3 -dof 1 2 3 4 5 6 reaction -timeSeries 1;
recorder Node -file ReaccionBase4.out -time -node 4 -dof 1 2 3 4 5 6 reaction -timeSeries 1;
puts "STATIC ANALYSIS...";
system UmfPack
numberer RCM
constraints Plain
test NormDispIncr 1.0e-5 10 0;
integrator LoadControl 0.1
algorithm Newton
analysis Static
analyze 10
loadConst -time 0.0



********************************************************************************************************************
#MODEL 3

puts " DIMENSIONS...."
wipe
model BasicBuilder -ndm 3 -ndf 6
set NumModos 1;
set pi 3.1416;
puts " NODES...."
node 1 0 0 0
node 2 3 0 0
node 3 3 0 3
node 4 0 0 3
node 5 0 3 0
node 6 3 3 0
node 7 3 3 3
node 8 0 3 3
puts " DOF...."
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
puts " GEOTRANSF..."
geomTransf PDelta 1 -3 0 0
geomTransf PDelta 2 -3 0 0
geomTransf PDelta 3 -3 0 0
geomTransf PDelta 4 -3 0 0
geomTransf PDelta 5 0 3 0
geomTransf PDelta 6 3 0 0
geomTransf PDelta 7 0 3 0
geomTransf PDelta 8 3 0 0
geomTransf PDelta 9 -3 0 0
geomTransf PDelta 10 -3 0 0
puts " MATERIAL..."
set E1 20389019158.0
set Fy1 35153481.31
uniaxialMaterial Steel02 1000 $Fy1 $E1 1e-9 20 0.925 0.15 0.0001 0.0001 0.0001 0.0001
#[expr {Dens_mat} * {Area_sec}]
set g 9.81
set Dens1 [expr 7849.05 * 0.005969026041820605]

puts " SECTION TAGS..."
set DL1 2000
puts " ELEMENTS..."
section Fiber $DL1 -GJ 423464.244 {
patch circ 1000 5 5 0.0 0.0 0.09 0.10 0.0 360.0
}
puts " ELEMENTS..."
element dispBeamColumn 3000 1 5 4 $DL1 1 -mass [expr $Dens1/$g] -integration Lobatto
element dispBeamColumn 3001 2 6 4 $DL1 2 -mass [expr $Dens1/$g] -integration Lobatto
element dispBeamColumn 3002 3 7 4 $DL1 3 -mass [expr $Dens1/$g] -integration Lobatto
element dispBeamColumn 3003 4 8 4 $DL1 4 -mass [expr $Dens1/$g] -integration Lobatto
element dispBeamColumn 3004 5 6 4 $DL1 5 -mass [expr $Dens1/$g] -integration Lobatto
element dispBeamColumn 3005 6 7 4 $DL1 6 -mass [expr $Dens1/$g] -integration Lobatto
element dispBeamColumn 3006 8 7 4 $DL1 7 -mass [expr $Dens1/$g] -integration Lobatto
element dispBeamColumn 3007 5 8 4 $DL1 8 -mass [expr $Dens1/$g] -integration Lobatto

element dispBeamColumn 3008 1 8 4 $DL1 9 -mass [expr $Dens1/$g] -integration Lobatto
element dispBeamColumn 3009 4 5 4 $DL1 9 -mass [expr $Dens1/$g] -integration Lobatto

puts " TIMESERIES... "
set timeLinear 1
timeSeries Linear $timeLinear
puts " LOAD PATTERN... "
pattern Plain 1 $timeLinear {
#(- - +) original
eleLoad -ele 3000 -type -beamUniform 0.0 0.0 -46.85118385355202
eleLoad -ele 3001 -type -beamUniform 0.0 0.0 -46.85118385355202
eleLoad -ele 3002 -type -beamUniform 0.0 0.0 -46.85118385355202
eleLoad -ele 3003 -type -beamUniform 0.0 0.0 -46.85118385355202
eleLoad -ele 3004 -type -beamUniform 0.0 -46.85118385355202 0.0
eleLoad -ele 3005 -type -beamUniform 46.85118385355202 0.0 0.0
eleLoad -ele 3006 -type -beamUniform 0.0 -46.85118385355202 0.0
eleLoad -ele 3007 -type -beamUniform 46.85118385355202 0.0 0.0

eleLoad -ele 3008 -type -beamUniform -33.12878980946432 0.0 -33.12878980946432
eleLoad -ele 3009 -type -beamUniform 33.12878980946432 0.0 -33.12878980946432
}
#
puts " STATIC OUTPUT...";
#records
recorder Node -file ReaccionBase1.out -time -node 1 -dof 1 2 3 4 5 6 reaction -timeSeries 1;
recorder Node -file ReaccionBase2.out -time -node 2 -dof 1 2 3 4 5 6 reaction -timeSeries 1;
recorder Node -file ReaccionBase3.out -time -node 3 -dof 1 2 3 4 5 6 reaction -timeSeries 1;
recorder Node -file ReaccionBase4.out -time -node 4 -dof 1 2 3 4 5 6 reaction -timeSeries 1;

puts "STATIC ANALYSIS...";
system UmfPack
numberer RCM
constraints Plain
test NormDispIncr 1.0e-5 10 0;
integrator LoadControl 0.1
algorithm Newton
analysis Static
analyze 10
loadConst -time 0.0

[sobeli]

Hi!,

I've just realize what's hapenning, in opensees is necessary to add an extra node to create de xbrace, I share this if is useful for someone.