• 지질공학
• /
• 제32권4호
• /
• pp.643-659
• /
• 2022

최근 지중저장기술(예, 온실가스 심지층 처분, 인공지열저류층 발전 등)이 활발히 수행됨에 따라, 유체 주입과 저장부지 안정성 사이의 역학적 관계에 관한 정량적 이해의 중요성이 인지되고 있다. 지중 유체 주입은 공극압 및 지중응력 교란과 지층의 역학적 불안정성을 야기할 수 있어, 유체 주입에 대한 다공탄성 수치 모형 구축이 요구된다. 본 연구에서는 순차적인 COMSOL-PyLith-COMSOL 유체 주입-유발지진 다공탄성 수치 모사를 수행한다. 유한요소 상용 소프트웨어인 COMSOL을 이용해 단층에 가해지는 쿨롱 파괴 응력(CFS) 변화를 시간에 따라 추적하였고, CFS 변화량이 임계값(예, 0.1 MPa)을 초과할 경우, 모형의 정보(기하구조, 물성 등)를 유한요소 오픈소스 코드인 PyLith로 이동시키는 알고리즘을 구축했다. PyLith는 단층의 미끄러짐을 모사하고, 미끌림에 의한 변위장을 획득한다. 이후 변위장을 COMSOL로 이동시켜 지진에 의한 응력 및 표면 변위를 계산한다. 수치 모사 결과, 주입 기간 중엔 주입정 근거리에서 큰 변화(공극압, CFS 변화 등)를 보였고, 주입 종류 후에는 잔류 응력이 원거리 영역으로 확산하는 양상이 나타났다. 이는 주입 종료 후 지속적인 모니터링의 필요성을 제안한다. 또한, 단층과 주입층 물성(예, 투수계수, Biot-Willis 계수)에 따른 CFS 변화량 비교는 주입정 위치 선정 시 주입층 및 주변 지층에 대한 물성 파악이 중요함을 의미한다. 단층 미끄러짐 양에 따른 표면 변위 및 이암층에 가해지는 편차응력은 다양한 단층 미끌림 시나리오 설정의 필요성을 지시한다.

• 한국강구조학회 논문집
• /
• 제9권2호통권31호
• /
• pp.171-179
• /
• 1997

고장력 볼트 이음부에서의 힘의 전달은 부재간의 압축력 및 볼트축력에 의하여 외력과 균형을 이룸으로서 형성되어지며, 부재간의 압축력, 도입된 축력의 변화 및 미끄러짐하중 등은 볼트 체결부 및 볼트의 유효강성에 지배된다. 이때 유효강성은 일반적으로 접합부의 유효단면적과 탄성계수의 곱으로 나타내어지지만 접합부재의 유효단면적이 어느 정도가 되는지는 판단하기가 쉽지가 않다. 이를 위해서 종래에는 이음부의 탄성변형에 관계하는 유효단면적을 등가중공원통형으로 가정하여 이에 대한 여러 가지 검토가 이루어졌다. 그러나 이러한 제안식들은 설계상의 목적과 복잡한 계산을 피하기 위하여 매우 단순화시킨 것으로 어느 식이 타당한 지는 단정하기가 어려우며 이에 대한 상세한 해석적인 검토 및 실험적인 검증이 요구된다. 따라서, 본 연구에서는 볼트 축력에 의한 접합부재의 유효단면적 산정과 이에 대한 검증을 위하여 고장력 볼트 이음 시험편에 대해 피로시험을 실시하고 반복하중 재하 후의 마찰면의 형상을 관찰하여 유효마찰 영역을 측정하였으며, 유한 요소법에 의해 수치 해석을 실시하여 내부 압축응력의 분포형상을 밝히고 접촉면상에 발생하는 응력 분포영역을 근사화하므로서 마찰면의 형상과 내부 압축응력의 분포영역과의 관계 그리고 유효단면적의 산정 방법을 고찰하였다.

• 한국지반공학회논문집
• /
• 제24권5호
• /
• pp.65-78
• /
• 2008

절토사면 및 천심도 터널의 경우 불연속면을 따른 블록의 미끄러짐 및 회전 등이 안정성에 큰 영향을 미친다. 절토사면이나 터널의 초기 설계 단계에서 안정성 검토를 위해 수치해석을 널리 사용하는데, 대부분의 상용화된 프로그램들은 연속체 해석에 기반을 두고 있다. 불연속면을 고려하여 연속체 해석을 수행할 수도 있으며, 절리를 포함한 암반자체의 강도를 감소시킴으로써 모사하는 것이 대표적인 방법이다. 또한, 연속체 해석방법인 Mohr-Coulomb 모델에 편재 절리모델 (ubiquitous joint model)을 적용함으로써 불연속면을 고려하는 방법이 있다. 절토사면의 수치해석시 주 입력치인 절리암반의 점착력, 마찰각, 탄성계수 등은 암반분류법(rock classification system)에 의한 경험식으로 도출할 수 있으며, 본 논문에서는 이 중 RMR 및 GSI 시스템에 의한 경험식을 통해 얻은 입력치를 실제 위험절토사면에 적용함으로써 두 시스템의 차이를 비교해 보고자 하였다. 또한 RMR 및 GSI 시스템으로 도출한 입력치에 편재절리모델을 혼합하여 적용함으로써 계측을 통해 얻은 변위량 및 변위양상과의 차이를 분석하였다. 분석 결과 GSI 시스템으로 도출한 입력치에 편재절리모델을 혼합하였을 경우의 수치해석결과가 실제현장의 변위량 및 변위양상과 가장 유사한 결과를 얻을 수 있음을 확인할 수 있었다.