核能高效率及穩定的發電模式，仍是全球供電系統中不可或缺的一部分，其經濟且便利的背後也伴隨著高度的風險。台灣與日本同處板塊交界帶，地震活動頻繁，核能電廠受震時的安全性問題變得愈加重要，核子反應爐內槽狀容器會因內部液體的晃動模式造成儲槽結構震盪放大，更容易破壞導致嚴重的後果。
本研究利用國家地震工程研究中心(NCREE)之大型雙軸向剪力盒，考量儲槽之流固耦合及土壤結構互制效應進行振動台試驗，並觀察儲槽之流體晃動模式造成之影響，以及探討於土壤顯著頻率下土層的加速度放大效應。
使用有限元素法分析軟體ABAQUS/Explicit來進行時間域的顯性動力分析，透過ABAQUS有限元素軟體建立數值模型，依據土壤模擬方式不同分為實體有限元素土壤模型及土壤彈簧模型，除了進行儲槽流固耦合CEL分析外，同時考量土壤結構互制效應來如實呈現試驗時儲槽受震的狀態。首先以ABAQUS之模態分析求解模型之顯著頻率與試驗進行比對，驗證模型之準確性，再進行重力平衡的擬靜態分析，將重力平衡完成之模型分別以接近流體顯著頻率及接近砂土顯著頻率的正弦波地震歷時輸入於模型底部，分析結構物之土層加速度放大效應以及流體晃動模式，並與振動台試驗結果進行比較。

The high-efficiency and stable power generation mode of nuclear energy is still an indispensable part of the global power supply system, and its economy and convenience are also accompanied by a high degree of risk. Taiwan and Japan are located at the same plate boundary zone, and earthquake activities are frequent. The safety of nuclear power plants when earthquakes become more important. The tank container in the nuclear reactor will be amplified by the shaking mode of the liquid inside, and it is easier to damage. In this reserch, the large-scale biaxial shear box of National Center for Research on Earthquakes Engineering (NCREE) was used to conduct shaking table tests in consideration of the fluid-solid coupling and soil structure interaction effects of the storage tank, and to observe the impact of the fluid sloshing mode of the storage tank.
The finite element method analysis software ABAQUS/Explicit is used for time-domain explicit dynamic analysis. According to the different methods of soil simulation, it is divided into physical finite element soil model and soil spring model. In addition to the fluid-solid coupling CEL analysis, the interaction effect of the soil structure is also considered, and the state of the storage tank under the seismic action during the test is truly presented. First, the significant frequency of the model is compared with the ABAQUS modal analysis solution to verify the accuracy of the model, and then a quasi-static analysis is performed. The quasi-static analysis model is input into sine waves close to the effective frequency of the fluid and the effective frequency of the sand. Analyze the soil acceleration amplification effect of the structure and fluid sloshing mode, and compare with the shaking table test results.

