接触解析が難しい理由と収束させるための実践対策

接触解析が難しい本当の理由と、収束させるための実践的な対策

接触解析の設定を正しく行わないと、計算時間が弾性体同士の解析の約4倍に膨れ上がる場合があります。

接触解析が難しいと言われる3つの根本的な理由

金属部品同士の組み合わせを扱う設計者の多くが、最初に壁にぶつかるのが接触解析です。「設定してみたが計算が終わらない」「エラーが出て何が悪いのかわからない」という声は現場でも絶えません。そもそも、なぜ接触解析はこれほど難しいとされているのでしょうか。

まず根本的に理解しておきたいのは、接触現象そのものが「状態が急激に変わる非線形現象」であるという点です。通常の線形解析では変形量と応力が一定の比例関係を持ちますが、接触解析ではパーツ同士の接触・分離・滑りといった状態が解析の進捗に応じて刻一刻と変化します。この「急激な状況変化」が計算を不安定にさせる根源です。

具体的には、主に3つの問題が絡み合っています。

1つ目は「接触剛性の設定問題」です。接触面には仮想的なバネ（接触剛性）が張られ、物体同士のすり抜けを防ぎます。このバネ剛性が高すぎると計算が発散し、低すぎると精度が落ちるという、絶妙なバランスを求められます。デフォルト設定のまま流して失敗するのは、まさにここが原因です。

2つ目は「摩擦係数の扱いの難しさ」です。実際の金属加工現場では鋼材と治具の間の摩擦係数が不明なことが多く、数値を仮定するしかありません。しかも摩擦係数が大きすぎると解の収束がさらに困難になるため、初期段階では摩擦係数を0（完全すべり）で試してから段階的に現実の値に近づけるアプローチが有効です。

3つ目は「多体接触におけるすべての接触面の定義」です。アセンブリモデルでは接触する可能性のある面を漏れなく定義しなければならず、定義の抜けがあると計算中に部品がすり抜けてしまいます。近年はCAEソフトが自動抽出してくれる機能を持つものも増えましたが、複雑な金属部品の組み合わせでは人間の目での確認も欠かせません。

難しいということですね。しかしこれらの原因を一つひとつ把握していれば、対策を立てることは十分に可能です。

SBDテクノロジー：有限要素法の脇役シリーズ Vol.13 – 厄介な接触問題（接触解析の本質的な難しさをわかりやすく解説）

接触解析の収束を妨げる「剛体移動」の正体と対処法

接触解析の現場で最も多い失敗のひとつが「剛体移動による発散」です。これは金属加工の設計者にとっても非常に身近なトラブルです。

剛体移動とは何でしょうか？ 静解析において、固定されていない物体に力を加えると、その物体は無限遠に飛んでいってしまいます。これが剛体移動です。通常の解析であれば拘束条件を設けることで防げますが、接触解析では「接触そのもので部品を支える」設定にするため、接触が成立するまでの瞬間に部品が浮いた状態になり、剛体移動が起きやすいのです。

たとえば、プレスブレーキで曲げ加工するときの型と板材の関係を思い浮かべてください。板材は型に押しつけられる前、何にも支えられていない宙ぶらりんの状態です。その状態に荷重をかけると、計算上は板材が無限に移動してしまいます。これがまさに剛体移動による発散です。

特にやっかいなのが「初期ギャップ」の問題です。CADモデル上では接触しているように見えるのに、メッシュを切ると有限要素モデル上ではわずかにギャップが生じているケースがあります。このギャップがあると計算開始時から接触が認識されず、剛体移動が発生してしまいます。1/10mm以下のギャップでも計算は平気で崩壊します。

対処法は複数あります。まず推奨されるのは、CAEソフトの「接触間調整」機能（Ansys Mechanicalなら「ギャップ食い込みをなくす」設定）を使って接触判定位置を自動オフセットする方法です。ただしオフセット量が大きすぎると形状と解析結果に矛盾が生じるため、あくまで小さなギャップへの対処として使います。

