please help me

[zdzdzdzd]

ArpackSolver::Error with _saupd info = -1
N must be positive.
WARNING DirectIntegrationAnalysis::eigen() - EigenSOE failed in solve()


node 1 255624 -15729.4 5682.7
node 2 250212 -7864.7 10805.3
node 3 244582 -15770.2 4825.2
node 4 244800 -15241.5 15928
node 5 244376 -16270.6 -5682.7
node 6 239793 -7759.57 10097.2
node 7 239488 -8135.3 -816.3
node 8 234600 -15807.1 4050.1
node 9 244800 2.76486e-10 15928
node 10 250212 7864.7 10805.3
node 11 234600 6.18456e-11 4050.1
node 12 229500 -7903.55 8781.2
node 13 234600 -14813.4 24918.9
node 14 239793 7759.57 10097.2
node 15 239488 8135.3 -816.3
node 16 255624 15729.4 5682.7
node 17 224400 -15356.6 13512.3
node 18 244582 15770.2 4825.2
node 19 244800 15241.5 15928
node 20 229500 7903.55 8781.2
node 21 244376 16270.6 -5682.7
node 22 224400 6.00267e-11 13512.3
node 23 224400 -7435.24 23720.4
node 24 234600 15807.1 4050.1
node 25 224400 -14414.2 33302
node 26 234600 14813.4 24918.9
node 27 224400 7435.24 23720.4
node 28 224400 7.00311e-11 33302
node 29 214200 -14937.8 22306.7
node 30 224400 15356.6 13512.3
node 31 224400 14414.2 33302
node 32 214200 14937.8 22306.7
node 33 214200 -14042.4 41110.1
node 34 204000 -14549 30470.7
node 35 204000 5.09317e-11 30470.7
node 36 214200 14042.4 41110.1
node 37 204000 -7054.94 39692.4
node 38 204000 14549 30470.7
node 39 204000 7054.94 39692.4
node 40 204000 -13696.5 48373.7
node 41 204000 7.27596e-11 48373.7
node 42 193800 -14188.6 38039.4
node 43 204000 13696.5 48373.7
node 44 193800 14188.6 38039.4
node 45 193800 -13375.2 55121.2
node 46 183600 -13855 45045
node 47 183600 2.00089e-11 45045
node 48 193800 13375.2 55121.2
node 49 183600 -6727.43 53447.9
node 50 183600 13855 45045
node 51 183600 6727.43 53447.9
node 52 183600 -13077.2 61379.1

fix 1 1 1 1 1 1 1
fix 5 1 1 1 1 1 1
fix 16 1 1 1 1 1 1
fix 21 1 1 1 1 1 1
fix 276 1 1 1 1 1 1
fix 280 1 1 1 1 1 1
fix 285 1 1 1 1 1 1
fix 286 1 1 1 1 1 1

# --------------------------------------------------------------------------------------------------------------
# M A S S E S
# --------------------------------------------------------------------------------------------------------------

# Mass Definition : mass $NodeTag $(ndf nodal mass values corresponding to each DOF)

