Problems in static analysis with steel 3d model

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

Moderators: silvia, selimgunay, Moderators

Post Reply
sobeli
Posts: 6
Joined: Thu Jul 27, 2023 10:38 am

Problems in static analysis with steel 3d model

Post by sobeli »

Hi everyone!

I'm dealing with a steel structure model in 3d, all the sections are circular fibers, the mass, and the loads are along the elements (Also I tried with nodal options). I ran the static analysis and all my Forces ar correct, but my moments in X are wird (in nodes 1 and 2 are low and in node 3 is too high). I've checked my local axis and forces several times. Another thing is that all the structure is completely welded so I do not release the connections, in the intersections I subdivide them with a node to avoid numerical problems. I share my code...

puts " DIMENSIONS...."
wipe
model BasicBuilder -ndm 3 -ndf 6
set NumModos 1;
set pi 3.1416;
puts " NODES...."
node 1 0.000 0.000 0.000
node 2 18.050 0.000 0.000
node 3 9.025 0.000 15.632
node 4 0.263 2.106 0.228
node 5 4.463 2.106 0.228
node 6 9.025 2.106 0.228
node 7 13.587 2.106 0.228
node 8 17.787 2.106 0.228
node 9 2.364 2.106 3.865
node 10 6.766 2.106 3.825
node 11 11.263 2.106 3.817
node 12 15.687 2.106 3.865
node 13 4.568 2.106 7.684
node 14 9.025 2.106 7.806
node 15 13.482 2.106 7.685
node 16 6.924 2.106 11.767
node 17 11.125 2.106 11.767
node 18 9.025 2.106 15.404
node 19 9.025 13.096 1.418
node 20 12.719 13.096 7.816
node 21 5.331 13.096 7.816
node 22 2.638 21.107 2.285
node 23 9.025 21.107 2.285
node 24 15.412 21.107 2.286
node 25 5.831 21.107 7.816
node 26 12.219 21.107 7.816
node 27 9.025 21.107 13.348
node 28 9.025 28.749 3.113
node 29 11.741 28.749 7.816
node 30 6.310 28.750 7.816
node 31 4.301 34.402 3.725
node 32 9.025 34.402 3.725
node 33 13.749 34.402 3.725
node 34 6.663 34.402 7.816
node 35 11.387 34.402 7.816
node 36 9.025 34.402 11.907
node 37 4.396 35.164 3.808
node 38 13.654 35.164 3.808
node 39 9.025 35.164 11.826
node 40 13.597 35.621 3.857
node 41 4.453 35.621 3.857
node 42 9.025 35.621 11.776
node 43 4.453 40.206 3.857
node 44 13.597 40.206 3.857
node 45 9.025 40.206 11.776

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

