初期貫通エラー

その通り、幾何学的に重なり合っている状態です。接触条件を定義した2つの部品（例えばボルトと穴）のメッシュが、解析開始時点（時刻0）で空間的に侵入し合っていることを指します。Abaqusのマニュアルでは「initial overclosure」、Ansysでは「initial penetration」と表記されます。接触面の法線方向で、主面が従面を貫通していると判断されます。

数値解法上の大きな問題です。ほとんどの接触アルゴリズムは、接触力

幾何学的形状と材料の変形を考慮する順序が逆です。現実のはめあいは、干渉量（例えばJIS B 0401のH7/p6）に従った幾何形状同士を、材料の弾性変形を利用して組み付けます。CAEでは、まず干渉のない（または微小クリアランスのある）メッシュ同士を配置し、解析ステップで干渉量を「与える」ことで、材料の変形挙動を「計算」します。初期状態で干渉させてしまうと、その変形プロセスを計算できず、いきなり非現実的な高応力が発生するのです。

基本的には、接触対に定義された「従面」の全ての積分点（ガウス点）またはノードについて、最も近い「主面」セグメントへの垂線距離を計算します。この距離が負（主面の法線ベクトルと逆方向）なら貫通と判定されます。Ansys Mechanicalのデフォルト設定では、この検出距離「Pinball Radius」は要素サイズに基づいて自動設定されますが、ユーザーが変更可能です。

それらは緩和策であって根本解決ではありません。「自動調整」（Abaqusの`ADJUST`、Ansysの`Adjust to Touch`）は、初期ステップで従面ノードを主面まで移動させ、貫通量をゼロにします。しかし、これはノード位置を「修正」するので、モデルの体積や質量が微妙に変化します。微小貫通（要素サイズの数%）の修正には有効ですが、大きな干渉を無理やり修正すると、不自然な初期応力（`INITIAL STRESS`）が発生し、後続の解析結果を歪めます。

部分的に正しく、危険でもあります。Ansysの「Softened Contact」やAbaqusの「Soft Contact」は、接触圧-オーバークロージャ曲線を線形または指数関数的にし、初期貫通による過大な反力を緩和します。しかし、これはあくまで「計算を続行させる」ための手法です。例えば、1mmの貫通をデフォルトの剛性で計算すると数GPaの圧力が発生しますが、ソフト接触で1/100に緩和しても、依然として数十MPaの非現実的な初期応力が生じた状態で解析が始まります。根本原因の修正にはなりません。

まずは「可視化」です。Ansysなら`Contact Tool`で`Initial Information`を、Abaqusなら`Interaction`モジュールで`Query`→`Contact constraints`を実行し、貫通しているノードとその量を色分け表示させます。次に、その位置の幾何形状（CAD）とメッシュを確認。CAD上でクリアランスがあるのにメッシュで貫通しているなら、メッシュ品質（歪んだ要素）が原因です。CAD自体が干渉しているなら、設計データの見直しが必要です。

要素サイズを小さくするのは一つの方法ですが、計算コストが増加します。より効果的なのは「メッシュ整合性」を高めることです。接触面同士でメッシュの密度とパターンを可能な限り一致させます。Ansysの`Match Control`やAbaqusの`Mesh Seeding`を使います。また、曲率の大きな部分ではメッシュを細かくし、要素のアスペクト比を1:10以内に収め、スキュー角を60度以下に抑えることで、メッシュの「見かけ上の貫通」を大幅に減らせます。

一般的な経験則として、局所的な要素サイズの1%から5%以下の貫通は、自動調整で修正しても結果への影響は無視できる場合が多いです。例えば、接触面の要素サイズが0.5mmなら、0.025mm以下の貫通は許容可能な範囲と考えることができます。ただし、これは線形静解析の場合で、非線形解析や動解析ではより厳しい基準が必要です。重要なのは、自動調整後に「初期応力が生じていないか」を必ず確認することです。

