接触解析 (Contact Analysis) — CAE用語解説
接触解析とは
接触解析ってどんなときに必要なんですか？ 普通の線形静解析と何が違うんですか？
部品同士が「触れたり離れたりする」現象を扱うときに必要だよ。例えばボルト締結体のフランジ面、歯車のかみ合い、プレス成形の金型とブランクなど。普通の静解析だと節点の拘束は最初から決まってるけど、接触解析では荷重の大きさや方向によって「どこが接触していてどこが離れているか」が変化する。これが境界非線形と呼ばれるもので、線形解析では本質的に扱えないんだ。
なるほど、境界条件が「解いてみないとわからない」ってことですね。具体的にどういう数式で接触を扱うんですか？
接触問題の本質は不等式拘束条件だ。法線方向のギャップを $g_N$ とすると、基本的な条件は3つ：
これをKarush-Kuhn-Tucker (KKT) 条件、あるいはHertz-Signorini-Moreau条件と呼ぶ。3番目の式がポイントで、「ギャップが開いていれば接触力ゼロ、接触力が出ていればギャップはゼロ」という相補性条件を表しているんだ。
定義と基本概念
FEMの設定で「master面」と「slave面」を選ぶじゃないですか。あれって何が違うんですか？ どっちをどっちにすればいいか迷います。
slave面の節点がmaster面に貫通しないように拘束するのが基本的な仕組みだ。だから一般的なルールは：
- 剛性が高い方をmaster（例：金型がmaster、ブランクがslave）
- メッシュが粗い方をmaster（細かいslave節点で接触を検知する方が精度が出る）
- 面積が大きい方をmaster
ただし最近のソルバー（Abaqus, Ansys等）ではsurface-to-surface定式化が主流で、主従の影響が小さくなってきている。それでも迷ったらこのルールに従っておくのが無難だよ。
接触タイプ（node-to-surface / surface-to-surface）
node-to-surface と surface-to-surface って、結局どっちを使えばいいんですか？
結論から言うと、特別な理由がなければ surface-to-surface (STS) を使うべきだ。理由を整理しよう：
Node-to-Surface (NTS):
- slave側の各節点がmaster面に投影される
- 計算は軽いが、接触圧分布がガタガタになりやすい
- master面がslave節点間を「すり抜ける」ことがある（ジッピング現象）
- 曲面が粗いメッシュだと精度が低下する
Surface-to-Surface (STS):
- 面全体で接触拘束を積分するので、接触圧が滑らかになる
- 主従依存性が小さい（どちらをmasterにしても結果がほぼ同じ）
- 計算コストはNTSより高いが、収束性は良い場合が多い
例えば自動車のボルト締結フランジ解析だと、NTSでは接触圧に節点ごとのチャタリングが出るけど、STSなら滑らかな分布が得られる。Abaqusのデフォルトも今はSTSだよ。
ペナルティ法と拡張ラグランジュ法
接触の設定で「ペナルティ法」と「拡張ラグランジュ法」を選べるんですが、何が違うんですか？
これは接触拘束 $g_N \geq 0$ をどうやって数値的に処理するかの話だ。
ペナルティ法は、貫通量に比例した反力を発生させる。仮想的な「接触バネ」を噛ませるイメージだ：
$$p_N = \epsilon_N \cdot \langle -g_N \rangle$$ここで $\epsilon_N$ はペナルティ剛性（接触スティフネス）、$\langle \cdot \rangle$ はMacauley括弧（負の値をゼロにカット）。メリットは追加の未知数が増えないので計算が速いこと。デメリットは $\epsilon_N$ が小さいと貫通が大きくなり、大きすぎると条件数が悪化して収束しなくなること。
ペナルティ剛性の「ちょうどいい値」を探すのが大変そうですね...。拡張ラグランジュ法だとその問題は解決するんですか？
その通り。拡張ラグランジュ法 (Augmented Lagrangian) は、ペナルティ項にラグランジュ乗数 $\lambda_N$ を加えて、反復的に貫通をゼロに追い込む方法だ：
各反復で $\lambda_N$ を更新して：
$$\lambda_N^{(k+1)} = \lambda_N^{(k)} + \epsilon_N \cdot g_N^{(k)}$$こうすると、$\epsilon_N$ が多少小さくても最終的な貫通がほぼゼロになる。Ansys Mechanicalのデフォルトはこの方式だね。ペナルティ法に比べて反復回数は増えるけど、接触圧の精度と$\epsilon_N$に対するロバスト性が格段に上がる。
使い分けの目安はこう：
- ペナルティ法：陽解法の衝突解析、大規模モデルで速度優先のとき
- 拡張ラグランジュ法：静解析・陰解法で接触圧精度が重要なとき
- 純粋ラグランジュ乗数法：貫通ゼロが絶対条件のとき（ただし剛性行列にゼロ対角が入るので直接法ソルバーが必要）
摩擦モデル
摩擦って単にクーロン摩擦の $\mu$ を入れるだけじゃダメなんですか？
基本はクーロン摩擦でOKだけど、数値的には工夫が必要だ。クーロン摩擦は「固着(stick)」と「すべり(slip)」の2状態を持つ：
問題は、固着からすべりへの遷移が不連続であること。これが収束を悪くする。実務では以下のような手法が使われる：
- 接線方向ペナルティ法：固着状態でもわずかなすべりを許容（弾性すべり）
- 指数型摩擦：すべり速度に依存して $\mu$ を滑らかに変化させる
- Coulomb + 粘性項：$\tau = \mu p_N + c \dot{u}_T$ で遷移を滑らかに
鋼同士で $\mu = 0.15$〜$0.3$、アルミ同士で $0.1$〜$0.35$、ゴムと金属で $0.5$〜$1.0$ あたりが目安。ただし実測値が一番信頼できるし、表面粗さ・潤滑油・温度で大きく変わるから注意が必要だ。
摩擦を入れると途端に収束しなくなるって先輩から聞いたんですが、何かコツはありますか？
実務の鉄板テクニックがあるよ：
- まず摩擦なし（$\mu=0$）で接触を安定させる。接触の開閉だけで収束するか確認
- 次に摩擦を徐々に導入する。ステップ1で $\mu=0$、ステップ2で $\mu=0.15$ のようにランプアップ
- ペナルティの接線剛性を法線の0.1倍程度にする。デフォルトが大きすぎることがある
例えばプレス成形のシミュレーションで、金型接触面に最初から $\mu=0.12$ を入れると収束しないのに、$\mu=0$ で成形途中まで流して、その後摩擦を入れるとスムーズに走る、みたいなケースは本当によくある。
収束性の課題と対策
接触解析の収束でよくハマるパターンって何ですか？ エラーメッセージを見てもどこが悪いかわからなくて...
実務で多い「接触が収束しない原因トップ5」を挙げるね：
- 接触面の法線方向が逆 — シェル要素で裏表が反転しているケースが一番多い。ソルバーが「離れている」と判定して接触が効かない
- 初期貫通 — CADのフィレットやクリアランスの関係で、メッシュ段階で既に面が交差している。Abaqusなら
*CONTACT CLEARANCE ASSIGNMENTで調整 - ペナルティ剛性の不適切な値 — 大きすぎると振動（チャタリング）、小さすぎると過大な貫通
- 荷重ステップが粗すぎる — 接触面が一気に開閉するとNewton-Raphson法が追従できない
- 剛体運動 — 接触前に部品がフリーだと剛性行列が特異になる。弱いバネや安定化で拘束するか、初期接触を確保する
チャタリングって具体的にどういう現象ですか？ どう対処すればいいんですか？
チャタリングは、ある節点が「接触→離れる→接触→離れる」を反復のたびに繰り返す現象だ。Newton-Raphson反復で接触状態が毎回切り替わるから、いつまで経っても収束しない。
対策としては：
- 接触安定化 (Contact stabilization)：離れた瞬間に微小なダッシュポット力を加えて急激な状態変化を緩和する
- 時間ステップのカットバック：自動タイムステップを有効にして、チャタリングが検出されたらステップを細かくする
- Over-closure tolerance：微小な貫通を許容する閾値を設定して、わずかな貫通で状態を切り替えないようにする
実務的には、Abaqusの *CONTACT CONTROLS で STABILIZE を使うか、Ansysなら CNOF / PINB を調整するのが定石だよ。
接触解析で「これだけは覚えておけ」ってポイントをまとめるとどうなりますか？
実務で最低限押さえるべきポイントは5つだ：
- 接触は境界非線形問題。線形解析のつもりで設定するとハマる
- surface-to-surface定式化を基本にする。node-to-surfaceは古い・精度が低い
- ペナルティ法は速いが貫通する、拡張ラグランジュ法は正確だが遅い。問題に応じて選ぶ
- 摩擦はまず無しで接触を安定させてから追加する
- 収束しないときは設定ミス（法線方向、初期貫通）を最初に疑う
接触解析は構造CAEの中でもトラブルが多い分野だけど、この5つを理解していれば8割の問題は自力で解決できるよ。
関連用語
- ペナルティ法 — 拘束条件をペナルティ項で近似する手法
- ラグランジュ乗数法 — 拘束条件を厳密に満たすための乗数法
- 摩擦 — 接触面間のすべり抵抗
- 境界非線形 — 境界条件自体が解に依存する非線形性
- Newton-Raphson法 — 非線形FEMの標準的な反復解法
- 接触応力 — 接触面に生じる面圧
- Hertz接触 — 弾性体の接触理論の古典解
- 収束 — 反復解法が解に近づくこと
CAE用語の正確な理解は、チーム内のコミュニケーションの基盤です。 — Project NovaSolverは実務者の学習支援も視野に入れています。
接触解析の実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
お問い合わせ（準備中）関連トピック
なった
詳しく
報告