탄소성 해석

시뮬레이션 구조해석 고체역학 수치해석 마지막 수정: 2026-07-19 04:12:40

1. 개요[편집]

탄소성 해석 Elasto-Plastic Analysis
분야	고체역학 × 비선형 구조해석
핵심 알고리즘	리턴 매핑 (탄성 예측 - 소성 보정)
수렴 열쇠	일관 접선 강성 (algorithmic tangent)
구조	전역 뉴턴 루프 × 국소 응력 적분 루프

탄소성 해석(elasto-plastic analysis)은 재료가 항복 전에는 탄성, 항복 후에는 소성 거동을 하는 이력 의존적 구성식을 유한요소 프레임에 얹어, 하중 이력을 따라 응력과 변형을 시간(또는 하중) 적분으로 추적하는 수치해석이다. 항복 조건·경화 법칙·유동법칙이 무엇인지는 소성 문서가 다루고, 이 문서는 그것을 코드가 실제로 어떻게 푸는가를 다룬다.

핵심 난점은 하나다. 응력과 변형률의 선형 관계에서는 변형률만 알면 응력이 즉시 나오지만, 탄소성에서는 같은 변형률에도 어떻게 도달했느냐에 따라 응력이 달라진다. 응력은 상태가 아니라 경로의 함수다. 그래서 탄소성 해석은 필연적으로 증분(incremental) 해석이 되고, 매 증분마다 적분점 하나하나에서 미분방정식을 적분하는 국소 문제가 전역 강성행렬 방정식 안에 중첩된다. 이 이중 구조를 이해하면 탄소성 해석의 90%는 이해한 셈이다.

2. 이중 루프 구조[편집]

상용 코드가 탄소성 문제를 푸는 뼈대는 어느 솔버든 거의 동일하다.

전역 루프 — 하중을 $N$ 개 증분으로 쪼개고, 각 증분에서 내력과 외력의 불평형 잔차 $\mathbf{R} = \mathbf{f}^{ext} - \mathbf{f}^{int}$ 를 뉴턴-랩슨법으로 0에 수렴시킨다.
국소 루프 — 전역 반복이 준 변형률 증분 $\Delta\boldsymbol{\varepsilon}$ 을 받아, 각 적분점마다 구성식을 적분해 갱신된 응력 $\boldsymbol{\sigma}_{n+1}$ 과 상태변수(누적 소성변형률, 후방응력)를 돌려준다. 동시에 접선 강성 $\partial\Delta\boldsymbol{\sigma}/\partial\Delta\boldsymbol{\varepsilon}$ 도 같이 돌려준다.

즉 전역 솔버는 “변형률 줄 테니 응력과 접선 내놔”라고 요구하는 소비자고, 재료 모델은 그 요구에 응답하는 서브루틴이다. Abaqus의 UMAT, LS-DYNA의 사용자 재료, ANSYS Mechanical의 USERMAT이 전부 이 인터페이스다. 재료 모델을 새로 만든다는 건 사실상 이 국소 루프를 짜는 일이다.¹

3. 방사 리턴 매핑[편집]

$J_2$ 소성(폰 미제스 + 등방경화)의 국소 적분은 놀랍도록 깔끔하게 닫힌 형태로 풀린다. 이것이 방사 리턴(radial return)이며, 탄소성 해석에서 가장 널리 구현된 알고리즘이다.

먼저 증분 전체가 탄성이라 가정한 시험 응력의 편차 성분 $\mathbf{s}^{tr}$ 을 계산한다. 이것이 항복면 밖이면 소성 보정이 필요하고, 그 크기는 다음 스칼라 방정식에서 나온다.

\|\mathbf{s}^{tr}\| - 2\mu\,\Delta\gamma - \sqrt{\tfrac{2}{3}}\,\sigma_Y(\bar{\varepsilon}^p_n + \sqrt{\tfrac{2}{3}}\Delta\gamma) = 0

여기서 $\mu$ 는 전단탄성계수, $\Delta\gamma$ 는 소성 승수 증분이다. 선형 등방경화( $\sigma_Y = \sigma_{Y0} + H\bar{\varepsilon}^p$ )라면 이 식은 1차식이 되어 반복 없이 바로 풀린다.

\Delta\gamma = \frac{\|\mathbf{s}^{tr}\| - \sqrt{\tfrac{2}{3}}\,\sigma_Y(\bar{\varepsilon}^p_n)}{2\mu + \tfrac{2}{3}H}

갱신된 편차응력은 시험 응력과 같은 방향이며 크기만 줄어든다. 즉 응력점이 항복면 중심을 향해 방사 방향으로 곧게 끌려 들어온다. “방사 리턴”이라는 이름이 여기서 왔다. 정수압 성분은 소성 비압축성 때문에 손대지 않는다. 이 방향 불변성 덕에 $J_2$ 소성은 국소 반복조차 필요 없는, 소성 알고리즘 중 가장 값싼 축에 속한다.²

비선형 경화(Voce·멱승법칙)나 이동경화가 섞이면 위 식이 비선형이 되어 적분점 내부에서 다시 뉴턴 반복을 돌려야 한다. 그래도 미지수는 스칼라 하나뿐이라 비용은 여전히 저렴하다. Drucker-Prager나 Hill48 이방성처럼 항복면이 방사 대칭이 아니면 이 우아함이 무너지고, 미지수 여러 개짜리 완전 암시적 리턴 매핑을 풀어야 한다.

4. 일관 접선 강성[편집]

탄소성 해석에서 초심자가 가장 자주 밟는 지뢰가 여기 있다. 리턴 매핑으로 응력은 정확히 구했는데 전역 뉴턴 반복이 한없이 기어가는 현상이다. 원인은 대개 접선 강성을 잘못 준 것이다.

