🔬 물리학

Influence of Kelvin-Voigt Viscoelasticity on Transient Pore Pressure Response of Cylindrical Cavities in Unsaturated Soil

이 논문은 켈빈-포이그트(Kelvin-Voigt) 점탄성 골격을 통합함으로써 불포화 토양 내 원통형 공동에 대한 해석적 해를 확장하며, 변위와 응력은 거의 영향을 받지 않는 반면 간극수압은 특히 펄스 하중 조건에서 체적 점탄성이 빠른 압축파에 대한 저역 통과 필터로 작용함에 따라 크게 감쇄된다는 것을 밝혀냈다.

원저자: Xiong Hao

게시일 2026-10-13

📖 4 분 읽기☕ 가벼운 읽기

CC BY 4.0

원저자: Xiong Hao

원본 논문은 CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/) 라이선스로 제공됩니다. ✨ 이것은 아래 논문에 대한 AI 생성 설명입니다. 저자가 작성하거나 승인한 것이 아닙니다. 기술적 정확성을 위해서는 원본 논문을 참조하세요. 전체 면책 조항 읽기

우리 발밑의 지표면은 결코 단순하고 단단한 하나의 덩어리가 아닙니다. 많은 곳에서, 특히 지표면 근처의 토양은 고체 입자와 그 사이의 아주 작은 틈새에 갇힌 물, 그리고 공기가 뒤섞인 복잡한 혼합물입니다. 이러한 종류의 흙이 말뚝 박기, 지하 폭발, 또는 터널의 급격한 굴착과 같은 갑작스러운 사건에 의해 교란될 때, 그 힘은 예측하기 어려운 방식으로 토양을 통해 파동처럼 퍼져 나갑니다. 엔지니어들은 오랫동안 토양의 고체 골격이 완벽하게 탄성적이라고 가정하는 수학적 모델, 즉 딱딱한 스프링처럼 압축된 후 즉각적이고 예측 가능하게 원래 상태로 되돌아온다고 간주하는 모델에 의존해 왔습니다. 그러나 실제의 토양, 특히 부드러운 점토나 유기질 실트(silt)는 다르게 행동합니다. 이들은 '점탄성(viscoelastic)'을 띠는데, 이는 재료가 단순히 강성(stiffness) 때문에 움직임에 저항하는 것이 아니라, 내부 마찰로 인해 움직임에 저항하며 시간이 지남에 따라 에너지를 흡수하고 더 느리게 반응하는 특성을 의미합니다. 이러한 시간 의존적 거동이 토양 내부의 물과 공기 압력을 어떻게 변화시키는지 이해하는 것은, 갑작스러운 충격을 견뎌낼 수 있는 안전한 지하 구조물을 설계하는 데 매우 중요합니다.

타이저우 대학교(Taizhou University)의 한 연구자는 기존 모델에 이러한 점성 거동을 포함하도록 업데이트함으로써 이 현상에 대한 우리의 이해를 정교하게 개선했습니다. 이 연구는 내부 폭발이나 급격한 굴착과 같은 조건을 시뮬레이션하기 위해, 갑작스러운 반경 방향 하중(radial load)을 받는 지중 원통형 구멍 내부에서 어떤 일이 발생하는지에 초점을 맞춥니다. 연구자는 토양이 변형되는 속도를 고려하는 켈빈-보이그트(Kelvin-Voigt) 모델이라는 특정 수학적 프레임워크를 적용하여, 공극 내 유체의 압력이 시간에 따라 어떻게 변하는지 계산할 수 있었습니다. 이 작업은 기존의 탄성 토양 솔루션을 기반으로 하지만, 토양의 골격이 즉각적으로 반응하지 않는다는 결정적인 교정 사항을 도입했습니다. 대신, 토양 입자 내부의 마찰이 댐퍼(damper) 역할을 하여 힘의 전달을 늦추고, 공극 내의 유체와 가스를 통과하는 압력파를 변화시킵니다.

이 연구의 가장 놀라운 발견은 토양 자체의 움직임과 고체 입자에 가해지는 응력은 기존의 더 단순한 모델이 예측한 것과 거의 정확히 일치하는 반면, 토양 내부의 물과 공기 압력은 현저하게 떨어진다는 점입니다. 토양을 점탄성으로 모델링할 경우, 하중이 가해지는 속도에 따라 급격한 하중에 의해 생성되는 피크 압력이 탄성 예측치보다 최대 43%까지 감소할 수 있습니다. 만약 하중이 계단 함수(step function)처럼 점진적으로 가해진다면 감소 폭은 약 2% 정도로 작습니다. 그러나 하중이 날카로운 삼각형 펄스(triangular pulse) 형태로 가해진다면 그 감소는 극적입니다. 이는 토양의 점성적 성질이 높은 유체 압력을 생성하는 급격한 압착에 저항하기 때문입니다. 고체 골격이 유체로 전달되기 전에 내부 마찰을 통해 에너지의 일부를 흡수함으로써, 마치 압력 스파이크를 억제하는 필터처럼 작동하는 것입니다.

