ナノスケールからマクロスケールへ：法線荷重下での粗面接触面積の進化

概要

粗面の法線荷重に対する接触面積の進化は、地震ダイナミクスから機械摩耗まで、基礎的かつ実用的に極めて重要である。本研究は、原子スケールとマクロスケールの間のギャップを、法線接触挙動で橋渡しする。実接触面積は、多数の離散的接触（クラスタ）アンサンブルで構成され、見かけ表面積よりはるかに小さいことが証明される。離散的接触クラスタの分布と相互作用が、接触固体メカニズム解明の鍵である。これを目的として、グリーン関数分子動力学（GFMD）を用い、接触クラスタが原子スケールからマクロスケールへどう進化するか、およびクラスタ間相互作用を研究した。クラスタ間相互作用がそれら形成に強い影響を持つことが判明。接触クラスタの形成と分布は、アスペリティモデル予測より遥かに複雑。相互作用無視は接触力過大評価につながる。実際の接触では、アスペリティ間相互作用により接触クラスタはより小さく、より離散的になる。法線荷重に対する正確な接触面積の性質理解は、摩擦研究への基盤必須。

背景

ほとんどのマクロスコピック表面は粗くフラクタルと見なされる。粗面間接触挙動は完全平滑面より遥かに複雑。実接触面積は多数の離散的接触領域（クラスタ）アンサンブルで形成され、見かけ表面積よりはるかに小さい。法線力と接触クラスタのサイズ、形状、分布が接触挙動解明鍵であり、摩擦研究へ必須。

1960年代以来、接触面積と荷重関係を求める多数モデル提案されている。中でもアスペリティモデルが最もシンプルで普及。初期応用でGreenwoodとWilliamsonは、アスペリティ同一半径・異なる高さ仮定で接触界面粗さ記述。以来アスペリティモデル主流、膨大文献蓄積。WhitehouseとArchardはアスペリティ先端曲率半径ランダム性考慮でG-Wモデル発展。Nayakはガウス粗さ解析にランダム過程理論導入、Bushらが粗面接触に適用。

アスペリティモデル基本仮定の一つは、アスペリティ間相互作用無視可能。これにより潜在接触アスペリティが表面幾何学から事前決定可能とされる。しかしこの仮定は接触力と接触面積の不正確推定招く可能性。接触クラスタ進化と相互作用取得のため、グリーン関数分子動力学（GFMD）用いフラクタル粗面研究。

本研究は原子スケールからマクロスケールへの法線接触挙動ギャップ橋渡し目的。アスペリティ相互作用考慮し、数値例で接触面積進化を原子からマクロスケールへ実証。以下、フラクタル表面生成、GFMDモデル、接触クラスタ検出アルゴリズム、数値実験設計手法簡述後、接触クラスタ形成・発展と界面挙動へ影響焦点化。

方法

フラクタル粗面生成

粗面接触挙動数値モデル研究のため表面生成必要。フラクタル表面生成アルゴリズム複数存在。本研究ではフーリエ変換法採用（図1参照）。フラクタル粗面幾何決定に4パラメータ必要：最大周波数（w_H）、最小周波数（w_L）、ハースト指数（H）、振幅標準偏差（P）。表面基本統計パラメータ、RMS高さ（√M₀）、RMS傾斜（√M₂）、RMS曲率（√M₄）が界面挙動鍵、M_iは表面スペクトル第iラジアルモーメント。表面統計パラメータM_iはプロファイル統計パラメータm_iと以下関係：M₀=m₀、M₂=2m₂、M₄=4/3m₄。アスペリティ密度n（表面山/谷）は以下式で決定：

n = 1/(6π√3) · (m₄/m₂) （式1）

また、表面総山/谷数Nは：

N = A₀ × n = A₀/(6π√3) · (m₄/m₂) （式2）

ここでA₀は見かけ表面積。自己アフィンフラクタル表面では、統計パラメータは入力パラメータ（w、H、P）と以下関係：

m_i = ∫_{w_L}^{w_H} ω^i Φ_φ(ω) dω = ∫_{w_L}^{w_H} ω^i B ω^{-(1+2H)} dω （式3）

Bは表面粗さ定数、Pと関係。式（3）と（2）よりフラクタル表面山/谷数は波長とハースト指数依存。フラクタル表面統計特性詳細議論は文献参照。

