스키 유발 눈판 사태의 3차원 탄성‑소성 시뮬레이션

초록

본 연구는 스키가 가하는 하중이 약층 위 눈판에 미치는 응력 변화를 3차원 유한요소 모델로 분석한다. 약층은 Mohr‑Coulomb‑Cap 모델을 이용해 탄성‑소성 거동을 구현했으며, 연속적인 연질‑연질‑연질 구조와 경‑연‑경 구조 두 가지 층형을 비교하였다. 결과는 경층이 ‘다리’ 역할을 하여 하중을 넓은 영역에 분산시키고, 약층의 전단·압축 응력이 감소함을 보여준다. 또한 전단 변형과 임계 균열 길이 사이의 연관성을 제시한다.

상세 분석

이 논문은 스키가 눈판에 가하는 동적 하중을 정적 하중으로 가정하고, 3차원 유한요소(FE) 해석을 통해 약층(weak‑layer) 위의 응력장 변화를 정량화한다. 약층의 소성 거동을 묘사하기 위해 Mohr‑Coulomb‑Cap(MCC) 모델을 채택했으며, 수치적 불안정을 방지하기 위해 Drucker‑Prager 근사화를 적용하였다. MCC 모델은 전단 파괴와 압축 캡을 동시에 고려함으로써 눈의 비선형 거동을 보다 현실적으로 재현한다.

시뮬레이션은 두 종류의 층구조를 사용하였다. 첫 번째는 동일한 연질(soft) 층이 3개 연속된 Soft‑Soft‑Soft 구조로, 전단강도와 탄성계수가 균일하게 분포한다. 두 번째는 상부와 하부에 경질(hard) 층이 삽입된 Hard‑Soft‑Hard 구조로, 경층이 ‘다리(bridge)’ 역할을 하여 하중을 넓은 범위에 전달한다. 각각의 구조에 대해 스키가 경사 방향으로 내려가는 경우와 경사와 평행하게 움직이는 경우를 모델링하였다.

결과적으로 Soft‑Soft‑Soft 구조에서는 스키 하중이 약층 바로 아래에서 -2330 Pa까지의 압축응력을 발생시키고, 약층 상부까지 -414 Pa 정도로 전달된다. 반면 Hard‑Soft‑Hard 구조에서는 경층이 하중을 흡수·분산시켜 약층 상부의 압축응력이 -108 Pa로 크게 감소한다. 이는 경층이 전단응력을 감소시키고, 약층의 파괴 가능성을 낮추는 ‘다리 효과’를 입증한다.

전단응력(τ_xz, τ_yz, τ_xy) 분석에서도 유사한 차이가 관찰되었다. Soft‑Soft‑Soft에서는 스키 양쪽 끝에서 집중된 전단응력이 나타나며, Hard‑Soft‑Hard에서는 전단응력이 더 넓게 퍼지고 최대값이 절반 이하로 감소한다. 이러한 전단응력 분포는 균열 전파 경로와 연관이 있다.

또한 논문은 전단응력과 체적 소성 변형률을 이용해 임계 균열 길이 a_c를 계산하였다. MCC 모델에 의해 얻어진 체적 소성 변형 영역이 넓을수록 a_c가 감소하여 균열이 더 쉽게 전파될 수 있음을 시사한다. 이는 기존의 스키 안정성 지수에 약층의 소성 거동을 포함시켜 보다 정밀한 위험 평가가 가능함을 의미한다.

수치해석 측면에서, 901 083 자유도(DoF)를 가진 모델을 16코어 클러스터에서 약 10분 내에 해결했으며, 메쉬는 하중 집중 영역에 고밀도로 배치하였다. 경계조건은 하부 고정, 나머지는 자유경계로 설정했으며, 중력 가속도는 시뮬레이션 중 제외하고 사후 보정식으로 추가하였다. 이러한 설정은 실제 눈판의 복잡한 변형을 단순화하면서도 핵심 물리 현상을 포착하는 데 충분했다.

전반적으로 이 연구는 눈판‑약층 시스템에 대한 3차원 탄성‑소성 해석을 최초로 수행했으며, 층 구조에 따른 하중 전달 메커니즘과 약층 파괴 위험을 정량적으로 제시한다. 결과는 스키어가 선택하는 주행 방향과 눈판의 층구조가 avalanche 위험에 미치는 영향을 이해하고, 현장 위험 평가 모델에 소성 거동을 통합하는 데 기여한다.