A Physically Constrained Energy-Based Excess Pore Pressure Model for Saturated Sands with an Ordinary Differential Equation Interpretation

기술 요약: 포화된 모래를 위한 물리적 제약이 적용된 에너지 기반 과잉 간극수압 모델

문제 정의
토양 액상화는 여전히 중요한 지반 공학적 위험 요소이며, 전통적으로 응력 기반, 변형률 기반, 간극수압 기반 및 에너지 기반 방법들을 통해 평가되어 왔다. 이 중 에너지 기반 방법(EBM)은 누적 소산 에너지와 과잉 간극수압 생성을 연관시킴으로써, 등가 사이클 변환의 불확실성 없이 불규칙한 지진 하중을 수용할 수 있다는 뚜렷한 장점을 제공한다. 가장 영향력 있는 EBM 정식화는 Berrill–Davis(1985)의 거듭제곱 법칙($r_u = \alpha W^\beta$)이다. 그러나 본 연구는 Berrill–Davis 모델의 두 가지 근본적인 한계를 식별하였다:

1. 무제한적 응답: 거듭제곱 법칙은 수학적으로 무제한이며, 과잉 간극수압비($r_u$)가 1을 초과할 수 있게 한다. 이는 $r_u$가 액상화 시 1.0에 점근해야 하는 포화된 모래의 물리적 현상과 맞지 않는 불가능한 결과이다.
2. 매개변수 상호 의존성: 보정 매개변수인 $\alpha$와 $\beta$는 흔히 독립적인 것으로 취급되지만, 경험적 연구에 따르면 이들은 서로 상관관계가 있다. 이러한 의존성은 보정을 복잡하게 만들고, 물리적 해석력을 저하시키며, 외삽을 방해한다.

방법론
저자들은 이러한 한계점들을 해결하기 위해 균질하고 깨끗한 포화 모래를 대상으로 재구성된 모델을 제안하며, 이는 다음과 같은 다단계 유도 과정을 거친다:

