Hi there,
I am doing a grillage analysis of a bridge where I have longitudinal and transverse members. I have only loaded the longitudinal members. When I use the globalForce element recorder I understand that the output is in the form Fx Fy Fz Mx My Mz for the i and j nodes of the element. When I output results for a few elements on one of the longitudinal beams, am I wrong in expecting that the moment(say My) at the j node of one element should be the same as the i node of the following element, as they are effectively the same node? My code is below(sorry there are loads of elements)
My global axis is x to the right, z up, y into the screen(when you are looking at the screen)
thanks
#Grillage model of a Two Span right slab Bridge
#10/02/11
#
#units N, m, kg
#
#
#SET UP ----------------------------------------------------------------------------
wipe; # clear memory of all past model definitions
model BasicBuilder -ndm 3 -ndf 6; # Define the model builder, ndm=#dimension, ndf=#dofs
set dataDir DataBridge3DwithColumns; # set up name of data directory (you can remove this)
file mkdir $dataDir; # create data directory
set GMdir "GMfiles"; # ground-motion file directory
#DEFINE GEOMETRY---------------------------------------------------------------------
# X Y Z
#Deck----------------
node 1 0. 0. 0.;
node 2 1.5 0. 0.;
node 3 3. 0. 0.;
node 4 4.5 0. 0.;
node 5 6. 0. 0.;
node 6 7.5 0. 0.;
node 7 9. 0. 0.;
node 8 10.5 0. 0.;
node 9 12. 0. 0.;
node 10 13.5 0. 0.;
node 11 15. 0. 0.;
node 12 16.5 0. 0.;
node 13 18. 0. 0.;
node 14 19.5 0. 0.;
node 15 21. 0. 0.;
node 16 22.5 0. 0.;
node 17 24. 0. 0.;
node 18 25.5 0. 0.;
node 19 27. 0. 0.;
node 20 28.5 0. 0.;
node 21 30. 0. 0.;
node 22 31.5 0. 0.;
node 23 33. 0. 0.;
node 24 34.5 0. 0.;
node 25 36. 0. 0.;
node 26 0. 1.26 0.;
node 27 1.5 1.26 0.;
node 28 3. 1.26 0.;
node 29 4.5 1.26 0.;
node 30 6. 1.26 0.;
node 31 7.5 1.26 0.;
node 32 9. 1.26 0.;
node 33 10.5 1.26 0.;
node 34 12. 1.26 0.;
node 35 13.5 1.26 0.;
node 36 15. 1.26 0.;
node 37 16.5 1.26 0.;
node 38 18. 1.26 0.;
node 39 19.5 1.26 0.;
node 40 21. 1.26 0.;
node 41 22.5 1.26 0.;
node 42 24. 1.26 0.;
node 43 25.5 1.26 0.;
node 44 27. 1.26 0.;
node 45 28.5 1.26 0.;
node 46 30. 1.26 0.;
node 47 31.5 1.26 0.;
node 48 33. 1.26 0.;
node 49 34.5 1.26 0.;
node 50 36. 1.26 0.;
node 51 0. 2.26 0.;
node 52 1.5 2.26 0.;
node 53 3. 2.26 0.;
node 54 4.5 2.26 0.;
node 55 6. 2.26 0.;
node 56 7.5 2.26 0.;
node 57 9. 2.26 0.;
node 58 10.5 2.26 0.;
node 59 12. 2.26 0.;
node 60 13.5 2.26 0.;
node 61 15. 2.26 0.;
node 62 16.5 2.26 0.;
node 63 18. 2.26 0.;
node 64 19.5 2.26 0.;
node 65 21. 2.26 0.;
node 66 22.5 2.26 0.;
node 67 24. 2.26 0.;
node 68 25.5 2.26 0.;
node 69 27. 2.26 0.;
node 70 28.5 2.26 0.;
node 71 30. 2.26 0.;
node 72 31.5 2.26 0.;
node 73 33. 2.26 0.;
node 74 34.5 2.26 0.;
node 75 36. 2.26 0.;
node 76 0. 3.46 0.;
node 77 1.5 3.46 0.;
node 78 3. 3.46 0.;
node 79 4.5 3.46 0.;
node 80 6. 3.46 0.;
node 81 7.5 3.46 0.;
node 82 9. 3.46 0.;
node 83 10.5 3.46 0.;
node 84 12. 3.46 0.;
node 85 13.5 3.46 0.;
node 86 15. 3.46 0.;
node 87 16.5 3.46 0.;
node 88 18. 3.46 0.;
node 89 19.5 3.46 0.;
node 90 21. 3.46 0.;
node 91 22.5 3.46 0.;
node 92 24. 3.46 0.;
node 93 25.5 3.46 0.;
node 94 27. 3.46 0.;
node 95 28.5 3.46 0.;
node 96 30. 3.46 0.;
node 97 31.5 3.46 0.;
node 98 33. 3.46 0.;
node 99 34.5 3.46 0.;
node 100 36. 3.46 0.;
node 101 0. 4.56 0.;
node 102 1.5 4.56 0.;
node 103 3. 4.56 0.;
node 104 4.5 4.56 0.;
node 105 6. 4.56 0.;
node 106 7.5 4.56 0.;
node 107 9. 4.56 0.;
node 108 10.5 4.56 0.;
node 109 12. 4.56 0.;
node 110 13.5 4.56 0.;
node 111 15. 4.56 0.;
node 112 16.5 4.56 0.;
node 113 18. 4.56 0.;
node 114 19.5 4.56 0.;
node 115 21. 4.56 0.;
node 116 22.5 4.56 0.;
node 117 24. 4.56 0.;
node 118 25.5 4.56 0.;
node 119 27. 4.56 0.;
node 120 28.5 4.56 0.;
node 121 30. 4.56 0.;
node 122 31.5 4.56 0.;
node 123 33. 4.56 0.;
node 124 34.5 4.56 0.;
node 125 36. 4.56 0.;
node 126 0. 5.66 0.;
node 127 1.5 5.66 0.;
node 128 3. 5.66 0.;
node 129 4.5 5.66 0.;
node 130 6. 5.66 0.;
node 131 7.5 5.66 0.;
node 132 9. 5.66 0.;
node 133 10.5 5.66 0.;
node 134 12. 5.66 0.;
node 135 13.5 5.66 0.;
node 136 15. 5.66 0.;
node 137 16.5 5.66 0.;
node 138 18. 5.66 0.;
node 139 19.5 5.66 0.;
node 140 21. 5.66 0.;
node 141 22.5 5.66 0.;
node 142 24. 5.66 0.;
node 143 25.5 5.66 0.;
node 144 27. 5.66 0.;
node 145 28.5 5.66 0.;
node 146 30. 5.66 0.;
node 147 31.5 5.66 0.;
node 148 33. 5.66 0.;
node 149 34.5 5.66 0.;
node 150 36. 5.66 0.;
node 151 0. 6.86 0.;
node 152 1.5 6.86 0.;
node 153 3. 6.86 0.;
node 154 4.5 6.86 0.;
node 155 6. 6.86 0.;
node 156 7.5 6.86 0.;
node 157 9. 6.86 0.;
node 158 10.5 6.86 0.;
node 159 12. 6.86 0.;
node 160 13.5 6.86 0.;
node 161 15. 6.86 0.;
node 162 16.5 6.86 0.;
node 163 18. 6.86 0.;
node 164 19.5 6.86 0.;
node 165 21. 6.86 0.;
node 166 22.5 6.86 0.;
node 167 24. 6.86 0.;
node 168 25.5 6.86 0.;
node 169 27. 6.86 0.;
node 170 28.5 6.86 0.;
node 171 30. 6.86 0.;
node 172 31.5 6.86 0.;
node 173 33. 6.86 0.;
node 174 34.5 6.86 0.;
node 175 36. 6.86 0.;
node 176 0. 7.86 0.;
node 177 1.5 7.86 0.;
node 178 3. 7.86 0.;
node 179 4.5 7.86 0.;
node 180 6. 7.86 0.;
node 181 7.5 7.86 0.;
node 182 9. 7.86 0.;
node 183 10.5 7.86 0.;
node 184 12. 7.86 0.;
node 185 13.5 7.86 0.;
node 186 15. 7.86 0.;
node 187 16.5 7.86 0.;
node 188 18. 7.86 0.;
node 189 19.5 7.86 0.;
node 190 21. 7.86 0.;
node 191 22.5 7.86 0.;
node 192 24. 7.86 0.;
node 193 25.5 7.86 0.;
node 194 27. 7.86 0.;
node 195 28.5 7.86 0.;
node 196 30. 7.86 0.;
node 197 31.5 7.86 0.;
node 198 33. 7.86 0.;
node 199 34.5 7.86 0.;
node 200 36. 7.86 0.;
node 201 0. 9.12 0.;
node 202 1.5 9.12 0.;
node 203 3. 9.12 0.;
node 204 4.5 9.12 0.;
node 205 6. 9.12 0.;
node 206 7.5 9.12 0.;
node 207 9. 9.12 0.;
node 208 10.5 9.12 0.;
node 209 12. 9.12 0.;
node 210 13.5 9.12 0.;
node 211 15. 9.12 0.
node 212 16.5 9.12 0.;
node 213 18. 9.12 0.;
node 214 19.5 9.12 0.;
node 215 21. 9.12 0.;
node 216 22.5 9.12 0.;
node 217 24. 9.12 0.;
node 218 25.5 9.12 0.;
node 219 27. 9.12 0.;
node 220 28.5 9.12 0.;
node 221 30. 9.12 0.;
node 222 31.5 9.12 0.;
node 223 33. 9.12 0.;
node 224 34.5 9.12 0.;
node 225 36. 9.12 0.;
#Columns------------
node 226 18. 2.26 -1.;
node 227 18. 6.86 -1.;
node 228 18. 2.26 -2.;
node 229 18. 6.86 -2.;
node 230 18. 2.26 -3.;
node 231 18. 6.86 -3. ;
node 232 18. 2.26 -4.;
node 233 18. 6.86 -4.;
#DEFINE BOUNDARY CONDITIONS-------------------------------------------------------------------------------
fix 26 1 1 1 1 1 1; #Left hand supports
fix 76 1 1 1 1 1 1;
fix 126 1 1 1 1 1 1;
fix 176 1 1 1 1 1 1;
fix 232 1 1 1 1 1 1; #Central supports
fix 233 1 1 1 1 1 1;
fix 50 1 1 1 1 1 1; #Right hand supports
fix 100 1 1 1 1 1 1;
fix 150 1 1 1 1 1 1;
fix 200 1 1 1 1 1 1;
#DEFINE MATERIAL PARAMETERS----------------------------------------------------------
set Ec 30000000000.; #Elastic modulus in N/m^2
set fc -40000000.; #Compressive strength in N/m^2
set v .2;
set Gc [expr $Ec/(2.*(1+$v))];
set dens 2400.; #density in Kg/m^3
set IDconcCore 1; # material ID tag -- confined core concrete
set IDconcCover 2; # material ID tag -- unconfined cover concrete
set IDreinf 3; # material ID tag -- reinforcement
#confined concrete
set Kfc 1.3; # ratio of confined to unconfined concrete strength
set fc1C [expr $Kfc*$fc]; # CONFINED concrete (mander model), maximum stress
set eps1C [expr 2.*$fc1C/$Ec]; # strain at maximum stress
set fc2C [expr 0.2*$fc1C]; # ultimate stress
set eps2C [expr 5*$eps1C]; # strain at ultimate stress
#unconfined concrete
set fc1U $fc; # UNCONFINED concrete (todeschini parabolic model), maximum stress
set eps1U -0.003; # strain at maximum strength of unconfined concrete
set fc2U [expr 0.2*$fc1U]; # ultimate stress
set eps2U -0.01; # strain at ultimate stress
set lambda 0.1; # ratio between unloading slope at $eps2 and initial slope $Ec
#tensile-strength properties
set ftC [expr -0.14*$fc1C]; # tensile strength +tension
set ftU [expr -0.14*$fc1U]; # tensile strength +tension
set Ets [expr $ftU/0.002]; # tension softening stiffness
#-----------
set Fy 460000000; # STEEL yield stress
set Es 200000000000; # modulus of steel
set Bs 0.01; # strain-hardening ratio
set R0 18; # control the transition from elastic to plastic branches
set cR1 0.925; # control the transition from elastic to plastic branches
set cR2 0.15; # control the transition from elastic to plastic branches
uniaxialMaterial Concrete02 $IDconcCore $fc1C $eps1C $fc2C $eps2C $lambda $ftC $Ets; # build core concrete (confined)
uniaxialMaterial Concrete02 $IDconcCover $fc1U $eps1U $fc2U $eps2U $lambda $ftU $Ets; # build cover concrete (unconfined)
uniaxialMaterial Steel02 $IDreinf $Fy $Es $Bs $R0 $cR1 $cR2; # build reinforcement material
##DEFINE SECTIONS-------------------------------------------------------------
#Slab Sections--------------------
set LongSec1Tag 1;
set LongSec2Tag 2;
set LongSec3Tag 3;
set LongSec4Tag 4;
set TransSec1Tag 5;
set TransSec2Tag 6;
set ALongSec1 [expr .87*.8];
set IzLongSec1 [expr 1./12.*.87*pow(.8,3)];
set IyLongSec1 [expr 1./12.*.8*pow(.87,3)];
set JLongSec1 [expr 1./6.*.87*pow(.8,3)];
set ALongSec2 [expr 1.13*.8];
set IzLongSec2 [expr 1./12.*1.13*pow(.8,3)];
set IyLongSec2 [expr 1./12.*.8*pow(1.13,3)];
set JLongSec2 [expr 1./6.*1.13*pow(.8,3)];
set ALongSec3 [expr 1.1*.8];
set IzLongSec3 [expr 1./12.*1.1*pow(.8,3)];
set IyLongSec3 [expr 1./12.*.8*pow(1.1,3)];
set JLongSec3 [expr 1./6.*1.1*pow(.8,3)];
set ALongSec4 [expr 1.15*.8];
set IzLongSec4 [expr 1./12.*1.15*pow(.8,3)];
set IyLongSec4 [expr 1./12.*.8*pow(1.15,3)];
set JLongSec4 [expr 1./6.*1.15*pow(.8,3)];
set ATransSec1 [expr 1.25*.8];
set IzTransSec1 [expr 1./12.*1.25*pow(.8,3)];
set IyTransSec1 [expr 1./12.*.8*pow(1.25,3)];
set JTransSec1 [expr 1./6.*1.25*pow(.8,3)];
set ATransSec2 [expr 1.5*.8];
set IzTransSec2 [expr 1./12.*1.5*pow(.8,3)];
set IyTransSec2 [expr 1./12.*.8*pow(1.5,3)];
set JTransSec2 [expr 1./6.*1.5*pow(.8,3)];
section Elastic $LongSec1Tag $Ec $ALongSec1 $IzLongSec1 $IyLongSec1 $Gc $JLongSec1;
section Elastic $LongSec2Tag $Ec $ALongSec2 $IzLongSec2 $IyLongSec2 $Gc $JLongSec2;
section Elastic $LongSec3Tag $Ec $ALongSec3 $IzLongSec3 $IyLongSec3 $Gc $JLongSec3;
section Elastic $LongSec4Tag $Ec $ALongSec4 $IzLongSec4 $IyLongSec4 $Gc $JLongSec4;
section Elastic $TransSec1Tag $Ec $ATransSec1 $IzTransSec1 $IyTransSec1 $Gc $JTransSec1;
section Elastic $TransSec2Tag $Ec $ATransSec2 $IzTransSec2 $IyTransSec2 $Gc $JTransSec2;
#Column Sections--------------------------------
set DSec 1.; # Column Diameter
set coverSec .1; # Column cover to reinforcing steel NA.
set numBarsSec 16; # number of uniformly-distributed longitudinal-reinforcement bars
set barAreaSec .001500; # area of longitudinal-reinforcement bars
set ColSecTag 7; # set tag for symmetric section
# Notes
# The center of the reinforcing bars are placed at the inner radius
# The core concrete ends at the inner radius (same as reinforcing bars)
# The reinforcing bars are all the same size
# The center of the section is at (0,0) in the local axis system
# Zero degrees is along section y-axis
set ri 0.0; # inner radius of the section, only for hollow sections
set ro [expr $DSec/2]; # overall (outer) radius of the section
set nfCoreR 8; # number of radial divisions in the core (number of "rings")
set nfCoreT 8; # number of theta divisions in the core (number of "wedges")
set nfCoverR 4; # number of radial divisions in the cover
set nfCoverT 8; # number of theta divisions in the cover
#Define the fiber section
section fiberSec $ColSecTag {
set rc [expr $ro-$coverSec]; # Core radius
patch circ $IDconcCore $nfCoreT $nfCoreR 0 0 $ri $rc 0 360; # Define the core patch
patch circ $IDconcCover $nfCoverT $nfCoverR 0 0 $rc $ro 0 360; # Define the cover patch
set theta [expr 360.0/$numBarsSec]; # Determine angle increment between bars
layer circ $IDreinf $numBarsSec $barAreaSec 0 0 $rc $theta 360; # Define the reinforcing layer
}
#assign torsional Stiffness for 3D Model
set SecTag2 99; # ID tag for torsional section behavior
set SecTorTag3D 8; # ID tag for combined behavior for 3D model
uniaxialMaterial Elastic $SecTorTag3D 1.0e10; # define elastic torsional stiffness
section Aggregator $SecTag2 $SecTorTag3D T -section $ColSecTag; # combine section properties
#define geometric transformation: performs a linear geometric transformation of beam stiffness and resisting force from the basic system to the global-coordinate system
set LongBeamTransfTag 1;
set TransBeamTransfTag 2;
set ColTransfTag 3;
set ColTransfType PDelta; # options, Linear PDelta Corotational
geomTransf $ColTransfType $ColTransfTag 0 -1 0 ; # only columns can have PDelta effects (gravity effects)
geomTransf Linear $LongBeamTransfTag 0 -1 0;
geomTransf Linear $TransBeamTransfTag 1 0 0;
#DEFINE NODAL MASSES---------------------------------------------------------
set Load1 [expr 1.25*.87*.8*$dens];
set Load2 [expr 1.5*.87*.8*$dens];
set Load3 [expr 1.25*1.13*.8*$dens];
set Load4 [expr 1.5*1.13*.8*$dens];
set Load5 [expr 1.25*1.1*.8*$dens];
set Load6 [expr 1.5*1.1*.8*$dens];
set Load7 [expr 1.25*1.15*.8*$dens];
set Load8 [expr 1.5*1.15*.8*$dens];
set LoadCol [expr 3.14*(pow($ro,2))*4.*$dens];
mass 1 1.0e-9 $Load1 1.0e-9 0. 0. 0.;
mass 2 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 3 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 4 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 5 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 6 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 7 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 8 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 9 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 10 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 11 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 12 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 13 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 14 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 15 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 16 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 17 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 18 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 19 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 20 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 21 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 22 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 23 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 24 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 25 1.0e-9 $Load1 1.0e-9 0. 0. 0.;
mass 26 1.0e-9 $Load3 1.0e-9 0. 0. 0.;
mass 27 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 28 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 29 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 30 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 31 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 32 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 33 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 34 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 35 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 36 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 37 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 38 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 39 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 40 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 41 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 42 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 43 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 44 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 45 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 46 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 47 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 48 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 49 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 50 1.0e-9 $Load3 1.0e-9 0. 0. 0.;
mass 51 1.0e-9 $Load5 1.0e-9 0. 0. 0.;
mass 52 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 53 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 54 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 55 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 56 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 57 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 58 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 59 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 60 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 61 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 62 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 63 1.0e-9 [expr $Load6+($LoadCol/8)] 1.0e-9 0. 0. 0.;
mass 64 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 65 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 66 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 67 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 68 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 69 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 70 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 71 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 72 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 73 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 74 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 75 1.0e-9 $Load5 1.0e-9 0. 0. 0.;
mass 76 1.0e-9 $Load7 1.0e-9 0. 0. 0.;
mass 77 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 78 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 79 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 80 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 81 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 82 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 83 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 84 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 85 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 86 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 87 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 88 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 89 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 90 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 91 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 92 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 93 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 94 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 95 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 96 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 97 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 98 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 99 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 100 1.0e-9 $Load7 1.0e-9 0. 0. 0.;
mass 101 1.0e-9 $Load5 1.0e-9 0. 0. 0.;
mass 102 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 103 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 104 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 105 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 106 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 107 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 108 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 109 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 110 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 111 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 112 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 113 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 114 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 115 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 116 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 117 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 118 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 119 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 120 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 121 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 122 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 123 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 124 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 125 1.0e-9 $Load5 1.0e-9 0. 0. 0.;
mass 126 1.0e-9 $Load7 1.0e-9 0. 0. 0.;
mass 127 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 128 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 129 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 130 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 131 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 132 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 133 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 134 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 135 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 136 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 137 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 138 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 139 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 140 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 141 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 142 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 143 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 144 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 145 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 146 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 147 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 148 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 149 1.0e-9 $Load8 1.0e-9 0. 0. 0.;
mass 150 1.0e-9 $Load7 1.0e-9 0. 0. 0.;
mass 151 1.0e-9 $Load5 1.0e-9 0. 0. 0.;
mass 152 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 153 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 154 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 155 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 156 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 157 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 158 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 159 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 160 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 161 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 162 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 163 1.0e-9 [expr $Load6+($LoadCol/8)] 1.0e-9 0. 0. 0.;
mass 164 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 165 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 166 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 167 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 168 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 169 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 170 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 171 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 172 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 173 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 174 1.0e-9 $Load6 1.0e-9 0. 0. 0.;
mass 175 1.0e-9 $Load5 1.0e-9 0. 0. 0.;
mass 176 1.0e-9 $Load3 1.0e-9 0. 0. 0.;
mass 177 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 178 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 179 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 180 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 181 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 182 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 183 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 184 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 185 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 186 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 187 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 188 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 189 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 190 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 191 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 192 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 193 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 194 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 195 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 196 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 197 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 198 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 199 1.0e-9 $Load4 1.0e-9 0. 0. 0.;
mass 200 1.0e-9 $Load3 1.0e-9 0. 0. 0.;
mass 201 1.0e-9 $Load1 1.0e-9 0. 0. 0.;
mass 202 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 203 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 204 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 205 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 206 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 207 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 208 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 209 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 210 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 211 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 212 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 213 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 214 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 215 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 216 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 217 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 218 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 219 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 220 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 221 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 222 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 223 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 224 1.0e-9 $Load2 1.0e-9 0. 0. 0.;
mass 225 1.0e-9 $Load1 1.0e-9 0. 0. 0.;
mass 226 1.0e-9 [expr $LoadCol/4] 1.0e-9 0. 0. 0.;
mass 227 1.0e-9 [expr $LoadCol/4] 1.0e-9 0. 0. 0.;
mass 228 1.0e-9 [expr $LoadCol/4] 1.0e-9 0. 0. 0.;
mass 229 1.0e-9 [expr $LoadCol/4] 1.0e-9 0. 0. 0.;
mass 230 1.0e-9 [expr $LoadCol/4] 1.0e-9 0. 0. 0.;
mass 231 1.0e-9 [expr $LoadCol/4] 1.0e-9 0. 0. 0.;
mass 232 1.0e-9 [expr $LoadCol/8] 1.0e-9 0. 0. 0.;
mass 233 1.0e-9 [expr $LoadCol/8] 1.0e-9 0. 0. 0.;
#element connectivity:
#longitudinal elements
set numIntgrPts 5;
element dispBeamColumn 1 1 2 $numIntgrPts $LongSec1Tag $LongBeamTransfTag; # self-explanatory when using variables
element dispBeamColumn 2 2 3 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 3 3 4 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 4 4 5 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 5 5 6 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 6 6 7 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 7 7 8 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 8 8 9 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 9 9 10 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 10 10 11 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 11 11 12 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 12 12 13 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 13 13 14 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 14 14 15 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 15 15 16 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 16 16 17 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 17 17 18 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 18 18 19 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 19 19 20 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 20 20 21 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 21 21 22 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 22 22 23 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 23 23 24 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 24 24 25 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 25 26 27 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 26 27 28 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 27 28 29 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 28 29 30 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 29 30 31 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 30 31 32 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 31 32 33 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 32 33 34 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 33 34 35 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 34 35 36 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 35 36 37 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 36 37 38 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 37 38 39 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 38 39 40 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 39 40 41 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 40 41 42 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 41 42 43 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 42 43 44 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 43 44 45 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 44 45 46 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 45 46 47 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 46 47 48 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 47 48 49 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 48 49 50 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 49 51 52 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 50 52 53 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 51 53 54 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 52 54 55 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 53 55 56 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 54 56 57 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 55 57 58 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 56 58 59 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 57 59 60 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 58 60 61 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 59 61 62 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 60 62 63 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 61 63 64 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 62 64 65 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 63 65 66 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 64 66 67 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 65 67 68 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 66 68 69 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 67 69 70 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 68 70 71 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 69 71 72 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 70 72 73 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 71 73 74 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 72 74 75 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 73 76 77 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 74 77 78 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 75 78 79 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 76 79 80 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 77 80 81 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 78 81 82 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 79 82 83 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 80 83 84 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 81 84 85 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 82 85 86 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 83 86 87 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 84 87 88 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 85 88 89 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 86 89 90 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 87 90 91 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 88 91 92 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 89 92 93 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 90 93 94 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 91 94 95 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 92 95 96 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 93 96 97 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 94 97 98 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 95 98 99 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 96 99 100 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 97 101 102 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 98 102 103 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 99 103 104 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 100 104 105 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 101 105 106 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 102 106 107 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 103 107 108 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 104 108 109 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 105 109 110 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 106 110 111 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 107 111 112 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 108 112 113 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 109 113 114 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 110 114 115 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 111 115 116 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 112 116 117 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 113 117 118 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 114 118 119 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 115 119 120 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 116 120 121 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 117 121 122 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 118 122 123 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 119 123 124 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 120 124 125 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 121 126 127 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 122 127 128 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 123 128 129 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 124 129 130 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 125 130 131 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 126 131 132 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 127 132 133 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 128 133 134 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 129 134 135 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 130 135 136 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 131 136 137 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 132 137 138 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 133 138 139 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 134 139 140 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 135 140 141 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 136 141 142 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 137 142 143 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 138 143 144 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 139 144 145 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 140 145 146 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 141 146 147 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 142 147 148 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 143 148 149 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 144 149 150 $numIntgrPts $LongSec4Tag $LongBeamTransfTag;
element dispBeamColumn 145 151 152 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 146 152 153 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 147 153 154 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 148 154 155 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 149 155 156 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 150 156 157 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 151 157 158 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 152 158 159 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 153 159 160 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 154 160 161 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 155 161 162 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 156 162 163 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 157 163 164 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 158 164 165 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 159 165 166 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 160 166 167 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 161 167 168 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 162 168 169 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 163 169 170 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 164 170 171 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 165 171 172 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 166 172 173 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 167 173 174 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 168 174 175 $numIntgrPts $LongSec3Tag $LongBeamTransfTag;
element dispBeamColumn 169 176 177 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 170 177 178 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 171 178 179 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 172 179 180 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 173 180 181 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 174 181 182 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 175 182 183 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 176 183 184 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 177 184 185 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 178 185 186 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 179 186 187 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 180 187 188 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 181 188 189 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 182 189 190 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 183 190 191 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 184 191 192 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 185 192 193 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 186 193 194 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 187 194 195 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 188 195 196 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 189 196 197 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 190 197 198 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 191 198 199 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 192 199 200 $numIntgrPts $LongSec2Tag $LongBeamTransfTag;
element dispBeamColumn 193 201 202 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 194 202 203 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 195 203 204 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 196 204 205 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 197 205 206 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 198 206 207 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 199 207 208 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 200 208 209 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 201 209 210 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 202 210 211 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 203 211 212 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 204 212 213 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 205 213 214 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 206 214 215 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 207 215 216 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 208 216 217 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 209 217 218 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 210 218 219 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 211 219 220 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 212 220 221 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 213 221 222 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 214 222 223 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 215 223 224 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
element dispBeamColumn 216 224 225 $numIntgrPts $LongSec1Tag $LongBeamTransfTag;
#transverse element
element dispBeamColumn 217 1 26 $numIntgrPts $TransSec1Tag $TransBeamTransfTag;
element dispBeamColumn 218 26 51 $numIntgrPts $TransSec1Tag $TransBeamTransfTag;
element dispBeamColumn 219 51 76 $numIntgrPts $TransSec1Tag $TransBeamTransfTag;
element dispBeamColumn 220 76 101 $numIntgrPts $TransSec1Tag $TransBeamTransfTag;
element dispBeamColumn 221 101 126 $numIntgrPts $TransSec1Tag $TransBeamTransfTag;
element dispBeamColumn 222 126 151 $numIntgrPts $TransSec1Tag $TransBeamTransfTag;
element dispBeamColumn 223 151 176 $numIntgrPts $TransSec1Tag $TransBeamTransfTag;
element dispBeamColumn 224 176 201 $numIntgrPts $TransSec1Tag $TransBeamTransfTag;
element dispBeamColumn 225 2 27 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 226 27 52 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 227 52 77 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 228 77 102 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 229 102 127 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 230 127 152 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 231 152 177 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 232 177 202 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 233 3 28 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 234 28 53 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 235 53 78 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 236 78 103 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 237 103 128 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 238 128 153 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 239 153 178 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 240 178 203 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 241 4 29 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 242 29 54 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 243 54 79 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 244 79 104 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 245 104 129 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 246 129 154 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 247 154 179 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 248 179 204 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 249 5 30 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 250 30 55 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 251 55 80 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 252 80 105 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 253 105 130 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 254 130 155 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 255 155 180 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 256 180 205 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 257 6 31 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 258 31 56 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 259 56 81 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 260 81 106 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 261 106 131 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 262 131 156 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 263 156 181 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 264 181 206 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 265 7 32 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 266 32 57 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 267 57 82 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 268 82 107 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 269 107 132 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 270 132 157 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 271 157 182 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 272 182 207 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 273 8 33 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 274 33 58 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 275 58 83 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 276 83 108 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 277 108 133 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 278 133 158 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 279 158 183 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 280 183 208 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 281 9 34 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 282 34 59 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 283 59 84 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 284 84 109 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 285 109 134 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 286 134 159 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 287 159 184 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 288 184 209 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 289 10 35 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 290 35 60 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 291 60 85 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 292 85 110 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 293 110 135 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 294 135 160 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 295 160 185 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 296 185 210 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 297 11 36 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 298 36 61 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 299 61 86 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 300 86 111 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 301 111 136 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 302 136 161 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 303 161 186 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 304 186 211 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 305 12 37 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 306 37 62 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 307 62 87 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 308 87 112 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 309 112 137 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 310 137 162 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 311 162 187 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 312 187 212 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 313 13 38 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 314 38 63 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 315 63 88 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 316 88 113 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 317 113 138 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 318 138 163 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 319 163 188 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 320 188 213 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 321 14 39 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 322 39 64 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 323 64 89 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 324 89 114 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 325 114 139 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 326 139 164 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 327 164 189 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 328 189 214 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 329 15 40 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 330 40 65 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 331 65 90 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 332 90 115 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 333 115 140 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 334 140 165 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 335 165 190 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 336 190 215 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 337 16 41 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 338 41 66 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 339 66 91 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 340 91 116 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 341 116 141 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 342 141 166 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 343 166 191 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 344 191 216 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 345 17 42 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 346 42 67 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 347 67 92 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 348 92 117 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 349 117 142 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 350 142 167 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 351 167 192 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 352 192 217 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 353 18 43 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 354 43 68 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 355 68 93 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 356 93 118 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 357 118 143 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 358 143 168 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 359 168 193 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 360 193 218 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 361 19 44 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 362 44 69 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 363 69 94 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 364 94 119 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 365 119 144 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 366 144 169 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 367 169 194 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 368 194 219 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 369 20 45 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 370 45 70 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 371 70 95 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 372 95 120 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 373 120 145 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 374 145 170 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 375 170 195 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 376 195 220 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 377 21 46 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 378 46 71 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 379 71 96 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 380 96 121 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 381 121 146 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 382 146 171 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 383 171 196 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 384 196 221 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 385 22 47 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 386 47 72 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 387 72 97 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 388 97 122 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 389 122 147 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 390 147 172 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 391 172 197 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 392 197 222 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 393 23 48 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 394 48 73 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 395 73 98 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 396 98 123 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 397 123 148 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 398 148 173 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 399 173 198 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 400 198 223 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 401 24 49 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 402 49 74 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 403 74 99 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 404 99 124 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 405 124 149 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 406 149 174 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 407 174 199 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 408 199 224 $numIntgrPts $TransSec2Tag $TransBeamTransfTag;
element dispBeamColumn 409 25 50 $numIntgrPts $TransSec1Tag $TransBeamTransfTag;
element dispBeamColumn 410 50 75 $numIntgrPts $TransSec1Tag $TransBeamTransfTag;
element dispBeamColumn 411 75 100 $numIntgrPts $TransSec1Tag $TransBeamTransfTag;
element dispBeamColumn 412 100 125 $numIntgrPts $TransSec1Tag $TransBeamTransfTag;
element dispBeamColumn 413 125 150 $numIntgrPts $TransSec1Tag $TransBeamTransfTag;
element dispBeamColumn 414 150 175 $numIntgrPts $TransSec1Tag $TransBeamTransfTag;
element dispBeamColumn 415 175 200 $numIntgrPts $TransSec1Tag $TransBeamTransfTag;
element dispBeamColumn 416 200 225 $numIntgrPts $TransSec1Tag $TransBeamTransfTag;
#Column Elements-------------------------------------------------------------------------
element dispBeamColumn 417 63 226 $numIntgrPts $SecTag2 $ColTransfTag;
element dispBeamColumn 418 163 227 $numIntgrPts $SecTag2 $ColTransfTag;
element dispBeamColumn 419 226 228 $numIntgrPts $SecTag2 $ColTransfTag;
element dispBeamColumn 420 227 229 $numIntgrPts $SecTag2 $ColTransfTag;
element dispBeamColumn 421 228 230 $numIntgrPts $SecTag2 $ColTransfTag;
element dispBeamColumn 422 229 231 $numIntgrPts $SecTag2 $ColTransfTag;
element dispBeamColumn 423 230 232 $numIntgrPts $SecTag2 $ColTransfTag;
element dispBeamColumn 424 231 233 $numIntgrPts $SecTag2 $ColTransfTag;
recorder Element -file $dataDir/Fcentre.out -time -ele 61 62 63 64 65 66 67 68 69 70 71 72 globalForce;
recorder Node -file $dataDir/Rsupport.out -time -node 26 76 126 176 232 233 50 100 150 200 -dof 3 reaction;
recorder Node -file $dataDir/Disp.out -time -node 63 163 -dof 1 2 3 disp;
#DEFINE GRAVITY---------------------------------------------------------------------
set g 9.81; #accel due to gravity
#set gravity loads
set UDL1 [expr .87*.8*$dens*$g];
set UDL2 [expr 1.13*.8*$dens*$g];
set UDL3 [expr 1.1*.8*$dens*$g];
set UDL4 [expr 1.15*.8*$dens*$g];
set UDLCol [expr 3.14*pow($ro,2)*$dens*$g];
#Gravity Load Case
pattern Plain 1 Linear {
eleLoad -ele 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 -type -beamUniform -$UDL1 0. 0.; #UDL on longitudinal beam 1
eleLoad -ele 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 -type -beamUniform -$UDL2 0. 0.;
eleLoad -ele 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 -type -beamUniform -$UDL3 0. 0.;
eleLoad -ele 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 -type -beamUniform -$UDL4 0. 0.;
eleLoad -ele 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 -type -beamUniform -$UDL3 0. 0.;
eleLoad -ele 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 -type -beamUniform -$UDL4 0. 0.;
eleLoad -ele 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 -type -beamUniform -$UDL3 0. 0.;
eleLoad -ele 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 -type -beamUniform -$UDL2 0. 0.;
eleLoad -ele 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 -type -beamUniform -$UDL1 0. 0.;
eleLoad -ele 417 418 419 420 421 422 423 424 -type -beamUniform 0. 0. $UDLCol;
}
#Start of Analysis------------------------------------------------------------------
#Gravity analysis
#Set up convergence criteria:
set Tol 1.0e-8; # convergence tolerance for test
set NumIter 10; # Set the number of iterations for convergence
test NormDispIncr $Tol $NumIter;
constraints Plain;
numberer RCM;
system BandGeneral;
algorithm Newton;
set LoadSteps 10;
set LoadIncr [expr 1.0/$LoadSteps];
integrator LoadControl $LoadIncr;
analysis Static;
analyze $LoadSteps; # applying gravity
#maintain constant gravity loads and reset time to zero
loadConst -time 0.0
puts "Model built"
Grillage analysis of Bridge
Moderators: silvia, selimgunay, Moderators
Re: Grillage analysis of Bridge
if the solution is correct, nodal equilibrium should be in effect.
Re: Grillage analysis of Bridge
Thank you very much. Do you have any ideas about what is likely to be wrong? I have checked and rechecked the model alot so any advice on what to look out for?
Re: Grillage analysis of Bridge
i don't have time to look at large models .. do a small model of less than 10 elements and tell me what node and elements you have problems with.
Re: Grillage analysis of Bridge
I think the problem is that I am reading the results incorrectly. I have included a small model below with 17 elements. It is a 4*1.8*.25 slab with 3 longitudinal and 4 transverse members. Sorry to bother you but could you run it and tell me what values the columns in the globalForce recorder correspond to? I understood from the manual that it was Fx Fy Fz Mx My Mz for the i and j nodes of each element but I think I may be wrong as nodal equilibrium is not in effect if that is the case.
I have recorded the forces in elements 4 5 6 (the central longitudinal member)
Thanks very much for your time.
#SET UP ----------------------------------------------------------------------------
wipe; # clear memory of all past model definitions
model BasicBuilder -ndm 3 -ndf 6; # Define the model builder, ndm=#dimension, ndf=#dofs
set dataDir DataSlab1; # set up name of data directory (you can remove this)
file mkdir $dataDir; # create data directory
#define geometry
node 1 0. 0. 0.;
node 2 1. 0. 0.;
node 3 2. 0. 0.;
node 4 3. 0. 0.;
node 5 0. .6 0.;
node 6 1. .6 0.;
node 7 2. .6 0.;
node 8 3. .6 0.;
node 9 0. 1.2 0.;
node 10 1. 1.2 0.;
node 11 2. 1.2 0.;
node 12 3. 1.2 0.;
#left hand support
fix 5 1 1 1 1 1 1;
#right hand support
fix 8 1 1 1 1 1 1;
#DEFINE MATERIAL PARAMETERS----------------------------------------------------------
set Ec 30000000000.; #Elastic modulus in N/m^2
set fc -40000000.; #Compressive strength in N/m^2
set v .2;
set Gc [expr $Ec/(2.*(1+$v))];
set dens 2400.; #density in Kg/m^3
#DEFINE SECTIONS------------------------------------------------------------------
set LongSecTag 1;
set TransSecTag 2;
#Longitudinal Section
set ALongSec [expr .6*.25];
set IzLongSec [expr 1./12.*.6*pow(.25,3)];
set IyLongSec [expr 1./12.*.25*pow(.6,3)];
set JLongSec [expr 1./6.*.6*pow(.25,3)];
#Transverse Section
set ATransSec [expr 1.*.25];
set IzTransSec [expr 1./12.*1.*pow(.25,3)];
set IyTransSec [expr 1./12.*.25*pow(1.,3)];
set JTransSec [expr 1./6.*1.*pow(.25,3)];
section Elastic $LongSecTag $Ec $ALongSec $IzLongSec $IyLongSec $Gc $JLongSec;
section Elastic $TransSecTag $Ec $ATransSec $IzTransSec $IyTransSec $Gc $JTransSec;
# DEFINE GEOMETRIC TRANSFORMATION-------------------------------------------------
set LongBeamTransfTag 1;
set TransBeamTransfTag 2;
geomTransf Linear $LongBeamTransfTag 0 -1 0;
geomTransf Linear $TransBeamTransfTag 1 0 0;
# DEFINE NODAL MASSES---------------------------------------------------------------
set Load [expr .6*1.*.25*$dens];
mass 1 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 2 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 3 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 4 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 5 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 6 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 7 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 8 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 9 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 10 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 11 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 12 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
#define elements
set numIntgrPts 5;
#longitudinal elements
element nonlinearBeamColumn 1 1 2 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 2 2 3 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 3 3 4 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 4 5 6 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 5 6 7 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 6 7 8 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 7 9 10 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 8 10 11 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 9 11 12 $numIntgrPts $LongSecTag $LongBeamTransfTag;
#transverse elements
element nonlinearBeamColumn 10 1 5 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 11 5 9 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 12 2 6 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 13 6 10 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 14 3 7 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 15 7 11 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 16 4 8 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 17 8 12 $numIntgrPts $TransSecTag $TransBeamTransfTag;
recorder Element -file $dataDir/Force.out -time -ele 4 5 6 globalForce;
recorder Node -file $dataDir/RSupport.out -time -node 5 8 -dof 3 reaction;
# DEFINE GRAVITY---------------------------------------------------
set g 9.81; #acceleration due to gravity
set UDL [expr $ALongSec*$dens*$g]; #longitudinal member self weight
# Gravity loadcase
pattern Plain 1 Linear {
eleLoad -ele 1 2 3 4 5 6 7 8 9 -type -beamUniform -$UDL 0. 0.; #Load the longitudinal members with self weight
}
# ANALYSIS--------------
# Set up convergence criteria:
set Tol 1.0e-8; # convergence tolerance for test
set NumIter 10; # Set the number of iterations for convergence
test NormDispIncr $Tol $NumIter;
constraints Plain;
numberer RCM;
system BandGeneral;
algorithm Newton;
set LoadSteps 10;
set LoadIncr [expr 1.0/$LoadSteps];
integrator LoadControl $LoadIncr;
analysis Static;
analyze $LoadSteps; # applying gravity
#maintain constant gravity loads and reset time to zero
loadConst -time 0.0
puts "Model built"
I have recorded the forces in elements 4 5 6 (the central longitudinal member)
Thanks very much for your time.
#SET UP ----------------------------------------------------------------------------
wipe; # clear memory of all past model definitions
model BasicBuilder -ndm 3 -ndf 6; # Define the model builder, ndm=#dimension, ndf=#dofs
set dataDir DataSlab1; # set up name of data directory (you can remove this)
file mkdir $dataDir; # create data directory
#define geometry
node 1 0. 0. 0.;
node 2 1. 0. 0.;
node 3 2. 0. 0.;
node 4 3. 0. 0.;
node 5 0. .6 0.;
node 6 1. .6 0.;
node 7 2. .6 0.;
node 8 3. .6 0.;
node 9 0. 1.2 0.;
node 10 1. 1.2 0.;
node 11 2. 1.2 0.;
node 12 3. 1.2 0.;
#left hand support
fix 5 1 1 1 1 1 1;
#right hand support
fix 8 1 1 1 1 1 1;
#DEFINE MATERIAL PARAMETERS----------------------------------------------------------
set Ec 30000000000.; #Elastic modulus in N/m^2
set fc -40000000.; #Compressive strength in N/m^2
set v .2;
set Gc [expr $Ec/(2.*(1+$v))];
set dens 2400.; #density in Kg/m^3
#DEFINE SECTIONS------------------------------------------------------------------
set LongSecTag 1;
set TransSecTag 2;
#Longitudinal Section
set ALongSec [expr .6*.25];
set IzLongSec [expr 1./12.*.6*pow(.25,3)];
set IyLongSec [expr 1./12.*.25*pow(.6,3)];
set JLongSec [expr 1./6.*.6*pow(.25,3)];
#Transverse Section
set ATransSec [expr 1.*.25];
set IzTransSec [expr 1./12.*1.*pow(.25,3)];
set IyTransSec [expr 1./12.*.25*pow(1.,3)];
set JTransSec [expr 1./6.*1.*pow(.25,3)];
section Elastic $LongSecTag $Ec $ALongSec $IzLongSec $IyLongSec $Gc $JLongSec;
section Elastic $TransSecTag $Ec $ATransSec $IzTransSec $IyTransSec $Gc $JTransSec;
# DEFINE GEOMETRIC TRANSFORMATION-------------------------------------------------
set LongBeamTransfTag 1;
set TransBeamTransfTag 2;
geomTransf Linear $LongBeamTransfTag 0 -1 0;
geomTransf Linear $TransBeamTransfTag 1 0 0;
# DEFINE NODAL MASSES---------------------------------------------------------------
set Load [expr .6*1.*.25*$dens];
mass 1 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 2 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 3 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 4 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 5 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 6 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 7 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 8 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 9 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 10 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 11 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
mass 12 1.0e-9 1.0e-9 $Load 1.0e-9 1.0e-9 1.0e-9;
#define elements
set numIntgrPts 5;
#longitudinal elements
element nonlinearBeamColumn 1 1 2 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 2 2 3 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 3 3 4 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 4 5 6 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 5 6 7 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 6 7 8 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 7 9 10 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 8 10 11 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 9 11 12 $numIntgrPts $LongSecTag $LongBeamTransfTag;
#transverse elements
element nonlinearBeamColumn 10 1 5 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 11 5 9 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 12 2 6 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 13 6 10 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 14 3 7 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 15 7 11 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 16 4 8 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 17 8 12 $numIntgrPts $TransSecTag $TransBeamTransfTag;
recorder Element -file $dataDir/Force.out -time -ele 4 5 6 globalForce;
recorder Node -file $dataDir/RSupport.out -time -node 5 8 -dof 3 reaction;
# DEFINE GRAVITY---------------------------------------------------
set g 9.81; #acceleration due to gravity
set UDL [expr $ALongSec*$dens*$g]; #longitudinal member self weight
# Gravity loadcase
pattern Plain 1 Linear {
eleLoad -ele 1 2 3 4 5 6 7 8 9 -type -beamUniform -$UDL 0. 0.; #Load the longitudinal members with self weight
}
# ANALYSIS--------------
# Set up convergence criteria:
set Tol 1.0e-8; # convergence tolerance for test
set NumIter 10; # Set the number of iterations for convergence
test NormDispIncr $Tol $NumIter;
constraints Plain;
numberer RCM;
system BandGeneral;
algorithm Newton;
set LoadSteps 10;
set LoadIncr [expr 1.0/$LoadSteps];
integrator LoadControl $LoadIncr;
analysis Static;
analyze $LoadSteps; # applying gravity
#maintain constant gravity loads and reset time to zero
loadConst -time 0.0
puts "Model built"
Re: Grillage analysis of Bridge
Hi again,
I am still having the problem of having different forces at the same nodes when I introduce transverse elements. I have pasted a more simple model with only 8 nodes. For this model the Mx values at node 2 are not matching. For the previous example however the My values didnt match either.
I am new to modelling bridges and wondering if anyone has any advice on modelling the superstructure. Would it be satisfactory to model the entire superstructure cross section as one section rather than breaking it up into a grillage?
Thanks very much for your time.
Regards
Mairead
#SET UP ----------------------------------------------------------------------------
wipe; # clear memory of all past model definitions
model BasicBuilder -ndm 3 -ndf 6; # Define the model builder, ndm=#dimension, ndf=#dofs
set dataDir DataSlab3; # set up name of data directory (you can remove this)
file mkdir $dataDir; # create data directory
# define geometry
node 1 0. 0. 0.;
node 2 1. 0. 0.;
node 3 2. 0. 0.;
node 4 3. 0. 0.;
node 5 0. .6 0.;
node 6 1. .6 0.;
node 7 2. .6 0.;
node 8 3. .6 0.;
#boundary conditions
fix 1 1 1 1 1 1 1;
fix 5 1 1 1 1 1 1;
fix 4 1 1 1 1 1 1;
fix 8 1 1 1 1 1 1;
#DEFINE MATERIAL PARAMETERS----------------------------------------------------------
set Ec 30000000000.; #Elastic modulus in N/m^2
set fc -40000000.; #Compressive strength in N/m^2
set v .2;
set Gc [expr $Ec/(2.*(1+$v))];
set dens 2400.; #density in Kg/m^3
#DEFINE SECTIONS------------------------------------------------------------------
set LongSecTag 1;
set TransSecTag 2;
#Longitudinal Sections
set ALongSec [expr .6*.25];
set IzLongSec [expr 1./12.*.6*pow(.25,3)];
set IyLongSec [expr 1./12.*.25*pow(.6,3)];
set JLongSec [expr 1./6.*.6*pow(.25,3)];
set ATransSec [expr 1.*.25];
set IzTransSec [expr 1./12.*1.*pow(.25,3)];
set IyTransSec [expr 1./12.*.25*pow(1.,3)];
set JTransSec [expr 1./6.*1.*pow(.25,3)];
section Elastic $LongSecTag $Ec $ALongSec $IzLongSec $IyLongSec $Gc $JLongSec;
section Elastic $TransSecTag $Ec $ATransSec $IzTransSec $IyTransSec $Gc $JTransSec;
# DEFINE GEOMETRIC TRANSFORMATION-------------------------------------------------
set LongBeamTransfTag 1;
set TransBeamTransfTag 2;
geomTransf Linear $LongBeamTransfTag 0 -1 0;
geomTransf Linear $TransBeamTransfTag 1 0 0;
# DEFINE NODAL MASSES---------------------------------------------------------------
set Load [expr .6*1.*.25*$dens];
mass 1 $Load 1.0e-9 1.0e-9 1.0e-9 1.0e-9 1.0e-9;
mass 2 $Load 1.0e-9 1.0e-9 1.0e-9 1.0e-9 1.0e-9;
mass 3 $Load 1.0e-9 1.0e-9 1.0e-9 1.0e-9 1.0e-9;
mass 4 $Load 1.0e-9 1.0e-9 1.0e-9 1.0e-9 1.0e-9;
mass 5 $Load 1.0e-9 1.0e-9 1.0e-9 1.0e-9 1.0e-9;
mass 6 $Load 1.0e-9 1.0e-9 1.0e-9 1.0e-9 1.0e-9;
mass 7 $Load 1.0e-9 1.0e-9 1.0e-9 1.0e-9 1.0e-9;
mass 8 $Load 1.0e-9 1.0e-9 1.0e-9 1.0e-9 1.0e-9;
# define elements
set numIntgrPts 5;
# longitudinal
element nonlinearBeamColumn 1 1 2 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 2 2 3 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 3 3 4 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 4 5 6 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 5 6 7 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 6 7 8 $numIntgrPts $LongSecTag $LongBeamTransfTag;
# transverse
element nonlinearBeamColumn 7 1 5 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 8 2 6 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 9 3 7 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 10 4 8 $numIntgrPts $TransSecTag $TransBeamTransfTag;
recorder Element -file $dataDir/Force.out -time -ele 1 2 globalForce;
recorder Node -file $dataDir/RSupport.out -time -node 1 5 4 8 -dof 3 reaction;
# DEFINE GRAVITY---------------------------------------------------
set g 9.81;
set UDL [expr $ALongSec*$dens*$g];
# Gravity loadcase
pattern Plain 1 Linear {
eleLoad -ele 1 2 3 4 5 6 -type -beamUniform -$UDL 0. 0.;
}
#ANALYSIS--------------
# Set up convergence criteria:
set Tol 1.0e-8; # convergence tolerance for test
set NumIter 10; # Set the number of iterations for convergence
test NormDispIncr $Tol $NumIter;
constraints Plain;
numberer RCM;
system BandGeneral;
algorithm Newton;
set LoadSteps 10;
set LoadIncr [expr 1.0/$LoadSteps];
integrator LoadControl $LoadIncr;
analysis Static;
analyze $LoadSteps; # applying gravity
#maintain constant gravity loads and reset time to zero
loadConst -time 0.0
puts "Model built"
I am still having the problem of having different forces at the same nodes when I introduce transverse elements. I have pasted a more simple model with only 8 nodes. For this model the Mx values at node 2 are not matching. For the previous example however the My values didnt match either.
I am new to modelling bridges and wondering if anyone has any advice on modelling the superstructure. Would it be satisfactory to model the entire superstructure cross section as one section rather than breaking it up into a grillage?
Thanks very much for your time.
Regards
Mairead
#SET UP ----------------------------------------------------------------------------
wipe; # clear memory of all past model definitions
model BasicBuilder -ndm 3 -ndf 6; # Define the model builder, ndm=#dimension, ndf=#dofs
set dataDir DataSlab3; # set up name of data directory (you can remove this)
file mkdir $dataDir; # create data directory
# define geometry
node 1 0. 0. 0.;
node 2 1. 0. 0.;
node 3 2. 0. 0.;
node 4 3. 0. 0.;
node 5 0. .6 0.;
node 6 1. .6 0.;
node 7 2. .6 0.;
node 8 3. .6 0.;
#boundary conditions
fix 1 1 1 1 1 1 1;
fix 5 1 1 1 1 1 1;
fix 4 1 1 1 1 1 1;
fix 8 1 1 1 1 1 1;
#DEFINE MATERIAL PARAMETERS----------------------------------------------------------
set Ec 30000000000.; #Elastic modulus in N/m^2
set fc -40000000.; #Compressive strength in N/m^2
set v .2;
set Gc [expr $Ec/(2.*(1+$v))];
set dens 2400.; #density in Kg/m^3
#DEFINE SECTIONS------------------------------------------------------------------
set LongSecTag 1;
set TransSecTag 2;
#Longitudinal Sections
set ALongSec [expr .6*.25];
set IzLongSec [expr 1./12.*.6*pow(.25,3)];
set IyLongSec [expr 1./12.*.25*pow(.6,3)];
set JLongSec [expr 1./6.*.6*pow(.25,3)];
set ATransSec [expr 1.*.25];
set IzTransSec [expr 1./12.*1.*pow(.25,3)];
set IyTransSec [expr 1./12.*.25*pow(1.,3)];
set JTransSec [expr 1./6.*1.*pow(.25,3)];
section Elastic $LongSecTag $Ec $ALongSec $IzLongSec $IyLongSec $Gc $JLongSec;
section Elastic $TransSecTag $Ec $ATransSec $IzTransSec $IyTransSec $Gc $JTransSec;
# DEFINE GEOMETRIC TRANSFORMATION-------------------------------------------------
set LongBeamTransfTag 1;
set TransBeamTransfTag 2;
geomTransf Linear $LongBeamTransfTag 0 -1 0;
geomTransf Linear $TransBeamTransfTag 1 0 0;
# DEFINE NODAL MASSES---------------------------------------------------------------
set Load [expr .6*1.*.25*$dens];
mass 1 $Load 1.0e-9 1.0e-9 1.0e-9 1.0e-9 1.0e-9;
mass 2 $Load 1.0e-9 1.0e-9 1.0e-9 1.0e-9 1.0e-9;
mass 3 $Load 1.0e-9 1.0e-9 1.0e-9 1.0e-9 1.0e-9;
mass 4 $Load 1.0e-9 1.0e-9 1.0e-9 1.0e-9 1.0e-9;
mass 5 $Load 1.0e-9 1.0e-9 1.0e-9 1.0e-9 1.0e-9;
mass 6 $Load 1.0e-9 1.0e-9 1.0e-9 1.0e-9 1.0e-9;
mass 7 $Load 1.0e-9 1.0e-9 1.0e-9 1.0e-9 1.0e-9;
mass 8 $Load 1.0e-9 1.0e-9 1.0e-9 1.0e-9 1.0e-9;
# define elements
set numIntgrPts 5;
# longitudinal
element nonlinearBeamColumn 1 1 2 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 2 2 3 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 3 3 4 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 4 5 6 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 5 6 7 $numIntgrPts $LongSecTag $LongBeamTransfTag;
element nonlinearBeamColumn 6 7 8 $numIntgrPts $LongSecTag $LongBeamTransfTag;
# transverse
element nonlinearBeamColumn 7 1 5 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 8 2 6 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 9 3 7 $numIntgrPts $TransSecTag $TransBeamTransfTag;
element nonlinearBeamColumn 10 4 8 $numIntgrPts $TransSecTag $TransBeamTransfTag;
recorder Element -file $dataDir/Force.out -time -ele 1 2 globalForce;
recorder Node -file $dataDir/RSupport.out -time -node 1 5 4 8 -dof 3 reaction;
# DEFINE GRAVITY---------------------------------------------------
set g 9.81;
set UDL [expr $ALongSec*$dens*$g];
# Gravity loadcase
pattern Plain 1 Linear {
eleLoad -ele 1 2 3 4 5 6 -type -beamUniform -$UDL 0. 0.;
}
#ANALYSIS--------------
# Set up convergence criteria:
set Tol 1.0e-8; # convergence tolerance for test
set NumIter 10; # Set the number of iterations for convergence
test NormDispIncr $Tol $NumIter;
constraints Plain;
numberer RCM;
system BandGeneral;
algorithm Newton;
set LoadSteps 10;
set LoadIncr [expr 1.0/$LoadSteps];
integrator LoadControl $LoadIncr;
analysis Static;
analyze $LoadSteps; # applying gravity
#maintain constant gravity loads and reset time to zero
loadConst -time 0.0
puts "Model built"
-
rasoulazadi
- Posts: 14
- Joined: Mon Apr 18, 2011 10:27 pm
- Location: tehran
Re: Grillage analysis of Bridge
hi dear vesna
i need to model a concrete bridge in opensees
if you have an example about that please send me
thanks
best regards
i need to model a concrete bridge in opensees
if you have an example about that please send me
thanks
best regards