Dynamic sliding and rolling friction models for linear viscoelastic contact pairs

기술 요약: 선형 점탄성 접촉 쌍을 위한 동적 미끄럼 및 구름 마찰 모델

문제 정의
서보 시스템에서 휠-레일 상호작용에 이르기까지 공학 응용 분야에서의 마찰 모델링은 전통적으로 미끄럼(sliding)과 구름(rolling) 영역을 구분해 왔다. 미끄럼 마찰은 일반적으로 유한 차원 상미분 방정식(ODE) 시스템으로 모델링되는데, 이는 접촉 계면에서 재료의 이송(advection)이 부재한다는 가정(예: 강체 평면 위를 미끄러지는 연성 블록)에 기초한다. 반대로, 구름 마찰은 Kalker가 개척한 바와 같이 접촉 패치를 통한 힘과 변형의 이송을 설명하기 위해 편미분 방정식(PDE)을 사용하여 정식화된다.

그러나 두 물체가 모두 변형 가능하거나 탄성 특성이 현저히 다른 경우, 특히 점탄성(타이어와 같은 고분자에서 흔히 나타남)이 존재하는 경우에는 이러한 이분법이 문제가 될 수 있다. 이러한 시나리오에서는 재료 입자와 내부 상태 변수가 접촉하는 물체의 류학적(rheological) 특성에 따라 서로 다른 속도로 이송되기 때문에, 미끄럼 마찰을 순수하게 집중된 형태(ODE 기반)로 모델링할 수 있다는 가정이 무너지게 된다. 기존 모델들은 이러한 역학을 통합하거나 점탄성 접촉 쌍의 운동학적 결합을 정확하게 처리하는 데 어려움을 겪는 경우가 많다.

방법론
본 논문은 선형 점탄성 류학 모델과 비선형 동적 마찰 모델을 결합한 통합된 이론적 프레임워크를 제안한다. 방법론은 다음과 같이 진행된다:

1. 류학 모델링: 접촉 계면은 "브러시 형태(bristle-like)" 요소들을 사용하여 모델링된다. 이 요소들의 점탄성 거동은 안정적인 선형 상태 공간 시스템으로 표현되는 일반화된 켈빈-보이트(Generalized Kelvin-Voigt, GKV) 모델을 사용하여 설명된다. 이를 통해 광범위한 주파수 대역에 걸친 복잡한 메모리 및 이력 현상(hysteresis)을 기술할 수 있다.
2. 점탄성-운동학적 유도: 저자들은 류학적 상태 공간 표현과 운동학적 관계를 결격시켜 지배 방정식 세트를 유도한다. 결정적으로, 이들은 재료 입자와 브러시 변형($z$)이 접촉 면적을 통해 이동하는 것을 고려하기 위해 표준 시간 미분 대신 물질(Eulerian) 미분(material derivatives)을 사용한다.
3. PDE 정식화: 구성 법칙을 운동학적 이송 방정식과 통합함으로써, 저자들은 접촉 영역 전체에 걸친 접선 방향 접촉력($f$)과 내부 상태 변수($\zeta$)의 진화를 지배하는 일련의 준선형 쌍곡형 편미분 방정식(PDE) 시스템을 유도한다.
4. 수학적 분석: 도출된 PDE 시스템은 쌍곡형 특성을 갖도록 분석된다. 저자들은 이 시스템이 리만 변수(Riemann variables)로 변환될 수 있음을 입증하며, 이를 통해 접촉 쌍의 탄성 분할에 의해 결정되는 특성 방향을 따라 방정식들을 디커플링(decoupling)할 수 있음을 보여준다.
5. 경계 및 초기 조건: 접촉 영역의 "선단(leading)", "중립(neutral)", "후단(trailing)" 엣지를 기반으로 적절한 경계 조건(BC)을 정식화하여 문제의 소산성(dissipativity)과 적정성(well-posedness)을 보장한다.

주요 기여
본 논문은 접촉 역학 및 마찰 모델링 분야에 세 가지 주요 기여를 한다:

