從建立應(yīng)力平衡方程、水連續(xù)性方程和彈塑性矩陣入手,使用Galerkin方法,將各控制方程分別在空間域和時間域進行離散,初步開發(fā)出了一個用于分析飽和-非飽和孔隙介質(zhì)中水-應(yīng)力耦合彈塑性問題的二維有限元程序.通過引入一組特定的與非飽和狀態(tài)有關(guān)的計算公式,對一個假定的非飽和土體中水-應(yīng)力耦合問題進行了數(shù)值計算,考察了不同時間土體中的位移、孔隙水壓力、有效主應(yīng)力、流速和塑性區(qū)的分布,定性驗證了該程序的正確性.