📄 other

Multiscale Calculation of Equivalent Elastic Modulus for Rock Mass Based on Perturbation Theory

본 논문은 암반의 등가 탄성 계수를 정확하게 결정하기 위해 초기 손상과 미세 구조적 불균질성을 포함하는 섭동 이론 기반의 다규모 계산 접근법을 제안하며, 실험실 및 현장 시험 결과와의 검증을 통해 기존 방식보다 우수한 정확도를 입증한다.

원저자: Lei Wen, Wu Yang, Huayun Yang

게시일 2026-07-24

📖 3 분 읽기☕ 가벼운 읽기

CC BY 4.0

원저자: Lei Wen, Wu Yang, Huayun Yang

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

거대한 돌무더기로 된 무너져가는 성이 압력을 받을 때 어떻게 버틸지 예측하려고 한다고 상상해 보십시오. 하지만 여기 반전이 있습니다. 이 성은 매끄러운 대리석 덩어리가 아닙니다. 그것은 매끄러운 것과 거친 것이 섞인 돌들이 시멘트 조각들로 연결되어 있고, 그 사이사이로 거미줄 같은 미세한 균열들이 숨겨져 있는 지저질한 돌무더기입니다. 이것이 바로 '암반(rock mass)' 공학의 실제 세계입니다. 엔지니어들이 터널, 광산, 또는 댐을 건설할 때, 실험실에서 테스트한 완벽하고 작은 정육면체 암석 하나만을 보고 거대한 산 전체가 똑같이 움직일 것이라고 가정할 수는 없습니다. 그것은 마치 건물 전체가 지진을 견딜 수 있을지 확인하기 위해 벽돌 한 장을 테스트하는 것과 같습니다. 문제는 암석이 매우 복잡하며, 어떻게 휘어지고 부서지는지를 변화시키는 미세한 결함과 균열로 가득 차 있다는 점입니다. 과학자들은 오랫동안 미시적인 균열의 세계와 거대한 암석 형성층의 세계 사이의 간극을 메우려 노력해 왔지만, 그들의 오래된 지도들은 종종 잘못된 길로 인도했습니다. 그들에게는 미시적인 세계와 거시적인 세계를 동시에 볼 수 있는 새로운 방법이 필요했습니다.

"섭동 이론(Perturbation Theory)에 기반한 암반의 등가 탄성 계수에 대한 다규모 계산(Multiscale Calculation of Equivalent Elastic Modulus for Rock Mass Based on Perturbation Theory)"이라는 제목의 이 논문은, 엔지니어들이 울퉁불퉁하고 균열이 있는 암석이 어떻게 찌그러지고 늘어나는지를 정확하게 볼 수 있게 해주는 새로운, 아주 똑똑한 안경과도 같습니다. 저자인 린 원(Lei Wen), 우 양(Wu Yang), 후이윈 양(Huayun Yang)은 암석의 초기 손상(이미 존재하는 미세한 균열들)을 출발점으로 삼는 '다규모 계산'이라는 영리한 수학적 기법을 개발했습니다. 그들은 암석이 완벽하다고 가정하는 대신, '섭동 이론'이라 불리는 방법—이는 기본적으로 "완벽한 시스템을 살짝 건드려 어떻게 흔들리는지 보자"는 방식입니다—을 사용하여 미세한 균열들이 암석의 강성에 어떤 영향을 미치는지 파악합니다. 그들은 이러한 미세한 균열들로 인해 발생하는 변형(squishing)을 미시 세계와 거시 세계를 연결하는 '다리(bridge)'라고 부릅니다.

연구팀은 단순히 컴퓨터 앞에만 앉아 있었던 것이 아니라, 실제 암석을 채취하기 위해 중국 남부의 깊은 금속 광산으로 직접 들어갔습니다. 그들은 암석을 서로 다른 힘(구체적으로 완벽한 암석 강도의 0.25, 0.4, 0.5, 0.65, 0.75배)으로 눌러서 의도적으로 손상을 입힘으로써 '초기 손상'을 만들었습니다. 그런 다음, 이 손상된 암석들을 두 가지 테스트, 즉 '브라질 분할 시험'(오렌지를 가르듯 측면에서 눌러 암석을 쪼개는 시험)과 표준 '일축 압축 시험'(위아래에서 짓누르는 시험)에 통과시켰습니다. 또한 그들은 암석 내부의 보이지 않는 균열을 관찰하기 위해 음파(P파 속도)와 CT 스 l 스캔을 사용했습니다.

