潰れた要素エラー

有限要素法では、要素の形状が極端に歪むと数値積分が正しく行えなくなり、ソルバーが計算を停止します。特に、要素の体積がゼロまたは負になる状態を「潰れた要素」と呼びます。例えば、四角形要素の4つの節点がほぼ一直線上に並んだり、四面体要素の体積が計算精度の限界を下回ったりする場合です。

計算上の話です。要素の節点座標からヤコビ行列の行列式、つまり要素の体積を計算します。変形が大きすぎて節点の順序が入れ替わると、この行列式が負の値になります。例えば、シェル要素が180度以上反転した「裏返り」状態が典型的です。Abaqus/Explicitのマニュアルでは、要素のアスペクト比が1000:1を超えるような極端な変形で発生すると記載されています。

要素剛性マトリックス

各ソルバーは内部で許容される要素の歪みの限界値を設定しています。例えば、Ansys Mechanicalの非線形ソルバーでは、デフォルトで要素の歪みが1.0（つまり100%）を超えると警告を出し、設定によっては計算を停止します。また、要素体積が初期体積の約1e-12倍以下になると、数値積分点での重みが実質的にゼロと見なされ、エラーとなります。

Lagrangian形式では限界があります。大変形解析で有効なのはALE（Arbitrary Lagrangian-Eulerian）法や、SPH（Smoothed Particle Hydrodynamics）のようなメッシュレス法です。例えば、金属成形シミュレーションのデファクトスタンダードであるAutoFormは、独自のインクリメンタルメッシング技術で、ツールに押し付けられて極端に歪む板金メッシュを随時再構成し、潰れエラーを回避しています。

まずは要素タイプの選択です。完全積分要素はロッキングを起こしやすく潰れエラーの原因になります。減積分要素（AbaqusのC3D8Rなど）の使用が基本です。また、初期メッシュ品質が悪いとすぐに潰れるので、アスペクト比は5以下、スキュー角は60度以下、ワーピング量は許容範囲内に収めるように前処理ソフト（ANSAやHyperMesh）でチェックしてください。

非常に重要です。現実の材料はある程度の圧縮に対して体積がほぼ不変（非圧縮性）ですが、ソルバーがこれを厳密に扱うと数値的に硬くなり、反発力が急増して潰れエラーを誘発します。これを防ぐため、Abaqusのハイパーエラスティック材料や多くの粘塑性モデルでは、体積弾性率をせん断弾性率の約100〜1000倍に設定し、準非圧縮性として扱います。安易にデフォルト値を使わないでください。

「ペナルティ法」の剛性が大きすぎると、接触面が衝突した瞬間に局所的に極めて大きな応力が発生し、要素を押しつぶすことがあります。Ansysでは「FKN」接触剛性係数を下げて（例えば0.01など）、より柔らかい接触に設定します。また、「自動二分法」を有効にし、増分サイズを自動的に小さくして、接触の瞬間の衝撃を和らげる設定が有効です。

メッセージと哲学が異なります。

はい、ソルバーごとに独自の拡張機能があります。

1. エラー要素の特定: ログファイルや結果ファイル（Abaqusの.odb、Ansysの.rst）を開き、最後の収束した増分ステップで、どの要素/部品が極端に変形しているかを可視化します。変形図に加え、要素歪み（Element Strain）のコンターを表示させると一目瞭然です。

必ずしもそうではありません。むしろ、メッシュを細かくしすぎると、小さな要素がより早く極限変形に達し、エラーが早期に発生する可能性があります。まずは、メッシュの粗い領域と細かい領域の移行をなだらかにする「メッシュバイアス」を調整してください。急激なメッシュサイズの変化が応力集中と大変形を引き起こしているケースが多いです。

その場合は、Lagrangianメッシュの限界を認め、別のアプローチに切り替える必要があります。選択肢は3つです。