⚡ electrical engineering

A discrete‑natured 3D DEM framework for spatial fracture analysis of unidirectional CFRP composites

본 연구는 단방향 CFRP 복합재의 횡방향 파괴를 효과적으로 모델링하기 위해 단위 셀을 연속체 요소로 취급하는 입방 충전 방식(cubic packing scheme)을 갖춘 새로운 3D 이산 요소법 프레임워크를 제시하며, 섬유 분포는 전역적 기계적 응답에 미치는 영향이 미미한 반면 섬유 부피 분율과 반경은 강성, 강도 및 손상 국부화에 결정적인 영향을 미친다는 것을 밝혀낸다.

원저자: Dechao Sun, Xiaoxuan Ding, Shuqi Zhang, Jianlin Liu, Zewen Gu, Xiaonan Hou

게시일 2026-09-11

📖 3 분 읽기☕ 가벼운 읽기

CC BY 4.0

원저자: Dechao Sun, Xiaoxuan Ding, Shuqi Zhang, Jianlin Liu, Zewen Gu, Xiaonan Hou

원본 논문은 CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/) 라이선스로 제공됩니다. ✨ 이것은 아래 논문에 대한 AI 생성 설명입니다. 저자가 작성하거나 승인한 것이 아닙니다. 기술적 정확성을 위해서는 원본 논문을 참조하세요. 전체 면책 조항 읽기

현대 공학은 탄소 섬유 강화 폴리머와 같이 강하면서도 가벼운 재료에 크게 의존하고 있습니다. 이러한 복합 재료는 얇고 단단한 탄소 실을 더 부드러운 플라스틱 접착제 안에 박아 넣어 만들어지며, 항공우주 및 자동차 분야에서 무거운 하중을 견딜 수 있는 구조를 형성합니다. 하지만 어떤 재료와 마찬가지로 이들에게도 약점이 있는데, 측면에서 잡아당길 때 섬유 사이의 플라스틱 접착제가 갈라질 수 있다는 점입니다. 횡방향 파괴(transverse fracture)라고 불리는 이러한 균열은 흔히 눈에 보이지 않는 아주 작은 결함에서 시작되어 점점 커지고 서로 연결되면서, 결국 구조 전체의 강도를 상실하게 만듭니다. 이러한 균열이 복잡한 섬유와 접착제의 미로 속에서 어떻게 형성되고 이동하는지 정확히 이해하는 것은 어렵습니다. 왜냐하면 손상이 일반적인 도구로는 쉽게 볼 수 없을 만큼 작은 규모에서 발생하며, 균열이 직선으로 이동하지 않기 때문입니다. 균열은 3차원 공간에서 뒤틀리고 꺾이며 섬유와 상호작용하므로, 부품이 언제 파손될지 예측하기가 어렵습니다.