mass 1 0 0 -190619 0 0 0
mass 1 0 0 -9390 0 0 0
mass 1 0 0 -190619 0 0 0
mass 1 0 0 -13047.7 0 0 0
mass 2 0 0 -13047.7 0 0 0
mass 3 0 0 -9390 0 0 0
mass 3 0 0 -8488.27 0 0 0
mass 3 0 0 -8920.41 0 0 0
mass 3 0 0 -9425.43 0 0 0
mass 4 0 0 -2540.45 0 0 0
mass 4 0 0 -20070.9 0 0 0
mass 4 0 0 -173893 0 0 0
mass 4 0 0 -190619 0 0 0
mass 4 0 0 -2540.45 0 0 0
mass 4 0 0 -190619 0 0 0
mass 4 0 0 -173893 0 0 0
mass 4 0 0 -9425.43 0 0 0
mass 4 0 0 -13047.7 0 0 0
mass 5 0 0 -176588 0 0 0
mass 5 0 0 -12844.3 0 0 0
mass 5 0 0 -176588 0 0 0
mass 5 0 0 -8920.41 0 0 0
mass 6 0 0 -5202.36 0 0 0
mass 6 0 0 -5202.36 0 0 0
mass 6 0 0 -2634.73 0 0 0
mass 6 0 0 -5016.21 0 0 0
mass 6 0 0 -2634.73 0 0 0
mass 6 0 0 -2540.45 0 0 0
mass 6 0 0 -5016.21 0 0 0
mass 6 0 0 -2540.45 0 0 0
mass 7 0 0 -12844.3 0 0 0
mass 8 0 0 -178098 0 0 0
mass 8 0 0 -178098 0 0 0
mass 8 0 0 -12513.4 0 0 0
mass 8 0 0 -176588 0 0 0
mass 8 0 0 -20336.7 0 0 0
mass 8 0 0 -176588 0 0 0
mass 8 0 0 -5202.36 0 0 0
mass 8 0 0 -5202.36 0 0 0
mass 8 0 0 -20070.9 0 0 0
mass 8 0 0 -8488.27 0 0 0
mass 8 0 0 -14424.7 0 0 0
mass 9 0 0 -12107.5 0 0 0
mass 9 0 0 -13047.7 0 0 0
mass 9 0 0 -5016.21 0 0 0
mass 9 0 0 -5016.21 0 0 0
mass 9 0 0 -5016.21 0 0 0
mass 9 0 0 -13047.7 0 0 0
mass 9 0 0 -12107.5 0 0 0
mass 9 0 0 -5016.21 0 0 0
mass 9 0 0 -13047.7 0 0 0
mass 9 0 0 -13047.7 0 0 0
mass 10 0 0 -13047.7 0 0 0
mass 11 0 0 -12844.3 0 0 0
mass 11 0 0 -2634.73 0 0 0
mass 11 0 0 -12844.3 0 0 0
mass 11 0 0 -2634.73 0 0 0
mass 11 0 0 -2634.73 0 0 0
mass 11 0 0 -2634.73 0 0 0
mass 11 0 0 -12513.4 0 0 0
mass 11 0 0 -12513.4 0 0 0
mass 12 0 0 -14424.7 0 0 0
mass 13 0 0 -20336.7 0 0 0
mass 13 0 0 -168695 0 0 0
mass 13 0 0 -168695 0 0 0
mass 13 0 0 -173893 0 0 0
mass 13 0 0 -173893 0 0 0
mass 13 0 0 -12107.5 0 0 0
mass 13 0 0 -12981.2 0 0 0
mass 13 0 0 -11948.3 0 0 0
mass 14 0 0 -2540.45 0 0 0
mass 14 0 0 -5016.21 0 0 0
mass 14 0 0 -2634.73 0 0 0
mass 14 0 0 -5016.21 0 0 0
mass 14 0 0 -5202.36 0 0 0
mass 14 0 0 -5202.36 0 0 0
mass 14 0 0 -2634.73 0 0 0
mass 14 0 0 -2540.45 0 0 0
mass 15 0 0 -12844.3 0 0 0
mass 16 0 0 -190619 0 0 0
mass 16 0 0 -190619 0 0 0
mass 16 0 0 -9390 0 0 0
mass 16 0 0 -13047.7 0 0 0
mass 17 0 0 -12981.2 0 0 0
mass 17 0 0 -19285.1 0 0 0
mass 17 0 0 -7831.18 0 0 0
mass 17 0 0 -178098 0 0 0
mass 17 0 0 -14506.4 0 0 0
mass 17 0 0 -178098 0 0 0
mass 17 0 0 -172341 0 0 0
mass 17 0 0 -172341 0 0 0
mass 17 0 0 -7831.18 0 0 0
mass 18 0 0 -9425.43 0 0 0
mass 18 0 0 -8488.27 0 0 0
mass 18 0 0 -9390 0 0 0
mass 18 0 0 -8920.41 0 0 0
mass 19 0 0 -173893 0 0 0
mass 19 0 0 -9425.43 0 0 0
mass 19 0 0 -173893 0 0 0
mass 19 0 0 -2540.45 0 0 0

