Mesoscale Modelling of Confined Split-Hopkinson Pressure Bar Tests on Concrete: Effects of Internal Damage and Strain Rates

이 연구는 유한요소법을 기반으로 한 메조스케일 모델링을 통해 콘크리트의 SHPB 실험을 시뮬레이션하여, 하중 램프 속도, 내부 마찰 및 구속 압력이 동적 증가 계수 (DIF) 와 변형률 속도 효과에 미치는 영향 및 그 미시적 손상 메커니즘을 규명하였다.

핵심

🏗️ 연구의 핵심: "콘크리트 속의 미니어처 도시"

콘크리트는 단순히 회색 덩어리가 아닙니다. 연구자들은 이를 세 가지로 이루어진 복잡한 도시로 봅니다.

1. 모르타르 (시멘트 반죽): 도시의 도로와 건물 사이를 채우는 흙.
2. 골재 (자갈): 도시의 거대한 건물들.
3. ITZ (접합부): 건물과 흙이 만나는 경계선 (가장 약한 고리).

이 연구는 이 '도시'를 컴퓨터 안에 3D 로 재현하고, **분리형 호프킨슨 압축 시험기 (SHPB)**라는 장비를 이용해 마치 폭탄을 터뜨리거나 기차로 들이받는 것처럼 아주 빠른 속도로 충격을 가해 보았습니다.

🚀 무엇을 발견했나요? (세 가지 주요 변수)

연구자들은 충격이 가해질 때 세 가지 요소를 바꿔가며 실험했습니다.

1. 충격을 주는 '속도' (로드 램프 속도)

• 비유: 갑자기 차를 멈추는 것 vs 천천히 브레이크를 밟는 것.
• 발견: 충격을 더욱 급격하게 (빠르게) 가할수록 콘크리트는 더 단단해집니다.
• 이유: 충격이 너무 빨라서 콘크리트 속의 '자갈'과 '흙'이 서로 미끄러지거나 부서질 시간이 없기 때문입니다. 마치 급하게 문을 닫으면 문이 더 단단하게 잠기는 것과 비슷합니다.

2. 내부의 '마찰력' (내부 마찰)

• 비유: 미끄러운 얼음 위를 걷는 것 vs 거친 모래 위를 걷는 것.
• 발견: 콘크리트 내부 입자들 사이의 마찰이 클수록 전체적인 강도는 높아지지만, 충격 속도에 따른 강도 증가 효과는 오히려 줄어듭니다.
• 이유: 마찰이 크면 입자들이 서로 꽉 잡혀서 움직이지 못합니다. 하지만 충격이 너무 빨라지면, 이 꽉 잡힌 상태가 오히려 '자갈'이 깨지는 것을 막아주지 못하고, 오히려 '흙 (모르타르)' 부분에서 더 쉽게 무너지게 만들어 속도가 빨라져도 강도가 크게 오르지 않게 합니다.

3. 옆에서 누르는 '압력' (구속 압력)

• 비유: 풍선을 손으로 꽉 쥐고 있는 상태 vs 그냥 놓아둔 상태.
• 발견: 옆에서 누르는 힘이 강할수록 콘크리트는 훨씬 더 강해집니다. 하지만 충격 속도가 빨라져도 강도가 오르는 정도는 느려집니다.
• 이유: 옆에서 누르는 힘 (구속) 이 있으면 콘크리트가 터져나갈 공간이 없습니다. 그래서 자갈이 깨지는 것은 막아주지만, 오히려 '흙 (모르타르)' 부분이 찌그러지면서 무너지는 경향이 생겨, 속도가 빨라져도 강도 증가폭이 작아집니다.

🔍 연구의 핵심 통찰: "속도계와 파손 지도"

이 논문이 가장 혁신적인 점은 왜 그런 현상이 일어나는지 그 '속'을 보여줬다는 것입니다.

• 내부 속도계 (Internal Strain Rate): 충격이 가해졌을 때, 콘크리트 속의 각 부분 (자갈, 흙) 이 얼마나 빠르게 변형되는지 측정했습니다.
• 파손 지도 (Local Damage): 어디가 먼저 깨지고, 어디가 더 많이 부서지는지 시각화했습니다.

결론적으로:
충격이 가해졌을 때, 콘크리트 속의 '흙 (모르타르)' 부분이 얼마나 빠르게 변형되고 부서지느냐가 전체 콘크리트의 강도 변화를 결정하는 열쇠였습니다.

