2009年8月7日莫拉克颱風侵襲台灣，短時間內之大量降雨引發高雄縣甲仙鄉後方獻肚山的崩塌，大量崩塌的岩塊將山坡下方之小林村掩埋並造成許多人員傷亡。面對此大型之崩塌，使用數值分析軟體進行模擬並與現地狀況比對為普遍之研究方法，但對於尚未發生崩塌之南投縣九份二山殘坡，在無實際結果比對下，如何利用數值分析軟體訂定出崩塌後之影響範圍，並提前對可能受災範圍進行防護應為目前較為重要之課題。
Distinct Element Method (DEM)為解決非連續體分析問題所發展出之數值分析方法，本研究使用以其為基礎所發展出之三維元素法分析軟體3DEC對大型邊坡進行模擬，探討獻肚山與九份二山地滑區在重力滑落下的崩塌歷程與崩塌範圍，並以不同摩擦角對邊坡進行敏感度分析。由獻肚山崩塌之模擬結果顯示，將590高地頂部及其東側之刮削情況納入考量並使用低摩擦角後所得出之模擬崩塌範圍會較接近真實情況，而九份二山殘坡崩塌之模擬結果顯示，當山崩發生時，崩落之岩塊會對龍南道路造成影響，使龍南道路上之行車安全受到危害，此時模擬出之崩塌範圍可做為當地安全評估上之參考，使災害發生時能減少對人民生命安全上的危害。

On 9 August 2009, the heavy rainfall brought by Typhoon Morakot induced the Hsien-du-shan rock avalanche. The falling rock buried the Hsaiolin village and caused lots of damage. Facing this kind of large-scale landslide, using numerical analysis to simulate the range of landslide and compare to the in-situ condition is the common research method. But for the residual slope of Chiu-fen-erh-shan, there is no result of slope failure can be compared. Under this situation, how to use the numerical method to simulate the range of landslide and then decrease the loss of disaster is an important issue for researchers currently.
Discrete element method (DEM) is a numerical method developed to solve the non-continuum analysis problem. In this thesis, the program 3DEC was used based on the theory of DEM to simulate the Hsien-du-shan and Chiu-fen-erh-shan landslide. And use different frictional coefficient to simulate the range of landslide and discuss the slope failure mechanism.
The result shows that considering the collapse of 590 height and using lower frictional coefficient to simulate the Hsien-du-shan landslide can obtain a better result. And the result of Chiu-fen-erh-shan residual slope landslide simulation indicates the Lung-nan path will be affected by the falling rock when rock avalanche occurs. At this time, the range of Chiu-fen-erh-shan residual slope landslide simulated by 3DEC can provide a safety assessment to local residents for decreasing the threat of disaster.

[1]中興顧問社，防災科技研究中心.(2015). 大規模崩塌地調查與數值模擬. http://dptrc.sinotech.org.tw/chinese/01_group/07_print.php?group=2&pid=265.

[2]王卫华, 李夕兵, 胡盛斌. (2005). 模型参数对3DEC动态建模的影响. 岩石力学与工程学报. 24, 4791-4797.

[3]巨能攀, 赵建军, 黄润秋, 邓辉, 王彦东. (2009). 基于 3DEC 的边坡块体稳定性分析. 辽宁工程技术大学学报: 自然科学版, 28(6), 925-928.

[4]吳秋雅. (2007). 九份二山山崩前後構造地形特徵及坡體破壞機制探討. 國立臺灣大學地質科學研究所碩士學位論文, 1-117.

[5]李錫堤, 董家鈞, 林銘郎. (2009). 小林村災變之地質背景探討. 地工技術. 122, 87-94.

[6]肖燃, 郑达, 巨能攀. (2014). 崩塌誘發其下部滑坡變形的分析研究. 工程地質學報. 22(3), 413-419.

[7]周正义, 曹平, 林杭. (2006). 3DEC 应用中节理岩体力学参数的选取. 西部探矿工程, 18(7), 163-165.

[8]唐昭榮, 胡植慶, 羅佳明, 林銘郎. (2009). 遽變式山崩之 PFC3D 模擬初探—以草嶺與小林村為例. 地工技術, (122), 143-152.