puts " GEOTRANSF..."
geomTransf Pdelta 1 0.0 0.0 1.0
geomTransf Pdelta 2 0.0 0.0 1.0
geomTransf Pdelta 3 0.0 0.0 1.0
geomTransf Pdelta 4 0.0 0.0 1.0
geomTransf Pdelta 5 0.0 0.0 1.0
geomTransf Pdelta 6 0.0 0.0 1.0
geomTransf Pdelta 7 0.0 0.0 -1.0
geomTransf Pdelta 8 0.0 0.0 1.0
geomTransf Pdelta 9 0.0 0.0 -1.0
geomTransf Pdelta 10 0.0 0.0 1.0
geomTransf Pdelta 11 0.0 0.0 1.0
geomTransf Pdelta 12 0.0 0.0 -1.0
geomTransf Pdelta 13 0.0 0.0 1.0
geomTransf Pdelta 14 0.0 0.0 -1.0
geomTransf Pdelta 15 0.0 0.0 1.0
geomTransf Pdelta 16 0.0 0.0 -1.0
geomTransf Pdelta 17 0.0 0.0 1.0
geomTransf Pdelta 18 0.0 0.0 -1.0
geomTransf Pdelta 19 0.0 0.0 1.0
geomTransf Pdelta 20 0.0 0.0 1.0
geomTransf Pdelta 21 0.0 0.0 -1.0
geomTransf Pdelta 22 0.0 0.0 -1.0
geomTransf Pdelta 23 0.0 0.0 1.0
geomTransf Pdelta 24 0.0 0.0 -1.0
geomTransf Pdelta 25 0.0 0.0 1.0
geomTransf Pdelta 26 0.0 0.0 -1.0
geomTransf Pdelta 27 0.0 0.0 1.0
geomTransf Pdelta 28 0.0 0.0 -1.0
geomTransf Pdelta 29 0.0 0.0 1.0
geomTransf Pdelta 30 0.0 0.0 1.0
geomTransf Pdelta 31 0.0 0.0 1.0
geomTransf Pdelta 32 0.0 0.0 1.0
geomTransf Pdelta 33 0.0 0.0 1.0
geomTransf Pdelta 34 0.0 0.0 -1.0
geomTransf Pdelta 35 0.0 0.0 1.0
geomTransf Pdelta 36 0.0 0.0 -1.0
geomTransf Pdelta 37 0.0 0.0 -1.0
geomTransf Pdelta 38 0.0 0.0 1.0
geomTransf Pdelta 39 0.0 0.0 1.0
geomTransf Pdelta 40 0.0 0.0 1.0
geomTransf Pdelta 41 0.0 0.0 1.0
geomTransf Pdelta 42 0.0 0.0 1.0
geomTransf Pdelta 43 0.0 0.0 -1.0
geomTransf Pdelta 44 0.0 0.0 1.0
geomTransf Pdelta 45 0.0 0.0 -1.0
geomTransf Pdelta 46 0.0 0.0 -1.0
geomTransf Pdelta 47 0.0 0.0 1.0
geomTransf Pdelta 48 0.0 0.0 1.0
geomTransf Pdelta 49 0.0 0.0 1.0
geomTransf Pdelta 50 0.0 0.0 1.0
geomTransf Pdelta 51 0.0 0.0 -1.0
geomTransf Pdelta 52 -0.228 0.0 0.263
geomTransf Pdelta 53 -0.228 0.0 -0.263
geomTransf Pdelta 54 1.0 0.0 0.0
geomTransf Pdelta 55 -2.057 0.0 2.375
geomTransf Pdelta 56 -2.058 0.0 -2.374
geomTransf Pdelta 57 1.0 0.0 0.0
geomTransf Pdelta 58 -1.44 0.0 1.663
geomTransf Pdelta 59 -1.439 0.0 -1.663
geomTransf Pdelta 60 1.0 0.0 0.0
geomTransf Pdelta 61 -0.082 0.0 0.094
geomTransf Pdelta 62 -0.082 0.0 -0.095
geomTransf Pdelta 63 1.0 0.0 0.0
geomTransf Pdelta 64 -0.049 0.0 0.057
geomTransf Pdelta 65 -0.049 0.0 -0.057
geomTransf Pdelta 66 1.0 0.0 0.0
geomTransf Pdelta 67 -7.588 0.0 5.068
geomTransf Pdelta 68 -1.19 0.0 -8.761
geomTransf Pdelta 69 7.588 0.0 3.693
geomTransf Pdelta 70 -0.867 0.0 -6.387
geomTransf Pdelta 71 5.529 0.0 2.693
geomTransf Pdelta 72 -5.532 0.0 3.694
geomTransf Pdelta 73 -1.19 0.0 8.762
geomTransf Pdelta 74 -7.588 0.0 -5.068
geomTransf Pdelta 75 7.588 0.0 -3.694
geomTransf Pdelta 76 5.531 0.0 -2.693
geomTransf Pdelta 77 -0.868 0.0 6.387
geomTransf Pdelta 78 -5.532 0.0 -3.693
geomTransf Pdelta 79 -0.827 0.0 6.387
geomTransf Pdelta 80 -5.529 0.0 -3.671
geomTransf Pdelta 81 5.532 0.0 -2.715
geomTransf Pdelta 82 4.090 0.0 -2.008
geomTransf Pdelta 83 -0.612 0.0 4.724
geomTransf Pdelta 84 -4.091 0.0 -2.715
geomTransf Pdelta 85 -5.531 0.0 3.671
geomTransf Pdelta 86 -0.827 0.0 -6.387
geomTransf Pdelta 87 5.532 0.0 2.715
geomTransf Pdelta 88 -0.612 0.0 -4.724
geomTransf Pdelta 89 4.090 0.0 2.008
geomTransf Pdelta 90 -4.091 0.0 2.715

