XFEM

1. 개요

XFEM eXtended Finite Element Method
한국어	확장 유한요소법
핵심 아이디어	불연속을 요소망 재구성 없이 근사공간에 확충
이론 기반	단위 분할(Partition of Unity)
주 응용	파괴역학 — 균열 전파, 상 경계, 공극

XFEM(eXtended Finite Element Method, 확장 유한요소법)은 균열·재료 경계·공극 같은 불연속(discontinuity)을 요소망을 그 불연속에 맞춰 재구성하지 않고도 표현하는 유한요소법의 확장 기법이다. 핵심 슬로건은 하나다. “격자는 그대로 두고, 근사 함수 공간에 불연속을 얹는다.” 균열이 요소 한복판을 뚫고 지나가도 격자를 다시 짜지 않는다.

이게 왜 혁명이었는지는 파괴역학 문서가 잘 설명한다. 전통적 FEM으로 균열을 다루려면 균열면을 따라 격자를 정확히 정렬하고, 균열이 자랄 때마다 매 스텝 재격자화(remeshing)를 해야 했다. 3차원 곡선 균열 전선에서는 이게 “메시 생성 코드와 영혼의 싸움”이었다. XFEM은 균열을 격자에서 해방시켜 이 싸움을 통째로 없앴다. 1999년 벨리치코(Belytschko)와 블랙(Black)이 제안한 이래 파괴역학 수치해석의 판을 바꿨다.¹

2. 단위 분할과 확충

XFEM의 수학적 토대는 단위 분할(Partition of Unity, PoU)이다. 표준 FEM 형상함수 $N_i(\mathbf{x})$ 는 임의의 점에서 모두 더하면 1이 된다.

$$\sum_i N_i(\mathbf{x}) = 1$$

이 성질 덕분에, 어떤 특별한 함수 $\psi(\mathbf{x})$ 를 형상함수에 곱해 근사공간에 끼워 넣으면, 그 함수가 요소 안에서 정확히 재현된다. XFEM은 이 원리로 표준 FEM 근사에 불연속·특이 함수를 “확충(enrichment)“한다. 확충된 변위장은 다음 형태를 갖는다.

$$\mathbf{u}^h(\mathbf{x}) = \sum_{i} N_i(\mathbf{x})\,\mathbf{u}_i + \sum_{j} N_j(\mathbf{x})\,\psi(\mathbf{x})\,\mathbf{a}_j$$

앞항은 우리가 아는 평범한 FEM 보간이고, 뒷항이 확충항이다. $\mathbf{a}_j$ 는 확충된 절점에만 추가로 붙는 확충 자유도(enriched DOF)로, 불연속의 세기를 담는다. 격자는 그대로인데 자유도만 국소적으로 늘어나는 셈이다. 어떤 절점을 확충할지, 어떤 $\psi$ 를 쓸지가 XFEM 설계의 전부다.

3. 헤비사이드와 팁 확충 함수

균열을 표현하려면 두 종류의 확충 함수가 필요하다. 균열이 만드는 물리가 위치마다 다르기 때문이다.

• 헤비사이드 확충(Heaviside enrichment) — 균열면을 완전히 가로지르는 요소에는 변위가 점프한다. 균열의 위쪽과 아래쪽이 따로 움직이는 것이다. 이 불연속을 부호 함수(±1)로 표현한다.

$$H(\mathbf{x}) = \begin{cases} +1, & \text{균열면 위쪽} \\ -1, & \text{균열면 아래쪽} \end{cases}$$

• 팁 확충(crack-tip / branch enrichment) — 균열 선단(tip)이 들어 있는 요소에는 파괴역학의 그 유명한 $1/\sqrt{r}$ 특이성이 있다. 선단 근방 점근해에서 유도된 네 개의 함수로 이 특이장을 정확히 재현한다.

$$\{\psi_k(r,\theta)\} = \sqrt{r}\left\{ \sin\tfrac{\theta}{2},\ \cos\tfrac{\theta}{2},\ \sin\tfrac{\theta}{2}\sin\theta,\ \cos\tfrac{\theta}{2}\sin\theta \right\}$$

첫 항 $\sqrt{r}\sin(\theta/2)$ 가 균열면에서 변위 점프를 만들면서 동시에 특이 변형률을 담당한다. 이 팁 확충 덕분에 선단에 특이요소(quarter-point element)를 심지 않고도 응력확대계수 $K$ 를 정확히 뽑을 수 있다. 균열면을 가로지르는 요소는 헤비사이드로, 선단 요소는 팁 함수로 확충하는 이 조합이 표준 균열 XFEM이다.²

