초탄성 재료 모델

1. 개요[편집]

초탄성 재료 모델(hyperelastic material model)은 응력이 변형률 에너지 밀도 함수 $W$ 의 도함수로 정의되는 재료 모델로, 고무·엘라스토머·연조직처럼 크게 늘어났다가 완전히 되돌아오는 비선형 탄성 거동을 기술한다. 응력-변형률 곡선이 직선인 응력과 변형률의 선형 탄성(후크의 법칙)과 달리, 수백 %의 변형에서도 에너지를 저장하고 완전히 복원하는 재료를 다룬다.

핵심은 경로 독립성이다. 초탄성체는 변형 경로가 어떻든 현재 변형 상태만으로 저장된 에너지가 결정되며, 하중을 제거하면 그 에너지를 100% 돌려준다. 에너지가 상태 함수라는 이 성질이 초탄성을 소성과 근본적으로 가른다 — 소성은 변형 경로에 따라 에너지가 소산되고 영구 변형이 남는다.¹ 따라서 초탄성은 비선형 구조해석의 대표적 재료 비선형 문제이면서도, 여전히 탄성(가역)이라는 독특한 자리에 있다.

2. 변형률 에너지 밀도 함수[편집]

초탄성 모델의 심장은 단위 부피당 저장 에너지, 즉 변형률 에너지 밀도 함수 $W$ 다. 모든 응력이 여기서 미분으로 나온다. 유한변형을 다루므로 변형은 변형 구배 텐서 $\mathbf{F} = \partial \mathbf{x}/\partial \mathbf{X}$ 로 기술하고, 강체 회전을 걸러내기 위해 오른쪽 코시-그린 텐서 $\mathbf{C} = \mathbf{F}^T\mathbf{F}$ 를 쓴다.

등방성 재료라면 $W$ 는 $\mathbf{C}$ 의 세 불변량 $I_1, I_2, I_3$ 의 함수로 쓸 수 있다.

$$I_1 = \mathrm{tr}(\mathbf{C}), \quad I_2 = \tfrac{1}{2}\!\left[(\mathrm{tr}\,\mathbf{C})^2 - \mathrm{tr}(\mathbf{C}^2)\right], \quad I_3 = \det(\mathbf{C}) = J^2$$

여기서 $J = \det\mathbf{F}$ 는 부피 변화율이다. 제2 피올라-키르히호프 응력은 $W$ 의 도함수로 정의된다.

$$\mathbf{S} = 2\,\frac{\partial W}{\partial \mathbf{C}}$$

이 한 줄이 초탄성의 전부다. 서로 다른 모델이란 결국 $W$ 를 어떤 함수 형태로 가정하느냐의 차이일 뿐이다. 고무는 사실상 비압축성($J \approx 1$)이므로, 실무에서는 $W$ 를 형상을 바꾸는 편차(deviatoric) 부분과 부피를 바꾸는 체적(volumetric) 부분으로 쪼개, 불변량을 $\bar{I}_1 = J^{-2/3} I_1$ 처럼 부피 보정해 다룬다.

3. 대표 모델들[편집]

$W$ 의 함수 형태에 따라 이름난 모델이 여럿이다. 무엇을 고르냐는 변형 범위와 확보한 실험 데이터가 결정한다.

• Neo-Hookean — 가장 단순. $W = C_1(\bar{I}_1 - 3)$. 상수 하나($C_1 = \mu/2$)로 끝나며, 고분자 사슬의 통계역학에서 유도된 물리적 뿌리가 있다. 중간 변형(~50%)까지는 잘 맞지만 큰 변형에서 실측 곡선의 상승을 못 따라간다.
• Mooney-Rivlin — $W = C_1(\bar{I}_1 - 3) + C_2(\bar{I}_2 - 3)$. $I_2$ 항을 더해 유연성을 키웠다. 2계수 형태가 고무 해석의 오랜 국룰이며, 중간~큰 변형을 무난히 잡는다.
• Ogden — 불변량 대신 주신장(principal stretch) $\lambda_i$ 를 직접 쓴다.

$$W = \sum_{p=1}^{N} \frac{\mu_p}{\alpha_p}\left(\lambda_1^{\alpha_p} + \lambda_2^{\alpha_p} + \lambda_3^{\alpha_p} - 3\right)$$

항을 늘리면($N=3$ 정도) 300% 이상의 극단적 변형까지 실측을 정밀하게 따라가, 타이어·씰 같은 고성능 해석의 표준으로 통한다. 대신 계수가 많아 실험 데이터 없이 쓰면 곡선이 물리적으로 폭주하기 쉽다.²

