Hi all. I have a problem with eigen command.
When I do eigenvalue analysis, why the eigenvalue change a little every time I run the program.
OpenSees > source dip.tcl
1st 1.329422e+003
3rd 1.329488e+003 2.944240e+004 3.479694e+004
4th 1.329237e+003 2.979972e+004 3.463235e+004 1.581617e+005
T1 0.172237698245
OpenSees > wipe
OpenSees > source dip.tcl
1st 1.328015e+003
3rd 1.327633e+003 3.024978e+004 3.490348e+004
4th 1.328162e+003 2.991004e+004 3.504907e+004 1.592121e+005
T1 0.172328914943
Thank you!
Eigenvalues are not constant
Moderators: silvia, selimgunay, Moderators
Hi fmk,
Thanks for your reply.
It is a model of bridge, whose cap beam and girder (I used section Elastic to simulate them) have a large stiffness. So I assign them with a large E(Young's Modulus) ,G(Shear Modulus) and J(torsional moment of inertia of section).
I found the problem. If the Yong's Modulus is assigned a very large number, such as 1.0e+20, the eigenvalue change a little every time. And when I reduced the number to 1.0e+10, the program works well. But the eigenvalue became larger.
So could you tell me if I want the cap beam and girder have a very large stiffness. How large the number is enough to assign to E, G, and J ?
Thank you in advance.
This is my script:
#################################################
#
#
# Available unit systems:
# 1) SI (kN, kPa, m, ton/m3)
# 2) English (inch, kip, kip-ft)
#
#################################################
# Unit system in this file: SI
# define UNITS ---------------------------------------------
set g 9.81; # gravitational acceleration
model basic -ndm 3 -ndf 6
#36_pile_left
node 1 0.00 0.00 0.00
node 2 0.00 0.00 -173.00
node 87 0.00 0.00 -43.25
node 88 0.00 0.00 -86.50
node 89 0.00 0.00 -129.75
node 3 46.49 99.69 0.00
node 4 46.49 99.69 -173.00
node 90 46.49 99.69 -43.25
node 91 46.49 99.69 -86.50
node 92 46.49 99.69 -129.75
node 5 92.98 199.39 0.00
node 6 92.98 199.39 -173.00
node 93 92.98 199.39 -43.25
node 94 92.98 199.39 -86.50
node 95 92.98 199.39 -129.75
node 7 139.46 299.08 0.00
node 8 139.46 299.08 -173.00
node 96 139.46 299.08 -43.25
node 97 139.46 299.08 -86.50
node 98 139.46 299.08 -129.75
node 9 185.95 398.78 0.00
node 10 185.95 398.78 -173.00
node 99 185.95 398.78 -43.25
node 100 185.95 398.78 -86.50
node 101 185.95 398.78 -129.75
node 11 232.44 498.47 0.00
node 12 232.44 498.47 -173.00
node 102 232.44 498.47 -43.25
node 103 232.44 498.47 -86.50
node 104 232.44 498.47 -129.75
node 13 278.93 598.16 0.00
node 14 278.93 598.16 -173.00
node 105 278.93 598.16 -43.25
node 106 278.93 598.16 -86.50
node 107 278.93 598.16 -129.75
node 15 325.42 697.86 0.00
node 16 325.42 697.86 -173.00
node 108 325.42 697.86 -43.25
node 109 325.42 697.86 -86.50
node 110 325.42 697.86 -129.75
node 17 371.90 797.55 0.00
node 18 371.90 797.55 -173.00
node 111 371.90 797.55 -43.25
node 112 371.90 797.55 -86.50
node 113 371.90 797.55 -129.75
node 19 418.39 897.24 0.00
node 20 418.39 897.24 -173.00
node 114 418.39 897.24 -43.25
node 115 418.39 897.24 -86.50
node 116 418.39 897.24 -129.75
#36_pile_right
node 21 1440.00 0.00 0.00
node 22 1440.00 0.00 -173.00
node 117 1440.00 0.00 -43.25
node 118 1440.00 0.00 -86.50
node 119 1440.00 0.00 -129.75
node 23 1486.49 99.69 0.00
node 24 1486.49 99.69 -173.00
node 120 1486.49 99.69 -43.25
node 121 1486.49 99.69 -86.50
node 122 1486.49 99.69 -129.75
node 25 1532.98 199.39 0.00
node 26 1532.98 199.39 -173.00
node 123 1532.98 199.39 -43.25
node 124 1532.98 199.39 -86.50
node 125 1532.98 199.39 -129.75
node 27 1579.46 299.08 0.00
node 28 1579.46 299.08 -173.00
node 126 1579.46 299.08 -43.25
node 127 1579.46 299.08 -86.50
node 128 1579.46 299.08 -129.75
node 29 1625.95 398.78 0.00
node 30 1625.95 398.78 -173.00
node 129 1625.95 398.78 -43.25
node 130 1625.95 398.78 -86.50
node 131 1625.95 398.78 -129.75
node 31 1672.44 498.47 0.00
node 32 1672.44 498.47 -173.00
node 132 1672.44 498.47 -43.25
node 133 1672.44 498.47 -86.50
node 134 1672.44 498.47 -129.75
node 33 1718.93 598.16 0.00
node 34 1718.93 598.16 -173.00
node 135 1718.93 598.16 -43.25
node 136 1718.93 598.16 -86.50
node 137 1718.93 598.16 -129.75
node 35 1765.42 697.86 0.00
node 36 1765.42 697.86 -173.00
node 138 1765.42 697.86 -43.25
node 139 1765.42 697.86 -86.50
node 140 1765.42 697.86 -129.75
node 37 1811.90 797.55 0.00
node 38 1811.90 797.55 -173.00
node 141 1811.90 797.55 -43.25
node 142 1811.90 797.55 -86.50
node 143 1811.90 797.55 -129.75
node 39 1858.39 897.24 0.00
node 40 1858.39 897.24 -173.00
node 144 1858.39 897.24 -43.25
node 145 1858.39 897.24 -86.50
node 146 1858.39 897.24 -129.75
#24_pile_left
node 41 -1408.00 0.00 0.00
node 42 -1408.00 0.00 -173.00
node 147 -1408.00 0.00 -43.25
node 148 -1408.00 0.00 -86.50
node 149 -1408.00 0.00 -129.75
node 43 -1366.16 89.72 0.00
node 44 -1366.16 89.72 -173.00
node 150 -1366.16 89.72 -43.25
node 151 -1366.16 89.72 -86.50
node 152 -1366.16 89.72 -129.75
node 45 -1324.32 179.45 0.00
node 46 -1324.32 179.45 -173.00
node 153 -1324.32 179.45 -43.25
node 154 -1324.32 179.45 -86.50
node 155 -1324.32 179.45 -129.75
node 47 -1282.48 269.17 0.00
node 48 -1282.48 269.17 -173.00
node 156 -1282.48 269.17 -43.25
node 157 -1282.48 269.17 -86.50
node 158 -1282.48 269.17 -129.75
node 49 -1240.64 358.90 0.00
node 50 -1240.64 358.90 -173.00
node 159 -1240.64 358.90 -43.25
node 160 -1240.64 358.90 -86.50
node 161 -1240.64 358.90 -129.75
node 51 -1198.80 448.62 0.00
node 52 -1198.80 448.62 -173.00
node 162 -1198.80 448.62 -43.25
node 163 -1198.80 448.62 -86.50
node 164 -1198.80 448.62 -129.75
node 53 -1156.96 538.35 0.00
node 54 -1156.96 538.35 -173.00
node 165 -1156.96 538.35 -43.25
node 166 -1156.96 538.35 -86.50
node 167 -1156.96 538.35 -129.75
node 55 -1115.13 628.07 0.00
node 56 -1115.13 628.07 -173.00
node 168 -1115.13 628.07 -43.25
node 169 -1115.13 628.07 -86.50
node 170 -1115.13 628.07 -129.75
node 57 -1073.29 717.80 0.00
node 58 -1073.29 717.80 -173.00
node 171 -1073.29 717.80 -43.25
node 172 -1073.29 717.80 -86.50
node 173 -1073.29 717.80 -129.75
node 59 -1031.45 807.52 0.00
node 60 -1031.45 807.52 -173.00
node 174 -1031.45 807.52 -43.25
node 175 -1031.45 807.52 -86.50
node 176 -1031.45 807.52 -129.75
node 61 -989.61 897.24 0.00
node 62 -989.61 897.24 -173.00
node 177 -989.61 897.24 -43.25
node 178 -989.61 897.24 -86.50
node 179 -989.61 897.24 -129.75
#24_pile_right
node 63 2848.00 0.00 0.00
node 64 2848.00 0.00 -173.00
node 180 2848.00 0.00 -43.25
node 181 2848.00 0.00 -86.50
node 182 2848.00 0.00 -129.75
node 65 2889.84 89.72 0.00
node 66 2889.84 89.72 -173.00
node 183 2889.84 89.72 -43.25
node 184 2889.84 89.72 -86.50
node 185 2889.84 89.72 -129.75
node 67 2931.68 179.45 0.00
node 68 2931.68 179.45 -173.00
node 186 2931.68 179.45 -43.25
node 187 2931.68 179.45 -86.50
node 188 2931.68 179.45 -129.75
node 69 2973.52 269.17 0.00
node 70 2973.52 269.17 -173.00
node 189 2973.52 269.17 -43.25
node 190 2973.52 269.17 -86.50
node 191 2973.52 269.17 -129.75
node 71 3015.36 358.90 0.00
node 72 3015.36 358.90 -173.00
node 192 3015.36 358.90 -43.25
node 193 3015.36 358.90 -86.50
node 194 3015.36 358.90 -129.75
node 73 3057.20 448.62 0.00
node 74 3057.20 448.62 -173.00
node 195 3057.20 448.62 -43.25
node 196 3057.20 448.62 -86.50
node 197 3057.20 448.62 -129.75
node 75 3099.04 538.35 0.00
node 76 3099.04 538.35 -173.00
node 198 3099.04 538.35 -43.25
node 199 3099.04 538.35 -86.50
node 200 3099.04 538.35 -129.75
node 77 3140.87 628.07 0.00
node 78 3140.87 628.07 -173.00
node 201 3140.87 628.07 -43.25
node 202 3140.87 628.07 -86.50
node 203 3140.87 628.07 -129.75
node 79 3182.71 717.80 0.00
node 80 3182.71 717.80 -173.00
node 204 3182.71 717.80 -43.25
node 205 3182.71 717.80 -86.50
node 206 3182.71 717.80 -129.75
node 81 3224.55 807.52 0.00
node 82 3224.55 807.52 -173.00
node 207 3224.55 807.52 -43.25
node 208 3224.55 807.52 -86.50
node 209 3224.55 807.52 -129.75
node 83 3266.39 897.24 0.00
node 84 3266.39 897.24 -173.00
node 210 3266.39 897.24 -43.25
node 211 3266.39 897.24 -86.50
node 212 3266.39 897.24 -129.75
#girder_main
node 85 -1408.00 398.78 0.00
node 86 3266.39 398.78 0.00
node 213 -1222.05 398.78 0.00
node 214 3033.95 398.78 0.00
#mass
mass 1 30.1 30.1 30.1 0 0 0
mass 3 30.1 30.1 30.1 0 0 0
mass 5 30.1 30.1 30.1 0 0 0
mass 7 30.1 30.1 30.1 0 0 0
mass 9 54.2 54.2 54.2 0 0 0
mass 11 30.1 30.1 30.1 0 0 0
mass 13 30.1 30.1 30.1 0 0 0
mass 15 30.1 30.1 30.1 0 0 0
mass 17 30.1 30.1 30.1 0 0 0
mass 19 30.1 30.1 30.1 0 0 0
mass 21 30.1 30.1 30.1 0 0 0
mass 23 30.1 30.1 30.1 0 0 0
mass 25 30.1 30.1 30.1 0 0 0
mass 27 30.1 30.1 30.1 0 0 0
mass 29 54.2 54.2 54.2 0 0 0
mass 31 30.1 30.1 30.1 0 0 0
mass 33 30.1 30.1 30.1 0 0 0
mass 35 30.1 30.1 30.1 0 0 0
mass 37 30.1 30.1 30.1 0 0 0
mass 39 30.1 30.1 30.1 0 0 0
mass 41 18.2 18.2 18.2 0 0 0
mass 43 18.2 18.2 18.2 0 0 0
mass 45 18.2 18.2 18.2 0 0 0
mass 47 18.2 18.2 18.2 0 0 0
mass 49 18.2 18.2 18.2 0 0 0
mass 51 18.2 18.2 18.2 0 0 0
mass 53 18.2 18.2 18.2 0 0 0
mass 55 18.2 18.2 18.2 0 0 0
mass 57 18.2 18.2 18.2 0 0 0
mass 59 18.2 18.2 18.2 0 0 0
mass 61 18.2 18.2 18.2 0 0 0
mass 63 18.2 18.2 18.2 0 0 0
mass 65 18.2 18.2 18.2 0 0 0
mass 67 18.2 18.2 18.2 0 0 0
mass 69 18.2 18.2 18.2 0 0 0
mass 71 18.2 18.2 18.2 0 0 0
mass 73 18.2 18.2 18.2 0 0 0
mass 75 18.2 18.2 18.2 0 0 0
mass 77 18.2 18.2 18.2 0 0 0
mass 79 18.2 18.2 18.2 0 0 0
mass 81 18.2 18.2 18.2 0 0 0
mass 83 18.2 18.2 18.2 0 0 0
mass 85 24.1 24.1 24.1 0 0 0
mass 86 24.1 24.1 24.1 0 0 0
mass 213 24.1 24.1 24.1 0 0 0
mass 214 24.1 24.1 24.1 0 0 0
#36_pile_left
fix 1 0 0 0 0 0 0
fix 2 1 1 1 1 1 1
fix 4 1 1 1 1 1 1
fix 6 1 1 1 1 1 1
fix 8 1 1 1 1 1 1
fix 10 1 1 1 1 1 1
fix 12 1 1 1 1 1 1
fix 14 1 1 1 1 1 1
fix 16 1 1 1 1 1 1
fix 18 1 1 1 1 1 1
fix 20 1 1 1 1 1 1
#36_pile_right
fix 22 1 1 1 1 1 1
fix 24 1 1 1 1 1 1
fix 26 1 1 1 1 1 1
fix 28 1 1 1 1 1 1
fix 30 1 1 1 1 1 1
fix 32 1 1 1 1 1 1
fix 34 1 1 1 1 1 1
fix 36 1 1 1 1 1 1
fix 38 1 1 1 1 1 1
fix 40 1 1 1 1 1 1
#24_pile_left
fix 42 1 1 1 1 1 1
fix 44 1 1 1 1 1 1
fix 46 1 1 1 1 1 1
fix 48 1 1 1 1 1 1
fix 50 1 1 1 1 1 1
fix 52 1 1 1 1 1 1
fix 54 1 1 1 1 1 1
fix 56 1 1 1 1 1 1
fix 58 1 1 1 1 1 1
fix 60 1 1 1 1 1 1
fix 62 1 1 1 1 1 1
#24_pile_right
fix 64 1 1 1 1 1 1
fix 66 1 1 1 1 1 1
fix 68 1 1 1 1 1 1
fix 70 1 1 1 1 1 1
fix 72 1 1 1 1 1 1
fix 74 1 1 1 1 1 1
fix 76 1 1 1 1 1 1
fix 78 1 1 1 1 1 1
fix 80 1 1 1 1 1 1
fix 82 1 1 1 1 1 1
fix 84 1 1 1 1 1 1
#girder_main
fix 85 1 1 1 1 1 1
fix 86 1 0 1 1 1 1
# Define materials for columns
uniaxialMaterial Concrete01 1001 -4.00 -0.002 0.0 -0.006
uniaxialMaterial Concrete01 1002 -5.00 -0.005 -3.5 -0.02
uniaxialMaterial Steel01 1003 60 30000 0.05
uniaxialMaterial Steel01 1004 60 30000 0.05
# 36'pile
section Fiber 10 {
patch circ 1001 16 4 0. 0. 0 17.25 0 360
patch circ 1003 16 1 0. 0. 17.25 18 0 360
layer circ 1004 10 0.724 0. 0. 14.41 0 360
}
# 24"pile
section Fiber 20 {
patch circ 1001 16 4 0. 0. 0 11.5 0 360
patch circ 1003 16 1 0. 0. 11.5 12 0 360
layer circ 1004 7 0.60 0. 0. 8.78 0 360
}
# cap beam
#section Elastic tag E A Iz Iy G J
section Elastic 3 1.0e+20 2520 7.56e+5 3.70e+5 1.0e+10 8.56e+5
section Elastic 4 1.0e+20 1680 5.04e+5 1.98e+5 1.0e+10 3.17e+5
#girder
section Elastic 5 1.0e+20 1215 5.05e+5 6.42e+5 1.0e+10 2.33e+5
uniaxialMaterial Elastic 2000 1e10; # define elastic torsional stiffness
# Attach torsion to the Fiber section
# tag uniTag uniCode secTag
section Aggregator 1 2000 T -section 10
section Aggregator 2 2000 T -section 20
geomTransf Linear 1001 0 -1 0; #pile
geomTransf Linear 1002 0 0 1; #main girder
geomTransf Linear 1003 0 0 1; #cap beam
set numIntgrPts 5
#36_pile_left leltag iN jN Nintpts Sec Trans
element dispBeamColumn 1 1 87 $numIntgrPts 1 1001
element dispBeamColumn 2 87 88 $numIntgrPts 1 1001
element dispBeamColumn 3 88 89 $numIntgrPts 1 1001
element dispBeamColumn 4 89 2 $numIntgrPts 1 1001
element dispBeamColumn 5 3 90 $numIntgrPts 1 1001
element dispBeamColumn 6 90 91 $numIntgrPts 1 1001
element dispBeamColumn 7 91 92 $numIntgrPts 1 1001
element dispBeamColumn 8 92 4 $numIntgrPts 1 1001
element dispBeamColumn 9 5 93 $numIntgrPts 1 1001
element dispBeamColumn 10 93 94 $numIntgrPts 1 1001
element dispBeamColumn 11 94 95 $numIntgrPts 1 1001
element dispBeamColumn 12 95 6 $numIntgrPts 1 1001
element dispBeamColumn 13 7 96 $numIntgrPts 1 1001
element dispBeamColumn 14 96 97 $numIntgrPts 1 1001
element dispBeamColumn 15 97 98 $numIntgrPts 1 1001
element dispBeamColumn 16 98 8 $numIntgrPts 1 1001
element dispBeamColumn 17 9 99 $numIntgrPts 1 1001
element dispBeamColumn 18 99 100 $numIntgrPts 1 1001
element dispBeamColumn 19 100 101 $numIntgrPts 1 1001
element dispBeamColumn 20 101 10 $numIntgrPts 1 1001
element dispBeamColumn 21 11 102 $numIntgrPts 1 1001
element dispBeamColumn 22 102 103 $numIntgrPts 1 1001
element dispBeamColumn 23 103 104 $numIntgrPts 1 1001
element dispBeamColumn 24 104 12 $numIntgrPts 1 1001
element dispBeamColumn 25 13 105 $numIntgrPts 1 1001
element dispBeamColumn 26 105 106 $numIntgrPts 1 1001
element dispBeamColumn 27 106 107 $numIntgrPts 1 1001
element dispBeamColumn 28 107 14 $numIntgrPts 1 1001
element dispBeamColumn 29 15 108 $numIntgrPts 1 1001
element dispBeamColumn 30 108 109 $numIntgrPts 1 1001
element dispBeamColumn 31 109 110 $numIntgrPts 1 1001
element dispBeamColumn 32 110 16 $numIntgrPts 1 1001
element dispBeamColumn 33 17 111 $numIntgrPts 1 1001
element dispBeamColumn 34 111 112 $numIntgrPts 1 1001
element dispBeamColumn 35 112 113 $numIntgrPts 1 1001
element dispBeamColumn 36 113 18 $numIntgrPts 1 1001
element dispBeamColumn 37 19 114 $numIntgrPts 1 1001
element dispBeamColumn 38 114 115 $numIntgrPts 1 1001
element dispBeamColumn 39 115 116 $numIntgrPts 1 1001
element dispBeamColumn 40 116 20 $numIntgrPts 1 1001
#36_pile_right
element dispBeamColumn 41 21 117 $numIntgrPts 1 1001
element dispBeamColumn 42 117 118 $numIntgrPts 1 1001
element dispBeamColumn 43 118 119 $numIntgrPts 1 1001
element dispBeamColumn 44 119 22 $numIntgrPts 1 1001
element dispBeamColumn 45 23 120 $numIntgrPts 1 1001
element dispBeamColumn 46 120 121 $numIntgrPts 1 1001
element dispBeamColumn 47 121 122 $numIntgrPts 1 1001
element dispBeamColumn 48 122 24 $numIntgrPts 1 1001
element dispBeamColumn 49 25 123 $numIntgrPts 1 1001
element dispBeamColumn 50 123 124 $numIntgrPts 1 1001
element dispBeamColumn 51 124 125 $numIntgrPts 1 1001
element dispBeamColumn 52 125 26 $numIntgrPts 1 1001
element dispBeamColumn 53 27 126 $numIntgrPts 1 1001
element dispBeamColumn 54 126 127 $numIntgrPts 1 1001
element dispBeamColumn 55 127 128 $numIntgrPts 1 1001
element dispBeamColumn 56 128 28 $numIntgrPts 1 1001
element dispBeamColumn 57 29 129 $numIntgrPts 1 1001
element dispBeamColumn 58 129 130 $numIntgrPts 1 1001
element dispBeamColumn 59 130 131 $numIntgrPts 1 1001
element dispBeamColumn 60 131 30 $numIntgrPts 1 1001
element dispBeamColumn 61 31 132 $numIntgrPts 1 1001
element dispBeamColumn 62 132 133 $numIntgrPts 1 1001
element dispBeamColumn 63 133 134 $numIntgrPts 1 1001
element dispBeamColumn 64 134 32 $numIntgrPts 1 1001
element dispBeamColumn 65 33 135 $numIntgrPts 1 1001
element dispBeamColumn 66 135 136 $numIntgrPts 1 1001
element dispBeamColumn 67 136 137 $numIntgrPts 1 1001
element dispBeamColumn 68 137 34 $numIntgrPts 1 1001
element dispBeamColumn 69 35 138 $numIntgrPts 1 1001
element dispBeamColumn 70 138 139 $numIntgrPts 1 1001
element dispBeamColumn 71 139 140 $numIntgrPts 1 1001
element dispBeamColumn 72 140 36 $numIntgrPts 1 1001
element dispBeamColumn 73 37 141 $numIntgrPts 1 1001
element dispBeamColumn 74 141 142 $numIntgrPts 1 1001
element dispBeamColumn 75 142 143 $numIntgrPts 1 1001
element dispBeamColumn 76 143 38 $numIntgrPts 1 1001
element dispBeamColumn 77 39 144 $numIntgrPts 1 1001
element dispBeamColumn 78 144 145 $numIntgrPts 1 1001
element dispBeamColumn 79 145 146 $numIntgrPts 1 1001
element dispBeamColumn 80 146 40 $numIntgrPts 1 1001
#24_pile_left
element dispBeamColumn 81 41 147 $numIntgrPts 2 1001
element dispBeamColumn 82 147 148 $numIntgrPts 2 1001
element dispBeamColumn 83 148 149 $numIntgrPts 2 1001
element dispBeamColumn 84 149 42 $numIntgrPts 2 1001
element dispBeamColumn 85 43 150 $numIntgrPts 2 1001
element dispBeamColumn 86 150 151 $numIntgrPts 2 1001
element dispBeamColumn 87 151 152 $numIntgrPts 2 1001
element dispBeamColumn 88 152 44 $numIntgrPts 2 1001
element dispBeamColumn 89 45 153 $numIntgrPts 2 1001
element dispBeamColumn 90 153 154 $numIntgrPts 2 1001
element dispBeamColumn 91 154 155 $numIntgrPts 2 1001
element dispBeamColumn 92 155 46 $numIntgrPts 2 1001
element dispBeamColumn 93 47 156 $numIntgrPts 2 1001
element dispBeamColumn 94 156 157 $numIntgrPts 2 1001
element dispBeamColumn 95 157 158 $numIntgrPts 2 1001
element dispBeamColumn 96 158 48 $numIntgrPts 2 1001
element dispBeamColumn 97 49 159 $numIntgrPts 2 1001
element dispBeamColumn 98 159 160 $numIntgrPts 2 1001
element dispBeamColumn 99 160 161 $numIntgrPts 2 1001
element dispBeamColumn 100 161 50 $numIntgrPts 2 1001
element dispBeamColumn 101 51 162 $numIntgrPts 2 1001
element dispBeamColumn 102 162 163 $numIntgrPts 2 1001
element dispBeamColumn 103 163 164 $numIntgrPts 2 1001
element dispBeamColumn 104 164 52 $numIntgrPts 2 1001
element dispBeamColumn 105 53 165 $numIntgrPts 2 1001
element dispBeamColumn 106 165 166 $numIntgrPts 2 1001
element dispBeamColumn 107 166 167 $numIntgrPts 2 1001
element dispBeamColumn 108 167 54 $numIntgrPts 2 1001
element dispBeamColumn 109 55 168 $numIntgrPts 2 1001
element dispBeamColumn 110 168 169 $numIntgrPts 2 1001
element dispBeamColumn 111 169 170 $numIntgrPts 2 1001
element dispBeamColumn 112 170 56 $numIntgrPts 2 1001
element dispBeamColumn 113 57 171 $numIntgrPts 2 1001
element dispBeamColumn 114 171 172 $numIntgrPts 2 1001
element dispBeamColumn 115 172 173 $numIntgrPts 2 1001
element dispBeamColumn 116 173 58 $numIntgrPts 2 1001
element dispBeamColumn 117 59 174 $numIntgrPts 2 1001
element dispBeamColumn 118 174 175 $numIntgrPts 2 1001
element dispBeamColumn 119 175 176 $numIntgrPts 2 1001
element dispBeamColumn 120 176 60 $numIntgrPts 2 1001
element dispBeamColumn 121 61 177 $numIntgrPts 2 1001
element dispBeamColumn 122 177 178 $numIntgrPts 2 1001
element dispBeamColumn 123 178 179 $numIntgrPts 2 1001
element dispBeamColumn 124 179 62 $numIntgrPts 2 1001
#24_pile_right
element dispBeamColumn 125 63 180 $numIntgrPts 2 1001
element dispBeamColumn 126 180 181 $numIntgrPts 2 1001
element dispBeamColumn 127 181 182 $numIntgrPts 2 1001
element dispBeamColumn 128 182 64 $numIntgrPts 2 1001
element dispBeamColumn 129 65 183 $numIntgrPts 2 1001
element dispBeamColumn 130 183 184 $numIntgrPts 2 1001
element dispBeamColumn 131 184 185 $numIntgrPts 2 1001
element dispBeamColumn 132 185 66 $numIntgrPts 2 1001
element dispBeamColumn 133 67 186 $numIntgrPts 2 1001
element dispBeamColumn 134 186 187 $numIntgrPts 2 1001
element dispBeamColumn 135 187 188 $numIntgrPts 2 1001
element dispBeamColumn 136 188 68 $numIntgrPts 2 1001
element dispBeamColumn 137 69 189 $numIntgrPts 2 1001
element dispBeamColumn 138 189 190 $numIntgrPts 2 1001
element dispBeamColumn 139 190 191 $numIntgrPts 2 1001
element dispBeamColumn 140 191 70 $numIntgrPts 2 1001
element dispBeamColumn 141 71 192 $numIntgrPts 2 1001
element dispBeamColumn 142 192 193 $numIntgrPts 2 1001
element dispBeamColumn 143 193 194 $numIntgrPts 2 1001
element dispBeamColumn 144 194 72 $numIntgrPts 2 1001
element dispBeamColumn 145 73 195 $numIntgrPts 2 1001
element dispBeamColumn 146 195 196 $numIntgrPts 2 1001
element dispBeamColumn 147 196 197 $numIntgrPts 2 1001
element dispBeamColumn 148 197 74 $numIntgrPts 2 1001
element dispBeamColumn 149 75 198 $numIntgrPts 2 1001
element dispBeamColumn 150 198 199 $numIntgrPts 2 1001
element dispBeamColumn 151 199 200 $numIntgrPts 2 1001
element dispBeamColumn 152 200 76 $numIntgrPts 2 1001
element dispBeamColumn 153 77 201 $numIntgrPts 2 1001
element dispBeamColumn 154 201 202 $numIntgrPts 2 1001
element dispBeamColumn 155 202 203 $numIntgrPts 2 1001
element dispBeamColumn 156 203 78 $numIntgrPts 2 1001
element dispBeamColumn 157 79 204 $numIntgrPts 2 1001
element dispBeamColumn 158 204 205 $numIntgrPts 2 1001
element dispBeamColumn 159 205 206 $numIntgrPts 2 1001
element dispBeamColumn 160 206 80 $numIntgrPts 2 1001
element dispBeamColumn 161 81 207 $numIntgrPts 2 1001
element dispBeamColumn 162 207 208 $numIntgrPts 2 1001
element dispBeamColumn 163 208 209 $numIntgrPts 2 1001
element dispBeamColumn 164 209 82 $numIntgrPts 2 1001
element dispBeamColumn 165 83 210 $numIntgrPts 2 1001
element dispBeamColumn 166 210 211 $numIntgrPts 2 1001
element dispBeamColumn 167 211 212 $numIntgrPts 2 1001
element dispBeamColumn 168 212 84 $numIntgrPts 2 1001
#girder_main
element dispBeamColumn 169 85 213 $numIntgrPts 5 1002
element dispBeamColumn 170 213 9 $numIntgrPts 5 1002
element dispBeamColumn 171 9 29 $numIntgrPts 5 1002
element dispBeamColumn 172 29 214 $numIntgrPts 5 1002
element dispBeamColumn 173 214 86 $numIntgrPts 5 1002
#girder_36L
element dispBeamColumn 174 1 3 $numIntgrPts 3 1003
element dispBeamColumn 175 3 5 $numIntgrPts 3 1003
element dispBeamColumn 176 5 7 $numIntgrPts 3 1003
element dispBeamColumn 177 7 9 $numIntgrPts 3 1003
element dispBeamColumn 178 9 11 $numIntgrPts 3 1003
element dispBeamColumn 179 11 13 $numIntgrPts 3 1003
element dispBeamColumn 180 13 15 $numIntgrPts 3 1003
element dispBeamColumn 181 15 17 $numIntgrPts 3 1003
element dispBeamColumn 182 17 19 $numIntgrPts 3 1003
#girder_36R
element dispBeamColumn 183 21 23 $numIntgrPts 3 1003
element dispBeamColumn 184 23 25 $numIntgrPts 3 1003
element dispBeamColumn 185 25 27 $numIntgrPts 3 1003
element dispBeamColumn 186 27 29 $numIntgrPts 3 1003
element dispBeamColumn 187 29 31 $numIntgrPts 3 1003
element dispBeamColumn 188 31 33 $numIntgrPts 3 1003
element dispBeamColumn 189 33 35 $numIntgrPts 3 1003
element dispBeamColumn 190 35 37 $numIntgrPts 3 1003
element dispBeamColumn 191 37 39 $numIntgrPts 3 1003
#girder_24L
element dispBeamColumn 192 41 43 $numIntgrPts 4 1003
element dispBeamColumn 193 43 45 $numIntgrPts 4 1003
element dispBeamColumn 194 45 47 $numIntgrPts 4 1003
element dispBeamColumn 195 47 49 $numIntgrPts 4 1003
element dispBeamColumn 196 49 213 $numIntgrPts 4 1003
element dispBeamColumn 197 213 51 $numIntgrPts 4 1003
element dispBeamColumn 198 51 53 $numIntgrPts 4 1003
element dispBeamColumn 199 53 55 $numIntgrPts 4 1003
element dispBeamColumn 200 55 57 $numIntgrPts 4 1003
element dispBeamColumn 201 57 59 $numIntgrPts 4 1003
element dispBeamColumn 202 59 61 $numIntgrPts 4 1003
#girder_24R
element dispBeamColumn 203 63 65 $numIntgrPts 4 1003
element dispBeamColumn 204 65 67 $numIntgrPts 4 1003
element dispBeamColumn 205 67 69 $numIntgrPts 4 1003
element dispBeamColumn 206 69 71 $numIntgrPts 4 1003
element dispBeamColumn 207 71 214 $numIntgrPts 4 1003
element dispBeamColumn 208 214 73 $numIntgrPts 4 1003
element dispBeamColumn 209 73 75 $numIntgrPts 4 1003
element dispBeamColumn 210 75 77 $numIntgrPts 4 1003
element dispBeamColumn 211 77 79 $numIntgrPts 4 1003
element dispBeamColumn 212 79 81 $numIntgrPts 4 1003
element dispBeamColumn 213 81 83 $numIntgrPts 4 1003
set first [eigen 1]
puts "1st $first"
set third [eigen 3]
puts "3rd $third"
set forth [eigen 4]
puts "4th $forth"
set omega [expr pow($first,0.5)]
set T1 [expr 2*3.14/$omega]
puts "T1 $T1"
Thanks for your reply.
It is a model of bridge, whose cap beam and girder (I used section Elastic to simulate them) have a large stiffness. So I assign them with a large E(Young's Modulus) ,G(Shear Modulus) and J(torsional moment of inertia of section).
I found the problem. If the Yong's Modulus is assigned a very large number, such as 1.0e+20, the eigenvalue change a little every time. And when I reduced the number to 1.0e+10, the program works well. But the eigenvalue became larger.
So could you tell me if I want the cap beam and girder have a very large stiffness. How large the number is enough to assign to E, G, and J ?
Thank you in advance.
This is my script:
#################################################
#
#
# Available unit systems:
# 1) SI (kN, kPa, m, ton/m3)
# 2) English (inch, kip, kip-ft)
#
#################################################
# Unit system in this file: SI
# define UNITS ---------------------------------------------
set g 9.81; # gravitational acceleration
model basic -ndm 3 -ndf 6
#36_pile_left
node 1 0.00 0.00 0.00
node 2 0.00 0.00 -173.00
node 87 0.00 0.00 -43.25
node 88 0.00 0.00 -86.50
node 89 0.00 0.00 -129.75
node 3 46.49 99.69 0.00
node 4 46.49 99.69 -173.00
node 90 46.49 99.69 -43.25
node 91 46.49 99.69 -86.50
node 92 46.49 99.69 -129.75
node 5 92.98 199.39 0.00
node 6 92.98 199.39 -173.00
node 93 92.98 199.39 -43.25
node 94 92.98 199.39 -86.50
node 95 92.98 199.39 -129.75
node 7 139.46 299.08 0.00
node 8 139.46 299.08 -173.00
node 96 139.46 299.08 -43.25
node 97 139.46 299.08 -86.50
node 98 139.46 299.08 -129.75
node 9 185.95 398.78 0.00
node 10 185.95 398.78 -173.00
node 99 185.95 398.78 -43.25
node 100 185.95 398.78 -86.50
node 101 185.95 398.78 -129.75
node 11 232.44 498.47 0.00
node 12 232.44 498.47 -173.00
node 102 232.44 498.47 -43.25
node 103 232.44 498.47 -86.50
node 104 232.44 498.47 -129.75
node 13 278.93 598.16 0.00
node 14 278.93 598.16 -173.00
node 105 278.93 598.16 -43.25
node 106 278.93 598.16 -86.50
node 107 278.93 598.16 -129.75
node 15 325.42 697.86 0.00
node 16 325.42 697.86 -173.00
node 108 325.42 697.86 -43.25
node 109 325.42 697.86 -86.50
node 110 325.42 697.86 -129.75
node 17 371.90 797.55 0.00
node 18 371.90 797.55 -173.00
node 111 371.90 797.55 -43.25
node 112 371.90 797.55 -86.50
node 113 371.90 797.55 -129.75
node 19 418.39 897.24 0.00
node 20 418.39 897.24 -173.00
node 114 418.39 897.24 -43.25
node 115 418.39 897.24 -86.50
node 116 418.39 897.24 -129.75
#36_pile_right
node 21 1440.00 0.00 0.00
node 22 1440.00 0.00 -173.00
node 117 1440.00 0.00 -43.25
node 118 1440.00 0.00 -86.50
node 119 1440.00 0.00 -129.75
node 23 1486.49 99.69 0.00
node 24 1486.49 99.69 -173.00
node 120 1486.49 99.69 -43.25
node 121 1486.49 99.69 -86.50
node 122 1486.49 99.69 -129.75
node 25 1532.98 199.39 0.00
node 26 1532.98 199.39 -173.00
node 123 1532.98 199.39 -43.25
node 124 1532.98 199.39 -86.50
node 125 1532.98 199.39 -129.75
node 27 1579.46 299.08 0.00
node 28 1579.46 299.08 -173.00
node 126 1579.46 299.08 -43.25
node 127 1579.46 299.08 -86.50
node 128 1579.46 299.08 -129.75
node 29 1625.95 398.78 0.00
node 30 1625.95 398.78 -173.00
node 129 1625.95 398.78 -43.25
node 130 1625.95 398.78 -86.50
node 131 1625.95 398.78 -129.75
node 31 1672.44 498.47 0.00
node 32 1672.44 498.47 -173.00
node 132 1672.44 498.47 -43.25
node 133 1672.44 498.47 -86.50
node 134 1672.44 498.47 -129.75
node 33 1718.93 598.16 0.00
node 34 1718.93 598.16 -173.00
node 135 1718.93 598.16 -43.25
node 136 1718.93 598.16 -86.50
node 137 1718.93 598.16 -129.75
node 35 1765.42 697.86 0.00
node 36 1765.42 697.86 -173.00
node 138 1765.42 697.86 -43.25
node 139 1765.42 697.86 -86.50
node 140 1765.42 697.86 -129.75
node 37 1811.90 797.55 0.00
node 38 1811.90 797.55 -173.00
node 141 1811.90 797.55 -43.25
node 142 1811.90 797.55 -86.50
node 143 1811.90 797.55 -129.75
node 39 1858.39 897.24 0.00
node 40 1858.39 897.24 -173.00
node 144 1858.39 897.24 -43.25
node 145 1858.39 897.24 -86.50
node 146 1858.39 897.24 -129.75
#24_pile_left
node 41 -1408.00 0.00 0.00
node 42 -1408.00 0.00 -173.00
node 147 -1408.00 0.00 -43.25
node 148 -1408.00 0.00 -86.50
node 149 -1408.00 0.00 -129.75
node 43 -1366.16 89.72 0.00
node 44 -1366.16 89.72 -173.00
node 150 -1366.16 89.72 -43.25
node 151 -1366.16 89.72 -86.50
node 152 -1366.16 89.72 -129.75
node 45 -1324.32 179.45 0.00
node 46 -1324.32 179.45 -173.00
node 153 -1324.32 179.45 -43.25
node 154 -1324.32 179.45 -86.50
node 155 -1324.32 179.45 -129.75
node 47 -1282.48 269.17 0.00
node 48 -1282.48 269.17 -173.00
node 156 -1282.48 269.17 -43.25
node 157 -1282.48 269.17 -86.50
node 158 -1282.48 269.17 -129.75
node 49 -1240.64 358.90 0.00
node 50 -1240.64 358.90 -173.00
node 159 -1240.64 358.90 -43.25
node 160 -1240.64 358.90 -86.50
node 161 -1240.64 358.90 -129.75
node 51 -1198.80 448.62 0.00
node 52 -1198.80 448.62 -173.00
node 162 -1198.80 448.62 -43.25
node 163 -1198.80 448.62 -86.50
node 164 -1198.80 448.62 -129.75
node 53 -1156.96 538.35 0.00
node 54 -1156.96 538.35 -173.00
node 165 -1156.96 538.35 -43.25
node 166 -1156.96 538.35 -86.50
node 167 -1156.96 538.35 -129.75
node 55 -1115.13 628.07 0.00
node 56 -1115.13 628.07 -173.00
node 168 -1115.13 628.07 -43.25
node 169 -1115.13 628.07 -86.50
node 170 -1115.13 628.07 -129.75
node 57 -1073.29 717.80 0.00
node 58 -1073.29 717.80 -173.00
node 171 -1073.29 717.80 -43.25
node 172 -1073.29 717.80 -86.50
node 173 -1073.29 717.80 -129.75
node 59 -1031.45 807.52 0.00
node 60 -1031.45 807.52 -173.00
node 174 -1031.45 807.52 -43.25
node 175 -1031.45 807.52 -86.50
node 176 -1031.45 807.52 -129.75
node 61 -989.61 897.24 0.00
node 62 -989.61 897.24 -173.00
node 177 -989.61 897.24 -43.25
node 178 -989.61 897.24 -86.50
node 179 -989.61 897.24 -129.75
#24_pile_right
node 63 2848.00 0.00 0.00
node 64 2848.00 0.00 -173.00
node 180 2848.00 0.00 -43.25
node 181 2848.00 0.00 -86.50
node 182 2848.00 0.00 -129.75
node 65 2889.84 89.72 0.00
node 66 2889.84 89.72 -173.00
node 183 2889.84 89.72 -43.25
node 184 2889.84 89.72 -86.50
node 185 2889.84 89.72 -129.75
node 67 2931.68 179.45 0.00
node 68 2931.68 179.45 -173.00
node 186 2931.68 179.45 -43.25
node 187 2931.68 179.45 -86.50
node 188 2931.68 179.45 -129.75
node 69 2973.52 269.17 0.00
node 70 2973.52 269.17 -173.00
node 189 2973.52 269.17 -43.25
node 190 2973.52 269.17 -86.50
node 191 2973.52 269.17 -129.75
node 71 3015.36 358.90 0.00
node 72 3015.36 358.90 -173.00
node 192 3015.36 358.90 -43.25
node 193 3015.36 358.90 -86.50
node 194 3015.36 358.90 -129.75
node 73 3057.20 448.62 0.00
node 74 3057.20 448.62 -173.00
node 195 3057.20 448.62 -43.25
node 196 3057.20 448.62 -86.50
node 197 3057.20 448.62 -129.75
node 75 3099.04 538.35 0.00
node 76 3099.04 538.35 -173.00
node 198 3099.04 538.35 -43.25
node 199 3099.04 538.35 -86.50
node 200 3099.04 538.35 -129.75
node 77 3140.87 628.07 0.00
node 78 3140.87 628.07 -173.00
node 201 3140.87 628.07 -43.25
node 202 3140.87 628.07 -86.50
node 203 3140.87 628.07 -129.75
node 79 3182.71 717.80 0.00
node 80 3182.71 717.80 -173.00
node 204 3182.71 717.80 -43.25
node 205 3182.71 717.80 -86.50
node 206 3182.71 717.80 -129.75
node 81 3224.55 807.52 0.00
node 82 3224.55 807.52 -173.00
node 207 3224.55 807.52 -43.25
node 208 3224.55 807.52 -86.50
node 209 3224.55 807.52 -129.75
node 83 3266.39 897.24 0.00
node 84 3266.39 897.24 -173.00
node 210 3266.39 897.24 -43.25
node 211 3266.39 897.24 -86.50
node 212 3266.39 897.24 -129.75
#girder_main
node 85 -1408.00 398.78 0.00
node 86 3266.39 398.78 0.00
node 213 -1222.05 398.78 0.00
node 214 3033.95 398.78 0.00
#mass
mass 1 30.1 30.1 30.1 0 0 0
mass 3 30.1 30.1 30.1 0 0 0
mass 5 30.1 30.1 30.1 0 0 0
mass 7 30.1 30.1 30.1 0 0 0
mass 9 54.2 54.2 54.2 0 0 0
mass 11 30.1 30.1 30.1 0 0 0
mass 13 30.1 30.1 30.1 0 0 0
mass 15 30.1 30.1 30.1 0 0 0
mass 17 30.1 30.1 30.1 0 0 0
mass 19 30.1 30.1 30.1 0 0 0
mass 21 30.1 30.1 30.1 0 0 0
mass 23 30.1 30.1 30.1 0 0 0
mass 25 30.1 30.1 30.1 0 0 0
mass 27 30.1 30.1 30.1 0 0 0
mass 29 54.2 54.2 54.2 0 0 0
mass 31 30.1 30.1 30.1 0 0 0
mass 33 30.1 30.1 30.1 0 0 0
mass 35 30.1 30.1 30.1 0 0 0
mass 37 30.1 30.1 30.1 0 0 0
mass 39 30.1 30.1 30.1 0 0 0
mass 41 18.2 18.2 18.2 0 0 0
mass 43 18.2 18.2 18.2 0 0 0
mass 45 18.2 18.2 18.2 0 0 0
mass 47 18.2 18.2 18.2 0 0 0
mass 49 18.2 18.2 18.2 0 0 0
mass 51 18.2 18.2 18.2 0 0 0
mass 53 18.2 18.2 18.2 0 0 0
mass 55 18.2 18.2 18.2 0 0 0
mass 57 18.2 18.2 18.2 0 0 0
mass 59 18.2 18.2 18.2 0 0 0
mass 61 18.2 18.2 18.2 0 0 0
mass 63 18.2 18.2 18.2 0 0 0
mass 65 18.2 18.2 18.2 0 0 0
mass 67 18.2 18.2 18.2 0 0 0
mass 69 18.2 18.2 18.2 0 0 0
mass 71 18.2 18.2 18.2 0 0 0
mass 73 18.2 18.2 18.2 0 0 0
mass 75 18.2 18.2 18.2 0 0 0
mass 77 18.2 18.2 18.2 0 0 0
mass 79 18.2 18.2 18.2 0 0 0
mass 81 18.2 18.2 18.2 0 0 0
mass 83 18.2 18.2 18.2 0 0 0
mass 85 24.1 24.1 24.1 0 0 0
mass 86 24.1 24.1 24.1 0 0 0
mass 213 24.1 24.1 24.1 0 0 0
mass 214 24.1 24.1 24.1 0 0 0
#36_pile_left
fix 1 0 0 0 0 0 0
fix 2 1 1 1 1 1 1
fix 4 1 1 1 1 1 1
fix 6 1 1 1 1 1 1
fix 8 1 1 1 1 1 1
fix 10 1 1 1 1 1 1
fix 12 1 1 1 1 1 1
fix 14 1 1 1 1 1 1
fix 16 1 1 1 1 1 1
fix 18 1 1 1 1 1 1
fix 20 1 1 1 1 1 1
#36_pile_right
fix 22 1 1 1 1 1 1
fix 24 1 1 1 1 1 1
fix 26 1 1 1 1 1 1
fix 28 1 1 1 1 1 1
fix 30 1 1 1 1 1 1
fix 32 1 1 1 1 1 1
fix 34 1 1 1 1 1 1
fix 36 1 1 1 1 1 1
fix 38 1 1 1 1 1 1
fix 40 1 1 1 1 1 1
#24_pile_left
fix 42 1 1 1 1 1 1
fix 44 1 1 1 1 1 1
fix 46 1 1 1 1 1 1
fix 48 1 1 1 1 1 1
fix 50 1 1 1 1 1 1
fix 52 1 1 1 1 1 1
fix 54 1 1 1 1 1 1
fix 56 1 1 1 1 1 1
fix 58 1 1 1 1 1 1
fix 60 1 1 1 1 1 1
fix 62 1 1 1 1 1 1
#24_pile_right
fix 64 1 1 1 1 1 1
fix 66 1 1 1 1 1 1
fix 68 1 1 1 1 1 1
fix 70 1 1 1 1 1 1
fix 72 1 1 1 1 1 1
fix 74 1 1 1 1 1 1
fix 76 1 1 1 1 1 1
fix 78 1 1 1 1 1 1
fix 80 1 1 1 1 1 1
fix 82 1 1 1 1 1 1
fix 84 1 1 1 1 1 1
#girder_main
fix 85 1 1 1 1 1 1
fix 86 1 0 1 1 1 1
# Define materials for columns
uniaxialMaterial Concrete01 1001 -4.00 -0.002 0.0 -0.006
uniaxialMaterial Concrete01 1002 -5.00 -0.005 -3.5 -0.02
uniaxialMaterial Steel01 1003 60 30000 0.05
uniaxialMaterial Steel01 1004 60 30000 0.05
# 36'pile
section Fiber 10 {
patch circ 1001 16 4 0. 0. 0 17.25 0 360
patch circ 1003 16 1 0. 0. 17.25 18 0 360
layer circ 1004 10 0.724 0. 0. 14.41 0 360
}
# 24"pile
section Fiber 20 {
patch circ 1001 16 4 0. 0. 0 11.5 0 360
patch circ 1003 16 1 0. 0. 11.5 12 0 360
layer circ 1004 7 0.60 0. 0. 8.78 0 360
}
# cap beam
#section Elastic tag E A Iz Iy G J
section Elastic 3 1.0e+20 2520 7.56e+5 3.70e+5 1.0e+10 8.56e+5
section Elastic 4 1.0e+20 1680 5.04e+5 1.98e+5 1.0e+10 3.17e+5
#girder
section Elastic 5 1.0e+20 1215 5.05e+5 6.42e+5 1.0e+10 2.33e+5
uniaxialMaterial Elastic 2000 1e10; # define elastic torsional stiffness
# Attach torsion to the Fiber section
# tag uniTag uniCode secTag
section Aggregator 1 2000 T -section 10
section Aggregator 2 2000 T -section 20
geomTransf Linear 1001 0 -1 0; #pile
geomTransf Linear 1002 0 0 1; #main girder
geomTransf Linear 1003 0 0 1; #cap beam
set numIntgrPts 5
#36_pile_left leltag iN jN Nintpts Sec Trans
element dispBeamColumn 1 1 87 $numIntgrPts 1 1001
element dispBeamColumn 2 87 88 $numIntgrPts 1 1001
element dispBeamColumn 3 88 89 $numIntgrPts 1 1001
element dispBeamColumn 4 89 2 $numIntgrPts 1 1001
element dispBeamColumn 5 3 90 $numIntgrPts 1 1001
element dispBeamColumn 6 90 91 $numIntgrPts 1 1001
element dispBeamColumn 7 91 92 $numIntgrPts 1 1001
element dispBeamColumn 8 92 4 $numIntgrPts 1 1001
element dispBeamColumn 9 5 93 $numIntgrPts 1 1001
element dispBeamColumn 10 93 94 $numIntgrPts 1 1001
element dispBeamColumn 11 94 95 $numIntgrPts 1 1001
element dispBeamColumn 12 95 6 $numIntgrPts 1 1001
element dispBeamColumn 13 7 96 $numIntgrPts 1 1001
element dispBeamColumn 14 96 97 $numIntgrPts 1 1001
element dispBeamColumn 15 97 98 $numIntgrPts 1 1001
element dispBeamColumn 16 98 8 $numIntgrPts 1 1001
element dispBeamColumn 17 9 99 $numIntgrPts 1 1001
element dispBeamColumn 18 99 100 $numIntgrPts 1 1001
element dispBeamColumn 19 100 101 $numIntgrPts 1 1001
element dispBeamColumn 20 101 10 $numIntgrPts 1 1001
element dispBeamColumn 21 11 102 $numIntgrPts 1 1001
element dispBeamColumn 22 102 103 $numIntgrPts 1 1001
element dispBeamColumn 23 103 104 $numIntgrPts 1 1001
element dispBeamColumn 24 104 12 $numIntgrPts 1 1001
element dispBeamColumn 25 13 105 $numIntgrPts 1 1001
element dispBeamColumn 26 105 106 $numIntgrPts 1 1001
element dispBeamColumn 27 106 107 $numIntgrPts 1 1001
element dispBeamColumn 28 107 14 $numIntgrPts 1 1001
element dispBeamColumn 29 15 108 $numIntgrPts 1 1001
element dispBeamColumn 30 108 109 $numIntgrPts 1 1001
element dispBeamColumn 31 109 110 $numIntgrPts 1 1001
element dispBeamColumn 32 110 16 $numIntgrPts 1 1001
element dispBeamColumn 33 17 111 $numIntgrPts 1 1001
element dispBeamColumn 34 111 112 $numIntgrPts 1 1001
element dispBeamColumn 35 112 113 $numIntgrPts 1 1001
element dispBeamColumn 36 113 18 $numIntgrPts 1 1001
element dispBeamColumn 37 19 114 $numIntgrPts 1 1001
element dispBeamColumn 38 114 115 $numIntgrPts 1 1001
element dispBeamColumn 39 115 116 $numIntgrPts 1 1001
element dispBeamColumn 40 116 20 $numIntgrPts 1 1001
#36_pile_right
element dispBeamColumn 41 21 117 $numIntgrPts 1 1001
element dispBeamColumn 42 117 118 $numIntgrPts 1 1001
element dispBeamColumn 43 118 119 $numIntgrPts 1 1001
element dispBeamColumn 44 119 22 $numIntgrPts 1 1001
element dispBeamColumn 45 23 120 $numIntgrPts 1 1001
element dispBeamColumn 46 120 121 $numIntgrPts 1 1001
element dispBeamColumn 47 121 122 $numIntgrPts 1 1001
element dispBeamColumn 48 122 24 $numIntgrPts 1 1001
element dispBeamColumn 49 25 123 $numIntgrPts 1 1001
element dispBeamColumn 50 123 124 $numIntgrPts 1 1001
element dispBeamColumn 51 124 125 $numIntgrPts 1 1001
element dispBeamColumn 52 125 26 $numIntgrPts 1 1001
element dispBeamColumn 53 27 126 $numIntgrPts 1 1001
element dispBeamColumn 54 126 127 $numIntgrPts 1 1001
element dispBeamColumn 55 127 128 $numIntgrPts 1 1001
element dispBeamColumn 56 128 28 $numIntgrPts 1 1001
element dispBeamColumn 57 29 129 $numIntgrPts 1 1001
element dispBeamColumn 58 129 130 $numIntgrPts 1 1001
element dispBeamColumn 59 130 131 $numIntgrPts 1 1001
element dispBeamColumn 60 131 30 $numIntgrPts 1 1001
element dispBeamColumn 61 31 132 $numIntgrPts 1 1001
element dispBeamColumn 62 132 133 $numIntgrPts 1 1001
element dispBeamColumn 63 133 134 $numIntgrPts 1 1001
element dispBeamColumn 64 134 32 $numIntgrPts 1 1001
element dispBeamColumn 65 33 135 $numIntgrPts 1 1001
element dispBeamColumn 66 135 136 $numIntgrPts 1 1001
element dispBeamColumn 67 136 137 $numIntgrPts 1 1001
element dispBeamColumn 68 137 34 $numIntgrPts 1 1001
element dispBeamColumn 69 35 138 $numIntgrPts 1 1001
element dispBeamColumn 70 138 139 $numIntgrPts 1 1001
element dispBeamColumn 71 139 140 $numIntgrPts 1 1001
element dispBeamColumn 72 140 36 $numIntgrPts 1 1001
element dispBeamColumn 73 37 141 $numIntgrPts 1 1001
element dispBeamColumn 74 141 142 $numIntgrPts 1 1001
element dispBeamColumn 75 142 143 $numIntgrPts 1 1001
element dispBeamColumn 76 143 38 $numIntgrPts 1 1001
element dispBeamColumn 77 39 144 $numIntgrPts 1 1001
element dispBeamColumn 78 144 145 $numIntgrPts 1 1001
element dispBeamColumn 79 145 146 $numIntgrPts 1 1001
element dispBeamColumn 80 146 40 $numIntgrPts 1 1001
#24_pile_left
element dispBeamColumn 81 41 147 $numIntgrPts 2 1001
element dispBeamColumn 82 147 148 $numIntgrPts 2 1001
element dispBeamColumn 83 148 149 $numIntgrPts 2 1001
element dispBeamColumn 84 149 42 $numIntgrPts 2 1001
element dispBeamColumn 85 43 150 $numIntgrPts 2 1001
element dispBeamColumn 86 150 151 $numIntgrPts 2 1001
element dispBeamColumn 87 151 152 $numIntgrPts 2 1001
element dispBeamColumn 88 152 44 $numIntgrPts 2 1001
element dispBeamColumn 89 45 153 $numIntgrPts 2 1001
element dispBeamColumn 90 153 154 $numIntgrPts 2 1001
element dispBeamColumn 91 154 155 $numIntgrPts 2 1001
element dispBeamColumn 92 155 46 $numIntgrPts 2 1001
element dispBeamColumn 93 47 156 $numIntgrPts 2 1001
element dispBeamColumn 94 156 157 $numIntgrPts 2 1001
element dispBeamColumn 95 157 158 $numIntgrPts 2 1001
element dispBeamColumn 96 158 48 $numIntgrPts 2 1001
element dispBeamColumn 97 49 159 $numIntgrPts 2 1001
element dispBeamColumn 98 159 160 $numIntgrPts 2 1001
element dispBeamColumn 99 160 161 $numIntgrPts 2 1001
element dispBeamColumn 100 161 50 $numIntgrPts 2 1001
element dispBeamColumn 101 51 162 $numIntgrPts 2 1001
element dispBeamColumn 102 162 163 $numIntgrPts 2 1001
element dispBeamColumn 103 163 164 $numIntgrPts 2 1001
element dispBeamColumn 104 164 52 $numIntgrPts 2 1001
element dispBeamColumn 105 53 165 $numIntgrPts 2 1001
element dispBeamColumn 106 165 166 $numIntgrPts 2 1001
element dispBeamColumn 107 166 167 $numIntgrPts 2 1001
element dispBeamColumn 108 167 54 $numIntgrPts 2 1001
element dispBeamColumn 109 55 168 $numIntgrPts 2 1001
element dispBeamColumn 110 168 169 $numIntgrPts 2 1001
element dispBeamColumn 111 169 170 $numIntgrPts 2 1001
element dispBeamColumn 112 170 56 $numIntgrPts 2 1001
element dispBeamColumn 113 57 171 $numIntgrPts 2 1001
element dispBeamColumn 114 171 172 $numIntgrPts 2 1001
element dispBeamColumn 115 172 173 $numIntgrPts 2 1001
element dispBeamColumn 116 173 58 $numIntgrPts 2 1001
element dispBeamColumn 117 59 174 $numIntgrPts 2 1001
element dispBeamColumn 118 174 175 $numIntgrPts 2 1001
element dispBeamColumn 119 175 176 $numIntgrPts 2 1001
element dispBeamColumn 120 176 60 $numIntgrPts 2 1001
element dispBeamColumn 121 61 177 $numIntgrPts 2 1001
element dispBeamColumn 122 177 178 $numIntgrPts 2 1001
element dispBeamColumn 123 178 179 $numIntgrPts 2 1001
element dispBeamColumn 124 179 62 $numIntgrPts 2 1001
#24_pile_right
element dispBeamColumn 125 63 180 $numIntgrPts 2 1001
element dispBeamColumn 126 180 181 $numIntgrPts 2 1001
element dispBeamColumn 127 181 182 $numIntgrPts 2 1001
element dispBeamColumn 128 182 64 $numIntgrPts 2 1001
element dispBeamColumn 129 65 183 $numIntgrPts 2 1001
element dispBeamColumn 130 183 184 $numIntgrPts 2 1001
element dispBeamColumn 131 184 185 $numIntgrPts 2 1001
element dispBeamColumn 132 185 66 $numIntgrPts 2 1001
element dispBeamColumn 133 67 186 $numIntgrPts 2 1001
element dispBeamColumn 134 186 187 $numIntgrPts 2 1001
element dispBeamColumn 135 187 188 $numIntgrPts 2 1001
element dispBeamColumn 136 188 68 $numIntgrPts 2 1001
element dispBeamColumn 137 69 189 $numIntgrPts 2 1001
element dispBeamColumn 138 189 190 $numIntgrPts 2 1001
element dispBeamColumn 139 190 191 $numIntgrPts 2 1001
element dispBeamColumn 140 191 70 $numIntgrPts 2 1001
element dispBeamColumn 141 71 192 $numIntgrPts 2 1001
element dispBeamColumn 142 192 193 $numIntgrPts 2 1001
element dispBeamColumn 143 193 194 $numIntgrPts 2 1001
element dispBeamColumn 144 194 72 $numIntgrPts 2 1001
element dispBeamColumn 145 73 195 $numIntgrPts 2 1001
element dispBeamColumn 146 195 196 $numIntgrPts 2 1001
element dispBeamColumn 147 196 197 $numIntgrPts 2 1001
element dispBeamColumn 148 197 74 $numIntgrPts 2 1001
element dispBeamColumn 149 75 198 $numIntgrPts 2 1001
element dispBeamColumn 150 198 199 $numIntgrPts 2 1001
element dispBeamColumn 151 199 200 $numIntgrPts 2 1001
element dispBeamColumn 152 200 76 $numIntgrPts 2 1001
element dispBeamColumn 153 77 201 $numIntgrPts 2 1001
element dispBeamColumn 154 201 202 $numIntgrPts 2 1001
element dispBeamColumn 155 202 203 $numIntgrPts 2 1001
element dispBeamColumn 156 203 78 $numIntgrPts 2 1001
element dispBeamColumn 157 79 204 $numIntgrPts 2 1001
element dispBeamColumn 158 204 205 $numIntgrPts 2 1001
element dispBeamColumn 159 205 206 $numIntgrPts 2 1001
element dispBeamColumn 160 206 80 $numIntgrPts 2 1001
element dispBeamColumn 161 81 207 $numIntgrPts 2 1001
element dispBeamColumn 162 207 208 $numIntgrPts 2 1001
element dispBeamColumn 163 208 209 $numIntgrPts 2 1001
element dispBeamColumn 164 209 82 $numIntgrPts 2 1001
element dispBeamColumn 165 83 210 $numIntgrPts 2 1001
element dispBeamColumn 166 210 211 $numIntgrPts 2 1001
element dispBeamColumn 167 211 212 $numIntgrPts 2 1001
element dispBeamColumn 168 212 84 $numIntgrPts 2 1001
#girder_main
element dispBeamColumn 169 85 213 $numIntgrPts 5 1002
element dispBeamColumn 170 213 9 $numIntgrPts 5 1002
element dispBeamColumn 171 9 29 $numIntgrPts 5 1002
element dispBeamColumn 172 29 214 $numIntgrPts 5 1002
element dispBeamColumn 173 214 86 $numIntgrPts 5 1002
#girder_36L
element dispBeamColumn 174 1 3 $numIntgrPts 3 1003
element dispBeamColumn 175 3 5 $numIntgrPts 3 1003
element dispBeamColumn 176 5 7 $numIntgrPts 3 1003
element dispBeamColumn 177 7 9 $numIntgrPts 3 1003
element dispBeamColumn 178 9 11 $numIntgrPts 3 1003
element dispBeamColumn 179 11 13 $numIntgrPts 3 1003
element dispBeamColumn 180 13 15 $numIntgrPts 3 1003
element dispBeamColumn 181 15 17 $numIntgrPts 3 1003
element dispBeamColumn 182 17 19 $numIntgrPts 3 1003
#girder_36R
element dispBeamColumn 183 21 23 $numIntgrPts 3 1003
element dispBeamColumn 184 23 25 $numIntgrPts 3 1003
element dispBeamColumn 185 25 27 $numIntgrPts 3 1003
element dispBeamColumn 186 27 29 $numIntgrPts 3 1003
element dispBeamColumn 187 29 31 $numIntgrPts 3 1003
element dispBeamColumn 188 31 33 $numIntgrPts 3 1003
element dispBeamColumn 189 33 35 $numIntgrPts 3 1003
element dispBeamColumn 190 35 37 $numIntgrPts 3 1003
element dispBeamColumn 191 37 39 $numIntgrPts 3 1003
#girder_24L
element dispBeamColumn 192 41 43 $numIntgrPts 4 1003
element dispBeamColumn 193 43 45 $numIntgrPts 4 1003
element dispBeamColumn 194 45 47 $numIntgrPts 4 1003
element dispBeamColumn 195 47 49 $numIntgrPts 4 1003
element dispBeamColumn 196 49 213 $numIntgrPts 4 1003
element dispBeamColumn 197 213 51 $numIntgrPts 4 1003
element dispBeamColumn 198 51 53 $numIntgrPts 4 1003
element dispBeamColumn 199 53 55 $numIntgrPts 4 1003
element dispBeamColumn 200 55 57 $numIntgrPts 4 1003
element dispBeamColumn 201 57 59 $numIntgrPts 4 1003
element dispBeamColumn 202 59 61 $numIntgrPts 4 1003
#girder_24R
element dispBeamColumn 203 63 65 $numIntgrPts 4 1003
element dispBeamColumn 204 65 67 $numIntgrPts 4 1003
element dispBeamColumn 205 67 69 $numIntgrPts 4 1003
element dispBeamColumn 206 69 71 $numIntgrPts 4 1003
element dispBeamColumn 207 71 214 $numIntgrPts 4 1003
element dispBeamColumn 208 214 73 $numIntgrPts 4 1003
element dispBeamColumn 209 73 75 $numIntgrPts 4 1003
element dispBeamColumn 210 75 77 $numIntgrPts 4 1003
element dispBeamColumn 211 77 79 $numIntgrPts 4 1003
element dispBeamColumn 212 79 81 $numIntgrPts 4 1003
element dispBeamColumn 213 81 83 $numIntgrPts 4 1003
set first [eigen 1]
puts "1st $first"
set third [eigen 3]
puts "3rd $third"
set forth [eigen 4]
puts "4th $forth"
set omega [expr pow($first,0.5)]
set T1 [expr 2*3.14/$omega]
puts "T1 $T1"
A good man!
Thanks fmk,
As you know it is a simple model. I have checked my script. I cannot find the problem and figure out why this happened.
I think it is must the same problem which caused there is not deformation in the ModeShape with DisplayPlane.tcl and DisplayModel3D.tcl.
Could you please give me some advice?
Looking forward to your reply!
As you know it is a simple model. I have checked my script. I cannot find the problem and figure out why this happened.
I think it is must the same problem which caused there is not deformation in the ModeShape with DisplayPlane.tcl and DisplayModel3D.tcl.
Could you please give me some advice?
Looking forward to your reply!