“Yttrium이 Mg‑Y 합금의 주된 텐션 트윈(TT1)·보조 트윈(TT2) 활성화와 국부 변형장에 미치는 영향 해부”

📝 Abstract

Twinning is a primary deformation mechanism in Mg alloys. This study focuses on tension twins during uniaxial compression of Mg-Y alloys, with three key aspects: the orientation specificity of twin grains, the relative evolution of CRSS with increasing Y content, and the local stress and strain evolution at twin sites. Experimental characterization and crystal plasticity modeling were performed. In Mg-7wt.%Y, TT2-{112-1} tension twins were observed in addition to the common TT1-{101-2} twins. Increasing Y suppressed TT1 formation while promoting TT2 activity. A previously unreported group of crystallographic orientations with a higher global Schmid factor for <c+a> slip was identified, which exhibited TT1 twinning with increasing compression strain. To elucidate Y effects on twin activity and local mechanical fields, both TT1 and TT2 tension twin modes were incorporated into PRISMS-Plasticity, an open-source, finite element-based crystal plasticity solver. Four binary Mg-Y alloys were modeled under compression, and statistical analysis was conducted to correlate initial orientations, stress-strain distributions, and twin activities as functions of Y concentration. The plasticity analysis revealed that increasing Y decreases the CRSS ratio of prismatic and pyramidal slip relative to TT1 twinning, while the slip-to-twin CRSS ratio for TT2 increases, thereby serving as a potential indicator of differential twin activity with Y addition in Mg alloys. Additionally, despite their small volume fraction, TT2 twin sites were predicted higher local strain accumulation locally, relative to the representative volume element and TT1 twins, suggesting their potential influence on localized phenomena such as recrystallization or twin nucleation. These findings provide insight into local mechanical behavior in Mg alloys and support alloy design for advanced engineering applications.

💡 Analysis

**

1. 연구 배경 및 의의

• Mg‑Y 합금의 산업적 가치 : 저밀도·고비중강도·우수한 탄성계수를 갖는 Mg 합금은 자동차·항공·에너지 분야에서 경량화와 친환경화에 핵심 소재로 부상하고 있다. Y와 같은 희토류(RE) 원소는 기본적인 베이식 텍스처를 약화시켜 가공성·연성(ductility)을 크게 향상시킨다.
• 트윈(Twinning)의 역할 : HCP 구조의 Mg에서는 c‑축 방향 변형을 수용하기 위해 텐션 트윈이 주요 메커니즘이다. 기존 연구는 주로 TT1 {101‑2} 트윈에 집중했으며, TT2 {112‑1} 트윈은 Mg‑RE 합금에서 드물게 보고된 바 있다.

2. 실험적 접근

3. 결정학적 플라스틱성 모델링 (CPFE)

• 프레임워크 : PRISMS‑Plasticity (오픈‑소스, 유한 요소 기반)
• 구현 내용
  • Slip 시스템: basal , prismatic , pyramidal , pyramidal <c + a> (총 4종)
  • Twin 시스템: TT1 {101‑2}와 TT2 {112‑1}을 pseudo‑slip으로 취급
  • CRSS 및 경화 파라미터는 실험 RVE(대표 부피 요소)와 매크로 흐름곡선, 트윈 부피비를 동시에 맞추는 방식으로 보정
• 핵심 결과
  1. CRSS 비 변화
     • Y 증가 → 프리즘·피라미드 전위 CRSS/TT1 CRSS 비 감소 (TT1 트윈이 상대적으로 쉬워짐)
     • 반대로 TT2 트윈에 대한 전위‑트윈 CRSS 비는 증가 (TT2 활성화가 전위에 비해 상대적으로 어려워짐)
  2. 국부 변형장
     • TT2 트윈 부위는 전체 부피 대비 작은 비율이지만, 동일한 전단 변형을 적은 부피로도 달성 → 해당 영역에서 국부 변형 집중(ε_local ≫ ε_RVE)
     • 이는 재결정화 핵, 추가 트윈 핵생성, 혹은 파단 전조 현상에 중요한 역할을 할 가능성 제시

