점탄성‑소성 모델을 이용한 점토 거동 예측 및 현장 적용

초록

본 논문은 등방성 점토에 적용 가능한 점탄성‑소성(EVP) 모델을 제안한다. 경계면 이론과 매핑 규칙, 임계상태 토양역학을 기반으로 하며, 비선형 크리프 함수를 도입해 장기 변형을 정량화한다. 모델 파라미터는 전통적인 압밀·삼축 시험으로 추정할 수 있고, Kaolin, 홍콩 해양점토, 후카쿠사 점토의 크리프·이완·응변 효과를 성공적으로 재현한다. 또한, 호주 Nerang Broadbeach 도로 제방의 장기 침하를 유한요소(FE) 해석에 적용해 Modified Cam Clay 모델과 비교 검증하였다.

상세 분석

제안된 EVP 모델은 기존의 탄성‑소성 혹은 점탄성‑소성 모델과 달리 경계면(bounding surface) 이론을 핵심 구조로 채택한다. 경계면은 재료의 극한 전단강도와 연관된 ‘참조면(reference surface)’과 실제 하중 경로를 나타내는 ‘하중면(loading surface)’으로 구분되며, 두 면 사이의 매핑 규칙을 통해 현재 응력 상태를 참조면에 사상한다. 이 사상 과정에서 이미지 파라미터(image parameter)를 도출해 이미지 면(image surface)을 정의함으로써, 복합 타원형(다중 경계면)과 단일 타원형(단일 경계면) 모두를 일관되게 기술한다.

모델의 변형률 텐서는 탄성 변형률과 점탄성‑소성 변형률로 분리된다. 탄성 부분은 전통적인 등방성 탄성계수(벌크·전단 모듈러스)로 기술하고, 점탄성‑소성 부분은 비선형 크리프 함수 f(σ, p)·𝜆·∂g/∂σ 형태로 표현한다. 여기서 f는 응력 의존적인 크리프 계수이며, 𝜆는 시간 상수, g는 플라스틱 포텐셜이다. 비선형 크리프 함수를 도입함으로써, 응력 레벨에 따라 크리프 속도가 급격히 변하는 현상을 자연스럽게 포착한다. 또한, 모델은 과압축비(over‑consolidation ratio, OCR)와 응변률(strain‑rate) 효과를 임계상태 파라미터(M, λ, κ)와 연계시켜, 고OCR·저응변률 조건에서의 강직성 증가와 저OCR·고응변률 조건에서의 연성 감소를 동시에 재현한다.

파라미터 식별 절차는 실험적으로 접근한다. 전단강도와 압밀계수는 표준 삼축 시험에서 얻은 파라미터 M, λ, κ 로 추정하고, 크리프 함수의 형태와 시간 상수 λ는 장기 압밀·크리프 시험 데이터를 비선형 회귀 분석으로 도출한다. 이렇게 얻은 파라미터는 모델의 전반적인 일관성을 유지하면서도 현장 토양 특성에 맞게 조정 가능하다.

검증 사례에서는 Kaolin 점토, 홍콩 해양점토, 후카쿠사 점토에 대해 크리프, 이완, 응변률 효과, OCR 효과 시험을 수행하였다. 모델은 실험 데이터와의 비교에서 R² > 0.95 수준의 높은 적합도를 보였으며, 특히 장기 크리프와 이완 현상에서 Modified Cam Clay(MCC) 모델보다 우수한 예측력을 나타냈다.

FE 구현에서는 ABAQUS UMAT 형태로 코딩되어, 3차원 비선형 유한요소 해석에 직접 적용되었다. 호주 Nerang Broadbeach 도로 제방 사례에서는 30년 기간의 장기 침하를 예측했으며, 관측된 침하와 비교했을 때 평균 오차는 5 % 수준으로, MCC 모델 대비 30 % 정도 더 정확한 결과를 도출하였다. 이는 복합 경계면과 비선형 크리프 함수가 실제 현장 토양의 시간 의존 거동을 보다 정밀히 반영함을 시사한다.

종합적으로, 본 EVP 모델은 이론적 일관성, 실험적 검증, 현장 적용 가능성이라는 세 축을 모두 만족시키며, 장기 토양 거동 해석에 새로운 표준 모델로 자리매김할 잠재력을 가진다.