mass 264 0 0 -12723 0 0 0
mass 264 0 0 -164122 0 0 0
mass 264 0 0 -13616.1 0 0 0
mass 265 0 0 -14424.7 0 0 0
mass 266 0 0 -14424.7 0 0 0
mass 267 0 0 -178098 0 0 0
mass 267 0 0 -20070.9 0 0 0
mass 267 0 0 -5202.36 0 0 0
mass 267 0 0 -176588 0 0 0
mass 267 0 0 -20336.7 0 0 0
mass 267 0 0 -14424.7 0 0 0
mass 267 0 0 -178098 0 0 0
mass 267 0 0 -176588 0 0 0
mass 267 0 0 -12513.4 0 0 0
mass 267 0 0 -5202.36 0 0 0
mass 267 0 0 -8488.27 0 0 0
mass 268 0 0 -2634.73 0 0 0
mass 268 0 0 -2634.73 0 0 0
mass 268 0 0 -2634.73 0 0 0
mass 268 0 0 -2634.73 0 0 0
mass 268 0 0 -12513.4 0 0 0
mass 268 0 0 -12844.3 0 0 0
mass 268 0 0 -12844.3 0 0 0
mass 268 0 0 -12513.4 0 0 0
mass 269 0 0 -20336.7 0 0 0
mass 269 0 0 -12107.5 0 0 0
mass 269 0 0 -168695 0 0 0
mass 269 0 0 -173893 0 0 0
mass 269 0 0 -173893 0 0 0
mass 269 0 0 -11948.3 0 0 0
mass 269 0 0 -168695 0 0 0
mass 269 0 0 -12981.2 0 0 0
mass 270 0 0 -176588 0 0 0
mass 270 0 0 -20336.7 0 0 0
mass 270 0 0 -20070.9 0 0 0
mass 270 0 0 -5202.36 0 0 0
mass 270 0 0 -176588 0 0 0
mass 270 0 0 -178098 0 0 0
mass 270 0 0 -178098 0 0 0
mass 270 0 0 -12513.4 0 0 0
mass 270 0 0 -14424.7 0 0 0
mass 270 0 0 -8488.27 0 0 0
mass 270 0 0 -5202.36 0 0 0
mass 271 0 0 -173893 0 0 0
mass 271 0 0 -168695 0 0 0
mass 271 0 0 -12107.5 0 0 0
mass 271 0 0 -11948.3 0 0 0
mass 271 0 0 -20336.7 0 0 0
mass 271 0 0 -168695 0 0 0
mass 271 0 0 -12981.2 0 0 0
mass 271 0 0 -173893 0 0 0
mass 272 0 0 -12844.3 0 0 0
mass 273 0 0 -5016.21 0 0 0
mass 273 0 0 -2634.73 0 0 0
mass 273 0 0 -2540.45 0 0 0
mass 273 0 0 -2634.73 0 0 0
mass 273 0 0 -5202.36 0 0 0
mass 273 0 0 -5202.36 0 0 0
mass 273 0 0 -5016.21 0 0 0
mass 273 0 0 -2540.45 0 0 0
mass 274 0 0 -12844.3 0 0 0
mass 275 0 0 -2634.73 0 0 0
mass 275 0 0 -5202.36 0 0 0
mass 275 0 0 -2634.73 0 0 0
mass 275 0 0 -2540.45 0 0 0
mass 275 0 0 -2540.45 0 0 0
mass 275 0 0 -5016.21 0 0 0
mass 275 0 0 -5202.36 0 0 0
mass 275 0 0 -5016.21 0 0 0
mass 276 0 0 -176588 0 0 0
mass 276 0 0 -8920.41 0 0 0
mass 276 0 0 -12844.3 0 0 0
mass 276 0 0 -176588 0 0 0
mass 277 0 0 -9425.43 0 0 0
mass 277 0 0 -8920.41 0 0 0
mass 277 0 0 -9390 0 0 0
mass 277 0 0 -8488.27 0 0 0
mass 278 0 0 -2540.45 0 0 0
mass 278 0 0 -9425.43 0 0 0
mass 278 0 0 -173893 0 0 0
mass 278 0 0 -190619 0 0 0
mass 278 0 0 -20070.9 0 0 0
mass 278 0 0 -190619 0 0 0
mass 278 0 0 -13047.7 0 0 0
mass 278 0 0 -2540.45 0 0 0
mass 278 0 0 -173893 0 0 0
mass 279 0 0 -5016.21 0 0 0
mass 279 0 0 -5016.21 0 0 0
mass 279 0 0 -13047.7 0 0 0
mass 279 0 0 -12107.5 0 0 0
mass 279 0 0 -5016.21 0 0 0
mass 279 0 0 -12107.5 0 0 0
mass 279 0 0 -13047.7 0 0 0
mass 279 0 0 -13047.7 0 0 0
mass 279 0 0 -13047.7 0 0 0
mass 279 0 0 -5016.21 0 0 0
mass 280 0 0 -176588 0 0 0
mass 280 0 0 -176588 0 0 0
mass 280 0 0 -8920.41 0 0 0
mass 280 0 0 -12844.3 0 0 0
mass 281 0 0 -8920.41 0 0 0
mass 281 0 0 -9390 0 0 0
mass 281 0 0 -9425.43 0 0 0
mass 281 0 0 -8488.27 0 0 0
mass 282 0 0 -13047.7 0 0 0
mass 282 0 0 -173893 0 0 0
mass 282 0 0 -20070.9 0 0 0
mass 282 0 0 -2540.45 0 0 0
mass 282 0 0 -2540.45 0 0 0
mass 282 0 0 -173893 0 0 0
mass 282 0 0 -9425.43 0 0 0
mass 282 0 0 -190619 0 0 0
mass 282 0 0 -190619 0 0 0
mass 283 0 0 -13047.7 0 0 0
mass 284 0 0 -13047.7 0 0 0
mass 285 0 0 -13047.7 0 0 0
mass 285 0 0 -190619 0 0 0
mass 285 0 0 -9390 0 0 0
mass 285 0 0 -190619 0 0 0
mass 286 0 0 -13047.7 0 0 0
mass 286 0 0 -9390 0 0 0
mass 286 0 0 -190619 0 0 0
mass 286 0 0 -190619 0 0 0