4. 과학적·공학적 시사점

1. 합금 설계 가이드라인

   • Y 함량을 조절함으로써 TT1·TT2 트윈 비율을 의도적으로 제어 가능 → 원하는 연성·강도 균형을 설계할 수 있다.
   • 특히, TT2 트윈이 제공하는 높은 전단 변형(γ≈0.616) 특성을 활용하면 국부 변형 집중을 이용한 미세구조 조정(재결정화·세분화) 전략을 구상할 수 있다.

2. 모델링 방법론의 확장성

   • PRISMS‑Plasticity에 TT2 트윈을 포함한 접근은 기존 CPFE가 다루지 못했던 다중 트윈 메커니즘을 정량적으로 예측하도록 확장한다.
   • 향후 다른 RE‑Mg 합금(예: Gd, Nd)이나 고온·고속 변형 상황에도 동일 프레임워크 적용이 가능하다.

3. 제한점 및 향후 연구 과제

   • 실험‑모델 일치성: CRSS 보정이 평균적인 메소스케일을 기반으로 하므로, 실제 시편 경계조건(마찰, 온도 구배 등)과는 차이가 존재한다.
   • TT2 트윈의 동역학: 현재는 pseudo‑slip 형태로 구현했으나, 실제 트윈 전파 속도·핵생성 메커니즘을 포함한 동적 모델링이 필요하다.
   • 다중 스케일 연계: 원자 수준의 GPFE 계산과 연계해 Y 원소가 트윈 경계 에너지에 미치는 정량적 영향을 파악하면, CRSS 비 변화의 물리적 근거를 더욱 명확히 할 수 있다.
   • 실용적 적용: 압출·압연 공정에서 Y 함량에 따른 텍스처 제어와 TT2 트윈 활용을 실제 부품 제조에 적용하기 위한 공정 파라미터 최적화 연구가 필요하다.

5. 결론 요약

• Y 함량이 증가함에 따라 Mg‑Y 합금에서 TT1 트윈 억제·TT2 트윈 촉진이라는 뚜렷한 전이 현상이 발생한다.
• CRSS 비의 변화는 TT1·TT2 트윈 활성화 차이를 설명하는 핵심 매개변수이며, TT2 트윈은 작은 부피에도 불구하고 국부 변형 집중을 유발한다.
• PRISMS‑Plasticity 기반의 다중 트윈 CPFE 모델은 실험 결과와 정량적으로 일치하며, 향후 Mg‑RE 합금 설계와 고성능 구조재 개발에 유용한 예측 도구가 될 것이다.

**

📄 Content

마그네슘 합금은 낮은 밀도, 높은 비강도, 높은 탄성계수와 같은 특성을 가지고 있어 운송, 항공우주, 에너지 분야에서 경제적이며 환경 친화적인 대안으로 각광받고 있다[1,2]. 이러한 강도를 더욱 향상시키기 위해 희토류 원소(RE)를 첨가한 Mg 합금이 개발되었다[3,4]. Y(이트륨)와 같은 RE 첨가원소는 제조 과정에서 발생하는 강한 basal 조직의 형성을 약화시켜 Mg의 성형성을 개선한다[5,6]. 또한 Y 첨가는 합금의 연성[7,8], 크리프 저항성[9,10], 부식 저항성[11]을 향상시키는 것으로 보고되었다.

Y를 Mg에 첨가했을 때 가장 두드러지는 효과는 피라미드 <c+a> 전위와 같은 비베이식 전위 활성화이지만[12‑14], Y가 쌍정(twinning) 거동에 미치는 영향 역시 큰 관심을 받고 있다.

