파괴역학

시뮬레이션 구조해석 고체역학 마지막 수정: 2026-07-09 04:11:27

파괴역학 Fracture Mechanics
분야	고체역학 × 재료공학 × 수치해석
핵심 파라미터	응력확대계수 K, 에너지 방출률 G, J-적분
두 갈래	LEFM(선형탄성), EPFM(탄소성)
수치 기법	유한요소법, XFEM, 코히시브 존

1. 개요[편집]

파괴역학(Fracture Mechanics)은 재료 내부에 이미 존재하는 균열(crack)이 언제, 어떻게 진전(propagation)하여 구조물을 파괴에 이르게 하는지를 정량적으로 다루는 학문이다. 전통적인 응력과 변형률 기반 설계가 “최대 응력이 항복강도를 넘지 않으면 안전하다”는 관점이라면, 파괴역학은 “결함은 어디에나 있다”는 전제에서 출발한다. 아무리 잘 만든 부재라도 미세한 균열이나 개재물(inclusion)은 반드시 존재하며, 그 결함 선단(crack tip)에서 응력이 이론적으로 무한대로 발산한다는 불편한 진실이 이 학문의 출발점이다.¹

이 관점의 실용적 무게는 결코 가볍지 않다. 타이타닉의 취성 파괴, 제2차 세계대전 리버티선의 선체 균열, 코멧 여객기의 창문 모서리에서 시작된 피로 파괴 — 굵직한 참사의 배후에는 대부분 파괴역학으로 설명되는 메커니즘이 있었다. 그래서 항공, 원자력, 압력용기 업계에서는 “손상 허용 설계(damage tolerance design)“라는 이름으로 이 학문을 법규 수준으로 강제한다.

2. 균열 선단 응력장과 응력확대계수[편집]

선형탄성 이론에서 균열 선단 주변의 응력장은 선단으로부터의 거리 $r$ 이 0으로 갈 때 $1/\sqrt{r}$ 속도로 발산한다. 이 특이성의 세기를 나타내는 단 하나의 파라미터가 응력확대계수(Stress Intensity Factor) $K$ 이다. 균열 선단 근방의 응력은 다음과 같이 쓸 수 있다.

\sigma_{ij}(r,\theta) = \frac{K}{\sqrt{2\pi r}}\, f_{ij}(\theta) + \cdots

여기서 $f_{ij}(\theta)$ 는 각도에만 의존하는 보편 함수다. 하중이 걸리는 방식에 따라 균열은 세 가지 모드로 나뉜다. 인장으로 벌어지는 모드 I(opening), 면내 전단인 모드 II(sliding), 면외 전단인 모드 III(tearing). 실무에서 압도적으로 중요한 것은 모드 I이며, 대표적으로 무한 평판 중앙에 길이 $2a$ 의 균열이 있고 원거리 응력 $\sigma$ 가 걸릴 때

K_{\mathrm{I}} = \sigma \sqrt{\pi a}

로 표현된다. 재료가 견딜 수 있는 $K$ 의 한계값을 파괴인성(fracture toughness) $K_{\mathrm{IC}}$ 라 하며, $K_{\mathrm{I}} \geq K_{\mathrm{IC}}$ 가 되는 순간 균열은 불안정하게 달려나간다.² $K_{\mathrm{IC}}$ 는 응력이나 형상이 아니라 재료 고유의 물성이라는 점이 핵심 — 그래서 이 값 하나로 서로 다른 형상의 부재를 비교할 수 있다.

3. 에너지 방출률과 그리피스[편집]

역사적으로 파괴역학의 진짜 시작점은 응력이 아니라 에너지였다. 1920년 A. A. 그리피스(Griffith)는 균열이 한 단위 면적만큼 진전할 때 방출되는 탄성 에너지가 새 표면을 만드는 데 드는 에너지보다 크면 균열이 자발적으로 성장한다는 에너지 균형 논리를 세웠다. 이 방출되는 에너지가 에너지 방출률(energy release rate) $G$ 이다. 선형탄성 조건에서 $G$ 와 $K$ 는 동전의 양면이며, 평면변형(plane strain) 상태에서