또한 이 연구는 모든 유형의 토양 강성이 이 효과에 동일하게 기여하는 것은 아니라는 점을 밝혀냈습니다. 압력 감소는 거의 전적으로 부피 변화에 대한 토양의 저항인 체적 점성(bulk viscosity)에 의해 주도됩니다. 형태 변화에 대한 저항인 전단 점성(shear viscosity)은 공극 내 압력이 토양이 비틀리거나 전단되는 것이 아니라 압착됨으로써 생성되기 때문에 무시할 수 있는 역할을 합니다. 이러한 구분은 엔지니어들에게 매우 중요한데, 이는 토양의 모든 가능한 변형 방식을 특성화할 필요 없이, 토양의 부피 관련 점성 특성을 측정하는 데에만 집중하면 정확하게 공극 압력을 예측할 수 있음을 의미하기 때문입니다. 더욱이, 이 효과는 토양이 단단하고 투수성이 높아 유체가 자유롭게 이동할 수 있는 경우, 그리고 토양이 물로 중등도 이상 포화된 상태일 때 가장 두드러집니다. 부드럽고 투수성이 낮은 점토의 경우 점성 교정 효과가 미미하며, 이는 기존의 탄성 모델이 해당 재료들에 대해서는 여전히 상당히 정확하다는 것을 시사합니다.

연구자는 이 발견을 검증하기 위해 고급 수치 해석 기법을 사용하여 토양의 거동을 지배하는 복잡한 방정식들을 풀었으며, 이를 통해 결과의 정밀도를 매우 높은 수준으로 확보했습니다. 분석 결과, 토양을 통해 에너지를 전달하는 주요 파동인 '빠른 압축파(fast compressional wave)'가 토양의 점성에 의해 가장 큰 영향을 받는 반면, 더 느린 확산파(diffusive waves)는 거의 변하지 않는다는 것이 확인되었습니다. 이러한 선택적 댐핑(damping) 현상은 왜 서로 다른 유형의 하중이 서로 다른 수준의 압력 감소를 유발하는지를 설명해 줍니다. 즉, 폭발과 같이 고주파 성분을 가진 하중이 느린 하중보다 더 효과적으로 필터링됩니다. 연구는 결론적으로 전통적인 탄성 모델이 지반의 움직임이나 토양 구조가 받는 응력을 예측하는 데에는 신뢰할 만하지만, 점탄성 지반에서의 위험한 유체 압력을 과대평가하는 경 경향이 있다고 밝힙니다. 사질 실트(sandy silts)와 같은 토양에서 발생하는 갑작스럽고 고에너지인 사건이 포함된 프로젝트의 경우, 이 점성 댐핑을 무시하면 지나치게 보수적인 설계를 초래할 수 있는 반면, 이를 고려하면 작용하는 힘에 대한 더 현실적인 그림을 제공할 수 있습니다.

기술 요약: 불포화 토양 내 원통형 공동의 과도적 간극수압 응답에 미치는 Kelvin-Voigt 점탄성 영향

문제 제기
불포화 토양 내 원통형 공동의 과도적 동적 응답에 대한 기존의 해석적 해법들은 선형 탄성 고체 골격(skeleton) 가정을 기반으로 한다. 그러나 연약 점토, 유기 실트, 이탄(peat)과 같은 천연 다공성 매질은 탄성 모델로는 포착할 수 없는 속도 의존적 변형을 나타낸다. 기존 연구들이 불포해 포로탄성(poroelasticity)을 다루거나(Li 등 [27]), 더 단순한 구성에서의 점탄성 파동 전파를 다룬 바 있으나(Guo 등 [30]), 점탄성 골격을 가진 3상(고체-액체-기체) 불포화 매질 내 원통형 공동의 과도적 응답에 대한 종합적인 해석적 해법은 부재했다. 본 논문은 탄성 기초 해법을 Kelvin-Voigt(KV) 점탄성으로 확장함으로써 이 공백을 메우며, 특히 골격 점성이 과도적 하중 하에서 간극수압, 변위 및 응력에 어떻게 영향을 미치는지 조사한다.

방법론
본 연구는 다음 프레임워크를 기반으로 한 해석적 접근 방식을 채택한다:

