[論文レビュー] The contact mechanics challenge: Problem definition

Martin H. Müser, Wolf B. Dapp|arXiv (Cornell University)|Dec 8, 2015

Adhesion, Friction, and Surface Interactions被引用数 12

ひとこと要約

本論文は、自己相似的粗さ、線形弾性、小傾斜近似、および短距離接着を特徴とする接触力学における標準化されたベンチマーク問題を定義する。理論家およびシミュレーターに対し、単なる荷重-接触面積比例を越えて、ギャップおよび応力分布、接触剛性、Tabor係数依存性といった高度な観測量を予測するよう挑戦し、提出後、参照データを共有することで、理論検証のための厳密な比較を可能にする。

ABSTRACT

We present a contact mechanics problem, which we consider to be representative for contacts between nominally flat surfaces. The main ingredients of the mathematically fully defined contact problem are: Self-affine roughness, linear elasticity, the small-slope approximation, and short-range adhesion between the frictionless surfaces. Surface energies, elastic contact modulus and computer-generated surface topographies are provided at www.lms.uni-saarland.de/contact-mechanics-challenge. To minimize the undesirable but frequent problem of unit conversion errors, we provide some benchmark results, such as the relative contact area as a function of load $a_{ m r}(L)$ between $0.1\%$ and $15\%$ relative contact. We call theorists and numericists alike to predict quantities that contain more information than $a_{ m r}(L)$ and provide information on how to submit predictions. Examples for quantities of interest are the mean gap or contact stiffness as a function of load as well as distributions of contact patch size, interfacial stress, and interfacial separation at a reference load. Numerically accurate reference results will be disseminated in subsequent work including an evaluation of the submitted results.

研究の動機と目的

• 名目的に平らで、接着性を有する粗い表面を代表する、標準化され、数学的に定義された接触力学問題を確立すること。
• 従来の理論-シミュレーション比較に起因する制限を克服し、単に荷重-接触面積比例をテストするのではなく、相対接触面積を超えるより豊かな観測量、例えばギャップおよび応力分布、接触剛性、Tabor係数依存性の予測を促すこと。
• 単位換算誤差を最小限に抑えるために、参考パラメータおよびベンチマーク結果を提供し、0.1% から 15% の接触における $ a_{\rm r}(L) $ を含む。
• 高精度なシミュレーションと理論モデルの厳密な検証を促進するため、構造的な提出プロトコルと共同著者機会を提供すること。

提案手法

• 表面の粗さは、Hurst指数 $ H $、ロールオーバー波長 $ \lambda_{\rm r} = 20\ \mu\text{m} $、短波長カットオフ $ \lambda_{\rm s} = 0.1\ \mu\text{m} $ で定義される自己相似的高さスペクトルを用いたフーリエ法により生成される。
• 線形弾性、摩擦なしの接触、およびハードウォール反発モデルを用いることで、接触点および界面分離の明確な定義がなされる。
• 弾性接触モジュラスと表面エネルギーは固定され、Tabor係数 $ \mu_{\rm T} $ は接着遷移を調べるための調整可能なパラメータとして用いられる。
• 周期的シミュレーションセル（サイズ $ 0.1\ \text{mm} \times 0.1\ \text{mm} $）を用い、高さプロファイルは平均二乗勾配 $ \bar{g} = 1 $ に正規化される。
• 短距離接着はLennard-Jonesに類似したポテンシャルを用いてモデル化され、空間分解能を細かくすることで連続体極限に近づける数値的解法が採用される。
• 参加者は、ギャップ、界面応力、接触パッチサイズの分布を含む9つの特定の観測量について予測を提出することが招待される。特に漸近的挙動が予想される場合には対数スケールが推奨される。

実験結果

リサーチクエスチョン

• RQ1 接着性を有する自己相似的粗い接触において、荷重が増加するに従い、接触パッチサイズの分布はどのように変化するか？
• RQ2 高荷重域、完全接触に近づく領域において、相対接触面積 $ a_{\rm r}(p) $ の関数的形態はどのようなものか？
• RQ3 Tabor係数 $ \mu_{\rm T} $ に応じて、ギャップおよび界面応力分布はどのように変化するか。特に平均ギャップの局所的極値付近での挙動に注目する。
• RQ4 平均ギャップまたは法線方向変位が $ \mu_{\rm T} $ の関数としてどのように振る舞うか。その極値はどこに現れるか？
• RQ5 例えば $ x = 50\ \mu\text{m} $ における断面に沿った応力およびギャッププロファイルは、微視的シミュレーションと均質化モデルとでどのように比較できるか？

主な発見

• 0.1% から 15% の接触範囲における相対接触面積 $ a_{\rm r}(L) $ がベンチマークとして提供され、提出における単位換算誤差が最小限に抑えられる。
• Tabor係数 $ \mu_{\rm T} $ に対する平均ギャップは局所的極値を示し、接着と弾性変形の相互作用に起因する非単調応答を示している。
• 接触パッチサイズ、ギャップ、界面応力の分布は、高精度なシミュレーションと理論の比較に不可欠な主要な観測量として特定された。
• 界面応力スペクトルおよび応力分布関数は、高応力域でべき乗則に類似した挙動を示し、明確化のため対数スケールの使用が推奨される。
• 荷重に応じた接触剛性および平均ギャップは、小さな分離領域で漸近的挙動を示し、対数軸を用いることで最もよく解明される。
• 提出プロセスには、成功した提出者に対して共同著者機会が含まれ、データは *Tribology Letters* に共有され、公開される。