Finite compressibility and strain hardening in elasto-plastic models of amorphous matter

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

비정질 고체 (금속 유리, 콜로이드 유리, 에멀전 등) 의 소성 변형 및 항복 거동은 국부적인 전단 변환 (Shear Transformations) 에 의해 제어되는 것으로 알려져 있습니다. 전단 변환이 발생하면 Eshelby 의 포함체 문제 (Inclusion problem) 에 따라 주변 물질에 응력 변화가 발생합니다.

• 핵심 문제: 기존 연구들은 대부분 비압축성 (incompressible) 인계를 가정하거나, Poisson 비 (또는 압축률 $K$와 전단률 $\mu$의 비율 $K/\mu$) 에 따른 유한 압축성 효과를 정량적으로 고려하지 않았습니다.
• 관찰된 현상: 실험적으로 금속 유리의 연성 (ductility) 은 $K/\mu$ 비율과 밀접한 관련이 있습니다 ($K/\mu$가 높을수록 연성이 좋음). 그러나 이 현상을 설명하는 미시적 메커니즘에 대한 정량적 이해는 부족했습니다.
• 연구 목적: 유한 압축성 (finite compressibility) 을 가진 탄성 - 소성 자동자 (elasto-plastic automaton) 모델을 통해, $K/\mu$ 비율이 비정질 물질의 순환 하중 (cyclic shear) 과 정방향 전단 (forward steady shear) 거동에 미치는 영향을 규명하는 것입니다.

2. 방법론 (Methodology)

• 모델: 공간을 동일한 전단 변환 사이트 (sites) 로 채우는 메조스코픽 탄성 - 소성 자동자 모델을 사용했습니다.
• 압축성 도입: 선형 탄성 방정식을 이산화 (discretization) 할 때, staggered grid (교차 격자) 방식을 도입하여 Eshelby 포함체 문제를 수치적으로 해결했습니다. 이를 통해 압축률 ($K$) 과 전단률 ($\mu$) 의 비율인 무차원 파라미터 $\alpha = 2K/(K+\mu)$를 모델에 명시적으로 포함시켰습니다.
  • $\alpha$는 $K/\mu=0$ (완전 압축성) 일 때 0 에서 $K/\mu \to \infty$ (비압축성) 일 때 2 로 단조 증가합니다.
• 시뮬레이션 조건:
  1. 순환 전단 (Cyclic Shear): 진폭 $\Gamma$를 가진 전후 반복 하중을 가하여, 탄성 영역 (ER), 가역 소성 영역 (RPR), 확산 영역 (DR) 사이의 전이를 관찰했습니다.
  2. 정방향 전단 (Forward Shear): 일정한 전단 하중을 가하여 정상 상태 유동 응력 ($\sigma_y$) 을 측정했습니다.
• 초기 조건: 무작위 소성 필드를 가진 급속 냉각 (rapid quench) 상태의 비정질 시스템을 초기 상태로 설정했습니다.

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

A. Eshelby 후방 응력 (Eshelby Back Stress, $\sigma_0$) 의 역할 규명

• 전단 변환이 발생한 후 포함체 (site) 에 작용하는 응력인 '후방 응력' $\sigma_0$는 $K/\mu$에 의존합니다.
• 결과: $K/\mu$가 증가할수록 (압축성이 감소할수록) $\sigma_0$는 더 큰 음의 값을 가집니다. 이는 주변 매트릭스가 변환된 영역을 더 강하게 억제한다는 것을 의미합니다.

B. 순환 전단 (Cyclic Shear) 거동

• 전이 진폭 ($\Gamma_0$) 의 변화: 탄성 영역 (ER) 에서 가역 소성 영역 (RPR) 으로 전이되는 임계 진폭 $\Gamma_0$는 $\sigma_0$와 다음과 같은 관계를 가집니다:
  $$\Gamma_0 = \frac{1 + \sigma_0}{2}$$
