平面ひずみ問題 — トラブルシューティングガイド
より充実した内容を plane-strain.html でご覧いただけます。
平面ひずみ解析のトラブル
平面ひずみ解析でよくあるトラブルを教えてください。
平面応力のトラブルに加えて、平面ひずみ特有の問題がある。
体積ロッキング
$\nu$ が大きい材料で変位が異常に小さく出ます。
体積ロッキングだ。平面ひずみの最も代表的なトラブル。
確認方法：
- 応力コンターにチェッカーボードパターン（市松模様）が出ていないか
- $\nu$ を小さくすると結果が大きく変わるか（$\nu = 0.3$ と $\nu = 0.499$ で比較）
対策：
1. 低減積分要素に切り替え — CPE4R（Abaqus）
2. ハイブリッド要素に切り替え — CPE4H / CPE8RH
3. 二次要素の低減積分 — CPE8R（最も安定）
4. $\nu = 0.5$ ではなく $\nu = 0.499$ — 完全な非圧縮を避ける（バルク弾性率を有限に保つ）
$\nu = 0.499$ と $\nu = 0.5$ で差がありますか？
数値的には大きな差がある。$\nu = 0.5$ では $1/(1-2\nu)$ が無限大になり、ソルバーが正常に動作しない。$\nu = 0.4999$ 程度にすれば実質的に非圧縮だがソルバーは安定する。ただしハイブリッド要素なら $\nu = 0.5$ でも問題ない。
初期地圧の設定ミス
地盤解析で掘削前の初期状態がおかしいです。
初期地圧の設定ミスは地盤解析で最も多いトラブルだ。
確認項目：
- $K_0$ の値は正しいか — 正規圧密粘土で $K_0 = 1-\sin\phi'$（Jaky式）
- 自重と初期応力が整合しているか — 初期ステップで大きな変位が出たら不整合
- 水圧の設定 — 全応力か有効応力かで $K_0$ の適用先が変わる
初期ステップで変位がゼロにならないのは問題ですか？
問題だ。正しく初期地圧が設定されていれば、初期ステップの変位はゼロ（すでに平衡状態）になるはず。変位が出るということは、自重による応力と設定した初期応力が一致していない。
Abaqusでは *GEOSTATIC ステップで初期平衡を自動的に確認できる。Plaxisでは $K_0$ procedure で初期状態を自動生成する。汎用FEMで手動設定する場合はこのチェックが特に重要だ。
掘削解析で応力が飛ぶ
掘削ステップで要素を除去した直後に応力が異常に大きくなります。
掘削面の応力解放が一度に行われると、特に弾塑性地盤ではオーバーシュートが起きる。
対策：
- 掘削を複数ステップに分割 — 1ステップで1 m掘削など
- 応力解放率を段階的に — Plaxisの $\beta$ パラメータで応力解放を段階化
- 増分数を増やす — 各掘削ステップの増分を細かく
結果の単位に注意
平面ひずみの結果の単位でよくあるミスは？
最も多いのは力の単位だ。平面ひずみは奥行き1単位あたりの解析だから：
- 反力は「kN/m」（奥行き1 mあたり）
- 集中荷重は「kN/m」として入力する（線荷重と同等）
- 実構造の全幅に換算するには奥行き長さを掛ける
「1 kNの集中荷重」を与えたつもりが、実は「1 kN/mの線荷重」だった…。
そう。平面ひずみの「集中荷重」は実際には奥行き方向に一様な線荷重だ。面荷重も同様で、入力値の単位が kN/m² なら結果は kN/m で出る。この換算を報告書で明記すること。
まとめ
平面ひずみのトラブル対処、整理します。
体積ロッキングと初期地圧が二大トラブルですね。
前者はFEMの数値的な問題、後者は地盤力学の物理的な問題。両方理解していないと平面ひずみの地盤解析は成り立たない。
平面ひずみ境界条件の設定ミス
平面ひずみ解析でz方向変位を全節点に固定する誤りは初心者に多い。正しくはεz=0を構成則レベルで強制し、変位境界条件は与えないのが原則だ。ANSYS PLANE182では要素タイプのKEYOPT(3)=2（Plane strain）を正しく設定しないと平面応力として計算される。この設定ミスによる応力差は最大でν/(1-ν)≈43%（ν=0.3の場合）になりうる。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——平面ひずみ問題の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告