# --------------------------------------------------------------------------------------------------------------
# F O R C E - B A S E D B E A M - C O L U M N E L E M E N T S
# --------------------------------------------------------------------------------------------------------------

# Geometric Transformation

geomTransf Linear 1 -1 0 0; # vertical
geomTransf Linear 2 0 0 1; # non-vertical
geomTransf PDelta 3 -1 0 0; # vertical
geomTransf PDelta 4 0 1 0; # non-vertical
geomTransf Corotational 5 -1 0 0; # vertical
geomTransf Corotational 6 0 1 0; # non-vertica

# Sections Definition used by forceBeamColumn Elements
# (飊ly if they have not already been defined on this model domain)

section Fiber 174 -GJ 1e10 {
uniaxialMaterial ConcreteCM 121 -70000 -0.0027 4.1e+07 7 1.035 4600 0.00011 7 10000
set MatTag 121
source "../Scripts/Fiber_Custom-10_1.tcl"
uniaxialMaterial Steel01 163 500000 2.1e+08 0.01 0 1 0 1
set MatTag 163
source "../Scripts/Fiber_Custom-10_2.tcl"
}
section Fiber 173 -GJ 1e10 {
set MatTag 121
source "../Scripts/Fiber_Custom-9_1.tcl"
set MatTag 163
source "../Scripts/Fiber_Custom-9_2.tcl"
}
section Fiber 170 -GJ 1e10 {
set MatTag 121
source "../Scripts/Fiber_Custom-6_1.tcl"
set MatTag 163
source "../Scripts/Fiber_Custom-6_2.tcl"
}
section Fiber 168 -GJ 1e10 {
set MatTag 121
source "../Scripts/Fiber_Custom-4_1.tcl"
set MatTag 163
source "../Scripts/Fiber_Custom-4_2.tcl"
}
section Fiber 166 -GJ 1e10 {
set MatTag 121
source "../Scripts/Fiber_Custom-2_1.tcl"
set MatTag 163
source "../Scripts/Fiber_Custom-2_2.tcl"
}
section Fiber 165 -GJ 1e10 {
set MatTag 121
source "../Scripts/Fiber_Custom-1_1.tcl"
set MatTag 163
source "../Scripts/Fiber_Custom-1_2.tcl"
}
section Fiber 172 -GJ 1e10 {
set MatTag 121
source "../Scripts/Fiber_Custom-8_1.tcl"
set MatTag 163
source "../Scripts/Fiber_Custom-8_2.tcl"
}
section Fiber 171 -GJ 1e10 {
set MatTag 121
source "../Scripts/Fiber_Custom-7_1.tcl"
set MatTag 163
source "../Scripts/Fiber_Custom-7_2.tcl"
}
section Fiber 169 -GJ 1e10 {
set MatTag 121
source "../Scripts/Fiber_Custom-5_1.tcl"
set MatTag 163
source "../Scripts/Fiber_Custom-5_2.tcl"
}
section Fiber 176 -GJ 1e10 {
uniaxialMaterial Steel01 164 370000 2.1e+08 0.01 0 1 0 1
set MatTag 164
source "../Scripts/Fiber_Custom-12_1.tcl"
}
section Fiber 178 -GJ 1e10 {
set MatTag 164
source "../Scripts/Fiber_Custom-13_1.tcl"
}
section Fiber 177 -GJ 1e10 {
set MatTag 164
source "../Scripts/Fiber_Custom-12a_1.tcl"
}
section Fiber 179 -GJ 1e10 {
set MatTag 164
source "../Scripts/Fiber_Custom-14_1.tcl"
}
section Fiber 175 -GJ 1e10 {
set MatTag 164
source "../Scripts/Fiber_Custom-11_1.tcl"
}
section Fiber 191 -GJ 1e10 {
set MatTag 164
source "../Scripts/Fiber_Custom-20a_1.tcl"
}
section Fiber 194 -GJ 1e10 {
set MatTag 164
source "../Scripts/Fiber_Custom-21a_1.tcl"
}
section Fiber 192 -GJ 1e10 {
set MatTag 164
source "../Scripts/Fiber_Custom-20b_1.tcl"
}
section Fiber 193 -GJ 1e10 {
set MatTag 164
source "../Scripts/Fiber_Custom-21_1.tcl"
}
section Fiber 190 -GJ 1e10 {
set MatTag 164
source "../Scripts/Fiber_Custom-20_1.tcl"
}
section Fiber 195 -GJ 1e10 {
set MatTag 164
source "../Scripts/Fiber_Custom-22_1.tcl"
}
section Fiber 181 -GJ 1e10 {
set MatTag 164
source "../Scripts/Fiber_Custom-16_1.tcl"
}
section Fiber 189 -GJ 1e10 {
set MatTag 164
source "../Scripts/Fiber_Custom-19a_1.tcl"
}
section Fiber 182 -GJ 1e10 {
set MatTag 164
source "../Scripts/Fiber_Custom-16a_1.tcl"
}


