마찰 전이에서 소산이 파열 속도에 미치는 영향

초록

본 논문은 프렌켈‑콘타브라 모델을 기반으로 마찰 전이 과정에 소산(dissipation) 효과를 도입한 새로운 이론을 제시한다. 정적 마찰에서 동적 마찰로 전이될 때 발생하는 파열 전파 속도는 전단‑정규응력 비만이 결정하고, 마찰계수와는 무관함을 보여준다. 반면 슬립 속도는 마찰에 민감하게 변한다. 이 모델은 지진부터 느린 슬립 이벤트까지 7 옥텟에 걸친 속도 범위를 포괄적으로 설명한다.

상세 분석

이 연구는 전통적인 프렌켈‑콘타브라(Frenkel‑Kontorova, FK) 모델을 거시 마찰 현상에 적용한 뒤, 소산(dissipation) 항을 추가함으로써 정적 마찰에서 동적 마찰로 전이되는 과정을 정량적으로 기술한다. FK 모델은 원자 혹은 입자들이 주기적인 포텐셜에 의해 배열되는 1차원 체인을 가정하고, 각 입자는 인접 입자와 탄성 결합을 가진다. 이를 거시 규모로 확대하면, 접촉면의 미세 구조가 유사한 비선형 탄성‑플라스틱 거동을 보인다고 볼 수 있다. 기존 모델에서는 마찰력(또는 마찰계수)만이 파열 전파 속도에 직접적인 영향을 미친다고 가정했으나, 실제 실험에서는 파열 전파 속도가 전단‑정규응력 비(τ/σ)와 강체 탄성 파라미터에 의해 주도되는 경우가 많았다.

논문은 소산 항을 도입해 운동 방정식에 댐핑 계수를 포함시켰다. 이때 중요한 수학적 결과는 파열 전파 속도 v_r이 댐핑(마찰) 항에 대해 독립적이라는 점이다. 구체적으로, 비선형 파동 방정식의 솔루션을 카우치(“kink”) 형태의 전이 파동으로 가정하고, 에너지 흐름을 분석하면 전파 속도는 다음과 같이 표현된다:

v_r = c · √(1 – (τ/σ)²)

여기서 c는 재료의 전단파 속도이며, τ/σ는 전단응력과 정규응력의 비율이다. 이 식은 마찰계수 μ가 어떠한 값이든 v_r에 영향을 주지 않음을 명시한다. 반면 슬립 속도 v_s는 에너지 손실을 나타내는 댐핑 계수 η에 비례하여 감소한다. 즉,

v_s ∝ 1/η · f(τ, σ, c)

이러한 결과는 직관에 반하는데, 마찰이 파열 전파 자체를 억제한다는 일반적인 기대와는 달리, 파열 전파는 응력 비와 탄성 파라미터에 의해 전적으로 결정된다는 점이다. 이는 지진학에서 관측되는 ‘느린 파열’ 현상과도 일맥상통한다. 실제 지진 파열 전파 속도는 수 km/s에서 수 km/일에 이르기까지 광범위한 변동을 보이며, 이 논문의 모델은 그 전체 범위를 하나의 통합 이론으로 포괄한다.

또한, 논문은 수치 시뮬레이션을 통해 파열 전파와 슬립 속도의 시간·공간 프로파일을 제시한다. 소산 파라미터를 변화시켜도 파열 전파 전선은 동일한 속도로 진행하지만, 슬립 영역의 폭과 전단 변형률은 크게 달라진다. 이는 실제 실험에서 관측되는 ‘파열 전선’과 ‘슬립 구역’이 서로 다른 물리적 메커니즘에 의해 제어된다는 점을 뒷받침한다.

결론적으로, 이 연구는 마찰 전이 현상을 두 개의 독립적인 동역학 변수—파열 전파 속도와 슬립 속도—로 분리함으로써, 기존의 단일 마찰계수 기반 모델의 한계를 극복한다. 특히, 전단‑정규응력 비가 파열 전파 속도를 결정한다는 점은 지진 위험 평가와 인공 구조물의 설계에 중요한 함의를 제공한다.