puts " MATERIAL..."
set E1 2.039e10
set Fy1 2.95e7
uniaxialMaterial Steel02 1000 $Fy1 $E1 1e-9 20 0.925 0.15 0.0001 0.0001 0.0001 0.0001
set E2 2.039e10
set Fy2 4.05e7
uniaxialMaterial Steel02 1001 $Fy2 $E2 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.1045]
set Dens2 [expr 7849.05 * 0.1045]
set Dens3 [expr 7849.05 * 0.1045]
set Dens4 [expr 7849.05 * 0.0352]
set Dens5 [expr 7849.05 * 0.0352]
set Dens6 [expr 7849.05 * 0.0352]
set Dens7 [expr 7849.05 * 0.0352]
set Dens8 [expr 7849.05 * 0.0352]
set Dens9 [expr 7849.05 * 0.0352]
set Dens10 [expr 7849.05 * 0.035]
set Dens11 [expr 7849.05 * 0.035]
set Dens12 [expr 7849.05 * 0.035]
set Dens13 [expr 7849.05 * 0.035]
set Dens14 [expr 7849.05 * 0.035]
set Dens15 [expr 7849.05 * 0.035]
set Dens16 [expr 7849.05 * 0.020]
set Dens17 [expr 7849.05 * 0.020]
set Dens18 [expr 7849.05 * 0.020]
set Dens19 [expr 7849.05 * 0.020]
set Dens20 [expr 7849.05 * 0.020]
set Dens21 [expr 7849.05 * 0.020]
set Dens22 [expr 7849.05 * 0.014]
set Dens23 [expr 7849.05 * 0.014]
set Dens24 [expr 7849.05 * 0.014]
set Dens25 [expr 7849.05 * 0.014]
set Dens26 [expr 7849.05 * 0.014]
set Dens27 [expr 7849.05 * 0.014]
set Dens28 [expr 7849.05 * 0.014]
set Dens29 [expr 7849.05 * 0.014]
set Dens30 [expr 7849.05 * 0.014]
set Dens31 [expr 7849.05 * 0.020]
set Dens32 [expr 7849.05 * 0.020]
set Dens33 [expr 7849.05 * 0.020]
set Dens34 [expr 7849.05 * 0.020]
set Dens35 [expr 7849.05 * 0.020]
set Dens36 [expr 7849.05 * 0.020]
set Dens37 [expr 7849.05 * 0.016]
set Dens38 [expr 7849.05 * 0.016]
set Dens39 [expr 7849.05 * 0.016]
set Dens40 [expr 7849.05 * 0.018]
set Dens41 [expr 7849.05 * 0.018]
set Dens42 [expr 7849.05 * 0.018]
set Dens43 [expr 7849.05 * 0.018]
set Dens44 [expr 7849.05 * 0.018]
set Dens45 [expr 7849.05 * 0.018]
set Dens46 [expr 7849.05 * 0.016]
set Dens47 [expr 7849.05 * 0.016]
set Dens48 [expr 7849.05 * 0.016]
set Dens49 [expr 7849.05 * 0.016]
set Dens50 [expr 7849.05 * 0.016]
set Dens51 [expr 7849.05 * 0.016]
set Dens52 [expr 7849.05 * 0.100]
set Dens53 [expr 7849.05 * 0.100]
set Dens54 [expr 7849.05 * 0.050]
set Dens55 [expr 7849.05 * 0.050]
set Dens56 [expr 7849.05 * 0.050]
set Dens57 [expr 7849.05 * 0.050]
set Dens58 [expr 7849.05 * 0.100]
set Dens59 [expr 7849.05 * 0.100]
set Dens60 [expr 7849.05 * 0.100]
set Dens61 [expr 7849.05 * 0.100]
set Dens62 [expr 7849.05 * 0.100]
set Dens63 [expr 7849.05 * 0.100]
set Dens64 [expr 7849.05 * 0.104]
set Dens65 [expr 7849.05 * 0.104]
set Dens66 [expr 7849.05 * 0.104]
set Dens67 [expr 7849.05 * 0.035]
set Dens68 [expr 7849.05 * 0.035]
set Dens69 [expr 7849.05 * 0.035]
set Dens70 [expr 7849.05 * 0.035]
set Dens71 [expr 7849.05 * 0.035]
set Dens72 [expr 7849.05 * 0.035]
set Dens73 [expr 7849.05 * 0.035]
set Dens74 [expr 7849.05 * 0.035]
set Dens75 [expr 7849.05 * 0.035]
set Dens76 [expr 7849.05 * 0.038]
set Dens77 [expr 7849.05 * 0.038]
set Dens78 [expr 7849.05 * 0.038]
set Dens79 [expr 7849.05 * 0.029]
set Dens80 [expr 7849.05 * 0.029]
set Dens81 [expr 7849.05 * 0.029]
set Dens82 [expr 7849.05 * 0.029]
set Dens83 [expr 7849.05 * 0.029]
set Dens84 [expr 7849.05 * 0.029]
set Dens85 [expr 7849.05 * 0.029]
set Dens86 [expr 7849.05 * 0.029]
set Dens87 [expr 7849.05 * 0.029]
set Dens88 [expr 7849.05 * 0.029]
set Dens89 [expr 7849.05 * 0.029]
set Dens90 [expr 7849.05 * 0.029]

puts " SECTION TAGS..."

