歪んだ要素エラー

良い質問だ。要素の歪み（Distortion）は、理想的な形状（例えば正方形や正四面体）からどれだけ要素が歪んでいるかを定量的に評価する指標だ。具体的には、要素の頂点座標から計算されるヤコビ行列の行列式や、アスペクト比、内角など、様々な幾何学的パラメータで評価される。例えば、四面体要素の歪み度合いを表す指標の一つに、次のようなものがある。

本質的な理由は2つある。第一に、数値積分の精度の低下だ。FEMではガウス積分点で応力やひずみを評価するが、歪んだ要素では積分点の物理的な位置が偏り、積分精度が大幅に落ちる。第二に、形状関数の導関数の計算が不安定になる。ヤコビ行列が特異に近づくと、その逆行列の計算で数値誤差が爆発し、剛性マトリクスが正しく計算できなくなる。例えば、ある四面体要素のヤコビ行列式が1.0e-10以下になると、倍精度計算でも有効桁数がほとんど失われ、ソルバーが発散する。

その通りだ。各ソルバーは内部的な閾値を持っている。例えば、Abaqus/Standardは四面体・六面体要素に対して「形状チェック」を実行し、デフォルトで歪み許容度を0.1〜1.0の範囲で警告やエラーを出す。一方、Ansys Mechanicalでは「Element Quality」という指標を多用し、0.0〜1.0で評価する。1が理想（正四面体や直方体）で、一般的に0.1以下は「Poor」と判定される。これらの閾値は、線形静解析と非線形・動解析では異なることが多い。動解析の方がより厳しい基準が要求される。

基本的には全要素に対するチェックだ。プリプロセッサ（例えばAnsys MeshingやAbaqus/CAE）がメッシュ生成時に、あるいはインポート時にバッチ処理で実行する。具体的なアルゴリズムとしては、各要素タイプ（四面体、六面体など）ごとに定義された「品質指標」を計算する。四面体なら先ほどの体積比や、各辺の長さの比（アスペクト比）、内角の最大・最小値などを計算する。六面体の場合はより複雑で、例えば「歪みパラメータ」

Jacobian Ratioは、要素内の複数地点（通常は積分点や頂点）で計算したヤコビ行列式の最小値と最大値の比だ。理想的な直方体要素では、全ての地点でヤコビ行列式は同じ値になるので比は1になる。要素が歪むと、地点によって値がばらつき、比が1から離れる。Ansysでは、二次要素（中間節点あり）の評価で特に重要だ。例えば、四面体二次要素（SOLID187）では、頂点と辺の中点、体積中心など、多数のサンプリング点でヤコビ行列式を計算する。許容値は解析タイプによるが、非線解析や接触を含む解析では30以上になると警告、100を超えるとエラーとなることが多い。

圧倒的に二次要素の方が敏感だ。その理由は形状関数にある。一次要素の形状関数は線形なので、要素内の変形勾配は一定だ。しかし二次要素では形状関数が二次式であり、中間節点の位置が要素内の変形マッピングを大きく左右する。中間節点が辺の中央から大きくずれていると、要素内部のマッピングが高度に非線形になり、ヤコビ行列式が場所によって極端に変化する。したがって、AnsysやAbaqusでは二次要素を使う場合、「中間節点を辺の中点に置く」オプション（Abaqusでは「Midside node parameter」）をデフォルトで有効にすることが推奨される。これを無効にすると、わずかな幾何形状の凹凸でも歪みエラーが多発する。

まず、エラーメッセージに記載されている要素IDを特定し、その要素がモデルのどの物理的な位置にあるかを可視化せよ。Ansysなら「Select > By Element Number」、Abaqus/CAEなら「Query > Element」で探せる。その要素が「なぜそこに生成されたか」を考える。鋭い角、非常に細長い面、2つの部品が接近している接触領域などが典型的な発生箇所だ。

幾何形状の修正が第一選択肢だ。CADでその鋭角なリブに0.5mmや1mmの面取り（フィレット）を追加する。現実の製品にも面取りはあることが多い。CAD修正が不可能なら、メッシュコントロールで対処する。Ansys Meshingでは「Face Meshing」や「Edge Sizing」を使って、そのリブに沿って要素を規則的に配置する。または、メッシュの種類を変える。四面体（Tet）メッシュで歪む場合、その領域だけをスイープメッシュで六面体（Hex）にできないか検討する。どうしてもダメな場合は、メッシュサイズを全局的に細かくする。細かくすると歪みは軽減されるが、計算コストが上がる。

必ずしも全て修正する必要はないが、戦略的な判断が必要だ。まず、それらの低品質要素が「高応力集中が予想される領域」にあるかどうかを確認する。ボルト穴の縁やノッチ根部など、結果が重要な領域に低品質要素がある場合は、たとえ警告レベルでも修正すべきだ。一方、拘束面や遠方の領域で、応力がほとんどかかっていない場所にあるなら、無視しても問題ないことが多い。ただし、非線形解析（材料非線形、接触）や固有値解析（モーダル解析）では、低品質要素が数値的不安定性を引き起こしやすいため、より厳格に対応する。

状況による。メッシュスムージングは、要素の内部節点を移動させて形状を改善する手法で、Abaqus/CAEやHyperMeshに実装されている。ただし、境界上の節点は移動しないので、根本的な形状問題は解決できない。また、過度にスムージングをかけると、要素が互いに貫通したり、本来の幾何形状から大きく逸脱するリスクがある。実務では、まず局所的なリメッシュやメッシュサイズ調整を行い、それでもダメな場合の最終手段として、軽度のスムージングを適用する。Ansys FluentなどのCFDソルバーでは、メッシュスムージングは一般的だが、FEMでは慎重に使うべき機能だ。

