第五章 結論與建議
5.2 建議
以下為本研究採用 OpenSees 之集中塑性模型(Concentrated Plasticity Model)執行側推分析，並提供參考與改良之建議：
1. 在 OpenSees 分析軟體上，使用者必須具備有一定的基本程式撰寫能 力，且必須清楚其指令的使用方式，而使用者在程式的模型建構上 自由度相當高，可以藉由各種概念建構模型並進行分析，最後在進 行其模型分析結果驗證，確認模型之可行性。
2. 在本研究之模型中皆未考慮結構本身自重，後續研究者可嘗試加入 自重與活載重等。本研究除了懸臂梁為3D 模型，其餘範例皆為2D，
後續研究者可嘗試創建3D 模型和加入樓板元件，並考慮剛性樓板之 設定。
3. 在本研究中有幾個範例，在分析結果中，容量曲線並未產生下降端，
後續研究者可以從 Analysis Object 物件中的 ConvergenceTest type object 與 SolutionAlgorithm object 嘗試改變其收斂方式和迭代算法。
參考文獻
1. 鐘立來、葉勇凱、簡文郁 等，「校舍結構耐震評估與補強技術手冊 第二版」，國家地震中心，NCREE-09-023，2009。
2. ACT40, Seismic evaluation and retrofit of concrete buildings. Report No. SSC 96-1, Applied Technology Council, 1996.
3. FEMA 273, NEHRP Guidelines for the seismic rehabilitation of buildings, Federal Emergency Management Agency, Washington, D.C., 1997.
4. Freeman, Nicoletti, J.P. , “Evaluations to existing buildings for seismic risk-A case study of Puget Sound Naval Shipyard”, Washington, Proceeding of 1st U.S. National Conference on Earthquake Engineering, EERI, Berkeley, 113-122, 1975.
5. Saiidi, M., M. A. Sozen., “Simple nonlinear seismic analysis of R/C structures”, Journal of Structural Division , 107(8), pp. 937-953, 1981.
6. Miranda, E., “Seismic evaluation and upgrading of existing buildings”, Ph. D. Dissertation, University of California at Berkeley, California, 1991.
7. 張凡，「高層建築結構的 Pushover 分析方法及應用研究」，碩士論 文，西安建築科技大學，中國，2005。
8. 朱杰江、呂西林、容柏生，「繁雜體型高層結構的推覆分析方法和 應用」，地震工程與工程震動期刊，23(2)：26-36，中國，2003。
9. 何象鏞，「含牆鋼筋混凝土結構側推分析之研究」，博士論文，國立 中央大學，2007。
10. Kilar and Fajfar, “Simple pushover analysis of asymmetric buildings”, Journal of Earthquake Engineering and Structural Dynamics, Vol. 26, pp. 223-249, 1997.
11. Elwood, K. J., Moehle, J. P., “Axial capacity model for shear damaged columns”, ACI Structural Journal, Vol. 102, NO. 4, 578-587, 2005.
12. Elwood, K. J., Moehle, J. P., “Drift capacity of reinforced concrete columns with light transverse reinforcement,” Earthquake Spectra, Vol.
21, No. 1, 71-89, 2005.
13. ACI Committee 318, Building code requirements for structural concrete (ACI 318-05) and commentary (ACI 318R-05). American Concrete Institute, Farmington Hills, MI, 2005.
14. Moehle, J. P.; Elwood, K. J.; and Sezen, H., “Gravity Load Collapse of Building Frames during Earthquakes,” S. M. Uzumeri Symposium:
Behavior and Design of Concrete Structures for Seismic Performance, SP-197, S. A. Sheikh and O. Bayrak, eds., American Concrete Institute, Farmington Hills, pp. 215-238, Mich., 2002.
15. ASCE 41-06, Seismic Rehabilitation of Existing Buildings, ASCE, 2006.
16. 楊元森、謝東儒，「地震工程模擬之 GPGPU 超多核心平行化研 究 」， 國 科 會 多 年 期 專 題 研 究 計 畫 期 中 研 究 成 果 ， NSC 99-2221-E-027-042，2010。
17. Lignos, D.G., and Krawinkler,“Deterioration modeling of steel components in support of collapse prediction of steel moment frames under earthquake loading”, Journal of Structural Engineering, ASCE,
18. Fenves. G. L., ”Annual Workshop on Open System for Earthquake Engineering Simulation”, Pacific Earthquake Engineering Research Center,UC Berkeley, 2001.
19.
http://opensees.berkeley.edu，柏克萊大學地震研究中心。
20. Zhong, “Model-Based Simulation of Reinforced-Concrete Plane Stress Structures, ”Ph.D. dissertation, Dept. of Civil Engineering, University of Huston, 2005.
21. 楊元森、許國志，「非線性動力結構分析之 GPGPU 平行化與效能 評估」，國家地震中心，NCREE-11-023，2011。
22. Elwood, J. P., K. J., Sezen, H., “Gravity Load Collapse of Building Frames during Earthquakes,” S. M. Uzumeri Symposium: Behavior and Design of Concrete Structures for Seismic Performance, SP-197, S. A. Sheikh and O. Bayrak, eds., American Concrete Institute, Farmington Hills, pp. 215-238, Mich., 2002.
圖表附錄
表2.1 RC 柱彎矩塑性鉸之參數[1]
表2.2 RC 柱剪力塑性鉸之參數[1]
表2.3 RC 梁彎矩塑性鉸之參數[15]
表2.4 RC 梁彎矩塑性鉸參數計算表[15]
表2.5 RC 剪力塑性鉸之參數[15]
表2.6 RC 梁剪力塑性鉸參數計算表[15]
表2.7 本研究之 RC 柱彎矩塑性鉸參數
Point Moment/SF Rotation/SF
A 0 0
B 1 0
C 1 0.0166
D 0 0.0263
E 0 0.2626
表2.8 本研究之 RC 梁彎矩塑性鉸參數
Point Moment/SF Rotation/SF
A 0 0
B 1 0
C 1 0.02
D 0.2 0.02
E 0.2 0.03
圖2. 1 側推分析之容量曲線[1]
圖2. 2 雙曲率柱破壞過程[1]
圖2. 3 柱彎矩塑性鉸性質與側向載重位移曲線[1]
圖2. 4軸向破壞時的變位角之修正係數 k'[1]
圖2. 5 柱剪力塑性鉸性質與側向載重位移曲線[1]
圖2. 6 ASCE 41-06建議之 RC 梁彎矩塑性鉸載重位移曲線[15]
圖2. 7 ASCE 41-06建議之 RC 梁剪力塑性鉸載重位移曲線[15]
圖3. 1 OpenSees 主要物件[1]
圖3. 2 OpenSees ModelBuilder Object[19]
圖3. 3 OpenSees Analysis Object[19]
圖3. 4 Bilin Material [17] [19]
圖3. 5 Element Classes Object[19]
圖3. 6 載重模式分類圖[19]
圖3. 7 OpenSees 分析物件
圖3. 8 集中塑性模型
圖3. 9 beam WithHinges command
(a)柱斷面 (b)梁斷面 圖3. 10 梁、柱斷面圖
圖3. 11 RC 柱彎矩塑性鉸性質
圖3. 12 RC 梁彎矩塑性鉸性質
圖3. 13 懸臂梁示意圖與節點編號
圖3. 14 OpenSees 懸臂梁容量曲線圖
圖3. 15 塑性鉸之彎矩與轉角關係圖
圖3. 16 1x1 Frame 示意圖
圖3. 17 1x1 Frame 節點編號與塑性鉸編號圖
圖3. 18 1x1 Frame 容量曲線圖
圖3. 19 1x1 Frame 塑性鉸之彎矩與側推位移圖
1 1
2
3 2 3 3
3
圖3. 20 1x1 Frame 塑性鉸之彎矩與轉角圖
圖3. 21 1x1 Frame 驗證之模型1
圖3. 22 1x1 Frame 驗證之模型2
圖4. 1 1x1 Frame 示意圖
圖4. 2 1x1 Frame 節點編號與桿件性質
圖4. 3 1x1 Frame ETABS 3D 立體圖
圖4. 4 ETABS 梁柱開裂勁度折減修改
圖4. 5 ETABS 剛域修改
圖4. 6 ETABS 梁 End Length Offsets
圖4. 7 ETABS 柱 End Length Offsets
圖4. 8 ETABS 柱之彎矩塑性鉸設定
圖4. 9 ETABS 梁之彎矩塑性鉸設定
圖4. 10 ETABS 塑性鉸分布位置
圖4. 11 ETABS Pushover Curve
圖4. 12 1x1 Frame OpenSees Pushover Curve
1
2 3
1 2 3
圖4. 13 ETABS Step3 變形圖
圖4. 14 柱底端編號 C1-1彎矩塑性鉸 側推位移與彎矩圖
1
1
圖4. 15 柱底端編號 C2-1彎矩塑性鉸 側推位移與彎矩圖
圖4. 16 ETABS Step5 變形圖
1
1
圖4. 17 柱頂端編號 C1-2彎矩塑性鉸 側推位移與彎矩圖
圖4. 18 柱頂端編號 C2-2彎矩塑性鉸 側推位移與彎矩圖
2
2
2
2
圖4. 19 OpenSees 柱彎矩塑性鉸 Moment-Rotation 圖
3 3
3 3
圖4. 20 ETABS 柱彎矩塑性鉸 Moment-Rotation 圖
圖4. 21 2x1 Frame 示意圖
圖4. 22 2x1 ETABS Pushover Curve
圖4. 23 2x1 OpenSees Pushover Curve
1
1 2
2 3
3
圖4. 24 ETABS Step3 變形圖
圖4. 25 編號 B1-1梁彎矩塑性鉸 側推位移與彎矩圖
1
1
圖4. 26 ETABS Step4 變形圖
圖4. 27 柱底端編號 C2-1彎矩塑性鉸 側推位移與彎矩圖
2
2
圖4. 28 ETABS Step5 變形圖
圖4. 29 柱底端編號 C1-1彎矩塑性鉸 側推位移與彎矩圖
3
3
圖4. 30 柱頂端編號 C2-2彎矩塑性鉸 側推位移與彎矩圖
圖4. 31 ETABS Step9 變形圖
3
3
圖4. 32 OpenSees 編號 C2-1之柱彎矩塑性鉸 Moment-Rotation 圖
圖4. 33 ETABS Step10 變形圖
圖4. 34 柱頂端編號 C1-2彎矩塑性鉸 側推位移與彎矩圖
圖4. 35 ETABS Step11 變形圖
圖4. 36 ETABS Step12 變形圖
圖4. 37 2x8 Frame 示意圖
圖4. 38 OpenSees 2x8 Frame Pushover Curve
圖4. 39 ETABS 2x8 Frame Pushover Curve
1
2
3
4
圖4. 40 ETABS 2x8 Frame Pushover Curve table
圖4. 41 2x8 Frame ETABS Step14 變形圖
圖4. 42 編號 B1-1梁正彎矩塑性鉸之側推位移與彎矩
圖4. 43 2x8 Frame ETABS Step15 變形圖
圖4. 44 編號 C2~C6-1、C8-1 柱彎矩塑性鉸之側推位移與彎矩圖
圖4. 45 2x8 Frame ETABS Step28 變形圖
圖4. 46 2x8 Frame ETABS Step31 變形圖
圖4. 47 編號 C2-1柱之彎矩塑性鉸 轉角與側推位移
圖4. 48 2x8 Frame ETABS Step32 變形圖
圖4. 49 2x8 Frame ETABS Step41 變形圖
圖4. 50 2x8 Frame ETABS Step49 變形圖
側推範例程式碼 1x1
圖 1 OpenSees 1x1模型 節點編號與桿件性質
1. wipe
2. #--- 3. # Concentrated -1X1-Frame 4. #--- 5. model BasicBuilder -ndm 2 -ndf 3 6. #--- 7. # Define Parameter 8. #--- 9. #--- 10. # Define Column Plastic Hinge
11. #--- 12. set Lc 360.0 ;
13. set bc 50.0 ; 14. set hc 30.0 ; 15. set fc_c 160.0 ; 16. set offset_c 60 ;
17. set Ec_c [expr 15000.0*sqrt($fc_c)] ;
54. set Ec_b [expr 15000.0*sqrt($fc_b)] ; 55. set E_b 204000000000
56. set Ab [expr $bb*$hb] ; 63. set My_pb1 1845540;
64. set My_nb1 -1845540;
65. set My_pb2 4407730;
66. set My_nb2 -4407730;
67. set LS_b 1000;
94. node 21 0.0 [expr 360.0-$offset_c]
95. node 34 400.0 0
96. node 43 400.0 [expr 360.0-$offset_c]
97. #--- 98. node 24 [expr 0.0+$offset_b] 360.0 99. node 42 [expr 400.0-$offset_b] 360.0 100.
101. node 2121 0.0 [expr 360.0-$offset_c]
102. node 4343 400.0 [expr 360.0-$offset_c]
103. node 2424 [expr 0.0+$offset_b] 360.0 104. node 4242 [expr 400.0-$offset_b] 360.0 105. #---
113. uniaxialMaterial Bilin 1 $E_c $u_pc $u_nc $My_pc $My_nc $LS_c $LC_c $LA_c $LK_c
$cS_c $cC_c $cA_c $cK_c $th_P_pc $th_P_nc $th_PC_pc $th_PC_nc $th_R_pc $th_R_nc
$th_U_pc $th_U_nc $DP_c $DN_c 114. #--- 115. # Material Beam 116. #---
117. uniaxialMaterial Bilin 2 $E_b $u_pb $u_nb $My_pb1 $My_nb1 $LS_b $LC_b $LA_b $LK_b
$cS_b $cC_b $cA_b $cK_b $th_P_pb $th_P_nb $th_PC_pb $th_PC_nb $th_R_pb $th_R_nb
$th_U_pb $th_U_nb $DP_b $DN_b
118. uniaxialMaterial Bilin 3 $E_b $u_pb $u_nb $My_pb2 $My_nb2 $LS_b $LC_b $LA_b $LK_b
$cS_b $cC_b $cA_b $cK_b $th_P_pb $th_P_nb $th_PC_pb $th_PC_nb $th_R_pb $th_R_nb
$th_U_pb $th_U_nb $DP_b $DN_b 119.
130. element elasticBeamColumn 1221 12 2121 $Ac $Ec_c $IZc_cr 1 131. element elasticBeamColumn 3443 34 4343 $Ac $Ec_c $IZc_cr 1 132. #rigid
133. element elasticBeamColumn 221 2 21 $Ac $E_c $IZc_cr 2 134. element elasticBeamColumn 443 4 43 $Ac $E_c $IZc_cr 2 135. #rigid
141. element elasticBeamColumn 2442 2424 4242 $Ab $Ec_b $IZb_cr 2 142. #rigid
143. element elasticBeamColumn 224 2 24 $Ab $E_b $IZb_cr 2 144. element elasticBeamColumn 442 4 42 $Ab $E_b $IZb_cr 2 145. #rigid
146. #--- 147. # Rotational Springs 148. #--- 149. element zeroLength 12 1 12 -mat 1 -dir 6 150. element zeroLength 21 21 2121 -mat 1 -dir 6 151. element zeroLength 24 24 2424 -mat 2 -dir 6 152. element zeroLength 34 3 34 -mat 1 -dir 6 153. element zeroLength 42 42 4242 -mat 3 -dir 6 154. element zeroLength 43 43 4343 -mat 1 -dir 6
155. #---
176. recorder Node -file Node2.out -time -node 2 -dof 1 2 3 disp 177. recorder Element -file MomentC.out -region 1 force 178. recorder Element -file RotationC.out -region 1 deformation 179. recorder Element -file MomentB.out -region 2 force 180. recorder Element -file RotationB.out -region 2 deformation 181. #---
198. loadConst -time 0.0
199. #--- 200. # Lateral-Load 201. #--- 202. pattern Plain 2 Linear {
203. load 2 60000 0 0 204. }
205. #--- 206. # Pushover-Analysis 207. #--- 208. constraints Plain
209. test EnergyIncr 1.0e-8 6 200
210. integrator DisplacementControl 2 1 0.05 211. algorithm Newton
212. numberer Plain 213. system BandGeneral 214. analysis Static 215. analyze 200 216.
217. #--- 218. # Display Data 219. #--- 220. print node
221. print ele
2x1
圖 2 OpenSees 2x1模型 節點編號與桿件性質
1. wipe
17. set Ec_c [expr 15000.0*sqrt($fc_c)] ; 18. set E_c 204000000000;
49. set Lb 400.0 ; 50. set bb 30.0 ; 51. set hb 60.0 ; 52. set fc_b 160.0 ; 53. set offset_b 15;
54. set Ec_b [expr 15000.0*sqrt($fc_b)] ; 55. set E_b 204000000000
56. set Ab [expr $bb*$hb] ; 63. set My_pb1 1845540;
64. set My_nb1 -1845540;
65. set LS_b 1000; 86. set My_pb2 4407730;
87. set My_nb2 -4407730;
88. #---
97. #--- 98. #colum plastic hinge
99. node 12 0.0 0 113. #beam plastic hinge
114. node 24 [expr 0.0+$offset_b] 360.0 115. node 42 [expr 400.0-$offset_b] 360.0 116. node 56 [expr 0.0+$offset_b] 720.0 117. node 65 [expr 400.0-$offset_b] 720.0 118. node 2424 [expr 0.0+$offset_b] 360.0 119. node 4242 [expr 400.0-$offset_b] 360.0 120. node 5656 [expr 0.0+$offset_b] 720.0 121. node 6565 [expr 400.0-$offset_b] 720.0 122.
131. uniaxialMaterial Bilin 1 $E_c $u_pc $u_nc $My_pc $My_nc $LS_c $LC_c $LA_c
$LK_c $cS_c $cC_c $cA_c $cK_c $th_P_pc $th_P_nc $th_PC_pc $th_PC_nc $th_R_pc
$th_R_nc $th_U_pc $th_U_nc $DP_c $DN_c 132. #--- 133. # Material Beam 134. #---
135. uniaxialMaterial Bilin 2 $E_b $u_pb $u_nb $My_pb1 $My_nb1 $LS_b $LC_b $LA_b
$LK_b $cS_b $cC_b $cA_b $cK_b $th_P_pb $th_P_nb $th_PC_pb $th_PC_nb $th_R_pb
$th_R_nb $th_U_pb $th_U_nb $DP_b $DN_b 136.
137. uniaxialMaterial Bilin 3 $E_b $u_pb $u_nb $My_pb2 $My_nb2 $LS_b $LC_b $LA_b
$LK_b $cS_b $cC_b $cA_b $cK_b $th_P_pb $th_P_nb $th_PC_pb $th_PC_nb $th_R_pb
$th_R_nb $th_U_pb $th_U_nb $DP_b $DN_b 138.
139. #---
149. element elasticBeamColumn 1221 12 2121 $Ac $Ec_c $IZc_cr 1 150. element elasticBeamColumn 3443 34 4343 $Ac $Ec_c $IZc_cr 1 151. element elasticBeamColumn 2552 25 5252 $Ac $Ec_c $IZc_cr 1 152. element elasticBeamColumn 4664 46 6464 $Ac $Ec_c $IZc_cr 1 153. #rigid
154. element elasticBeamColumn 221 2 21 $Ac $E_c $IZc_cr 2 155. element elasticBeamColumn 443 4 43 $Ac $E_c $IZc_cr 2 156.
157. element elasticBeamColumn 664 6 64 $Ac $E_c $IZc_cr 2 158. element elasticBeamColumn 552 5 52 $Ac $E_c $IZc_cr 2 159. #rigid
165. element elasticBeamColumn 2442 2424 4242 $Ab $Ec_b $IZb_cr 2 166. element elasticBeamColumn 5665 5656 6565 $Ab $Ec_b $IZb_cr 2 167. #rigid
168. element elasticBeamColumn 224 2 24 $Ab $E_b $IZb_cr 2 169. element elasticBeamColumn 442 4 42 $Ab $E_b $IZb_cr 2 170. element elasticBeamColumn 556 5 56 $Ab $E_b $IZb_cr 2 171. element elasticBeamColumn 665 6 65 $Ab $E_b $IZb_cr 2 172. #rigid
173. #--- 174. # Rotational Springs 175. #--- 176. element zeroLength 112 1 12 -mat 1 -dir 6 177. element zeroLength 212121 21 2121 -mat 1 -dir 6 178. element zeroLength 242424 24 2424 -mat 2 -dir 6 179. element zeroLength 334 3 34 -mat 1 -dir 6 180. element zeroLength 424242 42 4242 -mat 3 -dir 6 181. element zeroLength 434343 43 4343 -mat 1 -dir 6 182. element zeroLength 225 2 25 -mat 1 -dir 6 183. element zeroLength 525252 52 5252 -mat 1 -dir 6 184. element zeroLength 446 4 46 -mat 1 -dir 6 185. element zeroLength 646464 64 6464 -mat 1 -dir 6 186. element zeroLength 565656 56 5656 -mat 2 -dir 6
187. element zeroLength 656565 65 6565 -mat 3 -dir 6
219. region 1 -ele 112 212121 334 434343 225 525252 446 646464 220. #---
221. # Region Beam 222. #--- 223. region 2 -ele 242424 424242 565656 656565 224. #--- 225. # Recorder
226. #---
227. recorder Node -file Node5.out -time -node 5 -dof 1 2 3 disp 228. recorder Element -file MomentC.out -region 1 force 229. recorder Element -file RotationC.out -region 1 deformation 230. recorder Element -file MomentB.out -region 2 force 231. recorder Element -file RotationB.out -region 2 deformation 232. #---
233. # Gravity
234. #---
235. pattern Plain 1 Linear {
249. loadConst -time 0.0
250. #---
263. algorithm Newton 264. numberer Plain
2x8
17. set Ec_c [expr 15000.0*sqrt($fc_c)] ; 18. set E_c 204000000000;
47. # Define Beam Plastic Hinge 1
54. set Ec_b [expr 15000.0*sqrt($fc_b)] ; 55. set E_b 204000000000
56. set Ab [expr $bb*$hb] ; 63. set My_pb1 1845540;
64. set My_nb1 -1845540;
65. set LS_b 1000; 86. set My_pb2 4407730;
87. set My_nb2 -4407730;
88. #---
95. #plastic hinge 96. node 12 0 0
97. node 21 0 [expr 360-$offset_c]
98. node 2121 0 [expr 360-$offset_c]
99. node 23 0 360
100. node 32 0 [expr 720-$offset_c]
101. node 3232 0 [expr 720-$offset_c]
102.
103. #C2
104. node 4 400 0 105. node 5 400 360 106. node 6 400 720 107. #plastic hinge 108. node 45 400 0
109. node 54 400 [expr 360-$offset_c]
110. node 5454 400 [expr 360-$offset_c]
111. node 56 400 360
112. node 65 400 [expr 720-$offset_c]
113. node 6565 400 [expr 720-$offset_c]
114.
115. #C3
116. node 7 800 0 117. node 8 800 360 118. node 9 800 720 119. #plastic hinge 120. node 78 800 0
121. node 87 800 [expr 360-$offset_c]
122. node 8787 800 [expr 360-$offset_c]
123. node 89 800 360
124. node 98 800 [expr 720-$offset_c]
125. node 9898 800 [expr 720-$offset_c]
126.
127. #C4
128. node 10 1200 0 129. node 11 1200 360 130. node 121212 1200 720 131. #plastic hinge
132. node 1011 1200 0
133. node 1110 1200 [expr 360-$offset_c]
134. node 11102 1200 [expr 360-$offset_c]
135. node 1112 1200 360
136. node 1211 1200 [expr 720-$offset_c]
137. node 12112 1200 [expr 720-$offset_c]
138.
139. #C5
140. node 13 1600 0 141. node 14 1600 360 142. node 15 1600 720
143. #plastic hinge 144. node 1314 1600 0
145. node 1413 1600 [expr 360-$offset_c]
146. node 14132 1600 [expr 360-$offset_c]
147. node 1415 1600 360
148. node 1514 1600 [expr 720-$offset_c]
149. node 15142 1600 [expr 720-$offset_c]
150.
151. #C6
152. node 16 2000 0 153. node 17 2000 360 154. node 18 2000 720 155. #plastic hinge 156. node 1617 2000 0
157. node 1716 2000 [expr 360-$offset_c]
158. node 17162 2000 [expr 360-$offset_c]
159. node 1718 2000 360
160. node 1817 2000 [expr 720-$offset_c]
161. node 18172 2000 [expr 720-$offset_c]
162.
163. #C7
164. node 19 2400 0 165. node 20 2400 360 166. node 212121 2400 720 167. #plastic hinge
168. node 1920 2400 0
169. node 2019 2400 [expr 360-$offset_c]
170. node 20192 2400 [expr 360-$offset_c]
171. node 2021 2400 360
172. node 2120 2400 [expr 720-$offset_c]
173. node 21202 2400 [expr 720-$offset_c]
174.
175. #C8
176. node 22 2800 0 177. node 232323 2800 360 178. node 24 2800 720 179. #plastic hinge 180. node 2223 2800 0
181. node 2322 2800 [expr 360-$offset_c]
182. node 23222 2800 [expr 360-$offset_c]
183. node 2324 2800 360
184. node 2423 2800 [expr 720-$offset_c]
185. node 24232 2800 [expr 720-$offset_c]
186.
187. #C9
188. node 252525 3200 0 189. node 26 3200 360 190. node 27 3200 720
191. #plastic hinge 192. node 2526 3200 0
193. node 2625 3200 [expr 360-$offset_c]
194. node 26252 3200 [expr 360-$offset_c]
195. node 2627 3200 360
196. node 2726 3200 [expr 720-$offset_c]
197. node 27262 3200 [expr 720-$offset_c]
198. #--- 199. #B1 plastic hinge
200. node 25 [expr 0.0+$offset_b] 360.0 201. node 52 [expr 400.0-$offset_b] 360.0 202. node 2525 [expr 0.0+$offset_b] 360.0 203. node 5252 [expr 400.0-$offset_b] 360.0 204.
205. node 58 [expr 400.0+$offset_b] 360.0 206. node 85 [expr 800.0-$offset_b] 360.0 207. node 5858 [expr 400.0+$offset_b] 360.0 208. node 8585 [expr 800.0-$offset_b] 360.0 209.
210. node 811 [expr 800.0+$offset_b] 360.0 211. node 118 [expr 1200.0-$offset_b] 360.0 212. node 8112 [expr 800.0+$offset_b] 360.0 213. node 1182 [expr 1200.0-$offset_b] 360.0 214.
215. node 1114 [expr 1200.0+$offset_b] 360.0 216. node 1411 [expr 1600.0-$offset_b] 360.0 217. node 11142 [expr 1200.0+$offset_b] 360.0 218. node 14112 [expr 1600.0-$offset_b] 360.0 219.
220. node 1417 [expr 1600.0+$offset_b] 360.0 221. node 1714 [expr 2000.0-$offset_b] 360.0 222. node 14172 [expr 1600.0+$offset_b] 360.0 223. node 17142 [expr 2000.0-$offset_b] 360.0 224.
225. node 1720 [expr 2000.0+$offset_b] 360.0 226. node 2017 [expr 2400.0-$offset_b] 360.0 227. node 17202 [expr 2000.0+$offset_b] 360.0 228. node 20172 [expr 2400.0-$offset_b] 360.0 229.
230. node 2023 [expr 2400.0+$offset_b] 360.0 231. node 2320 [expr 2800.0-$offset_b] 360.0 232. node 20232 [expr 2400.0+$offset_b] 360.0 233. node 23202 [expr 2800.0-$offset_b] 360.0 234.
235. node 2326 [expr 2800.0+$offset_b] 360.0 236. node 2623 [expr 3200.0-$offset_b] 360.0 237. node 23262 [expr 2800.0+$offset_b] 360.0 238. node 26232 [expr 3200.0-$offset_b] 360.0
239.
240. #B2 plastic hinge
241. node 36 [expr 0.0+$offset_b] 720.0 242. node 63 [expr 400.0-$offset_b] 72.0 243. node 3636 [expr 0.0+$offset_b] 720.0 244. node 6363 [expr 400.0-$offset_b] 720.0 245.
246. node 69 [expr 400.0+$offset_b] 720.0 247. node 96 [expr 800.0-$offset_b] 720.0 248. node 6969 [expr 400.0+$offset_b] 720.0 249. node 9696 [expr 800.0-$offset_b] 720.0 250.
251. node 912 [expr 800.0+$offset_b] 720.0 252. node 129 [expr 1200.0-$offset_b] 720.0 253. node 9122 [expr 800.0+$offset_b] 720.0 254. node 1292 [expr 1200.0-$offset_b] 720.0 255.
256. node 1215 [expr 1200.0+$offset_b] 720.0 257. node 1512 [expr 1600.0-$offset_b] 720.0 258. node 12152 [expr 1200.0+$offset_b] 720.0 259. node 15122 [expr 1600.0-$offset_b] 720.0 260.
261. node 1518 [expr 1600.0+$offset_b] 720.0 262. node 1815 [expr 2000.0-$offset_b] 720.0 263. node 15182 [expr 1600.0+$offset_b] 720.0 264. node 18152 [expr 2000.0-$offset_b] 720.0 265.
266. node 1821 [expr 2000.0+$offset_b] 720.0 267. node 2118 [expr 2400.0-$offset_b] 720.0 268. node 18212 [expr 2000.0+$offset_b] 720.0 269. node 21182 [expr 2400.0-$offset_b] 720.0 270.
271. node 2124 [expr 2400.0+$offset_b] 720.0 272. node 2421 [expr 2800.0-$offset_b] 720.0 273. node 21242 [expr 2400.0+$offset_b] 720.0 274. node 24212 [expr 2800.0-$offset_b] 720.0 275.
276. node 2427 [expr 2800.0+$offset_b] 720.0 277. node 2724 [expr 3200.0-$offset_b] 720.0 278. node 24272 [expr 2800.0+$offset_b] 720.0 279. node 27242 [expr 3200.0-$offset_b] 720.0 280.
287. fix 10 1 1 1
296. uniaxialMaterial Bilin 1 $E_c $u_pc $u_nc $My_pc $My_nc $LS_c $LC_c $LA_c
$LK_c $cS_c $cC_c $cA_c $cK_c $th_P_pc $th_P_nc $th_PC_pc $th_PC_nc $th_R_pc
$th_R_nc $th_U_pc $th_U_nc $DP_c $DN_c 297. #--- 298. # Material Beam 299. #---
300. uniaxialMaterial Bilin 2 $E_b $u_pb $u_nb $My_pb1 $My_nb1 $LS_b $LC_b
$LA_b $LK_b $cS_b $cC_b $cA_b $cK_b $th_P_pb $th_P_nb $th_PC_pb $th_PC_nb
$th_R_pb $th_R_nb $th_U_pb $th_U_nb $DP_b $DN_b 301.
302. uniaxialMaterial Bilin 3 $E_b $u_pb $u_nb $My_pb2 $My_nb2 $LS_b $LC_b
$LA_b $LK_b $cS_b $cC_b $cA_b $cK_b $th_P_pb $th_P_nb $th_PC_pb $th_PC_nb
$th_R_pb $th_R_nb $th_U_pb $th_U_nb $DP_b $DN_b 303.
314. element elasticBeamColumn 1 12 2121 $Ac $Ec_c $IZc_cr 1 315. element elasticBeamColumn 10 23 3232 $Ac $Ec_c $IZc_cr 1 316. #rigid
317. element elasticBeamColumn 35 2 21 $Ac $E_c $IZc_cr 2 318. element elasticBeamColumn 44 3 32 $Ac $E_c $IZc_cr 2 319.
320. #C2 321. #Elastic
322. element elasticBeamColumn 2 45 5454 $Ac $Ec_c $IZc_cr 1 323. element elasticBeamColumn 11 56 6565 $Ac $Ec_c $IZc_cr 1 324. #rigid
325. element elasticBeamColumn 36 5 54 $Ac $E_c $IZc_cr 2 326. element elasticBeamColumn 45 6 65 $Ac $E_c $IZc_cr 2 327.
328. #C3
329. #Elastic
330. element elasticBeamColumn 3 78 8787 $Ac $Ec_c $IZc_cr 1 331. element elasticBeamColumn 12 89 9898 $Ac $Ec_c $IZc_cr 1 332. #rigid
333. element elasticBeamColumn 37 8 87 $Ac $E_c $IZc_cr 2 334. element elasticBeamColumn 46 9 98 $Ac $E_c $IZc_cr 2 335.
336. #C4 337. #Elastic
338. element elasticBeamColumn 4 1011 11102 $Ac $Ec_c $IZc_cr 1 339. element elasticBeamColumn 13 1112 12112 $Ac $Ec_c $IZc_cr 1 340. #rigid
341. element elasticBeamColumn 38 11 1110 $Ac $E_c $IZc_cr 2 342. element elasticBeamColumn 47 121212 1211 $Ac $E_c $IZc_cr 2 343.
344. #C5 345. #Elastic
346. element elasticBeamColumn 5 1314 14132 $Ac $Ec_c $IZc_cr 1 347. element elasticBeamColumn 14 1415 15142 $Ac $Ec_c $IZc_cr 1 348. #rigid
349. element elasticBeamColumn 39 14 1413 $Ac $E_c $IZc_cr 2 350. element elasticBeamColumn 48 15 1514 $Ac $E_c $IZc_cr 2 351.
352. #C6 353. #Elastic
354. element elasticBeamColumn 6 1617 17162 $Ac $Ec_c $IZc_cr 1 355. element elasticBeamColumn 15 1718 18172 $Ac $Ec_c $IZc_cr 1 356. #rigid
357. element elasticBeamColumn 40 17 1716 $Ac $E_c $IZc_cr 2 358. element elasticBeamColumn 49 18 1817 $Ac $E_c $IZc_cr 2 359.
360. #C7 361. #Elastic
362. element elasticBeamColumn 7 1920 20192 $Ac $Ec_c $IZc_cr 1 363. element elasticBeamColumn 16 2021 21202 $Ac $Ec_c $IZc_cr 1 364. #rigid
365. element elasticBeamColumn 41 20 2019 $Ac $E_c $IZc_cr 2 366. element elasticBeamColumn 50 212121 2120 $Ac $E_c $IZc_cr 2 367.
368. #C8 369. #Elastic
370. element elasticBeamColumn 8 2223 23222 $Ac $Ec_c $IZc_cr 1 371. element elasticBeamColumn 17 2324 24232 $Ac $Ec_c $IZc_cr 1 372. #rigid
373. element elasticBeamColumn 42 232323 2322 $Ac $E_c $IZc_cr 2 374. element elasticBeamColumn 51 24 2423 $Ac $E_c $IZc_cr 2 375.
376. #C9
377. #Elastic
378. element elasticBeamColumn 9 2526 26252 $Ac $Ec_c $IZc_cr 1 379. element elasticBeamColumn 18 2627 27262 $Ac $Ec_c $IZc_cr 1 380. #rigid
381. element elasticBeamColumn 43 26 2625 $Ac $E_c $IZc_cr 2 382. element elasticBeamColumn 52 27 2726 $Ac $E_c $IZc_cr 2 383.
389. element elasticBeamColumn 19 2525 5252 $Ab $Ec_b $IZb_cr 2 390. element elasticBeamColumn 20 5858 8585 $Ab $Ec_b $IZb_cr 2 391. element elasticBeamColumn 21 8112 1182 $Ab $Ec_b $IZb_cr 2 392. element elasticBeamColumn 22 11142 14112 $Ab $Ec_b $IZb_cr 2 393. element elasticBeamColumn 23 14172 17142 $Ab $Ec_b $IZb_cr 2 394. element elasticBeamColumn 24 17202 20172 $Ab $Ec_b $IZb_cr 2 395. element elasticBeamColumn 25 20232 23202 $Ab $Ec_b $IZb_cr 2 396. element elasticBeamColumn 26 23262 26232 $Ab $Ec_b $IZb_cr 2 397. #rigid
398. element elasticBeamColumn 53 2 25 $Ab $E_b $IZb_cr 2 399. element elasticBeamColumn 54 5 52 $Ab $E_b $IZb_cr 2 400.
401. element elasticBeamColumn 55 5 58 $Ab $E_b $IZb_cr 2 402. element elasticBeamColumn 56 8 85 $Ab $E_b $IZb_cr 2 403.
404. element elasticBeamColumn 57 8 811 $Ab $E_b $IZb_cr 2 405. element elasticBeamColumn 58 11 118 $Ab $E_b $IZb_cr 2 406.
407. element elasticBeamColumn 59 11 1114 $Ab $E_b $IZb_cr 2 408. element elasticBeamColumn 60 14 1411 $Ab $E_b $IZb_cr 2 409.
410. element elasticBeamColumn 61 14 1417 $Ab $E_b $IZb_cr 2 411. element elasticBeamColumn 62 17 1714 $Ab $E_b $IZb_cr 2 412.
413. element elasticBeamColumn 63 17 1720 $Ab $E_b $IZb_cr 2 414. element elasticBeamColumn 64 20 2017 $Ab $E_b $IZb_cr 2 415.
416. element elasticBeamColumn 65 20 2023 $Ab $E_b $IZb_cr 2 417. element elasticBeamColumn 66 232323 2320 $Ab $E_b $IZb_cr 2 418.
419. element elasticBeamColumn 67 232323 2326 $Ab $E_b $IZb_cr 2 420. element elasticBeamColumn 68 26 2623 $Ab $E_b $IZb_cr 2 421. #rigid
422.
423. #B2 424. #Elastic
425. element elasticBeamColumn 27 3636 6363 $Ab $Ec_b $IZb_cr 2 426. element elasticBeamColumn 28 6969 9696 $Ab $Ec_b $IZb_cr 2 427. element elasticBeamColumn 29 9122 1292 $Ab $Ec_b $IZb_cr 2 428. element elasticBeamColumn 30 12152 15122 $Ab $Ec_b $IZb_cr 2 429. element elasticBeamColumn 31 15182 18152 $Ab $Ec_b $IZb_cr 2 430. element elasticBeamColumn 32 18212 21182 $Ab $Ec_b $IZb_cr 2 431. element elasticBeamColumn 33 21242 24212 $Ab $Ec_b $IZb_cr 2 432. element elasticBeamColumn 34 24272 27242 $Ab $Ec_b $IZb_cr 2 433.
434. #rigid
435. element elasticBeamColumn 69 3 36 $Ab $E_b $IZb_cr 2
435. element elasticBeamColumn 69 3 36 $Ab $E_b $IZb_cr 2