現在位置首頁 > 博碩士論文 > 詳目
論文中文名稱:分離元素法山崩模型建置及震滑機制行為初探 [以論文名稱查詢館藏系統]
論文英文名稱:Landslide model set up and simulation using discrete element method [以論文名稱查詢館藏系統]
院校名稱:臺北科技大學
學院名稱:工程學院
系所名稱:土木與防災研究所
畢業學年度:99
出版年度:100
中文姓名:翁晨桓
英文姓名:Chan-Huan Weng
研究生學號:98428036
學位類別:碩士
語文別:中文
口試日期:2011-07-18
論文頁數:109
指導教授中文名:張國楨
口試委員中文名:林銘郎;王泰典
中文關鍵詞:分離元素法數值地形模型草嶺山崩集集地震
英文關鍵詞:Discrete element methodDigital terrain modelTsaoling landslideChi-Chi earthquake
論文中文摘要:本研究以1999年集集地震所引發之草嶺山崩事件為例,運用數值地形模型(DTM)資訊,結合顆粒流分離元素之PFC 3D程式分析山崩之行為及機制。為決定山崩塊體及山崩影響範圍,分別運用災前、災後DTM來分析,並用以規範崩塌區之上、下部地形面,定義其間所夾空間即為崩塌之區塊。首先,本研究嘗試數種不同方法來建置崩塌區塊體及模型顆粒叢集等數值模擬之準備,包含NO SHADOW指令法、霣降法、三角柱法、重力排除法、顆粒汰選法等,比較各方法之優缺點及模型可行性等。本研究結果指出顆粒汰選法為最佳方法,其方法為依模擬岩石試體製作流程預先製作一含蓋整個崩塌塊體之顆粒匯集區,再將顆粒予以排除或汰選,配合施予顆粒元素間鍵結作用,進而建立一接近實際外幾何之形態、厚度之崩塌區塊,滑動面亦由DTM資料來建立,符合真實地形特徵之地表模型。
本研究考慮草嶺強震站之地震加速度資料,將之轉換為給定模型之地表地動速度,並考慮岩層中存在孔隙水壓所造成之影響,另考量岩層因受震產生高速滑移,並引致岩層摩擦係數弱化之影響。考量以上因素後,最後以簡化方式一併列入模擬之不同參數及相關模擬條件。
由模擬結果數據資料建置之DTM與真實DTM比對顯示岩層弱化後摩擦係數為0.087時可得最接近真實之堆積型態,其堆積區高程差平均值為-3.2m,樣本標準差為17.5。山崩發生之最大滑動距離為3662 m,最大落距高度為632 m,最高速度為87.2 m/s,本速度峰值之時間序列和地震訊號所得結果吻合。
研究成果顯示以三維離散元素顆粒流數值方法可模擬山崩發生之完整動態歷程,相較於一般常用之極限平衡法及有限元素法,除可求得破壞面及應力應變得知穩定性外,該方法還可模擬山崩發生之大變形行為,包含裂縫、運動軌跡、速度變化、滑動距離,乃至三維的堆積型態,進一步掌握災害發生之影響範圍。研究成果指示了本模擬之成效,以及震滑相關現象之推論。
論文英文摘要:This research takes Tsaoling landslide caused by Chi Chi Earthquake as an example, adopts from DTM data and also integrates the PFC 3D program of Particle flow discrete element method to analyze the sling mechanism and the landslide behaviors.
To determine the affected area of landslide and to constraint the slid block, we use pre-disaster DTM and post-event DTM to build two topographic surfaces (upper topographic surface and lower topographic surface), thus the slid block between upper topographic surface and lower topographic surface could be defined.
First of all, in this study several different methods are developed and tried to build the slidblock and the particle clusters to confine the slidblock, including NO SHADOW command Method, Pluviation Method, Triangular Prism Method, Gravity Elimination Method、Particle Detection Method, etc. The above-mentioned methods are compared according to the availability and feasibility. The Particle Detection Method is composed of several procedures: Firstly, prepare one dense and cubic particle clusters and then transpose to the mass center of the slid block; Secondary, the particle cluster should be rescaled in order to cover the entire slid block; Third, detect the particlesif the position is situated within the slid blocks confined by the upper and lower surface of the DTM before and after landslide, otherwise, the particle is been eliminated from the cluster. After the procedure, all there maining particles are bonded according to the material properties and served as the model for further landslide simulation. Therefore an almost realistic geometry concerning the thickness and shapes of slid block could be established, and the sliding surface also could be built from DTM data to match the real terrain features of topography.
This study considers earthquake ground excitation by the strong motion acceleration data gather from the CHY080 seismostation situated just nearby the sliding area. Furthermore the existence of pore water pressure within the rock mass and the effect of friction coefficient attenuation caused by rock mass shaking and sliding are also integrated in this study.
Compare the topography results from the simulation data with the real DTM, it says that the gained accumulation patterns of the simulation result could be mostly similar to real situation (DTM after landslide) at the friction coefficient of 0.087. The average difference between the DTM and PFC results is -3.2m and the standard deviation is about 17.5. When landslides occur, the largest sliding distance is about 3662 m, the maximum drop height is about 632m and the highest sliding speed is 87.2 m/s. The timing of maximum sling velocity coincides with seismograph of CHY080.
The results shows that use the numerical method of three-dimensional discrete element particle flow can simulate the complete landslide process dynamically. Comparing with the limit equilibrium method and finite element method which are the most common to be used, the discrete element particle flow provides not only the states of instability but also the dynamical process of sliding. The result from this study states the effect of this simulation and states the inference of the related earthquake-sliding phenomena as well.
論文目次:誌謝 I
摘 要 III
ABSTACT V
目 錄 VII
表目錄 IX
圖目錄 X
第一章 緒論 1
1.1 研究動機 1
1.2 研究目的 2
1.3 研究流程 3
1.4 論文內容 5
第二章 文獻回顧 6
2.1草嶺崩塌區之活動與歷史概述 6
2.2草嶺之地形與地質 9
2.2.1 草嶺之地形 9
2.2.2 草嶺之地層 10
2.2.3 草嶺之地質構造 13
2.3 草嶺山崩現象與促崩因素 15
2.3.1 草嶺之山崩現象 15
2.3.2 草嶺之山崩促崩因素 16
2.4 分離元素法概述 17
2.5 分離元素法應用於大地工程之相關研究 18
第三章 數值模擬分析方法 20
3.1 地震作用模擬方法 20
3.2 岩層孔隙水壓作用模擬方法 26
3.3 PFC3D(PARTICLE FLOW CODE IN 3 DIMENSIONS)概述 29
3.4 PFC3D之假設 29
3.5 PFC3D之運算邏輯 30
3.6 PFC3D之微觀參數 31
3.6.1 顆粒元素參數 31
3.6.2 鍵結力參數 33
3.6.3 牆面元素參數 36
3.6.4 微觀參數檢核之步驟 37
第四章 模型建置方法 38
4.1區域選取與地表模型建置 38
4.2 NO SHADOW指令法 41
4.3 霣降法 43
4.4 三角柱法 44
4.5 重力排除法 45
4.6 顆粒汰選法 51
4.7 各模型建置方法之比較 54
第五章 山崩模擬結果分析與討論 55
5.1 模擬條件與假設 55
5.2 模型參數選定 56
5.3 微、巨觀參數校準 64
5.4 堆積型態及剖面比對與滑動歷程展示 72
5.4 監測顆粒佈設 82
5.5 地震作用裂縫監測與應力歷程 83
5.6 山崩滑動軌跡分析 87
5.7 山崩滑動速度與距離分析及地震訊號比對 91
5.8 地震力與孔隙水壓力之影響初探 99
第六章 結論與建議 103
6.1 結論 103
6.2 建議 105
參考文獻 106
論文參考文獻:【1】 中興工程顧問股份有限公司,「草嶺崩坍之處理對策評估」,經濟部水利處,1990。
【2】 何維信,航空攝影測量學,大中國圖書公司,台北,1995。
【3】 呂榮華,「以單顆粒勁度探討卵礫石層之抗剪強度」,碩士論文,國立中興大學土木工程研究所,台中,1999。
【4】 李宏輝,「砂岩力學行為之微觀機制-以個別元素法探討」,博士論文,國立台灣大學土木工程學系,台北,2008。
【5】 李後彥,「數值地形模型於山崩及後續地形演化之應用-以草嶺山崩為例」,碩士論文,國立台北科技大學土木與防災研究所,台北,2009。
【6】 李憲忠,「九二一集集地震三維震源過程與震波傳遞分析」,博士論文,國立中央大學地球物理研究所,桃園,2003。
【7】 李錫堤、洪如江、林銘郎、蔡龍珆,「草嶺崩塌地工程地質調查與穩定性評估」,中興工程顧問社,專題研究報告,1993。
【8】 李錫堤、林銘郎、吳禮浩、鄭俊昇等人,「草嶺大崩山區的地質調查及歷次大崩山滑動面的決定」,1994 年岩盤工程研討會,1994。
【9】 李錫堤、洪如江、林銘郎、蔡龍珆等人,「草嶺崩塌地工程地質調查與穏定性評估-專題研究報告」,經濟部水資局統一規劃委員會,1993。
【10】 林立欣,「離散元素模擬之平行策略」,碩士論文,國立台灣大學土木工程學研究所,台北,2005。
【11】 林昇鴻,「卵形顆粒介質運動之數值模擬」,碩士論文,國立中央大學土木工程研究所,桃園,2006。
【12】 林郁修,「分離元素法於岩石貫切破壞試驗之模擬分析」,碩士論文,國立台北科技大學土木與防災研究所,台北,2007。
【13】 洪如江、李錫堤、林美聆、林銘郎、鄭富書、陳正與,「天塹可以飛渡、崩山足以斷流(草嶺順向坡滑動) 」,地工技術第77 期,2000。
【14】 唐昭榮,「臺灣遽變式山崩傳送與堆積之顆粒流離散元素模擬」,博士論文,國立臺灣大學理學院地質科學研究所,台北,2010。
【15】 高贈智,「集集地震引致台灣中部山區邊坡臨界滑移量之分析」,碩士論文,國立台灣大學土木工程學研究所,台北,2004。
【16】 張來福,「車籠埔斷層周圍岩石力學特性之初探」,碩士論文,國立台灣大學土木工程學研究所,台北,2008。
【17】 張家偉,「應用分離元素法探討紅菜坪地區地滑演化」,碩士論文,國立台灣大學土木工程學系,台北,2007。
【18】 張家詮,「分離元素法於擬脆性岩材微觀破裂機制之初探」,碩士論文,國立台北科技大學土木與防災研究所,台北,2007。
【19】 連俊凱,「具節理強度岩體因落石作用撞擊碎裂與堆積行為」,碩士論文,國立台灣大學土木工程學系,台北,2010。
【20】 陳正旺,「車籠埔斷層周圍岩石力學特性之初探」,碩士論文,國立台灣大學土木工程學系,台北,2005。
【21】 陳春華,「非連續體變形分析法模擬集集地震引致草嶺邊坡崩塌模擬之研究」,碩士論文,國立成功大學土木工程研究所,台南,2008。
【22】 彭健豪,「九份二山地滑區滑動歷程與積行為之研究」,碩士論文,國立台灣大學土木工程學系,台北,2008。