フーリエ変換アルゴリズム典型例として、ハースト成分H=0.5、最大周波数w_L=1/(24σ)、最小周波数w_H=1/(256σ)、振幅標準偏差P=0.69、系サイズ512×512原子（初期間隔1.12σ）設定。これら入力から、RMS傾斜√M₂=0.077、RMS曲率√M₄=0.0077の表面生成。式（2）で総山/谷数150、数値カウントで山158、谷159、誤差5%以内、系サイズ統計的に許容。実際2048×2048原子まで拡大しても統計パラメータは小系と整合。

GFMDモデル

粒子間相互作用実験的捕捉極めて困難。最近分子動力学で粒子間相互作用シミュレート、接触/摩擦メカニズム分子起源解明目指す。しかし大規模MDシミュレーション計算コスト極大。したがってGFMD導入、高効率で表面シミュレート。GFMDは界面原子（ここでは2層）相互作用をMDシミュレート、非界面層（通常弾性挙動）をグリーン関数シミュレート。巨大原子系を界面2層原子へ縮約（図1参照）、計算コスト劇的削減。GFMD詳細議論は文献参照。GFMDモデルでレナード-ジョーンズ（LJ）ポテンシャル使用：

u(r) = 4ε[(σ/r)¹² − (σ/r)⁶] （式4）

εはポテンシャル井戸深さ、σは粒子間ポテンシャル零距離、rは粒子間距離。ε、σ、ε/σをエネルギー、距離、力単位とする。LJポテンシャルより、r=2¹ᐟ⁶σ≈1.12σで粒子間力零。r>1.12σで引力、r<1.12σで斥力。本研究接着考慮せず、カットオフ距離1.12σ設定。原子層結晶構造は面心立方（FCC）。対称性利用し界面層のみで表面幾何形成（図1参照）、平坦面下弾性ブロックをグリーン関数シミュレート。

接触クラスタ認識法

図1に見られる3スケール存在：（1）原子スケール、LJポテンシャルシミュレート；（2）アスペリティスケール、接触原子群効果；（3）マクロスケール、接触クラスタ群効果。接触クラスタサイズ、形状、位置、分布が分子挙動と界面特性の重要ブリッジ。ナノスケールで原子接触領域定義困難。ここでは接触原子を法線成分力f_z>0で定義。連結接触原子を接触クラスタ定義。ラベリング技法使用し接触クラスタ探索。高速化のため再帰探索回避修正アルゴリズム採用。アルゴリズムフローチャート図2、GFMDシミュレートから原子力データ抽出。アルゴリズム8キーステップ：

ステップ1. 行探索開始、新原子データ取得（行ごと原子探索）。
ステップ2. 原子接触判定。非接触ならステップ1戻る。接触なら次ステップ。
ステップ3. 同行前原子と比較。前原子も接触なら前原子所属クラスタへマージ、同番号ラベル。非接触なら新番号（前番号+1）ラベル。
ステップ4. 最終原子判定。非最終ならステップ1、最終なら列探索へ。
ステップ5. 列探索開始、新原子データ取得（列ごと原子探索）。
ステップ6. 原子接触判定。非接触ならステップ5、接触なら次。
ステップ7. 同列前原子と比較。前原子も接触かつ異クラスタ所属なら現クラスタを前原子所属クラスタへマージ、同番号ラベルし保存。非接触または同クラスタなら次。
ステップ8. 現在原子最終判定。非最終ならステップ5、最終なら探索完了。

数値実験設計

二粗面接触問題は、等価弾性率E*導入により、一複合剛性粗面と一平坦弾性面問題に簡略化可能：

1/E* = (1−ν₁²)/E₁ + (1−ν₂²)/E₂ （式5）

E₁、E₂は上下面弾性率。簡略化のため、剛性粗面と弾性平滑面接触とし、接触クラスタ形成・発展と力-面積挙動研究。以下、上記生成表面使用（上面剛性粗面E₁=∞、下面平滑弾性E₂=3ε/σ³）、ν₁=ν₂=0.5。系サイズ512×512原子（初期間隔1.12σ）、x-y平面周期境界。弾性ブロック深さ1024原子層（初期間隔1.12σ）。通常MDなら268,697,600原子；GFMDで524,288原子（2層）へ削減（図1参照）。剛性粗面（上）を平坦弾性面へ徐々押込み。剛性面荷重は変位制御。各変位ステップ0.01σ、GFMDアルゴリズムで原子位置更新、原子力収束基準L¹ノルム=0.01ε/σまで。最大反復50,000設定、無限ループ防止。

結果と考察

接触クラスタ分布と発展

