비선형 구조해석

구조해석 고체역학 마지막 수정: 2026-07-03 14:22:08

비선형 구조해석 Nonlinear Structural Analysis
분야	고체역학 × 수치해석
비선형의 종류	기하 / 재료 / 접촉
핵심 해법	증분-반복 (뉴턴-랩슨법)
수렴 무기	접선강성행렬, 호장법(arc-length)
주 기법	유한요소법

1. 개요[편집]

선형 해석은 한 번에 끝난다. 비선형 해석은 수렴할 때까지 기도해야 끝난다.

비선형 구조해석(nonlinear structural analysis)은 하중과 응답 사이의 관계가 비례하지 않는 — 즉 강성이 상수가 아닌 — 구조 문제를 다루는 해석이다. 선형 정적 해석이 $\mathbf{K}\mathbf{u} = \mathbf{f}$ 라는 한 방의 연립방정식으로 끝나는 반면, 비선형 해석에서는 강성 자체가 변위에 의존해서 $\mathbf{K}(\mathbf{u})\,\mathbf{u} = \mathbf{f}$ 가 되고, 이걸 풀려면 하중을 조금씩 나눠 걸며 매번 반복 계산을 돌려야 한다.¹

현실의 구조 대부분은 사실 비선형이다. 고무는 크게 늘어나고, 금속은 어느 지점부터 소성 변형하며, 부품끼리는 접촉하며 미끄러진다. 선형 해석은 이 모든 것을 “작은 변형, 착한 재료, 접촉 없음”이라고 가정한 편의점 버전일 뿐이다. 정밀함이 필요한 순간, 엔지니어는 결국 비선형의 늪으로 걸어 들어간다.

2. 세 가지 비선형[편집]

비선형성의 근원은 크게 세 가지로 분류되며, 실무에서는 이들이 종종 한꺼번에 몰려온다.

2.1. 기하 비선형 (대변형)[편집]

변형이 커서 원래 형상 기준의 선형 가정이 깨지는 경우다. 낚싯대가 활처럼 휘거나, 얇은 금속판이 크게 접히거나, 케이블이 늘어지는 상황. 변형이 커지면 하중의 방향과 작용점, 그리고 구조의 기하가 계속 바뀌기 때문에 평형 방정식을 변형된 형상에서 다시 세워야 한다. 좌굴 후 거동(좌굴 참조)도 대표적인 기하 비선형 문제다. 응력-변형률은 여전히 선형이어도, 기하 하나만으로 비선형이 된다.

2.2. 재료 비선형 (소성·초탄성)[편집]

재료의 응력-변형률 관계 자체가 직선이 아닌 경우다.

소성(plasticity) — 금속이 항복점을 넘으면 하중을 제거해도 변형이 남는다. 소성 유동 법칙과 경화 규칙(hardening rule)으로 이 비가역 거동을 모델링한다. 프레스 성형, 볼트 체결, 충돌 해석의 필수 요소.
초탄성(hyperelasticity) — 고무·엘라스토머는 수백 %까지 탄성적으로 늘어나며 응력-변형률이 극도로 비선형이다. Mooney-Rivlin, Ogden 같은 변형 에너지 밀도 함수로 기술한다. 타이어, 씰(seal), 방진 마운트 해석의 주역.

2.3. 접촉 비선형[편집]

두 물체가 닿았다 떨어졌다 하는 경계 조건 자체가 해에 의존하는 경우다. 어디가 닿을지, 얼마나 세게 눌릴지, 미끄러질지 붙어 있을지는 풀어보기 전에는 모른다. 접촉 상태가 매 반복마다 바뀌면서 강성이 불연속적으로 튀기 때문에, 세 비선형 중에서도 수렴이 가장 악명 높다.² 기어 물림, 조립체 간섭, 스냅핏 체결이 전형적인 접촉 문제다.

3. 증분-반복 해법[편집]

비선형 방정식 $\mathbf{K}(\mathbf{u})\,\mathbf{u} = \mathbf{f}$ 는 닫힌 형태로 풀 수 없으므로, 전체 하중을 여러 증분(increment/step)으로 쪼갠 뒤 각 증분에서 반복(iteration)으로 평형을 찾는다. 반복의 표준 도구가 뉴턴-랩슨법이다.

각 반복에서 우리가 없애려는 것은 내력과 외력의 차이, 즉 잔차(residual) $\mathbf{r}(\mathbf{u}) = \mathbf{f}_{\text{int}}(\mathbf{u}) - \mathbf{f}_{\text{ext}}$ 이다. 뉴턴-랩슨법은 잔차를 현재 변위 주변에서 선형화하여 변위 보정량 $\Delta\mathbf{u}$ 를 구한다.

