붕괴하는 모래의 미세역학 모델과 살아있는 퀵샌드 현상

초록

**
본 논문은 입자 규모에서 토양·모래의 거동을 재현하는 이산 입자법(Discrete Element Method)을 기반으로, 응집력을 갖는 미세 입자들의 접촉 동역학을 구현한다. ‘살아있는 퀵샌드’라 명명한 특수 토양의 붕괴·침투 현상을 시뮬레이션하고, 침투 물체의 밀도와 구동력에 따른 속도‑힘 관계가 (F\propto v^{2}) 형태의 1/2 차 전력법칙을 따른다는 것을 확인한다.

**

상세 분석

**
이 연구는 전통적인 접촉 동역학(Contact Dynamics) 프레임워크에 두 가지 핵심 요소를 추가함으로써 기존 DEM 모델의 한계를 극복한다. 첫째, 미세 입자 사이의 응집력을 파라미터화하여, 물리적으로는 정전기·수분에 의한 미세입자 결합을 모사한다. 이 응집력은 접촉 시 정상력에 일정한 ‘접착력’ (F_c) 를 부여하고, 파단 조건을 설정함으로써 입자 간 결합·파괴를 동적으로 제어한다. 둘째, ‘살아있는 퀵샌드’라 부르는 특수 토양은 초기에는 매우 낮은 체적밀도(≈0.4)와 높은 공극률을 가지지만, 외부 충격이나 하중이 가해지면 급격히 체적밀도가 증가하고 전단강도가 급격히 상승한다는 실험적 관찰을 모델에 반영한다. 이를 위해 입자 간 마찰계수와 전단강도 함수를 체적밀도에 의존하도록 설계했으며, 파괴 후 재구성 과정에서 입자 재배열을 허용한다.

시뮬레이션 설정은 원통형 용기 안에 3차원 구형 입자들을 무작위 배치하고, 상부에 일정한 질량을 가진 침투체(‘intruder’)를 놓는 방식이다. 침투체에 가해지는 구동력 (F_{drive}) 를 단계적으로 증가시키며, 침투 속도 (v) 를 측정한다. 결과는 두 가지 주요 패턴을 보인다. (1) 체적밀도가 낮을수록 침투 저항이 작아 침투 깊이가 크게 증가한다. (2) 일정 범위 이상의 구동력에서 (F_{drive}) 와 (v) 사이에 (F_{drive}\propto v^{2}) 관계가 성립한다. 이는 전통적인 유체역학에서 관찰되는 드래그 힘 (F_D = \frac{1}{2} C_d \rho A v^2) 와 유사하지만, 여기서는 입자 간 응집·파괴와 체적밀도 변화가 복합적으로 작용한다는 점이 차별점이다. 또한, 시뮬레이션에서 관찰된 ‘전이점’(critical force) 은 입자 네트워크가 전단강도 한계에 도달해 붕괴하기 시작하는 순간과 일치한다. 이 전이점은 입자 크기 분포, 응집력 강도, 마찰계수 등에 민감하게 반응하므로, 실제 현장 적용 시 토양 특성에 대한 정밀한 파라미터 추정이 필요함을 시사한다.

모델 검증을 위해 실험적 퀵샌드 데이터와 비교했으며, 침투 깊이와 속도‑힘 곡선이 정량적으로 일치함을 확인했다. 특히 ‘살아있는 퀵샌드’라 명명된 토양은 물리적 습윤 상태와 미세 입자 응집력이 동시에 작용해, 일반적인 건조 퀵샌드와는 다른 비선형 저항 특성을 보인다. 이러한 결과는 토양 공학, 지반 안정성 평가, 그리고 로봇·탐사 장비 설계에 있어, 급격한 토양 붕괴 현상을 예측하고 회피 전략을 수립하는 데 유용한 이론적·수치적 기반을 제공한다.

**