입자 회전 제한이 토양 임계 상태에 미치는 영향

초록

본 연구는 입자 수준에서 회전을 방해하는 마찰계수를 도입해 단순 전단 시험을 수행한 분자동역학 시뮬레이션을 통해 토양의 임계 상태에 미치는 영향을 조사하였다. 회전 저항이 클수록 전단 강도가 증가하고, 구조적 이방성 및 힘 사슬 분포가 변함을 확인하였다. 결과는 회전 변수를 토양의 본 constitutive 모델에 포함시켜야 함을 시사한다.

상세 분석

본 논문은 입자 회전 자유도가 토양 입계 거동에 미치는 역할을 정량화하고자, 입자 간 접촉에 회전 마찰(토크) 한계를 부여하는 새로운 마찰계수 μ_r 을 도입하였다. 기존의 전단 시뮬레이션에서는 전단 강도와 체적 변형률이 입자 간 전단 마찰계수 μ_t 에만 의존한다고 가정했지만, 실제 입자들은 회전 저항을 갖는 경우가 많아 전단 전이 현상에 영향을 줄 수 있다. 연구팀은 다중 크기 분포를 가진 구형 입자 집합을 3‑차원 박스 안에 채우고, 주기적 경계조건 하에 단순 전단(단위 전단 변형률 γ̇ = 10⁻⁴ s⁻¹)을 적용하였다. 접촉 모델은 Hertz‑Mindlin 법칙에 기반해 정상·전단력과 토크를 동시에 계산했으며, 회전 마찰계수 μ_r 가 클수록 입자 회전이 제한되어 접촉면에서 발생하는 토크가 증가한다. 시뮬레이션 파라미터는 입자 직경 d = 0.5 mm, 평균 밀도 ρ = 2600 kg/m³, 초기 부피비 e₀ ≈ 0.75 이며, μ_r 값을 0, 0.2, 0.5, 0.8 네 단계로 변화시켰다.

결과 분석에서는 전단 응력 σ_xy 대 전단 변형 γ 곡선을 통해 임계 상태에 도달하는 시점을 확인하였다. μ_r 가 0인 경우 전형적인 임계 상태 전이(전단 강도 q_c ≈ 30 kPa, 부피비 e_c ≈ 0.68)를 보였지만, μ_r 가 0.8에 이르면 전단 강도가 약 45 kPa까지 상승하고 부피비는 e_c ≈ 0.62로 더 촘촘히 압축되었다. 이는 회전 저항이 입자 간 전단 전이를 억제하고, 보다 강인한 힘 사슬을 형성하게 함을 의미한다.

구조적 이방성은 Fabric Tensor F_ij 를 통해 정량화했으며, μ_r 증가에 따라 F_12 성분이 크게 상승해 전단 방향으로의 입자 배열이 강화되는 것을 확인했다. 또한, 접촉력 분포는 정규화된 정상력 P_n 과 전단력 P_t 의 확률밀도함수(PDF)에서 μ_r 증가 시 꼬리 부분이 두드러지게 늘어나, 강한 접촉이 더 많이 발생함을 보여준다. 이는 전단 저항이 높은 접촉이 전체 네트워크를 지배하게 되어 전단 강도가 상승하는 메커니즘과 일치한다.

마지막으로, 입자 회전 제한이 미세구조에 미치는 영향을 토양 공학의 본질적 모델에 반영하기 위해, 기존의 elasto‑plastic 모델에 회전 전이 변수 θ (예: 평균 회전각 제한)와 그에 대응하는 경화 법칙을 도입하는 방안을 제시한다. 이러한 확장은 전단‑압축 시험, 지반‑구조 상호작용, 그리고 비균질 매질의 동적 거동을 보다 정확히 예측할 수 있게 할 것으로 기대된다.