接触オーバークロージャ（初期貫通）

オーバークロージャ（初期貫通）は、解析開始時点（t=0）で、接触条件を定義した2つの部品が、すでに互いにめり込んでいる状態を指します。例えば、公差設計で0.1mmの隙間を想定した軸と穴を、誤って0.05mmの干渉状態でCADモデリングしてしまうと、これが初期貫通になります。物理的には、部品が重なって存在するという、現実にはあり得ない初期状態を意味します。

多くの接触アルゴリズム、特にペナルティ法では、大きな問題になります。接触力は「貫通量×ペナルティ剛性」で計算されます。初期貫通量が大きいと、最初の反復計算で非常に大きな仮想的な反力が発生し、

それは接触検出アルゴリズムに依存します。Ansys Mechanicalの「ノード対面」検出では、スレーブ節点がマスター面に貫通しているかで判定します。一方、「面対面」検出では、ガウス積分点での距離を評価するため、節点が貫通していなくても要素中央で干渉していると検出されることがあります。後者の方が一般的に精度が高く、特に二次要素を使う時に重要です。初期貫通量は、この検出点での最短距離が負の値として計算されます。

主要なソルバーはいくつかの自動補正オプションを提供しています。例えば、Abaqus/Standardの「*CONTACT INTERFERENCE」では、指定したオーバークロージャ量を初期荷重条件として扱い、最初のインクリメントで徐々に解放します。Ansysの「Automatic Bisection」は、最初のサブステップをより小さなステップに分割して貫通を解消しようと試みます。しかし、これらは「応力のない初期状態」という仮定を崩すため、アセンブリ応力が生じ、結果に影響を与える可能性があります。

理論的には、ペナルティ剛性

はい、以下の順で確認することを推奨します。
1. **CAD段階**: アセンブリの干渉チェック（SolidWorksの「干渉検出」やCATIAの「Clash」機能）を実行し、意図しない干渉（例えば0.001mm以上）がないか確認。
2. **メッシュ段階**: 接触面のメッシュサイズが粗すぎないか確認。粗いメッシュはジグザグ形状になり、見かけ上の貫通を生む。
3. **プリ処理段階**: ソフトウェアの初期接触状態可視化ツールを使う（Ansysの「Contact Tool」、Abaqusの「Initial Contact Status」）。「Near Field」設定値（デフォルトはメッシュサイズの倍数）を適切に調整し、接触ペアが正しく認識されているか確認。
4. **接触設定**: 「Adjust to Touch」や「Slave Adjustment」などのオプションで、微小な隙間/干渉を自動補正する機能を試す。

まず、その貫通が「物理的な干渉」なのか「メッシュによる見かけ上の干渉」なのかを区別します。可視化ツールで貫通量を定量化します。0.01mm以下の微小な貫通で、かつアセンブリ上で隙間が意図されている場合は、ソフトウェアの調整機能（Ansysの「Pinball Region」内での「Adjust to Touch」）で対応可能です。しかし、0.1mm以上の明確な干渉や、圧入配合などを模擬する意図的な干渉の場合は、CADモデルを修正するか、解析設定で「Interference Fit」として明示的に定義する必要があります。安易な自動調整は初期応力を生み、疲労解析などでは重大な誤差になります。

そのために接触ツールがあります。例えば、Ansys Mechanicalの「Contact Tool」で「Initial Information」を表示し、「Status」で「Overclosed」のみをフィルタリングします。Abaqus/CAEでは、「Interaction」モジュールで「Contact Controls」の「Interference Check」を実行します。より高度な方法として、貫通量をコンター図で表示させ、色で深刻度を判断します。多くの場合、問題は1〜2箇所の局部に集中しています。それらを修正した後、再度全体チェックを実行するのが効率的です。

各ソフトの哲学が反映されています。
**Ansys Mechanical**: デフォルトの「Program Controlled」では、微小な貫通に対して自動的に「Adjust to Touch」を適用しようとします。ユーザーは「Interface Treatment」で「Adjust to Touch」「Add Offset」「Ramped Effects」などを細かく設定可能。「MPC」や「Beam」接触では初期貫通を許容しません。
**Abaqus/Standard**: より厳格で、初期貫通は警告またはエラーとして報告されます。明示的な対策として「*CONTACT INTERFERENCE」を用い、オーバークロージャ量を指定して解放します。「*CONTACT CONTROLS」で「AUTOMATIC TOLERANCE」を設定することも可能です。
**COMSOL Multiphysics**: 「接触ペナルティ法」の設定で「初期ギャップ/干渉」を直接入力するフィールドがあります。初期貫通を「負のギャップ」として定義し、それを解消するように解析を開始できます。柔軟ですが、ユーザーが物理的に正しい値を知っている必要があります。

オープンソースソルバーでも基本的な機能はありますが、GUIでの可視化や自動修正が弱く、ユーザーの知識がより求められます。**CalculiX**（Abaqusと互換性のある入力形式）では、`*CONTACT PAIR` に `INTERFERENCE` パラメータを設定できます。**Code_Aster**（Salome-Meca環境）では、`DEFI_CONTACT` コマンドで `DEPL_INIT_MAIT` や `DEPL_INIT_ESCL` を用いて初期位置補正を定義します。しかし、商用ソフトのような直感的な「初期状態診断ツール」は乏しいため、コマンドやPythonスクリプトを書いて貫通量を事前計算するなど、手間がかかることが多いです。

以下の3点を即座に確認してください。
1. **接触面の定義距離（Pinball Radius, Search Distance）**: これが小さすぎると、本来接触すべき面が認識されず、別の面と干渉していると誤検出されることがあります。デフォルト値（多くの場合「Program Controlled」）から始め、必要に応じてメッシュサイズの2〜3倍に広げてみます。
2. **部品の初期位置**: アセンブリ内での部品の配置がCADと一致しているか。特にインポート時の単位系ミス（mmとmの取り違い）は、1000倍の貫通を生みます。
3. **メッシュ品質**: 接触面付近に極端に歪んだ要素（アスペクト比>20）がないか。これにより節点位置がおかしくなることがあります。

その場合、手動での位置調整を検討します。具体的な手順は：
1. **局所メッシュ細分化**: 貫通している接触面周辺のみメッシュを細かくし、形状を正確に表現。
2. **仮想部品の挿入**: 貫通が生じている2部品の間に、厚さ0.01mmなどの仮想的な「シム」要素を挿入し、物理的な隙間を作る。解析後、その要素の影響を評価する。
3. **初期ステップの分割**: 最初のロードステップを10や100などのより小さなサブステップに強制分割し、ソルバーがゆっくりと接触状態を解決できるようにする（Abaqusの「Stabilization」、Ansysの「Auto Time Stepping」の積極的利用）。
4. **接触アルゴリズムの変更**: ペナルティ法からAugmented Lagrange法やMPC法に切り替える。これらは初期貫通に対する耐性が異なります。

ほぼ間違いなく信用できない結果になります。主に以下の問題が発生します：
1. **非物理的な高応力**: 最初の反復で巨大な接触力が発生し、局部に現実とはかけ離れた高応力（例えば降伏応力の10倍）が生じる。
2. **エネルギーの不整合**: 偽りのひずみエネルギーが発生し、特に動解析では系のエネルギー保存則が破綻する。
3. **収束不良と不安定性**: ソルバーが振動し、収束せずに停止するか、誤った平衡状態に落ち着く。
「結果が得られた」という出力自体が危険です。常に初期接触状態のレポートを確認し、貫通量が許容範囲（通常、メッシュサイズの1%以下）内であることを確認することが、信頼性のある接触解析の絶対条件です。