고에너지 재료의 미세역학 기반 열탄성 특성 예측

초록

본 연구는 고체 폭약인 PBX 9501과 같은 고에너지 재료의 입자 복합체에서, 실험이 어려운 열탄성 계수를 미세역학 모델을 이용해 예측한다. 높은 입자 부피비(>90 %)와 큰 탄성계수 대비(5,000–10,000)를 고려해, 엄격한 경계값, 평균장 이론, 일반화 셀 방법(Generalized Method of Cells, GMC) 및 재귀 셀 방법(Recursive Cell Method, RCM)을 적용하고, 유한요소법(FEM) 및 실험 데이터와 비교 검증한다. 결과적으로 RCM이 계산 효율성과 정확도 면에서 가장 유망함을 확인하였다.

상세 분석

본 논문은 고에너지 물질, 특히 폴리머 결합 폭약(PBX) 9501의 열탄성 거동을 미세역학적으로 해석하는 데 초점을 맞추었다. PBX는 폭약 결정 입자가 고무성 바인더에 90 % 이상 부피비로 분산된 두 성분 복합체이며, 입자와 바인더 간 탄성계수 차이가 5,000~10,000배에 달한다는 특징을 가진다. 이러한 극단적인 물성 대비와 높은 입자 함량은 전통적인 평균장 이론이나 단순 혼합 규칙으로는 정확한 유효 탄성계수를 예측하기 어렵게 만든다.

먼저 저자는 Hill, Hashin‑Shtrikman, Voigt‑Reuss 등 고전적인 엄격 경계값을 정리하고, PBX 9501의 실제 물성(입자와 바인더의 Young’s modulus, Poisson 비, 열팽창계수 등)을 대입해 이론적 한계를 산출하였다. 경계값은 실험값을 포괄하지만, 그 폭이 너무 넓어 설계에 직접 활용하기엔 한계가 있다.

다음으로 평균장 기반의 분석 모델, 즉 Mori‑Tanaka와 Self‑Consistent Scheme을 적용하였다. 이들 모델은 입자 형태를 구형으로 가정하고, 상호작용을 평균화하지만, 고부피비와 높은 강성 대비에서는 입자 간 접촉 및 응력 집중 현상을 충분히 반영하지 못한다는 점을 지적한다. 실제로 계산된 유효 Young’s modulus는 실험값보다 15~30 % 낮게 나타났다.

연구의 핵심은 두 가지 수치 동질화 기법, Generalized Method of Cells(GMC)와 Recursive Cell Method(RCM)이다. GMC는 복합체를 규칙적인 셀 네트워크로 분할하고, 각 셀 내에서 평균장과 변형률을 연계해 전체 응답을 도출한다. 이 방법은 셀 수가 늘어날수록 정확도가 향상되지만, 고부피비와 비균질 입자 배열을 재현하려면 셀 크기가 매우 작아져 계산 비용이 급증한다.

RCM은 이러한 문제를 해결하기 위해 다중 스케일링 접근을 채택한다. 먼저 작은 대표 부피셀(RVE)을 상세 FEM으로 해석해 유효 탄성 텐서를 얻고, 이를 상위 레벨 셀의 재료 특성으로 재귀적으로 전달한다. 각 단계에서 필요한 자유도는 크게 감소하므로 전체 계산 시간이 GMC 대비 5~10배 단축된다. 또한, 입자 간 접촉을 직접 모델링하므로 응력 집중과 국부 변형을 보다 정확히 포착한다.

검증 결과, RCM이 예측한 유효 Young’s modulus와 열팽창계수는 실험값과 3 % 이내 차이를 보였으며, FEM 전체 모델(수십만 자유도)과 비교했을 때 오차는 12 %에 불과했다. 반면, GMC는 동일한 셀 해상도에서 812 % 정도의 오차를 나타냈다. 열탄성 특성에서도 RCM은 온도 의존성을 포함한 복합 응답을 안정적으로 재현했으며, 고온(≈80 °C)에서의 탄성 저하를 정확히 예측하였다.

이러한 결과는 고에너지 복합재의 설계 단계에서 실험적 위험을 최소화하면서도 신뢰성 있는 물성 데이터를 제공할 수 있음을 의미한다. 특히, RCM은 고부피비·고강성 대비 복합재의 다중 물리적 거동(기계·열·동적)을 통합적으로 해석할 수 있는 효율적인 플랫폼으로 활용 가능하다. 다만, 입자 형태가 비구형이거나, 바인더가 비선형 점탄성인 경우 추가적인 모델 확장이 필요하며, RCM의 재귀 단계에서 경계 조건 선택이 결과에 민감하게 작용한다는 점도 주의해야 한다.