고에너지 폭발에 의한 원통 파편화 시뮬레이션 물질점법의 가능성

초록

본 논문은 고에너지 물질이 내부에 존재하는 원통형 구조가 폭발 가스로 팽창하면서 파편화되는 현상을 물질점법(MPM)으로 시뮬레이션한다. 금속의 소성 변형은 가법적 변형률 텐서 분해와 방사형 복귀 알고리즘을 이용한 저탄성‑소성 응력 업데이트로 모델링하고, 다양한 온도·속도 의존 흐름법칙과 항복조건을 검토한다. 파손은 공극률·손상·분기 조건으로 판단하며, 고속 충격 실험 데이터와 비교해 모델을 검증한다. 결과는 MPM이 대변형·고속 유동‑구조 상호작용 문제에 유망함을 보여준다.

상세 분석

본 연구는 전통적인 유한요소법(FEM)이 대변형, 격자 얽힘, 접촉 처리 복잡성 및 유동‑구조 상호작용(FSI) 해석에서 겪는 한계를 극복하기 위해 물질점법(MPM)을 선택한 점이 핵심이다. MPM은 라그라지안 입자를 에이전트로 두고, 이들을 고정된 오일러 격자와 교차시켜 물리량을 전송함으로써 격자 변형 없이 대변형을 자연스럽게 처리한다. 논문에서는 금속의 소성 거동을 가법적 변형률 텐서 분해(additive decomposition)와 방사형 복귀(radial return) 기법을 결합한 저탄성‑소성(hypoelastic‑plastic) 응력 업데이트 모델로 구현하였다. 이 접근법은 변형률 텐서의 대칭성을 유지하면서도 소성 흐름법칙을 명시적으로 적용할 수 있다.

다양한 플라스틱 흐름법칙—예를 들어 온도·속도 의존 Johnson‑Cook 모델, Zerilli‑Armstrong 모델, 그리고 고온 고속 변형에 적합한 Steinberg‑Guinan 모델—을 시험하여 각 모델이 파편화 과정에서 에너지 흡수와 온도 상승에 미치는 영향을 비교하였다. 항복조건 역시 von Mises와 Drucker‑Prager를 포함해 다중 항복면을 고려했으며, 이는 고압 가스와 금속 사이의 복합 응력 상태를 보다 정확히 포착한다.

파손 판단 기준으로는 입자 수준에서의 공극률(porosity) 증가, 손상 변수(damage) 누적, 그리고 국부적인 응력 텐서의 분기(bifurcation) 조건을 도입하였다. 특히 공극률 기반 파손 모델은 가스 팽창에 의해 발생하는 급격한 압축 해제와 연계된 미세균열 성장 메커니즘을 재현한다.

수치 검증은 고속 충격 실험 데이터와 비교했으며, 충격 파라미터(속도, 압력)와 파편화 패턴에서 MPM 시뮬레이션이 실험과 높은 일치를 보였다. 이는 격자 독립적인 입자 기반 접근법이 복잡한 충격‑파편화 현상을 정량적으로 예측할 수 있음을 의미한다. 최종적으로 연구진은 MPM이 대변형, 고속 유동‑구조 상호작용, 그리고 다중 물질 접촉을 포함한 문제에 대해 높은 정확도와 계산 효율성을 제공한다는 결론을 내렸다.