• 충격을 급격히 주면 흙이 빠르게 변형되면서 전체가 더 단단해집니다.
• 하지만 옆에서 누르거나 마찰이 크면, 흙의 변형 속도가 느려져 전체 강도 증가폭이 줄어듭니다.

💡 이 연구가 우리에게 주는 메시지

이 연구는 단순히 "콘크리트가 강하다"는 사실을 넘어, "어떤 조건에서, 왜, 어떻게 부서지는지" 그 미세한 메커니즘을 밝혀냈습니다.

• 실제 적용: 지진, 폭발, 낙석, 선박 충돌 등 예측 불가능한 재해에 대비한 더 안전한 교량, 터널, 방호벽을 설계하는 데 큰 도움이 됩니다.
• 미래: 이제 우리는 콘크리트의 '속'을 컴퓨터로 미리 시뮬레이션하여, 재해가 발생하기 전에 가장 약한 부분을 찾아내고 보강할 수 있게 되었습니다.

한 줄 요약:

"콘크리트는 충격이 가해질 때, 속의 작은 자갈과 흙이 서로 어떻게 부딪히고 부서지느냐에 따라 그 강도가 결정됩니다. 이 연구는 그 미세한 싸움의 과정을 컴퓨터로 재현하여, 더 튼튼한 건축물을 만드는 비법을 찾아냈습니다."

기술 요약

1. 연구 배경 및 문제 제기 (Problem)

• 배경: 콘크리트는 교량, 터널, 해상 풍력 발전기 등 다양한 인프라에 사용되며, 지진, 낙석, 쓰나미, 선박 충돌 등 극한 동적 하중을 견뎌야 합니다. 이러한 동적 하중 하에서 콘크리트의 거동을 이해하는 것은 구조물 설계에 필수적입니다.
• 문제점:
  • 콘크리트는 모르터, 골재, 계면 전이층 (ITZ) 으로 구성된 이질적 재료로, 거시적 실험만으로는 동적 파괴의 미세한 메커니즘을 명확히 규명하기 어렵습니다.
  • 기존 실험 (SHPB 등) 은 하중의 가해지는 속도 (램프 속도), 내부 마찰, 구속 압력 (Confining pressure) 이 동적 강도 증가 계수 (DIF) 에 미치는 영향을 정량화하는 데 한계가 있습니다.
  • 특히, 하중 파형의 모양 (램프 속도) 이나 구속 조건 변화에 따른 미시적 변형률 분포와 국부 손상 진화의 상관관계에 대한 체계적인 연구가 부족합니다.

2. 연구 방법론 (Methodology)

이 연구는 유한 요소법 (FEM) 을 기반으로 한 메조스케일 (Mesoscale) 모델링을 통해 SHPB 시험을 시뮬레이션하였습니다.

• 메조스케일 모델 생성:
  • 콘크리트를 3 상 (모르터, 골재, ITZ) 복합재로 모델링했습니다.
  • 현실적인 골재 형상: 볼로노이 테셀레이션 (Voronoi tessellation) 과 푸리에 기술자 (Fourier descriptor) 기반 알고리즘을 결합하여 구형이 아닌 불규칙하고 현실적인 골재 형상을 생성했습니다.
  • 메쉬: 사면체 요소 (약 1mm) 를 사용하여 모르터와 골재를 모델링하고, 0 두께의 코히시브 인터페이스 요소 (CIE) 를 사용하여 ITZ 를 표현했습니다.
• 구성 모델 (Constitutive Models):
  • 모르터 및 골재: 콘크리트 손상 소성 (CDP) 모델을 사용하며, FIB 모델 코드 2010 에 기반하여 압축 및 인장 강도가 변형률 속도에 따라 증가하도록 정의했습니다.
  • ITZ: 코히시브 인터페이스 요소 (CIE) 를 사용하여 인장 - 분리 법칙 (Traction-separation law) 을 적용했습니다.
  • 동적 하중: SHPB 시험을 모사하기 위해 스타이커 (striker) 가 입사봉을 타격하여 생성된 응력파를 입력 조건으로 사용했습니다.
• 변수 분석:
  • 하중 램프 속도 ($V_L$): 파형의 상승 속도를 변화시켜 분석.
  • 내부 마찰 계수 ($f_i$): 재료 간 접촉면의 마찰 계수를 변화시켜 분석.
  • 구속 압력 ($\sigma_c$): 측면 구속 압력을 적용하여 분석.
