A multi-phase-field model for fiber-reinforced composite laminates based on puck failure theory

이 논문은 Puck 파괴 이론을 기반으로 두 개의 위상장 변수와 메쉬 오버레이 기법을 활용하여 섬유 강화 복합재 적층판의 섬유 및 섬유간 파괴 모드를 정성적·정량적으로 예측하는 2 차원 다중 위상장 모델을 제안하고, 다양한 실험 결과와 비교하여 그 유효성을 입증합니다.

핵심

이 논문은 **"복합재료 (Fiber-Reinforced Composites) 가 어떻게 부서지는지 컴퓨터로 정확하게 예측하는 새로운 방법"**을 소개하고 있습니다.

복합재료는 탄소섬유와 수지를 섞어 만든 재질로, 비행기 (보잉 787, 에어버스 A350) 나 고성능 스포츠 용품에 많이 쓰입니다. 가볍고 튼튼하지만, 어떤 방향으로 힘을 가하면 어떻게 깨지는지 예측하기 매우 어렵습니다.

이 논문은 그 어려운 문제를 해결하기 위해 **'다중 위상장 (Multi-Phase-Field)'**이라는 기술을 사용했습니다. 이를 일상적인 언어와 비유로 설명해 드리겠습니다.

---

1. 핵심 아이디어: "두 개의 감시 카메라"

기존의 컴퓨터 시뮬레이션은 복합재료가 깨질 때, "섬유가 끊어지는 것"과 "섬유 사이의 수지 (접착제) 가 찢어지는 것"을 구분하기가 힘들었습니다. 마치 한 대의 카메라로 두 가지 다른 사건을 동시에 찍으려다 흐릿하게 보이는 것과 비슷합니다.

이 연구팀은 **두 개의 독립된 '감시 카메라 (위상장 변수)'**를 설치했습니다.

• 카메라 A (섬유 감시): 섬유가 끊어지는 순간을 포착합니다.
• 카메라 B (수지 감시): 섬유 사이의 접착제가 찢어지는 순간을 포착합니다.

이 두 카메라가 서로의 상황을 공유하며, **"어느 부분이 먼저 약해졌는지, 그리고 그로 인해 다른 부분이 어떻게 영향을 받는지"**를 정교하게 분석합니다.

2. 방법론: "투명한 오버레이 (Mesh Overlay)"

복합재료는 여러 겹의 천 (적층판) 을 쌓아 만듭니다. 각 층마다 섬유 방향이 다릅니다 (예: 0 도, 45 도, 90 도).
기존 방식은 이 모든 층을 3 차원 입체로 모델링해야 해서 컴퓨터가 매우 무거워졌습니다.

이 연구팀은 "투명한 오버레이 (겹쳐서 보기)" 방식을 썼습니다.

• 비유: 얇은 투명 시트 (층) 여러 장을 하나씩 겹쳐서 쌓는다고 상상해 보세요.
• 각 시트에는 서로 다른 방향의 섬유가 그려져 있습니다.
• 컴퓨터는 이 시트들이 하나의 공간에 겹쳐져 있지만, 각 시트의 상태 (손상 여부) 는 따로 계산합니다.
• 이렇게 하면 3 차원 모델보다 훨씬 가볍고 빠르게, 하지만 정확도는 유지하면서 층별 손상을 계산할 수 있습니다.

3. 규칙: "푸크 (Puck) 의 위험 신호"

손상이 언제 시작될지 결정하는 기준은 **'푸크 (Puck) 이론'**이라는 공학적 규칙을 따릅니다.

• 비유: 비행기 날개에 붙은 **'경고등'**과 같습니다.
• 이 경고등은 섬유가 당겨질 때 (인장) 와 눌릴 때 (압축) 에 따라, 그리고 섬유 사이가 당겨질 때와 눌릴 때에 따라 서로 다른 조건에서 켜집니다.
• 이 연구팀은 이 복잡한 경고등 시스템 (푸크 이론) 을 컴퓨터 시뮬레이션에 완벽하게 접목시켜, "이제부터는 섬유가 끊어질 거야", **"아니, 지금부터는 수지가 찢어질 거야"**라고 미리 예측하게 만들었습니다.

4. 검증: "실제 실험과의 대결"

이론만으로는 부족합니다. 연구팀은 실제 실험 데이터와 비교했습니다.

• 구멍 뚫린 판을 당기는 실험 (Open Hole Tension): 구멍 주변에서 어떻게 갈라지는지.
• 양쪽 끝이 잘린 판을 당기는 실험 (Double Notched Tension): 잘린 끝에서 균열이 어떻게 퍼지는지.
• 압축 실험: 누르면 어떻게 부러지는지.

결과적으로, 이 새로운 모델은 **실제 실험에서 관찰된 균열의 모양 (Qualitative)**과 **얼마나 많은 힘을 견디는지 (Quantitative)**를 매우 잘 예측했습니다. 특히 섬유가 끊어지는 과정과 수지가 찢어지는 과정이 서로 어떻게 경쟁하고 협력하며 파괴를 일으키는지 보여줍니다.

