圧密解析 — CAE用語解説

本質的な違いは「時間」です。砂地盤の即時沈下は主に弾性変形ですが、粘性土の圧密沈下は、間隙水が絞り出される「排水」プロセスに支配されます。例えば、厚さ10mの粘土層に盛土を載せると、最終沈下量が30cmだとしても、その90%に達するのに数年から数十年かかることがあります。この時間遅れを予測するのが圧密解析です。

最も重要なのは「圧密係数（Coefficient of Consolidation）」

テルツァーギの一次元圧密理論が基礎です。過剰間隙水圧

その通りです。地盤中の全間隙水圧は「静水圧 + 過剰間隙水圧」です。荷重がかかると過剰分が発生し、これが排水とともに消散していきます。初期条件として、例えば瞬間載荷では載荷直後の過剰間隙水圧増分は載荷応力と等しくなります。境界条件は、砂層など排水面では

はい。それが「有効応力解析」または「u-p（変位-圧力）定式化」の核心です。未知数は節点変位

圧密解析では無条件安定な陰解法（後退差分法）が必須です。時間刻みは、現象の時間スケールと要素サイズから決まる「制限」があります。大雑把な目安として、

両方の手法があります。「モノリシック解法」は全未知数を一度に解き、安定ですが計算コストが高い。「分割解法」、特に「UW（Updated Weight）法」は、変位と圧力を交互に解き、効率的です。例えば、Abaqus/Standardの圧密解析では、デフォルトでモノリシック直接法ソルバーを使用しますが、大規模モデルでは反復ソルバーと前処理を組み合わせた手法も選択できます。重要なのは、分割解法を使う場合、小さな時間刻みを設定して「分割誤差」を抑えることです。

最低限必要なのは「土質定数」と「圧密定数」です。土質定数は、単位体積重量、ポアソン比、弾性係数または圧縮指数

界面の扱いが沈下予測の精度を左右します。砂層は透水性が極めて高いため、圧密解析では通常「完全排水層」とみなします。つまり、砂層内の間隙水圧は常に静水圧（過剰間隙水圧=0）と仮定します。FEMモデルでは、砂層要素に非常に大きな透水係数（例：1.0×10⁻² m/s）を設定するか、あるいは砂層の節点に直接排水境界条件（p=0）を適用する方法があります。粘土層と砂層の界面節点は、排水境界条件が適用される点になります。

通常、現地盤は長い年月をかけて自重による圧密は完了している、つまり「正常圧密状態」にあると仮定します。したがって、初期有効応力は土被り重量から計算されたK₀状態（静止土圧状態）で設定します。重要なのは、この初期有効応力と初期間隙水圧（静水圧）の和が全応力となり、これが初期の平衡状態を構成することです。ソフトウェアでは「地盤生成（Geostatic）」ステップでこの平衡状態を正確に計算させます。ここで収束しない場合は、パラメータや境界条件の見直しが必要です。

以下の順で検証します。1. **入力定数**：特に

核となるu-p定式化は同じですが、対象規模や前処理・後処理の思想が大きく異なります。**Abaqus/Standard**は汎用FEMソルバーで、極めて複雑な構成則や多物理場連成が可能ですが、地盤モデルの前処理は手間がかかります。**Plaxis 2D/3D**は地盤構造物解析に特化しており、自動メッシュ生成や階段状載荷の設定、施工段階の定義が直感的です。**GeoStudio SIGMA/W**はシンプルなUIと比較的安価なライセンスが特徴で、一次元・二次元の標準的な圧密解析に向いています。

これはPlaxisの重要な概念です。「非排水」は総応力解析で、間隙水圧の変化を考慮しません。「非排水(U)」は有効応力解析の一種で、間隙水圧の発生・消散を計算しますが、解析ステップの間は「透水係数=0」と仮定する、つまりそのステップ中は排水が起こらない設定です。例えば、盛土施工中（数日間）は「非排水(U)」ステップで過剰間隙水圧の発生を計算し、施工後の放置期間は「排水」ステップに切り替えて圧密過程を計算する、といった使い方をします。

「沈下量」と「沈下の時間経過」で影響が分かれます。弾性モデルでは、沈下量は応力に比例しますが、実際の粘土は圧縮指数

計算エンジンとしては可能です。例えば**Code_Aster**には飽和多孔質媒質の力学モデル（THMの一部）が実装されています。しかし、実用上の大きな壁は「前処理」と「パラメータ設定の煩雑さ」です。商用ソフトのように地盤層を直感的に定義し、自動で初期地盤応力場を生成する機能は乏しいです。研究目的で完全な制御性を求める場合や、独自の構成則を実装したい場合は有力な選択肢ですが、実務設計で効率性と確実性を求めるなら、PlaxisやGeoStudioなどの専用ソフトが圧倒的に優れています。

最も多い原因は「初期条件と境界条件の不整合」です。具体的には、1. 設定した初期間隙水圧（通常は静水圧分布）と、底面の排水/非排水条件が矛盾している。2. 側方境界の水平変位拘束と、自重による初期水平応力（K₀状態）が釣り合っていない。対策として、まずは全ての境界をローラー支承（水平変位固定、鉛直変位自由）で実行し、鉛直方向の平衡が取れるか確認します。その後、実際の境界条件に戻していきます。また、材料の剛性が低すぎる（非常に柔らかい粘土）場合も収束が悪くなります。

これは「マンドゥル（Mandel）-クライエル（Cryer）効果」という、二次元・三次元圧密に特有の現象で、実は非物理的ではありません。載荷直後、排水面付近の要素では間隙水圧がすぐに消散し有効応力が増加して収縮しますが、領域内部の要素では間隙水圧が上昇し、平均有効応力が一時的に「減少」することがあります。これに伴い、ポアソン比の効果で内部領域が一時的に膨張します。一次元解析では起こりません。結果がおかしいと判断する前に、間隙水圧と有効応力の時空間分布を詳細に確認してください。多くの商用ソフトの検証例題にもこの現象は登場します。

この症状は、過剰間隙水圧がほとんど消散した後（圧密度U>90%）で発生しがちです。残りの微小な圧密を追うためにソルバーが刻みを極端に小さくし（例えば1秒以下）、収束判定が難しくなるためです。対策は二つ。1. **解析終了条件の見直し**：実務上、圧密度95%や過剰間隙水圧が初期値の5%以下になった時点で解析を打ち切る。2. **最小時間刻みの設定**：ソフトウェアのオプションで、例えば「最小刻み=0.1日」などを設定し、それ以上小さくしないようにする。圧密終盤の微小な変化は、対数時間軸で外挿して推定するのが現実的です。