三维边界元素法在异向性裂缝岩石之应用.pdf

2006 岩盤工程研討會論文集，台南 2006 Taiwan Rock Engineering Symposium, Tainan June 26-27 三維邊界元素法在異向性裂縫岩石之應用 1 1 1,2 陳家豪 陳昭旭 楊振成 1 國立成功大學 資源工程學系 2經濟部工業局新營工業區服務中心 摘要 本文主要提出一套適用於處理三維異向性裂縫岩石之線彈性破壞力學問題的相關分析。 採用三維邊界元素法(3D-BEM)結合異向性彈性理論來處理三維異向性裂縫岩石試體的邊界 及內部點位移及應力分析。並運用裂縫開口位移(crack opening displacement, COD)準則來求取 裂縫前緣(crack front)三種常見破壞模式的應力強度因子(stress intensity factors, SIFs)K 、K 及 I II KIII ，此為控制裂縫延伸與破壞的重要參數。於研究中亦以具有明顯層理的大理岩為試驗材 料，因此假設其為橫向等向性岩石，以此種岩石之彈性參數來探求大理岩在受到外加載重及 內部徹體力作用下之數值分析應用。 關鍵字：異向性裂縫岩石、三維邊界元素法、裂縫前緣、應力強度因子。 一、前言 在1958年Irwin提出了能量與應力強度K之關係，也就是當材料受到外來載重造成裂縫尖 端的應力分佈達到臨界狀況後，使得裂縫發生延伸破壞時，可用臨界應力強度因子(critical stress intensity factor) Kc 來描述此臨界狀況之材料性質。因此，破裂韌度(fracture toughness)即 被定義為應力強度因子(stress intensity factor, SIF)的臨界值。所以，如何決定裂縫尖端附近之 應力強度因子，已成為研究線彈性破壞力學的學者之重要的目標，以致目前已研擬出許多方 法來決定應力強度因子。然而，為了能夠探知裂縫延伸的破壞機制，所以在試樣上會先給予 預裂缺口，並在試驗中觀察此缺口之生長過程 ，藉以推斷實際情況之裂縫發展，一般引用三 個應力強度因子K 、K 及K ，來表示模式I(開裂 ) 、II(滑裂)及III(撕裂) 等裂縫受力延伸的三種 I II III 破壞模式。 在察看早期學者於岩石破壞力學的研究中，為了探求各種破壞模式而苦心研究各種試驗 方法，絕大部分學者皆在求取模態I的應力強度因子 KI ，再者即為模態II及混合模態 I-II ， 鮮少有模態 III的發展，這是因為前人皆是採用二維的分析方法來處理問題，由於受到維度的 - 61 - 陳家豪、陳昭旭、楊振成 三維邊界元素法在異向性裂縫岩石之應用 61-72 限制，而使得在處理力學問題時不能自由的展現所求。有鑑於此，本研究遂開發三維的分析 處理方法來改善既有之缺憾。然而，在面對一些複雜的環境條件或材料性質，解析解並不易 求得，因而必須藉助數值的方法來近似求解，一般常用的方法如有限差分法(finite difference method, FDM) 、有限元素法(finite element method, FEM)及邊界元素法 (boundary element method, BEM