element forceBeamColumn 1 282 270 3 178 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 2 270 271 3 177 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 3 271 263 3 176 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 4 263 264 3 177 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 5 264 254 3 176 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 6 254 256 3 177 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 7 256 248 3 176 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 8 248 252 3 179 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 9 252 242 3 176 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 10 242 244 3 176 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 11 244 236 3 175 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 12 236 240 3 176 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 13 240 230 3 175 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 14 230 232 3 176 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 15 232 223 3 176 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 16 223 228 3 179 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 17 228 218 3 176 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 18 218 220 3 175 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 19 220 213 3 175 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 20 213 216 3 175 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 21 216 208 3 175 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 22 208 210 3 175 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 23 210 203 3 176 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 807 72 71 3 185 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 808 75 74 3 187 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 809 85 84 3 187 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 810 82 81 3 186 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 811 92 91 3 185 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 812 95 94 3 187 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 813 132 131 3 185 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 814 122 121 3 187 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 815 112 111 3 185 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 816 102 101 3 187 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 817 135 134 3 189 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 818 125 124 3 188 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 819 115 114 3 187 2 -iter 10 1.00e-06 -mass 0
element forceBeamColumn 820 105 104 3 188 2 -iter 10 1.00e-06 -mass 0

# --------------------------------------------------------------------------------------------------------------
#
# D O M A I N C O M M O N S
#
# --------------------------------------------------------------------------------------------------------------