쌍정 형성은 c축을 따라 변형을 수용함으로써 Mg의 기계적 거동에 크게 영향을 미친다. 그림 1은 두 종류의 장축 쌍정, TT1 {101̅2}<101̅1>와 TT2 {112̅1}<112̅1>의 개략도를 보여준다. 두 쌍정 모두 변형 중 축 방향의 인장 변형을 수용한다. HCP 금속에서 TT2 쌍정은 전단 방향으로 원자 재배열이 추가로 필요하지 않고 모든 격자점이 올바른 쌍정 위치로 이동할 수 있는 유일한 쌍정 모드이다[15]. TT1 쌍정의 경우, 미소회전 각도 분포는 <112̅0> 회전축을 중심으로 84.8°(±5°)의 특성 피크를 보인다[16]. TT2 {112̅1} 쌍정은 <101̅0> 회전축을 중심으로 35.1°(±5°)의 피크를 나타낸다[16]. TT1의 쌍정 전단 변형은 γ₁ = (a/c)·√(3/2) 로 표현되며, 여기서 R = a/c이다. TT2의 전단 변형은 γ₂ = (c/a)·√(3/2) 로 주어진다[17]. Mg의 c/a 비율이 1.624이므로, TT1의 특성 전단 변형은 약 0.129, TT2는 약 0.616이다. 따라서 TT2 쌍정이 수용하는 전단 변형은 TT1보다 거의 5배에 달하며, 이는 변형 중 변형 수용 및 쌍정 진화에 큰 영향을 미칠 수 있다.

TT1 쌍정은 Mg 합금에서 흔히 관찰되며, 그 핵생, 전파, 그리고 쌍정 간 상호작용을 이해하기 위한 광범위한 연구가 수행되어 왔다.

반면 TT2 쌍정은 Co[18], Ti 합금[19], Zr 합금[20,21] 등 다른 HCP 금속에서는 보고되었지만 Mg 합금에서는 상대적으로 드물다. 그럼에도 불구하고 WE54[22], WE43[23], Mg‑Gd‑Y‑Zr[24], Mg‑17 wt.% Gd[25]와 같은 다양한 Mg‑RE 합금에서 TT2 쌍정이 보고된 바 있다.

특히 Y를 첨가한 경우, Mg‑9 wt.% Y[26]와 Mg‑10 wt.% Y[27]에서 TT2 형성이 관찰되었다.