G = \frac{K_{\mathrm{I}}^2}{E'}, \qquad E' = \frac{E}{1-\nu^2}

의 관계로 정확히 연결된다. 파괴 조건은 $G \geq G_c$ 로도 동등하게 쓸 수 있다. 그리피스의 논리는 원래 유리 같은 순수 취성 재료에만 잘 맞았고, 금속처럼 선단에서 소성 변형이 크게 일어나는 재료로 확장한 사람이 어윈(Irwin)이다. 그가 소성 영역의 에너지 소산까지 $G_c$ 에 흡수시키면서 파괴역학은 비로소 공학 재료 전반에 쓸 수 있는 도구가 되었다.

4. LEFM vs EPFM — J-적분의 등장[편집]

선형탄성 파괴역학(LEFM)은 균열 선단의 소성 영역이 균열 길이나 부재 두께에 비해 충분히 작을 때(small-scale yielding)만 유효하다. 그런데 연성이 큰 금속이나 얇은 판재에서는 선단 소성역이 무시할 수 없을 만큼 커진다. 이때 $K$ 는 물리적 의미를 잃고, 탄소성 파괴역학(EPFM)으로 넘어가야 한다.

EPFM의 왕관은 라이스(Rice)가 1968년에 제안한 J-적분(J-integral)이다. J-적분은 균열 선단을 감싸는 임의의 경로를 따라 계산하는 선적분인데,

J = \int_{\Gamma} \left( W\, dy - \mathbf{t} \cdot \frac{\partial \mathbf{u}}{\partial x}\, ds \right)

가 그 정의다( $W$ 는 변형에너지 밀도, $\mathbf{t}$ 는 응력 벡터). J-적분의 마법 같은 성질은 경로 독립성이다 — 어떤 경로를 택하든 같은 값이 나온다. 덕분에 균열 선단에서 응력이 발산하는 지저분한 영역을 피해 멀찍이 떨어진 경로에서 계산할 수 있고, 이것이 수치해석에서 J-적분이 사랑받는 결정적 이유다. 선형탄성 극한에서는 $J = G$ 로 자연스럽게 환원되므로, J-적분은 $G$ 의 탄소성 일반화라고 봐도 무방하다.

5. FEM에서 균열 모델링[편집]

균열을 유한요소법으로 다루는 것은 오랫동안 골칫거리였다. 균열은 형상 자체가 불연속이라, 전통적 FEM에서는 균열면을 따라 격자를 정확히 정렬하고 선단에는 $1/\sqrt{r}$ 특이성을 표현하는 특이요소(quarter-point element)를 심어야 했다. 문제는 균열이 자라면 매 스텝마다 격자를 다시 짜야 한다는 것 — 3차원에서는 거의 재앙에 가깝다.³

이 고통을 덜어준 두 가지 현대적 접근이 있다.

XFEM(확장 유한요소법, eXtended FEM): 균열이 격자를 뚫고 지나가도록 허용한다. 단위 분할(partition of unity) 개념을 이용해 불연속 함수와 선단 특이 함수를 요소의 형상함수에 추가로 얹는 방식으로, 균열이 요소 경계와 무관하게 지나갈 수 있다. 매 스텝 재격자화의 지옥에서 해방되는 것이 최대 미덕.
코히시브 존 모델(Cohesive Zone Model, CZM): 균열 선단 앞에 가상의 “응집 영역”을 두고, 그 안에서 견인력(traction)과 분리 변위(separation) 사이의 구성 관계로 균열 진전을 기술한다. 응력 특이성을 원천적으로 제거하고 파괴를 에너지 소산 문제로 바꾸므로, 복합재 해석의 층간 박리(delamination)나 접착 조인트 해석의 국룰로 자리 잡았다.

여기에 최근에는 균열을 날카로운 불연속이 아니라 매끄러운 손상장(damage field)으로 번지게 표현하는 위상장 파괴 모델(phase-field fracture)도 각광받는다. 균열이 알아서 갈라지고(branching) 합쳐지는 복잡한 경로를 재격자화 없이 잡아낸다는 점에서 매력적이지만, 그 대가로 격자를 잔인하게 갈아 넣어야 한다는 반전이 있다.

6. 관련 문서[편집]

“응력이 무한대”라는 말에 물리학도는 경기를 일으키겠지만, 이것은 연속체 선형탄성 이론의 이상화가 만든 수학적 산물이다. 실제 재료에서는 선단에서 소성 변형이나 원자 결합 파괴가 일어나 응력을 유한하게 잘라낸다. 파괴역학은 그 무한대를 직접 다루는 대신, 발산의 “세기”인 $K$ 를 물리량으로 삼는 영리한 우회로다. ↩
여기서 “불안정하게”가 무섭다. 유리컵에 금이 가는 순간을 떠올려보라. $K_{\mathrm{IC}}$ 를 넘긴 취성 파괴는 음속에 가까운 속도로 전파되며, 엔지니어가 “어어?” 하는 순간 이미 끝나 있다. ↩
3차원 균열 전선(crack front)은 곡선이고, 진전하면서 형상이 바뀐다. 이걸 격자로 따라가겠다고 마음먹으면 메시 생성 코드와 영혼의 싸움을 벌이게 된다. XFEM이 논문 인용 수 폭발을 일으킨 데는 다 이유가 있다. ↩