Cam-Clayモデルで非排水せん断解析が発散する原因は何ですか？

臨界状態線を超えた軟化域で応力-ひずみ関係が非凸となり、ヤコビアン（接線剛性マトリックス）が負定値になるためです。これによりNewton-Raphson法の反復計算が振動・発散し、解析が不安定になります。

Cam-Clayモデルの収束性を改善する方法はありますか？

SloanとRandolph（1980年代）が提案した「サブインクリメント法」が有効です。1ステップの応力積分を100〜200の微小ステップに分割して計算精度を高めることで、安定性を確保できます。現在のAbaqusなど主要CAEソフトも同様の手法を採用しています。

AbaqusでCam-Clayモデルを使う際の注意点は？

非排水せん断解析では、材料の軟化による数値的不安定性に注意が必要です。Abaqusは内部でサブインクリメント法を適用して安定化を図っていますが、収束しない場合は要素サイズや増分ステップの見直し、また材料パラメータ（特に臨界状態線）の設定を再確認してください。

サブインクリメント法とは何ですか？

Cam-Clayモデルなどの弾塑性解析で数値的不安定を防ぐ応力積分手法です。1つの解析ステップを100〜200などの微小なサブステップに分割し、各サブステップで塑性ひずみや応力を逐次更新することで、高精度な積分を実現し、Newton-Raphson法の収束性を大幅に改善します。

Cam-Clayモデル — トラブルシューティングガイド

カテゴリ: 構造解析 | 2026-02-20

この記事は統合版に移行しました
より充実した内容を cam-clay.html でご覧いただけます。

CAE visualization for cam clay troubleshoot - technical simulation diagram

Cam-Clayのトラブル

🎓

初期状態で破壊する → 初期応力が降伏面の外にある。$K_0$ と $p_0$ の整合を確認

沈下が過大/過小 → $\lambda$ と $\kappa$ の値を確認。圧密試験のデータと比較

排水/非排水の区別 → 排水: 有効応力で計算。非排水: 間隙水圧を考慮

Coffee Break よもやま話

収束しない圧密解析の落とし穴

Cam-Clayで非排水せん断解析を行うと、臨界状態線を超えた軟化域でヤコビアンが負定値になりやすく、Newton-Raphson法が振動・発散する。1980年代にSloanとRandolphが「サブインクリメント法」を提案し、応力積分を100〜200の微小ステップに分割する手法で安定性を確保した。現在のAbaqusもこの考え方を踏襲している。

トラブル解決の考え方

「解析が合わない」と思ったら

まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
最小再現ケースを作る——Cam-Clayモデルの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う

Cam-Clayモデル — トラブルシューティングガイド

Cam-Clayのトラブル

収束しない圧密解析の落とし穴

トラブル解決の考え方

「解析が合わない」と思ったら

関連記事

関連トピック