変位過大エラー

いい質問だ。単なる警告ではない。これはソルバーが、ある増分ステップで計算した変位が、要素の物理的サイズを超えるような非現実的な値になったことを示す。例えば、辺長10mmの要素が、1ステップで15mmの変位を計算してしまった場合だ。これは材料の構成則（応力-ひずみ関係）が、その大きな変形に対して適切に定義されていない、あるいは境界条件や荷重が現実的でないことを強く示唆している。

線形弾性モデル

材料非線形性を考慮しても、ソルバーが解を探す「増分ステップ」の過程で、一時的に非物理的な状態を通ることがある。例えば、Abaqus/Standardのデフォルト設定では、ある増分ステップで平衡反復が収束しない場合、自動的に増分サイズを25%に切り下げて再試行する。しかし、問題の設定そのものが悪い（例えば、拘束不足で剛体運動が可能）場合、増分サイズをどれだけ小さくしても、解が発散してこのエラーに至る。

主要なソルバーは、予測子-修正子法に基づくNewton-Raphson法を用いている。この反復計算の各サイクルで、変位増分の大きさを監視している。Ansys Mechanicalのデフォルトでは、ある反復で計算された変位増分の最大値が、現在の要素サイズの20倍を超えると「Large Deflection」警告を出し、さらに大きくなるとエラーで停止する。この「20倍」という値は、ソルバー設定で変更可能だが、物理的に無意味な計算を続行させないための安全装置だ。

全くの別物だ。収束判定は、平衡方程式

可能性はあるが、根本解決にはならない。確かに、Ansysでは「Initial Substeps」を1000など大きな値に、Abaqusでは「Initial Increment」を1e-6など小さな値に設定すれば、第一歩目の変位増分は小さく抑えられる。しかし、問題の根源が「剛体運動モードの存在」や「材料モデルの不安定性」にある場合、2ステップ目、3ステップ目で再びエラーが発生する。計算コストが膨大になるだけで、対症療法に過ぎない。

まずは「拘束不足」を疑え。特に、静解析でこのエラーが出た場合、最も多い原因だ。部品が宙に浮いており、剛体運動（並進・回転）を抑制する拘束が足りていない。具体的には、ISO 10791-7などの機械の剛性試験規格を参考に、現実の治具を模した拘束条件を適用する。3次元モデルでは、最低6つの自由度（3並進、3回転）を過不足なく拘束する必要がある。

次は「単位系の不一致」だ。これは初学者に非常に多い。材料データ（ヤング率：GPa）、幾何（mm）、荷重（N）の単位が整合しているか。例えば、ヤング率を210 GPa (=210,000 N/mm²)、幾何をmm、荷重を1 Nとした場合、変位はµmオーダーになる。もしヤング率を210 Paと誤入力していれば、変位は極端に大きくなり、エラーの原因となる。常にSI単位系（m, kg, s, N, Pa）か、mm単位系（mm, tonne, s, N, MPa）のいずれかに統一すること。

大いにある。接触条件が適切に定義されていないと、部品が互いに貫通し、その反動で異常な変位が発生する。Ansysの「Frictional」接触で「Normal Stiffness Factor」がデフォルトの1.0だと、硬い部品同士の接触で数値的不安定性を起こすことがある。10.0や0.1に調整してみる。また、Abaqusの「Hard Contact」は圧力-過剰侵入量関係の勾配が無限大なので、「Softened Contact」に変更して現実的な接触剛性を定義するのも手だ。

出る。線形静解析は連立一次方程式

Ansysでは、「Solver pivot warnings」や「An excessive deformation (13.4mm) of element 45217 has occurred」といった具体的なメッセージが`.err`ファイルに出力される。Abaqusでは、`.msg`ファイルに「***ERROR: TOO MANY ATTEMPTS MADE FOR THIS INCREMENT」と表示され、その前の行に「THE SOLUTION APPEARS TO BE DIVERGING. CONVERGENCE IS JUDGED UNLIKELY.」や「EXCESSIVE DISPLACEMENT AT NODE...」といった詳細が記される。Abaqusのメッセージの方が原因を特定しやすい構成だ。

あるが、根本治療ではない。Ansysの「Nonlinear Controls」では、「Line Search」をオンにすると、変位増分の大きさを自動調整して発散を抑制できる。Abaqus/Standardでは、ステップ定義で「Stabilization」タブから粘性減衰を導入できる。これは仮想的なダンパーを追加して剛体運動を抑制する方法で、特に接触問題の初期安定化に有効だが、エネルギーが追加されるため結果の解釈には注意が必要だ。COMSOLでは、「定式化」を「完全ラグランジュ」から「更新ラグランジュ」に変えることで大変形の追従性を改善できる場合がある。

基本的な数値アルゴリズムは同じなので、同様の問題は発生する。MSC Nastranの非線形ソルバーSOL 400では、「FATAL ERROR 4080 (GFLAS2)」として「EXCESSIVE DISPLACEMENT」が報告される。そのマニュアルには「拘束不足または材料不安定性をチェックせよ」と明確に書かれている。RADIOSS（Altair製品）では、エラーメッセージと共に、変位が大きくなっている特定のノードや要素のIDを出力し、可視化ツールで直接その部位を強調表示できる機能があり、デバッグがしやすい。

典型的なケースだ。ボルト軸方向には拘束されていても、ボルトの「ねじ山」とナット（または雌ねじ）の接触が適切に定義されていない可能性が高い。単純な摩擦接触だけでは、回転方向の剛性が不足し、ボルトが「空転」するような数値的不安定性を生む。対策は二つ。第一に、Ansysの「Bolt Pretension」やAbaqusの「Bolt Load」といった専用機能を使い、軸力と締め付けトルクの関係を適切にモデル化する。第二に、ねじ山の詳細モデルを作成する代わりに、「Beam Connection」や「Spring」要素でねじ部の軸方向・せん断方向剛性をJIS B 1082やVDI 2230規格に基づいて定義する。

これもよくある。超弾性体モデル（例えばMooney-RivlinやOgdenモデル）は、試験データからフィッティングされたひずみエネルギー密度関数

メッシュが粗すぎたために生じた「ロッキング」や「アワーグラッシング」と呼ばれる数値的困難が緩和されたからだ。特に、体積非圧縮性を持つ材料（ゴム、金属塑性）や、せん断変形が支配的な領域で、一次要素（四面体、六面体とも）を使うと、過剰剛性（ロッキング）により、物理的には柔らかい変形モードを表現できず、代償として別の自由度で異常な変位が発生することがある。メッシュを細かくする、または二次要素（中間節点付き）に変更することで、要素の変形自由度が増え、真の変形モードをより正確に表現できるようになり、発散が回避される。

あるが、強く推奨しない。Ansysの「Solver Pivot Warning」を無視する設定や、Abaqusで「*CONTROLS, PARAMETERS=FIELD」を編集して発散判定の閾値を緩めることは技術的には可能だ。しかし、これは「エンジンチェックランプを切って走り続ける」ような行為である。得られた結果は物理的に無意味で、設計判断に用いてはならない。CAEの目的は現実を模倣することであり、数値的不安定性を無理やり押し通すことではない。エラーは問題設定の不備を教えてくれる貴重なシグナルとして捉え、モデルと物理的理解を修正する機会とすべきだ。