退化要素エラー

良い質問だ。典型的なのは、四面体要素の1つの辺が潰れて、4つの節点がほぼ同一平面上に並んでしまうケースだ。例えば、理想的な四面体の体積は

自動メッシャーのアルゴリズムに起因する。特に、複雑なCADジオメトリを「テトラメッシュ」で埋めようとすると、狭い隙間や鋭いエッジの近くで要素が押しつぶされる。また、メッシュの「スイープ」操作で、掃引経路が急激に変化する部分でも発生しやすい。Ansys Meshingの「Patch Conforming」アルゴリズムはこの問題に敏感で、ジオメトリの細かい特徴を忠実に追従しようとするあまり、微小な面の近くに極端に扁平な要素を生成してしまう。

数値積分の精度が致命的に低下するからだ。FEMでは要素内の応力やひずみを計算するためにガウス積分点を使用する。退化した要素では、形状関数のヤコビアン行列

両方の段階でチェックされるが、基準と深刻度が異なる。前処理メッシャー（例：Altair HyperMeshの「Quality Index」チェック）は、形状比、歪み、最小角、ヤコビアンなど多角的な指標で「低品質」要素を検出し、警告や自動修復を試みる。一方、ソルバー（例：MSC Nastran）は剛性マトリクスを組み立てる際に要素剛性マトリクスの計算を試み、数値的に特異（singular）と判断した場合に致命的エラー（FATAL MSG 3009など）を出力して解析を停止させる。

ソルバー内部の要素ライブラリルーチンで行われる。例えば、等パラメトリック要素では、各積分点でヤコビアン行列式

大きな違いがある。二次要素（中間節点あり）は形状関数が高次であるため、辺が曲がって表現できる代わりに、中間節点の位置が不適切だと更容易にヤコビアンが負になる「inverted element」を生じやすい。一方、一次要素はヤコビアンが節点位置だけで決まるので、検出は比較的単純だが、一度退化すると積分点自体の値がおかしくなる。このため、非線形解析で大きな変形が生じ、要素が潰れる過程では、二次要素の方が早い段階でエラーを吐く傾向がある。

まず、エラーメッセージに記載されている要素IDを特定せよ。ほとんどのソルバー（例：Ansys Mechanicalの「Solution Information」）はエラーを起こした要素の番号を出力する。その要素を可視化し、どのようなジオメトリ特徴（鋭い角、非常に小さい面、薄肉部など）の近くにあるかを確認する。これが根本原因の特定の第一歩だ。

ある。これを「ジオメトリの理想化」または「特徴の抑制」と呼ぶ。例えば、Siemens NX Simcenter 3DやCATIAのCAEモジュールでは、「Virtual Topology」機能を使って、解析に不要な小さなエッジや面を隣接する大きな面にマージできる。具体的には、直径1mm以下のボスやリブ、0.5mm以下のフィレットは、構造強度に寄与しないことが多く、抑制対象の候補だ。ただし、応力集中が予想される重要な箇所のフィレットは残す必要がある。

幾つかある。第一に、メッシュサイズを調整する。問題の小さな面の近くで「ローカルメッシュサイズ」を設定し、その面の寸法（例えば0.5mm）よりも少し大きいサイズ（例えば1.0mm）でメッシュを切ると、メッシャーがその面を無視して通過する可能性がある。第二に、メッシュタイプを変える。四面体（テトラ）から六面体主体（主にヘキサ）のメッシュに変更する。スイープメッシュやマップドメッシュは、要素のアスペクト比を制御しやすく、退化を生じにくい。ただし、形状が複雑だと適用できない場合もある。

業界の経験則がある。線形静解析の場合、最小ヤコビアン比（理想的な形状に対する比）が0.5以上あればまずまず、0.7以上あれば良好とされる。しかし、非線形解析や接触解析ではより厳しい基準が必要で、0.8以上を目標とする。歪み（Skewness）については、四面体で0.9以下、六面体で0.8以下が一般的な目標値だ。これらの値は、NASAの「Structural Analysis Reference Guide」や自動車業界の内部ベストプラクティスにも記載されている。

