A variational multi-phase model for elastoplastic materials with microstructure evolution

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

• 배경: 다양한 재료 (토양, 비정질 금속, 유리질 실리카 등) 에서 전단대 (shear bands) 와 같은 미세구조 형성이 상변태 (phase transformations) 와 밀접하게 연관되어 발생합니다. 예를 들어, 고압 하에서 실리카의 밀도 변화나 비정질 금속의 국소적 상 분리 현상 등이 관찰됩니다.
• 문제점: 이러한 미세구조의 진화와 상변태를 연속체 역학 관점에서 설명하는 기존 모델들은 한계가 있습니다. 특히, 상변태가 순간적으로 일어나는지, 아니면 시간에 따라 연속적으로 진화하는지에 대한 수학적 형식화가 명확하지 않으며, 확률적 과정과 변분 원리를 통합하여 미세구조의 시간 연속적 진화를 설명하는 모델이 부족했습니다.
• 목표: 본 연구는 선형 운동 경화 (linear kinematic hardening) 를 갖는 탄소성 재료를 대상으로, 상변태와 연관된 미세구조 진화를 설명할 수 있는 일반적인 변분 다상 모델 (Variational Multi-phase Model) 을 제안하는 것입니다.

2. 방법론 (Methodology)

본 연구는 변분 원리 (Variational Principles) 에 기반한 열역학적 프레임워크를 사용하여 모델을 구성했습니다.

• 기본 가정:
  • 무한소 변형 이론 (Infinitesimal strain theory) 적용.
  • 등온 과정 가정.
  • 재료는 $k$개의 가능한 상 (phases) 을 가질 수 있으며, 각 상의 부피 분율 ($\lambda_i$) 로 표현됨.
• 에너지 함수 및 소산 (Dissipation):
  • 헬름홀츠 자유 에너지 밀도 ($\Psi$): 각 상의 탄성 에너지와 선형 운동 경화 항을 포함하며, 전체 에너지는 각 상의 부피 분율에 따른 가중 평균으로 정의됨.
  • 소산 포텐셜 ($\Delta$):
    1. 소성 변형 ($\Delta_{plast}$): 각 상 내의 소성 변형률 속도에 비례.
    2. 상변태 ($\Delta_{trans}$): 상 간의 전이를 설명하기 위해 소산 거리 (Dissipation Distance) 개념 도입. 이는 시간 내에서의 즉각적인 상 전이를 확률 행렬 (Transition Matrix, $g_{ij}$) 을 통해 모델링함.
    3. 규제 항 ($\Delta_{reg}$): 수학적 처리의 용이성과 계산적 안정성을 위해 점성 규제 (Viscous regularization) 항 추가.
• 미세구조 진화 및 완화 에너지 (Relaxed Energy):
  • 영 측정 (Young Measures): 새로운 상의 부피 분율을 나타내기 위해 영 측정을 사용하여 시간 연속적인 미세구조 진화를 표현.
  • 완화된 에너지 ($\Psi_{rel}$): 전체 에너지를 최소화하여 유효 강성 텐서 ($C_{eff}$) 와 유효 소성 변형률 ($\varepsilon^p_{eff}$) 을 도출. 이는 Reuß 하한선 (homogenous Reuß lower bound) 형식과 유사함.
• 진화 방정식 유도:
  • 라그랑지안 ($L = \dot{\Psi} + \Delta$) 을 내부 변수 ($\dot{\varepsilon}^p_i$) 와 전이율 ($g_{ij}$) 에 대해 최소화하여 미분 포함식 (Differential Inclusions) 을 유도.
  • 이를 통해 소성 변형률 진화 방정식과 상 전이율 방정식을 명시적으로 도출.
  • 전이율 방정식은 두 상 간의 에너지 차이와 소성 변형률 차이, 그리고 소산 거리를 기반으로 하여, 한 번의 시간 단계에서 양방향 전이가 동시에 일어나지 않도록 보장함.