[9]陳建宏. (2011). 小林村獻肚山深層崩塌地質構造及地質材料特性之研究. 國立成功大學土木工程學系碩士學位論文, 1-142.

[10]鄭文棠, 徐衛亞, 童富果, 石安池. (2007). 複雜邊坡三維地質可視化和數值模型構建. 岩石力學與工程學報, 26(8), 1634-1644.

[11]魏群. (1991). 離散單元法的基本原理數值方法及程序. 科學出版社.
[12]闕衍智. (2013). 以三維分離元素法模擬大型邊坡受地震力作用下影響之研究. 成功大學土木工程學系碩士學位論文, 1-123.

[13]羅佳明, 林彥享, 魏倫瑋, 黃韋凱, 鄭添耀, 林銘郎, 李錫堤. (2012). 遽變式山崩調查與動態模擬之研究. 中國土木水利工程學會期刊, 24(3), 231-243.

[14]Cundall, P. A., & Strack, O. D. (1979). A discrete numerical model for granular assemblies. Geotechnique, 29(1), 47-65.

[15]Huang, C. S., Chen, M. M., & Hsu, M. I. (2002). A preliminary report on the Chiufenershan landslide triggered by the 921 Chichi earthquake in Nantou, central Taiwan. Terrestrial, Atmospheric and Oceanic Sciences, 13(3), 387-395.

[16]Itasca Consulting Group. (2007). 3DEC - Theory and background. Itasca Consulting Group.

[17]Itasca Consulting Group. (2007). 3DEC - Command Reference. Itasca Consulting Group.

[18]Itasca Consulting Group. (2007). 3DEC - User’s Guide. Itasca Consulting Group.

[19]Jing, L., & Stephansson, O. (2007). Fundamentals of Discrete Element Methods for Rock Engineering: Theory and Applications. Elsevier, 85.

[20]Kulatilake, P. H. S. W., Ucpirti, H., Wang, S., Radberg, G., & Stephansson, O. (1992). Use of the distinct element method to perform stress analysis in rock with non-persistent joints and to study the effect of joint geometry parameters on the strength and deformability of rock masses. Rock Mechanics and Rock Engineering, 25(4), 253-274.

[21]Kuo, C. Y., Tai, Y. C., Chen, C. C., Chang, K. J., Siau, A. Y., Dong, J. J., Lee, C. T. (2011). The landslide stage of the Hsiaolin catastrophe: Simulation and validation. Journal of Geophysical Research: Earth Surface (2003–2012), 116(F4).

[22]Lo, C. M., Lin, M. L., Tang, C. L., & Hu, J. C. (2011). A kinematic model of the Hsiaolin landslide calibrated to the morphology of the landslide deposit. Engineering Geology, 123(1), 22-39.

[23]Scheidegger, A. E. (1973). On the prediction of the reach and velocity of catastrophic landslides. Rock mechanics, 5(4), 231-236.

[24]Shou, K. J., & Wang, C. F. (2003). Analysis of the Chiufengershan landslide triggered by the 1999 Chi-Chi earthquake in Taiwan. Engineering Geology, 68(3), 237-250.

[25]Tsou, C. Y., Feng, Z. Y., Chigira, M. (2011). Catastrophic landslide induced by typhoon Morakot, Shiaolin, Taiwan. Geomorphology, 127(3), 166-178.

[26]Wu, J. H., Chen, J. H., Lu, C. W. (2013). Investigation of the Hsien-du-Shan rock avalanche caused by typhoon Morakot in 2009 at Kaohsiung county, Taiwan. International Journal of Rock Mechanics and Mining Sciences, 60, 148-159.

[27]Wu, J. H., Lin, J. S., & Chen, C. S. (2009). Dynamic discrete analysis of an earthquake-induced large-scale landslide. International Journal of Rock Mechanics and Mining Sciences, 46(2), 397-407.

[28]Yano K, Shimamoto T, Oohashi K, Dong JJ, Lee CT(2009) Ultra-low friction of shale and clayey fault gouge at high velocities: implication for Jiufenger-shan landslide induced by 1999 Chi–Chi earthquake.The next generation of research on earthquake-induced landslides: an international conference incommemoration of 10th anniversary of the Chi–Chi Earthquake, 402–406.