もう一つの有効な方法が、力荷重を強制変位に切り替える手法です。強制変位を使うと物体の位置が常に拘束されるため、剛体移動が原理的に発生しません。接触設定が正しくできているか確認するための予備解析としても最適な方法です。設計者が最初の検証段階で使うべきテクニックです。

収束しない場合の最終手段として、弱いバネで部品を支える方法もあります。ただしバネ剛性が高すぎると解析結果に影響が出るため、バネ端部の反力を必ず確認することが条件です。

サイバネットシステム：接触解析のトラブルシューティング（剛体移動・接触剛性・エッジ接触の対処法を詳解）

接触解析で絶対に確認すべきメッシュと境界条件の設定ポイント

接触解析を失敗させる原因として、剛体移動と同じくらい多いのがメッシュ品質と境界条件の設定ミスです。ここを正しく理解しておくかどうかで、解析の成功率が大きく変わります。

まずメッシュについてです。接触面のメッシュが粗すぎると、部品同士の接触認識が不安定になります。特に注意が必要なのは、角（エッジ）や頂点が当たる接触です。Ansys Mechanicalなどのツールでは、接触判定点はデフォルトで「ガウス点」と呼ばれるコーナー節点の内側にあります。そのためエッジが当たる場合、ガウス点では認識できずに過大な食い込みが発生し、計算が発散します。

これはいいことですね……と言いたいところですが、実際には見落としやすい落とし穴です。

対策として、エッジや頂点同士が当たる接触では判定点を「節点」に切り替える設定が有効です。ただし節点に変更すると判定点の数が減少するため（例：隣接2要素の場合、ガウス点8点→節点6点に減少）、精度面で注意が必要です。可能であればフィレットをつけて面対面での接触に形状変更することを優先すべきです。

次に境界条件の設定についてです。接触解析では「接触タイプ」の選択も重要なポイントです。

| 接触タイプ | 特徴 | 使用場面の例 |
|---|---|---|
| 固着（Bonded） | 相対変位なし、線形解析可 | 溶接・接着部位 |
| 摩擦なし（Frictionless） | 法線方向のみ力伝達 | 軸受面・スライド面 |
| 摩擦あり（Frictional） | 静摩擦・動摩擦を考慮 | ボルト締結・プレス加工 |
| ラフ接触（Rough） | 完全固着・摩擦係数相当∞ | 完全に滑らない接触面 |

摩擦係数の大きな接触（たとえば0.4以上）では収束がさらに困難になります。摩擦が大きい金属接触面を解析する際には、まず摩擦なし設定で収束するかを確認し、収束できたら段階的に摩擦係数を上げていく手順が安全です。

また、境界条件として「固定支持」を誤って設定すると、部品が動かないよう完全に拘束されてしまい、接触による変形挙動が正しく評価できなくなります。回転や摺動を伴う金属部品の解析では、特に自由度の設定に気を配る必要があります。自由度の設定が条件です。

サブステップ数も見逃せません。接触解析は状態変化が急激なため、サブステップが不足すると接触が認識される前に部品がすり抜けてしまいます。落下・衝突・プレスの初期段階では最小サブステップ数を意識的に細かく設定することが不可欠です。

MONO塾：CAE解析によくある失敗事例（メッシュ品質・境界条件ミスの具体例と改善方法）

接触解析の計算時間を大幅に短縮する「剛体モデル活用」という視点

ここからは、上位記事にはあまり書かれていない独自の視点をお伝えします。「接触解析は難しい」という文脈で語られることが少ないのが、「剛体モデルを使うと計算コストが約1/4に削減できる」という事実です。

接触解析では、相互作用する2つの部品の両方を弾性体として扱うのが基本です。しかし、剛性差が大きい2つの部品が接触する場合、剛性の高い方を「剛体」として定義することで、計算コストを大幅に削減できます。

具体的なイメージとしては、圧延ロール（鋼製）と加工材料（軟鋼・アルミ）の組み合わせが典型例です。実際の圧延解析では、剛性が格段に高いロール側を剛体と定義することで、弾性体同士の解析と比べて計算コストが約1/4に低減されたというデータがあります。時間に換算すると、4時間かかっていた計算が1時間で終わるイメージです。

