接触チャタリング（状態振動）

接触状態が「離れる→接触する」を短時間で繰り返す不安定現象だ。例えば、10ミリ秒の解析タイムステップで、ある節点が「接触→非接触→接触」と1ステップ内で状態を変える。物理的には、接触面に微小な隙間や弾性変形があり、ソルバーが平衡状態を見つけられずに振動する。ギアの歯面やクラッチ板で実際に発生する「バズノイズ」の原因の一つでもある。

接触条件の不等式制約と運動方程式の連成だ。ペナルティ法では、接触力

良い指摘だ。数値的に発生する「偽のエネルギー」だ。物理的には減衰があるが、数値モデルで減衰が不十分だったり、ペナルティ剛性が高すぎると、侵入した変形状態が持つ弾性エネルギー

一般的には、条件が悪いとどちらも起こるが、陽解法（LS-DYNAのExplicitやAbaqus/Explicit）の方が顕著に出やすい。陽解法はタイムステップが小さいため、状態の切り替わりを細かく「見てしまう」。例えば、10^-6秒ステップで接触・非接触が毎ステップ変われば、それは1MHzの高周波振動として観測される。陰解法（Abaqus/StandardやANSYS Mechanical）はニュートン・ラフソン法の反復中に接触状態を固定する「パーシスタンス」戦略をとるので、多少は抑制される。

非常に大きい。サーチ距離（サーチボックス）を大きくしすぎると、本来接触しないはずの面を候補に含め、不必要な接触力が発生する起点になる。逆に小さすぎると、実際に接触する瞬間を捉えられず、大きな侵入を許してしまう。ANSYSでは「pinball radius」、Abaqusでは「contact tracking」の設定だ。実務では部品間クリアランスの2〜3倍を初期値とし、調整する。

ケースによる。自動調整（Abaqusの「AUTOMATIC」、ANSYSの「program controlled」）は便利だが、局所的に硬い要素があると、そこだけ極端に高い剛性が設定され、チャタリングの原因になる。手動設定では、接触面周辺の最も柔らかい部材のヤング率と要素サイズから概算する。例えば、ゴム（E=5MPa）と鋼（E=210GPa）の接触なら、ゴム側の剛性を基準に

まず、接触面の初期隙間または初期侵入をチェックせよ。CADアセンブリでは「きっちり合わせ」だが、現実にはクリアランスや面圧による変形がある。ANSYSの「Contact Tool」やAbaqusの「Initial Status」で、初期ギャップが0.1mm以下か、あるいは0.01mm程度の微小侵入がある状態から始めるのが現実的だ。完全にぴったり（ゼロギャップ）は、数値的に不安定になりやすい。

接触面同士のメッシュサイズとパターンを可能な限り一致させること。片方が粗く、片方が細かいと、細かいメッシュの節点が粗いメッシュの要素面を「叩く」ような挙動になり、チャタリングを誘発する。また、二次要素（中間節点あり）は接触面の幾何学追従性は良いが、接触圧力の計算で振動を起こしやすいことがある。最初は一次要素で安定性を確認し、必要なら二次要素に移行するという手順が確実だ。

Abaqus/Standardなら「STATUS」ファイルの接触状態変化（closed to open, open to closed）の回数だ。1インクリメント中に同じ接触対で何度も状態が変わっていれば、ほぼ間違いなくチャタリングしている。ANSYSの「Solver Output」では、接触要素の実効ペナルティ剛性が極端に変動していないか確認する。また、反力や接触圧力の時系列データをプロットし、高周波のノイズ（例えば、平均10Nの反力が±5Nで振動）が出ていないか目視確認する。

各社とも「接触減衰」または「安定化」機能を用意している。ANSYS Mechanicalは「Normal Damping Factor」と「Stabilization」タブにある「Damping Factor」で、剛性マトリックスに人工減衰を加える。デフォルトは1e-4程度だが、0.01まで上げることもある。Abaqusは「Contact Controls」で「Stabilization」を選択し、減衰係数またはエネルギー比を指定する。LS-DYNAの陽解法では「CONTROL_CONTACT」カードの「SOFT」オプションや「VDA」オプションが、ペナルティ剛性を調整して衝撃を和らげる。

その通りだ。だからこそ、エネルギー監視が必須になる。Abaqusの「ALLSD」やANSYSの「Artificial Strain Energy」をモニターし、総ひずみエネルギーの1%未満に収まっているか確認する。例えば、総エネルギーが100Jの解析で、安定化による散逸エネルギーが5Jもあれば、それは過剰だ。あくまで数値的不安定を抑える「微量の糊」として使う。LS-DYNAのマニュアルにも「SOFT=2は最終手段」と書かれている。

商用ソフトのような高度な自動安定化機能は少ない。CalculiXでは「*CONTACT PAIR」に「STABILIZE」パラメータを追加できるが、基本的にはペナルティ剛性（

まず「接触剛性を下げる」だ。Abaqusなら「Contact Controls」の「Penalty Stiffness Scale Factor」を0.1にする。ANSYSなら「Normal Stiffness Factor」を「Program Controlled」から「Manual」に切り替え、値を0.1倍から始める。これで収束すれば、徐々に1.0に近づけて精度を上げる。同時に、最小タイムステップ/インクリメントサイズを小さくする（例えば1e-5秒や1e-3ステップ）。これで状態変化をソルバーがより細かく追跡できる。

接触面の「滑り」と「粘着」の切り替わりによる摩擦チャタリングだ。特に静摩擦係数と動摩擦係数を同じ値に設定していないか確認する。ISO 12156-1のような潤滑試験規格でも、静摩擦係数は動摩擦係数の1.2〜1.5倍だ。ANSYSやAbaqusで「Static/Kinetic Friction Ratio」を1.2に設定する。また、摩擦係数そのものが高すぎる（例えば鋼-鋼で0.5以上）場合、現実的には潤滑があるので0.1〜0.15に下げてみる。

最終手段は二つ。1. 接触アルゴリズムそのものを変える。ペナルティ法からラグランジュ乗数法（Abaqusの「Hard Contact」+「Lagrange Multiplier」）や増分ラグランジュ法に変更する。拘束条件を厳密に満たすのでチャタリングは起きにくいが、計算コストが上がり、接触離脱時の扱いが難しくなる。2. 物理的に微小な減衰を追加する。Abaqusの「Material Damping」でRayleigh減衰のβ係数に、最低周波数の逆数（例えば100Hzなら0.0016）程度の値を設定する。これらは「現象を変えてしまう」可能性があるため、結果の検証が必須だ。

エネルギー収支と反力の時定数を見る。まず、人工減衰や安定化による散逸エネルギーが全エネルギーの1%以下か。次に、反力の振動周期を計算する。接触面の局部剛性と質量から概算される固有周期