최근 TT2 쌍정 형성에 대한 연구 관심이 높아지면서, 최첨단 실험 기법과 모델링 방법이 다수 적용되고 있다. 예를 들어, TEM/STEM 분석을 통해 Chen 등은 WE43에서 TT2 쌍정 변이들의 쌍정‑쌍정 상호작용을 조사했으며[23], Zhang 등은 Mg‑Gd‑Y‑Zr 합금에서 높은 변형률에서의 TT2 활성화를 연구하였다[24]. Li 등은 용액 처리된 Mg‑Gd‑Y 합금에서 TT2 쌍정이 TT1보다 먼저 활성화되는 현상을 관찰했지만, 추가 노화 처리된 합금에서는 순서가 뒤바뀌었다[28]. Gengor 등은 ab‑initio 시뮬레이션을 이용해 HCP 재료의 TT2 변이들에 대한 일반화된 평면 결함 에너지(GPFE)를 계산하고, TT2 쌍정에 대해 전단면에 수직인 비영(非零) 셔플 왜곡이 존재함을 강조하였다[29]. 또한 탄성 강성 계수 C₄₄와 쌍정계면 에너지 사이에 비선형 양의 상관관계가 있음을 보고하였다. 이처럼 최근 연구는 주로 열처리 혹은 변형 조건에 따른 쌍정 에너지 혹은 활성화 메커니즘을 이해하는 데 초점을 맞추고 있다. 그러나 Mg 합금에서 TT2 쌍정 활동과 그 변형 중 영향을 정량적으로 분석한 연구는 아직 부족한 실정이다. 따라서 본 연구에서는 상세한 실험적 특성화와 결정 플라스틱성 분석을 결합하여 TT2 쌍정 활동과 미세기계장 진화 사이의 상관관계를 규명하고자 한다. 현재까지 Mg‑Y 합금에 결정 플라스틱성(CP) 모델을 적용한 연구는 드물다. 예를 들어, Li 등은 Mg‑0.8 wt.% Y에 대해 CPFE 분석을 수행해 높은 변형률에서 비베이식 전위 활동이 증가하고 쌍정 활동은 감소한다는 결과를 제시하였다[30]. Huang 등은 Mg‑2 wt.% Y 합금에서 비정상적인 TT1 쌍정 활동을 CPFE로 해석했으며[31], Su 등은 Mg‑2 wt.% Y 합금에 대해 나노인덴테이션과 CPFE를 결합해 변형 중 결정방향 의존성을 분석하였다[32]. Kula 등은 Y 함량이 증가함에 따라 활성 변형 모드들의 임계 전단응력(CRSS)이 단조롭게 증가한다는 것을 보여주었다[33]. 이들 연구는 모두 Y 함량이 비교적 낮은 합금을 대상으로 했으며, 따라서 TT2 쌍정 모드는 CPFE 시뮬레이션에 포함되지 않았다. 최근 Siska 등은 타원형 쌍정 형상을 가진 단결정 내 응력 분포에 대한 통찰을 제공했지만[34], 다결정 변형은 개별 입자의 배향 및 주변 입자와의 상호작용에 의해 추가적인 영향을 받을 수 있다. 따라서 본 연구에서는 Mg‑Y 합금에서 TT1·TT2 두 종류의 쌍정이 동시에 활성화되는 다결정 변형을 체계적으로 통계 분석하고, 결정 플라스틱성 모델과 실험적 쌍정 특성화를 결합한다.

본 연구에 사용된 두 종류의 Mg‑Y 합금은 Y 함량이 각각 1 wt.%(0.28 at.%)와 7 wt.%(2.02 at.%)인 Mg‑1Y와 Mg‑7Y이다. 이 재료들은 캐나다 CanmetMATERIALS에서 압출된 막대로 제공되었으며, 합금 조성 및 압출 공정은 이전 논문[38]에 상세히 기술되어 있다. 직경 15 mm의 압출 막대를 절단한 뒤, 400 °C에서 24 h 열처리하고 수냉하였다.

열처리 후, 직경 10 mm·높이 12 mm인 원통형 시편을 전기 방전 가공(EDM)으로 가공하였다. 원통형 시편의 축방향은 압출 방향과 일치하며, 이는 압축 시험의 축과도 동일하다. 압축 시험은 100 kN 로드셀을 장착한 Instron 장비를 사용해 변위 속도 0.03 mm/min으로 수행하였다. 각 합금에 대해 3 %, 5 %, 8 %의 총 변형률을 적용하였다. 시편 양 끝은 800 SiC 사포로 손으로 연마하고, 마찰을 감소시키기 위해 두 면 사이에 테플론 테이프를 부착하였다. 여기서 변형률은 시편 전체의 평균 변형률을 의미한다.

압축 후, 시편을 하중축에 평행하게 절단하여 미세구조를 관찰하였다. 금속 시편은 Mg 합금용 표준 연마·폴리싱 절차에 따라 연마했으며, 최종 단계에서는 1 µm 다이아몬드 페이스트를 사용하였다. 연마 윤활제로는 물 대신 Buehler MetaDi 액체를 사용하였다. 에칭은 에탄올 60 mL, 물 20 mL, 글레이셜 아세트산 15 mL, 질산 5 mL를 혼합한 용액으로 약 15 s 동안 진행하였다.

결정학적 방향성을 확인하기 위해 전자후방산란(EBSD) 스캔을 수행하였다. Tescan Mira 3 전자현미경에 EDAX Velocity 카메라를 장착하고, 30 kV, 작업거리 약 20 mm, 스텝 크기 500 nm로 667 µm × 500 µm 영역을 스캔하였다. 각 압축 조건마다 시료 중심부에서 4개의 EBSD 스캔을 획득해 통계적 비교를 가능하게 하였다. 수집된 EBSD 원시 데이터는 OIM 소프트웨어 또는 오픈소스 MATLAB 툴박스인 MTEX[40]를 이용해 분석하였다.