これは使えそうです。

剛体モデルを使うメリットはコスト削減だけではありません。剛体ボディ側には変形が発生しないため、接触剛性マトリクスがシンプルになり、収束性も向上します。「難しい」と言われる接触解析の難易度そのものが下がるわけです。

注意点もあります。剛体として設定したボディには変形量や応力といった結果を評価できません。つまり、工具や型側の強度を評価したい場合には剛体化できないということです。評価対象の部品は弾性体のまま残し、評価不要な治具・型・ロールなどを剛体化する設計が基本です。

大規模なアセンブリモデルでも同様の発想が使えます。応力評価が必要なパーツ以外をすべて剛体化してしまうことで、部品点数が多いモデルでも現実的な計算時間で解析を完了できます。オーバーヘッドバルブエンジンの動的解析でカム以外の全部品を剛体化し、計算時間を大幅に短縮した事例もあります。

まず剛体モデルを試す、というアプローチを接触解析の第一段階に組み込むことで、「収束しない→時間だけが過ぎる」という最悪のパターンを避けられます。

サイバネットシステム：剛体-弾性体接触解析（圧延加工解析での計算コスト1/4削減の事例を紹介）

接触解析の結果を正しく読むためのエラーメッセージ活用法

計算が発散したとき、多くの設計者は「また収束しなかった」と途方に暮れます。しかし実は、ソルバーが出力するエラーメッセージには解決のヒントが詰まっています。エラーを正しく読む習慣を持つだけで、対処時間を数時間単位で節約できます。

代表的なエラーメッセージとその意味を整理します。

- 「DOF Limit Exceeded」：自由度の値が上限を超えた。ほぼ剛体移動が原因です。初期接触状態を確認し、拘束条件の見直しや固定忘れがないかをチェックします。

- 「Error in Element Formulation」：要素が極短時間に過大な変形をした場合に発生します。物性値・境界条件・単位系の設定ミスを疑います。単位をmmとNで設定すべきところをm/kNで入力していた、というケースは意外に多いです。

- 「Unconverged Solution」：ニュートン-ラプソン法の残差が収束しなかった。接触部位での残差が大きい場合は接触剛性（垂直剛性係数）の見直しで改善できることが多いです。

- 「BEGIN BISECTION NUMBER ○○」：計算中のバイセクション（サブステップを半分にして収束を試みる動作）の記録です。多発している場合は初期サブステップ数が不足しています。

エラーメッセージは英語で出力されることがほとんどですが、発生タイミングと前後のログを合わせて読むと原因箇所が絞り込めます。発散直前のサブステップの結果を保存しておく設定（Ansys Mechanicalなら途中結果出力を有効化）をしておくと、発散した時点の変位状態を目視確認でき、原因特定が格段に楽になります。

収束判定グラフも活用しましょう。ニュートン-ラプソン残差グラフで、どのステップで残差が下がらなくなっているかを視覚的に確認できます。接触開始直後に発散する場合は接触剛性、解析途中で発散する場合は剛体移動や摩擦の影響が疑われます。

また、解析前に「初期接触診断ツール」を使うことを強く推奨します。Ansys Mechanicalでは初期接触状態を可視化できるツールが実装されており、接触すべき箇所が接触として認識されているか、ギャップや食い込みが存在するかを解析実行前に確認できます。これを習慣にするだけで、無駄な発散を大幅に減らせます。初期確認が原則です。

解析後も確認は必要です。食い込み量コンターを出力し、接触面で発生している食い込みが許容範囲内かをチェックします。食い込みが大きすぎる場合は接触剛性ファクターを上げる方向で精度を改善します。

接触解析は「エラーが出て終わり」ではなく、「エラーを読んで原因を特定し、次の改善策を打つ」というサイクルを繰り返す解析です。エラーメッセージは敵ではなく、収束へのナビゲーションです。それが大事です。

村田ソフトウェア（Femtetヘルプ）：接触解析が収束しない場合の対処法（設定パラメータと原因ごとの対策を解説）

商品名