• 압축성의 영향: $K/\mu$가 증가하면 $\sigma_0$가 더 음수가 되어 $\Gamma_0$가 증가합니다. 즉, 압축성이 낮은 (비압축성에 가까운) 시스템일수록 더 큰 변형 진폭까지 순수 탄성 거동을 유지하며, 효과적으로 '단단해집니다' (harder).
• 변형 경화 (Hardening) 전이:
  • $\Gamma < \Gamma_0$: 등방성 경화 (Isotropic hardening) 거동을 보입니다.
  • $\Gamma > \Gamma_0$: 운동학적 경화 (Kinematic hardening) 거동으로 전이됩니다. 이는 $\Gamma_0$에서 정확히 발생하며, 기존에 임의의 경화 규칙 없이도 이러한 복잡한 경화 거동이 모델에서 자연스럽게 나타남을 보여줍니다.

C. 정방향 전단 (Forward Shear) 거동

• 정상 상태 유동 응력 ($\sigma_y$): $K/\mu$가 증가함에 따라 정상 상태 유동 응력 $\sigma_y$는 감소합니다. 즉, 압축성이 낮은 시스템이 더 낮은 응력에서 소성 유동을 시작합니다.
• 응력 분포의 붕괴: 정상 상태에서의 국부 응력 분포 $P(\sigma)$는 $\sigma_0$에서 특징적인 어깨 (shoulder) 를 보입니다. 이를 $\sigma' = (\sigma - \sigma_0)/(1 - \sigma_0)$로 재스케일링하면 모든 $K/\mu$ 값에 대해 분포가 완벽하게 중첩 (collapse) 됩니다. 이를 통해 $\sigma_y$가 $\sigma_0$에 의존하는 정량적 관계를 유도했습니다.

D. 응력 재분배 (Stress Redistribution) 의 메커니즘

• $K/\mu$가 증가하면 국부 항복 사건 이후의 응력 재분배 진폭은 커지지만, 사건으로 인한 순 응력 완화 (net stress relaxation) 양은 변하지 않습니다.
• 이는 평균장 이론 (Mean-field theories) 에서 응력 재분배 진폭을 임의의 매개변수로 취급하는 것과 대조적으로, 본 모델에서는 이 진폭이 탄성 상호작용 (Eshelby 해) 과 탄성 계수에서 자연스럽게 유도됨을 보여줍니다.

4. 의의 및 결론 (Significance)

1. 미시적 메커니즘의 정량적 설명: 비정질 합금의 연성과 $K/\mu$ 비율 간의 경험적 관계를, 'Eshelby 후방 응력'의 변화와 '응력 재분배 진폭'의 증가라는 물리적 메커니즘으로 정량적으로 설명했습니다.
2. 경화 거동의 자발적 발생: 별도의 미세 구조적 경화 규칙 (ad-hoc hardening rules) 을 추가하지 않고도, 순환 하중 진폭에 따라 등방성 경화에서 운동학적 경화로 전이되는 복잡한 거동이 모델에서 자연스럽게 나타남을 증명했습니다.
3. 일반화 가능성: 금속 유리뿐만 아니라 콜로이드 유리, 에멀전 등 $K/\mu$가 조성, 노화, 열처리 등에 따라 변할 수 있는 다양한 비정질 시스템의 거동을 이해하는 데 중요한 통찰을 제공합니다.
4. 이론적 기여: 평균장 이론에서 종종 임의의 파라미터로 처리되던 '응력 재분배 진폭'이 실제 탄성 계수 ($K, \mu$) 에 의해 결정됨을 보여주어, 비정질 물질의 소성 이론을 보다 물리적으로 정립하는 데 기여했습니다.

요약하자면, 이 연구는 비정질 물질의 압축성 ($K/\mu$) 이 Eshelby 후방 응력을 조절하여, 순환 하중에서의 탄성 한계 ($\Gamma_0$) 와 정방향 전단에서의 항복 응력 ($\sigma_y$) 을 결정하며, 동시에 복잡한 변형 경화 거동을 유도한다는 핵심 결론을 도출했습니다.