【23】 黃鑑水、何信昌、劉桓吉等人,「台灣中部草嶺地區之地質與山崩」,經濟部地質調所彙刊,第二號,1983。
【24】 葉怡成,測量學-21世紀觀點,東華書局,2006,台北。
【25】 廖豐益,「以數值分析探討粒狀材料之破壞行為」,碩士論文,國立台灣工業技術學院營建工程技術研究所,台北,1994。
【26】 蕭安佑,「洪氾及山崩模擬-以金山及小林村為例」,碩士論文,國立台北科技大學土木與防災研究所,台北,2010。
【27】 簡士堯,「低凝聚力材料趾部受侵蝕崩塌滑動模式研究」,碩士論文,國立台灣大學土木工程學系,台北,2009。
【28】 Chen, R.-F., K.-J. Chang, J. Angelier, Y.-C. Chan, B. Deffontaines, C.-T. Lee, and M.-L. Lin (2006), Topographical changes revealed by high-resolution airbone LiDAR data: The 1999 Tsaoling landslide induced by the Chi-Chi earthquake, Eng. Geol., 88, 160-172.
【29】 Chen, R.-F., Y.-C. Chan, J. Angelier, J.-C. Hu, C. Huang, K.-J. Chang, and T.-Y. Shih (2005), Large earthquake-triggered landslides and mountain belt erosion: The Tsaoling case, Taiwan, Comptes Rendus Géoscience, Acad. des Sci. Paris, 337, 1164-1172.
【30】 Cundall, P. A. and O. D. L. Strack, “A Discrete Numerical Model for
【31】 Cundall, P.A. (1971), A computer model for simulating progressive large scale movements in blocky rock systems. Proc. of the Sym, Int. Soc. Rock Mech., 1 No. II-8.
Granular Assemblies,” Geotechnique, Vol. 29, No. 1, pp. 47-65 (1979).
【32】 Hung, J.-J. (2000), Chi-Chi earthquake induced landslides in Taiwan, Earthquake Eng. Eng. Seismol., 2 (2), 25-33.
【33】 Itasca (2006), User’s Manual, PFC 3D, Itasca Consulting Group, Inc., Minneapolis, Minnesota.
【34】 Kuo, C.-Y., Y.-C. Tai, F. Bouchut, A. Mangeney, M. Pelanti, R.-F. Chen, and K.-J. Chang (2009), Simulation of Tsaoling landslide, Taiwan, based on Saint Venant equations over general topography, Eng. Geol., 104, 181-189.
【35】 Potyondy, D.O., and P.A. Cundall (2004), A bonded-particle model for rock, Int. J. Rock Mech. Min. Sci., 41, 1239-1364.
【36】 Togo, T., T. Shimamoto, S.L. Ma, and T. Hirose (2009), High-velocity friction of faults: A review and implication for landslide studies. The Next Generation of Research on Earthquske-induces landslides, An Int. Conf. Comm. 10th Anni. of the Chi-Chi Earthquake, 205-216.
【37】 USGS (2004), Landslide Types and Processes, Complied by Highland, L., USGS Fact Sheet 2004-3072.
論文全文使用權限:不同意授權