유승경
(Dept. of Civil Engineering, Myongji College)
,
김주봉
(School of Civil and Environmental Engineering, Chung-Ang Univ.)
,
한중근
(School of Civil and Environmental Engineering, Urban Design and Study, Chung-Ang Univ.)
,
홍기권
(Institute of Technology Research and Development, Korea Engineering & Construction)
,
윤중만
(Dept. of Construction Information System, ShinAnsan Univ.)
,
이강일
(Dept. of Civil Engineering, Deajin Univ.)
지반 공동 발생으로 인한 지반함몰 사고를 미연에 방지하기 위해서는 발생된 공동 주변 지반의 이완영역에 대한 역학적 특성 평가와 더불어 이에 상응하는 적합한 보강 대책을 신속히 수립하는 것이 필요하다. 본 논문에서는 지중 공동 주변의 이완영역 발생에 대한 메커니즘 분석을 위해 유한요소 수치해석을 실시하였다. 수치해석에서는 지중 공동과 이완영역을 모사하기 위해 강제변위법을 적용하였으며, 수치해석 결과는 실내모형실험을 수행한 기존 연구 결과를 이용하여 검증하였다. 유한요소 수치해석 결과로부터 토사 유실로 인한 공동 주변 지반의 간극비 분포 특성을 파악하였으며, 전단응력 감소 특성을 분석하여 이완영역의 범위를 정량적으로 제시하였다.
지반 공동 발생으로 인한 지반함몰 사고를 미연에 방지하기 위해서는 발생된 공동 주변 지반의 이완영역에 대한 역학적 특성 평가와 더불어 이에 상응하는 적합한 보강 대책을 신속히 수립하는 것이 필요하다. 본 논문에서는 지중 공동 주변의 이완영역 발생에 대한 메커니즘 분석을 위해 유한요소 수치해석을 실시하였다. 수치해석에서는 지중 공동과 이완영역을 모사하기 위해 강제변위법을 적용하였으며, 수치해석 결과는 실내모형실험을 수행한 기존 연구 결과를 이용하여 검증하였다. 유한요소 수치해석 결과로부터 토사 유실로 인한 공동 주변 지반의 간극비 분포 특성을 파악하였으며, 전단응력 감소 특성을 분석하여 이완영역의 범위를 정량적으로 제시하였다.
In order to prevent the ground subsidence accidents caused by the occurrence of underground cavity, it is necessary to evaluate the mechanical characteristics in the relaxation zone of the underground cavity. Also, the relaxation zone including underground cavity be appropriately reinforced. This pa...
In order to prevent the ground subsidence accidents caused by the occurrence of underground cavity, it is necessary to evaluate the mechanical characteristics in the relaxation zone of the underground cavity. Also, the relaxation zone including underground cavity be appropriately reinforced. This paper described analysis results based on finite element method that was conducted to analyze the mechanism for occurrence of the relaxation zone around the underground cavity. The finite element analysis applied in forced displacement was carried out to simulate the underground cavity and relaxation zone, and then there were compared with previous research results. The analysis results showed that the void distribution of soil around the underground cavity has figured out. As a result, the area of the relaxation zone could be quantitatively presented by reduction characteristics of the shear stress.
In order to prevent the ground subsidence accidents caused by the occurrence of underground cavity, it is necessary to evaluate the mechanical characteristics in the relaxation zone of the underground cavity. Also, the relaxation zone including underground cavity be appropriately reinforced. This paper described analysis results based on finite element method that was conducted to analyze the mechanism for occurrence of the relaxation zone around the underground cavity. The finite element analysis applied in forced displacement was carried out to simulate the underground cavity and relaxation zone, and then there were compared with previous research results. The analysis results showed that the void distribution of soil around the underground cavity has figured out. As a result, the area of the relaxation zone could be quantitatively presented by reduction characteristics of the shear stress.
* AI 자동 식별 결과로 적합하지 않은 문장이 있을 수 있으니, 이용에 유의하시기 바랍니다.
문제 정의
본 논문에서는 하수관 상부 손상 시, 토사유출에 의한 지중 공동 및 이완영역의 발생 메커니즘을 분석하기 위하여 유한요소 수치해석을 실시하였다. 수치해석에서는 지중의 공동과 함께 이완영역을 모사하기 위해 강제변위법을 적용하였다.
본 연구에서는 지중 공동 주변 지반의 역학거동을 분석하기 위해 유한요소해석을 실시하였다. 이 때, 본 수치해석 결과의 신뢰성 검증을 위해 실내모형실험과 개별요소 수치해석을 수행한 선행연구(Kim et al.
본 연구에서는 하수관 상부 손상 시, 토사유출에 의한 지중 공동 및 이완영역의 발생 메커니즘을 분석하기 위하여, 강제변위법을 활용한 유한요소해석을 수행하여 공동 및 이완영역을 모사하였으며, 기존 연구결과와의 비교를 통해 이를 검증하였다. 또한 공동 주변 지반의 전단응력감소비 분포 특성을 분석하여 공동 및 이완영역의 범위를 정량적으로 제시하였으며, 연구결과를 요약하면 다음과 같다.
상기에서는 탄성 모델을 이용하여 지반의 공동과 이완 영역을 모사함에 있어 유한요소법의 적용성을 평가하였다. 그러나 보다 엄밀히 실제 지반의 거동을 분석하기 위해서는 탄소성 이론에 근거하여 해석함이 바람직하며, 더욱이 지중에서의 공동 발생은 지하수위 조건에 영향을 받기 때문에 이러한 제반 조건을 고려한 해석이 필요하다.
제안 방법
강제변위를 적용한 유한요소해석 결과가 공동 및 이완영역 모사를 적절하게 나타낼 수 있는지 평가하기 위하여, 선행 연구(Kim et al., 2017) 및 본 해석결과에 의한 간극비 분포도를 비교하였다. Fig.
또한 트랩도어를 통한 토사 유출 현상을 안정적으로 모사하기 위해 모형지반 하단부에 강재 요소를 추가 적용하였으며, 강재 중앙부의 Geometry를 제거하여 트랩도어 개방을 모사하였다. 그리고 Geometry 제거로 인해 노출된 영역의 모형지반 절점에 강제변위를 적용하여 토사유출을 모사하였으며, 해석 결과로부터 지반의 공동부와 이완영역을 분석하였다.
그러나 보다 엄밀히 실제 지반의 거동을 분석하기 위해서는 탄소성 이론에 근거하여 해석함이 바람직하며, 더욱이 지중에서의 공동 발생은 지하수위 조건에 영향을 받기 때문에 이러한 제반 조건을 고려한 해석이 필요하다. 따라서 실제 지반에서의 공동 생성 및 주변 지반의 역학거동을 보다 엄밀히 분석하기 위하여, 지하수위 조건에 따라 수행된 기존의 실험적 연구(Kuwano et al., 2010)를 참고하여 탄소성 모델을 적용한 유한요소해석을 추가적으로 실시하였다.
본 연구에서는 하수관 상부 손상 시, 토사유출에 의한 지중 공동 및 이완영역의 발생 메커니즘을 분석하기 위하여, 강제변위법을 활용한 유한요소해석을 수행하여 공동 및 이완영역을 모사하였으며, 기존 연구결과와의 비교를 통해 이를 검증하였다. 또한 공동 주변 지반의 전단응력감소비 분포 특성을 분석하여 공동 및 이완영역의 범위를 정량적으로 제시하였으며, 연구결과를 요약하면 다음과 같다.
1). 또한 트랩도어를 통한 토사 유출 현상을 안정적으로 모사하기 위해 모형지반 하단부에 강재 요소를 추가 적용하였으며, 강재 중앙부의 Geometry를 제거하여 트랩도어 개방을 모사하였다. 그리고 Geometry 제거로 인해 노출된 영역의 모형지반 절점에 강제변위를 적용하여 토사유출을 모사하였으며, 해석 결과로부터 지반의 공동부와 이완영역을 분석하였다.
탄성거동 해석에서 사용되는 2가지 파라미터는 유효 영률(E′)과 유효 포아송 비(ν′)이며 식 (1)와 같이 지반의 변형을 나타낸다. 또한, 탄성 거동 모델을 적용한 지반의 물성은 개별요소 해석시 사용된 값으로 결정하였다(Table 1).
상기의 기존 연구결과를 바탕으로 동일한 조건에서의 유한요소해석을 수행하였다. 수치해석에서는 상기의 실내모형실험과 동일한 경계조건으로 지반을 모델링하였으며, Fig.
수치해석에서는 상기의 실내모형실험과 동일한 경계조건으로 지반을 모델링하였으며, Fig. 6에 나타낸 바와 같이, 지하수위 조건을 단계(0cm → 10cm → 20cm → 40cm)별로 적용하고, 각 지하수위 단계에서 강제변위를 적용하여 공동의 생성과 확장을 모사하였다.
(2010)이 수행한 실내모형실험 결과를 이용해 비교·검증 하였다. 유한요소 수치해석 결과로부터 토사유실로 인한 공동 주변 지반의 간극비 분포 특성을 파악하였으며, 공동 주변 지반의 전단응력 감소비 분포 특성을 분석하여 이완영역의 범위를 정량적으로 제시하였다.
본 연구에서는 지중 공동 주변 지반의 역학거동을 분석하기 위해 유한요소해석을 실시하였다. 이 때, 본 수치해석 결과의 신뢰성 검증을 위해 실내모형실험과 개별요소 수치해석을 수행한 선행연구(Kim et al., 2017)와 동일한 규모와 경계조건으로 모형지반 및 트랩도어를 모델링하였다(Fig. 1). 또한 트랩도어를 통한 토사 유출 현상을 안정적으로 모사하기 위해 모형지반 하단부에 강재 요소를 추가 적용하였으며, 강재 중앙부의 Geometry를 제거하여 트랩도어 개방을 모사하였다.
지하수위 조건에 따라 공동을 모사한 모형실험결과는 공동의 규모 및 형상을 확인할 수는 있지만, 공동 생성에 의해 주변 지반의 응력 재분배로 인한 이완영역 평가는 어려운 일이다. 이에 유한요소해석 결과를 이용하여 공동 주변 지반의 이완영역을 평가하고자, 전단응력 감소비를 통해 주변 지반에 나타나는 이완영역의 응력거동을 분석하였다. 이 때, 전단응력 감소비는 초기 전단응력과 간극비증가로 인해 감소한 전단응력의 비로 정의하였다.
4에서 보는 바와 같이, Toyoura 모래를 이용하여 폭 160cm, 높이 110cm의 모형지반을 제작하였으며. 하수관 파손에 따른 인위적 공동 발생 유도를 위해 토조하단 중심부에 폭 5mm의 트랩도어를 형성시켰다. 이 때, 지하수위는 10cm, 20cm, 40cm의 순으로 상승시키며 지중 공동을 생성하였다.
대상 데이터
(2010)는 지하수위 변화에 따라 토사 유출로 인해 발생하는 공동 및 지반거동을 평가하기 위하여 공동 발생에 대한 실내모형실험을 수행하였으며, 지하수 위 변화에 따른 공동의 형상과 규모를 분석하였다. 모형실험은 Fig. 4에서 보는 바와 같이, Toyoura 모래를 이용하여 폭 160cm, 높이 110cm의 모형지반을 제작하였으며. 하수관 파손에 따른 인위적 공동 발생 유도를 위해 토조하단 중심부에 폭 5mm의 트랩도어를 형성시켰다.
데이터처리
강제변위법을 활용한 지중의 공동과 이완영역은 기존에 필자들이 수행한 연구 결과(Kim et al., 2017)를 활용하여 설정하였으며, 수치해석 결과는 Kuwano et al. (2010)이 수행한 실내모형실험 결과를 이용해 비교·검증 하였다.
이론/모형
수치해석에서 모형 지반의 구성모델은 선행연구의 개별요소 해석에서 적용한 탄성거동을 모사하는 Linear Elastic 모델을 적용하였다. 탄성거동 해석에서 사용되는 2가지 파라미터는 유효 영률(E′)과 유효 포아송 비(ν′)이며 식 (1)와 같이 지반의 변형을 나타낸다.
본 논문에서는 하수관 상부 손상 시, 토사유출에 의한 지중 공동 및 이완영역의 발생 메커니즘을 분석하기 위하여 유한요소 수치해석을 실시하였다. 수치해석에서는 지중의 공동과 함께 이완영역을 모사하기 위해 강제변위법을 적용하였다. 강제변위법을 활용한 지중의 공동과 이완영역은 기존에 필자들이 수행한 연구 결과(Kim et al.
6에 나타낸 바와 같이, 지하수위 조건을 단계(0cm → 10cm → 20cm → 40cm)별로 적용하고, 각 지하수위 단계에서 강제변위를 적용하여 공동의 생성과 확장을 모사하였다. 이 때, 지반 모델은 탄소성 거동을 모사하기 위해 Hardening soil 모델을 적용하였다. Hardening soil 모델은 기존 탄소성 모델에 비하여 소성거동을 세밀히 모사할 수 있다고 알려져 있으며, 흙의 다일러턴시 효과를 고려하고 있으므로 공동 생성에 따른 주변 지반의 거동을 보다 엄밀히 나타낼 수 있다.
유한요소해석을 통한 지중 공동 생성에 대한 모사는 mesh가 연속적인 유한요소법의 특성으로 기인해 직접적인 모사가 불가능하다. 이에 유한요소법을 이용한 공동 모사방법은 임의의 변위를 적용하여 공동을 모사하는 방법이 적합할 수 있기 때문에, 강제변위법을 적용하여 공동과 이완영역을 모사하였다. Fig.
성능/효과
(1) 모형실험 및 개별요소법 기반의 수치해석결과에 의한 간극비 분포도를 유한요소해석 결과와 비교함으로써 유한요소해석의 한계점을 확인하였지만, 공동의 형상 및 규모와 공동 주변지반의 이완영역의 확장정도를 바탕으로, 유한요소해석을 활용한 공동 모사는 충분히 가능한 것으로 판단되었다.
(2) 기존 모형실험 결과를 대상으로 탄소성 유한요소해석을 수행하여 지하수위 조건에 따른 공동과 이완영역을 모사한 결과, 공동생성이 진전됨에 따라 실제 지반에서 생성되는 공동과 이완영역의 규모와 형상이 모형실험 결과와 유사해지는 것을 확인할 수 있었다.
(3) 수치해석결과로부터 강제변위 적용으로 인해 지반의 전단응력이 80% 이상 감소되는 영역은 지반 공동으로 추정할 수 있었으며, 그 경계를 벗어나면서 전단강도 감소비가 점차 감소하며 확장되는 영역을 이완영역으로 추정할 수 있었다.
(2010) 연구결과를 요약하면, Fig. 5(a)에서 보는 바와 같이, 지하수위를 지반 하부로부터 10cm 상승시킨 경우에 토사 유출로 인해 지중에서의 공동이 발생하기 시작하였으며, 공동의 높이는 지하수위보다 상부에서 형성된 것을 확인하였다. 또한 공동의 수평방향으로의 폭은 높이에 비하여 좁게 나타났고, 이후의 지반거동은 토사유출이 더 이상 발생하지 않은 정지상태인 것으로 나타났다.
Fig. 5(b)와 같이 지하수위를 20cm까지 상승시킨 경우, 공동의 높이는 지하수위 10cm일 때와 유사하였지만, 그 폭이 확장되는 것으로 나타났으며, 지하수위를 40cm까지 상승시켰을 경우에는 공동의 높이가 지하수위 높이만큼 연직으로 확장됨과 동시에, 폭 또한 2배 이상으로 확장되는 것을 확인하였다(Fig. 5(c)).
11에서 그 경계를 결정할 수 있었다. 따라서 금번 수행한 유한요소해석의 결과 분석에서도 지반 공동과 이완 영역을 각각 0.52와 0.11의 간극비를 경계로 결정하였다. 유한요소해석에서 모사한 간극비 분포도(Fig.
또한 이 경계를 벗어난 부분에서 전단강도 감소비가 점차 감소하며 확장되는 형상을 볼 수 있으며 이러한 부분을 이완영역으로 추정할 수 있다. 따라서 본 수치해석의 결과를 통해 지반의 공동과 이완영역의 경계를 추정할 수 있으며, 이완영역의 지반 강도 저하 정도를 정량적으로 산정할 수 있다고 판단된다.
5cm보다는 다소 낮게 나타났다. 또한 공동의 폭은 16.5cm로 실험결과인 10.7cm보다 다소 넓게 형성되는 것으로 확인되었다. 지하수위가 20cm인 경우에는 공동의 높이와 폭은 각각 24.
7에 나타내었다. 본 수치해석에서 지하수위 조건(10cm, 20cm 및 40cm)에 따른 최대 강제변위량은 각각 20cm, 32cm 및 70cm인 것으로 나타남을 알 수 있다.
, 2017)에서 설정한 육안 공동 영역(Cavity relative to total area)과 동일한 개념으로써, 전체 모형지반의 면적 대비 분석된 공동의 면적비로 정의하였다. 유한요소해석을 통한 공동비는 모형 실험에 비하여 지하수위 조건에 따라 각각 약 33%, 31% 및 13%가 증가한 것으로 확인되었으며, 이는 공동생성이 진전됨에 따라 유한요소해석에 의한 공동모사가 실제 지반에서 생성되는 공동과 유사해지는 것을 의미한다. 즉, 지반의 공동이 발생되는 초기단계에서는 유한요소해석 결과가 공동비를 다소 과다하게 평가하는 것으로 나타났으나, 지하수위 상승과 함께 지속적인 토사 유출로 인하여 공동이 확장되는 상태에 대해서는 유한요소해석으로 공동을 적절히 모사할 수 있다고 판단된다.
또한 이완영역에 대한 간극비 분포 양상은 기존 연구결과와 다소 차이를 보이고 있지만, 전반적인 형상과 규모는 유사하게 나타나고 있음을 알 수 있다. 이러한 결과로부터 본 연구에서 수행된 강제변위법에 의한 지반 공동과 이완영역의 모사는 그 형상과 규모, 확장성 등을 고려해 볼 때 적절히 수행되어졌다고 판단할 수 있다.
유한요소해석을 통한 공동비는 모형 실험에 비하여 지하수위 조건에 따라 각각 약 33%, 31% 및 13%가 증가한 것으로 확인되었으며, 이는 공동생성이 진전됨에 따라 유한요소해석에 의한 공동모사가 실제 지반에서 생성되는 공동과 유사해지는 것을 의미한다. 즉, 지반의 공동이 발생되는 초기단계에서는 유한요소해석 결과가 공동비를 다소 과다하게 평가하는 것으로 나타났으나, 지하수위 상승과 함께 지속적인 토사 유출로 인하여 공동이 확장되는 상태에 대해서는 유한요소해석으로 공동을 적절히 모사할 수 있다고 판단된다.
지하수 위가 10cm 및 20cm인 경우에는 이완영역이 유사한 것으로 확인되었지만, 지하수위가 40cm일 때에는 이완영역이 크게 증가하는 것으로 나타났다. 즉, 지하수위가 낮은 경우보다 높게 유지되는 경우에 이완영역의 확장이 더욱 크게 발생하는 것으로 분석되었다. 이는 파손된 하수관 상부에 지하수위가 높게 형성되는 경우에 이완영역의 규모가 증가할 수 있음을 시사한다.
, 2017). 지하수 위가 10cm 및 20cm인 경우에는 이완영역이 유사한 것으로 확인되었지만, 지하수위가 40cm일 때에는 이완영역이 크게 증가하는 것으로 나타났다. 즉, 지하수위가 낮은 경우보다 높게 유지되는 경우에 이완영역의 확장이 더욱 크게 발생하는 것으로 분석되었다.
후속연구
(4) 유한요소 수치해석을 통한 전단응력 감소비를 활용하여 공동 및 주변 지반의 이완영역 경계를 추정할 수 있으며, 이완영역의 지반 강도 저하 정도를 정량적으로 산정하여 공동발생에 따른 이완영역 예측 및 그에 대한 합리적인 보강 방안 수립이 가능할 것으로 판단된다.
질의응답
핵심어
질문
논문에서 추출한 답변
도심지 지반함몰의 원인은?
최근 국내외적으로 사회적 이슈가 되고 있는 도심지 지반함몰은 기존 매설되어 있는 상하수도관의 손상이나 지중구조물 시공 중 수행되는 굴착이나 되메움 시 시공관리 부실로 인한 토사유출이 중요 원인으로 분석되고 있다(Park Park, 2014; Choi et al., 2016).
최근 5년간 서울시에서 발생한 지반함몰 건수는?
, 2016). 특히 국내에서는 인구밀도가 높고 다수의 노후한 상하수도관이나 지하매설물 등이 분포하는 서울시에서 매년 크고 작은 지반함몰 발생이 증가하고 있으며, 최근 5년간 약 3200여건에 달하는 지반함몰 발생이 보고되고 있는 실정이다(한국지질자원연구원, 2014). 또한 지반함몰 발생 건수 중, 하수관 손상에 의한 지반함몰이 전체 발생 건수의 약 85%(2,636건)를 차지하고 있으며 이는 지중매설관의 손상이 지반함몰 발생빈도와 매우 밀접한 연관성을 갖는 것으로 분석되고 있다(Choi et al.
지반함몰 사고를 예방하기 위해서는 지반 내 공동 방지 대책 마련과 함께 이를 위한 공동 발생의 원인 및 메커니즘 규명에 대한 연구가 필요한 이유는?
특히 국내에서는 인구밀도가 높고 다수의 노후한 상하수도관이나 지하매설물 등이 분포하는 서울시에서 매년 크고 작은 지반함몰 발생이 증가하고 있으며, 최근 5년간 약 3200여건에 달하는 지반함몰 발생이 보고되고 있는 실정이다(한국지질자원연구원, 2014). 또한 지반함몰 발생 건수 중, 하수관 손상에 의한 지반함몰이 전체 발생 건수의 약 85%(2,636건)를 차지하고 있으며 이는 지중매설관의 손상이 지반함몰 발생빈도와 매우 밀접한 연관성을 갖는 것으로 분석되고 있다(Choi et al., 2016).
참고문헌 (10)
Brady, B. H. G. and Brown, E. T. (1985), Rock mechanics for underground mining, George Allen and Unwin., pp.527.
Choi S. K., Back S. I., An J. B., and Kwon T. I. (2016), "Geotechnical investigation on causes and mitigation of ground subsidence during underground structure construction", Journal of Korean Tunnelling and Underground Space Assoication, Vol.18, No.2, pp.143-154.
Kim, J. B., You, S. K., Han, J. G., Hong, G. G. and Park, J. B. (2017), "A Study on Simulation of Cavity and Relaxation Zone Using Laboratory Model Test and Discrete Element Method", J. Korean Geosynthetics Society, Vol.16, No.2, pp.11-21.
Kim, N. Y. and Umm, T. W. (2013), "A case study on the Chimney collapse of tunnel under construction", Proceeding of Korean geo-environmental society, Seoul, pp.43-53.
Korea institute of geoscience and mineral resources (2014), Research on causes and policy suggestions by sinkhole type, Research report, pp.18-39.
Kuwano, R., Sato, M., and Sera, R. (2010), "Study on the detection of underground cavity and ground loosening for the prevention of ground cave-in accident", Japanese Geotechnical Journal Vol.5, No.2, pp.219-229.
Lee, S. H., Lee H. L., and Song, K. I. (2016), "The effect of formation of spherical underground cavity on ground surface settlement : Numerical analysis using 3D DEM", Journal of Korean Tunnelling and Underground Space Association, Vol.18, No.2, pp.129-142.
Park, I. J. and Park, S. H. (2014), "Cause analysis and counterplan for sinkhole", Magazine of Korean Society of Hazard Mitigation, Vol.14, No.5, pp.12-17.
Song, K. I. and Yoon, J. S. (2015), Tunnel deformation mechanism, CIR, Seoul.
Suchowerska, A. M., Merifield, R. S., Carter, J. P., and Clausen, J. (2012), "Prediction of underground cavity roof collapse using the Hoek-Brown failure criterion", Computers and Geotechnics, Vol.44, pp.325-342.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.