• 방사성폐기물학회지
• /
• 제10권4호
• /
• pp.281-294
• /
• 2012

현재까지 개발된 고준위폐기물 심지층처분장의 열-수리-역학적 복합거동 해석을 위한 전산 코드의 현황을 조사하고, 문헌에 보고된 각 코드에 의한 계산치와 현장실험 측정치의 비교 결과를 이용하여, 기존 전산 코드들의 신뢰도를 분석하였다. 개발된 전산코드들은 완충재가 없는 처분장에서는 붕괴열에 따른 암반의 열-수리-역학적 거동을 비교적 잘 모사하였으나, 포화 경암층에 위치한 완충재가 존재하는 처분장의 공학적방벽시스템 내에서 일어나는 열-수리-역학적 복합거동의 예측은 만족스럽지 못하였다. 현재 제안된 열-수리-역학적 복합거동 해석모델을 고준위폐기물 처분장 공학적방벽시스템의 거동 해석에 적용하기 위해서는 완충재 내의 수분함량 및 전 압력 분포를 보다 정교하게 모사할 수 있도록 수학적 모델의 개선이 필요하다.

• 지질공학
• /
• 제11권3호
• /
• pp.255-267
• /
• 2001

지하 50m의 화강암반내 단층지역에 위치한 지하 방사성폐기물 처분장 인접지역에서의 열, 수리, 및 역학적 연성거동을 비교하고 분석하였다. 해석에는 2차원 해석코드인 UDEC을 사용하였다. 해석모델은 화강암반, 처분공내의 압축 벤토나이트로 둘러싸인 PWR 사용후 핵연료 처분용기, 및 처분동굴내에 채워진 혼합 벤토나이트를 포함한다. 수리-역학적, 열-역학적, 및 열-수리-역학적 연성거동을 비교 및 분석하였다. PWR 사용후 핵연료내의 방사성 물질로부터 나오는 시간의존 방사성 붕괴열이 처분장 및 인접지역에 미치는 영향을 분석하였다. 수리해석에는 steady state flow 알고리즘을 사용하였다.

• 방사성폐기물학회지
• /
• 제4권2호
• /
• pp.117-131
• /
• 2006

한국의 고준위폐기물 기준 처분 시스템의 공학적 방벽에서의 T-H-M(Thermo-Hydro-Mechanical) 거동 실증을 위한 KENTEX(KAERI Engineering-scale T-H-M Experiment for Engineered Barrier System)실험 장치를 대상으로 열-수리-역학 연동현상 해석을 하여 온도, 포화도 및 응력의 변화를 예측하였다. 그리고 이들 변수와 열-수리-역학의 연동현상에 사용된 세물성법칙인 탄성물성법칙, 공극탄성 물성법칙 및 공극탄성-소성 물성법칙과의 관계를 분석하였다. 열-수리-역학 연동현상을 계산하는 데는 상용 유한요소 코드인 ABAQUS를 사용하였다. 열 계산에서 벤토나이트 내 온도는 히터 가열 후 초기에는 급격히 증가하다가 얼마의 시간이 경과한 후에는 거의 일정한 값에 도달하였다. 이 도달시간은 약 37.5일로 반경방향의 모든 지점(H=0.68m 일때)에서 정상상태에 도달한 것을 알 수 있었다. 즉, 히터와 벤토나이트 경계면에서는 $90^{\circ}C$, 벤토나이트와 외부 셀 경계면에서는 약 $70^{\circ}C$를 유지하였다. 열-수리-역학 연동현상 계산에서 시간에 따른 벤토나이트 포화도는 탄성 물성법칙, 공극탄성 물성법칙 및 공극탄성-소성 물성법칙의 세 경우 모두 거의 차이가 없었다. 열-수리-역학 계산 결과와 수리-역학 계산 결과의 비교에서 온도의 증가는 탄성 물성법칙 및 공극탄성 물성법칙 각각에 대해 시간이 경과함에 따라 포화도가 증가함을 초래해 포화가 빨리 진행됨을 알 수 있었다. 특히 히터에 가까운 쪽에서는물이 침투하고 있는 쪽 보다 포화도 증가가 큰 것으로 나타나 벤토나이트가 물로 포화되기 전의초기상태가 온도의 영향을 많이 받는 것을 알 수 있었다. 또한 응력은 세 물성 법칙 모두 시간의 경과에 따라 증가하는 경향을 보이나 탄성 물성법칙의 경우가 다른 두 경우보다 현저한 변화를 보이는데 이는 변형율이 탄성한계를 넘어서도 계속 작용하여 공극비 변화를 고려한 다른 두 물성법칙과 차이가 있음을 나타내고 있다. 그러나 공극탄성 물성법칙 및 공극탄성-소성 물성법칙의 경우에 열-수리-역학 계산 결과와 수리-역학 계산 결과를 비교하면 시간이 경과함에 따라 응력은 증가하지만 온도의 변화에 따른 서로의 응력의 차이는 작은 것을 알 수 있다. 즉 온도변화의 영향보다는 시간에 따른 포화도 변화의 영향이 더 큰 것으로 생각된다. 따라서 벤토나이트의 열-수리-역학 연동현상 해석에서 벤토나이트는 온도의 증가로 포화가 빨라지고, 포화도 증가는 응력을 증가시키는 결과를 보이므로 공극비, 열팽창 및 팽윤압 등의 영향을 받고 있는 것으로 이해된다. 그래서 벤토나이트의 열-수리-역학 연동현상 해석에서 벤토나이트는 공극비, 열팽창 및 팽윤압 등의 영향을 받으므로 탄성과 소성을 동시에 고려할 수 있는 물성법칙을 선택하는 것이 바람직하다.