# --------------------------------------------------------------------------------------------------------------
# R E C O R D E R S
# --------------------------------------------------------------------------------------------------------------

recorder Node -file Node_displacements.out -time -nodeRange 1 286 -dof 1 2 3 disp
recorder Node -file Node_rotations.out -time -nodeRange 1 286 -dof 4 5 6 disp
recorder Node -file Node_forceReactions.out -time -nodeRange 1 286 -dof 1 2 3 reaction
recorder Node -file Node_momentReactions.out -time -nodeRange 1 286 -dof 4 5 6 reaction



puts "Running interval 1\n"

pattern Plain 100 Linear {
load 1 0 0 -1868066.2 0 0 0
load 1 0 0 -92022 0 0 0
load 1 0 0 -1868066.2 0 0 0
load 1 0 0 -127867.46 0 0 0
load 2 0 0 -127867.46 0 0 0
load 3 0 0 -92022 0 0 0
load 3 0 0 -83185.046 0 0 0
load 3 0 0 -87420.018 0 0 0
load 3 0 0 -92369.214 0 0 0
load 4 0 0 -24896.41 0 0 0
load 4 0 0 -196694.82 0 0 0
load 4 0 0 -1704151.4 0 0 0
load 4 0 0 -1868066.2 0 0 0
load 4 0 0 -24896.41 0 0 0
load 4 0 0 -1868066.2 0 0 0
load 4 0 0 -1704151.4 0 0 0
load 4 0 0 -92369.214 0 0 0
load 4 0 0 -127867.46 0 0 0
load 5 0 0 -1730562.4 0 0 0
load 5 0 0 -125874.14 0 0 0
load 5 0 0 -1730562.4 0 0 0
load 5 0 0 -87420.018 0 0 0
load 6 0 0 -50983.128 0 0 0
load 6 0 0 -50983.128 0 0 0
load 6 0 0 -25820.354 0 0 0
load 6 0 0 -49158.858 0 0 0
load 6 0 0 -25820.354 0 0 0
load 6 0 0 -24896.41 0 0 0
load 286 0 0 -92022 0 0 0
load 286 0 0 -1868066.2 0 0 0
load 286 0 0 -1868066.2 0 0 0
}

# recording the initial status

record

# Analysis options

system BandGeneral
numberer RCM
constraints Plain
integrator LoadControl 0.01
algorithm Linear
analysis Static
set committedSteps 1
for {set i 1} { $i <= 100 } {incr i 1} {
set AnalOk [analyze 1]
if {$AnalOk !=0} {
break
} else {
set committedSteps [expr $committedSteps+1]
}
}
if {$AnalOk == 0} {
puts "\nAnalysis completed SUCCESSFULLY"
puts "Committed steps : $committedSteps\n"
} else {
puts "\nAnalysis FAILED"
puts "Committed steps : $committedSteps\n"
}

# all previously defined patterns are constant for so on.

loadConst -time 0.0

# --------------------------------------------------------------------------------------------------------------

set hour 0.0
set minute 0.0
set second 0.0
set time_end [clock seconds]
set analysisTime [expr $time_end-$time_start]

puts "Analysis finished at [clock format $time_end -format %H:%M:%S]\n"

if {$analysisTime<60} {
if {$analysisTime==0} {
puts "Analysis time : less than one second"
} elseif {$analysisTime==1} {
puts "Analysis time : 1 second"
} else {
puts "Analysis time : $analysisTime seconds"
}

} elseif {$analysisTime<3600} {
set minutes [expr $analysisTime/60]
set seconds [expr $analysisTime%60]

if {$minutes==1} {
puts -nonewline "Analysis time : 1 minute"
} else {
puts -nonewline "Analysis time : $minutes minutes"
}

if {$seconds==0} {
puts ""
} elseif {$seconds==1} {
puts " and 1 second"
} else {
puts " and $seconds seconds"
}

} else {
set hours [expr $analysisTime/3600]
set minutes [expr ($analysisTime%3600)/60]
set seconds [expr ($analysisTime%3600)%60]

if {$hours==1} {
puts -nonewline "Analysis time : 1 hour"
} else {
puts -nonewline "Analysis time : $hours hours"
}

if {$minutes==0} {
} elseif {$minute==1} {
puts -nonewline ", 1 minute"
} else {
puts -nonewline ", $minutes minutes"
}

if {$seconds==0} {
puts ""
} elseif {$second==1} {
puts " and 1 second"
} else {
puts " and $seconds seconds"
}
}

