에너지 지형 기반 탄성소성 모델로 본 비정질 고체의 주기적 전단 거동

초록

본 연구는 메조스케일 블록을 에너지 지형상의 “메소스테이트”로 정의하고, 각 블록의 탄성·소성 전이를 확률적으로 선택하는 엘라스토플라스틱 모델을 구축한다. 유한요소법(FEM)으로 블록 간 탄성 상호작용을 계산하고, 다양한 초기 어닐링 상태와 전단 진폭에 대해 균일 전단 및 주기 전단을 시뮬레이션한다. 모델은 어닐링에 따른 취성‑연성 전이, 바우싱거 효과, 주기 전단에서의 임계 에너지와 전단 밴드 형성, 그리고 파괴 시간의 발산 등 실험·분자동역학 결과와 일치하는 현상을 재현한다. 또한 파라미터(밀도 상태, 플라스틱 변형 선택) 의 민감도가 결과에 큰 영향을 미침을 보여준다.

상세 분석

이 논문은 비정질 고체의 거동을 “에너지 지형 기반 엘라스토플라스틱 모델(EPM)”로 접근한다. 핵심 아이디어는 각 메조스케일 부피(메소블록)를 여러 메소스테이트(에너지 최소점) 중 하나에 할당하고, 그 에너지 E₀와 무응력 플라스틱 변형 γ₀에 따라 선형 탄성 포텐셜 E = E₀ + μ/2(γ‑γ₀)²를 부여한다. 메소스테이트는 Gaussian 형태의 밀도 상태(Ω(E₀), 평균 –0.5, 표준편차 0.1)에서 무작위로 선택되며, 깊은 에너지 최소일수록 변형 안정 범위가 넓어지는 물리적 근거를 반영한다. 전단이 가해져 현재 γ가 안정 구간을 벗어나면 “소성 전이”가 발생하고, 새로운 메소스테이트와 플라스틱 변형을 선택한다. 플라스틱 변형 선택 규칙은(1) 균일, (2) SGR(현재 γ에 고정), (3) 최대 변위(현재 γ와 가장 멀리) 세 가지를 제시하고, 본 연구에서는 3번을 주로 사용한다.

블록 간 탄성 상호작용은 유한요소법(FEM)으로 구현한다. FEM은 한 번의 단위 플라스틱 변형에 대한 응답 텐서를 사전 계산하고, 실제 변형장에 대해 선형 중첩을 적용함으로써 계산 비용을 크게 낮춘다. 이는 기존의 무한 매질 커널 방식과 달리 비균질 탄성 상수와 경계조건을 손쉽게 포함시킬 수 있다.

시뮬레이션은 초기 어닐링 정도를 부모 온도 Tₚ로 조절한다. 높은 Tₚ(덜 어닐링)일수록 메소스테이트 분포가 넓고, 초기 플라스틱 변형의 표준편차가 커진다. 이렇게 준비된 시료에 대해(1) 균일 전단, (2) 주기 전단을 수행한다.

① 균일 전단 결과는 어닐링에 따라 취성‑연성 전이가 나타난다. 잘 어닐링된 시료는 전단 초기에 응력 급증 후 전단 밴드가 형성되어 파괴가 일어나며, 전단 밴드의 성장은 √γ 스케일을 보인다. 반면 덜 어닐링된 시료는 전단 전역에 걸쳐 플라스틱 활동이 퍼져 연성 거동을 보이고, 최종적으로 응력은 0.737, 평균 E₀는 –0.525로 수렴한다. 바우싱거 효과도 재현되는데, 전방 전단 후 역전단을 하면 역방향 강성이 감소한다. 이는 전방 전단 동안 낮은 안정 구간을 가진 메소스테이트가 소멸하고, 역방향에 대한 안정 구간이 상대적으로 넓어지기 때문이다.

② 주기 전단에서는 “임계 에너지” 개념이 부각된다. 포아리(포리) 어닐링된 시료(Tₚ≈3.16)는 전단 진폭 γ_max이 임계값 γ_yield≈0.43 이하이면 매 사이클마다 에너지가 감소해 흡수 상태(absorbing state)로 수렴한다. γ_max이 임계값을 초과하면 초기에는 에너지 감소가 이어지다 급격히 상승하면서 전단 밴드가 형성된다. 밴드는 이후 시스템 전체를 순회하며, 결국 모든 초기 정보를 소멸시킨다. 잘 어닐링된 시료는 더 큰 γ_max에서야 전단 밴드가 나타나며, 임계 진폭이 어닐링 정도에 따라 이동한다.

또한, 전단 밴드 형성 전후의 플라스틱 이벤트율(yield rate)은 비단조적이다. 초기에는 이벤트가 감소하다가 밴드가 형성되면 급증하고, 이후 일정한 수준으로 유지된다. 이는 실험·분자동역학 결과와 일치한다. 파괴까지 필요한 사이클 수 τ_f는 γ_max‑γ_yield에 대한 거듭 제곱근 법칙을 보이며, 포아리와 잘 어닐링된 경우 각각 β≈‑1.0, ‑2.2의 서로 다른 지수를 갖는다.

마지막으로 모델 파라미터(밀도 상태 형태, 플라스틱 변형 선택 규칙 등)에 대한 민감도 분석을 수행했으며, 특히 플라스틱 변형 선택이 “최대 변위”일 때 전단 밴드가 장기간 고정되는 트렌칭 현상이 나타난다. 이는 실제 유리에서 관찰되는 장기 잔류 전단 밴드와 유사할 수도 있지만, 추가 검증이 필요하다.

전반적으로 이 연구는 에너지 지형을 기반으로 한 메소스테이트 모델이 비정질 고체의 복잡한 전단 거동을 정량적으로 설명할 수 있음을 보여준다. FEM과 결합된 EPM은 다양한 초기 조건과 전단 프로토콜을 손쉽게 다룰 수 있는 강력한 프레임워크이며, 향후 실험 데이터와의 직접 비교 및 파라미터 최적화를 통해 비정질 물질 설계에 활용될 가능성이 크다.