メッシュ修復の哲学と自動化の度合いに違いが顕著だ。Ansys Workbenchは「Mesh Metrics」で詳細な品質指標を可視化し、「Patch Independent」メッシュ法という、小さなジオメトリ特徴をある程度無視してメッシュを生成するオプションを提供する。Abaqus/CAEのメッシャーはジオメトリに忠実だが、「Mesh Edit」ツールキットで個々の要素を手動で移動・分割する機能が強力で、熟練者がピンポイントで修復するのに向く。COMSOL Multiphysicsは「Free Tetrahedral」メッシュに「極端に細長い要素を避ける」アルゴリズムを内蔵しており、事前にある程度の予防がなされるが、一度発生した場合の診断情報はやや簡素だ。

その通りだ。オープンソースのメッシャー（例：Netgen、Gmsh）はアルゴリズム的には優秀だが、商用ソフトのような頑健な「ジオメトリヒーリング」機能が乏しい。そのため、少し汚いSTEPやIGESファイルを読み込むと、直接的に退化要素が生成されやすい。また、CalculiXなどのソルバーはエラーメッセージが簡素で「negative jacobian」としか出ず、問題の要素を特定するのに手間がかかる。商用ソフトの強みは、このような「汚い現実」のデータに対する前処理の耐性と、豊富な診断情報にある。

二つの理由がある。第一に、それらは「メッシュ品質管理」を最優先したUIと機能を備えている。例えば、ANSAの「Morphing」ボックスは、メッシュ節点を直接引きずって形状を修正しながら、リアルタイムで品質指標（ヤコビアン、アスペクト比など）の変化をモニターできる。第二に、業界特有のベストプラクティスが組み込まれている。自動車業界向けの「セミオートマチック中面メッシュ」機能は、薄板部品を効率的に中面化し、そこで発生しがちな退化要素を未然に防ぐアルゴリズムを内包している。

原因は類似しているが、数値的な状態が異なる。「Negative Jacobian」は、要素の節点順序がおかしくて実際に行列式の値が負になっている状態を示す。これは、大きな変形解析で要素が裏返った（inverted）時や、メッシュ生成時の節点番号付けの不具合で発生する。「Zero or Negative Jacobian」は、負であることに加えて、ゼロに近すぎる状態も含む。後者の方がより一般的で、単に要素が潰れた場合にも該当する。いずれにせよ、ソルバーはその要素の剛性マトリクスを正しく計算できないと判断してエラーを出す。

最終手段として幾つかあるが、結果の信頼性は低下するので注意が必要だ。第一に、ソルバー設定を変える。Abaqus/Explicitでは、デフォルトの「Element deletion」設定を変更して、潰れた要素を削除して解析を続行させる方法がある。第二に、要素タイプを変える。Abaqus/StandardでC3D4（線形テトラ）がダメなら、C3D4H（線形テトラ、ハイブリッド圧力定式化）に変えると、体積ロックが緩和され、多少の形状歪みに耐えられる場合がある。第三に、メッシュを極端に粗くする。問題の局部だけを選択し、メッシュサイズを2倍や3倍にすると、その小さな特徴をメッシュが飛び越えてしまう可能性がある。

これは「大変形による要素の破壊」が原因だ。初期状態では良好な形状比だった要素が、解析中に大きな変形や座屈を起こし、物理的に押しつぶされてしまう。特に、ゴム材料の超弾性解析や、金属の大塑性変形解析で頻発する。対策は二つ。一つは、メッシュ適応（リメッシュ）機能を使うこと。Abaqusの「ALE adaptive meshing」やLS-DYNAの「Adaptive Remeshing」は、変形が大きい領域のメッシュを解析中に自動的に再生成する。もう一つは、問題領域の初期メッシュを意図的に細かくし、変形に伴う形状歪みを分散させることだ。

間接的な原因になり得る。特に「ペナルティ法」による接触では、接触面に垂直な方向に非常に大きな仮想バネが挿入される。この時、接触面のメッシュが不適切（例えば、従属面のメッシュが粗すぎる）だと、その反力が局所的な節点に集中し、その周辺の要素を数値的に歪ませ、ヤコビアンを悪化させる可能性がある。対策としては、接触面のメッシュサイズを均一にし、従属面は主面より細かいメッシュにする（経験則では1.2倍から2倍細かい）ことが推奨される。また、Ansysでは「Normal Lagrange」法に切り替えることで、ペナルティ法による過剰な反力集中を緩和できる場合がある。