이 퍼즐을 풀기 위해, 중국 석유 대학교의 첸차오 순(Dechao Sun)과 그의 팀은 랭커스터 대학교의 샤오난 허우(Xiaonan Hou)와 협력하여 컴퓨터상에서 이러한 재료를 시뮬레이션하는 새로운 방법을 개발했습니다. 연구진은 재료를 매끄럽고 단단한 하나의 덩어리로 취급하는 대신, 서로 맞닿고 밀어내는 방대한 양의 미세한 개별 입자들로 분해했습니다. 이 방식은 이산 요소법(discrete element method)이라고 불리며, 컴퓨터가 마치 수백만 개의 미세한 구성 블록으로 만들어진 것처럼 재료가 어떻게 작동하는지 관찰할 수 있게 해줍니다. 연구팀은 이 복합 재료의 정육면체에 대한 특정한 3차원 모델을 만들고, 그 안에 가상의 탄소 섬유와 에폭시 수지를 채워 넣었습니다. 그런 다음 그들은 이 입자들 사이의 미세한 연결 강도를 재료의 전체 강도와 연결하는 새로운 수학적 규칙을 고안했습니다. 이 방법의 핵심적인 혁신은, 컴퓨터에 미리 설정된 숫자를 강제로 입력하는 대신 입자들의 움직임을 관찰함으로써 재료가 압축될 때 옆으로 부풀어 오르는 성질인 포아송 비(Poisson's ratio)를 자동으로 생성한다는 점입니다.

연구진이 이 새로운 모델을 테스트했을 때, 그것이 재료의 실제 세계에서의 거동을 정확하게 재현할 수 있음을 발견했습니다. 시뮬레이션은 재료가 늘어나고 파괴되는 방식이 기존의 실험실 실험 및 이전의 2차원 컴퓨터 모델들과 일치함을 보여주었습니다. 그러나 3차원 관점은 기존 모델들이 놓쳤던 것을 밝혀냈습니다. 즉, 균열이 빵을 칼로 자르듯 한꺼번에 재료를 가로지르는 것이 아니라는 점입니다. 대신, 균열 전선(crack front)은 불균일하게 전진합니다. 어떤 지점에서는 균열이 빠르게 앞으로 나아가는 반면, 다른 지점에서는 뒤처지기도 하며, 재료의 두께를 따라 구불구불하고 뒤틀린 경로를 만들어냅니다. 연구진이 '균열 전선 곡절성(crack-front tortuosity)'이라 부르는 이 불균일한 진행은 오직 완전한 3차원 관점만이 포착할 수 있는 중요한 세부 사항입니다.

연구팀은 이 도구를 사용하여 섬유의 다양한 배열이 재료의 강도에 어떤 영향을 미치는지 조사했습니다. 그들은 플라스틱 접착제 내부의 섬유가 무작위로 배치되는 것이 재료의 내구성에 변화를 주는지 테스트했습니다. 연구 결과, 균열이 지나가는 정확한 경로는 섬유가 위치한 곳에 크게 좌우되었지만, 재료의 전반적인 강도와 강성은 배열에 관계없이 거의 동일하게 유지된다는 것을 발견했습니다. 이는 복합 재료의 일반적인 성능을 예측하는 데 있어, 섬유의 특정 무작위 패턴이 이전에 생각했던 것보다 덜 중요하다는 것을 시사합니다. 그러나 재료 속에 채워진 섬유의 양은 상당한 차이를 만들었습니다. 섬유 함량을 50%에서 68%로 높이면 재료는 더 단단하고 강해졌지만, 동시에 더 취성(brittle)이 강해져 더 낮은 신장률에서 파괴되었습니다. 연구진은 섬유가 너무 빽빽하게 채워지면 스트레스를 흡수할 접착제가 충분하지 않아 급격한 균열이 발생한다는 것을 관찰했습니다.

마지막으로, 연구진은 섬유의 양은 일정하게 유지한 채 섬유 자체의 크기를 살펴보았습니다. 그들은 섬유의 반지름, 즉 두께가 재료가 응력을 다루는 방식에 영향을 미친다는 것을 발견했습니다. 3.5 마이크로미터의 섬유 반경이 하중 전달과 접착제가 늘어나는 능력 사이에서 최적의 균형을 제공하는 가장 효율적인 구성임이 입증되었습니다. 시뮬레이션은 섬유의 크기가 섬유 간의 간격을 결정하며, 이 간격이 스트레스가 집중되는 곳과 균열이 시작되는 곳을 제어한다는 것을 보여주었습니다. 가상의 입자들 사이의 힘을 매핑함으로써, 연구팀은 하중이 접착제에서 섬유로 어떻게 이동하는지, 그리고 재료가 최종적으로 어디서 무너지는지를 정확히 볼 수 있었습니다. 이 새로운 3차원 프레임워크는 엔지니어들이 복합 재료의 숨겨진 역학을 이해할 수 있는 강력한 도구를 제공하며, 이러한 복잡한 구조가 어떻게 실패하는지에 대한 더 명확한 그림을 제시하여 더 오래 지속되도록 설계하는 데 도움을 줍니다.

기술 요약: 단방향 CFRP 복합재의 공간적 파괴 분석을 위한 이산적 특성의 3D DEM 프레임워크

문제 정의
탄소섬유 강화 폴리머(CFRP)는 항공우주 및 자동차 분야에서 매우 중요하지만, 특히 횡방향 균열(transverse cracking)과 관련하여 제한된 손상 허용 성능을 보입니다. 이러한 실패 모드는 미세 균열을 생성하여 층 두께를 통해 전파되며, 이는 강성 저하와 구조적 결함으로 이어집니다. 연속체 손상 역학(CDM), 응집 구역 모델(CZM), 확장 유한 요소법(XFEM)과 같은 기존의 수치적 접근 방식들이 이 분야를 발전시켜 왔으나, 이들은 종종 미리 정의된 균열 경로에 의존하거나 다수의 교차하는 3차원 균열의 수렴성 및 복잡성을 처리하는 데 어려움을 겪습니다. 또한, 기존의 이산 요소법(DEM) 연구는 주로 2D 단면 분석에 국한되어 있어, 두께 방향의 균열 연속성, 3D 손상 연결성, 그리고 불균질한 미세 구조에서 발생하는 공간적 응력 재분배를 포착하지 못합니다. 따라서 등가 이방성 물성이나 균질화된 구조에 의존하지 않고 단방향 CFRP의 횡방향 균열을 시뮬레이션할 수 있는 명시적인 3D DEM 모델에 대한 공백이 존재합니다.

방법론
저자들은 탄소섬유 강화 에폭시 복합재의 횡방향 파괴를 조사하기 위해 새로운 3도(3D) 이산 요소법(DEM) 프레임워크를 제안합니다. 본 방법론은 다음의 기둥을 기반으로 구축되었습니다:

이산 모델링: 3D 입방 충전 방식(cubic packing scheme)을 사용하여 49 µm 크기의 정육면체 대표 체적 요소(RVE)를 구성하였습니다. 미세 구조는 에폭시 매트릭스 내에 삽입된 38개의 원형 섬유(반지름 3.5 µm)로 구성되며, 상호작용은 6개의 직교 방향을 따르는 중심력 스프링(central-force springs)을 통해 모델링됩니다.
이산적 특성 교정(Discrete-Natured Calibration): 결합(bond)을 법선( $k_n$ $k_{n}$ ) 및 접선( $k_s$ $k_{s}$ ) 강성을 동시에 교정해야 하는 연속체 요소로 취급하는 전통적인 결합 입자 모델과 달리, 본 연구는 혁신적인 에너지 등가 관계를 도입합니다. 이 모델은 각 단위 셀을 연속체 요소로 취급하며, 거시적 포아송 비(Poisson's ratio)는 규정되는 것이 아니라 셀 변형에 의해 자동으로 생성됩니다.
- 프레임워크는 미시적 법선 강도( $k_n$ )와 거시적 영률( $E$ ) 및 포아송 비( $\nu$ )를 연결하는 명시적인 해석적 관계(식 41 및 42)를 도출합니다.
- 이 접근 방식은 미시적 접촉력과 거시적 연속체 변형 사이의 에너지 등가를 충족하기 위해 접선 강도( $k_s$ )를 교정 과정에서 배제하고 오직 법선 강도에만 의존함으로써, $k_s$ 의 필요성을 제거합니다.
검증: 교정된 모델은 실험 데이터 및 2D DEM 결과와 비교하여 검증됩니다. 시뮬레이션은 준정적 변위 제어 인장 하중 조건에서 수행되어 손상 진화, 균열 전파 및 거시적 응력-변형률 응답을 관찰합니다.

주요 결과
본 연구는 무작위 섬유 분포, 섬유 부피 분율( $V_f$ ), 섬유 반지름( $r$ )의 효과에 대한 파라미터 조사를 수행합니다:

모델 검증: 3D DEM은 순수 에폭시 및 CFRP RVE의 거시적 기계적 거동을 성공적으로 재현하였으며, 시뮬레이션된 영률과 포아송 비는 목표값에서 2% 미만의 편차를 보였습니다. 모델은 두께 방향의 균열 연속성과 균열 전면의 굴곡성(tortuosity)을 포착하였으며, 이를 통해 균열 전파가 2D 모델의 평면 궤적과 달리 섬유 길이를 따라 비동시적으로 발생함을 보여주었습니다.
섬유 분포: 무작위 섬유 분포는 국부적인 균열 개시 지점과 최종 균열 경로(예: 조밀하게 밀집된 섬유 사이의 변동)를 크게 변화시키지만, 전역적인 기계적 응답에는 미미한 영향을 미칩니다. 통계적으로 동등한 RVE들은 거의 동일한 응력-변형률 거동과 영률 진화를 나타냅니다.
섬유 부피 분율 ( $V_f$ ): $V_f$ 를 증가시키면(50%에서 68%로) 횡방향 강성과 피크 강도는 향상되지만, 파단 변형률은 감소하여 더 취성적인 포스트 피크(post-peak) 응답을 유도합니다. 높은 $V_f$ 는 수지 부족을 유발하여 국부적인 응력 집중과 조기 매트릭스 주도 손상을 촉진합니다.
섬유 반지름: 고정된 부피 분율( $V_f = 60\%$ )에서 섬유 반지름은 계면 응력 전달과 매트릭스 변형 수용 사이의 균형을 결정합니다. 섬유 반지름 $r = 3.5$ µm인 구성이 $r = 3.2$ µm 및 $r = 3.8$ µm와 비교했을 때 가장 효율적인 하중 지지 능력과 가장 높은 강성을 입증했습니다.

의의 및 주장
본 논문은 제안된 3D DEM 프레임워크가 단방향 CFRP를 위한 다척도 점진적 손상 분석의 효과적인 도구가 될 수 있다고 주장합니다. 그 주요 의의는 다음과 같습니다:

이산적 특성 교정: DEM의 불연속적 특성을 존중하면서도 포아송 비를 자연스럽게 생성하는, 미시적 파라미터에서 거시적 물성으로의 엄격한 에너지 기반 매핑을 제공합니다.
공간적 파괴 분석: 두께 방향의 균열 관통, 균열 전면의 굴곡성, 손상 합체(coalescence)를 포함한 복잡한 3D 파괴 메커니즘을 성공적으로 포착하며, 이는 2D 단순화로는 달성할 수 없는 것입니다.
메커니즘적 통찰: 섬유 분포가 국부적인 파괴 패턴을 결정하는 반면, 섬유 부피 분율이 전역적인 기계적 성능을 제어하는 지배적인 요인임을 밝혀냈습니다.

저자들은 현재의 모델이 단방향 적층판의 횡방향 파괴를 효과적으로 다루고 있으나, 향유 다방향 섬유 적층 구조, 정밀한 섬유-매트릭스 계면 구성 방정식, 그리고 다척도 손상 결합 전략을 통합하기 위한 향후 연구가 필요하다고 결론짓습니다.

유사한 논문