set DL1 2000
set H10 2001
set H11 2002
set H12 2003
set H14 2004
set H15 2005
set H16 2006
set H20 2007
set H21 2008
set H30 2009
set H31 2010
set H40 2011
set LG1 2012
set LG3 2013
set LG4 2014
set PL4 2015
set V10 2016
set V11 2017
set V15 2018
set V20 2019
set V21 2020

puts " ELEMENTS..."

section Fiber $DL1 -GJ 1.58e8 {
patch circ 1000 3 8 0.0 0.0 0.419 0.457 0.0 360.0
}
#
section Fiber $H10 -GJ 2.42e7 {
patch circ 1000 3 8 0.0 0.0 0.286 0.305 0.0 360.0
}
#
section Fiber $H11 -GJ 2.42e7 {
patch circ 1000 3 8 0.0 0.0 0.286 0.305 0.0 360.0
}
#
section Fiber $H12 -GJ 2.42e7 {
patch circ 1000 3 8 0.0 0.0 0.286 0.305 0.0 360.0
}
#
section Fiber $H14 -GJ 2.42e7 {
patch circ 1000 3 8 0.0 0.0 0.286 0.305 0.0 360.0
}
#
section Fiber $H15 -GJ 9.72e6 {
patch circ 1000 3 8 0.0 0.0 0.241 0.254 0.0 360.0
}
#
section Fiber $H16 -GJ 3.24e6 {
patch circ 1000 3 8 0.0 0.0 0.164 0.178 0.0 360.0
}
#
section Fiber $H20 -GJ 9.72e6 {
patch circ 1000 3 8 0.0 0.0 0.241 0.254 0.0 360.0
}
#
section Fiber $H21 -GJ 4.87e6 {
patch circ 1000 3 8 0.0 0.0 0.19 0.203 0.0 360.0
}
#
section Fiber $H30 -GJ 7.01e6 {
patch circ 1000 3 8 0.0 0.0 0.215 0.228 0.0 360.0
}
#
section Fiber $H31 -GJ 4.87e6 {
patch circ 1000 3 8 0.0 0.0 0.19 0.203 0.0 360.0
}
#
section Fiber $H40 -GJ 4.87e6 {
patch circ 1000 3 8 0.0 0.0 0.19 0.203 0.0 360.0
}
#
section Fiber $LG1 -GJ 1.96e8 {
patch circ 1001 3 8 0.0 0.0 0.482 0.5145 0.0 360.0
}
#
section Fiber $LG3 -GJ 1.02e8 {
patch circ 1000 3 8 0.0 0.0 0.498 0.5145 0.0 360.0
}
#
section Fiber $LG4 -GJ 1.96e8 {
patch circ 1001 3 8 0.0 0.0 0.482 0.5145 0.0 360.0
}
#
section Fiber $PL4 -GJ 1.58e8 {
patch circ 1000 3 8 0.0 0.0 0.419 0.457 0.0 360.0
}
#
section Fiber $V10 -GJ 2.42e7 {
patch circ 1000 3 8 0.0 0.0 0.286 0.305 0.0 360.0
}
#
section Fiber $V11 -GJ 2.42e7 {
patch circ 1000 3 8 0.0 0.0 0.286 0.305 0.0 360.0
}
#
section Fiber $V15 -GJ 2.65e7 {
patch circ 1000 3 8 0.0 0.0 0.284 0.305 0.0 360.0
}
#
section Fiber $V20 -GJ 1.37e7 {
patch circ 1000 3 8 0.0 0.0 0.235 0.254 0.0 360.0
}
#
section Fiber $V21 -GJ 1.37e7 {
patch circ 1000 3 8 0.0 0.0 0.235 0.254 0.0 360.0
}
#

puts " ELEMENTS..."