アスペリティモデルはアスペリティ球形/楕円形仮定、接触アスペリティ間相互作用不考慮。本研究では上記生成表面からアスペリティモデル用アスペリティ抽出。アスペリティモデルでは潜在接触アスペリティが表面幾何から事前決定可能、表面山/谷が高さにより接触クラスタ形成。しかし現実ではアスペリティ形状不規則、通常複数隣接アスペリティが大きな一つへマージ（図2参照）。最初6独立アスペリティ観察、接触力増加で最終的に一つ大接触クラスタへマージ（図3）。アスペリティ間距離十分遠く相互影響なし仮定が不正確結果導く可能性示唆。

図4でクラスタ数が接触面積増加に伴い先増加後減少、表面アスペリティ数は常に増加。これは図3説明マージ効果起因。

クラスタマージ効果はアスペリティモデルとGFMDモデル両方で観察。しかし同接触面積でGFMDモデルクラスタ数がアスペリティモデルより遥か多い（図5）。GFMDクラスタ数ほぼアスペリティモデル2倍（図5）。主理由はアスペリティモデルがアスペリティ間相互作用不考慮。GFMDでは接触クラスタ相互影響。クラスタ生成変位場が表面全域連続。大アンサンブルクラスタ変位が弾性面上に新幾何生成、新接触クラスタ形成影響。したがって接触クラスタ形成は剛性粗面高さ基準のみでなく、平滑弾性面変形にも影響される。図6でアスペリティモデルとGFMDモデル各接触面積でのクラスタ分布確認可能。接触面積5%でアスペリティモデル17、GFMD34クラスタ；10%で24と52。GFMDクラスタより離散的。GFMDで平均クラスタサイズ小さいが、ほとんど山/谷と一致（図7）。アスペリティモデルは山または谷どちらか潜在アスペリティ（接触側依存）。しかし図8で接触面積増加に伴い山谷両方接触クラスタ形成可能。小接触面積では谷が大半、10%以上で山も増加形成。

接触面積-荷重関係

法線荷重下力-面積関係は接触挙動必須。先行モデルではアスペリティ通常球形/楕円形仮定。しかし実際クラスタ遥か複雑。ここでは3モデル力-面積関係比較：（1）GFMDモデル；（2）アスペリティモデル（AM）、GFMD使用前表面から直接抽出アスペリティ平坦弾性面押込み（予期せぬクラスタ形成回避）；（3）Greenwood-Williamsonモデル（G-W）、アスペリティ等価球変換。球半径：

1/R = 8/3 (m₄/π)¹ᐟ² （式6）

GFMDと抽出アスペリティモデルで界面総力はGFMD抽出各クラスタ力合計。G-Wモデルでは各アスペリティ力にヘルツ理論使用（GFMD同材料特性）、総力F：

F = Σ f_i = Σ (4/3)E*R¹ᐟ²(d−z_i)³ᐟ² （式7）

z_iはアスペリティ高さ、dは剛性面印加変位、f_iはヘルツ理論基づくアスペリティ接触力。

図9で3モデル力-面積関係比較、線形関係示す。GFMD総力がアスペリティモデルとG-Wモデルより遥か小さい。アスペリティモデルFはGFMD予測1.80倍、G-Wモデル1.54倍。これは接触クラスタRMS傾斜で説明可能。法線荷重はRMS傾斜比例、L∝√M₂。GFMDでは多数クラスタ構成、貫入浅く、アスペリティモデルよりクラスタRMS傾斜小さい。図10で3モデルクラスタRMS傾斜比較。GFMDクラスタRMS傾斜は表面RMS傾斜0.077未満、他2モデルクラスタRMS傾斜は表面RMS傾斜超過。

結論

原子スケールからマクロスケールへの接触面積進化探求のため、GFMDモデルでフラクタル粗面接触問題研究。原子接触を力>0存在で定義、同一系で3長さスケール研究。アスペリティ間相互作用が接触クラスタ形成必須と判明。一部クラスタ十分近接し大クラスタへマージ。実接触領域は弾性平滑面変形により表面幾何予測より遥か複雑。大半クラスタ位置は表面山/谷と一致。しかしクラスタサイズ小さく、形成は表面アスペリティ高さ決定せず。接触面積増加で山谷両方クラスタ形成可能。GFMDでは力遥か小さいがクラスタ数遥か多い。GFMDクラスタRMS傾斜アスペリティモデルより小さい、アスペリティモデル高圧力予測理由説明。実接触面積は表面幾何単純予測不可。法線荷重下実際接触面積は摩擦後続研究へ重要。