극한 동적 하중 하의 비정질 열가소성 폴리머 거동 통합 이론

초록

본 논문은 비정질 열가소성 폴리머, 특히 PMMA의 충격·스팔·용융·분해 등 극한 동적 하중 상황을 포괄적으로 기술하는 비선형 연속체 모델을 제시한다. 열역학적 내부 상태 변수와 전이(order) 파라미터를 결합해 열탄성·점탄·점소성·취성·연성 파괴·유리전이·충격분해를 동시에 다루며, 위상장 모델을 통해 균열·공극 성장까지 묘사한다. 모델은 고압 Hugoniot, 방출파 속도, 저압 점탄 파동 및 스팔 강도 예측에서 실험과 좋은 일치를 보이며, 다양한 온도·속도 조건에서의 물성 변화를 통합적으로 설명한다.

상세 분석

이 연구는 비정질 열가소성 폴리머의 복합적인 동적 거동을 하나의 연속체 이론으로 통합하려는 시도이다. 가장 큰 특징은 열역학적 내부 상태 변수(internal state variables, ISV)를 이용해 분자 사슬의 얽힘, 자유 부피, 체인 회전 등 미시구조 변화를 매크로 응답에 연결한 점이다. ISV는 점탄성(Viscoelasticity)과 점소성(Viscoplasticity) 각각에 대해 텐서형 변수와 스칼라 변수로 구분되며, 변형률 속도와 온도에 대한 민감도를 자연스럽게 포함한다.

변형도(Deformation gradient)는 총 열탄성 성분(F)과 잔류 플라스틱 성분(P)으로 분해되는데, 여기서 열탄성 성분은 열팽창과 점탄성을 동시에 포함한다. 기존의 복합 분해(F=Fe·Fv·Fp 등)와 달리, 현재 모델은 열탄성·점탄을 하나의 텐서로 취급함으로써 수식상의 복잡성을 크게 낮추면서도 모든 비가역 변형을 속도 형태(velocity gradient)에서 일관되게 기술한다.

열탄성 부분은 로그형 열탄성 변형 측정(logarithmic strain)과 온도·압력 의존성을 갖는 자유 에너지 함수를 사용해 고압 EOS를 자연스럽게 도출한다. 점탄성은 선형 점탄성 이론을 일반화한 내부 변수 기반의 Prony‑like 표현이 아니라, Holzapfel‑Simo 프레임워크에 기반한 텐서형 내적 변수를 통해 비선형 점탄성을 구현한다. 이는 고주파·고속 충격 상황에서 관측되는 이완 시간의 비선형 변화를 정확히 포착한다.

점소성은 등방성 및 이방성 경화·연화 메커니즘을 동시에 포함한다. 스칼라 변수는 등방성 경화(예: 온도·압력에 의한 강도 상승)를, 텐서 변수는 플라스틱 스트레치 텐서를 통해 전단‑축 방향 경화·연화를 기술한다. 플라스틱 흐름 규칙은 온도·압력 의존성을 갖는 연속적인 흐름 함수와, 유리전이·충격분해 임계값을 초과할 때 급격히 변하는 전이 함수를 결합한다.

파괴 메커니즘은 두 종류의 전이 파라미터(phase‑field order parameters)로 구분된다. 하나는 연성 파괴(전단 항복에 의한 국부 소성)용, 다른 하나는 취성 파괴(크레이징·공극 성장)용이다. 전이 파라미터는 스칼라 형태이며, 변형·온도·압력에 따라 연속적이거나 급격히 변하는 Kinetics를 갖는다. 파괴 에너지 밀도에 공극 부피 비율을 직접 연결함으로써, 크레이징에 따른 체적 팽창을 자연스럽게 모델링한다.

유리전이와 충격분해는 각각 별도의 스칼라 전이 파라미터로 다루며, 온도·압력·변형률 속도에 대한 임계 함수를 통해 물성 변화를 스위칭한다. 특히 충격분해는 화학 반응을 가정하지 않고, 압축에 의해 발생하는 급격한 체적 감소와 내부 에너지 급증을 반영한다.

모델 검증을 위해 PMMA에 적용했으며, 고압 Hugoniot(압력 0–120 GPa)에서 실험 데이터와 거의 일치하는 압축 곡선을 얻었다. 방출파 속도 역시 실험값을 재현했으며, 저압 점탄 파동 해석을 통해 이완 시간(τ)과 점탄성 모듈러스의 온도·압력 의존성을 정량화했다. 스팔 강도 예측에서는 초기 온도와 충격 속도에 따른 파괴 강도 변화를 정확히 포착했으며, 특히 온도가 유리전이 근처에서 급격히 감소하는 현상을 재현했다.

이론적 강점은 (1) 열·역학·점탄·점소성·파괴를 하나의 열역학적 프레임워크에 통합, (2) 내부 변수와 위상장 모델을 통한 물성 전이의 연속·불연속 제어, (3) 고압·고온·고속 조건에서도 수치적으로 안정적인 전개 가능성이다. 한계점으로는 파라미터 식별이 복잡하고, 다중 스케일(분자‑거시) 연결이 아직 정량적 검증이 부족하다는 점, 그리고 3‑D 전이 파라미터의 경계 조건 설정이 실험 데이터에 크게 의존한다는 점을 들 수 있다. 향후 연구에서는 파라미터 자동 추정, 실험적 마이크로‑스케일 검증, 그리고 복합 구조(다층·다재료) 적용을 목표로 할 필요가 있다.