現在位置首頁 > 博碩士論文 > 詳目
  • 同意授權
論文中文名稱:應用分離元素法探討順向坡失穩歷程-以國道三號崩塌事件為例 [以論文名稱查詢館藏系統]
論文英文名稱:Dip-slope failure assessment based on Discrete numerical simulation - example from Formosa highway [以論文名稱查詢館藏系統]
院校名稱:臺北科技大學
學院名稱:工程學院
系所名稱:土木與防災研究所
畢業學年度:100
出版年度:101
中文姓名:張賀翔
英文姓名:Ho-Shyang Chang
研究生學號:99428036
學位類別:碩士
語文別:中文
口試日期:2012-07-20
論文頁數:99
指導教授中文名:張國楨
口試委員中文名:王泰典;葉恩肇
中文關鍵詞:分離元素法數值地形模型國道三號順向坡
英文關鍵詞:discrete element methoddigital terrain modelFormosa highwaydip-slope failure
論文中文摘要:鑒於台灣板塊運動頻繁且劇烈,造成山多平原少,對於交通的拓展十分不便,迫於無奈,經常需要將道路建於邊坡上或貫穿邊坡,而當中很容易遭遇順向坡,因此在施工中,順向坡坡腳遭破壞,邊坡就失去原本穩定性,如果新建物無法負荷原本坡腳之載重,往往會發生重大災害。民國九十九年四月二十五日,國道三號3.1公里處發生地滑事件,造成四人罹難。對於地滑地點在設計規畫階段皆符合規範設計,但仍發生重大的公安災害令人十分震驚。
本研究運用數值地形模型(DTM)具有詳細三維座標的特點,並結合顆粒流離散元素分析軟體PFC^3D,基於其特性足以模擬大變形進行模擬,並分析地滑的行為及機制。首先,本研究利用災前、災後的航空照片分別建置成數值地形模型(DTM) ,由於本例不同以往之處為地滑塊體與滑動面並未分離,導致滑動面難以取得,因此,本研究利用上部露出之岩層滑動面進行內差,進而取得完整滑動面。之後,利用災前DTM與有完整滑動面之災後DTM來分離出滑動塊體,經計算後,地滑塊體約為十五萬立方公尺之土石。在PFC^3D建置地滑塊體方面,首先,生成由顆粒所匯集而成之矩形塊體,等向應力及孔隙率皆於此步驟先行設置完成,並配合本團隊所研發之程式,將災前、災後DTM所夾空間完全涵蓋矩形塊體,並由程式進行排除即得地滑塊體之顆粒。此外,程式可將矩形塊體任意縮放比例,達到套疊任一崩塌區之快速建模目的。滑動模型由災前DTM所轉換之地形面與顆粒所代表之地滑區塊所組成。
本研究進一步由地質圖判斷,將地滑塊體分為因含水弱化的砂頁岩互層及上部較完整的砂岩兩部分區塊,並在前端另設置一區顆粒來模擬地錨作用。考慮以上因素後,將參數一併輸入,完成模型準備。在模型變數上,本研究選擇滑動面弱化之摩擦角由10度弱化至2度,地錨剩餘效率70%、50%以及30%,將滑動面弱化程度及地錨折減率,以組合之方式設置於地滑模型來模擬真實情況下滑動面弱化情況,模擬結果為當滑動面弱化至4度、地錨剩餘效率為30%時與真實情況最為符合現地堆積狀況。山崩歷時約17秒,到達國道三號邊緣時間為4.36秒,掩埋南下路段時間為6.96秒,掩埋北上路段時間為10.39秒,崩落過程所產生之最高速度為15.34m/s,最遠水平滑動距離為102.84 m,最大落距為29.43 m。
使用數值地形模型(DTM)搭配顆粒流離散元素分析軟體PFC^3D其優點為將真實地形面完整呈現於模型上,以更符合災害發生實際狀況,並可模擬地滑發生之完整動態歷程。相較於一般常用之有限元素分析軟體,僅可求得破壞當時應力應變得知穩定性外,本方法還可將破壞後塊體持續發生之大變形行為,如運動軌跡、速度變化、滑動距離,以及停止運動後堆積型態,進一步掌握災害發生之影響範圍。
論文英文摘要:It is because that frequent and intense of the plate movement in Taiwan make it so inconvenient to construct transportation system, and thus transportation is restricted to develop on the dip-slope which would impair the toe of slope leading to instability during construction. Accordingly, if the new-built construction cannot bear the weight of toe of slope, tremendous disaster will occur.Four people passed away on April 25th 2010 as a result of a landslide in Formosa Freeway. It was astonishing that under legitimate regulation the design of the construction still caused great casualty of public security.
This research applied the Digital Terrain Model (DTM) featuring specific three dimensional coordinate, combines Particle Flow Code (PFC3D) which could simulate great deformation and analyze the landslide formation. First of all, builds up DTM based on the aerial photos before and after landslide. The different point from other cases is that there is no separation between slide block and slide slope result in difficult to acquire its slide slope. Using slide slop on upper part of exposed strata to obtain completed slide slop through interpolation method.Second, via calculation, the slide block is approximately 150 thousand cubic meters. In the PFC3D part, the topographic surface, based on pre-disaster and post-event DTM, serve as sliding surface of the numerical simulations models. To begin with, particle clusters production that the isotropic stress and porosity should take into account. To ensure the production range should cover the entire slid mass from the source area. Besides that, to transpose the particle clusters to the numerical model where the source is. Moreover, keep the particles which the positions are situated within the slide block confined by the upper and lower surface of the DTM before and after landslide, or eliminate from the clusters.
According to the geological map, the particles of slide block can be divided into two parts: 1) interbedded sandstone and shale which soften by water 2) upper part of sandstone. Furthermore, set up a layer of particles to simulate ground anchor. In the end, the model is ready for simulation. The advantages of DTM collocate PFC3d are that real terrain can be represented on the model, and can be simulated the complete landslide process dynamically.
Comparing with other finite element analysis software that only provides state of instability but by using this method we can know the dynamical process of sliding include trajectory, velocity change, sliding distance and also accumulation patterns after stopping and know the effect areas from the disaster event.
論文目次:摘 要 i
ABSTRCT iii
誌謝 v
目錄 vi
表目錄 viii
圖目錄 ix
第一章 緒論 1
1.1 研究動機與目的 1
1.2 研究方法與範圍 2
1.3 論文架構與內容 3
第二章 文獻回顧 5
2.1 國道三號之歷史概述與災害 5
2.1.1 國道三號之歷史 5
2.1.2 國道三號之災害 6
2.1.3 國道三號之歷史災害紀錄 8
2.2 國道三號之地形與地質 10
2.2.1 國道三號附近區域之地質 10
2.2.2 國道三號附近地區地形特徵 13
2.2.3 國道三號之山崩因素探討 13
2.3 離散元素法概述 15
2.4 PFC應用於大地工程之相關研究 16
第三章 數值模擬分析方法 17
3.1 PFC3D之概述 17
3.2 PFC3D之基本假設 18
3.3 PFC3D之運算邏輯 18
3.4 PFC3D之微觀參數 19
3.4.1 顆粒元素參數 19
3.4.2鍵結力參數 21
3.4.3 牆面元素參數 24
3.4.4 微觀參數檢核之步驟 25
第四章 模型建置方法 26
4.1 數值地形模型建置與邊界設定 26
4.2 滑動面前處理 28
4.2.1 應用數值地形模型生成正射影像 28
4.2.2 滑動面決定 29
4.3 PFC3D地表模型建置 31
4.4 PFC3D模型之崩塌塊體 34
4.4.1顆粒匯集之矩形塊體 35
4.4.2矩形塊體之應用 39
4.4.3顆粒排除原理 40
4.5 岩層模擬法 43
4.6 地錨模擬方法 45
4.6.1 模擬地錨 46
4.6.2 模擬地錨位置 49
4.7完整模型建置 49
第五章 模擬結果分析與討論 53
5.1模擬條件與假設 53
5.2 模型微觀參數設定 54
5.3 微、巨觀參數校正 61
5.4 弱化分析及堆積型態比對分析 65
5.4.1參數弱化分析 65
5.4.2堆積型態比對 66
5.5 裂縫分析 72
5.6 監測顆粒選定 74
5.6.1 滑動軌跡分析 75
5.6.2滑動速度、距離與落距分析 79
5.7 災害歷時分析 86
第六章 結論與建議 89
6.1 結論 89
6.2 建議 90
參考文獻 91
附錄A 94
論文參考文獻:[1] 何維信,航空攝影測量學,台北,大中國圖書公司,1995。
[2] 呂榮華,以單顆粒勁度探討卵礫石層之抗剪強度,碩士論文,國立中興大學土木工程研究所,台中,1999。
[3] 李宏輝,砂岩力學行為之微觀機制-以個別元素法探討,博士論文,國立台灣大學土木工程學系,台北,2008。
[4] 林立欣,離散元素模擬之平行策略,碩士論文,國立台灣大學土木工程學研究所,台北,2005。
[5] 林郁修,分離元素法於岩石貫切破壞試驗之模擬分析,碩士論文,國立台北科技大學土木與防災研究所,台北,2007。
[6] 唐昭榮,臺灣遽變式山崩傳送與堆積之顆粒流離散元素模擬,博士論文,國立臺灣大學理學院地質科學研究所,台北,2010。
[7] 張家偉,應用分離元素法探討紅菜坪地區地滑演化,碩士論文,國立台灣大學土木工程學系,台北,2007。
[8] 張家詮,分離元素法於擬脆性岩材微觀破裂機制之初探,碩士論文,國立台北科技大學土木與防災研究所,台北,2007。
[9] 廖豐益,以數值分析探討粒狀材料之破壞行為,碩士論文,國立台灣工業技術學院營建工程技術研究所,台北,1994。
[10] 粘為東,以PFC2D模擬砂土直剪實驗中之剪動帶及應用之研究,碩士論文,國立臺灣大學工學院土木工程學系,台北,2010。
[11] 江明軒,FLAC耦合PFC2D應用於雙楔刀貫切破壞之初探,碩士論文,國立台北科技大學土木與防災研究所,台北,2008。
[12] 吳曉琦,利用數值模擬探討覆瓦狀構造的力學機制,碩士論文,國立臺灣大學工學院土木工程學系,台北,2008。
[13] 翁晨桓,分離元素法山崩模型建置及震滑機制行為初探,碩士論文,國立台北科技大學土木與防災研究所,台北,2011。
[14] 羅若瑜,床面剪應力及顆粒彈跳引致岩石磨蝕之微觀數值模擬,碩士論文,國立交通大學土木工程學系,新竹,2011。
[15] 廖洪鈞、石原研而等人,國道3號3.1公里崩塌事件原因調查工作總結報告,社團法人中華民國大地工程學會,台北,2011。
[16] 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, pp.160-172.
[17] 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 Geoscience, Acad. des Sci. Paris, 337, pp.1164-1172.
[18] Cundall, P. A. and O. D. L. Strack, “A Discrete Numerical Model for Granular Assemblies, “Geotechnique, Vol. 29, No. 1, pp. 47-65 (1979).
[19] 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.
[20] Hung, J.-J. (2000), Chi-Chi earthquake induced landslides in Taiwan, Earthquake Eng. Eng. Seismol., 2 (2), pp.25-33.
[21] Itasca (2006), User’s Manual, PFC 3D, Itasca Consulting Group, Inc., Minneapolis, Minnesota.
[22] Chia-Ming Lo, Ming-Lang Lin, Chao-Lung Tang, Jyr-Ching Hu (2011),A kinematic model of the Hsiaolin landslide calibrated to the morphology of thelandslide deposit, Engineering Geology., 123,pp. 22-39.
[23] 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, pp.181-189.
[24] Potyondy, D.O., and P.A. Cundall (2004), A bonded-particle model for rock, Int. J. Rock Mech. Min. Sci., 41, pp.1239-1364.
[25] 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.
[26] USGS (2004), Landslide Types and Processes, Complied by Highland, L., USGS Fact Sheet 2004-3072.
[27] Jeoungseok Yoon (2007), Application of experimental design and optimization to PFC model calibration in uniaxial compression simulation, Int. J. Rock Mech. Min. Sci., 44, 871-889.
[28] Yuannian Wang, Fulvio Tonon (2009), Modeling Lac du Bonnet granite using a discrete element model, Int. J. Rock Mech. Min. Sci., 46, 1124-1135.
論文全文使用權限:同意授權於2017-08-22起公開