5. 요약 및 의의

이 논문은 **"복합재료의 파괴를 예측하는 새로운 지도"**를 만들었습니다.

• 기존의 문제: 복잡한 3D 모델은 너무 무겁고, 섬유와 수지 손상을 구분하기 어려웠다.
• 이 연구의 해결책: 두 개의 감시 카메라 (위상장) 와 투명한 오버레이 기법을 써서, 층별로 정교하게 분석하면서도 계산은 가볍게 했다.
• 결과: 비행기나 자동차 설계자가 "여기가 약할 것이다"라고 미리 알 수 있게 도와주어, 더 안전하고 가벼운 구조물을 설계하는 데 기여할 것입니다.

한 줄 요약:

"복합재료의 층층이 쌓인 구조 속에서, 섬유가 끊어지는 것과 접착제가 찢어지는 것을 두 개의 눈으로 따로따로 지켜보며, 가볍고 정확한 파괴 예측을 가능하게 한 혁신적인 시뮬레이션 기술입니다."

기술 요약

논문 요약: Puck 파괴 이론 기반의 섬유 강화 복합재 적층판용 다상계장 (Multi-Phase-Field) 모델

1. 연구 배경 및 문제 제기 (Problem)

• 복합재료의 복잡성: 섬유 강화 복합재 (FRC) 는 항공우주 및 자동차 등 경량화 분야에서 널리 사용되지만, 복잡한 파괴 메커니즘 (섬유 파단, 기지재 균열, 층간 박리 등) 으로 인해 파괴 예측이 어렵습니다.
• 기존 수치 해석의 한계:
  • 기존 유한요소법 (FEM) 기반의 손상 역학 (CDM) 모델은 균열의 전파를 모사하는 데 한계가 있으며, 균열 경로가 미리 정의되어야 하는 경우가 많습니다.
  • 기존 위상계장 (Phase-Field, PF) 방법들은 주로 단일 적층판이나 단순한 교차 적층 (Cross-ply) 구조에 국한되어 있었으며, 다양한 적층 순서 (Ply sequence) 를 가진 적층판의 면내 (In-plane) 균열 전파를 정확히 예측하는 데는 미흡했습니다.
  • 특히, 인장 및 압축 하중 하에서의 파괴 개시 (Crack nucleation) 를 설명하기 위해 Puck 과 같은 현상론적 파괴 기준을 PF 모델에 통합한 연구는 부족했습니다.
• 주요 목표: 다양한 적층 순서를 가진 FRC 적층판의 면내 손상을 예측하기 위해, Puck 파괴 이론을 기반으로 한 2 차원 다상계장 (Multi-Phase-Field) 모델을 개발하고, 이를 통해 섬유 파단과 기지재 (Inter-fiber) 파괴를 독립적으로 트리거하고 상호작용을 분석하는 것입니다.

2. 제안된 방법론 (Methodology)

가. 다상계장 모델 (Multi-Phase-Field Formulation)

• 독립된 위상계장 변수: 섬유 파괴 ($d_f$) 와 기지재 (Inter-fiber) 파괴 ($d_m$) 를 각각 독립적인 위상계장 변수로 정의하여 두 파괴 모드의 상호작용을 모사합니다.
• 에너지 분해: 시스템의 총 에너지를 섬유 관련 에너지 ($\Psi_f$) 와 기지재 관련 에너지 ($\Psi_m$) 로 분해하며, 각각에 대해 별도의 특성 길이 척도 (Characteristic length scale, $\ell_f, \ell_m$) 와 구조 텐서 (Structural tensor) 를 도입하여 파괴 모드를 정밀하게 기술합니다.
• 열역학적 일관성: 클라우지우스 - 더헴 (Clausius-Duhem) 부등식을 만족하도록 변분 형식 (Variational formulation) 을 유도하며, 파괴의 비가역성을 보장하기 위해 KKT 조건을 적용합니다.

나. Puck 파괴 이론의 통합

• 파괴 개시 기준: 위상계장의 구동력 (Driving force) 을 결정하기 위해 Puck 파괴 이론을 도입합니다.
  • 섬유 파괴: 인장 ($F^T_f$) 과 압축 ($F^C_f$) 모드를 구분합니다.
  • 기지재 파괴: 전단 및 인장/압축 응력 상태에 따라 모드 A(횡방향 인장), 모드 B(횡방향 압축 + 전단), 모드 C(우세한 횡방향 압축) 로 구분합니다.
• 구동력 (Driving Force): Puck 기준을 만족할 때의 에너지 노출 계수를 기반으로 위상계장의 진화를 유도하는 구동력 ($H_f, H_m$) 을 정의합니다.

다. 메시 오버레이 기법 (Mesh Overlay Method)