• 터널과지하공간
• /
• 제21권1호
• /
• pp.66-81
• /
• 2011

암반 내 균열 생성, 진전, 파괴 등과 같은 지하 암반의 역학적 거동과 이들 균열을 통한 지하수 유동 및 온도 변화에 기인한 열응력이 역학적 거동에 미치는 상호작용을 모델링하기 위한 수치해석코드를 개발하였다. 개발된 수치해석코드에서는 기존의 2차원 FRACOD(Shen & Stephasson, 1993)에 열-역학 및 수리-역학 상호 거동을 모델링하기 위한 해석모듈을 개발하여 추가하였다. 열-역학 연계를 위해서는 가상열원법과 시간전진기법을 도입하였으며, 수리-역학 연계에서는 양해법에 의한 반복기법을 적용하였다. 수치해석결과와 해석해와의 비교를 통해 개발된 해석모듈의 유용성을 검증하고 해석사례를 통해 온도변화에 따른 열균열 발생 및 수압파쇄과정에서의 균열 진전 양상과 같은 절리암반 복합거동을 잘 표현할 수 있음을 확인하였다.

• 한국지반환경공학회 논문집
• /
• 제10권2호
• /
• pp.61-68
• /
• 2009

본 논문은 수리학적 현장 주입시험 기간 동안 고려될 필요가 있는 수리역학적 상호작용에 따른 균열암반매질의 수리학적 거동에 대한 수치적 연구이다. 이러한 주입시험은 굴착정 내에 설치된 개구간(open hole section)을 따라 높은 압력을 가진 유량을 주입하며, 이를 통해 굴착정을 가로지르는 균열로부터 유량을 측정하는 것이다. 시간에 따라 변화되는 유량측정결과는 수리역학적 상호작용에 대한 분석 및 예측을 위해 개발된 수치모형의 현장 적용성을 평가하기 위해 사용되었다. 유량측정결과 전도성이 있는 균열들은 상호의존적인 균열망을 형성하며, 이로 인해 균열망을 구성하는 개별적인 균열요소들은 독립적인 시스템으로서 적절하게 모형화 될 수 없음을 보여주었다. 또한 간극수압이 굴착정 주위에 작용하는 최소주응력을 초과할 때 새로운 유체유입영역이 발생되며, 이러한 최소주응력보다 훨씬 큰 간극수압은 굴착정 주위의 균열들에 의해 유지될 수 있다는 것을 보여주었다. 본 연구에서 이러한 특성들이 자연 상태의 균열망의 구성형태에 따라 어떻게 영향을 미치게 되는지 이산 균열망 모형을 통해 수치적으로 분석되었다.