교과서의 연속체 탄소성 접선(continuum tangent)은 $\Delta t \to 0$ 극한의 미분이다. 하지만 코드가 실제로 계산하는 응력은 유한한 증분에 대해 리턴 매핑이 뱉은 값이다. 그러므로 뉴턴법이 요구하는 정확한 야코비안은 그 알고리즘 자체를 미분한 것이어야 한다. 이것을 일관 접선 강성(consistent/algorithmic tangent)이라 부른다. $J_2$ 방사 리턴의 경우 잘 알려진 닫힌 형태가 있다.

\mathbb{C}^{alg} = \kappa\,\mathbf{1}\otimes\mathbf{1} + 2\mu\theta\left(\mathbb{I}^{dev}\right) - 2\mu\bar{\theta}\,\mathbf{n}\otimes\mathbf{n}

\theta = 1 - \frac{2\mu\,\Delta\gamma}{\|\mathbf{s}^{tr}\|},\qquad \bar{\theta} = \frac{1}{1 + H/(3\mu)} - (1-\theta)

연속체 접선과의 차이는 $\theta$ 항 하나뿐이지만, 그 하나가 뉴턴-랩슨법의 2차 수렴을 살리느냐 죽이느냐를 가른다. 일관 접선을 쓰면 증분당 반복이 3~5회로 끝나고, 연속체 접선을 쓰면 20회를 넘기다 결국 증분 절반 나기(cutback)에 들어간다. UMAT을 짜본 사람이 “응력은 맞는데 왜 안 수렴하죠”라고 묻는다면 답은 십중팔구 DDSDDE다.³

참고로 명시적 동해석은 강성행렬을 아예 안 만들므로 이 고민이 없다. 성형·충돌 해석이 명시적으로 도망가는 실질적 이유 중 하나다. 대신 시간 스텝이 안정 조건에 묶인다.

5. 유한변형과 요소 선택[편집]

금속 성형처럼 변형률이 수십 %를 넘어가면 미소변형 가정이 깨지고, 응력 갱신에 객관성(objectivity) 문제가 생긴다. 강체 회전만 시켜도 성분이 변하는 응력 텐서를 그대로 적분하면 회전만으로 가짜 응력이 생기기 때문이다. 그래서 Jaumann·Green-Naghdi 같은 객관 응력률을 쓰거나, 회전을 걷어낸 변형 구배 기반으로 정식화한 뒤 리턴 매핑을 그대로 재사용한다. 놀랍게도 대수 변형률(logarithmic strain) 공간에서 정식화하면 미소변형용 리턴 매핑 코드를 거의 그대로 쓸 수 있어, 초탄성 재료 모델의 유한변형 틀과 소성을 자연스럽게 결합할 수 있다.

요소 선택도 결과를 좌우한다. 소성 변형은 부피를 보존하므로, 완전적분 저차 요소는 체적 잠김(volumetric locking)에 걸려 실제보다 몇 배 뻣뻣한 답을 낸다. 대책은 감차적분(+모래시계 제어), B-bar/F-bar, 혼합 정식화 등이다. 소성 해석 결과가 실험보다 유난히 딱딱하게 나온다면 재료 계수를 의심하기 전에 요소 종류부터 보는 게 순서다.

한편 완전소성에 가까운 재료가 국부적으로 연화되면 해가 격자 크기에 의존하는 변형률 국부화(strain localization)가 나타난다. 격자를 조밀하게 할수록 파단 밴드가 얇아지며 흡수 에너지가 0으로 수렴하는 병리적 현상인데, 비국소 모델이나 손상역학 기반 정규화로 달래야 한다. 검증 및 확인에서 격자 수렴성 시험이 유독 안 되는 문제가 있다면 이걸 의심하자.

6. 실무 체크리스트[편집]

증분 크기. 한 증분에서 항복면을 크게 넘어가면 적분 오차가 누적된다. 응력 경로가 꺾이는 구간에서는 증분을 잘게.
경화 곡선 입력. 대부분의 코드는 진응력-진소성변형률을 요구한다. 공칭 곡선을 그대로 넣는 실수가 초심자 사고의 절반이다.
경화 데이터 범위. 입력 곡선의 마지막 점 너머는 코드가 평평하게(완전소성) 외삽한다. 변형률이 데이터 범위를 넘었는지 항상 확인.
잔여 응력 상태. 잔류응력이 중요한 문제(숏피닝, 용접 해석)는 제하(unloading) 단계까지 반드시 해석해야 의미가 있다.
극한하중. 구조가 붕괴 근처로 가면 하중 제어 뉴턴법이 발산한다. 변위 제어나 호장법으로 전환.

7. 관련 문서[편집]

UMAT은 구조해석 업계의 통과의례다. 인터페이스가 요구하는 건 딱 두 개 — 갱신된 응력 배열 STRESS와 접선 강성 행렬 DDSDDE. 그런데 이 두 줄짜리 요구사항 앞에서 대학원생 학기 하나가 증발한다. ↩
이 우아함은 폰 미제스 항복면이 편차 응력 공간에서 완벽한 구(球)라는 기하학적 사실에서 나온다. 구면으로의 최단 투영은 언제나 중심 방향이므로 방향이 보존된다. 트레스카였으면 모서리에서 방향이 정의 안 돼 난리가 났을 것이다. 수학적 편의가 또 한 번 표준을 결정한 사례. ↩
접선이 틀려도 답 자체는 맞는다는 점이 이 버그를 지독하게 만든다. 잔차가 0으로 갔다면 응력은 옳다. 틀린 접선은 “틀린 답”이 아니라 “느린 답”으로 나타나므로, 검증 단계에서 안 잡히고 대형 모델에서 계산 시간 폭발로 뒤늦게 드러난다. ↩