그들이 발견한 것은, 그들의 새로운 수학적 방법이 기존의 표준 컴퓨터 시뮬레이션보다 현실을 훨씬 더 잘 예측한다는 사실입니다. 그들이 결과를 실제 실험 결과와 비교했을 때, 그들의 '다규모' 접근 방식은 실제 암석과 거의 완벽하게 일치하는 곡선을 그려냈습니다. 기존의 방법들은 정점 강도(부서지기 전까지 얼마나 세게 밀 수 있는지)를 추측하는 데는 괜찮았지만, 곡선의 형태는 틀리게 그렸습니다. 새로운 방법은 인장(tensile) 하중-변위 곡선(분할 시험 결과)이 실제와 같이 특정한 '엄격한 하부 볼록(strictly lower convex)' 형태로 굽어진다는 것을 보여주었으나, 기존 시뮬레이션은 이 미묘한 차이를 놓쳤습니다. 반면, 압축(compression) 응력-변형률 곡선(짓누르는 시험 결과)은 다른 '볼록(convex)'한 형태를 따릅니다. 그들은 심지어 이 방법을 납-아연 광산이 붕괴되었던 실제 건설 현장에 적용했습니다. 그들은 '그라우팅 고결체(grouting solidified bodies)', 즉 시멘트 슬러리로 다시 붙여진 부서진 암석의 강도를 계산하는 데 이 방법을 사용했습니다. 이 수학 모델은 암석이 파괴점에 도달하기 전의 거동을 정확하게 예측했습니다.

하지만 저자들은 자신들의 방법이 마법이 아니라는 점을 주의 깊게 언급합니다. 그들의 방법은 암석이 부서지기 전의 거동은 완벽하게 잡아냈지만, 암석이 완전히 무너지기 시작하면서 발생하는 '변형 연화(strain-softening)' 현상을 정확히 예측하는 데는 어려움을 겪었습니다. 이는 그들의 모델이 암석이 완전히 파괴될 때 일어나는 지저분한 변형 연화 과정을 아직 반영하지 못하기 때문입니다. 그러나 암반이 굴복하기 전, 암반이 얼마나 단단하고 강한지를 파악하는 가장 중요한 부분에 있어서, 미시 세계와 거시 세계 사이의 이 새로운 '다리'는 기존의 도구들보다 더 효과적입니다. 이는 암석의 초기 흉터를 인정하고 이 섭동 수학을 사용함으로써, 엔지니어들이 암반이 버틸 수 있을지 아니면 붕괴할지를 훨씬 더 명확하게 파악할 수 있으며, 이를 통해 우리의 터널과 광산을 더 안전하게 만들 수 있음을 시사합니다.

기술 요약: 섭동 이론에 기반한 암반의 등가 탄성 계수 다중 스케일 계산

문제 제기
등가 탄성 계수를 정확하게 결정하는 것은 암반 공학의 안전 및 안정성 평가에 있어 매우 중요하다. 현재의 방법론들은 상당한 한계에 직면해 있다. 건실한 암석 시편에 대한 실험실 테스트는 공학적 암반에 내재된 거시적 및 미시적 절리와 균열을 반영하지 못하며, 현장 테스트는 제한된 공간적 범위로 인해 대표성이 떨어지는 경우가 많다. 또한, 실험실 결과를 현장 조건으로 스케일링하기 위해 사용되는 전통적인 저감법은 신뢰할 수 없는 저감 계수에 의존하며, 암반 파괴 과정 중의 역학적 응답 과정을 포착할 수 없다. 따라서 건실한 암석 역학과 손상된 암반의 불균질한 거동 사이의 간극을 메울 수 있는 체계적인 다중 스케일 결합 접근 방식이 필요하다는 인식이 존재한다.

방법론
저자들은 암반의 초기 손상을 고려하기 위해 설계된, 섭동 이론(perturbation theory)에 기반한 다중 스계 계산 접근법을 제안한다. 이 방법론은 암반을 내부 균열의 다중 스케일 확장 과정을 반영하는 작은 매개변수( $\delta$ )를 가진 섭동 시스템으로 취급한다.

주요 방법론적 단계는 다음과 같다:

이론적 프레임워크: 중간 규모의 불균질성(meso-heterogeneity)에 의해 발생하는 변형률 텐서를 거시적 척도와 미시적 척도를 연결하는 "가교"로 활용한다. 변위장은 점근적 전개(asymptotic expansion)를 거치며, 변수를 거시적 좌표( $x$ )와 미시적 근방( $y$ )으로 분리한다.
지배 방정식: 미세 탄성 손상 계수( $D^\delta$ )를 포함하는 구성 방정식과 평형 미분 방정식을 결합하여 지배 방정식을 유도한다. 초기 손상은 초기 탄성 손상 텐서로 표현된다.
균질화(Homogenization): 서로 다른 물리적 차원 스케일 하에서 운동 방정식을 해결함으로써, 저자들은 균질화된 거시 등가 탄성 계수( $E_{mnkl}^\delta$ )를 도출한다. 이는 변분 원리에 기초하여 유한 요소법을 사용하여 텐서 고유함수( $\chi$ )를 구하는 과정을 포함한다.
수치적 구현: 접근 방식은 손상 계수의 반복 계산을 위한 MATLAB과 유한 요소 모델링을 위한 COMSOL을 사용하여 구현되었다. 모델은 진화하는 손상 상태를 반영하기 위해 각 단계에서 하중( $Q$ )과 변위( $u$ )를 업데이트한다.
검증: 본 접근법은 두 가지 주요 사례를 통해 검증되었다:
1. 실험실 검증: 중국 남부의 심부 금속 광산에서 채취한 석회암 시편을 대상으로 일축 압축 시험 및 브라질리언 인장 시험(Brazil split tensile tests)을 수행하였다. 시편은 다양한 정도( $0.25\sigma_{cs}$ 에서 $0.75\sigma_{cs}$ )로 사전 손상을 입혔으며, P파 속도와 CT 스캐닝을 통해 특성을 규명하였다.
2. 현장 검증: 본 방법을 납-아연 광산 스톱(stope)의 그라우팅 고결체의 적용하였다. 이때 그라우팅된 재료는 불규칙한 절리 균열을 가진 암반으로 취급되었다.

