적응형 준연속체 기법으로 폴리머 네트워크의 파괴를 모델링하다

초록

본 연구는 네트워크 구조를 가진 재료(예: 폴리머, 겔)의 다중 스케일 파괴 현상을 효율적으로 모델링하기 위한 적응형 준연속체(QC) 수치 기법을 제안한다. 균열 선단 등 국부적 고변형률 영역에서는 이산적인 폴리머 체인을 명시적으로 모델링하고, 변형이 균일한 영역에서는 계산적 균질화를 통해 연속체로 근사한다. 동적 적응성을 통해 이 두 해상도 간 전환이 가능하며, 이를 통해 파괴 에너지에 대한 사전 가정 없이 대규모 시스템에서도 정확한 균열 전파 예측이 가능하다.

상세 분석

이 논문이 제안하는 방법론의 기술적 핵심은 기존 원자尺度 시뮬레이션용으로 개발된 준연속체(Quasi-Continuum, QC) 방법을 비선형 탄성 거동과 대변형을 수반하는 불규칙 폴리머 네트워크로 확장한 데 있다. 기존 QC 방법의 철학—고관심 영역의 이산적 모델링과 저관심 영역의 효율적 균질화—을 유지하면서, 폴리머 네트워크의 물리적 특성에 맞춰 두 가지 중요한 발전을 이루었다.

첫째, **미시-거시 에너지 일관성 조건(Hill-Mandel 조건)**을 기반으로 한 새로운 균질화 체계를 도입했다. 이는 미시尺度(개별 폴리머 체인)에서의 가상일(virtual work) 변화량과 거시尺度(균질화된 연속체)에서의 가상일 변화량이 동일해야 한다는 원리이다. 이를 통해 균질화된 1차 피올라-키르히호프 응력 텐서와 물성 탄성 텐서를 유도하며, 이 텐서들은 기저 네트워크 구조의 이방성을 그대로 반영한다. 이 접근법은 작은 변형과 선형 재료를 가정한 기존 균질화 연구를 넘어, 폴리머 체인의 본질적인 비선형성(예: Worm-Like Chain 모델)과 유한 변형을 포괄한다.

둘째, 동적 적응성(Dynamic Adaptivity) 알고리즘을 통한 효율성 극대화이다. 균열이 전파됨에 따라 고변형률 영역이 이동하면, 해당 영역을 중심으로 한 새로운 ‘완전 해상도(fully-resolved)’ 영역이 자동으로 정의되고, 변형이 완화된 영역은 다시 균질화 영역으로 전환된다. 이는 계산 자원을 항상 가장 필요한 곳에 집중시켜, 순수 이산 모델링의 계산 비용 문제를 해결하면서도 국부적 파괴 메커니즘의 정확한 포착을 가능하게 한다.

이 방법의 주요 장점은 파괴 역학에 대한 사전 지식의 불필요성에 있다. 위상 장 이론(Phase Field) 등 많은 연속체 파괴 모델이 요구하는 ‘규칙화 길이尺度(regularization length scale)‘나 임계 파괴 에너지와 같은 현상론적 매개변수를 필요로 하지 않는다. 대신, 파괴는 미시尺度 네트워크 토폴로지(가교 밀도, 결합 유형)와 개별 체인의 비선형 응답에서 자연스럽게 발생한다. 따라서 이 방법은 네트워크 설계(예: 이상적 공유결합 네트워크)와 그 파괴 저항성 간의 명시적인 관계를 탐구하는 데 매우 유용한 도구가 될 수 있다.