1. 데이터 분석 및 분해: Ottawa F65 모래(LEAP-2017)와 Akpınar 모래의 반복 삼축 시험 데이터를 사용하여, 정규화된 간극수압비($r_u$)와 정규화된 소산 에너지($W_N$) 사이의 관계를 로그-로그 스케일에서 분석한다. 저자들은 데이터가 단일 직선을 따르는 것이 아니라 최소 세 개의 선형 구간으로 구성되어 있음을 관찰한다.
2. 다항식 근사: 데이터를 상수항이 없는 6차 다항식으로 근사한다. 이 함수는 주요 성장 항($W_N^b$)과 매개 함수 $g(W_N)$으로 분해된다.
3. 테일러 전개 및 단순화: 매개 함수를 추가로 분석한다. 지수 성분에 대한 테일러 전개를 통해 3차 다항식이 데이터 추세를 적절히 표현함을 밝힌다. 이를 통해 후보 모델 형태인 $r_u = W_N^b (a_0 + a_1 W_N + a_2 W_N^2 + a_3 W_N^3)$를 도출한다.
4. 상미분 방정식(ODE) 해석: 수학적 구조를 인식하여, 이 모델을 비동차 코시-오일러(Cauchy–Euler) 상미분 방정식의 해로 재해석한다. 좌변은 시스템의 내부 응답을 나타내며, 우변은 외부 하중 조건을 나타낸다.
5. 물리적 제약: 비물리적인 변동과 매개변수 간의 상호 의존성을 제거하기 위해, 액상화 지점($W_N = 1, r_u = 1$)에서의 특정 경계 조건을 모델에 적용한다:
   • $W_N = 1$일 때 $r_u = 1$.
   • $W_N = 1$에서 수평 접선($r_u' = 0$)을 가져, 곡선이 한계치에 도달할 때 평탄해지도록 보장한다.
   • 하중 과정 전반에 걸쳐 단조 증가($r_u' > 0$).
   • 액상화 시 비동차 항이 소멸하는 것에서 유도된 계수 간의 특정 관계($a_2 = -3a_3$).

이러한 제약 조건들은 매개 다항식의 계수들을 단일 매개변수인 $b$에 대한 식으로 결정하며, 모델을 $b$에 의해서만 지배되는 형태로 축소시킨다.

주요 기여

• 물리적 제약이 적용된 정식화: 제안된 모델은 $0 \le r_u \le 1$을 본질적으로 만족하도록 설계되어, Berrill–Davis 모델에서 흔히 발생하는 비물리적인 과잉 간극수압 과대 예측을 방지한다.
• 단일 보정 매개변수: 매개 계수들을 성장 지수 $b$와 분석적으로 연결함으로써, 기존 정식화에서 발견된 상호 의존성 문제를 해결하고 여러 매개변수를 독립적으로 보정할 필요를 없앴다.
• ODE 프레임워크: 본 연구는 경험적인 간극수압 관계를 비동차 코시-오일러 상미분 방정식의 해로 해석하는 새로운 관점을 제공한다. 이 프레임워크는 응답을 내부 재료 특성을 지배하는 동차 부분(초기 단계에서 지배적)과 외부 하중 조건을 지배하는 특수 부분으로 구분한다.
• 검증: 모델은 다양한 상대 밀도(느슨함에서 매우 조밀함까지) 범위에 걸친 Ottawa F65 및 Akpınar 모래의 고품질 반복 삼축 시험 데이터에 대해 $R^2 \ge 0.99$를 달 Achieve 하며 검증되었다.

결과

• 단순화된 모델(식 20)은 Ottawa F65 및 Akpınar 모래의 측정된 간극수압 이력을 성공적으로 재현한다.
• 매개변수 $b$는 상대 밀도에 따라 변하는 것으로 나타났다 (예: 상대 밀도 51%일 때 $b \approx 0.36$, 80–91%일 때 $b \approx 0.43$).
• 모델은 액상화에 접근함에 따라 발생하는 초기 급격한 상승과 이후의 평탄화 과정을 포함한 비선형적 간극수압 진화를 포착한다.
• 기존 Berrill–Davis 모델과의 비교를 통해, 제안된 정식화가 특히 완전한 액상화에 도달하는 종단 조건 근처에서 원래 모델이 보이는 큰 산포와 무제한적 성장을 극복하고, 간극수압 진화를 더 일관되게 설명함을 입증하였다.
• 단순화된 모델을 나머지 LEAP-2017 데이터셋에 적용한 결과, 국부 최댓값을 피하기 위한 조건인 $b \ge 0.27$을 만족하는 경우 다양한 상대 밀도와 반복 응력비에 걸쳐 간극수압 생성을 시뮬레이션할 수 있음이 확인되었다.

의의
본 논문은 이 연구가 액상화 연구를 위한 균질하고 깨끗한 모래에 대해 클래식한 Berrill–Davis 거듭제곱 법칙을 대체할 수 있는 실행 가능하고, 간결하며, 물리적으로 일관된 대안을 제공한다고 주장한다. 응답을 제한하고 ODE 기반 유도를 통해 매개변수 간의 상호 의존성을 제거함으로써, 모델은 높은 정확도를 유지하면서도 보정을 단순화한다. 저자들은 이 정식화가 경험적임에도 불구하고, 엄격한 물리적 제약 하의 회귀 분석으로부터 유도되었기에 액상화의 근본 메커니즘을 포착할 수 있다고 강조한다. 즉, 내부 입자 특성에 의해 지배되는 초기 단계(동차해)와 재료가 기계적 기억을 잃고 완전한 액상화에 도달하기 위해 상당한 에너지를 필요로 하는 단계(특수해)를 포착한다. 연구의 범위는 균질하고 깨끗한 포화 모래로 명시적으로 제한된다.