지배 방정식: 문제는 Van Genuchten 토양-수분 특성 곡선을 포함하는 3상 다공성 매질의 연속체 이론(Wei와 Muraleetharan [15])을 사용하여 정식화된다. 시스템은 $u_s$ - $w_f$ - $w_a$ (고체, 액체, 기체 상대 변위) 정식화를 사용하여 무차원 항으로 기술된다.
점탄성 확장: Kelvin-Voigt 구성 모델이 토양 골격에 적용된다. 라플라스 영역에서의 탄성-점탄성 대응 원리(elastic-viscoelastic correspondence principle)를 사용하여, 탄성 계수는 주파수 의존적 복소 계수인 $\lambda^*(s) = \lambda^*(1 + t^*_r s)$ 및 $\mu^*(s) = 1 + t^*_r s$ 로 대체된다 (여기서 $t^*_r$ 은 무차원 이완 시간이다).
해의 구조: 대응 원리를 통해 기존 탄성 해의 구조를 유지할 수 있다. 헬름홀츠 분해(Helmholtz decomposition)를 통해 결합된 시스템을 비결합된 헬름홀츠형 방정식으로 축소한다. 일반해는 세 가지 압축 파동 모드(P1, P2, P3)에 대해 수정된 베셀 함수( $K_n$ )를 사용하여 표현된다.
수치적 구현: 경계값 문제(투수 및 불투수 조건)는 라플라스 영역에서 해결된다. 시간 영역의 해는 복소 평면상의 수직 경로와 가우스 구적법을 사용하는 Den Iseger의 수치적 역 라플라스 변환 알고리즘을 통해 얻어진다. 점탄성 계수가 복소수 값을 가지기 때문에 이 방법이 기존 탄성 연구에서 사용된 Stehfest 알고리즘보다 선택되었다.
검증: 해법은 탄성 한계( $t^*_r = 0$ ), 포화 한계( $S_r \to 1$ ), 그리고 독립적인 주파수 영역 분산 분석(Guo 등 [30])을 통해 검증된다.

주요 기여 및 결과
본 논문은 벤치마크 해석 해법을 제공하며 골격 점탄성의 구체적인 효과를 정량화한다:

간극수압 감쇠: 가장 중요한 발견은 골격 점탄성이 간극수압을 강력하게 감쇠시킨다는 것이다. 하중 유형과 이완 시간에 따라, 피크 간극수압은 $t^*_r = 2.0$ 일 때 단계 하중(step loading)의 경우 2%, 삼각형 펄스 하중(triangular pulse loading)의 경우 43%까지 감소한다.
변위 및 응력의 불변성: 간극수압과 대조적으로, 반경 방향 변위와 원주 방향 응력은 골격 점탄성에 의해 거의 변화되지 않는다. 상대적 변화는 $10^{-6}$ 미만이며, 이는 점탄성 지반에서도 총 응력 및 변위장을 예측하는 데 탄성 해법이 여전히 정확함을 시사한다.
체적 vs 전단 기여: 분해 분석 결과, 간극수압 감쇠는 거의 전적으로 체적(bulk) 점탄성( $\lambda$ )에 의해 주도된다. 전단(deviatoric) 점탄성( $\mu$ )의 기여는 무시할 만한데, 이는 간극수압 생성이 전단 변형이 아닌 체적 변형에 결합되어 있기 때문이다.
주파수 영역 메커니즘: 점탄성 효과는 저역 통과 필터(low-pass filter) 역할을 한다. 점탄성 대비 탄성 간극수압의 비를 나타내는 전달 함수 $H(s)$ 는 주파수가 증가함에 따라 단조 감소한다. 이는 왜 고주파 과도 하중(예: 삼각형 펄스 또는 폭발)이 저주파 지속 하중(예: 단계 함수)보다 더 큰 감쇠를 경험하는지를 설명한다.
파동 모드 선택성: 파동 모드 분석을 통해 빠른 P1 압축파(골격 지배)가 점탄성 댐핑의 주요 대상임을 확인하였으며, 이 파동의 감쇠 길이는 10~30배 감소한다. 느린 P2 및 P3 확산 모드는 골격 계수보다는 수송 계수에 의해 지배되므로 거의 영향을 받지 않는다.
매개변수 민절도: 점탄성 보정의 크기는 단단하고( $\lambda^*$ 높음), 투수성이 높으며( $b^*_f$ 낮음), 중등도 이상의 포화도( $S_r > 0.6$ )를 가진 토양에서 가장 크다. 연약하고 투수성이 낮은 점토에서는 보정이 미미하다(< 2%).

의의 및 주장
본 논문은 점탄성 지반에 대한 엄밀한 해석적 벤치마크를 제공하고, 점탄성 지반에서 탄성 모델이 가질 수 있는 간극수압 과대 예측을 교정한다는 점에서 그 주요 의의가 있다고 주장한다.

실무적 시사점: 점탄성 지반에서 간극수압 예측 시 탄성 해법은 비보수적(non-conservative)이며, 특히 과도적인 폭발 하중 하에서 그러하다. 과대 예측 정도는 상당할 수 있으나(최대 43%), 변위 및 총 응력 예측은 신뢰할 수 있다.
방법론적 효율성: 점탄성으로의 확장이 복잡한 경계값 문제를 재유도할 필요 없이, 라플라스 영역에서의 대수적 매개변수 치환을 통해 달성될 수 있음을 입증하였다.
실험적 가이드: 결과에 따르면 간극수압 응답을 추정하기 위해 체적 점탄성 매개변수의 실험적 특성 파악만으로 충분하며, 전단 점성은 간극수압 응답에 미치는 영향이 미미하다.

본 연구는 KV 모델이 속도 의존적 강성(stiffening)을 포착할 수 있지만, 향후 연구에서는 응력 이완을 포착하기 위해 표준 선형 고체(standard linear solid) 또는 분수 미분 모델로 확장할 수 있으며, 탄성 터널 라이닝의 강성 효과를 포함할 수 있다고 결론짓는다.

유사한 논문