4. 레벨셋으로 균열 추적

균열이 격자와 무관하게 지나간다면, 그 균열이 지금 어디에 있는지를 격자 없이 기술해야 한다. 이때 쓰는 것이 레벨셋 방법(level set method)이다. 균열 기하를 두 개의 부호 거리 함수로 암시적으로 표현한다.

• $\phi(\mathbf{x})$ — 균열면까지의 부호 거리. $\phi = 0$ 이 균열면이다.
• $\psi(\mathbf{x})$ — 균열 선단까지의 거리. $\psi \le 0$ 이면서 $\phi = 0$ 인 점이 균열이 실제로 존재하는 영역이다.

각 절점에서 두 레벨셋 값의 부호만 보면, 그 요소가 균열에 완전히 잘렸는지(헤비사이드 확충 대상), 선단을 품고 있는지(팁 확충 대상)를 자동으로 판정할 수 있다. 균열이 자라면 레벨셋 함수만 갱신하면 되고, 격자는 손대지 않는다. 불연속의 기하를 격자에서 완전히 분리한 이 조합(XFEM + 레벨셋)이 균열 전파 해석을 우아하게 만든 핵심이다.

5. 응용과 한계

XFEM은 균열만의 도구가 아니다. PoU 확충이라는 발상 자체가 범용이라, 확충 함수만 바꾸면 온갖 불연속을 다룬다.

• 균열 전파 — 원조이자 최대 응용. 재격자화 없이 균열 성장 경로를 추적한다.
• 재료 경계·개재물 — 변위는 연속이되 변형률이 꺾이는 약한 불연속(weak discontinuity)은 기울기 점프 함수로 확충한다. 복합재·다결정 재료 해석에 쓴다.
• 공극·보이드, 상 경계 — 유동의 자유표면이나 이상(two-phase) 경계도 같은 틀로 다룬다.

물론 공짜는 아니다. XFEM의 그늘도 뚜렷하다.

• 수치 적분 — 확충된 요소 안에 불연속이 지나가면 표준 가우스 적분이 안 맞는다. 요소를 불연속 기준으로 잘게 나눠(subtriangulation) 적분해야 해서 구현이 번거롭다.
• 조건수 악화 — 균열이 요소 경계에 아슬아슬하게 붙으면 확충 자유도가 거의 종속이 되어 강성행렬의 조건수가 폭발한다. 안정화 기법이나 전처리기가 필요하다.
• 혼합 자유도 — 확충 절점과 표준 절점이 섞인 경계(blending element)에서 정확도가 떨어지는 문제가 오래 연구 주제였다.

그래서 균열 경로가 복잡하게 갈라지고 합쳐지는 문제에서는 파괴역학 문서가 언급한 위상장 파괴(phase-field) 모델이 대안으로 부상했다. XFEM이 “격자를 아끼고 자유도를 쓰는” 쪽이라면, 위상장은 “자유도를 아끼고 격자를 갈아 넣는” 반대 철학이다. 둘 다 재격자화 지옥을 피하려는 같은 목적지를 향한 다른 길이다.³

6. 관련 문서

• 파괴역학 · 피로 해석
• 유한요소법 · 강성행렬 · 조건수
• 전처리기 · 비선형 구조해석
• 응력과 변형률 · 소성
• 복합재 해석 · 구조해석
• 잔류응력 · Abaqus

7. Footnotes

1. XFEM 논문의 인용 수 폭발은 파괴역학 연구자들이 재격자화에서 얼마나 고통받고 있었는지의 역설적 증거다. “이제 균열마다 메시를 다시 안 짜도 된다”는 한 문장이 수천 편의 후속 논문을 낳았다. 편의는 인용을 부른다. ↩

2. 팁 확충 함수 네 개 중 첫 번째만 균열면에서 실제로 불연속이고 나머지 셋은 연속이다. 그런데도 넷을 다 쓰는 이유는, 그래야 선단 근방 응력장의 각도 의존성 $f_{ij}(\theta)$ 까지 정확히 재현되어 $K$ 값이 정밀해지기 때문이다. 하나로 줄이면 편해지지만 정확도를 잃는 전형적 트레이드오프다. ↩

3. “격자와 물리를 분리한다”는 XFEM의 정신은 사실 무요소법(meshfree)이나 등기하해석(IGA)과 같은 시대적 흐름의 일부다. 다들 “격자가 물리를 지배하는” 전통 FEM의 속박에서 벗어나려 몸부림친 결과물이다. XFEM은 그중 가장 실용적으로 살아남은 타협안이다. ↩