A Finite Element Process Modeling Framework for Predicting Residual Stresses and Microcracking in Ceramic Matrix Composites

기술 요약: 세라믹 매트릭스 복합재의 잔류 응력 및 미세 균열 예측을 위한 유한 요소 공정 모델링 프레임워크

문제 정의
세라믹 매트릭스 복합재(CMC)는 1200–1600 °C 이상의 온도에서도 탁월한 열 저항성과 구조적 무결성을 보유하고 있어 차세대 극초음속 응용 분야에 매우 중요하다. 그러나 주요 제조 방법인 고분자 함침 및 열분해(PIP) 공정은 미세 균열 및 기공과 같은 상당한 결함 형성을 야기한다. 이러한 결함은 완전한 치밀화를 달성하기 위해 비용이 많이 드는 다수의 재함침 사이클을 필요로 한다. 근본 원인은 복잡하게 결합된 열화학 반응과 열분해 과정 중 전세라믹 고분자(PCP) 내에서 발생하는 기계적 특성의 동시 진화에 있다. 이러한 공정은 실험적으로 정량화하기 어려운 거대한 잔류 응력을 생성하는데, 이는 극심한 공정 온도와 반응 역학, 질량 손실, 치밀화, 체적 수축 및 가스 수송이 동시에 발생하는 특성 때문이다. 따라서 이러한 결합된 메커니즘을 시뮬레이션하고 결함 형성을 예측할 수 있는 물리 기반의 계산 프레임워크가 절실히 필요하다.

방법론
본 연구는 CMC의 미시적 시뮬레이션에 초점을 맞춘, 디지털 제조를 위한 물리 기반 유한 요소(FE) 프레임워크를 소개한다. 이 접근 방식은 경화(C-staging)를 이후의 열분해 단계와 직접 연결하는 새로운 2단계 공정을 활용한다.

• 지배 물리 법칙: 프레임워크는 세 가지 결합된 방정식을 해결한다: 즉각적인 선형 탄성 구성 방정식, 발열 항을 포함한 과도 열 방정식, 그리고 경화 및 열분해를 위한 역학 모델이다. 총 변형률은 기계적, 열적, 수축 성분으로 분해된다. 수축은 경화 단계에서는 경화도에 대한 선형 함수로, 열분해 단계에서는 질량 손실 및 치밀화로부터 유도된 체적 수율의 함수로 모델링된다.
• 손상 모델링: In-situ 미세 균열을 포착하기 위해 Bažant과 Oh의 크랙 밴드 이론에 기반한 점진적 손상 모델이 통합되었다. 손상 개시는 최대 주응력 기준을 따르며, 피크 이후의 거동은 변형 연화(strain softening)를 동반한 선형 트랙션-분리(linear traction-separation)를 통해 모델링된다. 스칼라 손상 매개 변수($D$)가 컴플라이언스 텐서를 저하시킨다.
• 구현: 지배 방정식은 사용자 정의 서브루틴(UMAT 및 UMATHT)을 사용하여 Abaqus FE 환경에 구현되었다. 프레임워크는 2단계 분석을 채택한다: 경화 분석이 초기 상태(잔류 응력, 손상, 재료 특성)를 설정하며, 이는 열분해 분석을 위한 초기 조건으로 임포트된다.
• 재료 모델링: 특정 시스템(SMP-10 매트릭스 및 AS4 탄소 섬유)에 대한 제한된 실험 데이터로 인해, 본 연구는 대치 재료 기술(stand-in material description)을 구축한다. 경화 특성은 유사 시스템(phenolic SC-1008 및 epoxy RIM R135-H1366)으로부터 도출되었으며, 열분해 역학 및 수축은 Key 등이 제안한 SMP-10의 확립된 모델을 기반으로 한다. 열분해 중의 재료 특성은 변화하는 세라믹 밀도의 함수로 진화한다고 가정된다.
• 기하 구조: 프레임워크는 주기적 경계 조건을 적용하여, SMP-10 매트릭스 내에 20개의 무작위로 분산된 AS4 탄소 섬유(체적 분율 45%)를 포함하는 대표 체적 요소(RVE)에서 시연된다.

주요 기여

