Numerical Study of the Mechanical Response and Failure Evolution of Layered Cemented Tailings Backfill under Biaxial Confinement

기술 요약: 이축 구속 조건 하에서의 층상 시멘트 테일링 뒷채움재의 기계적 응답 및 파괴 진화에 관한 수치 연구

문제 정의
광산 운영이 심부로 확장됨에 따라, 뒷채움재는 더 높은 지중 응력과 광산 교란에 직면하게 되며, 이에 따라 견고한 지지력과 장기적 안정성이 요구된다. 단계별 채움은 고함량 시멘트 하부층 뒤에 저함량 시멘트 층을 배치하는 방식이 일반적이며, 이는 상당한 기계적 성질 차이를 가진 층상 구조를 형성한다. 이러한 층 사이의 계면은 종종 결합력이 약하고 공극률이 높아, 균열 개시의 잠재적인 응력 집중 구역으로 작용한다. 이러한 재료의 일축 거동은 연구된 바 있으나, 시편 준비와 저강도 재료의 변형 측정의 어려움으로 인해 삼축 또는 이축 연구는 제한적이다. 또한, 기존의 수치 연구는 종종 확률적 변동성, 입자 충진(packing)의 차이, 그리고 구속 조건 하에서의 특정 층 순서 효과를 고려하지 못하는 경우가 많다. 본 연구는 고해 resolución 수치 접근 방식을 사용하여 이축 구속 하에서의 층상 시멘트 테일링 뒷채움재의 기계적 응답 및 미시적 결합 파괴(bond-break) 진화를 이해할 필요성을 다룬다.

연구 방법론
본 연구는 이산 요소 이론(discrete-element theory)에 기반한 입자 유동 코드(PFC2D)를 채택하였다. 방법론은 세 단계로 구성된다:

1. 기존 데이터에 대한 보정(Calibration): 새로운 일축 압축 실험은 수행되지 않았다. 대신, 저자들은 균질한 시멘트 대 테일링 비율이 1:4 및 1:8인 경우의 황 등(Huang et al., 2022)의 기존 일축 압축 데이터를 활용하였다. 약 38,700개의 입자를 포함하는 고해상도 PFC2D 모델은 특정 무작위 실현값(S10001)을 사용하여 보정되었다. 보정은 1:4 및 1:8 재료의 일축 압축 강도(UCS)와 피크 변형률을 일치시키는 데 중점을 두었다.
2. 독립적 검증: 매개변수가 고정된 후, 모델을 추가 튜닝 없이 검증하기 위해 두 개의 독립적인 무작위 실현값(S20001 및 S30001)을 생성하였다. 이는 모델이 확률적 틀 내에서 실험적 평균값을 재현할 수 있음을 보장하기 위함이다.
3. 이축 구속 시뮬레이션: 보정된 매개변수는 엄격하게 대칭을 이루는 수평 층상 조립체(조합 I: U1:8/L1:4; 조합 II: U1:4/L1:8)에 적용되었다. 시뮬레이션은 목표 구속 압력 0.2, 0.4, 0.6 MPa에서 수행되었다. 연구에서는 변동성을 정량화하기 위해 각 구성마다 세 개의 무작위 시드(random seed)를 사용하였다.
   • 계면 가정: 독립적인 계면 전단 데이터가 없었으므로, 수평 층 경계를 가로지르는 접촉점들은 더 약한 1:8 재료의 특성을 갖도록 할당되었다.
   • 하중: 이축 하중이 가해졌으며 측압 제어가 활성화되었다. 관성 효과를 평가하기 위해 하나의 케이스에 대해 속도 민감도 체크를 수행하였으나, 저자들은 수치적 변형률이 실험적 속도보다 현저히 빠르다는 점을 명시적으로 언급하였다.

주요 기여

• 고해상도 보정: 본 연구는 특정 일축 벤치마크에 따라 보정된 고해상도 PFC2D 모델을 구축하였으며, 독립적인 무작위 실현을 통해 예측 능력을 검증하였다.
• 엄격한 대칭 프로토콜: 동일한 입자 접촉 및 벽 조건을 가진 엄격하게 대칭되는 조립체(CI 및 CII)를 생성함으로써, 이 연구는 확률적 변동성과 경계 아티팩트로부터 층 순서의 효과를 분리하였다.
• 미시적 결합 파괴 분석: 본 연구는 구속 조건 하에서의 결합 파괴 진화(인장 vs 전단)에 대한 상세한 통계적 분석을 제공하며, 특히 심각한 후기 단계 파괴의 복잡성을 피하기 위해 초기 피크 이후 단계에 초점을 맞추었다.
• 변동성 정량화: 세 개의 무작위 시드를 사용함으로써 단일 실현 연구보다 더 견고한 통계적 근거를 제공하는 평균값과 표준 편차를 보고할 수 있었다.