element dispBeamColumn 3000 40 44 10 $DL1 1 -mass [expr $Dens1/$g] -integration NewtonCotes
element dispBeamColumn 3001 41 43 10 $DL1 2 -mass [expr $Dens2/$g] -integration NewtonCotes
element dispBeamColumn 3002 42 45 10 $DL1 3 -mass [expr $Dens3/$g] -integration NewtonCotes
element dispBeamColumn 3003 4 9 10 $H10 4 -mass [expr $Dens4/$g] -integration NewtonCotes
element dispBeamColumn 3004 5 6 10 $H10 5 -mass [expr $Dens5/$g] -integration NewtonCotes
element dispBeamColumn 3005 6 7 10 $H10 6 -mass [expr $Dens6/$g] -integration NewtonCotes
element dispBeamColumn 3006 8 12 10 $H10 7 -mass [expr $Dens7/$g] -integration NewtonCotes
element dispBeamColumn 3007 9 13 10 $H10 8 -mass [expr $Dens8/$g] -integration NewtonCotes
element dispBeamColumn 3008 12 15 10 $H10 9 -mass [expr $Dens9/$g] -integration NewtonCotes
element dispBeamColumn 3009 4 5 10 $H11 10 -mass [expr $Dens10/$g] -integration NewtonCotes
element dispBeamColumn 3010 7 8 10 $H11 11 -mass [expr $Dens11/$g] -integration NewtonCotes
element dispBeamColumn 3011 17 18 10 $H12 12 -mass [expr $Dens12/$g] -integration NewtonCotes
element dispBeamColumn 3012 13 16 10 $H14 13 -mass [expr $Dens13/$g] -integration NewtonCotes
element dispBeamColumn 3013 15 17 10 $H14 14 -mass [expr $Dens14/$g] -integration NewtonCotes
element dispBeamColumn 3014 16 18 10 $H14 15 -mass [expr $Dens15/$g] -integration NewtonCotes
element dispBeamColumn 3015 6 10 10 $H15 16 -mass [expr $Dens16/$g] -integration NewtonCotes
element dispBeamColumn 3016 6 11 10 $H15 17 -mass [expr $Dens17/$g] -integration NewtonCotes
element dispBeamColumn 3017 10 13 10 $H15 18 -mass [expr $Dens18/$g] -integration NewtonCotes
element dispBeamColumn 3018 11 15 10 $H15 19 -mass [expr $Dens19/$g] -integration NewtonCotes
element dispBeamColumn 3019 13 14 10 $H15 20 -mass [expr $Dens20/$g] -integration NewtonCotes
element dispBeamColumn 3020 15 14 10 $H15 21 -mass [expr $Dens21/$g] -integration NewtonCotes
element dispBeamColumn 3021 5 9 10 $H16 22 -mass [expr $Dens22/$g] -integration NewtonCotes
element dispBeamColumn 3022 5 10 10 $H16 23 -mass [expr $Dens23/$g] -integration NewtonCotes
element dispBeamColumn 3023 7 11 10 $H16 24 -mass [expr $Dens24/$g] -integration NewtonCotes
element dispBeamColumn 3024 7 12 10 $H16 25 -mass [expr $Dens25/$g] -integration NewtonCotes
element dispBeamColumn 3025 10 9 10 $H16 26 -mass [expr $Dens26/$g] -integration NewtonCotes
element dispBeamColumn 3026 11 12 10 $H16 27 -mass [expr $Dens27/$g] -integration NewtonCotes
element dispBeamColumn 3027 14 16 10 $H16 28 -mass [expr $Dens28/$g] -integration NewtonCotes
element dispBeamColumn 3028 14 17 10 $H16 29 -mass [expr $Dens29/$g] -integration NewtonCotes
element dispBeamColumn 3029 16 17 10 $H16 30 -mass [expr $Dens30/$g] -integration NewtonCotes
element dispBeamColumn 3030 22 23 10 $H20 31 -mass [expr $Dens31/$g] -integration NewtonCotes
element dispBeamColumn 3031 22 25 10 $H20 32 -mass [expr $Dens32/$g] -integration NewtonCotes
element dispBeamColumn 3032 23 24 10 $H20 33 -mass [expr $Dens33/$g] -integration NewtonCotes
element dispBeamColumn 3033 24 26 10 $H20 34 -mass [expr $Dens34/$g] -integration NewtonCotes
element dispBeamColumn 3034 25 27 10 $H20 35 -mass [expr $Dens35/$g] -integration NewtonCotes
element dispBeamColumn 3035 26 27 10 $H20 36 -mass [expr $Dens36/$g] -integration NewtonCotes
element dispBeamColumn 3036 23 25 10 $H21 37 -mass [expr $Dens37/$g] -integration NewtonCotes
element dispBeamColumn 3037 23 26 10 $H21 38 -mass [expr $Dens38/$g] -integration NewtonCotes
element dispBeamColumn 3038 25 26 10 $H21 39 -mass [expr $Dens39/$g] -integration NewtonCotes
element dispBeamColumn 3039 31 32 10 $H30 40 -mass [expr $Dens40/$g] -integration NewtonCotes
element dispBeamColumn 3040 31 34 10 $H30 41 -mass [expr $Dens41/$g] -integration NewtonCotes
element dispBeamColumn 3041 32 33 10 $H30 42 -mass [expr $Dens42/$g] -integration NewtonCotes
element dispBeamColumn 3042 33 35 10 $H30 43 -mass [expr $Dens43/$g] -integration NewtonCotes
element dispBeamColumn 3043 34 36 10 $H30 44 -mass [expr $Dens44/$g] -integration NewtonCotes
element dispBeamColumn 3044 35 36 10 $H30 45 -mass [expr $Dens45/$g] -integration NewtonCotes
element dispBeamColumn 3045 32 34 10 $H31 46 -mass [expr $Dens46/$g] -integration NewtonCotes
element dispBeamColumn 3046 32 35 10 $H31 47 -mass [expr $Dens47/$g] -integration NewtonCotes
element dispBeamColumn 3047 34 35 10 $H31 48 -mass [expr $Dens48/$g] -integration NewtonCotes
element dispBeamColumn 3048 43 44 10 $H40 49 -mass [expr $Dens49/$g] -integration NewtonCotes
element dispBeamColumn 3049 43 45 10 $H40 50 -mass [expr $Dens50/$g] -integration NewtonCotes
element dispBeamColumn 3050 44 45 10 $H40 51 -mass [expr $Dens51/$g] -integration NewtonCotes
element dispBeamColumn 3051 1 4 10 $LG1 52 -mass [expr $Dens52/$g] -integration NewtonCotes
element dispBeamColumn 3052 2 8 10 $LG1 53 -mass [expr $Dens53/$g] -integration NewtonCotes
element dispBeamColumn 3053 3 18 10 $LG3 54 -mass [expr $Dens54/$g] -integration NewtonCotes
element dispBeamColumn 3054 4 22 10 $LG3 55 -mass [expr $Dens55/$g] -integration NewtonCotes
element dispBeamColumn 3055 8 24 10 $LG3 56 -mass [expr $Dens56/$g] -integration NewtonCotes
element dispBeamColumn 3056 18 27 10 $LG3 57 -mass [expr $Dens57/$g] -integration NewtonCotes
element dispBeamColumn 3057 22 31 10 $LG4 58 -mass [expr $Dens58/$g] -integration NewtonCotes
element dispBeamColumn 3058 24 33 10 $LG4 59 -mass [expr $Dens59/$g] -integration NewtonCotes
element dispBeamColumn 3059 27 36 10 $LG4 60 -mass [expr $Dens60/$g] -integration NewtonCotes
element dispBeamColumn 3060 31 37 10 $LG4 61 -mass [expr $Dens61/$g] -integration NewtonCotes
element dispBeamColumn 3061 33 38 10 $LG4 62 -mass [expr $Dens62/$g] -integration NewtonCotes
element dispBeamColumn 3062 36 39 10 $LG4 63 -mass [expr $Dens63/$g] -integration NewtonCotes
element dispBeamColumn 3063 37 41 10 $PL4 64 -mass [expr $Dens64/$g] -integration NewtonCotes
element dispBeamColumn 3064 38 40 10 $PL4 65 -mass [expr $Dens65/$g] -integration NewtonCotes
element dispBeamColumn 3065 39 42 10 $PL4 66 -mass [expr $Dens66/$g] -integration NewtonCotes
element dispBeamColumn 3066 4 21 10 $V10 67 -mass [expr $Dens67/$g] -integration NewtonCotes
element dispBeamColumn 3067 8 19 10 $V10 68 -mass [expr $Dens68/$g] -integration NewtonCotes
element dispBeamColumn 3068 18 20 10 $V10 69 -mass [expr $Dens69/$g] -integration NewtonCotes
element dispBeamColumn 3069 19 22 10 $V10 70 -mass [expr $Dens70/$g] -integration NewtonCotes
element dispBeamColumn 3070 20 24 10 $V10 71 -mass [expr $Dens71/$g] -integration NewtonCotes
element dispBeamColumn 3071 21 27 10 $V10 72 -mass [expr $Dens72/$g] -integration NewtonCotes
element dispBeamColumn 3072 4 19 10 $V11 73 -mass [expr $Dens73/$g] -integration NewtonCotes
element dispBeamColumn 3073 8 20 10 $V11 74 -mass [expr $Dens74/$g] -integration NewtonCotes
element dispBeamColumn 3074 18 21 10 $V11 75 -mass [expr $Dens75/$g] -integration NewtonCotes
element dispBeamColumn 3075 21 22 10 $V15 76 -mass [expr $Dens76/$g] -integration NewtonCotes
element dispBeamColumn 3076 19 24 10 $V15 77 -mass [expr $Dens77/$g] -integration NewtonCotes
element dispBeamColumn 3077 20 27 10 $V15 78 -mass [expr $Dens78/$g] -integration NewtonCotes
element dispBeamColumn 3078 22 28 10 $V20 79 -mass [expr $Dens79/$g] -integration NewtonCotes
element dispBeamColumn 3079 24 29 10 $V20 80 -mass [expr $Dens80/$g] -integration NewtonCotes
element dispBeamColumn 3080 27 30 10 $V20 81 -mass [expr $Dens81/$g] -integration NewtonCotes
element dispBeamColumn 3081 30 31 10 $V20 82 -mass [expr $Dens82/$g] -integration NewtonCotes
element dispBeamColumn 3082 28 33 10 $V20 83 -mass [expr $Dens83/$g] -integration NewtonCotes
element dispBeamColumn 3083 29 36 10 $V20 84 -mass [expr $Dens84/$g] -integration NewtonCotes
element dispBeamColumn 3084 22 30 10 $V21 85 -mass [expr $Dens85/$g] -integration NewtonCotes
element dispBeamColumn 3085 24 28 10 $V21 86 -mass [expr $Dens86/$g] -integration NewtonCotes
element dispBeamColumn 3086 27 29 10 $V21 87 -mass [expr $Dens87/$g] -integration NewtonCotes
element dispBeamColumn 3087 28 31 10 $V21 88 -mass [expr $Dens88/$g] -integration NewtonCotes
element dispBeamColumn 3088 29 33 10 $V21 89 -mass [expr $Dens89/$g] -integration NewtonCotes
element dispBeamColumn 3089 30 36 10 $V21 90 -mass [expr $Dens90/$g] -integration NewtonCotes