주요 기여

다중 스케일 섭동 모델: 본 논문은 중간 규모의 불균질성에서 발생하는 변형률 텐서가 미시적 척도와 거시적 척도를 연결하는 메커니즘 역할을 하여, 초기 손상을 고려한 등가 탄성 계수 계산을 가능하게 하는 이론적 프레임워크를 구축하였다.
동적 손상 진화: 탄성 계수를 상수로 취급하는 기존 모델과 달리, 이 접근법은 하중 과정 중의 등가 탄성 계수의 동적 발달과 거시적 및 미시적 손상의 진화를 모두 반영한다.
정확도 향 향상: 본 연구는 제안된 다중 스케일 계산이 기존의 유한 요소 시뮬레이션 방법보다 실제 실험실 및 현장 테스트 데이터에 더 근접한 결과를 산출함을 입증하였다.

결과

실험실 결과:
- P파 속도( $V_p$ )는 초기 손상이 증가함에 따라 감소하는 것으로 나타났으며, 이는 손상 지표로서의 유용성을 확인시켜 준다.
- CT 스캔 결과, 높은 초기 손상은 높은 공극률 및 균열 밀도와 상관관계가 있음이 밝혀졌다.
- 브라질리언 인장 시험: 하중-변위 곡선은 엄격하게 낮은 볼록 함수(lower convex function) 형태를 보였다. 다중 스케일 계산 결과는 특히 피크 이전 영역(pre-peak region)에서 기존 시뮬레이션보다 실험 곡선과 더 밀접하게 일치하였다.
- 일축 압축 시험: 응력-변형률 곡선은 단조 증가하는 오목 함수 형태를 보였다. 다중 스계 접근법은 등가 탄성 계수를 정확하게 예측하였으며, 그 값은 기존 시뮬레이션보다 테스트 결과에 더 근접하였다. 손상 진화는 초기에 점진적으로 증가하다가 급격히 발달한 후 안정화되는 양상을 보였다.
현장 결과:
- 그라우팅 고결체의 현장 압축 시험에서, 다중 스케일 접근법은 최대 일축 압축 응력까지의 응력-변형률 곡선을 정확하게 예측하였다.
- 그라우팅 고결체에 대해 계산된 등가 탄성 계수는 현장 테스트 데이터와 일치하였으며, 이를 통해 복잡한 공학적 시나리오(파쇄된 광석 및 시멘트 슬러리 골재 포함)에 대한 방법론의 적용 가능성을 검증하였다.
관찰된 한계점: 본 연구는 피크 이전 결과는 정확했으나, 그라우팅 고결체의 피크 이후(post-peak) 응력-변형률 곡선은 현장 결과와 차이를 보였다. 이러한 불일치는 모델에 심부 취성 암석과는 다른 재료의 파괴 후 거동을 포착하는 데 필요한 변형 연화(strain-softening) 매개변수가 결여되어 있기 때문인 것으로 분석되었다.

의의 및 주장
본 논문은 제안된 다중 스케일 계산 접근법이 암반 공학의 복잡한 문제를 해결하기 위한 견고한 이론적 방법을 제공한다고 주장한다. 초기 손상을 고려하고 섭동 이론을 활용함으로써, 이 방법은 미시 구조의 불균질성과 거시적 역학 거동 사이의 간극을 성공적으로 메운다. 저자들은 이 접근법이 암반의 등가 탄성 계수를 결정하는 데 있어 기존의 시뮬레이션 방법보다 더 정확한 대안을 제공한다고 단언한다. 또한, 거시적 척도와 미시적 척도 사이의 적절한 가교를 식별할 수 있다면, 이 방법론을 다른 불균질 재료의 역학적 매개변수를 결정하는 데에도 확장 적용할 수 있을 것이라고 제안한다. 본 연구는 특히 초기 손상과 스케일 효과가 중요한 요소인 공학 프로젝트에서 암반 안정성 평가를 개선하기 위한 참고 자료로 제시된다.

유사한 논문