Dynamic sliding and rolling friction models for linear viscoelastic contact pairs

技術要約：線形粘弾性接触対における動的滑りおよび転がり摩擦モデル

問題提起
サーボシステムから車輪とレールの相互作用に至るまで、工学応用における摩擦モデリングでは、伝統的に滑りと転がりが区別されてきた。滑り摩擦は、通常、有限次元の常微分方程式（ODE）系としてモデル化されるが、これは多くの場合、接触界面において材料の移流が存在しないという仮定（例：剛平面上を滑る軟らかいブロック）に基づいている。対照的に、転がり摩擦は、Kalkerが先駆的に示したように、接触パッチを通じた力と変形の移流を考慮するために偏微分方程式（PDE）を用いて定式化される。

しかし、両方の物体が変形可能であるか、あるいは弾性特性が著しく異なる（タイヤのようなポリマーに共通する）接触対を扱う場合、この二分法は問題となる。このようなシナリオでは、接触する物体のレオロジー特性に応じて材料粒子や内部状態変数が異なる速度で移流するため、滑り摩擦を純粋に集約されたもの（ODEベース）として扱うという仮定が崩れる。既存のモデルは、こうしたダイナミクスを統一したり、粘弾性接触対における運動学的結合を正しく処理したりすることに苦慮することが多い。

手法
本論文は、線形粘弾性レオロジーと非線形動的摩擦モデルを組み合わせた統一的な理論フレームワークを提案している。その手法は以下の通りである。

1. レオロジー・モデリング： 接触界面は「ブリッスル（毛）状」要素を用いてモデル化される。これらの要素の粘弾性挙動は、安定な線形状態空間システムとして表される一般化ケルビン・フォークト（GKV）モデルによって記述される。これにより、幅広い周波数領域にわたる複雑なメモリ効果とヒステリシス効果を記述することが可能となる。
2. 粘弾性・運動学的導出： 著者らは、レオロジーの状態空間表現と運動学的関係を結合させることで、支配方程式のセットを導出している。決定的なのは、標準的な時間微分を、接触領域を通じた内部状態変数（$\zeta$）およびブリッスル変形（$z$）の輸送を考慮するために、物質（オイラー）微分に置き換えている点である。
3. PDE定式化： 構成則と運動学的輸送方程式を統合することにより、著者らは、接線方向の接触力（$f$）および内部状態変数（$\zeta$）の進化を支配する、半線形双曲型偏微分方程式（PDE）系を導出する。
4. 数学的解析： 得られたPDE系は、その双曲型の性質を証明するために解析される。著者らは、このシステムをリーマン変数へと変換し、接触対の弾性分配によって決定される特性方向に沿って方程式をデカップリング（分離）できることを示している。
5. 境界条件および初期条件： 接触領域の「リーディング（前縁）」、「ニュートラル（中間）」、および「トレイリング（後縁）」のエッジに基づいた適切な境界条件（BC）が定式化され、問題の散逸性とウェル・ポーズド性（適切解の存在性）が保証されている。

主な貢献
本研究は、接触力学および摩擦モデリングの分野において、主に以下の3つの貢献を行う。

• 分散フレームワークへの体系的な拡張： 本研究は、動的粘弾性摩擦モデルを集約された（ODE形式の）定式化から、分散型（PDE形式の）フレームワークへと拡張するものである。基材がコンプライアント（柔軟）であるか、あるいは両方の物体が変形可能である場合、滑り接触であっても、支配方程式は自然にPDEとなることを示している。
• 粘弾性・運動学的方程式の導出： 著者らは、滑りおよび転がり接触の両方に対して、摩擦力、材料変形、および内部状態変数の結合した進化を支配する特定のクラスの方程式を導出している。これらの方程式は、2つの物体の異なる輸送速度を明示的に考慮している。
• 滑りと転がりのダイナミクスの統一： 本理論は、滑りと転がりの摩擦が共通の数学的構造を共有していることを明らかにしている。すなわち、両者は線形または半線形双曲型PDE系によって記述される。ODEベースの滑りモデルとPDEベースの転がりモデルの区別は、運動の種類に固有のものではなく、レオロジーに依存するものであることが示された。

結果
提案された理論は、FrBDモデルをGKV要素と結合させた拡張版を用いて、代表的な滑りおよび転がり接触の例を通じて検証されている。

• 滑り接触： 粘弾性基材上を滑るブロックの場合、モデルは基材が剛体である極限（移流がない場合）においてのみ、古典的なODE挙動を回収する。基材が粘弾性である場合、モデルは空間的に変化する摩擦力とブリッスル変形を予測する。結果として、摩擦力の大きさは、両者の剛性比およびレオロジーモデルの散逸分岐の数に依存することが示される。高次のレオロジーモデルは空間的な緩和をもたらし、接触長に沿って摩擦力を減衰させる。
• 転がり接触： 転がり接触の場合、材料の移流が転がり運動に固有であるため、モデルは剛性比に関わらず完全なPDEとしての性質を保持する。解析により、システムの双曲性が維持されていることが確認され、特性速度は主要な弾性相互作用モードによって決定される。
• 過渡現象および定常状態の挙動： シミュレーションにより、システムは粘弾性材料の緩和時間に依存して明確な過渡挙動を示すことが示された。定常状態の滑りにおいて、摩擦力は摩擦係数と垂直荷重の積によって厳密に制限され、剛性比が剛体側に有利な場合にその上限に近づく。

意義と主張
本論文は、そのフレームワークが、広範なクラスの粘弾性摩擦システムに適用可能な一般的な理論的設定を提供すると主張している。主な意義は、滑りと転がりの摩擦モデル間の概念的な曖昧さを解消することにある。両現象が同じクラスの双曲型PDEによって記述できることを示すことにより、本研究は、従来のODE領域とPDE領域への分離は、物理的な本質的な違いではなく、簡略化された仮定（具体的には、剛体基材またはデカップリングされた運動学の仮定）によるアーティファクト（副産物）であることを示唆している。

著者らは、本アプローチが、弾性的に性質の異なるペアにおける材料粒子の運動学を正しく扱うものであることを強調している。これは、従来のブラシ型モデルにおいて、滑りとスピンの定義に関して概念的な困難を生じさせていたシナリオである。本理論は、実用的なアプリケーション（タイヤなど）に見られる粘弾性材料をモデル化するために必要な複雑さを組み込んだ、従来の弾性接触対に関する研究の厳密な拡張として提示されている。結論として、分散型定式化は単なる代替案ではなく、粘弾性接触を正確にモデル化するための必要な記述であり、滑りと転がりのダイナミクスの間の溝を埋める統一的な視点を提供するものである。