기술 요약: 극한 동적 하중에 노출된 철근 콘크리트 구조물의 손상 수준 예측을 위한 머신러닝 알고리즘 개발

1. 문제 정의

테러 공격의 빈도 증가와 잠재적인 산업 사고로 인해 철근 콘크리트(RC) 구조물의 취약성이 높아지고 있다. 폭발 하중에 의한 손상을 평가하는 전통적인 방법인 실험적 테스트와 수치 시뮬레이션은 상당한 한계에 직면해 있다. 실험적 분석은 비용이 많이 들고 시간이 오래 걸리며 심각한 안전 위험을 초래한다. 반대로, 수치 시뮬레이션(예: LS-DYNA를 이용한 유한 요소 해석)은 상세한 통찰력을 제공하지만, 방대한 계산 자원과 전문 지식을 필요로 하므로 신속한 예비 설계나 대규모 취약성 스크리닝에는 부적합하다.

폭발 하중을 받는 RC 기둥의 손상 수준을 예측하기 위한 더 단순하고 빠르며 정확한 방법이 절실히 필요하다. 본 연구는 전통적인 시뮬레이션 방법의 데이터 기반 대안으로서 머신러닝(ML) 알고리즘의 효용성을 조사함으로써 이 간극을 해결하고자 한다.

2. 연구 방법론

2.1 데이터 생성 및 수치 시뮬레이션

본 연구는 물리적 실험 대신 수치 시뮬레이션을 통해 생성된 포괄적인 데이터셋을 활용하였다.

• 소프트웨어: LS-DYNA를 사용하여 폭발 하중 하에서의 RC 기둥의 동적 응답을 모델링하였다.
• 모델 세부 사항: 모델은 특정 종방향 및 횡방향 철근을 가진 직사각형 RC 기둥(단면 500×700 mm, 높이 4400 mm)을 나타낸다.
• 재료 모델:
  • 콘크리트: 비선형 거동, 손상 축적 및 변형률 효과를 포착하기 위해 MAT_CONCRETE_DAMAGE_REL3 (MAT 72Rel3)를 사용하여 모델링하였다.
  • 강철: MAT_PIECEWISE_LINEAR_PLASTICITY (MAT 24)를 사용하여 모델링하였다.
  • 변형률 속도: 폭발 이벤트에서 전형적으로 나타나는 높은 변형률 속도를 고려하기 위해 CEB-FIB 모델 코드에 기반한 동적 증가 계수(DIF)를 콘크리트와 강철 모두에 적용하였다.
• 하중: 폭발 하중은 피크 과압(peak overpressure)과 양의 위상 지속 시간으로 정의된 단순화된 삼각형 압력-시간 이력을 사용하여 기둥 표면에 시뮬레이션되었다.
• 검증: 수치 모델은 Baylot와 Benvis (2007)의 실험 데이터와 비교 검증되었으며, 중간 높이 처짐(mid-height deflection) 및 잔류 변위에서 밀접한 일치성을 보였다.

2.2 손상 평가

손상 지수($D$)는 잔류 축하중 지지 능력을 기준으로 다음과 같이 정의되었다:
$$D = 1 - \frac{P_{residual}}{P_d}$$
여기서 $P_{residual}$은 손상된 기둥의 능력이고 $P_d$는 손상되지 않은 기둥의 능력이다. 손상 수준은 네 가지 범주로 분류되었다: 저(Low) ($0 < D < 0.2$), 중(Medium) ($0.2 < D < 0.5$), 고(High) ($0.5 < D < 0.8$), 붕괴(Collapse) ($0.8 < D < 1$).

2.3 머신러닝 프레임워크

유한 요소 모델에서 도출된 데이터셋을 바탕으로 다섯 가지 서로 다른 ML 알고리즘을 학습 및 테스트하였다:

1. 랜덤 포레스트 (Random Forest, RF)
2. 라이트 그래디언트 부스팅 머신 (Light Gradient Boosting Machine, LightGBM)
3. 카테고리컬 부스팅 (Categorical Boosting, CatBoost)
4. 엑스트림 그래디언트 부스팅 (Extreme Gradient Boosting, XGBoost)
5. 서포트 벡터 머신 (Support Vector Machine, SVM)

입력 특성 (10개 파라미터):

• 기하학적 특성: 기둥 깊이($D$), 폭($w$), 높이($H$).
• 재료 특성: 콘크리트 압축 강도($f_c$), 종방향 철근 항복 강도($f_y$), 횡방향 철근 항복 강도($f_{yt}$), 종방향 철근비($\rho$), 횡방향 철근비($\rho_s$).
• 폭발 특성: 폭약 중량($W$), 스케일 거리($Z$).

3. 주요 기여

• 비교 분석: 본 연구는 폭발 하중에 의한 RC 기둥의 손상 예측이라는 비선형적이고 임계값이 뚜렷한 문제를 해결하기 위해 다섯 가지 최첨단 ML 알고리즘을 엄격하게 비교하였다.
• 최적 알고리즘 식별: 이 특정 응용 분야에서 RF, XGBoost, LightGBM 및 SVM보다 정확도와 신뢰성 면에서 뛰어난 성능을 보이는 CatBoost를 최적의 모델로 식으로 식별하였다.
• 특성 중요도 정량화: 기하학적, 재료적, 폭발적 파라미터가 손상 심각도에 미치는 상대적 영향력을 정량화하여 설계 최적화를 위한 데이터 기반 통찰력을 제공하였다.
• 민감도 분석: 입력 변수(예: 폭약 중량, 기둥 깊이)와 결과적인 손상 지수 사이의 비선형 관계를 매핑하여 구조적 거동의 임계점을 식별하였다.

