다상 변형플라스틱 재료의 미세구조 진화 변분 모델

초록

**
본 논문은 선형 운동경화와 무한소 변형 이론을 기반으로, 상전이와 미세구조 변화를 동시에 기술하는 변분 모델을 제시한다. 즉시 전이를 나타내는 ‘소산 거리’와 부피분율을 Young 측정으로 연속적으로 표현하여, 두 차원 FEM 시뮬레이션을 통해 모델을 검증한다.

**

상세 분석

**
이 연구는 기존의 단일상 변형플라스틱 이론을 다상 시스템에 확장함으로써, 상전이와 미세구조 형성을 하나의 통합된 프레임워크 안에 포함시켰다. 핵심은 두 가지 변분 원리, 즉 자유에너지 최소화와 소산 퍼텐셜 최소화를 동시에 적용한 점이다. 자유에너지 Ψ는 각 상의 탄성·경화 에너지와 상별 상수 c_i 를 가중 평균한 형태로 정의되며, 소산 퍼텐셜 Δ는 플라스틱 변형, 상전이, 점성 정규화 세 부분으로 분리된다. 특히 상전이 항 Δ_trans 에서는 ‘소산 거리’ D_i 를 도입해 즉시 전이를 경로 독립적으로 기술하고, 전이 매트릭스 g_ij 로 전이율을 확률적·비대칭적으로 표현한다.

Young 측정의 도입은 미세구조의 연속적인 진화를 수학적으로 정당화한다. 내부 변수 z 의 분포를 λ_i·δ(z−z_i) 형태의 이산 측정으로 나타내어, 부피분율 λ_i 와 상별 플라스틱 변형 ε^p_i 를 동시에 최적화한다. 이를 통해 유효 강성 텐서 C_eff 와 유효 플라스틱 변형 ε^p_eff 를 도출하고, 전체 시스템의 이완 자유에너지 Ψ_rel 를 구한다.

변분 라그랑지안 L 은 Ψ_rel·̇ + Δ_reg + Δ_trans + Δ_plast 로 구성되며, 변분 조건 0∈∂L/∂·ε^p_i , 0∈∂L/∂g_ij 로부터 플라스틱 흐름 법칙과 상전이 유동 법칙을 얻는다. 플라스틱 흐름은 전형적인 비선형 쿠흔‑터커 조건을 만족하고, 상전이 유동은 두 상 사이의 자유에너지 차와 전이 저항 r_i·‖ε^p_j−ε^p_i‖ 에 의해 구동된다. 전이율 g_ij 은 양의 부분만 허용하도록 설계돼, 동시에 두 방향 전이가 일어나지 않도록 보장한다.

새로운 상이 생성될 때의 초기 플라스틱 변형 ε^p_j,ini 은 해당 전이의 수율 함수 ϕ_ij 를 최대화하는 값으로 정의한다. 이는 열역학적 가능성에 기반한 물리적 직관을 반영한다.

수치 구현에서는 시간 적분 알고리즘을 제시하고, 1‑D 재료점 실험과 2‑D 원형 구멍이 있는 판의 압축 시험을 통해 모델을 검증한다. 결과는 전이 시점에서 부피분율이 급격히 변하고, 플라스틱 변형과 응력‑변형 곡선이 상전이와 강하게 연관됨을 보여준다. 또한 점성 파라미터 η_1, η_2 가 전이 속도와 수렴성을 조절한다는 점이 확인되었다.

전반적으로 이 논문은 변분 원리를 다상 변형플라스틱 문제에 성공적으로 적용함으로써, 상전이와 미세구조 진화를 동시에 예측할 수 있는 강력한 수학·수치적 도구를 제공한다. 향후 복합재, 금속 유리, 비정질 실리카 등 다양한 재료의 복합 변형 거동을 해석하는 데 활용될 가능성이 크다.

**