• 검증: 생성된 모델의 결과 (DIF 대 변형률 속도 곡선) 를 기존 실험 데이터와 비교하여 모델의 신뢰성을 검증했습니다.

3. 주요 기여 및 결과 (Key Contributions & Results)

A. 하중 램프 속도 (Loading Ramp Rate) 의 영향

• 거시적 결과: 하중 램프 속도가 증가하면 동적 강도 증가 계수 (DIF) 가 증가하며, 특히 고 변형률 속도 영역에서 DIF 증가 폭이 더 커집니다 (DIF-변형률 속도 곡선의 기울기 $\beta$ 증가).
• 미시적 기작:
  • 높은 램프 속도는 시료 내부의 국부 변형률 속도 분포를 더 넓고 높은 값으로 확장시킵니다.
  • 모르터와 골재의 손상 (Damage) 이 급격히 증가하여 전체 강도 상승을 유도합니다.
  • 결론: 높은 램프 속도는 변형률 속도 효과를 증폭시켜 DIF 를 더욱 크게 만듭니다.

B. 내부 마찰 (Internal Friction) 의 영향

• 거시적 결과: 내부 마찰 계수가 증가하면 정적 강도는 증가하지만, 변형률 속도에 따른 DIF 증가 폭 (기울기 $\beta$) 은 감소합니다.
• 미시적 기작:
  • 높은 마찰은 국부 변형률을 억제하여 전체적인 변형률 속도를 낮춥니다.
  • ITZ 와 골재 간의 손상 경쟁 메커니즘이 발생하지만, DIF 기울기 감소의 주된 원인은 모르터 (Mortar) 의 손상 증가가 둔화되기 때문입니다.
  • 결론: 높은 내부 마찰은 변형률 속도 효과를 약화시킵니다.

C. 구속 압력 (Confining Pressure) 의 영향

• 거시적 결과: 구속 압력이 증가하면 콘크리트의 동적 강도 (DIF) 는 전반적으로 증가하지만, 변형률 속도에 따른 DIF 증가율은 감소합니다.
• 미시적 기작:
  • 구속 압력은 골재의 손상을 크게 증가시키지만, 모르터의 손상 증가율은 변형률 속도에 덜 민감해집니다.
  • 측면 구속으로 인해 변형이 제한되어 내부 변형률 속도 분포가 좁아집니다.
  • 결론: 구속 압력은 강도를 높이지만, 변형률 속도 의존성 (Strain-rate sensitivity) 을 약화시킵니다.

D. 미시적 지표의 상관관계

• 연구는 **내부 변형률 속도의 통계적 평균값 ($M_{\dot{\varepsilon}_{eq}}$)**과 정규화된 DIF 사이에 강한 양의 상관관계 ($R=0.94$) 가 있음을 발견했습니다.
• 이는 거시적인 동적 거동을 설명하는 데 있어 미시적인 변형률 이질성 (Heterogeneity) 이 핵심 지표임을 입증했습니다.

4. 연구의 의의 및 결론 (Significance & Conclusion)

• 실험적 한계 극복: 실험적으로 접근하기 어려운 내부 마찰, 구속 압력, 하중 파형의 미세한 변화가 콘크리트 동적 거동에 미치는 영향을 메조스케일 시뮬레이션을 통해 정량적으로 규명했습니다.
• 기작 규명: DIF 의 변화가 단순히 전체 변형률 속도뿐만 아니라, 모르터, 골재, ITZ 의 국부 손상 분포와 내부 변형률 속도의 공간적 분포에 의해 결정됨을 보여주었습니다. 특히 모르터의 손상 거동이 DIF-변형률 속도 관계의 기울기를 결정하는 주요 인자임을 밝혔습니다.
• 실용적 가치: 복잡한 동적 하중 조건 (지진, 충격 등) 하에서 콘크리트 구조물의 안전성을 평가하고, 더 정확한 동적 구성 관계식 (Constitutive relationships) 을 개발하는 데 이론적 기반을 제공합니다.
• 향후 과제: 단일 요인의 영향을 분석했으나, 향후 이러한 요인들 간의 결합 효과 (Coupled effects) 와 지속 가능한 콘크리트 재료 (재활용 골재 등) 에 대한 연구가 필요하다고 제안했습니다.

이 논문은 콘크리트의 동적 파괴 메커니즘을 거시적 현상과 미시적 물리적 과정 (손상 및 변형률) 을 연결하여 심층적으로 이해하는 중요한 통찰을 제공합니다.