摩擦なし接触の不安定性

良い質問だ。物理的には、接触面に摩擦が存在しない（

離散化された線形方程式系

その通り、接触チャタリングだ。時間積分（過渡解析）では、あるタイムステップで接触が成立し、次のステップでわずかに離れる、を高速で繰り返す。これは、接線方向の運動エネルギーを吸収するメカニズム（摩擦）が存在しないため、わずかな摂動でも平衡状態を維持できないからだ。Abaqus/Explicitで摩擦係数0.0の接触を設定すると、しばしば「CONTACT FORCE」の出力がステップごとに激しく正負に振動するのが観測される。

大きな違いがある。ペナルティ法は接触面にバネを置くようなものなので、接線方向にも微小なバネ剛性

代表的なのは「人工的な接線剛性の付加」と「減衰の導入」だ。例えば、MSC NastranのSOL 400では、`BSTIF` パラメータで接線方向の「ダミー剛性」を指定できる。また、Abaqus/Standardの「*CONTACT STABILIZATION」では、接触面に粘性減衰力を導入し、非物理的な剛体運動を抑制する。この減衰係数は、初期過渡期のみ作用させ、定常状態には影響しないように設計されている。

一般的には、陽解法（Abaqus/Explicit, LS-DYNA）の方が顕著に影響を受ける。なぜなら、陽解法は条件付き安定なため、接触による急激な剛性変化やチャタリングが安定限界タイムステップ

まず第一に、「本当に摩擦係数=0が必要か」を物理的に再考せよ。多くの実務では、0.01〜0.05の微小な摩擦を設定しても結果に大きな影響はない。どうしても0が必要なら、以下のワークフローを推奨する：1) 静解析から始め、接触安定化を有効にする（Ansysなら`Stabilization`、Abaqusなら`*CONTACT STABILIZATION`）。2) 非常にゆっくりとした擬似静的な過渡解析で挙動を確認する。3) 接触反力（特に接線方向）の履歴をモニターし、非物理的な振動がないかチェックする。4) メッシュ依存性を確認。粗いメッシュほど不安定になりやすい。

その懸念は正しい。Ansys Mechanicalでは、安定化エネルギーと全体のひずみエネルギーの比を`Solution Information`で監視せよ。経験則では、この比が1〜5%を超えないように調整する。Abaqusでは、`*CONTACT STABILIZATION`の`DAMPING FACTOR`をデフォルト（1E-4程度）から始め、収束する最小の値に設定する。最終的には、安定化をオフにした解析（微小摩擦あり）と結果を比較し、主要な応力や変位が例えば5%以内で一致することを確認するのが確実な検証法だ。

ある。第一に、接触面のメッシュを整合させる（matched mesh）ことが有効な場合が多い。これにより接線方向の「すべり」が要素エッジに沿って起こり、数値的なノイズが減る。第二に、要素形状はできるだけ正方形や正三角形に近づける。細長い要素（アスペクト比>10）は接線方向の変形モードを不正確に評価し、不安定性を助長する。第三に、Ansysの`CONTA174`やAbaqusの`*CONTACT PAIR`では、サーフェス平滑化（Surface Smoothing）オプションを有効にすると、接触力の計算が滑らかになり、安定性が向上することがある。

デフォルトの戦略が異なる。Ansys Mechanicalは、摩擦なし接触に対して自動的に微小な「摩擦ダンピング」を導入する傾向がある。詳細は`Solution Information`の「Stabilization」項を確認せよ。Abaqus/Standardは、より厳格で、不安定性を検知するとエラーメッセージ`***ERROR: TOO MANY ATTEMPTS`を出して停止する傾向が強い。COMSOL Multiphysicsの「接触」物理場では、デフォルトで「増強ラグランジュ法」が使われ、接線方向にペナルティパラメータを導入するオプションがあり、比較的ロバストだが、ユーザーが意識して設定する必要がある。

以下が主要なものだ。

その通り、ユーザー自身が数値的処置をより意識する必要がある。CalculiX（CCX）では、デフォルトの接触アルゴリズムはペナルティ法だが、摩擦なしの場合の安定化機能は限定的だ。`*CONTACT PAIR`で`STABILIZE`パラメータはない。代わりに、`*DAMPING`カードで全体または局部に数値減衰を加えるか、現実的な微小摩擦を設定するのが現実的だ。Code_Asterでは、`DEFI_CONTACT`コマンドで`ALGO_CONT='PENALISATION'`（ペナルティ法）を選択し、`COEF_MULT_IMPE`（ペナルティ係数増加倍率）を調整して安定性と精度のバランスを取る必要がある。ドキュメント`U4.44.01`に詳細がある。

高い確率で関連している。これらの警告は剛性マトリクスの特異性を示す。まず`/STATUS, CINT`コマンドで接触状態を確認せよ。摩擦なし接触が疑われる場合、対策は段階的だ：1) `CNCHECK`で接触ペアの初期貫通や法線方向を確認・修正。2) 接触設定で`KEYOPT(2)=1`（MPCアルゴリズムの使用）に変更してみる。3) ソルバー制御で`EQSLV, , ,1E-8`のようにピボット許容値を緩和する（最終手段）。根本的には、接触安定化を有効にするか、微小摩擦を導入することを検討せよ。

これは摩擦なし接触によるチャタリングの典型的な症状だ。ALLAE（人工ひずみエネルギー）は安定化や時間積分に伴う数値的エネルギーで、これが内部エネルギー（ALLIE）の5〜10%を超えると問題ありだ。対策：1) `*CONTACT`定義に`*CONTACT CONTROLS, STABILIZE`を追加。2) 接触面のスケーリング係数（`*CONTACT PAIR, ADJUST=位置調整値`）を使用し、初期貫通を防ぐ。3) 全体の体積粘性スケーリング（`*BULK VISCOSITY`）をデフォルトから増やす（例：Q1=1.5, Q2=1.0）。4) 根本的には、接触面のメッシュを細かくし、`*SECTION CONTROLS`で`TYPE=ENHANCED`を指定して接触検出精度を上げる。

有力な原因ではあるが、断定する前に他の可能性を潰すべきだ。手順：1) **モデル単純化**：同じ接触条件を単一要素や2Dモデルで再現し、現象を分離する。2) **拘束不足の確認**：摩擦なし接触面の接線方向と、モデル全体の剛体運動モードが重なっていないか。例えば円筒面の接線方向は、軸方向の並進や回転と一致する可能性がある。3) **接触法線方向の確認**：法線方向が一貫していない（例えば自己接触で反転する）と、接線方向の定義も不安定になる。4) これらのチェック後も問題が続くなら、摩擦係数を0.001に変更してみる。それで収束すれば、摩擦0が原因とほぼ断定できる。その後、安定化オプションの導入を検討すればよい。

それは接触圧力とせん断応力の関係による。クーロンの摩擦則