3. 주요 기여 (Key Contributions)

1. 변분 기반의 다상 모델 통합: 확률적 전이 행렬 (Stochastic probability matrix) 과 소산 거리 (Dissipation distance) 개념을 결합하여, 탄소성 재료 내에서의 상변태와 미세구조 진화를 통합적으로 기술하는 새로운 변분 모델을 제시함.
2. 시간 연속적 미세구조 표현: 영 측정 (Young measures) 을 도입하여 이산적인 상 전이를 시간적으로 연속적인 미세구조 진화 과정으로 모델링함.
3. 새로운 상의 생성 메커니즘: 열역학적으로 가능한 경우 새로운 상이 즉시 생성된다고 가정하고, 초기 소성 변형 값을 결정하는 수학적 조건을 제시함.
4. 수치적 검증: 유한요소법 (FEM) 을 이용한 2 차원 벤치마크 테스트를 통해 모델의 타당성을 입증함.

4. 결과 (Results)

• 재료 점 거동 (Material Point Behavior):
  • 1 차원 인장 - 압축 사이클 테스트에서, 항복 응력에 도달하자마자 새로운 상 ($\lambda_1, \lambda_2$) 의 부피 분율이 즉시 증가하는 것을 확인.
  • 하중이 제거되는 구간에서는 진화가 멈추고, 하중이 다시 가해지면 기존 상이 더 발전함.
  • 쇼크다운 (Shakedown) 거동: 하중이 반복됨에 따라 상 분율이 일정 값에 수렴하여 더 이상 새로운 상이 생성되지 않고 순수 탄성 거동만 나타나는 현상 관찰.
• 2 차원 유한요소 벤치마크 테스트:
  • 원형 구멍이 있는 정사각형 판에 압축 하중을 가하는 시뮬레이션 수행.
  • 상변태 시작 위치: 구멍 주변 (응력 집중 영역) 에서 상변태가 시작되어 시료 전체로 확산됨.
  • 수렴성: 3 가지 다른 메쉬 (Coarse, Fine, Superfine) 를 사용하여 계산한 결과, 응력 - 변형률 거동 및 상 진화 곡선이 메쉬 정련에 따라 수렴하는 것을 확인하여 모델의 수치적 안정성 입증.
  • 최종적으로 약 54% ($\lambda_3 \approx 0.54$) 의 부피가 상변태를 겪는 것으로 나타남.

5. 의의 및 결론 (Significance and Conclusion)

• 이론적 의의: 복잡한 미세구조 진화와 상변태를 제어할 수 있는 전이 행렬 (Transition Matrix) 을 통해 재료의 거동을 예측할 수 있음을 보임. 이는 기존 모델들이 간과했던 미세구조의 역학적 진화 과정을 정량적으로 설명할 수 있는 틀을 제공함.
• 실용적 의의: 비정질 금속, 유리질 실리카, 입상 물질 등 다양한 재료의 국소적 변형 및 파괴 메커니즘을 이해하는 데 기여할 수 있음.
• 향후 과제:
  • 현재 모델은 등온 상황을 가정하고 있으나, 향후 온도 의존성 및 열적 효과를 포함할 수 있음.
  • 전이 행렬을 지수 함수 형태의 내부 상호작용에 의존하도록 확장하여 더 정교한 물리적 현상 (예: [16, 17] 에서 제안된 접근법) 을 모사할 수 있음.
  • 라미네이트 구조나 기타 복잡한 미세구조 형상을 모델에 포함할 수 있는 가능성 제시.

요약하자면, 이 논문은 탄소성 재료의 미세구조 진화를 설명하기 위해 변분 원리, 소산 거리, 그리고 영 측정을 결합한 강력한 수학적 모델을 제시하였으며, 이를 통해 상변태와 미세구조 형성의 역학적 메커니즘을 정량적으로 예측할 수 있음을 수치 시뮬레이션을 통해 입증했습니다.