이질 재료 계면에서의 동적 파열·단층 개방·근접 입자 운동

초록

본 연구는 전단 지배 하에 서로 다른 탄성 물성을 가진 두 고체 사이의 평면 계면을 2‑D 격자 입자 모델(LPM)로 모사하여, 파열 전파와 동시에 발생하는 단층 개방 현상 및 근접 입자들의 비대칭 운동을 조사한다. 연성(연질) 물질 쪽으로 전파되는 파열은 느린 레일리 파속에 근접한 속도로 전파되며, 국소적인 계면 분리가 동반된다. 정상 방향 입자 속도가 전단 방향보다 크게 나타나며, 연성 측에서의 운동이 경질 측보다 현저히 강하다. 이러한 결과는 실험적 폼 고무 관찰과 Weertman‑dislocation 이론, Schallamach 파동과 일치한다.

상세 분석

본 논문은 이질 재료 계면을 따라 전단‑주도 파열이 어떻게 전파되는지를 2‑D 격자 입자 모델(LPM)로 정량화한다. 모델은 삼각 격자 배열의 입자들로 구성되며, 입자 간 상호작용은 수정된 Lennard‑Jones 포텐셜을 사용해 탄성 상수와 절단 길이를 직접적으로 매핑한다. 두 블록은 동일한 밀도(ρ)지만 서로 다른 전단 탄성계수(k)와 라메 상수(μ)를 갖으며, 인터페이스 강성 K는 두 물질의 탄성 상수 차이에 대한 mismatch 함수로 정의된다(식 5). 이때 K는 0.5 < k₁/k₂ < 1 구간에서 최소가 되며, 실제 단층에서 관측되는 물성 대비(전단파속 차이 0.7–1.0)와 일치한다.

전단 하중은 외부 경계에 일정한 전단 응력 τ∞와 압축 응력 σ∞를 가함으로써 구현하고, 마찰은 Coulomb 법칙(τ = f σ)으로 기술한다. 초기 파열은 좌측 가장자리의 거칠기를 통해 인위적으로 촉발된다. 시뮬레이션 결과, 파열 전파 속도는 연성 물질(소프트 블록)의 레일리 파속 V_R에 근접하며, 파열 전파 중에 국소적인 인터페이스 분리(단층 개방)가 발생한다. 이 개방은 두 측면의 정상 변위 차이로 나타나며, 정상 응력 σ가 순간적으로 0이 되면서 전단 응력 τ도 급격히 감소한다.

입자 운동을 자세히 살펴보면, 정상 방향 속도와 가속도가 전단 방향보다 수배 크게 나타난다. 특히 연성 측에서는 정상 변위가 경질 측의 약 2배에 달하며, 이는 Weertman이 제시한 “material contrast‑induced normal stress change”와 일치한다. 물성 대비가 40 % 이상이면 정상 변위 파형이 Schallamach 파동과 유사한 형태를 보이며, 이는 고무 실험에서 보고된 detachment wave와 동일한 메커니즘으로 해석된다.

또한, 파열 전파가 진행될수록 입자 운동은 펄스 형태로 점점 뾰�게 되고, 진폭이 증가한다. 이는 Andrews와 Ben‑Zion(1997)이 제시한 “traveling wave radiation”에 의한 파열 강화 메커니즘과 부합한다. 파열 과정에서 방출되는 표면파는 합성 지진계측기(seismogram)에서 큰 정상 및 전단 성분을 보이며, 방출된 지진 에너지는 슬립·개방·치유(정지) 단계에서의 입자 운동에 의해 공급된다. 따라서 전단 응력 감소는 부분적인 응력 강하(partial stress drop)로 해석된다.

이러한 결과는 기존의 고전적 마찰‑전단 모델이 설명하기 어려운 “열 흐름 역설”을 완화시키는 메커니즘을 제공한다. 정상 진동에 의한 유효 정상 응력 감소는 마찰열 발생을 억제하고, P‑파 방출을 증가시켜 관측된 비정상적인 P‑파 에너지와도 일치한다.