支持力係数

反復回数1000回でも、結果は全く変わりませんでした。

ｃ＝10kN/m2で固定した場合の結果は、以下の通り。

値はOKです。やはり、粘りがないですね。施工時の観測では困るでしょうね。弾性域は200kN/m2程度でしょうか？

これを含めた検証結果より、改変コードは正しく動作していると言えるでしょう。

ところで、支持力係数は道示ではグラフからの読み取りになっており、式が示されていません。古い土研資料が根拠のようですが、持っていないですね。
http://jsce.jp/pro/node/1212
今回は、建築基礎構造設計指針の Nc, Nq 式を使いました。

支持力係数の算定式はいくつかあります。式やグラフは、いろんな参考書に載っています。
技報堂出版「土の力学」は10年以上前の講習会で購入した教科書ですが、そこに係数のグラフが載っています。当時、そこに書き込みをしていました。

「Terzaghi は理論的に少し違う。Caquot・Kerisel が正解。」

極限支持力算出式で、根入れ分の土の荷重は見るが、強度を見ない所が理論的でない（安全側）といったことは覚えています。しかし、支持力係数の「理論的に少し違う」は覚えていないですね。何だったのでしょう？

理工図書「土質力学-全訂新版-」でも式が以下のように書かれています。

「上述の解析は本来、有効応力法であるクーロンの破壊基準に基づき、またプラントルの支持力破壊の手法を踏襲しているので、支持力の解析も有効応力法であるが、第4章で述べたような全応力法による見かけのせん断強さを用いても差し支えないことは、土圧や斜面の安定の場合と同様である。それらの結果に差があるのは現在ではやむを得ない。」

これは1991年の4版なので、ちょうど20年前ですね。20年経ってもあまり変わっていないです。