NAFEMS 接触ベンチマーク — V&V結果総括
ベンチマーク概要
先生、接触解析のベンチマークって他のV&Vとは難しさが違うんですか？
理論と物理
接触問題の基本概念
NAFEMSの接触ベンチマークで検証する「接触」とは、物理的にどう定義されるんですか？単に部品が触れ合うこととは違いますか？
厳密には「接触」は境界条件の一種です。二つの物体が互いに貫入しないという「非貫入条件」と、接触面で圧力が常に圧縮（引張りは生じない）という「圧力条件」、さらに接触面での摩擦の有無で定義されます。NAFEMSベンチマークでは、この条件を数値的にどう満たすかが焦点です。
非貫入条件を数式で書くとどうなりますか？また、それをソフトウェアが「厳密に」満たすのは難しいのですか？
シンプルな例では、接触点のギャップを
摩擦がある場合の条件はもっと複雑なんですか？NAFEMSではどんな摩擦モデルを扱っているんでしょう。
その通りです。摩擦があると、接触面のせん断力と滑り速度・変位の関係を定義する必要があります。NAFEMSのベンチマーク問題では、古典的なクーロン摩擦モデル
数値解法と実装
接触アルゴリズムの違い
教授がおっしゃったペナルティ法とラグランジュ乗数法の根本的な違いは何ですか？どちらが「正解」に近いんですか？
決定的な違いは、非貫入条件の扱い方です。ペナルティ法は小さな仮想的なバネで貫入を防ぎ、貫入量に比例した反力を発生させます。計算は比較的安定ですが、原理的に微小な貫入が許容されます。一方、ラグランジュ乗数法は追加の自由度（乗数）を導入して条件を厳密に満たそうとしますが、方程式が大きくなり、接触状態の判定が難しくなります。NAFEMSの結果を見ると、厳密解に近いのは多くの場合、Augmented Lagrangian法（両者を組み合わせた手法）を用いた高精度設定です。
メッシュの細かさは接触の結果にどれくらい影響しますか？要素サイズを半分にすると誤差も半分になりますか？
線形弾性問題のような単純な収束とは異なります。接触点や接触境界はメッシュに依存して離散化されるため、特に粗いメッシュでは結果が大きくぶれることがあります。NAFEMSの「LE10: Hertzian Contact」問題では、目標とする接触圧力分布を得るには、接触領域の少なくとも8要素以上が必要であると報告されています。誤差は単純に線形では減らず、アルゴリズムとメッシュの適合性が鍵です。
接触探索（サーチ）のアルゴリズムで、バケット法やモーテリー法と聞きますが、これらは精度ではなく計算速度のためですか？
主に計算速度のためですが、間接的に精度にも影響します。すべての節点と面の組み合わせをチェックする全探索（
実践ガイド
ベンチマークを模した解析設定
実際にAnsysやAbaqusで接触解析をする時、NAFEMSの知見からまず最初にチェックすべき設定は何ですか？
まずは「接触剛性（ペナルティ剛性）」または「ラグランジュ乗数の許容貫入量」の設定です。Ansysでは`FKN`、Abaqusでは「接触剛性スケールファクタ」に相当します。NAFEMSの結果を参照すると、デフォルト値（例えばAbaqusの1.0）では粗いメッシュで過大な貫入が生じる場合があり、問題に応じて0.1から10の間で調整が必要です。ただし、高くしすぎると数値的な悪条件（マトリクスが硬くなる）を招きます。
接触面の定義で、「マスター/スレーブ」や「面対面」「節点対面」といった選択肢がありますが、ベンチマークではどれが推奨されていますか？
NAFEMSのガイドラインでは、一般に「面対面」接触を推奨しています。これは接触圧力の分布をより滑らかに計算できるためです。「節点対面」は計算コストが低いですが、特にスレーブ面が粗い場合に「接触ロッキング」と呼ばれる過剰剛性化を引き起こすリスクがあります。マスター面はより剛直、またはメッシュが粗い面に設定するのが原則です。
収束が難しい時、荷重を一度に載せるのではなく、段階的に載荷する（ステップを分ける）のは有効ですか？
非常に有効で、これはNAFEMSの多くの問題でも採用されている基本戦略です。特に初期接触（ギャップや過干渉がある場合）では、まず小さな荷重で接触状態を安定させ、その後本荷重を載せます。Abaqus/Standardでは「ステップ」機能、Ansys Mechanicalでは「分析設定」の「段階的荷重」または「自動ステップ」を活用します。1ステップで全荷重を適用すると、接触状態の急激な変化により収束しないことが多々あります。
ソフトウェア比較
主要ソルバーの接触実装特性
Ansys Mechanical、Abaqus/Standard、COMSOLで同じNAFEMSベンチマークを解いた場合、結果に顕著な差は出るものですか？
デフォルト設定では差が出ますが、設定を適切にチューニングすればいずれも参照解に十分近づけられます。特徴を言うと、Abaqus/Standardの接触アルゴリズム（デフォルトはAugmented Lagrangian）は非常にロバストで産業界で広く信頼されています。Ansys Mechanicalはペナルティ法がデフォルトですが、詳細設定で貫入量を厳密に制御できます。COMSOLは「接触ペア」の設定が直感的ですが、大変形・大滑りを伴う問題では他の二つに比べて設定がやや難易度が高い印象です。
「ロバスト」とは具体的にどういうことですか？Abaqusが特別に優れている点は？
「ロバスト」とは、メッシュ品質や初期ギャップの設定が多少悪くても、とりあえず収束して物理的に妥当な答えを返す傾向が強い、ということです。Abaqus/Standardは、接触状態の判定と方程式の更新を非常に注意深く行うアルゴリズムを持っており、自動的にステップサイズを調整する能力に長けています。ただし、その分、AnsysのExplicit Dynamics (LS-DYNA) のような明示解法に比べると計算時間はかかる傾向があります。
無償/低価格のCAEソフト（CalculiX、Code_Asterなど）では、接触解析の信頼性はどう評価されていますか？
CalculiX（Abaqus風の入力）やCode_Asterは、NAFEMSベンチマークの多くを公式にパスしており、コアなアルゴリズムは高価格ソフトと同等の能力を持っています。しかし、前処理（メッシュ生成）や後処理、そして何より「収束しない場合の自動的な対処能力」やユーザーサポートが商用ソフトに比べて限られます。つまり、ユーザー自身が理論とソルバーの詳細を理解して設定する必要が高く、初心者にはハードルが高いと言えます。
トラブルシューティング
収束不良と結果の異常値
接触解析で「収束しない」というエラーが頻発します。最初に疑うべき原因と、具体的な対策を教えてください。
まず疑うのは「初期過干渉」です。部品が最初からめり込んでいる状態だと、ソルバーは大きな反力を計算し、発散します。対策1: モデルを確認し、クリアランスを設ける。対策2: ほとんどのソフトに「初期接触調整」機能があり（Abaqusの`*CONTACT INTERFERENCE`、Ansysの「接触ツール」）、初期段階で過干渉を解消するように設定します。NAFEMSのガイドでも、初期幾何定義の重要性が強く指摘されています。
収束はしたのですが、接触圧力の分布が非常にギザギザで、物理的に不自然です。これはメッシュが粗いからですか？
メッシュが粗い可能性は高いですが、もう一つの主要因は「接触条件の離散化」です。節点対面接触を使用している場合、この現象（圧力振動）が起きやすいです。対策としては、1) 面対面接触に変更する、2) メッシュを細かくする（特に接触領域）、3) 接触圧力の結果を「節点平滑化」する（後処理機能）があります。NAFEMS LE11問題では、この圧力振動を抑えるためのメッシュ要件が具体的に議論されています。
摩擦を入れると急に収束しなくなります。摩擦係数を0.1のような小さな値にしてもダメな場合、どうすれば？
摩擦は強い非線形性を加えます。対策は段階的です。第一に、まず摩擦無し（
反力の総和が外力とつり合っていない（例えば数%ずれている）場合、これは許容誤差ですか？それとも設定ミスのサインですか？
接触解析では、数%の不つり合いは設定ミスのサインである可能性が高いです。特に、接触領域が荷重伝達経路から外れている（接触が成立していない）、またはペナルティ剛性が低すぎて大きな貫入が生じ、実質的な「遊び」ができている場合です。まず、接触状態の可視化結果を確認し、意図した面で接触が確立しているかチェックしてください。NAFEMSの検証では、力のつり合い誤差は通常0.5%以内に収めることが目標とされます。
なった
詳しく
報告