以数值分析法求解版桩下游抗隆起之安全因素.docVIP

PAGE 1 PAGE 9 以數值分析法求解版樁下游抵抗隆起之安全因素 楊全成 徐佳仕 正修技術學院土木工程系 摘 要 擋土版樁土層之滲流量及版樁下游抵抗隆起之穩定分析為大地工程中相當重要課題之一，本文首先利用有限差分法，設計了一套交談式電腦程式(以下以APAS稱之)，可求解擋土版樁滲流引致之各節點過剩水頭。本文另以Excel應用軟體來計算引致隆起或砂湧區塊之平均過剩水頭以求解安全因素。 本文之APAS電腦程式可以依序輸入：(1)土層總厚度、(2)數值分析時土層每層分割厚度、(3)上升水流或下降水流、(4)滲流區域之上下游水頭差、(5)上游邊界座標及(6)下游邊界座標等數據。經反覆運算後，即可得到各節點之過剩水頭值，如此就可以進一步根據數值來繪製等勢能線及流線，並計算滲流量及抵抗隆起之安全因素。 本APAS電腦程式，可更進一步提供大地工程師設計安全之版樁埋土深度，以抵抗隆起之發生。APAS應用到某一版樁工程案例，獲得如下之結果：版樁入土深度越深，抵抗隆起之安全因素越高；上下游總水頭差越大，抵抗隆起之安全因素越低，且當土層厚10m，上下游水頭差5m，需求之安全因素1.5以上時，最低之版樁埋土深度為2.7m。 目前電腦速度已達Pentium級以上，以疊代法分析求解版樁滲流區各點之過剩水頭，所需之時間僅數秒至數分鐘即可達成；本APAS電腦程式實為求解版樁滲流區各點之過剩水頭快速有效之工具。未來如進一步設計出直接求解版樁下游抵抗隆起之安全因素，及自動繪出流線網並計算滲流量之電腦程式，以解決工程之問題，是吾人甚為期盼的。 關鍵詞：數值分析法，等勢能線，隆起，滲流量。 一、前 言 在大地工程學門上諸多問題可藉由電腦作為分析之工具，電腦可將龐大之計算程序，利用有限元素法等數值分析程式以求得其解[3]。這些應用在作邊坡穩定分析及土壤互制等方面用得甚多。更重要的是，電腦可用來求得不適宜用人力求其數值解的偏微分方程式解。困難基礎的地層[2]，複雜的壓密理論分析及滲透分析問題[4,5,6,7,8,9,10]，可歸屬於此類。偏微分方程式問題可採用有限差分法（finite difference method）[8]，或採用Taylor級數展開[4]來加以求得數值解。欲對此種方法的瞭解，故本文首先對滲流問題之有限差分法所使用公式之推導加以介紹，隨後再應用該有限差分公式，配合工程情況、土層條件與邊界條件來做對談式之電腦程式設計，最後並對電腦程式之應用結果作一介紹與討論。 土壤之滲流會引致版樁下游隆起之發生，而抵抗隆起安全因素之穩定分析為大地工程中相當重要之課題，本文以Gauss-Seidel法來設計一套求解擋土版樁滲流水頭損失之電腦程式。經由數值分析所得各節點之勢能值，可提供進一步之等勢能線之繪製，而相對之流線亦可繪製而得。此外根據Terzaghi分析版樁下游抵抗隆起之安全因素規定，本文提出一數值分析法來求解。經由一實例分析顯示：版樁入土深度越深，抵抗隆起之安全因素越高；上下游總水頭差越大，抵抗隆起之安全因素越低。 二、滲流理論之有限差分解析法 二度空間之土壤滲流情況，其水頭分佈可以用Laplace二次偏微分方程式表示之：[4] （公式1） 偏微分方程式的數值解若干大地技術問題非常適合用有限差分法求解。例如: 滲流分析、飽和粘土的壓密，以及用波動程式分析打樁情形。 有限差分法之近似解：[8] 用有限差分法之形式，來代表一般微分方程式或偏微分方程式，其微分的意義係用有限差分近似值來表示。一個函數y＝f(x)位於xi的一次微分近似值為： 可寫成 中點差分 （公式2） 增量Δx，代表沿x軸，或位在xi－1，xi，xi＋1等各點間的有限距離。註腳字母表示與指定xi值相對應的yi值。註腳字按x正值方向依次遞增。 圖一 位在xi之一次微分的中點差分近似值 圖一表示公式2的一次微分的有限差分近值情形。由圖一可以看出，位於xi點，函數y=f(x)的微分，即為位在xi點處曲線的斜率，可就曲線上兩點間的斜率作為近似值。此兩點分別用有限距離位於xi的兩側。公式2代表一次微分的中點差分近似值。xi的一次微分，亦可只考慮在xi前方的一點，或在xi後方的一點，求其近似值。此等近似值分別為一次微分的向前差分與向後差分近似值。用代數式表示時則為： 向前差分 （公式3） 向後差分 （公式4） 公式3與公式4跟公式2比起，對xi的微分的近似值，其準確性略