• 터널과지하공간
• /
• 제33권4호
• /
• pp.265-280
• /
• 2023

터널누수는 주변지반의 응력 및 간극수압을 변화시켜 터널 안정성 및 지반변형에 영향을 미칠 수 있는 결함요소이다. 장기간 또는 큰 규모의 누수발생은 터널 라이닝의 불안정성 및 지표침하와 같은 터널 구조물 및 주변지반 환경에 손상을 일으킬 수 있다. 본 연구에서는 누수발생 시 터널의 구조 안정성 및 지반거동에 미치는 영향을 수치해석적으로 분석하였다. 고려된 터널은 내부로 주변 지하수의 유입을 허용하지 않는 비배수 조건으로 가정하였고 터널 완공 후 라이닝에서 누수가 발생하는 것으로 설정하였다. 누수로 인한 터널 구조물 및 지반의 거동을 모사하기 위해 수리역학 연계해석이 수행되었으며 파이썬으로 개발된 TOUGH-FLAC 시뮬레이터가 사용되었다. 누수 발생량과 누수위치를 변화시켜 수치모사가 수행되었으며 수리역학 해석을 위한 연계항들이 복합거동 결과에 미치는 영향을 조사하였다.

• 터널과지하공간
• /
• 제17권2호
• /
• pp.109-118
• /
• 2007

굴착손상영역(EDZ)은 굴착으로 인해 현지 암반이 역학적으로 손상을 입게 되어 응력상태, 변위상태, 암반의 안정성, 지하수의 흐름상태 등에 변화가 일어나는 영역을 의미한다. EDZ의 역학적 특성과 관련한 많은 연구들이 수행되었지만, EDZ에서의 지하수 유동 특성에 관한 연구는 아직 부족한 수준이다. 본 연구에서는 굴착으로 인한 수리-역학 상호작용(coupling)에 의해 굴착면 주변의 수리적 간극값이 변하는 영역을 수리적 굴착 thstkdduddudr이라 정의하고, 이를 3차원 분리단열망(DFN)에 적용시켜 보았다. 이를 통해 수리적 간극변화가 3차원 불연속 망에서의 전반적인 지하수 유동에 미치는 영향을 파악하였다. 또한 3차원 DFN 지하수 유동 해석 시 주로 이용되는 수두 조건과 유량 조건의 적용성을 고찰하였다. 해석 결과 수리-역학적 상호거동에 의해 발생하는 굴착면 주변의 수리적 간극변화는 터널 내부로 유입되는 유량에 큰 영향을 미치는 것으로 나타났다. 또한 DFN 해석 시 다양한 경계조건에 따른 상이한 결과를 토대로 보다 합리적인 경계조건 설정에 대한 방향을 제시하였다. 마지막으로 실제 현장에서 수리해석을 실시한 자료를 바탕으로 수리적 간극 값의 변화를 고려할 때와 고려하지 않을 때의 유입유량 차이를 통해, 3차원 지하공동의 지하수 유동해석 시 수리적 간극 값의 변화를 고려하는 것이 보다 더 보수적인 결과를 나타내는 것을 확인하였다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2021년도 학술발표회
• /
• pp.238-238
• /
• 2021