4. 결과

4.1 모델 성능

성능은 $R^2$, 평균 제곱 오차(MSE), 평균 절대 오차(MAE), 설명 분산 점수(EVS), 평균 절대 백분율 오차(MAPE), 평균 제곱 로그 오차(MSLE)를 사용하여 평가되었다.

• CatBoost (Tier 1): $R^2$ 0.7125, MAE 1,061.09, MAPE **17.63%**를 달성하며 가장 높은 성능을 보였다. 이는 정확도와 신뢰성의 최적의 균형을 입증했다.
• Random Forest (Tier 2): 좋은 성능($R^2 = 0.6274$)을 보였으며, 가장 높은 안정성(가장 낮은 $R^2$ 표준 편차)을 보여 견고한 대안임을 나타냈다.
• XGBoost (Tier 2): 경쟁력 있는 평균 정확도($R^2 = 0.6109$)를 달성했으나, 높은 변동성(표준 편차 0.2313)을 보여 데이터 분할에 민감함을 시사했다.
• LightGBM (Tier 3): 잎 단위 성장(leaf-wise growth) 전략으로 인해 국부적 패턴에 과적합되어 한계적인 성능($R^2 = 0.5043$)을 보였다.
• SVM (Tier 4): $R^2$가 -0.0170으로 나타나 완전히 실패하였으며, 이는 폭발 하중에 내재된 복잡한 비선형성과 불연속적 전이(예: 국부적 손상에서 붕괴로의 전환)를 포착하지 못함을 의미한다.

4.2 특성 중요도

CatBoost 모델을 기반으로 입력 파라미터를 영향력 순으로 나열하였다:

1. 구조 요소의 깊이 ($D$): 38.45% (지배적 요인).
2. 폭약 중량 (Charge Weight): 27.25%.
3. 종방향 철근의 항복 강도 ($f_y$): 14.55%.
4. 스케일 거리 ($Z$): 6.42%.
5. 횡방향 철근비 ($\rho_s$): 3.04%.
6. 횡방향 철근의 항복 강도 ($f_{yt}$): 2.71%.
7. 콘크리트 압축 강도 ($f_c$): 2.66%.
8. 기둥 높이 ($H$): 2.21%.
9. 기둥 폭 ($w$): 1.78%.
10. 종방향 철근비 ($\rho$): 0.93%.

4.3 민감도 및 임계값

• 깊이: 550 mm와 650 mm 사이에 뚜렷한 임계점이 존재한다. 600 mm를 초과하여 깊이를 늘리는 것은 손상 감소 측면에서 수익 체감의 법칙이 적용된다.
• 폭약 중량: 손상은 폭약 중량에 따라 로그 함수적으로 증가하며, 구조물이 완전 붕괴에 가까워질수록 정체된다.
• 콘크리트 강도: 변형률 효과(DIF)가 정적 강도 차이를 정상화하기 때문에, 45 MPa를 초과하여 $f_c$를 높이는 것은 추가적인 이득이 거의 없다.
• 횡방향 철근: 연구 범위 내에서, 최소 설계 기준을 충족한다는 전제하에 횡방향 철근의 변화는 잔류 축하중 용량 손상 지표에 미치는 영향이 미미했다.

5. 의의 및 주장

본 논문은 개발된 ML 모델, 특히 CatBoost가 폭발 하중 하에서의 RC 기둥의 구조적 손상 및 파괴 모드를 예측하는 신속하고 정확한 도구를 제공한다고 주장한다.

• 실용적 유용성: 적절한 안전 계수가 적용된다는 전제하에, 이 모델들은 예비 설계 및 취약성 스크리닝을 위해 계산 비용이 많이 드는 유한 요소 시뮬레이션을 대체할 수 있다.
• 설계 지침: 본 연구는 데이터에 기반한 구체적인 설계 권고 사항을 제공한다:
  • 폭파 저항성을 위해 콘크리트 강도를 높이거나 횡방향 철근을 추가하는 것보다 기둥 깊이를 늘리는 것을 우선시하라.
  • 변형률 효과를 활용하기 위해 고강도 종방향 철근 ($f_y \ge 500$ MPa)을 사용하라.
  • 폭파 저항성을 위해 콘크리트 강도를 45 MPa 이상으로 높이는 것은 한계 이익이 낮으므로 피하라.
  • 부지 계획 시 스케일 거리(이격 거리)를 주요 보호 조치로 고려하라.

저자들은 개발된 모델이 사용된 파라미터 범위와 수치 시뮬레이션 데이터에 국한된다는 점을 명시하며, 이 모델들이 효율적인 데이터 기반 폭파 저항 설계로 나아가는 중요한 단계라고 결론짓는다. 또한, 향후 연구에는 실험 데이터의 포함과 다른 구조 요소로의 확장이 필요함을 제언한다.