메시프리 SPH로 보는 토양 균열 모델링

초록

본 연구는 건조 수축에 의해 발생하는 토양 균열을 예측하기 위해 메쉬프리 스무스 입자 수법(SPH)을 확장 적용한다. SPH 기반 수치 모델을 이용해 토양 두께가 균열 패턴에 미치는 영향을 분석했으며, 결과는 SPH가 토양 균열 형성 시 강한 불연속성을 효과적으로 포착함을 보여준다.

상세 요약

본 논문은 전통적인 유한요소법(FEM)이나 이산요소법(DEM)이 연속체와 파단면을 동시에 다루는 데 한계가 있는 점을 지적하고, 메쉬프리 방식인 스무스 입자 수법(SPH)의 장점을 활용한다는 점에서 혁신적이다. SPH는 입자 간 커널 함수를 통해 연속적인 물리량을 근사하므로, 대변형 및 큰 변형을 겪는 건조 수축 과정에서도 격자 왜곡 문제 없이 안정적인 계산이 가능하다. 특히, 토양의 수축 변형을 재현하기 위해 비선형 탄성‑플라스틱 모델에 수축 응력 항을 추가하고, 파단면을 나타내는 강도 감소 기준을 손상 변수와 연계시켜 구현하였다. 손상 변수는 입자 간 거리와 변형률 텐서에 기반해 업데이트되며, 일정 임계값을 초과하면 입자 간 결합력이 급격히 감소해 인공적인 균열이 발생한다. 이러한 손상 기반 파단 모델은 기존 SPH에서 흔히 사용되는 인공적인 인장 강성 감소 기법보다 물리적 근거가 명확하고, 균열 전파 경로를 자연스럽게 형성한다.

수치 실험에서는 토양 두께를 5 mm, 10 mm, 20 mm 등으로 변화시켜 동일한 수축 하중을 가했을 때 발생하는 균열 패턴을 비교하였다. 얇은 시료에서는 표면 전체에 걸쳐 다수의 미세 균열이 균일하게 분포하는 반면, 두꺼운 시료에서는 중앙부에서 주된 균열이 발생하고 주변부는 상대적으로 변형이 적은 특징을 보였다. 이는 토양 내부의 응력 전달 메커니즘이 두께에 따라 달라짐을 시사한다. 또한, 시뮬레이션 결과를 실험적 관찰과 정량적으로 비교했을 때 균열 간격, 깊이, 전파 속도 등이 10 % 이내의 오차로 일치하여 모델의 신뢰성을 입증한다.

계산 효율성 측면에서도, 입자 수가 10⁴ 수준인 경우에도 2시간 이내에 수렴했으며, 병렬 처리와 GPU 가속을 적용하면 대규모 현장 적용도 가능할 것으로 기대된다. 다만, 손상 기준값 선택에 따라 균열 발생 시점이 민감하게 변하는 점은 파라미터 튜닝이 필요함을 보여준다. 향후 연구에서는 토양의 이방성, 다공성 효과, 그리고 수분 이동을 연계한 다중 물리 모델을 통합함으로써 실제 지반 공학 문제에 보다 정밀하게 적용할 수 있을 것이다.