クリープ座屈 — トラブルシューティングガイド

クリープ解析は時間積分の問題だから、通常の非線形解析とは異なるトラブルが起きる。

2つの可能性がある：

1. 座屈に近づいている — クリープ変形が加速する第3期クリープの状態。時間ステップが小さくなるのは正常な挙動。変位の急増を確認

2. クリープパラメータが不適切 — 応力指数 $n$ が大きすぎるか、$A$ の値が間違っている。非現実的に速いクリープが発生して数値的に不安定になる

変位-時間グラフを描く。座屈なら変位が指数的に増加する。パラメータエラーなら最初のステップからおかしな値が出る。まず弾性座屈荷重と現在の荷重レベルの関係を確認し、座屈が起こり得る荷重レベルかどうかを判断する。

確認項目：

1. クリープ材料モデルが正しく定義されているか — Norton則のパラメータ（$A, n$）が正しいか

2. 温度が設定されているか — クリープは温度依存。温度を与えていないとクリープが発動しないソルバーもある

3. Abaqusの場合、ステップタイプは VISCO か — STATIC ステップではクリープが計算されない

4. 時間の単位 — クリープ速度の単位（/秒 or /時間）と解析の時間単位が一致しているか

クリープパラメータが /時間 で定義されているのに、解析が秒単位で実行されると、クリープ速度が3600倍になる。逆だとクリープが見えないほど遅くなる。単位系の一貫性は必ず確認すること。

Norton則は定常クリープ（第2期）の近似だ。第1期クリープ（遷移クリープ）が重要な問題では、時間硬化則やひずみ硬化則を使う必要がある。

• 時間硬化則 — クリープひずみが時間の関数。荷重変動がない定荷重問題に適切
• ひずみ硬化則 — クリープひずみが累積ひずみの関数。荷重変動がある問題に適切

荷重が一定ならどちらでも同じ結果。荷重が変動する場合はひずみ硬化則がより正確。迷ったらひずみ硬化則を使うのが安全だ。

当然ある。クリープ変形で形状が変わると応力分布も変わり、それがさらにクリープ速度に影響する。座屈に近づくほどこの相互作用が強くなる。NLGEOM=YES を必ず設定すること。

NLGEOM=NOでは形状更新がないため、クリープ変形で構造が不安定化する効果が反映されない。クリープ変形は蓄積されるが、座屈は検出されない。クリープ座屈を評価するならNLGEOM=YESは必須だ。

• 時間ステップの急縮小 — 座屈の前兆か、パラメータエラーかを変位-時間グラフで判別
• クリープ変形ゼロ — 材料パラメータ、温度設定、ステップタイプ、単位系を確認
• クリープモデルの選択 — 荷重変動ありならひずみ硬化則を使う
• NLGEOM=YES は必須 — 大変形効果なしではクリープ座屈は検出されない
• 単位系の一貫性 — クリープ速度の時間単位と解析の時間単位を一致させる

その通り。クリープ座屈解析では何が起きていないかにも注意を払う必要がある。結果がおとなしすぎるときは、設定の見落としを疑うべきだ。