• 2 차원 모델의 효율성: 3 차원 모델의 높은 계산 비용을 줄이기 위해 2 차원 모델을 사용합니다.
• 적층판 구현: 각 적층 (Ply) 을 별도의 메시 (Dummy mesh) 로 정의하고, 모든 적층이 동일한 노드 (Nodes) 를 공유하도록 메시 오버레이 기법을 적용합니다.
  • 결합: 모든 적층의 변위는 강하게 결합 (Coupled) 되지만, 응력과 손상은 각 적층의 섬유 방향 ($\theta$) 에 따라 독립적으로 계산됩니다.
  • 제한점: 이 방법은 층간 박리 (Delamination) 를 명시적으로 모사하지는 못하지만, 면내 손상 전파를 효율적으로 계산할 수 있습니다.

라. 수치 구현

• Abaqus 의 사용자 정의 요소 (UEL) 를 통해 구현되었으며, 뉴턴 - 랩슨 (Newton-Raphson) 반복법을 사용하여 비선형 시스템을 해결합니다.

3. 주요 결과 (Results)

논문은 4 가지 벤치마크 예시를 통해 제안된 모델의 예측 능력을 검증했습니다.

1. 커폰 시험 (Coupon Tests):

   • 다양한 적층 순서 ([(45/-45)8]s, [(60/-60)8]s 등) 에 대한 인장 및 압축 시험을 수행했습니다.
   • 결과: 실험 결과와 힘 - 변위 곡선 및 최대 하중을 정량적으로 잘 일치시켰습니다. 특히 Puck 기준에 따른 파괴 개시 후 균열 전파가 즉시 발생하여 최대 하중이 특성 길이 척도에 민감하지 않음을 확인했습니다.

2. 개구부 인장 (Open Hole Tension, OHT):

   • [0/90]s 적층: 하중 방향과 수직인 섬유 방향에서 기지재 손상이 먼저 발생하고, 이후 섬유 파단이 일어나는 과정을 정확히 재현했습니다.
   • [45/-45]s 적층: 비대칭적인 균열 전파 (Crack shielding) 현상을 잘 포착하여 실험적 균열 패턴과 정성적으로 높은 일치도를 보였습니다.

3. 압축 시험 (Compact Tension, CT):

   • 블록 적층 (Blocked ply) 시편을 대상으로 실험 결과와 비교했습니다.
   • 결과: 힘 - 변위 곡선의 초기 선형 구간과 파괴 패턴에서 실험과 잘 일치했으나, 실험에서 관찰된 층간 박리로 인한 비선형성은 모델에서 완전히 재현되지 않았습니다 (모델의 2 차원적 한계).

4. 이중 노치 인장 (Double-Notched Tension):

   • 교차 적층 (Cross-ply) 및 등방성 (Quasi-isotropic) 적층판을 대상으로 하였습니다.
   • 결과: 섬유 지배적 파괴와 기지재 지배적 파괴 간의 경쟁 관계와 상호작용을 성공적으로 시뮬레이션했습니다. 하중 증가에 따른 하위 임계 손상 (Sub-critical damage) 의 발달과 균열 전파 경로를 실험 결과와 정성적으로 매우 잘 매칭했습니다.

4. 주요 기여도 (Key Contributions)

1. Puck 이론과 위상계장의 통합: 복합재료의 파괴 개시 (Strength) 와 전파를 동시에 설명하기 위해 Puck 파괴 이론을 위상계장 프레임워크에 성공적으로 통합했습니다.
2. 다중 적층 순서 처리: 메시 오버레이 기법을 통해 다양한 적층 순서 (Ply sequence) 를 가진 적층판의 면내 손상을 효율적으로 모델링할 수 있는 2 차원 프레임워크를 제시했습니다.
3. 독립적 파괴 모드 모델링: 섬유 파단과 기지재 파괴를 별도의 위상계장 변수와 구조 텐서로 분리하여, 두 파괴 모드 간의 복잡한 상호작용을 정밀하게 포착했습니다.
4. 검증된 예측 능력: 다양한 하중 조건 (인장/압축) 과 기하학적 구조 (개구부, 노치) 에서 실험 결과와 정성적/정량적으로 높은 일치도를 보이는 것을 입증했습니다.

5. 의의 및 결론 (Significance)

• 계산 효율성: 3 차원 모델에 비해 계산 비용을 크게 절감하면서도 복잡한 적층판의 파괴 거동을 예측할 수 있는 효율적인 도구를 제공합니다.
• 설계 지원: 복합재료 구조물의 손상 허용 설계 (Damage-tolerant design) 에 있어 파괴 메커니즘을 이해하고 최적화하는 데 유용한 통찰을 제공합니다.
• 한계 및 향후 과제: 현재 모델은 층간 박리 (Delamination) 를 고려하지 않으며, 소성 변형이나 소성 경화/연화 거동을 포함하지 않습니다. 또한, 적층 순서에 따른 층간 효과는 완전히 포착되지 않을 수 있습니다.
• 데이터 공개: 논문은 모든 예제에 대한 코드와 데이터를 공개하여 연구의 투명성을 높이고, 다른 연구자들이 모델을 확장할 수 있도록 지원하고 있습니다.

이 연구는 복합재료의 파괴 역학을 모델링하는 데 있어 위상계장 방법의 적용 범위를 확장하고, Puck 이론과의 결합을 통해 실제 공학적 문제 해결에 기여한다는 점에서 중요한 의의를 가집니다.