손상역학

구조해석 고체역학 마지막 수정: 2026-07-22 04:34:19

1. 개요[편집]

손상역학(Continuum Damage Mechanics, CDM)은 재료 내부에 흩어져 발생하는 미세균열·공동(void)의 누적을 하나의 연속체 장 변수 **손상 변수 $D$ **로 뭉뚱그려 표현하고, 이 $D$ 가 커질수록 재료의 강성과 강도가 점진적으로 떨어지는 과정을 기술하는 고체역학 이론이다. $D=0$ 이면 처녀(intact) 상태, $D=1$ 이면 완전 파단 상태로 정의하며, 그 사이의 연속적인 값으로 “재료가 얼마나 골병들었는가”를 정량화한다.

핵심 아이디어는 프랑스의 카샤노프(L. Kachanov)가 1958년에 크리프 파단을 설명하려고 도입한 것으로¹, 이후 라보트노프(Rabotnov), 르메트르(J. Lemaitre), 샤보슈(J.-L. Chaboche)가 열역학적으로 다듬어 오늘날의 골격을 완성했다. 한마디로 “균열을 하나하나 그리기 귀찮으니 재료가 삭은 정도를 필드로 바르자”는 발상이다.

2. 파괴역학과 무엇이 다른가[편집]

이 문서를 읽는 사람이 가장 헷갈려 하는 지점이 바로 파괴역학과의 차이다. 결론부터 말하면 다루는 대상의 스케일과 표현 방식이 정반대다.

파괴역학(Fracture Mechanics)은 **이미 존재하는 하나의 뾰족한 균열 선단(crack tip)**을 명시적으로 놓고, 그 주변의 특이(singular) 응력장을 응력확대계수 $K$ 나 J-적분 같은 파라미터로 정량화한다. 균열은 기하학적으로 뚜렷한 불연속(discontinuity)이며, 관심사는 “이 균열이 언제 어떻게 진전하느냐”다.
손상역학은 균열 선단을 그리지 않는다. 대신 재료 전체에 퍼진 미세손상을 $D(\mathbf{x})$ 라는 매끈한 연속 필드로 번져 바른다(smear). 불연속 대신 강성 저하라는 연속적 물리량으로 파손을 표현한다는 점이 결정적 차이다.

비유하자면 파괴역학은 “찢어진 종이의 찢김 끝을 확대경으로 보는 것”이고, 손상역학은 “종이 전체가 얼마나 삭았는지 색으로 칠하는 것”이다. 실제 해석에서는 CDM으로 손상이 국부적으로 집중되어 $D \to 1$ 에 도달하면, 그 지점을 매크로 균열의 발생·전파로 해석해 두 이론을 이어 붙이기도 한다.²

3. 유효응력과 변형률 등가 원리[편집]

CDM의 수학적 엔진은 유효응력(effective stress) 개념이다. 손상으로 실제 하중을 견디는 유효 단면적이 줄었다고 보고, 명목응력 $\sigma$ 를 남은 면적으로 나눠 진짜 재료가 느끼는 응력을 정의한다.

\tilde{\sigma} = \frac{\sigma}{1 - D}

여기에 르메트르의 변형률 등가 원리(strain equivalence)를 얹는다. “손상된 재료의 변형 거동은, 명목응력을 유효응력 $\tilde\sigma$ 로 바꾼 처녀 재료의 구성식과 똑같다”는 가정이다. 이 원리를 1차원 선형 탄성에 적용하면 손상된 재료의 유효 탄성계수가 곧바로 나온다.

\varepsilon = \frac{\tilde{\sigma}}{E} = \frac{\sigma}{(1-D)E}, \qquad E_{\text{eff}} = (1-D)\,E

즉 $D$ 는 강성이 얼마나 깎였는지를 직접 나타낸다. 실험적으로 $D$ 를 재는 고전적 방법도 여기서 나온다. 하중을 걸었다 풀며 반복 재하 시 탄성 계수가 처음 대비 얼마나 떨어졌는지를 측정하면 $D = 1 - E_{\text{eff}}/E$ 로 손상량을 역산할 수 있다.³ 스칼라 $D$ 는 등방 손상 가정이고, 방향에 따라 손상이 다르면 2계·4계 손상 텐서로 확장한다. 다만 텐서 손상은 파라미터 지옥이라 실무에서는 웬만하면 스칼라로 버틴다.

4. 소성과의 연성 (탄성-손상-소성)[편집]

금속의 연성 파괴는 소성 변형과 손상이 손을 잡고 진행된다. 그래서 CDM은 보통 소성 구성식에 손상을 결합한 탄소성-손상(elasto-plastic damage) 형태로 쓰인다. 열역학적으로는 손상을 하나의 내부 상태변수로 놓고, 소성 변형이 축적될수록 손상이 자라도록 진화식(evolution law)을 세운다. 르메트르 모델의 대표적인 손상 진화식은 다음과 같다.

