メッシュ検証 — CAE用語解説

いいえ、それだけではありません。メッシュ検証の本質は、離散化された領域が元の連続体の物理をどれだけ正確に表現できるかを評価することです。例えば、歪んだ要素（高アスペクト比要素）では、剛性マトリクスが悪条件になり、応力計算が不正確になります。具体的には、アスペクト比が20を超えると、せん断ロッキングのリスクが顕著に高まります。

その通りです。定量的な指標があります。例えば、四面体要素の「歪み（Distortion）」は、理想的な正四面体からの幾何学的なずれを0から1の値で評価します。Ansys Meshingでは「Element Quality」として、体積と辺長の比に基づく指標を採用しています。値が1に近いほど理想的で、0.1以下は多くのソルバーで警告の対象です。

鋭角が小さすぎると、数値積分点が極端に近づき、ヤコビアン行列式がゼロに近づくからです。これにより、要素の剛性が過大評価され、局所的に非常に硬い挙動を示す「剛性マトリクスの特異性」が生じます。Abaqusのドキュメントでは、三角形シェル要素の内角が10度を下回ると、結果の精度に重大な影響を与えると明記されています。

主な指標と関連誤差は以下の通りです。

ヤコビアンは、要素の形状関数を実空間座標で微分する際の変換行列の行列式です。等パラメトリック要素では、ガウス積分点でこの値が正であることが必須条件です。負のヤコビアンは、要素が「裏返っている」ことを意味し、剛性マトリクスの計算が完全に破綻します。式で書くと、あるガウス点でのヤコビアン

要素の形状と使用される数値積分法が異なるからです。四面体（テトラ）要素は一次要素でも比較的ロバストですが、歪みに敏感です。一方、六面体（ヘキサ）要素は高次要素で高い精度を発揮しますが、非直交性やスキューが生じやすく、これらが誤差の主要因となります。例えば、Autodesk CFDのヘキサドミナントメッシュでは、直交性が0.3以上であることが推奨されています。

実務では以下のステップが一般的です。

あります。学術論文では非常に厳しい値が使われがちですが、実務では計算コストと精度のトレードオフを考慮します。一例として、自動車業界でよく使われる目安は：

「要素品質は良いが、メッシュの急激な変化（サイズ比）が大きい領域」です。隣接する要素のサイズ比が3倍を超えると、特に過渡解析やCFDで数値振動を引き起こす可能性が高まります。もう一つは「接触面のメッシュ不一致」です。接触を定義する二つの面で、メッシュサイズやパターンが大きく異なると、接触探索が失敗したり、接触圧力が振動したりします。このチェックは自動化されていないことが多いので、目視確認が必須です。

大きく異なります。

プリプロセッサに依存します。CalculiXのフロントエンドであるCalculiX GraphiX（CGX）には、メッシュ統計を表示する機能があります。OpenFOAMでは、メッシュ生成ツール`snappyHexMesh`や`blockMesh`の実行時に、コンソールに最小/最大直交性、非直交性、スキュー度などの統計が出力されます。また、`checkMesh`というユーティリティが最も重要で、$ checkMesh -allGeometry -allTopologyと実行すると、網羅的な幾何学・位相学的チェックが行われ、問題のある要素の数と位置を詳細にレポートします。

両者とも具体的な改善策を提示します。

これは致命的なエラーです。即座に以下の手順を踏みます。

まず第一に直交性（Orthogonality）とスキュー角（Skewness）を疑います。非直交性が高いと、拡散項の離散化に用いる補間が不正確になり、数値拡散が増加または振動を引き起こします。次に、急激なメッシュサイズ変化がないか確認します。隣接セルの体積比が10:1を超えるような領域は、移流項の計算が不安定になります。OpenFOAMの`checkMesh`で報告される「non-orthogonality」の最大値が70度を超えている場合は、明らかな原因です。

これは古典的な「特異点」の問題です。鋭い凹角（ゼロ半径のフィレット）など、理論上応力が無限大になる幾何学的特異点が存在すると、メッシュを細かくするほど計算される応力は増加し続けます。メッシュの品質ではなく、幾何学そのものが問題です。対策は二つ：