その場合は、初期状態では「接触」ではなく「微小なクリアランス」（例えば0.01mm）を設けたモデルを作成します。その後、解析ステップで「干渉量」を指定します。Abaqusなら`Contact Interference`、Ansysなら`Joint`の`Preload`や`Contact`の`Offset`機能を使用します。これにより、材料が変形しながら嵌め込まれる物理的なプロセスを、数値的に安定して追跡できます。初期貫通させて「はめあい応力」を再現するのは、絶対に避けるべき方法です。

大きく異なります。Ansys Mechanicalはデフォルトで「Adjust to Touch」がオフで、貫通があると警告またはエラーを出します。ユーザーが明示的にオンにする必要があります。Abaqus/Standardは、微小貫通に対して自動的に調整を試みる傾向があり、警告メッセージ(`***WARNING: OVERCLOSURE`)を出しながら計算を続行することが多いです。COMSOL Multiphysicsの「接触」ペアは、初期ギャップ/貫通を検出する機能があり、「補助スイープ」オプションで初期貫通ノードをサーフェス上に投影できます。

原理は同じですが、GUIでの可視化・修正ツールが貧弱なため、難易度は上がります。CalculiX（PrePoMaxフロントエンド）やCode_Aster（Salome-Meca環境）でも、接触ペナルティ法やラグランジュ乗数法は実装されています。問題は、貫通箇所を特定するのに、入力ファイル（.inp, .comm）の警告メッセージを読み解き、ノード番号と座標から位置を逆算する必要がある点です。商用ソフトのような直感的なコンタクトツールはほぼないため、メッシュ作成段階での予防が極めて重要になります。

Altair RADIOSSなどのExplicitソルバーを中心としたコードは、接触検出アルゴリズムそのものが異なります。タイムステップが非常に小さいため、初期貫通による「巨大な力」が一挙に発生せず、数ステップかけて「押しのけられる」挙動を示すことがあります。また、「サーフェス・オフセット」機能を使って、接触面をわずかにオフセットし、見かけ上の初期ギャップを作る手法が一般的です。ただし、これも質量や慣性の変化を招くため、動解析では注意深い検証が必要です。

主に2つの原因が考えられます。1つは、調整量が大きすぎて不自然な初期変形・応力が生じ、その状態が平衡状態ではなく、解析開始直後に「跳ね返る」ように変形するためです。もう1つは、接触面の定義が不適切で、調整により一部のノードは修正されたが、隣接するノードが新たに貫通する「調整の波及効果」が起きているためです。対策としては、まず調整量をログで確認し（Abaqusの.msgファイル、AnsysのSolution Information）、それが要素サイズに対して無視できない量（例：10%以上）でないか確認します。

「接触対を一時的に無効化して、1つずつ有効化しながら確認する」が基本です。Ansysなら`Contact Group`を作成して個別にオン/オフ。Abaqusなら`Model Tree`で`Interaction`をコメントアウト。それでも見つからない場合、CADの干渉チェック機能（SolidWorksの「干渉チェック」、CATIAの「Clash」）を事前に実行しておくことが最も効果的です。CAE前処理段階で、CADレベルでの干渉を排除しておけば、多くの初期貫通エラーは未然に防げます。

シェル要素は厚みを持たない「中間面」で定義されるからです。現実の板厚t=2mmの2枚の板が接触する場合、CADでは面同士の距離が2mmですが、シェルモデルでは両方の中間面が同じ位置に来てしまい、必然的に貫通状態になります。この場合は、接触オフセット（`Offset`）の設定が必須です。Ansysでは`Contact Surface Treatment`で`Add Offset`を指定し、値を`Radial`（半径方向）か`Uniform`で与えます。Abaqusでは`Contact Initialization`で`Specify tolerance:`を設定します。デフォルトは「無し」なので、必ず設定を確認してください。

「CAEモデルは、現実の物理状態をそのまま再現するものではない」という原則を理解することです。CAEは「計算可能な理想化されたモデル」です。初期状態は、力が加わっていない自然な状態（＝接触面同士が干渉せず、微小な隙間があるか、ちょうど接している状態）から始めなければなりません。全ての拘束や荷重、干渉は、解析ステップを通じて「適用する」ものです。この基本を守らずにCADの組み立て状態をそのままメッシュ化すると、高確率で初期貫通エラーに遭遇します。予防に勝る対策はありません。