이방성 고체 계면에서 동적 반면 슬라이딩 안정성 분석

본 논문은 평면 대칭성을 갖는 두 이방성 선형 탄성체 사이에서 발생하는 반면(anti‑plane) 미끄럼 현상의 동적 안정성을 체계적으로 분석한다. 연구는 다음과 같은 흐름으로 전개된다. 1. **문제 설정 및 배경** - 기존의 단순 스프링‑블록 모델과 레이트‑앤드‑스테이트 마찰법칙을 검토하고, 이들이 속도 약화(velocity‑weakening) 조건에서도 안정성을 보일 수 있음을 소개한다. - 동질·등방성 경우의 반면 슬라이딩에서 파수 k 형태의 교란에 대한 임계 파수 |k|₍cr₎와 중성 모드 위상 속도 c 를 기존 문헌(예: Rice & Ruina, 1983; Ranjith & Rice, 2001)과 비교한다. - 이질 재료(비동질)와 이방성 효과가 정상응력 변화를 유발해 불안정성을 촉진한다는 점을 강조한다. 2. **이방성 탄성체의 모델링** - x₁‑x₂ 평면을 대칭면으로 가정함으로써 변위는 오직 u₃만 존재하고, 정상응력 σ₂₃는 슬립에 의해 변하지 않는다. 이는 전통적인 이방성 문제에서 발생하는 복합 변위·응력 결합을 회피한다. - 각 물질은 C₄₄, C₄₅, C₅₅ 세 개의 독립 탄성 상수와 밀도 ρ 로 기술된다. 3. **동역학 방정식 도출** - 라플라스 변환(p)과 푸리에 변환(k)를 적용해 2차 편미분 방정식 C₅₅∂²u/∂x₁² + 2C₄₅∂²u/∂x₁∂x₂ + C₄₄∂²u/∂x₂² = ρp²u 를 얻는다. - 유한성을 만족하는 해를 찾기 위해 α 를 도입하고, α²C₄₄ + 2ikC₄₅α − (ρp² + C₅₅k²)=0 를 푼다. - 여기서 유효 전단계수 μ ≡ C₄₄C₅₅ − C₄₅² 와 전단파 속도 c₁ ≡ μ/√(C₄₄ρ) 를 정의한다. 4. **계면 조건 및 전단응력‑슬립 연계** - 물질 #1과 #2 각각에 대해 전단응력 교란 T̂ 와 변위 교란 Û 사이의 관계를 T̂ = −μ|k|(1 + p²/(k²c₁²))Û, T̂′ = +μ′|k|(1 + p²/(k²c₁′²))Û′ 로 정리한다. - 계면에서 전단응력 연속성과 슬립 정의 D(k,t)=U⁺−U⁻ 를 이용해, 슬립 속도 교란 V−V₀ = ∂D/∂t 와 전단응력 교란 사이에 하나의 동역학 연산자를 얻는다. 5. **레이트‑앤드‑스테이트 마찰법칙 적용** - τ = τ₀ + aσ₀ ln(V/V₀) + bσ₀ ln(V₀θ/L) 와 진화식 dθ/dt = 1 − Vθ/L 를 선형화한다. - 선형화 결과는 마찰 파라미터 a, b, L, σ₀ 가 포함된 식 p·