puts " TIMESERIES... "
set timeLinear 1
timeSeries Linear $timeLinear
puts " LOAD PATTERN... "
pattern Plain 1 $timeLinear {
eleLoad -ele 3000 -type -beamUniform 0.0 0.0 -820.832
eleLoad -ele 3001 -type -beamUniform 0.0 0.0 -820.832
eleLoad -ele 3002 -type -beamUniform 0.0 0.0 -820.832
eleLoad -ele 3003 -type -beamUniform -276.890 0.0 0.0
eleLoad -ele 3004 -type -beamUniform -276.890 0.0 0.0
eleLoad -ele 3005 -type -beamUniform -276.890 0.0 0.0
eleLoad -ele 3006 -type -beamUniform -276.890 0.0 0.0
eleLoad -ele 3007 -type -beamUniform -276.890 0.0 0.0
eleLoad -ele 3008 -type -beamUniform -276.890 0.0 0.0
eleLoad -ele 3009 -type -beamUniform -276.890 0.0 0.0
eleLoad -ele 3010 -type -beamUniform -276.890 0.0 0.0
eleLoad -ele 3011 -type -beamUniform -276.890 0.0 0.0
eleLoad -ele 3012 -type -beamUniform -276.890 0.0 0.0
eleLoad -ele 3013 -type -beamUniform -276.890 0.0 0.0
eleLoad -ele 3014 -type -beamUniform -276.890 0.0 0.0
eleLoad -ele 3015 -type -beamUniform -158.677 0.0 0.0
eleLoad -ele 3016 -type -beamUniform -158.677 0.0 0.0
eleLoad -ele 3017 -type -beamUniform -158.677 0.0 0.0
eleLoad -ele 3018 -type -beamUniform -158.677 0.0 0.0
eleLoad -ele 3019 -type -beamUniform -158.677 0.0 0.0
eleLoad -ele 3020 -type -beamUniform -158.677 0.0 0.0
eleLoad -ele 3021 -type -beamUniform -109.952 0.0 0.0
eleLoad -ele 3022 -type -beamUniform -109.952 0.0 0.0
eleLoad -ele 3023 -type -beamUniform -109.952 0.0 0.0
eleLoad -ele 3024 -type -beamUniform -109.952 0.0 0.0
eleLoad -ele 3025 -type -beamUniform -109.952 0.0 0.0
eleLoad -ele 3026 -type -beamUniform -109.952 0.0 0.0
eleLoad -ele 3027 -type -beamUniform -109.952 0.0 0.0
eleLoad -ele 3028 -type -beamUniform -109.952 0.0 0.0
eleLoad -ele 3029 -type -beamUniform -109.952 0.0 0.0
eleLoad -ele 3030 -type -beamUniform -158.677 0.0 0.0
eleLoad -ele 3031 -type -beamUniform -158.677 0.0 0.0
eleLoad -ele 3032 -type -beamUniform -158.677 0.0 0.0
eleLoad -ele 3033 -type -beamUniform -158.677 0.0 0.0
eleLoad -ele 3034 -type -beamUniform -158.677 0.0 0.0
eleLoad -ele 3035 -type -beamUniform -158.677 0.0 0.0
eleLoad -ele 3036 -type -beamUniform -125.980 0.0 0.0
eleLoad -ele 3037 -type -beamUniform -125.980 0.0 0.0
eleLoad -ele 3038 -type -beamUniform -125.980 0.0 0.0
eleLoad -ele 3039 -type -beamUniform -142.328 0.0 0.0
eleLoad -ele 3040 -type -beamUniform -142.328 0.0 0.0
eleLoad -ele 3041 -type -beamUniform -142.328 0.0 0.0
eleLoad -ele 3042 -type -beamUniform -142.328 0.0 0.0
eleLoad -ele 3043 -type -beamUniform -142.328 0.0 0.0
eleLoad -ele 3044 -type -beamUniform -142.328 0.0 0.0
eleLoad -ele 3045 -type -beamUniform -125.980 0.0 0.0
eleLoad -ele 3046 -type -beamUniform -125.980 0.0 0.0
eleLoad -ele 3047 -type -beamUniform -125.980 0.0 0.0
eleLoad -ele 3048 -type -beamUniform -125.980 0.0 0.0
eleLoad -ele 3049 -type -beamUniform -125.980 0.0 0.0
eleLoad -ele 3050 -type -beamUniform -125.980 0.0 0.0
eleLoad -ele 3051 -type -beamUniform -128.282 0.0 -776.175
eleLoad -ele 3052 -type -beamUniform -128.282 0.0 -776.175
eleLoad -ele 3053 -type -beamUniform -43.0172 0.0 -397.343
eleLoad -ele 3054 -type -beamUniform -65.2021 0.0 -394.310
eleLoad -ele 3055 -type -beamUniform -65.2153 0.0 -394.308
eleLoad -ele 3056 -type -beamUniform -42.9947 0.0 -397.345
eleLoad -ele 3057 -type -beamUniform -128.423 0.0 -776.152
eleLoad -ele 3058 -type -beamUniform -128.386 0.0 -776.158
eleLoad -ele 3059 -type -beamUniform -84.7718 0.0 -782.124
eleLoad -ele 3060 -type -beamUniform -128.491 0.0 -776.141
eleLoad -ele 3061 -type -beamUniform -128.491 0.0 -776.141
eleLoad -ele 3062 -type -beamUniform -83.1576 0.0 -782.297
eleLoad -ele 3063 -type -beamUniform -133.219 0.0 -809.950
eleLoad -ele 3064 -type -beamUniform -133.219 0.0 -809.950
eleLoad -ele 3065 -type -beamUniform -89.2739 0.0 -815.963
eleLoad -ele 3066 -type -beamUniform -176.877 0.0 -213.032
eleLoad -ele 3067 -type -beamUniform -173.575 0.0 -215.731
eleLoad -ele 3068 -type -beamUniform -168.641 0.0 -219.609
eleLoad -ele 3069 -type -beamUniform -173.575 0.0 -215.731
eleLoad -ele 3070 -type -beamUniform -168.625 0.0 -219.621
eleLoad -ele 3071 -type -beamUniform -176.886 0.0 -213.024
eleLoad -ele 3072 -type -beamUniform -173.575 0.0 -215.731
eleLoad -ele 3073 -type -beamUniform -176.877 0.0 -213.032
eleLoad -ele 3074 -type -beamUniform -168.641 0.0 -219.609
eleLoad -ele 3075 -type -beamUniform -185.762 0.0 -241.905
eleLoad -ele 3076 -type -beamUniform -191.199 0.0 -237.630
eleLoad -ele 3077 -type -beamUniform -194.844 0.0 -234.651
eleLoad -ele 3078 -type -beamUniform -147.641 0.0 -175.185
eleLoad -ele 3079 -type -beamUniform -150.233 0.0 -172.967
eleLoad -ele 3080 -type -beamUniform -143.799 0.0 -178.351
eleLoad -ele 3081 -type -beamUniform -143.811 0.0 -178.342
eleLoad -ele 3082 -type -beamUniform -147.628 0.0 -175.196
eleLoad -ele 3083 -type -beamUniform -150.242 0.0 -172.959
eleLoad -ele 3084 -type -beamUniform -150.239 0.0 -172.961
eleLoad -ele 3085 -type -beamUniform -147.639 0.0 -175.186
eleLoad -ele 3086 -type -beamUniform -143.817 0.0 -178.337
eleLoad -ele 3087 -type -beamUniform -147.628 0.0 -175.196
eleLoad -ele 3088 -type -beamUniform -143.788 0.0 -178.361
eleLoad -ele 3089 -type -beamUniform -150.247 0.0 -172.955
}
#
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;


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

Return to “OpenSees.exe Users”