Mg‑Zn‑Ca 합금 정적 재결정 속도 해석: PRISMS 통합 시뮬레이션과 실험 데이터 융합

초록

본 연구는 Mg‑3Zn‑0.1Ca 합금의 정적 재결정 과정을 PRISMS 기반 결정소성‑위상장 모델을 결합해 시뮬레이션하고, 20 % 변형 후 350 °C에서의 재결정 분율 변화를 실험과 비교한다. 시뮬레이션은 평균 입계 이동도와 플라스틱 작업이 저장에너지로 전환되는 비율을 조정함으로써 시간 스케일 팩터를 맞추며, 두 파라미터 중 하나가 알려지면 다른 파라미터를 역산할 수 있음을 보여준다. 저온(≤300 °C)에서는 실험에서 재결정 속도가 급격히 감소하는 현상이 관찰되며, 이는 용질‑드래그, 전단‑결합 이동, 입계 에너지 이방성 등 추가 물리 메커니즘을 모델에 포함시켜야 함을 시사한다.

상세 분석

이 논문은 정적 재결정 현상을 정량적으로 예측하기 위해 두 가지 핵심 모델을 통합한다. 첫 번째는 PRISMS‑Plasticity를 이용한 결정소성 유한요소 해석(CPFE)으로, 변형 과정에서 발생하는 전단·전단계와 트윈을 의사‑슬립 메커니즘으로 구현하고, 플라스틱 작업량 Wₚ를 계산한다. 여기서 플라스틱 작업 중 저장에너지로 전환되는 비율 β는 1 %~15 % 사이의 불확실성을 갖는다. 두 번째는 PRISMS‑PF를 이용한 위상장(Phase‑Field) 모델로, 각 입계는 오더 파라미터 ηᵢ로 표현되며, 저장에너지 차이 ΔEₛ가 구동력으로 작용한다. 입계 이동속도 v는 v = M·ΔEₛ (곡률 구동력 무시) 형태로 기술되며, 여기서 M은 평균 입계 이동도이다.

실험적으로는 Mg‑3Zn‑0.1Ca(ZX30) 합금을 5 %, 10 %, 20 %의 진응 변형 후 350 °C, 300 °C, 250 °C 등 다양한 온도에서 일정 시간 어닐링하고, EBSD‑GOS 맵을 통해 재결정 분율 X(t)를 측정한다. 시뮬레이션 결과는 350 °C, 20 % 변형 조건에서 시간 스케일 팩터 τ = t·(M_ref/M)·(β/β_ref) 로 맞출 경우 실험 데이터와 거의 일치한다. 이는 저장에너지 ΔEₛ가 충분히 큰 경우 입계 이동이 저장에너지 구동에 의해 지배된다는 가정을 검증한다.

하지만 저온에서는 실험에서 재결정이 후기 단계에서 급격히 둔화되는 현상이 관찰되는데, 시뮬레이션은 이러한 둔화를 재현하지 못한다. 저자들은 이를 설명하기 위해 다음과 같은 물리적 메커니즘을 제시한다. (1) 용질‑드래그: Zn·Ca 원소가 입계에 농축되어 이동도를 감소시킴. (2) 전단‑결합 입계 이동(shear‑coupled migration): 입계 이동이 순수한 곡률 구동이 아니라 전단 변형과 결합될 때 속도가 온도에 민감하게 변함. (3) 입계 에너지·이동도 이방성: 특정 미소축 방향·면에 따라 M이 크게 달라져 평균값만으로는 후기 성장 속도를 설명하기 어려움. (4) 재결정 씨앗 생성 메커니즘의 온도 의존성: 저온에서는 새로운 씨앗이 형성되는 속도가 감소해 전체 성장률이 제한됨.

또한, 논문은 β와 M 사이의 상호보완 관계를 정량적으로 제시한다. β가 독립적으로 측정 가능할 경우(예: 인‑시틸레이션 실험) M을 역산할 수 있으며, 반대로 M을 독립적으로 측정(예: X‑ray tomography)하면 β를 추정할 수 있다. 이는 재결정 모델링에서 가장 큰 불확실성인 “플라스틱 작업 → 저장에너지 전환 효율”을 실험적으로 보정할 수 있는 새로운 전략을 제공한다.

마지막으로, 저자들은 향후 모델 개선 방향으로 (i) 원자 수준 시뮬레이션을 통한 입계 이동도·에너지의 미세구조 의존성 파라미터화, (ii) 용질 확산·드래그 모델을 위상장에 결합, (iii) 전단‑결합 이동을 기술하는 추가 항을 도입하여 저온 재결정 둔화를 재현하는 방안을 제시한다. 이러한 확장은 Mg‑Zn‑Ca와 같은 저비용 경량 합금의 공정 설계와 텍스처 제어에 직접적인 영향을 미칠 것으로 기대된다.