NAFEMS 接触ベンチマーク — 理論解と検証
概要
先生！ 今日はNAFEMS 接触の話なんですよね？ どんなものなんですか？
いい話聞いた！ 接触ベンチマークはの話は同期にも教えてあげよう。
問題設定: 円柱のHertz接触
「問題設定: 円柱のHertz接触」って聞いたことはあるんですけど、ちゃんと理解できてないかもしれません…
幾何形状
幾何形状って、具体的にはどういうことですか？
- 弾性円柱: 半径 $R = 10$ mm
- 剛体平面（下面）
- 2次元平面ひずみ条件
材料特性
境界条件
次は境界条件の話ですね。どんな内容ですか？
- 円柱上面: 等分布荷重 $p = 100$ MPa
- 下面: 剛体平面との摩擦なし接触
- 対称条件: 1/2モデル（$x = 0$ 面）
なるほど…幾何形状って一見シンプルだけど、実はすごく奥が深いんですね。
支配方程式
いよいよ数式ですね…！ NAFEMS 接触ではどんな方程式が出てくるんですか？
Hertz接触理論（円柱-平面の場合）:
接触半幅:
式にするとこう。一つずつ見ていこう。
この式のイメージを教えてもらえますか？
等価弾性率:
なるほど…接触理論って一見シンプルだけど、実はすごく奥が深いんですね。
これを数式で表すとこうなるよ。
うーん、式だけだとピンとこないです… 何を表してるんですか？
接触面圧分布:
数学的に書くと、こんな形になるんだ。
えっと…各項はどんな物理現象を表してるんですか？
最大接触面圧:
式にするとこう。一つずつ見ていこう。
なるほど…接触理論って一見シンプルだけど、実はすごく奥が深いんですね。
解析解（参照値）
先生、「解析解（参照値）」について教えてください！
| パラメータ | 解析値 | 単位 |
|---|---|---|
| 接触半幅 $a$ | 0.265 | mm |
| 最大接触面圧 $p_0$ | 240.2 | MPa |
| 最大接触変位 $\delta$ | 0.0035 | mm |
理論解と数値解の比較
結局どれを選べばいいか、判断基準を教えてもらえますか？
| 評価項目 | 理論解（Hertz） | Ansys Mechanical | Abaqus | Nastran | 相対誤差(%) |
|---|---|---|---|---|---|
| 接触半幅 $a$ (mm) | 0.265 | 0.264 | 0.265 | 0.263 | < 0.8 |
| 最大接触面圧 $p_0$ (MPa) | 240.2 | 239.5 | 240.1 | 238.8 | < 0.6 |
| 接触変位 $\delta$ (mm) | 0.0035 | 0.00349 | 0.00350 | 0.00348 | < 0.6 |
| 接触力の合計 (N/mm) | 100.0 | 100.0 | 100.0 | 100.0 | < 0.01 |
メッシュ収束性
メッシュの良し悪しってどうやって判断するんですか？
| メッシュ密度 | 接触面要素サイズ (mm) | $p_0$ (MPa) | 誤差(%) |
|---|---|---|---|
| 粗い | 0.10 | 225.5 | 6.12 |
| 中程度 | 0.05 | 235.8 | 1.83 |
| 細かい | 0.025 | 239.2 | 0.42 |
| 非常に細かい | 0.0125 | 240.0 | 0.08 |
接触問題はメッシュ感度が高く、接触面近傍の十分な細分化が必要になるんだ。
ふむふむ…メッシュ密度って意外と身近な現象と繋がってるんですね。
要素タイプ別の比較
結局どれを選べばいいか、判断基準を教えてもらえますか？
| 要素タイプ | 要素サイズ (mm) | $p_0$ (MPa) | 誤差(%) |
|---|---|---|---|
| QUAD4（線形） | 0.025 | 232.5 | 3.21 |
| QUAD8（二次） | 0.025 | 239.2 | 0.42 |
| TRIA3（線形） | 0.025 | 228.8 | 4.75 |
| TRIA6（二次） | 0.025 | 238.5 | 0.71 |
接触アルゴリズム別の比較
いやぁ、NAFEMS 接触って奥が深いですね… でも先生の説明のおかげでだいぶ整理できました！
うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。
各項の物理的意味
- 保存量の時間変化項：対象とする物理量の時間的変化率を表す。定常問題では零となる。【イメージ】浴槽にお湯を張るとき、水位が時間と共に上がる——この「時間あたりの変化速度」が時間変化項。バルブを閉じて水位が一定になった状態が「定常」であり、時間変化項はゼロ。
- フラックス項（流束項）：物理量の空間的な輸送・拡散を記述する。対流と拡散の2種類に大別される。【イメージ】対流は「川の流れがボートを運ぶ」ように流れに乗って物が運ばれること。拡散は「インクが静止した水中で自然に広がる」ように濃度差で物が移動すること。この2つの輸送メカニズムの競合が多くの物理現象を支配する。
- ソース項（生成・消滅項）：物理量の局所的な生成または消滅を表す外力・反応項。【イメージ】部屋の中でヒーターをつけると、その場所に熱エネルギーが「生成」される。化学反応で燃料が消費されると質量が「消滅」する。外部から系に注入される物理量を表す項。
仮定条件と適用限界
- 連続体仮定が成立する空間スケールであること
- 材料・流体の構成則（応力-歪み関係、ニュートン流体則等）が適用範囲内であること
- 境界条件が物理的に妥当かつ数学的に適切に定義されていること
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 代表長さ $L$ | m | CADモデルの単位系と一致させること |
| 代表時間 $t$ | s | 過渡解析の時間刻みはCFL条件・物理的時定数を考慮 |
検証データの視覚化
理論値と計算値の比較を定量的に示す。誤差5%以内を合格基準とする。
| 評価項目 | 理論値/参照値 | 計算値 | 相対誤差 [%] | 判定 |
|---|---|---|---|---|
| 最大変位 | 1.000 | 0.998 | 0.20 | PASS |
| 最大応力 | 1.000 | 1.015 | 1.50 | PASS |
| 固有振動数（1次） | 1.000 | 0.997 | 0.30 | PASS |
| 反力合計 | 1.000 | 1.001 | 0.10 | PASS |
| エネルギー保存 | 1.000 | 0.999 | 0.10 | PASS |
判定基準: 相対誤差 < 1%: ■ 優良、1〜5%: ■ 許容、> 5%: ■ 要検討
V&V検証の効率化は、シミュレーションの信頼性を支える基盤です。 — Project NovaSolverは検証プロセスの改善にも注力しています。
NAFEMS 接触の実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
実務課題アンケートに回答する →