\mathbf{K}_T(\mathbf{u}_i)\,\Delta\mathbf{u} = -\,\mathbf{r}(\mathbf{u}_i), \qquad \mathbf{u}_{i+1} = \mathbf{u}_i + \Delta\mathbf{u}

여기서 $\mathbf{K}_T$ 는 접선강성행렬(tangent stiffness matrix)로, 잔차를 변위로 미분한 것이다.

\mathbf{K}_T = \frac{\partial \mathbf{r}}{\partial \mathbf{u}} = \frac{\partial \mathbf{f}_{\text{int}}}{\partial \mathbf{u}}

기하학적으로 접선강성은 하중-변위 곡선의 현재 기울기다. 잔차가 충분히 작아질 때까지 이 과정을 반복하면 해당 증분이 수렴하고, 다음 증분으로 넘어간다. 뉴턴-랩슨법은 해 근처에서 2차 수렴(quadratic convergence)을 보여 매우 빠르지만, 매 반복마다 접선강성행렬을 새로 만들고 분해해야 해서 비싸다. 이 비용을 아끼려고 접선강성을 몇 번 재사용하는 수정 뉴턴법(modified Newton)도 널리 쓰인다.

4. 수렴 문제와 호장법[편집]

비선형 해석의 절반은 물리이고 절반은 수렴과의 전쟁이다. 잔차가 얌전히 줄지 않고 발산하거나 진동하는 일이 다반사다. 흔한 처방은 하중 증분을 잘게 쪼개고(자동 시간 증분), 접촉을 부드럽게 켜고, 초기 결함으로 대칭 분기를 깨는 것 등이다.

특히 골치 아픈 것이 하중-변위 곡선에 극점(limit point) 이 나타나는 경우다. 좌굴 스냅스루(snap-through)처럼 하중이 증가하다 감소로 꺾이는 지점에서는, 하중을 계속 증가시키는 일반 하중 제어(load control)가 물리적으로 존재하지 않는 해를 찾다가 발산한다.

이때 등장하는 구원투수가 호장법(arc-length method, Riks/Crisfield 법)이다. 하중과 변위를 각각 독립으로 제어하는 대신, 하중-변위 공간에서 “한 걸음의 호(弧) 길이”를 제약 조건으로 걸어 곡선을 따라간다. 덕분에 하중이 줄어드는 구간, 심지어 변위가 되돌아가는 구간까지 곡선 전체를 안정적으로 추적할 수 있다. 후좌굴 거동(좌굴)이나 스냅스루가 있는 쉘 구조를 풀 때 사실상 표준 기법이다.

5. 실무에서의 비선형 해석[편집]

자동차 — 충돌 안전 해석은 대변형·소성·접촉이 밀리초 안에 한꺼번에 일어나는 비선형의 종합선물세트다. 이 경우엔 접선강성 분해 대신 외연적(explicit) 시간 적분을 쓰는 것이 보통이다.
성형·제조 — 프레스 성형, 단조, 판재 굽힘은 전형적인 대변형+소성 문제로, 스프링백(제거 후 되돌아옴) 예측이 핵심 관심사다.
씰·개스킷·타이어 — 초탄성 재료와 접촉이 결합된 고전적 난제. 수렴만 시켜도 절반은 성공이라는 말이 나온다.

비선형 해석의 결과는 선형보다 훨씬 현실적이지만, 그 대가로 해석 시간은 수십~수백 배로 뛰고 수렴 설정에 노하우가 필요하다. 그래서 실무 격언은 이렇게 말한다. “선형으로 될 일을 비선형으로 풀지 마라. 하지만 비선형이 필요한 걸 선형으로 우기지도 마라.”³

6. 관련 문서[편집]

그래서 비선형 해석 로그 파일에는 “increment”, “iteration”, “residual” 세 단어가 수만 번 반복된다. 이 로그를 스크롤하며 잔차가 줄어드는지 노려보는 것이 비선형 해석 엔지니어의 주요 업무이자 취미이자 스트레스의 근원이다. ↩
접촉이 안 붙으면 강성 행렬이 특이(singular)해져 “0으로 나누기” 사고가 나고, 갑자기 붙으면 강성이 폭증해 해가 튄다. 그래서 페널티(penalty)나 라그랑주 승수(Lagrange multiplier)로 접촉을 부드럽게 다루는데, 페널티 강성 하나 잘못 잡으면 밤새 안 수렴한다. ↩
전자를 어기면 회사 클러스터를 일주일 태우고도 선형 해석이면 5분에 끝났을 답을 얻는다. 후자를 어기면 5분 만에 답을 얻지만 그 답이 틀렸다. 어느 쪽이 더 나쁜지는 마감이 언제냐에 달렸다. ↩