ペナルティ法による接触定式化 — トラブルシューティングガイド
この記事は統合版に移行しました
より充実した内容を penalty-contact.html でご覧いただけます。
より充実した内容を penalty-contact.html でご覧いただけます。
ペナルティ法による接触定式化 — トラブルシューティングガイド
接触解析のトラブル
接触解析でよくあるトラブルを教えてください。
収束しない
接触解析の収束困難は最も一般的なトラブル。対策：
1. 初期増分を小さく（0.1→0.01→0.001）
2. 接触安定化を使用（*CONTACT STABILIZATION）
3. 初期ギャップを除去（最初から接触状態に）
4. 摩擦を段階的に導入（$\mu=0$→$\mu=0.1$→$\mu=0.3$）
5. ペナルティ剛性を下げる（貫通は増えるが収束しやすい）
貫通量が大きい
ペナルティ剛性が小さすぎる。$k_p$を上げるか、Augmented Lagrangian法に切り替え。
接触圧にチェッカーボードパターン
1次要素（TET4, HEX8完全積分）の接触面で圧力が振動する。C3D10M（Abaqus）に切り替えるか、二次要素の低減積分を使う。
スレーブがマスターに食い込む
マスター/スレーブの選択が逆。「硬い側＝マスター」の原則を確認。
まとめ
Coffee Break よもやま話
エネルギー誤差と砂時計
ペナルティ接触解析での典型的な失敗として「hourglass（砂時計）エネルギーの過大」がある。低次縮減積分要素とペナルティ接触を組み合わせると、接触力がhourglass変形を刺激し数値的に発散するケースがある。2008年のLivermore Software社のベンチマーク報告では、ERODED_SINGLE_SURFACE接触でhourglassエネルギーが全変形エネルギーの15%を超えると結果信頼性が著しく低下することが示され、hourglassコントロールの強化か接触定式化の変更が勧告された。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ペナルティ法による接触定式化の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
この記事の評価
ご回答ありがとうございます！
参考に
なった
なった
もっと
詳しく
詳しく
誤りを
報告
報告