1. 통합 2단계 프레임워크: 본 논문은 경화와 열분해를 직접 연결하여, 경화 중에 축적된 잔류 응력 상태와 손상이 열분해 시뮬레이션의 초기 조건 역할을 할 수 있도록 하는 직접적인 연결 고리를 제시한다. 이는 열분해 단계를 별도로 처리하거나 섬유가 명시적으로 모델링되지 않은 메조스케일에 집중했던 기존 연구들과 대조된다.
2. 점진적 손상 통합: 점진적 손상 모델의 포함은 두 공정 단계 전체에서 in-situ 미세 균열을 시뮬레이션할 수 있게 하여, 응력 진화와 재료 저하 사이의 상호작용을 포착한다.
3. 물리적으로 일관된 가정: 완전한 재료 데이터가 부재한 상황에서, 저자들은 SMP-10 특성을 위한 견고한 프록시 모델을 개발하였으며, 밀도 변화와 검증된 역학 모델로부터 특성 진화를 도출함으로써 물리적 일관성을 확보하였다.
4. 미시적 해상도: 모델은 미시적 규모에서 섬유-매트릭스 계면을 명시적으로 분해하여, 더 큰 규모의 모델에서 흔히 균질화되는 응력 집중 현상과 손상 국부화를 관찰할 수 있게 한다.

결과
CMC 미세 구조 시뮬레이션 결과는 다음과 같다:

• 경화 단계: 경화 중 발생한 잔류 응력은 상대적으로 낮았다(램프 구간 중 인장 응력 약 7 MPa 피크, 냉각 후 상온에서 약 65 MPa). 손상 개시는 미미했으며, 체적 평균 손상 매개 변수($D$)는 경화 종료 시점에 0.015 미만을 유지했다. 응력 집중은 주로 섬유가 근접한 곳에서 관찰되었다.
• 열분해 단계: 매트릭스의 상당한 미세 균열은 열분해 유도 수축에 의해 발생했다. 손상은 열적 램프 초기(약 285°C, 급격한 분해가 시작되는 시점)에 시작되었다. 체적 평균 손상 매개 변수는 급격히 상승하여 최종 1200°C 목표치의 절반 지점에 도달하기 전에 0.7을 초 exceeded했다.
• 응력 및 수축: 피크 인장 응력은 열적 램프 초기(약 300–400°C)에 발생하였으며, 점진적 손상에 의해 추가적인 응력 축적이 제한되기 전 약 55 MPa에 도달했다. 시뮬레이션은 약 8%의 최종 체적 수축을 예측했다.
• 특성 진화: 최종 탄성 계수와 인장 강도는 누적된 손상으로 인해 완전히 치밀한 SiC의 이론적 벌크 특성에 비해 크게 저하되었다(각각 ~160 GPa 및 1.4 GPa). 결과적인 미세 구조는 문헌 관찰 내용과 일치하는 파편화 형태(fragmentation morphologies)를 나타냈다.

의의 및 주장
저자들은 본 연구를 CMC의 결함 형성을 예측하기 위해 경화와 열분해를 연결하는 개념 증명(proof-of-concept)으로서의 계산 프레임워크로 규정한다. 본 연구는 공정 조건(열 사이클, PIP 파라미터)이 잔류 응력, 수축 및 미세 균열에 어떻게 영향을 미치는지 체계적으로 조사할 수 있는 "예측적, 공정 정보 기반 모델링"의 토대를 구축했다고 주장한다.

논문은 제한된 재료 데이터로 인해 현재의 결과가 프록시 기술과 가정에 의존하고 있음을 겸허히 언급한다. 그러나 본 프레임워크가 결합된 결함 형성 메커니즘을 이해하기 위한 필수적인 플랫폼을 제공한다고 단언한다. 저자들은 더 정확한 재료 특성 데이터가 확보되고 열분해 가스 수송과 같은 추가적인 물리 현상이 통합됨에 따라, 이 접근 방식이 CMC의 공정 후 성능을 개선하기 위한 제조 전략 설계를 가이드할 잠재력이 있다고 결론짓는다. 본 연구는 즉각적으로 모든 최적화 문제를 해결한다고 주장하는 것이 아니라, 그러한 조사를 시작하기 위한 필수적인 계산 도구 세트를 제공하는 데 목적이 있다.