열충격에 취약한 취성재료를 위한 완전 위상장 파괴 모델

초록

본 논문은 열충격으로 인한 온도 구배와 그에 따른 기계적 응력이 취성재료에 균열을 일으키는 현상을, 재료 강도·탄성·파괴인성도를 독립적으로 지정할 수 있는 완전 위상장 파괴 모델을 통해 정량적으로 예측한다. 모델은 Drucker‑Prager 강도면과 공간적 강도 변동을 포함하고, 열‑기계 연성 방정식을 순차적 스키마로 풀어 실험적 파열 패턴(유리판 퀜칭, 적외선 가열 세라믹 디스크, 급속 전력 펄스 핵연료 펠릿)을 성공적으로 재현한다.

상세 분석

이 연구는 기존 위상장 파괴 모델이 갖는 ‘그리피스 전역 에너지 기준’만으로는 균열 핵생을 설명하기 어려운 점을 극복하기 위해, 재료의 탄성계수, 파괴인성도(Gc), 그리고 강도(압축·인장 강도)를 완전히 독립적인 매개변수로 취급하는 ‘완전 모델(complete model)’을 제안한다. 강도면으로는 Drucker‑Prager식을 채택하고, 이를 위상장 방정식에 외부 미소힘 ce 형태로 삽입함으로써, 응력 상태가 강도면을 초과할 때 즉시 v̇<0 조건을 만족하도록 하였다. 또한, 실제 재료가 보이는 강도 변동성을 반영하기 위해 공간적 상관성을 갖는 무작위 강도장(translation model 혹은 random mosaic)을 도입, 이는 균열 경로와 분기 현상에 중요한 역할을 한다.

열‑기계 결합은 온도장 T와 변위장 u를 각각 열전도 방정식과 선형 탄성 방정식으로 기술하고, 위상장 v는 열에 직접적으로는 영향을 받지 않지만, 온도에 의한 열팽창 εT가 기계적 변형 E_m에 기여함으로써 간접적으로 파괴 구동력을 변화시킨다. 수치해석에서는 4‑노드 2D 사각요소와 8‑노드 3D 육면체 요소를 사용하고, 열·기계·위상장을 순차적으로 해결하는 staggered scheme을 적용하였다. 위상장의 비가역성(v̇≥0일 때 v̇=0) 제약은 Primal‑Dual Active Set 방법으로 구현했으며, 격자 의존성을 보정하기 위해 δ에 대한 메쉬 보정식(식 28)을 도입해 실제 파괴인성도와 일치하도록 조정하였다.

세 가지 실험 사례를 통해 모델의 범용성을 검증한다. 첫 번째는 대형 사전 균열이 존재하는 유리판의 급냉(quenching) 실험으로, 온도 구배가 큰 경우에도 강도 파라미터가 균열 전파 속도와 경로에 미치는 영향을 정량화하였다. 두 번째는 적외선 가열에 의해 열응력이 거의 균일하게 발생하는 세라믹 디스크 실험으로, 노치가 있는 경우 직선 균열이, 무노치인 경우 분기가 발생함을 강도장의 공간적 변동성과 Drucker‑Prager 면의 형태가 정확히 재현함을 보여준다. 마지막으로 핵연료 펠릿에 급속 전력 펄스를 가했을 때, 낮은 펄스에서는 전혀 파손되지 않지만 높은 펄스에서는 파편이 발생하는 현상을, 강도와 파괴인성도의 상호작용으로 설명한다. 특히, 파괴가 전적으로 그리피스 에너지에 의해 지배되는 것이 아니라, 강도 한계가 먼저 도달하는 경우가 있음을 강조한다.

전반적으로, 이 완전 위상장 모델은 기존 변분 위상장 모델이 갖는 ‘핵심 파괴 메커니즘(그리피스)’과 ‘강도 기반 핵생(압축·인장)’을 동시에 포괄함으로써, 열충격과 같은 급격한 비정상 환경에서도 취성재료의 파괴를 예측하는 데 있어 가장 포괄적이고 실용적인 수단을 제공한다.