테르자리 압밀 방정식 정확 해법

초록

본 논문은 테르자리 1차원 압밀 방정식에 대한 기존 테일러 근사 해법을 대체할 수 있는 정확 해를 제시한다. d’Alembert, Fourier, Laplace 방정식의 구조적 유사성을 이용해 파동‑감쇠 형태의 해 u(z,t)=u₀ e^{‑kz} cos(ωt‑kz) 를 도출하고, k와 ω를 기본 압밀 계수 c_v 와 시간 t 에 연결한다. 실험적 적용 사례를 통해 이 해가 테일러 해와 동일한 소산 시간(t₉₅, t₁₀₀)을 예측함을 확인한다.

상세 분석

논문은 먼저 d’Alembert 방정식이 2차 공간 미분과 2차 시간 미분을 포함한 파동 방정식임을 강조하고, Fourier 방정식이 열전달을 기술하지만 시간 차수가 하나 낮아 정적 상황에서는 Laplace 방정식으로 환원된다는 점을 지적한다. 이러한 수학적 구조의 차이를 이용해 테르자리 압밀 방정식(∂²u/∂z² = (1/c_v)∂u/∂t)을 d’Alembert 형태와 결합하면, 시간에 따라 감쇠되는 파동 해가 가능함을 보인다. 구체적으로 u(z,t)=u₀ e^{‑kz} cos(ωt‑kz) 형태를 제안하고, 편미분을 수행해 방정식에 대입하면 k² = 2c_v/ t, ω = k·c_v 라는 관계가 도출된다. 여기서 k는 깊이‑시간 의존적인 ‘압밀 변수’이며, 실험적으로 측정되는 c_v만으로 완전히 정의된다. 논문은 이 해가 기존 테일러 해의 차원 없는 파라미터 Z와 Tv를 필요로 하지 않으며, 오히려 압밀 과정의 물리적 메커니즘—즉, 과잉 간극압이 깊이에 따라 지수적으로 감소하고, 시간에 따라 진동‑감쇠 형태로 소멸한다—를 직접적으로 설명한다고 주장한다. 실용적 적용에서는 20 m 두께 점토층에 정하중 100 kPa를 가하고 c_v=0.02 m²/day인 경우, 해의 파라미터 ω와 k를 조정해 2 ~ 20 년 구간의 ‘소산 곡선’을 생성한다. 이 곡선들을 통해 t₁₀₀(완전 압밀 시점)을 17.5 ~ 18 년으로 추정했으며, 이는 테일러 해가 예측한 t₉₅≈15.5 년과 일치한다. 또한, 역으로 실험 데이터(Δu‑t 곡선)에서 t₁₀₀을 이용해 c_v를 추정할 수 있음을 보여준다. 논문은 이러한 접근이 1차원 이론의 한계를 보완하고, 3차원 압밀 해석으로 확장될 가능성을 제시한다. 다만, 해의 도출 과정에서 가정된 초기 하중 적용이 순간적이며, 비선형 토양 거동이나 다중 배수 조건을 고려하지 않은 점은 향후 연구 과제로 남는다.