puts "bridge gravity analysis completed"
wipe Analysis
source dynamic100wave01.tcl





constraints Plain
numberer RCM
system BandGeneral
test NormDispIncr 1.e-5 500
algorithm ModifiedNewton -initial
integrator Newmark 0.5 0.25
analysis Transient
puts "11111111111111"
#用瑞利阻尼求系数
set xDamp 0.05;
set MpropSwitch 1.0;
set KcurrSwitch 0.0;#当前迭代步刚度矩阵系数=0
set KcommSwitch 0.0;#上一个时步的迭代系数=0
set KinitSwitch 1.0;
puts "22222222222"
set nEigenI 1; # mode 1
set nEigenJ 2; # mode 2
puts "33333333333"
set lambdaN [eigen [expr $nEigenJ]]; # eigenvalue analysis for nEigenJ modes
set lambdaI [lindex $lambdaN [expr $nEigenI-1]]; # eigenvalue mode i
set lambdaJ [lindex $lambdaN [expr $nEigenJ-1]]; # eigenvalue mode j
puts "44444444444"
set omegaI [expr pow($lambdaI,0.5)];
set omegaJ [expr pow($lambdaJ,0.5)];
set alphaM [expr $MpropSwitch*$xDamp*(2*$omegaI*$omegaJ)/($omegaI+$omegaJ)]; # M-prop. damping; D = alphaM*M（质量矩阵系数）
set betaKcurr [expr $KcurrSwitch*2.*$xDamp/($omegaI+$omegaJ)]; # current-K; +beatKcurr*KCurrent，等于0
set betaKcomm [expr $KcommSwitch*2.*$xDamp/($omegaI+$omegaJ)]; # last-committed K; +betaKcomm*KlastCommitt，等于0
set betaKinit [expr $KinitSwitch*2.*$xDamp/($omegaI+$omegaJ)]; # initial-K; +beatKinit*Kini
rayleigh $alphaM $betaKcurr $betaKinit $betaKcomm; # RAYLEIGH damping，只有质量刚度矩阵系数，还有第三个刚度矩阵系数，其他都是0.
# puts "$alphaM $betaKcurr $betaKinit $betaKcomm"
# rayleigh 0.3054899518002805 0.0 0.002709476467775522 0.0

set T {}
foreach lam $lambdaN {
lappend T [expr 6.283185/sqrt($lam)];#开平方数squrt
puts "T: $T"
}


set GMdirection 1; # ground-motion direction，地震输入的方向，GM=ground-motion，即地震激励
set GMfile "la01" ; # ground-motion filenames
set GMfact [expr $PGA/461.255];
set dt [expr 0.02]; # time-step Dt for lateral analysis
set TmaxAnalysis [expr 53.48]; # 54 maximum duration of ground-motion analysis -- should be 50
# set g 9.8
set IDloadTag 400; # for uniformSupport excitation
set Wave GMfiles/$GMfile.txt; # set variable holding new filename (PEER files have .at2/dt2 extension)
set GMfatt [expr $g*$GMfact]; # data in input file is in g Unifts -- ACCELERATION TH（体系加速度放大系数为，峰值加速度与地面最大加速度之比）
set AccelSeries "Series -dt $dt -filePath $Wave -factor $GMfatt"; # time series information,=Pathseries
pattern UniformExcitation $IDloadTag $GMdirection -accel $AccelSeries ; # create Unifform excitation
set Nsteps [expr int($TmaxAnalysis/$dt)];
set ok [analyze $Nsteps $dt]; # actually perform analysis; returns ok=0 if analysis was successful

puts "bridge transient analysis completed"

[mhscott]

All of your masses are negative and only defined in the Z direction.

https://portwooddigital.com/2020/11/20/mass-and-weight/

Also, why do you define mass (and load) multiple times for each node?