A stochastic model for elastoplastic contact of rough surfaces incorporating scale-dependent hardness

기술 요약: 규모 의존적 경도를 포함한 거친 표면의 탄소성 접촉에 관한 확률론적 모델

문제 정의
자연적 및 인공적 표면은 고유한 거칠기를 가지고 있어 접촉 돌기(asperity)에서의 응력 집중을 유발하며, 이는 종종 소성 변형을 일으킨다. 이러한 탄소성 접촉을 모델링하는 것은 두 가지 경쟁적인 요인으로 인해 복잡해진다: 다중 스케일에 걸친 자기 유사적(self-affine) 표면 토포그래피와 소성 변형의 "크기 효과(size effect)"이다. 후자는 재료의 경도가 불변의 상수가 아니라, 기하학적 필수 전위(Geometrically Necessary Dislocations, GNDs)의 축적으로 인해 관찰 규모(또는 배율)가 감소함에 따라 증가함을 의미한다. Persson의 확산 방정식이나 다중 돌기 접근법에 기반한 기존 모델들은 종종 일정한 경도를 가정하거나, 높은 해상도를 유지하면서 규모 의존적 경도를 효율적으로 통합하는 데 어려움을 겪는다. 또한, 기존의 공식들은 (확산 기반 접근법의 삼각 행렬과 같이) 복잡한 선형 방정식 시스템을 해결해야 하는 경우가 많아 계산 집약적일 수 있다.

방법론
저자들은 규모 의존적 경도를 갖는 탄소성 접촉 문제를 해결하기 위해 복합 채프먼-콜모고로프(compounded Chapman-Kolmogorov, C-K) 방정식에 기반한 새로운 확률론적 프레임워크를 개발하였다. 핵심 가정 및 단계는 다음과 같다:

1. 마르코프 과정 가정: 배율($\zeta$)에 따른 무작위 접촉 압력의 변화를 마르코프 과정으로 취급한다.
2. 적분 방정식 공식화: 저자들은 다음의 진화를 기술하는 세 가지 적분 방정식 시스템을 유도한다:
   • 탄성 접촉 압력의 확률 밀도 함수(PDF), $P_0(p, \zeta)$.
   • 상대적 소성 접촉 면적, $A^*_{pl}(\zeta)$.
   • 상대적 비접촉 면적, $A^*_{nc}(\zeta)$.