• 분포형 프레임워크로의 체계적 확장: 본 연구는 동적 점탄성 마찰 모델을 집중형(ODE) 정식화에서 분포형(PDE) 프레임워크로 확장한다. 기판이 컴플라이언트(compliant)하거나 두 물체가 모두 변형 가능한 경우, 미끄럼 접촉 중에도 지배 방정식이 자연스럽게 PDE가 됨을 입증한다.
• 점탄성-운동학적 방정식의 유도: 저자들은 두 물체의 서로 다른 이송 속도를 명시적으로 고려하여, 미끄럼 및 구름 접촉 쌍 모두에 대한 마찰력, 재료 변형 및 내부 상태 변수의 결합된 진화를 지배하는 특정 클래스의 방정식을 유도한다.
• 미끄림 및 구름 역학의 통합: 이 이론은 미끄럼 마찰과 구름 마찰이 동일한 수학적 구조를 공유함을 밝혀낸다. 즉, 두 현상 모두 선형 또는 준선형 쌍곡형 PDE 시스템으로 설명된다. ODE 기반 미끄럼 모델과 PDE 기반 구름 모델의 차이는 운동의 유형에 따른 본질적인 차이가 아니라 류학(rheology)에 의존하는 것임을 보여준다.

결과
제안된 이론은 FrBD 모델의 확장과 GKV 요소를 결합하여 대표적인 미끄럼 및 구름 접촉 사례를 통해 입증된다:

• 미끄림 접촉: 점탄성 기판 위를 미끄러지는 블록의 경우, 모델은 기판이 강체인 극한(이송이 없는 경우)에서만 고전적인 ODE 거동을 회복한다. 기판이 점탄성인 경우, 모델은 공간적으로 변화하는 마찰력과 브러시 변형을 예측한다. 결과에 따르면 마찰력의 크기는 두 물체 사이의 강성비와 류학 모델의 소산 분지(dissipative branches) 수에 따라 달라진다. 고차 류학 모델은 공간적 완화(spatial relaxation)를 도입하여 접촉 길이를 따라 마찰력을 감쇄시킨다.
• 구름 접촉: 구름 접촉의 경우, 재료의 이송이 구름 운동에 내재되어 있으므로 강성비와 상관없이 PDE 특성을 온전히 유지한다. 분석 결과, 시스템의 쌍곡형 특성이 유지되며, 특성 속도는 주요 탄성 상호작용 모드에 의해 결정됨을 확인하였다.
• 과도 및 정상 상태 거동: 시뮬레이션은 점탄성 재료의 완화 시간(relaxation times)에 따라 시스템이 뚜렷한 과도 현상을 보임을 입증한다. 정상 상태 미끄럼에서 마찰력은 마찰 계수와 수직 하중의 곱에 의해 엄격히 제한되며, 강성비가 강체 쪽으로 유리할 때 이 경계값에 접근한다.

의의 및 주장
본 논문은 자신의 프레임워크가 광범위한 클래스의 점탄성 마찰 시스템에 적용 가능한 일반적인 이론적 환경을 제공한다고 주장한다. 주요 의의는 미끄럼 마찰과 구름 마찰 모델 사이의 개념적 모호성을 해결하는 데 있다. 두 현상 모두 동일한 클래스의 쌍곡형 PDE로 설명될 수 있음을 보여줌으로써, 본 연구는 미끄럼과 구름 문제를 ODE와 PDE 영역으로 나누는 전통적인 구분이 근본적인 물리적 차이가 아니라 단순화된 가정(특히 강체 기판 또는 디커플링된 운동학 가정)에 의한 산물임을 시사한다.

저자들은 자신의 접근 방식이 탄성적으로 상이한 쌍(elastically dissimilar pairs)에서 발생하는 재료 입자의 운동학을 정확하게 처리하며, 이는 기존의 브러시 유형 모델들이 슬립(slip)과 스핀(spin) 정의와 관련하여 겪었던 개념적 어려움을 해결한다고 강조한다. 본 이론은 이전의 탄성 접촉 쌍 연구에 대한 엄밀한 확장이며, 이제는 타이어와 같은 실제 응용 분야에서 발견되는 점탄성 재료를 모델링하기 위한 필요한 복잡성을 포함하고 있다. 결론적으로, 본 논문은 분포형 정식화가 단순히 대안적인 방법이 아니라, 점탄성 접촉의 정확한 모델링을 위해 필수적인 묘사이며, 미끄럼과 구름 역학 사이의 간극을 메우는 통합적 관점을 제공한다고 밝힌다.