[1] P. K.Malhotra, “Earthquake induced sloshing in tanks with insufficient freeboard,” Struct. Eng. Int. J. Int. Assoc. Bridg. Struct. Eng., vol. 16, no. 3, pp. 222–225, 2006, doi: 10.2749/101686606778026466.
[2] D.Watkins andM.Tobolski, “Modal Combination of Seismically Induced,” 2013.
[3] S.Lee, B.Kim, andY. J.Lee, “Seismic Fragility Analysis of Steel Liquid Storage Tanks Using Earthquake Ground Motions Recorded in Korea,” Math. Probl. Eng., vol. 2019, 2019, doi: 10.1155/2019/6190159.
[4] F.Paolacci, H. N.Phan, D.Corritore, S.Alessandri, O. S.Bursi, andM. S.Reza, “Seismic fragility analysis of steel storage tanks,” COMPDYN 2015 - 5th ECCOMAS Themat. Conf. Comput. Methods Struct. Dyn. Earthq. Eng., no. May, pp. 2054–2065, 2015, doi: 10.7712/120115.3522.1040.
[5] M.Sivý andM.Musil, “Seismic resistance of storage tanks containing liquid in accordance with principles of Eurocode 8 standard,” Stroj. Cas., vol. 66, no. 2, pp. 79–88, 2016, doi: 10.1515/scjme-2016-0021.
[6] K.Bandyopadhyay, A.Cornell, C.Costantino, R.Kennedy, C.Miller, andA.Veletsos, “Seismic Design and Evaluation Guidelines for the Department of Energy High-Level Waste Storage Tanks and Appurtenances,” vol. 52361, 1995.
[7] A. C. I.Committee, “Seismic Design of Liquid-Containing ( ACI 350 . 3-06 ),” Main, 2007.
[8] T.Ibata, I.Nakachi, K.Ishida, andJ.Yokozawa, “Damage to storage tanks caused by the 2011 Tohoku earthquake and tsunami and proposal for structural assessment method for cylindrical storage tanks,” IGT Int. Liq. Nat. Gas Conf. Proc., vol. 2, pp. 820–837, 2013.
[9] 許皓程, “馬鞍山核電廠圍阻體土壤結構互制之動力分析,” 國立成功大學, 2018.
[10] G.-S. L.Po-Chi Lo, “Analysis of Soil-structure Interaction By Sub-structure Method,” 2003.
[11] 葉鼎盛, “離岸風機結構與樁土互制之有限元素自振分析,” 國立成功大學, 2015.
[12] K.Tokimatsu, H.Suzuki, andM.Sato, “Effects of inertial and kinematic interaction on seismic behavior of pile with embedded foundation,” Soil Dyn. Earthq. Eng., vol. 25, no. 7–10, pp. 753–762, 2005, doi: 10.1016/j.soildyn.2004.11.018.
[13] 翁作新、陳家漢、彭立先、李偉誠, “大型振動台剪力盒土壤液化試驗(II)－大型砂試體之準備與振動台初期試驗,” 2003.
[14] M.Ashiquzzaman andK. J.Hong, “Simplified Model of Soil-Structure Interaction for Seismically Isolated Containment Buildings in Nuclear Power Plant,” Structures, vol. 10, pp. 209–218, 2017, doi: 10.1016/j.istruc.2016.09.014.
[15] 陳正興、柯永彥、許尚逸, “核能一廠乾式貯存設施結構地震之安全審查與確認分析,” 2006.
[16] R.Rajapakse, “Laterally loaded piles,” Pile Des. Constr. Rules Thumb, pp. 235–240, 2016, doi: 10.1016/b978-0-12-804202-1.00017-6.
[17] API, “Recommended Practice for Planning , Designing and Constructing Fixed Offshore Platforms — Working Stress Design,” Api Recomm. Pract., vol. 24-WSD, no. December 2000, p. 242, 2007, [Online]. Available: http://scholar.google.com/scholar?hl=en&btnG=Search&q=intitle:Recommended+Practice+for+planning,+designing+and+constructing+fixed+Offshore+Platforms+-+Working+stress+design#0.
[18] 宋学官, “流固耦合分析基礎,” in 流固耦合分析与工程实例, 水利水电出版社, 2012, pp. 1–13.
[19] S.B., “Coupled Eulerian Lagrange (CEL) Analysis with ABAQUS Eulerian and Lagrangian Approach.” .
[20] J.Jablonski, P.Carlucci, R.Thyagarajan, B.Nandi, andJ.Arata, “Simulating Underbelly Blast Events using Abaqus/Explicit -CEL,” no. December, p. 14, 2012.
[21] I.Smojver andD.Ivančević, “Coupled Euler Lagrangian approach using Abaqus/explicit in the bird strike aircraft damage analysis,” 2010 SIMULIA Cust. Conf., pp. 1–14, 2010.
[22] G.Qiu, S.Henke, andJ.Grabe, “Applications of Coupled Eulerian-Lagrangian Method to Geotechnical Problems with Large Deformations,” SIMULIA Cust. Conf., no. 2001, pp. 1–16, 2009.
[23] Dassault Systèmes Simulia, “Abaqus Analysis User’s Guide, vol4,” in ABAQUS 6.14 Analysis User’s Guide, vol. IV, 2014.
[24] 詹琍尹, “非韌性鋼筋混凝土結構受近斷層地震力之三維有限元素非線性分析,” 國立成功大學, 2019.
[25] Dassault Systèmes Simulia, “Abaqus Analysis User’s Guide, vol2,” in ABAQUS 6.14 Analysis User’s Guide, vol. II, 2014.
[26] 鄭鈞, “以實例說明隱式與顯式求解法 - YouTube,” Simutech Solution Corp. .
[27] I.Chowdhury andS. P.Dasgupta, “Computation of Rayleigh damping coefficients for large systems,” Electron. J. Geotech. Eng., vol. 8 C, 2003.
[28] W.Liu andX. J.Yang, “Damage evolution with growing cyclic creep and life prediction of MDYB-3 PMMA,” Fatigue Fract. Eng. Mater. Struct., vol. 36, no. 6, pp. 483–491, 2013, doi: 10.1111/ffe.12017.
[29] E.Engineering, “樁基礎沖刷橋梁模型之 振動台試驗研究 (i),” no. I.
[30] S.SAWADA, “A Simplified Equation to Approximate Natural Period of Layered Ground on the Elastic Bedrock for Seismic Design of Structures,” 13th World Conf. Earthq. Eng., no. 1100, p. Paper No. 1100, 2004.
[31] W. F.Robert, T. M.Alan, andJ. P.Philip, “Introduction to fluid mechanics, Fifth edition,” Introduction to Fluid Mechanics, Fifth Edition. pp. 1–745, 2015.
[32] J. D.Tippmann, S. C.Prasad, andP. N.Shah, “2-D Tank Sloshing Using the Coupled Eulerian- LaGrangian ( CEL ) Capability of Abaqus / Explicit,” Simulia Cust. Conf., pp. 1–11, 2009.
[33] E.Engineering, “單跨樁基礎橋梁模型之 振動台試驗研究.”
[34] Dassault Systèmes Simulia, “Abaqus Analysis User’s Guide, vol5,” ABAQUS 6.14 Anal. User’s Guid., vol. V, p. 946, 2014.
[35] X.Chen, J.Duan, andY.Li, “Mass proportional damping in nonlinear time-history analysis,” no. Ic3me, pp. 567–571, 2015, doi: 10.2991/ic3me-15.2015.112.
[36] A. G.Brady, “STUDIES OF RESPONSE TO EARTHQUAKE GROUND MOTION,” vol. 2004, 2004.
[37] P. E.John D. Stevenson, Ph.D., Seismic analysis of safety-related nuclear structures. American Society of Civil Engineers, 2017.
[38] F.Hammelmuller andC.Zehetner, “Increasing numerical efficiency in coupled eulerian-lagrangian metal forming simulations,” Proc. 8th Int. Conf. Comput. Plast. - Fundam. Appl. COMPLAS 2015, pp. 727–733, 2015.