결과

• 모델 정확도: 독립적 검증 실현값(S20001, S30001)은 발표된 실험 평균 대비 -5.14%(1:4) 및 -2.88%(1:8)의 오차를 보이는 평균 UCS 값을 나타내어, 보정된 매개변수에 대한 모델의 타당성을 확인하였다.
• 응력-변형률 응답: 이축 구속 하에서 피크 편차 응력은 목표 구속 압력에 따라 단조롭게 증가하지 않았다. CI 및 CII 구성 모두에서, 피크 응력은 구속 압력이 0.2에서 0.6 MPa로 증가함에 따라 일반적으로 감소하거나 안정된 상태를 유지하였다.
• 층 순서 효과: 층 순서를 역전시키는 것(CI vs CII)은 재현 가능한 강도 이점을 생성하지 않았다. 대칭 모델 조건에서 미러링된 구성 간의 피크 응력 차이는 무작위 시드와 구속 수준에 따라 부호와 크기가 변하였으며, 이는 특정 층 순서에 내재된 기계적 우월성이 없음을 시사한다.
• 결합 파괴 진화:
  • 구속이 증가함에 따라, 일치된 초기 피크 이후 단계(피크 응력의 80%)에서의 평균 총 결합 파괴 횟수가 크게 감소하였다 (예: CI의 경우 ~8,956에서 ~3,762로 감소).
  • 인장 결합 파괴가 파괴 과정을 지배하였으며, 전체 파괴의 97.8%에서 99.2%를 차지하였다. 전단 파괴는 미미했다.
  • 계면으로 라벨링된 파괴는 전체 파괴의 작고 배타적이지 않은 부분집합이었으며, 이는 초기 피크 이후 단계에서 파괴가 계면에만 국한되어 발생하지 않았음을 나타낸다.
• 속도 민감도: 대표적인 체크 결과, 수치적 변형률을 0.50에서 0.10 s⁻¹로 줄였을 때 피크 응력(<0.24%)과 변형률(<0.46%)의 변화가 미미하여, 피크 응답이 테스트된 수치 범위 내에서 관성력에 의해 지배되지 않음을 시사하였다. 그러나 이것이 실험적 하중 속도와의 동등성을 확립한 것은 아니다.

의의 및 주장
본 논문은 연구 결과를 모델의 가정 범위 내에서의 수치적 예측 및 기술적 비교로 겸허하게 규정한다. 저자들은 다음과 같이 명시하였다:

• 본 연구는 고해상도 수치 보정 및 이축 구속 예측을 통해 기존에 발표된 실험적 연구를 확장한다.
• 구속 증가에 따른 결합 파괴 횟수의 감소 관찰은 측면 구속이 인장 결합 열림을 억제하는 것과 일치하지만, 독립적인 구속 실험 데이터가 없기 때문에 이것이 검증된 물리적 메커니즘은 아니다.
• 재현 가능한 강도 이점이 없다는 결과는 층 순서 단독으로 대칭 조건 하에서 용량을 결정한다는 가정에 의문을 제기하며, 상당한 효과를 위해서는 측정되지 않은 비대칭성(예: 밀도 구배 또는 양생 변동)이 필요함을 시사한다.
• 결과는 세 가지 실현으로부터 도출된 기술적 추세이며, 보편적인 구성 법칙이나 일반적인 계면 제어 전단 메커니즘을 확립하는 것으로 해석되어서는 안 된다.

요약하자면, 본 연구는 층상 뒷채움재를 분석하기 위한 엄격한 수치적 프레임워크를 제공하며, 구속이 결합 파괴를 억제할 수는 있지만 반드시 피크 편차 응력을 단조롭게 증가시키는 것은 아니며, 특정 대칭 모델 조건 하에서는 층 순서가 본질적으로 강도를 결정하지 않는다는 것을 보여준다.