ユーザーインターフェースとデフォルトの厳しさに違いがある。Ansys Mechanicalは「Mesh Metric」として「Element Quality」「Skewness」「Orthogonal Quality」など多角的な指標を提供し、ヒストグラムで分布を見せるのが特徴だ。Abaqus/CAEは「Mesh」モジュールの「Verify」タブで「Analysis checks」を実行し、「Distortion」の度合いを要素ごとに色分け表示する。COMSOL Multiphysicsは「Mesh」ノードの「Statistics」で「Element quality」を表示し、0〜1の値で示す。デフォルトの警告閾値は、Abaqusがやや保守的（厳しい）で、Ansysは標準的、COMSOLはユーザー設定に委ねる部分が大きい印象だ。

大きく違う。Abaqus/Standardは歪んだ要素を検出すると、デフォルトでは解析を停止しエラーメッセージを出力する。しかし、キーワード「*DIAGNOSTICS, DISTORTION CONTROL=YES」を入力ファイルに追加すると、歪んだ要素を自動的に削除（あるいは修正）して計算を続行しようとする「歪み制御」機能が働く。Ansys Mechanicalでは、歪みがひどい要素があると、ソルバーが剛性マトリクスの生成に失敗し、「A general failure occurred in the solver」のような一般的なエラーで止まることが多く、要素特定までが少し手間になる。COMSOLは「歪んだメッシュ要素」に対する自動リメッシュ機能が強力で、適応型メッシュ細分化と組み合わせやすい。

オープンソースソルバーでは、プリプロセッサ（メッシャー）とソルバーの連携が商用ソフトほど密でないため、ユーザー自身が責任を持ってメッシュ品質を管理する必要がより高い。例えば、CalculiX（プリプロセッサはCGXなど）では、メッシュ生成時に品質チェックコマンドを実行するが、その基準は比較的単純だ。Code_AsterはSalome-Meca環境と連携するが、Salomeのメッシュモジュールで提供される品質指標（NICO、Quality）を解釈する知識が要求される。共通して言えるのは、エラーメッセージが商用ソフトより簡素で、「歪んだ要素」というより「ゼロまたは負の体積要素」といった直接的な表現でエラーが出ることが多い点だ。

これはほぼ間違いなく、CADデータの「幾何学的欠陥」が原因だ。STLファイルはファセットデータなので、面同士が完全に接続されていない「ギャップ」や、重なっている「オーバーラップ」がある。STEP/IGESファイルでも、非常に短いエッジ（長さ0.001mmなど）、微小面、サーフェス間のわずかなずれ（ギャップ）が含まれていることがある。メッシャーはこれらの不完全な幾何形状を忠実になぞろうとして、極端に細長い三角形や四角形を作り出し、結果として歪んだ要素が生成される。対策は、プリプロセッサの「幾何修復」機能（Ansysの「SCDM」、Abaqusの「Geometry Edit」）を使って、微小な特徴を削除したり、ギャップを縫合したりすることだ。

これは「幾何学的非線形性」に起因する典型的な問題だ。材料が大きく変形するにつれて、要素自体も引き伸ばされ、押しつぶされ、著しく歪む。初期状態で良くても、変形後の形状でヤコビ行列式が負（要素が裏返し）になることがある。対策はいくつかある。第一に、要素公式を見直す。Abaqusではハイブリッド要素（例：C3D8H）は体積変化の大きい材料に有効だ。第二に、メッシュを細かくする。変形をより細かい要素で追跡できる。第三に、ソルバー設定を変更する。Ansysでは大変形オプション（NLGEOM,ON）をONにし、さらに「Mesh Morphing」や「Rezoning」（Abaqusでは「Adaptive Meshing」）といった、変形途中でメッシュを再構成する機能を利用する。これらは計算コストが高いが、解析を継続させるためには必要だ。

接触計算は、要素の表面法線や侵入量の計算に要素形状が直接関わるからだ。歪んだ要素では表面法線の方向が急激に変化したり、接触判定のための積分点配置が偏ったりする。その結果、接触力の計算が不安定になり、収束しないか、あるいは異常に大きな接触力が発生して、さらに要素を歪ませる悪循環に陥る。具体的な対策としては、(1) 接触面には可能な限り六面体要素や、品質の高い四面体要素を使用する、(2) 接触面のメッシュサイズをそろえる（サイズ比が10:1以上だと問題が起きやすい）、(3) 初期接触設定を「Adjust to Touch」などにして、初期貫入を解消する、(4) 接触アルゴリズムを「Augmented Lagrange」から「Pure Penalty」に変えてみる（場合による）、などがある。

「歪み」以外のメッシュ品質指標を確認せよ。特に「アスペクト比」と「スキュー角」だ。細長い要素（アスペクト比が大きい要素）は、一方向の変形剛性が過大評価され、応力が真値からずれる。また、「スキュー角」が大きい（要素の辺が直交から大きくずれている）要素では、せん断変形の計算に誤差が生じる。Ansysの「Mesh Metric」でこれらを一括チェックできる。もう一点、結果評価の方法も重要だ。要素内部の応力はガウス積分点で計算され、節点に外挿される。歪みの大きい要素では、この外挿処理自体が誤差を増幅する。したがって、歪みが大きい領域の「要素解」ではなく「節点解（平均化済み）」を見ても、見かけ上は滑らかだが、それは誤った平均化の結果かもしれない。根本的解決は、その領域のメッシュ品質を向上させるしかない。