\dot{D} = \frac{\dot{p}}{1 - D}\left(\frac{Y}{S}\right)^{s}

여기서 $p$ 는 누적 소성 변형률, $Y$ 는 손상을 밀어붙이는 열역학적 구동력(손상 에너지 방출률), $S$ 와 $s$ 는 재료 상수다. 소성이 흐를 때만( $\dot p > 0$ ) 손상이 자라고, 응력 삼축비(triaxiality)가 높을수록 $Y$ 가 커져 손상이 빨라진다. 이 구동력이 정수압 응력 상태에 민감하다는 점이 연성 파괴가 응력 삼축비에 크게 좌우되는 실험 사실과 잘 맞는다.

손상이 강성을 깎으면 응력이 떨어지고, 이 연화(softening) 거동이 비선형 구조해석에서 해를 불안정하게 만든다. 특히 국소 CDM은 그물망을 촘촘히 할수록 손상이 한 요소 폭으로 쪼그라드는 격자 의존성(mesh dependency) 문제를 일으켜, 이를 막으려고 비국소(nonlocal)·구배(gradient) 손상 모델이나 특성 길이 도입이 필요하다.⁴

5. 크리프·피로에서의 활용[편집]

CDM의 뿌리가 크리프 해석이었던 만큼, 고온 장시간 하중을 받는 부품의 크리프 파단 수명 예측에 카샤노프-라보트노프 손상 모델이 지금도 표준으로 쓰인다. 응력이 걸린 채 시간이 흐르면 $D$ 가 서서히 자라다가 파단 직전 급격히 폭주해 파단 시각을 준다.

피로 해석과의 연결도 중요하다. 반복 하중 한 사이클마다 손상이 조금씩 축적된다고 보는 연속 손상 피로 모델이 있는데, 르메트르류 CDM을 사이클 단위로 적분하면 유명한 마이너 선형 누적 법칙(Miner’s rule, $\sum n_i/N_i = 1$ )이 손상역학의 한 특수 케이스로 유도된다. 즉 “손상 $D$ 가 1에 도달하면 파단”이라는 CDM의 관점이 피로 수명 예측의 누적 손상 개념과 그대로 포개진다.⁵ 비선형 손상 누적을 쓰면 하중 순서 효과(load sequence effect)까지 잡아낼 수 있어, 단순 마이너 법칙보다 정교한 피로 예측이 가능하다.

6. 한계와 현실[편집]

파라미터 동정(identification)이 고역이다. $S$ , $s$ , 임계 손상 $D_c$ 같은 상수를 뽑으려면 반복 재하 시험 데이터가 필요한데, 재료마다 다시 갈아 넣어야 한다.
연화 구간의 격자 의존성 때문에 “일단 돌려” 정신으로 국소 모델을 쓰면 결과가 격자 따라 춤춘다. 비국소·구배 정규화는 필수에 가깝다.
그럼에도 CDM은 균열을 일일이 안 그리고도 파손 개시부터 진전까지 하나의 연속체 틀로 다룰 수 있어, 적층 제조 시뮬레이션이나 복합재 해석의 층간분리(delamination) 예측 등에서 꾸준히 쓰인다.

7. 관련 문서[편집]

카샤노프가 원래 쓴 변수는 손상 $D$ 가 아니라 “연속도(continuity)” $\psi = 1 - D$ 였다. 즉 원조는 재료가 “얼마나 멀쩡한가”를 셌고, 우리는 “얼마나 골병들었는가”를 센다. 반잔의 물을 반대로 본 셈. ↩
그래서 실무자들 사이에선 “CDM은 균열 그리기 싫은 사람이 쓰는 파괴역학”이라는 농담이 돈다. 물론 CDM 하는 사람 앞에서 하면 킹받아 한다. ↩
이 강성 저하 측정법은 이론은 우아한데 실험은 짜증난다. 재하-제하를 수십 번 반복하며 미세한 기울기 변화를 재야 해서, 노이즈에 손상값이 통째로 묻히기 일쑤다. ↩
국소 손상 모델의 격자 의존성은 유한요소 초심자를 절망에 빠뜨리는 국룰 함정이다. “격자를 촘촘히 했더니 강도가 0으로 수렴한다”는 괴담의 정체가 바로 이것. ↩
마이너 법칙이 CDM의 특수 케이스라는 사실을 알고 나면, 피로 실무에서 그렇게 애지중지하던 $\sum n_i/N_i = 1$ 이 사실은 손상역학의 가장 게으른 버전이었다는 걸 깨닫고 묘한 배신감을 느끼게 된다. ↩