在講述本部分內容之前，先回顧一下關于土體密度的相關基礎知識：

上述公式分別為土體天然密度、干密度、含水率、飽和密度、浮密度以及相互之間的關系。

巖土工程中通常認為地下水位以下土體處于飽和狀態，即涉及到飽和滲流的問題，典型的就是飽和土的固結問題?？讐和ǔＳ袃煞N表示形式，即總孔壓或超孔壓。超孔壓指的是超出靜水壓力的那部分。abaqus/standard中的流體滲透/應力耦合分析中可以基于總孔壓，也可基于超孔壓進行分析。當模型的重力載荷采用gravity分布載荷進行定義時，則基于總孔壓進行分析；若模型重力通過施加體力來實現，則采用的是超孔壓。分析中關心的是載荷引起的超孔壓分布及消散，當采用總孔壓進行分析，則需通過*initial conditions，type=pore pressure來定義初始總孔壓的分布，其中材料的密度應為干密度。流體滲透/應力耦合分析問題中采用的單元具備孔壓自由度，其單元類型標識符以字母P結尾表示孔壓單元。

飽和土體的流固耦合問題通常采用瞬態分析，進而求解孔壓、沉降隨時間的變化過程。abaqus/standard采用非對稱的剛度矩陣存儲和求解方法（當采用MC模型進行分析時同樣采用該方法）。

固結計算中一般采用自動時間步長，在該分析步中通常設置UTOL（增量步中允許的孔壓變化最大值）和CETOL（若材料模型中不包含蠕變材料特性，則不需設置此項）

進行流固耦合分析過程中，需定義滲透系數，滲透系數的單位應該與分析中的單位相對應，如m/d（天），type（類型）選擇各向同性，液體重度（間隙流體比重）通常為9.8，同時還需定義孔隙比；隨著土體的固結沉降，孔隙比在深度方向上是變化的，通常為了簡化計算，孔隙比會設置為常數，當然在軟件中也可以設置沿深度變化。嚴格上來說，由于孔隙比的變化，干密度的取值也是變化的。（相關子程序中可以自行定義，這里不再贅述）

子選項中的數值通常在非飽和土體中進行設置

在進行流固耦合分析中，通常需另外設置排水邊界，即在土體表面設置零孔壓面，區域的選擇是結點（node)集合而不是面集合（surf。在非飽和土體中排水面的設置與此不同。

巖土工程分析中，Geostatic分析步通常都是巖土分析的第一步，在該步中，土體受到重力作用，所設置的初始應力應和重力相平衡，且不產生任何位移。然而，對于復雜的模型，模型中存在過大的非線性，可能造成Geostatic步的迭代失敗，需要調整初始應力狀態，因此，精確的設置初始應力并不容易。在這里，簡單講述一種設置初始應力的方法。通常認為土體中的初始水平有效應力只和z坐標有關和水平位置無關，因此在分析中，水平有效應力由豎向有效應力乘以水平土壓力系數得到。相關的關鍵字語句為：

*initial conditions ,type=stress,geostatic;

關鍵字stress表明定義的初始應力，而關鍵字geostatic表明定義隨深度線性變化的初始應力，相應的數據行依次為單元號或單元集合名稱，起點應力值，對應坐標，終點應力值，對應坐標，x方向的水平土壓力系數，y方向的水平土壓力系數。cae中設置如圖所示：