Hi everybody,
I have modeled a 3D 10-story steel building using forcebeamcolumn elements. the building has 3 bays in every direction ( x and Y). I want to perform an eigenvalue analysis. but, when the code runs the following message appears:
"Opensees has stopped working"
what is the problem?
Here is my code:
wipe;
model basic -ndm 3 -ndf 6
file mkdir data
# Create nodes
# ------------
node 1 0. 0. 0.
node 2 5. 0. 0.
node 3 10. 0. 0.
node 4 15. 0. 0.
node 5 0. 5. 0.
node 6 5. 5. 0.
node 7 10. 5. 0.
node 8 15. 5. 0.
node 9 0. 10. 0.
node 10 5. 10. 0.
node 11 10. 10. 0.
node 12 15. 10. 0.
node 13 0. 15. 0.
node 14 5. 15. 0.
node 15 10. 15. 0.
node 16 15. 15. 0.
node 17 0. 0. 3.2
node 18 5. 0. 3.2
node 19 10. 0. 3.2
node 20 15. 0. 3.2
node 21 0. 5. 3.2
node 22 5. 5. 3.2
node 23 10. 5. 3.2
node 24 15. 5. 3.2
node 25 0. 10. 3.2
node 26 5. 10. 3.2
node 27 10. 10. 3.2
node 28 15. 10. 3.2
node 29 0. 15. 3.2
node 30 5. 15. 3.2
node 31 10. 15. 3.2
node 32 15. 15. 3.2
node 33 0. 0. 6.4
node 34 5. 0. 6.4
node 35 10. 0. 6.4
node 36 15. 0. 6.4
node 37 0. 5. 6.4
node 38 5. 5. 6.4
node 39 10. 5. 6.4
node 40 15. 5. 6.4
node 41 0. 10. 6.4
node 42 5. 10. 6.4
node 43 10. 10. 6.4
node 44 15. 10. 6.4
node 45 0. 15. 6.4
node 46 5. 15. 6.4
node 47 10. 15. 6.4
node 48 15. 15. 6.4
node 49 0. 0. 9.6
node 50 5. 0. 9.6
node 51 10. 0. 9.6
node 52 15. 0. 9.6
node 53 0. 5. 9.6
node 54 5. 5. 9.6
node 55 10. 5. 9.6
node 56 15. 5. 9.6
node 57 0. 10. 9.6
node 58 5. 10. 9.6
node 59 10. 10. 9.6
node 60 15. 10. 9.6
node 61 0. 15. 9.6
node 62 5. 15. 9.6
node 63 10. 15. 9.6
node 64 15. 15. 9.6
node 65 0. 0. 12.8
node 66 5. 0. 12.8
node 67 10. 0. 12.8
node 68 15. 0. 12.8
node 69 0. 5. 12.8
node 70 5. 5. 12.8
node 71 10. 5. 12.8
node 72 15. 5. 12.8
node 73 0. 10. 12.8
node 74 5. 10. 12.8
node 75 10. 10. 12.8
node 76 15. 10. 12.8
node 77 0. 15. 12.8
node 78 5. 15. 12.8
node 79 10. 15. 12.8
node 80 15. 15. 12.8
node 81 0. 0. 16.0
node 82 5. 0. 16.0
node 83 10. 0. 16.0
node 84 15. 0. 16.0
node 85 0. 5. 16.0
node 86 5. 5. 16.0
node 87 10. 5. 16.0
node 88 15. 5. 16.0
node 89 0. 10. 16.0
node 90 5. 10. 16.0
node 91 10. 10. 16.0
node 92 15. 10. 16.0
node 93 0. 15. 16.0
node 94 5. 15. 16.0
node 95 10. 15. 16.0
node 96 15. 15. 16.0
node 97 0. 0. 19.2
node 98 5. 0. 19.2
node 99 10. 0. 19.2
node 100 15. 0. 19.2
node 101 0. 5. 19.2
node 102 5. 5. 19.2
node 103 10. 5. 19.2
node 104 15. 5. 19.2
node 105 0. 10. 19.2
node 106 5. 10. 19.2
node 107 10. 10. 19.2
node 108 15. 10. 19.2
node 109 0. 15. 19.2
node 110 5. 15. 19.2
node 111 10. 15. 19.2
node 112 15. 15. 19.2
node 113 0. 0. 22.4
node 114 5. 0. 22.4
node 115 10. 0. 22.4
node 116 15. 0. 22.4
node 117 0. 5. 22.4
node 118 5. 5. 22.4
node 119 10. 5. 22.4
node 120 15. 5. 22.4
node 121 0. 10. 22.4
node 122 5. 10. 22.4
node 123 10. 10. 22.4
node 124 15. 10. 22.4
node 125 0. 15. 22.4
node 126 5. 15. 22.4
node 127 10. 15. 22.4
node 128 15. 15. 22.4
node 129 0. 0. 25.6
node 130 5. 0. 25.6
node 131 10. 0. 25.6
node 132 15. 0. 25.6
node 133 0. 5. 25.6
node 134 5. 5. 25.6
node 135 10. 5. 25.6
node 136 15. 5. 25.6
node 137 0. 10. 25.6
node 138 5. 10. 25.6
node 139 10. 10. 25.6
node 140 15. 10. 25.6
node 141 0. 15. 25.6
node 142 5. 15. 25.6
node 143 10. 15. 25.6
node 144 15. 15. 25.6
node 145 0. 0. 28.8
node 146 5. 0. 28.8
node 147 10. 0. 28.8
node 148 15. 0. 28.8
node 149 0. 5. 28.8
node 150 5. 5. 28.8
node 151 10. 5. 28.8
node 152 15. 5. 28.8
node 153 0. 10. 28.8
node 154 5. 10. 28.8
node 155 10. 10. 28.8
node 156 15. 10. 28.8
node 157 0. 15. 28.8
node 158 5. 15. 28.8
node 159 10. 15. 28.8
node 160 15. 15. 28.8
node 161 0. 0. 32.
node 162 5. 0. 32.
node 163 10. 0. 32.
node 164 15. 0. 32.
node 165 0. 5. 32.
node 166 5. 5. 32.
node 167 10. 5. 32.
node 168 15. 5. 32.
node 169 0. 10. 32.
node 170 5. 10. 32.
node 171 10. 10. 32.
node 172 15. 10. 32.
node 173 0. 15. 32.
node 174 5. 15. 32.
node 175 10. 15. 32.
node 176 15. 15. 32.
node 200 7.5 7.5 3.2
node 201 7.5 7.5 6.4
node 202 7.5 7.5 9.6
node 203 7.5 7.5 12.8
node 204 7.5 7.5 16.
node 205 7.5 7.5 19.2
node 206 7.5 7.5 22.4
node 207 7.5 7.5 25.6
node 208 7.5 7.5 28.8
node 209 7.5 7.5 32.
for {set x 17} {$x<=176} {incr x 1} {
mass $x 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
}
mass 200 17329. 17329. 1.e-9 1.e-9 1.e-9 [expr (1./12)*(170000./9.81)*(pow(15,2)+pow(15,2))]
mass 201 17329. 17329. 1.e-9 1.e-9 1.e-9 [expr (1./12)*(170000./9.81)*(pow(15,2)+pow(15,2))]
mass 202 17329. 17329. 1.e-9 1.e-9 1.e-9 [expr (1./12)*(170000./9.81)*(pow(15,2)+pow(15,2))]
mass 203 17329. 17329. 1.e-9 1.e-9 1.e-9 [expr (1./12)*(170000./9.81)*(pow(15,2)+pow(15,2))]
mass 204 17329. 17329. 1.e-9 1.e-9 1.e-9 [expr (1./12)*(170000./9.81)*(pow(15,2)+pow(15,2))]
mass 205 17329. 17329. 1.e-9 1.e-9 1.e-9 [expr (1./12)*(170000./9.81)*(pow(15,2)+pow(15,2))]
mass 206 17329. 17329. 1.e-9 1.e-9 1.e-9 [expr (1./12)*(170000./9.81)*(pow(15,2)+pow(15,2))]
mass 207 17329. 17329. 1.e-9 1.e-9 1.e-9 [expr (1./12)*(170000./9.81)*(pow(15,2)+pow(15,2))]
mass 208 17329. 17329. 1.e-9 1.e-9 1.e-9 [expr (1./12)*(170000./9.81)*(pow(15,2)+pow(15,2))]
mass 209 17329. 17329. 1.e-9 1.e-9 1.e-9 [expr (1./12)*(170000./9.81)*(pow(15,2)+pow(15,2))]
fix 1 1 1 1 1 1 1
fix 2 1 1 1 1 1 1
fix 3 1 1 1 1 1 1
fix 4 1 1 1 1 1 1
fix 5 1 1 1 1 1 1
fix 6 1 1 1 1 1 1
fix 7 1 1 1 1 1 1
fix 8 1 1 1 1 1 1
fix 9 1 1 1 1 1 1
fix 10 1 1 1 1 1 1
fix 11 1 1 1 1 1 1
fix 12 1 1 1 1 1 1
fix 13 1 1 1 1 1 1
fix 14 1 1 1 1 1 1
fix 15 1 1 1 1 1 1
fix 16 1 1 1 1 1 1
fix 200 0 0 1 1 1 0
fix 201 0 0 1 1 1 0
fix 202 0 0 1 1 1 0
fix 203 0 0 1 1 1 0
fix 204 0 0 1 1 1 0
fix 205 0 0 1 1 1 0
fix 206 0 0 1 1 1 0
fix 207 0 0 1 1 1 0
fix 208 0 0 1 1 1 0
fix 209 0 0 1 1 1 0
rigidDiaphragm 3 200 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32
rigidDiaphragm 3 201 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48
rigidDiaphragm 3 202 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64
rigidDiaphragm 3 203 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80
rigidDiaphragm 3 204 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96
rigidDiaphragm 3 205 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112
rigidDiaphragm 3 206 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128
rigidDiaphragm 3 207 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144
rigidDiaphragm 3 208 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160
rigidDiaphragm 3 209 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176
set ColumnTransfTag 1
geomTransf PDelta $ColumnTransfTag 1 0 0
set BeamTransfTag 2
geomTransf Linear $BeamTransfTag 0 -1 0
set GirderTransfTag 3
geomTransf Linear $GirderTransfTag 1 0 0
uniaxialMaterial Steel01 1 24000000.0 21000000000.0 0.005
# SC5(column 45X3)
# set some paramaters
set h 0.45
set t .03
set y [expr $h/2.0]
section Fiber 1 {
# create top,bot,left,right
patch quad 1 30 4 [expr -$y] [expr $y] [expr -$y] [expr $y-$t] [expr $y] [expr $y-$t] [expr $y] [expr $y]
patch quad 1 30 4 [expr -$y] [expr -$y+$t] [expr -$y] [expr -$y] [expr $y] [expr -$y] [expr $y] [expr -$y+$t]
patch quad 1 4 25 [expr -$y] [expr $y-$t] [expr -$y] [expr -$y+$t] [expr -$y+$t] [expr -$y+$t] [expr -$y+$t] [expr $y-$t]
patch quad 1 4 25 [expr $y-$t] [expr $y-$t] [expr $y-$t] [expr -$y+$t] [expr $y] [expr -$y+$t] [expr $y] [expr $y-$t]
}
# SC4(column 40X2.5)
# set some paramaters
set h 0.4
set t .025
set y [expr $h/2.0]
section Fiber 2 {
# create top,bot,left,right
patch quad 1 30 4 [expr -$y] [expr $y] [expr -$y] [expr $y-$t] [expr $y] [expr $y-$t] [expr $y] [expr $y]
patch quad 1 30 4 [expr -$y] [expr -$y+$t] [expr -$y] [expr -$y] [expr $y] [expr -$y] [expr $y] [expr -$y+$t]
patch quad 1 4 25 [expr -$y] [expr $y-$t] [expr -$y] [expr -$y+$t] [expr -$y+$t] [expr -$y+$t] [expr -$y+$t] [expr $y-$t]
patch quad 1 4 25 [expr $y-$t] [expr $y-$t] [expr $y-$t] [expr -$y+$t] [expr $y] [expr -$y+$t] [expr $y] [expr $y-$t]
}
# SC3(column 35X2.5)
# set some paramaters
set h 0.35
set t .025
set y [expr $h/2.0]
section Fiber 3 {
# create top,bot,left,right
patch quad 1 30 4 [expr -$y] [expr $y] [expr -$y] [expr $y-$t] [expr $y] [expr $y-$t] [expr $y] [expr $y]
patch quad 1 30 4 [expr -$y] [expr -$y+$t] [expr -$y] [expr -$y] [expr $y] [expr -$y] [expr $y] [expr -$y+$t]
patch quad 1 4 25 [expr -$y] [expr $y-$t] [expr -$y] [expr -$y+$t] [expr -$y+$t] [expr -$y+$t] [expr -$y+$t] [expr $y-$t]
patch quad 1 4 25 [expr $y-$t] [expr $y-$t] [expr $y-$t] [expr -$y+$t] [expr $y] [expr -$y+$t] [expr $y] [expr $y-$t]
}
# SC2(column 30X2)
# set some paramaters
set h 0.3
set t .02
set y [expr $h/2.0]
section Fiber 4 {
# create top,bot,left,right
patch quad 1 30 4 [expr -$y] [expr $y] [expr -$y] [expr $y-$t] [expr $y] [expr $y-$t] [expr $y] [expr $y]
patch quad 1 30 4 [expr -$y] [expr -$y+$t] [expr -$y] [expr -$y] [expr $y] [expr -$y] [expr $y] [expr -$y+$t]
patch quad 1 4 25 [expr -$y] [expr $y-$t] [expr -$y] [expr -$y+$t] [expr -$y+$t] [expr -$y+$t] [expr -$y+$t] [expr $y-$t]
patch quad 1 4 25 [expr $y-$t] [expr $y-$t] [expr $y-$t] [expr -$y+$t] [expr $y] [expr -$y+$t] [expr $y] [expr $y-$t]
}
# SC1(column 25X1.5)
# set some paramaters
set h 0.25
set t .015
set y [expr $h/2.0]
section Fiber 5 {
# create top,bot,left,right
patch quad 1 30 4 [expr -$y] [expr $y] [expr -$y] [expr $y-$t] [expr $y] [expr $y-$t] [expr $y] [expr $y]
patch quad 1 30 4 [expr -$y] [expr -$y+$t] [expr -$y] [expr -$y] [expr $y] [expr -$y] [expr $y] [expr -$y+$t]
patch quad 1 4 25 [expr -$y] [expr $y-$t] [expr -$y] [expr -$y+$t] [expr -$y+$t] [expr -$y+$t] [expr -$y+$t] [expr $y-$t]
patch quad 1 4 25 [expr $y-$t] [expr $y-$t] [expr $y-$t] [expr -$y+$t] [expr $y] [expr -$y+$t] [expr $y] [expr $y-$t]
}
# SB5(beam 40 1 22.5 2)
# set some paramaters
set ht 0.4
set tw .01
set bf .225
set tf .02
set y1 [expr $ht/2.0]
set y2 [expr $bf/2.0]
set y3 [expr $y1-$tf]
set y4 [expr $tw/2.0]
section Fiber 7 {
# create bot flange,web, top flange
patch quad 1 4 20 [expr -$y1] [expr $y2] [expr -$y1] [expr -$y2] [expr -$y3] [expr -$y2] [expr -$y3] [expr $y2]
patch quad 1 30 2 [expr -$y3] [expr $y4] [expr -$y3] [expr -$y4] [expr $y3] [expr -$y4] [expr $y3] [expr $y4]
patch quad 1 4 20 [expr $y3] [expr $y2] [expr $y3] [expr -$y2] [expr $y1] [expr -$y2] [expr $y1] [expr $y2]
}
# SB4(beam 35 .8 22.5 2)
# set some paramaters
set ht 0.35
set tw .008
set bf .225
set tf .02
set y1 [expr $ht/2.0]
set y2 [expr $bf/2.0]
set y3 [expr $y1-$tf]
set y4 [expr $tw/2.0]
section Fiber 8 {
# create bot flange,web, top flange
patch quad 1 4 20 [expr -$y1] [expr $y2] [expr -$y1] [expr -$y2] [expr -$y3] [expr -$y2] [expr -$y3] [expr $y2]
patch quad 1 30 2 [expr -$y3] [expr $y4] [expr -$y3] [expr -$y4] [expr $y3] [expr -$y4] [expr $y3] [expr $y4]
patch quad 1 4 20 [expr $y3] [expr $y2] [expr $y3] [expr -$y2] [expr $y1] [expr -$y2] [expr $y1] [expr $y2]
}
# SB3(beam 30 0.8 20 1.5)
# set some paramaters
set ht 0.3
set tw .008
set bf .2
set tf .015
set y1 [expr $ht/2.0]
set y2 [expr $bf/2.0]
set y3 [expr $y1-$tf]
set y4 [expr $tw/2.0]
section Fiber 9 {
# create bot flange,web, top flange
patch quad 1 4 20 [expr -$y1] [expr $y2] [expr -$y1] [expr -$y2] [expr -$y3] [expr -$y2] [expr -$y3] [expr $y2]
patch quad 1 30 2 [expr -$y3] [expr $y4] [expr -$y3] [expr -$y4] [expr $y3] [expr -$y4] [expr $y3] [expr $y4]
patch quad 1 4 20 [expr $y3] [expr $y2] [expr $y3] [expr -$y2] [expr $y1] [expr -$y2] [expr $y1] [expr $y2]
}
# SB2(beam 30 0.8 15 1.5)
# set some paramaters
set ht 0.3
set tw .008
set bf .15
set tf .015
set y1 [expr $ht/2.0]
set y2 [expr $bf/2.0]
set y3 [expr $y1-$tf]
set y4 [expr $tw/2.0]
section Fiber 10 {
# create bot flange,web, top flange
patch quad 1 4 20 [expr -$y1] [expr $y2] [expr -$y1] [expr -$y2] [expr -$y3] [expr -$y2] [expr -$y3] [expr $y2]
patch quad 1 30 2 [expr -$y3] [expr $y4] [expr -$y3] [expr -$y4] [expr $y3] [expr -$y4] [expr $y3] [expr $y4]
patch quad 1 4 20 [expr $y3] [expr $y2] [expr $y3] [expr -$y2] [expr $y1] [expr -$y2] [expr $y1] [expr $y2]
}
# SB1 (beam 25 0.6 17.5 1.5)
# set some paramaters
set ht 0.25
set tw .006
set bf .175
set tf .015
set y1 [expr $ht/2.0]
set y2 [expr $bf/2.0]
set y3 [expr $y1-$tf]
set y4 [expr $tw/2.0]
section Fiber 11 {
# create bot flange,web, top flange
patch quad 1 4 20 [expr -$y1] [expr $y2] [expr -$y1] [expr -$y2] [expr -$y3] [expr -$y2] [expr -$y3] [expr $y2]
patch quad 1 30 2 [expr -$y3] [expr $y4] [expr -$y3] [expr -$y4] [expr $y3] [expr -$y4] [expr $y3] [expr $y4]
patch quad 1 4 20 [expr $y3] [expr $y2] [expr $y3] [expr -$y2] [expr $y1] [expr -$y2] [expr $y1] [expr $y2]
}
set e 21000000000.0
set nu 0.3
set G [expr $e/(2.*(1.+$nu))]
section Elastic 10 21000000000.0 .0504 .00148932 $G 0.536
section Elastic 20 21000000000.0 .0375 .0008828125 $G 0.533
section Elastic 30 21000000000.0 .0325 .0005755208 $G 0.538
section Elastic 40 21000000000.0 .0224 .0002941867 $G 0.536
section Elastic 50 21000000000.0 .0141 .0001303075 $G 0.532
section Elastic 70 21000000000.0 .0126 .00036408 $G 0.317
section Elastic 80 21000000000.0 .01148 .00026518567 $G 0.243
section Elastic 90 21000000000.0 .00816 .000135072 $G 0.294
section Elastic 100 21000000000.0 .00666 .0001045845 $G 0.360
section Elastic 110 21000000000.0 .00657 .00007790525 $G 0.228
# Define column elements
# ---------------------------------
# Geometry of column elements
#story 1 to 4
for {set x 0} {$x<=3} {incr x 1} {
# tag ndI ndJ nsecs secID transfTag
element forceBeamColumn [expr 1+16*$x] [expr 1+16*$x] [expr 1+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 2+16*$x] [expr 2+16*$x] [expr 2+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 3+16*$x] [expr 3+16*$x] [expr 3+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 4+16*$x] [expr 4+16*$x] [expr 4+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 5+16*$x] [expr 5+16*$x] [expr 5+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 6+16*$x] [expr 6+16*$x] [expr 6+16*($x+1)] $ColumnTransfTag "HingeRadau 1 .45 1 .45 10"
element forceBeamColumn [expr 7+16*$x] [expr 7+16*$x] [expr 7+16*($x+1)] $ColumnTransfTag "HingeRadau 1 .45 1 .45 10"
element forceBeamColumn [expr 8+16*$x] [expr 8+16*$x] [expr 8+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 9+16*$x] [expr 9+16*$x] [expr 9+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 10+16*$x] [expr 10+16*$x] [expr 10+16*($x+1)] $ColumnTransfTag "HingeRadau 1 .45 1 .45 10"
element forceBeamColumn [expr 11+16*$x] [expr 11+16*$x] [expr 11+16*($x+1)] $ColumnTransfTag "HingeRadau 1 .45 1 .45 10"
element forceBeamColumn [expr 12+16*$x] [expr 12+16*$x] [expr 12+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 13+16*$x] [expr 13+16*$x] [expr 13+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 14+16*$x] [expr 14+16*$x] [expr 14+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 15+16*$x] [expr 15+16*$x] [expr 15+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 16+16*$x] [expr 16+16*$x] [expr 16+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
}
#story 5
for {set x 4} {$x<=4} {incr x 1} {
# tag ndI ndJ nsecs secID transfTag
element forceBeamColumn [expr 1+16*$x] [expr 1+16*$x] [expr 1+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 2+16*$x] [expr 2+16*$x] [expr 2+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 3+16*$x] [expr 3+16*$x] [expr 3+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 4+16*$x] [expr 4+16*$x] [expr 4+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 5+16*$x] [expr 5+16*$x] [expr 5+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 6+16*$x] [expr 6+16*$x] [expr 6+16*($x+1)] $ColumnTransfTag "HingeRadau 1 .45 1 .45 10"
element forceBeamColumn [expr 7+16*$x] [expr 7+16*$x] [expr 7+16*($x+1)] $ColumnTransfTag "HingeRadau 1 .45 1 .45 10"
element forceBeamColumn [expr 8+16*$x] [expr 8+16*$x] [expr 8+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 9+16*$x] [expr 9+16*$x] [expr 9+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 10+16*$x] [expr 10+16*$x] [expr 10+16*($x+1)] $ColumnTransfTag "HingeRadau 1 .45 1 .45 10"
element forceBeamColumn [expr 11+16*$x] [expr 11+16*$x] [expr 11+16*($x+1)] $ColumnTransfTag "HingeRadau 1 .45 1 .45 10"
element forceBeamColumn [expr 12+16*$x] [expr 12+16*$x] [expr 12+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 13+16*$x] [expr 13+16*$x] [expr 13+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 14+16*$x] [expr 14+16*$x] [expr 14+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 15+16*$x] [expr 15+16*$x] [expr 15+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 16+16*$x] [expr 16+16*$x] [expr 16+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
}
#story 6
for {set x 5} {$x<=5} {incr x 1} {
# tag ndI ndJ nsecs secID transfTag
element forceBeamColumn [expr 1+16*$x] [expr 1+16*$x] [expr 1+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 2+16*$x] [expr 2+16*$x] [expr 2+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 3+16*$x] [expr 3+16*$x] [expr 3+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 4+16*$x] [expr 4+16*$x] [expr 4+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 5+16*$x] [expr 5+16*$x] [expr 5+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 6+16*$x] [expr 6+16*$x] [expr 6+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 7+16*$x] [expr 7+16*$x] [expr 7+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 8+16*$x] [expr 8+16*$x] [expr 8+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 9+16*$x] [expr 9+16*$x] [expr 9+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 10+16*$x] [expr 10+16*$x] [expr 10+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 11+16*$x] [expr 11+16*$x] [expr 11+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 12+16*$x] [expr 12+16*$x] [expr 12+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 13+16*$x] [expr 13+16*$x] [expr 13+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 14+16*$x] [expr 14+16*$x] [expr 14+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 15+16*$x] [expr 15+16*$x] [expr 15+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 16+16*$x] [expr 16+16*$x] [expr 16+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
}
#story 7
for {set x 6} {$x<=6} {incr x 1} {
# tag ndI ndJ nsecs secID transfTag
element forceBeamColumn [expr 1+16*$x] [expr 1+16*$x] [expr 1+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 2+16*$x] [expr 2+16*$x] [expr 2+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 3+16*$x] [expr 3+16*$x] [expr 3+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 4+16*$x] [expr 4+16*$x] [expr 4+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 5+16*$x] [expr 5+16*$x] [expr 5+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 6+16*$x] [expr 6+16*$x] [expr 6+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 7+16*$x] [expr 7+16*$x] [expr 7+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 8+16*$x] [expr 8+16*$x] [expr 8+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 9+16*$x] [expr 9+16*$x] [expr 9+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 10+16*$x] [expr 10+16*$x] [expr 10+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 11+16*$x] [expr 11+16*$x] [expr 11+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 12+16*$x] [expr 12+16*$x] [expr 12+16*($x+1)] $ColumnTransfTag "HingeRadau 2 .40 2 .40 20"
element forceBeamColumn [expr 13+16*$x] [expr 13+16*$x] [expr 13+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 14+16*$x] [expr 14+16*$x] [expr 14+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 15+16*$x] [expr 15+16*$x] [expr 15+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 16+16*$x] [expr 16+16*$x] [expr 16+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
}
#story 8
for {set x 7} {$x<=7} {incr x 1} {
# tag ndI ndJ nsecs secID transfTag
element forceBeamColumn [expr 1+16*$x] [expr 1+16*$x] [expr 1+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 2+16*$x] [expr 2+16*$x] [expr 2+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 3+16*$x] [expr 3+16*$x] [expr 3+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 4+16*$x] [expr 4+16*$x] [expr 4+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 5+16*$x] [expr 5+16*$x] [expr 5+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 6+16*$x] [expr 6+16*$x] [expr 6+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 7+16*$x] [expr 7+16*$x] [expr 7+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 8+16*$x] [expr 8+16*$x] [expr 8+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 9+16*$x] [expr 9+16*$x] [expr 9+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 10+16*$x] [expr 10+16*$x] [expr 10+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 11+16*$x] [expr 11+16*$x] [expr 11+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 12+16*$x] [expr 12+16*$x] [expr 12+16*($x+1)] $ColumnTransfTag "HingeRadau 3 .35 3 .35 30"
element forceBeamColumn [expr 13+16*$x] [expr 13+16*$x] [expr 13+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 14+16*$x] [expr 14+16*$x] [expr 14+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 15+16*$x] [expr 15+16*$x] [expr 15+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 16+16*$x] [expr 16+16*$x] [expr 16+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
}
#story 9 & 10
for {set x 8} {$x<=9} {incr x 1} {
# tag ndI ndJ nsecs secID transfTag
element forceBeamColumn [expr 1+16*$x] [expr 1+16*$x] [expr 1+16*($x+1)] $ColumnTransfTag "HingeRadau 5 .25 5 .25 50"
element forceBeamColumn [expr 2+16*$x] [expr 2+16*$x] [expr 2+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 3+16*$x] [expr 3+16*$x] [expr 3+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 4+16*$x] [expr 4+16*$x] [expr 4+16*($x+1)] $ColumnTransfTag "HingeRadau 5 .25 5 .25 50"
element forceBeamColumn [expr 5+16*$x] [expr 5+16*$x] [expr 5+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 6+16*$x] [expr 6+16*$x] [expr 6+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 7+16*$x] [expr 7+16*$x] [expr 7+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 8+16*$x] [expr 8+16*$x] [expr 8+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 9+16*$x] [expr 9+16*$x] [expr 9+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 10+16*$x] [expr 10+16*$x] [expr 10+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 11+16*$x] [expr 11+16*$x] [expr 11+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 12+16*$x] [expr 12+16*$x] [expr 12+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 13+16*$x] [expr 13+16*$x] [expr 13+16*($x+1)] $ColumnTransfTag "HingeRadau 5 .25 5 .25 50"
element forceBeamColumn [expr 14+16*$x] [expr 14+16*$x] [expr 14+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 15+16*$x] [expr 15+16*$x] [expr 15+16*($x+1)] $ColumnTransfTag "HingeRadau 4 .30 4 .30 40"
element forceBeamColumn [expr 16+16*$x] [expr 16+16*$x] [expr 16+16*($x+1)] $ColumnTransfTag "HingeRadau 5 .25 5 .25 50"
}
#Beams
#1-5
for {set x 0} {$x<=4} {incr x 1} {
# tag ndI ndJ nsecs secID transfTag
element forceBeamColumn [expr 161+12*$x] [expr 17+16*$x] [expr 18+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 162+12*$x] [expr 18+16*$x] [expr 19+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 163+12*$x] [expr 19+16*$x] [expr 20+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 164+12*$x] [expr 21+16*$x] [expr 22+16*$x] $BeamTransfTag "HingeRadau 8 .35 8 .35 80"
element forceBeamColumn [expr 165+12*$x] [expr 22+16*$x] [expr 23+16*$x] $BeamTransfTag "HingeRadau 8 .35 8 .35 80"
element forceBeamColumn [expr 166+12*$x] [expr 23+16*$x] [expr 24+16*$x] $BeamTransfTag "HingeRadau 8 .35 8 .35 80"
element forceBeamColumn [expr 167+12*$x] [expr 25+16*$x] [expr 26+16*$x] $BeamTransfTag "HingeRadau 8 .35 8 .35 80"
element forceBeamColumn [expr 168+12*$x] [expr 26+16*$x] [expr 27+16*$x] $BeamTransfTag "HingeRadau 8 .35 8 .35 80"
element forceBeamColumn [expr 169+12*$x] [expr 27+16*$x] [expr 28+16*$x] $BeamTransfTag "HingeRadau 8 .35 8 .35 80"
element forceBeamColumn [expr 170+12*$x] [expr 29+16*$x] [expr 30+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 171+12*$x] [expr 30+16*$x] [expr 31+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 172+12*$x] [expr 31+16*$x] [expr 32+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
}
#6-7
for {set x 5} {$x<=6} {incr x 1} {
# tag ndI ndJ nsecs secID transfTag
element forceBeamColumn [expr 161+12*$x] [expr 17+16*$x] [expr 18+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 162+12*$x] [expr 18+16*$x] [expr 19+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 163+12*$x] [expr 19+16*$x] [expr 20+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 164+12*$x] [expr 21+16*$x] [expr 22+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 165+12*$x] [expr 22+16*$x] [expr 23+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 166+12*$x] [expr 23+16*$x] [expr 24+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 167+12*$x] [expr 25+16*$x] [expr 26+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 168+12*$x] [expr 26+16*$x] [expr 27+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 169+12*$x] [expr 27+16*$x] [expr 28+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 170+12*$x] [expr 29+16*$x] [expr 30+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 171+12*$x] [expr 30+16*$x] [expr 31+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 172+12*$x] [expr 31+16*$x] [expr 32+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
}
#8
for {set x 7} {$x<=7} {incr x 1} {
# tag ndI ndJ nsecs secID transfTag
element forceBeamColumn [expr 161+12*$x] [expr 17+16*$x] [expr 18+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
element forceBeamColumn [expr 162+12*$x] [expr 18+16*$x] [expr 19+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
element forceBeamColumn [expr 163+12*$x] [expr 19+16*$x] [expr 20+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
element forceBeamColumn [expr 164+12*$x] [expr 21+16*$x] [expr 22+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 165+12*$x] [expr 22+16*$x] [expr 23+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 166+12*$x] [expr 23+16*$x] [expr 24+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 167+12*$x] [expr 25+16*$x] [expr 26+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 168+12*$x] [expr 26+16*$x] [expr 27+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 169+12*$x] [expr 27+16*$x] [expr 28+16*$x] $BeamTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 170+12*$x] [expr 29+16*$x] [expr 30+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
element forceBeamColumn [expr 171+12*$x] [expr 30+16*$x] [expr 31+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
element forceBeamColumn [expr 172+12*$x] [expr 31+16*$x] [expr 32+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
}
#9-10
for {set x 8} {$x<=9} {incr x 1} {
# tag ndI ndJ nsecs secID transfTag
element forceBeamColumn [expr 161+12*$x] [expr 17+16*$x] [expr 18+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
element forceBeamColumn [expr 162+12*$x] [expr 18+16*$x] [expr 19+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
element forceBeamColumn [expr 163+12*$x] [expr 19+16*$x] [expr 20+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
element forceBeamColumn [expr 164+12*$x] [expr 21+16*$x] [expr 22+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
element forceBeamColumn [expr 165+12*$x] [expr 22+16*$x] [expr 23+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
element forceBeamColumn [expr 166+12*$x] [expr 23+16*$x] [expr 24+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
element forceBeamColumn [expr 167+12*$x] [expr 25+16*$x] [expr 26+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
element forceBeamColumn [expr 168+12*$x] [expr 26+16*$x] [expr 27+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
element forceBeamColumn [expr 169+12*$x] [expr 27+16*$x] [expr 28+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
element forceBeamColumn [expr 170+12*$x] [expr 29+16*$x] [expr 30+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
element forceBeamColumn [expr 171+12*$x] [expr 30+16*$x] [expr 31+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
element forceBeamColumn [expr 172+12*$x] [expr 31+16*$x] [expr 32+16*$x] $BeamTransfTag "HingeRadau 11 .25 11 .25 110"
}
#Girders
#1-8
for {set x 0} {$x<=7} {incr x 1} {
# tag ndI ndJ nsecs secID transfTag
element forceBeamColumn [expr 281+12*$x] [expr 17+16*$x] [expr 21+16*$x] $GirderTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 282+12*$x] [expr 18+16*$x] [expr 22+16*$x] $GirderTransfTag "HingeRadau 7 .40 7 .40 70"
element forceBeamColumn [expr 283+12*$x] [expr 19+16*$x] [expr 23+16*$x] $GirderTransfTag "HingeRadau 7 .40 7 .40 70"
element forceBeamColumn [expr 284+12*$x] [expr 20+16*$x] [expr 24+16*$x] $GirderTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 285+12*$x] [expr 21+16*$x] [expr 25+16*$x] $GirderTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 286+12*$x] [expr 22+16*$x] [expr 26+16*$x] $GirderTransfTag "HingeRadau 7 .40 7 .40 70"
element forceBeamColumn [expr 287+12*$x] [expr 23+16*$x] [expr 27+16*$x] $GirderTransfTag "HingeRadau 7 .40 7 .40 70"
element forceBeamColumn [expr 288+12*$x] [expr 24+16*$x] [expr 28+16*$x] $GirderTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 289+12*$x] [expr 25+16*$x] [expr 29+16*$x] $GirderTransfTag "HingeRadau 9 .30 9 .30 90"
element forceBeamColumn [expr 290+12*$x] [expr 26+16*$x] [expr 30+16*$x] $GirderTransfTag "HingeRadau 7 .40 7 .40 70"
element forceBeamColumn [expr 291+12*$x] [expr 27+16*$x] [expr 31+16*$x] $GirderTransfTag "HingeRadau 7 .40 7 .40 70"
element forceBeamColumn [expr 292+12*$x] [expr 28+16*$x] [expr 32+16*$x] $GirderTransfTag "HingeRadau 9 .30 9 .30 90"
}
# 9,10
for {set x 8} {$x<=9} {incr x 1} {
# tag ndI ndJ nsecs secID transfTag
element forceBeamColumn [expr 281+12*$x] [expr 17+16*$x] [expr 21+16*$x] $GirderTransfTag "HingeRadau 10 .30 10 .30 100"
element forceBeamColumn [expr 282+12*$x] [expr 18+16*$x] [expr 22+16*$x] $GirderTransfTag "HingeRadau 8 .35 8 .35 80"
element forceBeamColumn [expr 283+12*$x] [expr 19+16*$x] [expr 23+16*$x] $GirderTransfTag "HingeRadau 8 .35 8 .35 80"
element forceBeamColumn [expr 284+12*$x] [expr 20+16*$x] [expr 24+16*$x] $GirderTransfTag "HingeRadau 10 .30 10 .30 100"
element forceBeamColumn [expr 285+12*$x] [expr 21+16*$x] [expr 25+16*$x] $GirderTransfTag "HingeRadau 10 .30 10 .30 100"
element forceBeamColumn [expr 286+12*$x] [expr 22+16*$x] [expr 26+16*$x] $GirderTransfTag "HingeRadau 8 .35 8 .35 80"
element forceBeamColumn [expr 287+12*$x] [expr 23+16*$x] [expr 27+16*$x] $GirderTransfTag "HingeRadau 8 .35 8 .35 80"
element forceBeamColumn [expr 288+12*$x] [expr 24+16*$x] [expr 28+16*$x] $GirderTransfTag "HingeRadau 10 .30 10 .30 100"
element forceBeamColumn [expr 289+12*$x] [expr 25+16*$x] [expr 29+16*$x] $GirderTransfTag "HingeRadau 10 .30 10 .30 100"
element forceBeamColumn [expr 290+12*$x] [expr 26+16*$x] [expr 30+16*$x] $GirderTransfTag "HingeRadau 8 .35 8 .35 80"
element forceBeamColumn [expr 291+12*$x] [expr 27+16*$x] [expr 31+16*$x] $GirderTransfTag "HingeRadau 8 .35 8 .35 80"
element forceBeamColumn [expr 292+12*$x] [expr 28+16*$x] [expr 32+16*$x] $GirderTransfTag "HingeRadau 10 .30 10 .30 100"
}
# Define gravity loads
# -------------------------------------
set p [expr -4650.]; # (1.1*650+1.1*200)*5=4675 kg/m
# Create a Plain load pattern with a Linear TimeSeries
pattern Plain 1 Linear {
# Create distributed loads
for {set x 281} {$x<=400} {incr x 1} {
eleLoad -ele $x -type -beamUniform $p
}
}
# Create the system of equation, a sparse solver with partial pivoting
system BandGeneral
# Create the constraint handler, the transformation method
constraints Transformation
# Create the DOF numberer, the reverse Cuthill-McKee algorithm
numberer RCM
# Create the convergence test, the norm of the residual with a tolerance of
# 1e-12 and a max number of iterations of 10
test NormDispIncr 1.0e-12 10 3
# Create the solution algorithm, a Newton-Raphson algorithm
algorithm Newton
# Create the integration scheme, the LoadControl scheme using steps of 0.1
integrator LoadControl 0.1
# Create the analysis object
analysis Static
analyze 10
loadConst -time 0.0;
#recorder Node -file eigen1.out -node 13 -dof 1 2 6 "eigen1"
#recorder Node -file eigen2.out -node 13 -dof 1 2 6 "eigen2"
#recorder Node -file eigen3.out -node 13 -dof 1 2 6 "eigen3"
set eig [eigen 3]
set pi [expr 4*atan(1.)]
puts "First Mode Period is:[expr 2*$pi/pow([lindex $eig 0],0.5)]"
puts "Second Mode Period is:[expr 2*$pi/pow([lindex $eig 1],0.5)]"
puts "Third Mode Period is:[expr 2*$pi/pow([lindex $eig 2],0.5)]"
integrator LoadControl 0 1 0 0
test EnergyIncr 1.0e-10 100 0
algorithm Newton
numberer RCM
constraints Transformation
system ProfileSPD
analysis Static
analyze 1
Opensees has stopped working!!!
Moderators: silvia, selimgunay, Moderators
Re: Opensees has stopped working!!!
i will update the code .. the problem is the Elastic Scetion definition .. this is a 3d problem, you are providing 2d properties
Re: Opensees has stopped working!!!
thanks dear fmk for your answer