토목 구조물에서 수문 (Gate)은 다양한 시설에서 설계되어 적용되고 있다. 이런 수문은 용도나규모, 사용재료, 구동방식에 의해 다양하게 분류되는데, 본 글에서 다룰 Girder 타입의 Sluice Gate 는 가장 많은 분야에서 포괄적으로 적용되는 수문 형식중 하나이다. 일반적으로 2m×2m 미만의 소규모 수문에서는 수문설계에 수리계산이나 응력검토를 요하지 않으나, 높이 5m 가 넘는 대형 수문의 경우 비체의 무게가 수십톤에 달하기 때문에 정수역학이나 동수역학과 같은 수리학적 설계가 반드시 필요하다. 수문설계는 크게 세 가지 해석을 다루는데, 첫 번째는 정수역학 또는 동수역학적 거동에 대한 수문 비체의 구조해석이고 두 번째는 수문의 상승 및 하강을 위한 권양장치 및 수문비체, 구동장치 등의 계산이며 마지막으로는 수문의 형상과 수류형상에 따른 수문의 진동해석 이다. 본 글에서는 두 번째에 해당하는 수문의 상승 및 하강에 영향을 미치는 수리학적 요소들을 통해 가장 합리적인 수문형상설계에 대해 논하고자 한다. 특히 국내에서는 Girder 타입의 수문설계에서 수문하단부 sill 부근의 형상을 쐐기 형태로 일률적으로 적용하고 있는데, 이를 개선하여 유체역학적으로나 경제적으로 보다 유리한 설계안을 본 글을 통해 제안하고자 한다.

• 터널과지하공간
• /
• 제29권1호
• /
• pp.12-29
• /
• 2019

최근의 이산화탄소 심지층 처분, 인공지열저류층 형성 및 발전, 원유회수증진 사업 등에는 고압의 유체를 주입하는 과정이 수반되고 이들 사업의 안전하고 성공적인 수행을 위해서는 대상 부지 절리 암반의 투수 특성 및 주입압에 의한 역학적 변형에 기인한 변화를 정확하게 평가하는 것이 중요하다. 본 고에서는 절리 암반 수리-역학 특성 파악을 위한 해석적 및 실험적 평가 방법에 대해 검토하였다. 먼저 유체 주입 전 절리 암반 초기 투과 특성 및 수리-역학 특성 평가에 고려해야 할 기술요소를 분석하고 현장 시추공 실험을 통해 이들 수리-역학 특성을 직접 측정하기 위한 최근의 SIMFIP 실험장치의 특징 및 활용방안에 대해서도 검토하였다.

• 한국터널지하공간학회 논문집
• /
• 제10권4호
• /
• pp.397-404
• /
• 2008

본 연구는 인천국제공항철도 신설 공사구간 중 서울지하철 2호선 하부에서 TRcM 공법으로 통과하는 구간의 수리 및 굴착 안정성에 대해 검토한 내용이다. 본 검토 대상 구간에서는 터널 굴착과정에서 하루 최대 86.4 ton의 예상치 못한 과다 용출수가 발생하였으며, 과다 용출수의 발생에 대한 원인 분석을 위해 관측정에서 계측된 지하수위 저하거동과 가장 적합한 수리해석 파라미터를 찾기 위한 수리 역해석을 실시하였다. 본 수리 역해석에서는 TRcM 구간에 접한 홍익대학교 정거장 굴착시의 측정된 자료중 TRcM 구간 굴착전까지의 지하수위 계측결과를 사용하였다. 암종별 수리전도계수와 서울지하철 2호선의 배수경계조건, 과업구간의 강우강도 등을 변화시켜가며 수리 역해석을 수행하였다. 그 결과 계측된 자료와 유사한 최종 지하수위 및 지하수위 저하 경향을 보이는 수리해석 조건을 얻을 수 있었다. 다음, 지반의 굴착에 의한 역학적 영향을 분석하기 위해 지하수를 고려하지 않은 역학적 해석을 수행하였다. 순수 역학적 해석 결과, 지표침하 0.85 mm, TRcM 구간 천단침하 1.32 mm 발생하는 것으로 나타났으며, 복합거동 해석결과 지표침하 1.2 mm, TRcM 구간 최대 천단침하 1.72 mm으로 나타나 지하수위 저하에 의해 추가 변위가 발생하나 침하량이 크지 않아 안정성에는 영향이 없는 것으로 분석되었다.