河道护坡稳定分析模式之发展与探讨.pdfVIP

第八屆海峽兩岸水利科技交流研討會，2004年 12月 河道護坡穩定分析模式之發展與探討 姜世偉 蔡東霖、楊錦釧 交通大學土木所 交通大學防災工程研究中心 摘要 本研究旨在發展一河道護坡穩定分析模式，模式乃根據多孔彈性介質理論(poro-elastic theory) 以及摩爾-庫倫(Mohr-Coulomb)破壞準則，計算飽和護坡之穩態破壞指標，控制方程式乃耦合考慮 穩態地下水流方程式與土體平衡方程式，並且可結合河道水理模式，以反應河道對護坡穩定的影 響。在數值處理上，利用有限解析法配合對角線卡氏座標系統，有限解析法乃是利用局部解析解 (local analytic solution)建構離散方程式 (discritization equations) ，較有限差分法(finite difference method)準確。此外，有限解析法更可以在卡氏座標下，簡單地處理不規則模擬邊界之計算問題， 所以較有限元素法(finite element method)容易建構。最後，將模式應用於一假想案例，模擬分析不 同條件對河道護坡穩定之影響。 一、前言 台灣河川特性多元，洪枯懸殊，沖淤頻繁，對於不同自然特性之河川及河段， 諸如河床坡度、底床質、灘瀨及基流量等自然條件，應妥慎研選適合之工法。在台灣 水文、地文條件極為特殊之情形下，第一線設計工程人員常有安全性不足之疑慮，因 而常裹足不前，致影響推展之進程。故有必要針對適合台灣河川自然條件之護岸工法 詳加分析評估，瞭解相關工法之適用條件及最佳配置，訂定一套合理且嚴謹之設計分 析模式。 一般常用於分析邊坡穩定的方法主要有極限平衡法與數值分析方法，極限平衡法 262 假設邊坡未破壞時土體為剛體且無變形，破壞發生時則為完全剪力破壞，沿破壞面均 達塑性，且破壞面上的每一點安全係數均相同，即未考慮材料之彈性變形，並以最可 能發生破壞面對應之安全係數來定義邊坡穩定性，另外未能考量土壤之張力強度，因 其應用上較為簡易，故為工程師所廣泛使用，常用的模式如 STABL 。而數值方法主 要考慮材料的應力-應變性質，其分析理論較為繁複並考慮到較多的影響因素，且由 於土體材料性質不易掌握，而限制了其在邊坡穩定分析上的實用性，常用的模式如 FLAC 。但上述兩個邊坡穩定分析模式由於並非專為河川護坡穩定分析目的所發展， 因此並無法有效考慮近岸水流與土體互制之效應。河川水位會影響地下水流之孔隙水 壓，同理地下水滲流量亦會影響河川水流，孔隙水壓改變會造成土體有效應力變化， 進而導致河岸護坡之破壞。因此河道護坡穩定與河川水流、地下水流與土體有效應力 之相互關係，乃是一值得探討的課題。 有限解析法(finite analytic method)由佛羅里達大學工學院院長陳景仁(Chen, C.J.) 教授於 1980 年所提出，已利用於計算流體動力學(computational fluid dynamics)以及 地下水質量傳輸(groundwater mass transport) 多年，獲得不錯之模擬結果。有限解析法 乃是利用局部解析解(local analytic solution)建構離散方程式 (discritization equations) ， 因此不但具有無條件穩定收斂(unconditionally stable)之特點，且其準確度可能高於有 限差分法且與有限元素法相當(Hwang, 1985 ；Chen, 1995) 。此外，有限解析法更可在 對角卡氏座標(diagonal Cartesian coordinates)系統下，處理不規則模擬邊界之計算問 題，所以有限解析法與有限差分法一樣，較有限元素法容易建構。 二、數值模式 1. 控制方程式 二維地下水流控制方程式符合達西定律(Darcy’s Law)以及土體為等向(isotropic) 之假設下可表示成