Mg‑1Y와 Mg‑7Y의 초기 입자 구조와 조직은 각각 그림 2(a), (b)와 같다. 두 합금 모두 400 °C, 24 h 열처리 후 완전 재결정된 구조를 보였다. 입자 크기는 Mg‑1Y가 51 ± 18 µm, Mg‑7Y가 32 ± 11 µm로 약간 차이를 보였다. 두 합금 모두 압출 막대의 미세구조에서 유래한 약한 basal 섬유 조직을 유지했으며, pole figure 분석에 따르면 <0001>‑c축이 방사 방향과 일치하고 압출 방향에 수직인 경향을 보였다.

결정 플라스틱성 유한요소(CPFE) 시뮬레이션은 오픈소스 PRISMS‑Plasticity 프레임워크[36,41]를 이용해 수행하였다. PRISMS‑Plasticity는 피로 시뮬레이션[42,43], 조직 진화[44,45], 입자 크기가 흐름 거동에 미치는 영향[46] 등 다양한 분야에서 검증된 도구이다. 유한 변형 프레임워크와 구현 상세는 [36]을 참고한다.

본 연구에서는 변형 시 속도 의존성을 무시한 비속도 의존성 구성을 사용하였다. 아래에서는 해당 구성을 간략히 요약한다.

변형 구배 F는 탄성 구배 Fᵉ와 소성 구배 Fᵖ의 곱으로 표현된다.

F = Fᵉ Fᵖ

변형된 현재 좌표계에서의 속도 구배 L은

L = Ḟ F⁻¹

위 식을 이용하면

L = Ḟᵉ Fᵉ⁻¹ + Fᵉ Ḟᵖ Fᵖ⁻¹ Fᵉ⁻¹

즉, 전체 속도 구배 텐서는 탄성 부분 Lᵉ와 소성 부분 Lᵖ로 가법적으로 분해된다.

L = Lᵉ + Lᵖ

여기서 Lᵉ = Ḟᵉ Fᵉ⁻¹, Lᵖ = Fᵉ Ḟᵖ Fᵖ⁻¹ Fᵉ⁻¹이며, Lᵖ는 중간 좌표계에서의 소성 속도 구배 Lᵖ = Ḟᵖ Fᵖ⁻¹ 로 정의된다. Lᵖ는 다수 전위계에 대한 전단률의 중첩으로 표현된다. 본 연구에서는 쌍정을 의사 전위계(pseudo‑slip system)로 간주하여, 전위와 쌍정에 의한 전체 전단률을 Lᵖ에 포함시켰다.

Lᵖ = ∑{α=1}^{N_s} γ̇_α S_α + ∑{k=1}^{N_t} γ̇_k S_k

여기서 γ̇_α는 전위계 α에 대한 전단률, S_α는 중간 좌표계에서의 슈미트 텐서이며, N_s는 전위계 수, N_t는 쌍정계 수이다. m_α와 n_α는 각각 전위계 α의 전단 방향과 전위면 법선 단위벡터이며, 쌍정계 k에 대해서는 m_k와 n_k가 각각 쌍정 전단 방향과 쌍정면 법선을 나타낸다.

비속도 의존성 모델에서는 전위계 α에 대한 전단 응력 τ_α가 임계 전단 응력 s_α에 도달하면 전위가 발생한다. 즉, 전위 조건은

τ_α ≥ s_α

전위가 활성화되지 않을 경우 τ_α < s_α 혹은 |τ_α| = s_α이며, 전단 누적은 위 조건을 만족할 때만 허용된다.

전

이 글은 AI가 자동 번역 및 요약한 내용입니다.