Terzaghiの式による比較

Ga3d改変コードの検証の続きです。

Terzaghiの式では、Nc, Nq, Nγ の3つの支持力係数を使用します。このうち、理論解が求められるのは前２つですので、γ=0として（Nγの影響を取り除いて）比較しました。

ピーク強度 c=10kN/m2、φ=30°kN/m2、残留強度 c=１kN/m2、φ=25°kN/m2 です。
結果は以下の通り。

値は合っていますが、いきなり破壊するんですね。粘土の方が降伏してからも、文字通り粘りが出るんですね。

ただこれは「ある収束値での修正NRの繰り返しが○○回以上で発散」という人的な設定値によるものも大きいと思います。今は500回なので、1000回ならもう少し伸びるでしょうか？あるいはcを落とさずにφだけ落としたらどうなるでしょう？破壊時におかしな値になっているのも気になります。

もう少し、やってみましょう。