• Yeoh — $\bar{I}_1$ 의 3차 다항식. $I_2$ 데이터가 부실할 때(단축 인장만 있을 때) 넓은 변형 범위를 그럭저럭 커버해 실무에서 인기.

4. 어디에 쓰나[편집]

초탄성이 필수인 자리는 “많이 늘어나고 완전히 되돌아오는” 모든 것이다.

• 엘라스토머 부품 — 타이어, O-링/씰, 방진 마운트, 벨트. 자동차·기계 산업의 고무 부품 해석은 사실상 초탄성 없이는 성립하지 않는다.
• 연조직(soft tissue) — 혈관, 피부, 근육, 연골. 생체역학에서 조직은 대개 초탄성으로 모델링되며, 섬유 방향성을 반영한 이방성 초탄성(Holzapfel 모델 등)으로 확장된다.
• 폼·젤·소프트 로봇 — 큰 변형이 기능 그 자체인 물질들. 게임·그래픽스의 천 시뮬레이션이나 소프트바디도 성긴 의미의 초탄성 에너지를 최소화하는 방식으로 굴러간다.

5. 수치해석에서의 골칫거리[편집]

초탄성 해석은 유한요소법으로 풀지만, 선형 탄성과는 차원이 다른 난이도를 자랑한다.

• 비압축성과 체적 잠김. $J \approx 1$ 이라는 거의 완전한 비압축 조건은 표준 요소를 **체적 잠김(volumetric locking)**으로 폭주시킨다. 응력과 변형률에서 $\nu \to 0.5$ 일 때 체적탄성계수가 발산하던 그 문제의 극단판이다. 압력을 독립 미지수로 두는 **혼합정식화(mixed u-P)**나 감차적분으로 우회한다.
• 강한 비선형과 수렴. $W$ 가 변형에 대해 비선형이라 뉴턴-랩슨법 반복이 필요하고, 하중을 한 번에 걸면 발산한다. 하중을 잘게 쪼개는 증분법이 필수이며, 좌굴·스냅을 동반하면 호장법까지 동원해야 한다.
• 재료 계수 피팅. 모델을 골랐어도 계수는 실험에서 나온다. 단축 인장·이축 인장·순수 전단 시험 곡선에 최소자승법으로 $W$ 를 맞추는데, 한 종류 시험만으로 피팅하면 다른 변형 모드에서 엉뚱한 응력을 내는 함정이 유명하다. 여러 변형 모드를 동시에 피팅하는 게 안전하다.
• 안정성(Drucker 조건). 계수를 잘못 잡으면 특정 변형 구간에서 접선 강성이 음수가 되어 재료가 물리적으로 불안정해진다. 상용 솔버가 “material stability check failed” 경고를 뱉는 이유가 대개 이것이다.

결국 초탄성 해석의 절반은 재료 모델을 얼마나 잘 고르고 계수를 얼마나 정직하게 피팅했느냐로 결정된다. “일단 Neo-Hookean으로 돌려”는 첫 시도로는 훌륭하지만, 큰 변형에서 실측과 안 맞으면 모델부터 의심하는 게 순서다.³

6. 관련 문서[편집]

• 응력과 변형률 · 고체역학
• 비선형 구조해석 · 소성
• 유한요소법 · 강성행렬
• 뉴턴-랩슨법 · 호장법
• 천 시뮬레이션 · 최소자승법

7. Footnotes[편집]

1. “초탄성(hyper-elastic)“의 hyper는 “매우 탄성적”이라는 뜻이 아니라, 응력이 퍼텐셜 함수($W$)의 도함수로 나온다는 수학적 구조를 가리킨다. 보존력이 퍼텐셜의 그래디언트인 것과 똑같은 구조라, “그린 탄성(Green-elastic)“이라고도 부른다. 이름에 낚여 “엄청 잘 늘어나는 재료”로만 이해하면 절반만 아는 것. ↩

2. Ogden 지수 $\alpha_p$ 는 정수가 아니어도 되고 음수여도 된다는 게 이 모델의 유연함이자 저주다. 자유도가 높아 실측을 잘 따라가지만, 데이터 없이 계수를 손으로 넣으면 곡선이 어느 순간 하늘로 솟거나 음의 강성으로 뒤집힌다. 자유도는 언제나 양날의 검이다. ↩

3. 고무 해석 초심자가 반드시 겪는 통과의례가 있다. 단축 인장 데이터만으로 Mooney-Rivlin 계수를 피팅해 놓고, 정작 부품이 이축으로 부풀 때 응력이 두 배씩 틀리는 것을 발견하는 순간. 교훈은 하나 — 재료는 네가 시험한 변형 모드에서만 정직하다. ↩