3. 규모 의존적 경도: 모델은 경도 함수 $H(\zeta) = H_0 \zeta^n$을 포함한다(여기서 $n$은 경도 지수). 이를 통해 관찰 규모가 작아짐(배율 증가)에 따라 탄성 압력의 상한값이 증가할 수 있게 한다.
4. 역류(Backflow) 메커니즘: 이는 일정한 경도를 가진 모델에서 발견되는 "재진입 불가(no re-entry)" 가정을 완화하는 중요한 이론적 진전이다. 경도가 규모에 따라 증가할 때, 이전에 소성 접촉 상태($p = H(\zeta')$)였던 표면 지점들이 더 높은 배율에서 탄성 접촉 영역이나 심지어 비접촉 영역으로 다시 돌아갈 수 있다. 이 확률의 "역류"는 적분 방정식 내에서 수학적으로 포착되며, 이는 본 모델을 기존의 일정 경도 공식과 구별 짓는 특징이다.
5. 수치적 구현: 배율 축은 로그 스케일을 사용하여 이산화된다. 솔루션은 흡수 및 비흡수 경계 조건을 갖는 전이 확률 커널을 활용하여 배율 단계를 따라 순차적으로 진행된다.

주요 결과
본 연구는 다양한 조건 하에서의 모델 동작을 보여주는 수치 해법과 수렴성 테스트를 제시한다:

• 수렴성: 모델은 이산화 지점($n_\zeta$)이 증가함에 따라 우수한 수렴성을 보이며, 선형 탄성 접촉에 대한 Persson의 이론 및 일정한 경도에 대한 기존의 폐쇄형 해(closed-form solutions)와 일치하는 결과를 나타낸다.
• 접촉 상태의 진화:
  • 경도 지수 ($n$): $n$이 증가하면(강한 크기 효과) 접촉 거동이 완전 소성에서 선형 탄성 쪽으로 이동한다. $n$이 충분히 크면, 완전 명목 접촉이 달성될 때까지 접촉은 탄성 상태를 유지한다.
  • 허스트 지수 ($H$): 더 높은 $H$ 값(더 매끄러운 표면) 역시 표면 높이의 분산과 압력 재분배율을 감소시킴으로써 탄성 거동을 촉진한다.
  • 대칭성 파괴: 탄성 압력 PDF가 $H_0/2$를 중심으로 대칭인 일정한 경도 모델과 달리, 규모 의전적 모델은 이 대칭성을 깨뜨린다. 역류 메커니즘으로 인해 현재의 경도 한계 $H(\zeta)$ 근처에 확률이 축적된다.
• 면적-하중 관계: 모델은 선형 탄성에서 탄소성 거동, 그리고 최종적으로 완전 소성으로 이어지는 매끄러운 전이를 예측한다. 저하중에서의 면적-하중 관계 기울기($\kappa$)는 경도 지수 $n$이 증가함에 따라 감소하며, 완전 소성 한계에서 선형 탄성 한계로 전이된다.
• 토포그래피 항복 매개변수: 저자들은 탄성 압력의 재분배와 경도 증가 사이의 경쟁을 정량화하는 새로운 토포그래피 항복 매개변수 $w(\zeta)$를 제안한다. 이 새로운 공식은 기존의 매개변수와 달리 영률(Young's modulus), 차단 파수(cut-off wavenumbers), 그리고 PSD 상수를 명시적으로 포함한다.
• 접촉 상태 도표: 미스핏 매개변수($M_{le}$ 및 $M_{fp}$)를 사용하여 순수 탄성 및 순수 소성 한계로부터의 편차를 정량화함으로써, $(H, n)$ 공간에서의 접촉 상태 도표를 생성한다. 이 도표는 선형 탄성, 탄소성, 그리고 완전 소성 영역을 구분한다. 결과는 일반적인 거친 표면($H \approx 0.8$)의 경우, 작은 크기 효과($n > 0$)만으로도 일정한 경도를 가정하는 모델의 예측과 달리 지배적인 탄성 거동을 유발할 수 있음을 시사한다.

의의 및 주장
본 논문은 접촉 역학 분야에 다음과 같은 중요한 기여를 한다고 주장한다:

1. 새로운 프레임워크: 거친 표면 접촉 분석에 복합 채프먼-콜모고로프 방정식을 적용한 선구적인 연구로서, 확산 기반 또는 다중 돌기 합산 방식과는 구별되는 확률 과정 프레임워크를 제공한다.
2. 계산 효율성: 적분 방정식 접근법은 확산 기반 모델에서 요구되는 대규모 선형 방정식 시스템(삼각 행렬)을 풀 필요가 없어, 잠재적으로 더 효율적인 계산 과정을 제공한다.
3. 크기 효과에 대한 물리적 통찰: 모델은 규모 의존적 경도가 어떻게 접촉 역학을 변화시키는지, 특히 소성 상태에서 탄성 상태로의 확률 "역류"가 어떻게 더 큰 규모(낮은 배율)에서 탄성 거동을 유도하는지에 대한 메커니즘을 제공한다.
4. 실용적 유용성: 유도된 접촉 상태 도표는 유한한 기계적 및 기하학적 특성을 바탕으로 접촉 영역(탄성 vs 소성)을 신속하게 식별할 수 있는 도구를 제공하여, 트라이볼로지 시스템의 설계 및 분석을 돕는다.
5. 광범위한 적용 가능성: 저자들은 인터페이스 특성의 규모에 따른 진화를 특징짓는 적분 방정식이 지진 역학, 전기 접촉, 접촉 전기화(contact electrification)와 같이 다중 스케일 거칠함이 포함된 다른 다학제적 분야에도 귀중한 통찰을 제공할 수 있다고 제안한다.

본 연구는 모델이 이전 연구들(Archard, Gao-Bower, Jackson-Streator)과 "돌기 간의 접촉(protuberance-on-protuberance)"이라는 기초적인 개념을 공유하고 있지만, 특정 기하학적 표현(규모 의존적 분산을 갖는 가우시안 거친 표면)과 C-K 방정를 통한 규모 의존적 경도의 포함이 독특한 지배 방정식과 물리적 예측을 이끌어낸다는 점을 강조하며 결론을 맺는다.