구조해석

1. 개요

구조해석 Structural Analysis
분야	고체역학 × 수치해석 × 기계·토목·항공 설계
표준 기법	유한요소법 (사실상 단일 체제)
지배 방정식	정적: Ku = f / 동적: Mü + Cu̇ + Ku = f(t)
기원	1956년 터너·클라우 등의 항공기 날개 강성 해석
해석 유형	선형 정적 · 모드 · 좌굴 · 비선형 · 피로 · 충돌
검증 기준	패치 테스트 · 메시 독립성 · ASME V&V 10
대표 코드	Abaqus · ANSYS Mechanical · NASTRAN · LS-DYNA

구조해석(structural analysis)은 고체 및 구조물이 하중을 받을 때 발생하는 응력(stress), 변형률(strain), 변위(displacement)를 수치적으로 예측하는 해석 분야이다. 다리가 무너지지 않을지, 부품이 부러지지 않을지, 흔들리다 공진하지 않을지를 시제품을 부수기 전에 컴퓨터로 미리 확인하는 것이 목적이다. 전산유체역학과 함께 CAE의 양대 기둥으로 꼽히며, 역사는 오히려 CFD보다 길다.

흔히 유한요소법의 출생신고로 인용되는 문서가 1956년 터너·클라우·마틴·톱의 항공기 날개 논문이고, 수학 쪽 뿌리는 그보다 앞선 1943년 쿠란트의 삼각형 분할 변분해법까지 올라간다. 즉 이 분야는 “항공기 날개가 얼마나 휘는지 알아야 했던 사람들”이 만든 수치해석이다.

2. 유한요소법이 표준인 이유

구조해석의 수치 기법은 사실상 유한요소법(FEM) 단일 체제다. 1950년대 항공기 구조 해석에서 출발한 태생부터가 고체역학이었고, 이후로도 이 분야와 함께 성장했다. 기술적인 이유도 명확하다.

• 약형식(weak form) 정식화가 복잡한 형상과 다양한 경계 조건을 자연스럽게 수용한다.
• 변위 기반 정식화에서 선형 탄성 문제의 강성행렬은 대칭 양정치(symmetric positive definite)가 되어, 촐레스키 분해처럼 효율적이고 안정적인 선형 솔버를 쓸 수 있다.
• 요소 차수를 높이는 것(p-refinement)만으로 정확도를 올릴 수 있어, 응력 집중부 해상도 확보에 유리하다.

유체 쪽에서 유한체적법이 보존성을 무기로 왕좌를 차지한 것과 대조적인 지형이다.

3. 해석 유형

• 선형 정적 해석 — 하중과 변형이 비례한다고 가정하는 가장 기본적인 해석. 대부분의 설계 검토는 여기서 시작한다.
• 모드 해석(고유치 해석) — 구조물의 고유진동수와 모드 형상을 구한다. 가진 주파수와 고유진동수가 겹치면 공진이므로, 회전기계나 차량 부품 설계에서 필수다.¹
• 좌굴 해석 — 압축 하중을 받는 세장 구조가 급격히 붕괴하는 임계 하중을 고유값 문제로 추정한다.
• 비선형 해석 — 대변형(기하 비선형), 소성·초탄성(재료 비선형), 접촉(경계 비선형)을 다룬다. 셋 다 켜면 수렴은 기도의 영역이 된다.
• 피로 해석 — 반복 하중에 의한 누적 손상과 수명을 평가한다. S-N 선도와 레인플로 계수법이 기본 도구다.
• 충돌·낙하 해석 — 외연적(explicit) 시간 적분으로 밀리초 단위의 충격 현상을 푼다. 자동차 충돌 안전 해석이 대표 사례다.

4. 지배 방정식

이산화를 마친 선형 정적 문제는 결국 하나의 연립방정식으로 귀결된다.

$$\mathbf{K}\mathbf{u} = \mathbf{f}$$

여기서 $\mathbf{K}$는 강성행렬, $\mathbf{u}$는 절점 변위 벡터, $\mathbf{f}$는 하중 벡터다. CFD 엔지니어의 눈에는 허무할 정도로 깔끔한 형태지만, 비선형 문제가 되면 강성이 변위에 의존하여

$$\mathbf{K}(\mathbf{u})\,\mathbf{u} = \mathbf{f}$$

가 되고, 뉴턴-랩슨법 반복으로 접선 강성을 갱신해 가며 푼다. 동적 문제는 질량 행렬 $\mathbf{M}$과 감쇠 행렬 $\mathbf{C}$가 추가된 운동 방정식

$$\mathbf{M}\ddot{\mathbf{u}} + \mathbf{C}\dot{\mathbf{u}} + \mathbf{K}\mathbf{u} = \mathbf{f}(t)$$

을 시간 적분한다. 내연적(implicit) 적분의 표준은 1959년 뉴마크가 제시한 뉴마크법($\gamma = 1/2$, $\beta = 1/4$의 평균 가속도 법이 무조건 안정)이고, 충돌 해석의 외연적 적분은 CFL 조건과 같은 원리의 안정 시간 간격 제한을 받아 가장 작은 요소가 전체 계산 시간을 지배한다.

5. 외팔보로 하는 직접 검산

구조해석 코드를 새로 붙였을 때 가장 먼저 때려 보는 것이 해석해가 있는 외팔보다. 아래는 강재 외팔보($E = 210$ GPa, $\rho = 7850$ kg/m³, $\nu = 0.3$, 길이 $L = 1$ m, 직사각 단면 $20 \times 10$ mm, 자유단 하중 $P = 10$ N)를 2절점 에르미트 보요소로 직접 풀어 본 결과다. 단면 2차 모멘트는 $I = bh^3/12 = 1.6667 \times 10^{-9}$ m⁴, 단면적은 $A = 2 \times 10^{-4}$ m².

오일러-베르누이 해석해는 $\delta = PL^3/(3EI) = 9.5238$ mm 이고, 1차 고유진동수는 자유진동 해의 초월방정식

$$f_1 = \frac{\beta_1^2}{2\pi}\sqrt{\frac{EI}{\rho A L^4}}, \qquad \cos(\beta_1 L)\cosh(\beta_1 L) + 1 = 0$$

의 첫 근 $\beta_1 L = 1.875104$ 를 써서 $f_1 = 8.3552$ Hz 다. FEM 결과는 이렇게 수렴한다.

요소 수	자유단 처짐 (mm)	해석해 상대차	1차 고유진동수 (Hz)	상대오차
1	9.523810	1.8e-16	8.3949	4.8e-03
2	9.523810	9.1e-16	8.3592	4.8e-04
4	9.523810	8.0e-15	8.3554	3.3e-05
8	9.523810	2.1e-13	8.3552	2.1e-06
16	9.523810	9.4e-13	8.3552	1.3e-07

읽을 거리가 세 개 있다.

• 정적 처짐은 요소 1개로도 기계 정밀도까지 맞는다. 3차 에르미트 형상함수가 집중하중 외팔보의 정확해(3차 다항식)를 그 안에 품고 있기 때문이다. 격자를 늘려도 더 좋아질 게 없다 — 이걸 모르고 “격자 수렴 곡선이 평평하다”고 당황하는 일이 흔하다.
• 고유진동수는 격자를 2배로 늘릴 때마다 오차가 약 16배 줄어든다($h^4$). 처짐과 달리 모드 형상은 다항식이 아니라 쌍곡선·삼각함수 조합이라 요소가 필요하다.
• 일관질량(consistent mass) FEM의 고유진동수는 항상 위에서 수렴한다. 표의 8.3949 → 8.3592 → 8.3554 가 단조 감소인 게 그것이다. 리츠 근사가 변위장을 제약해 구조를 실제보다 딱딱하게 만들기 때문. 반대로 집중질량(lumped mass)을 쓰면 아래로 내려갈 수 있어, “고유진동수가 해석해보다 높게 나왔다”는 그 자체로는 버그가 아니다.

참고로 같은 보의 오일러 좌굴 임계하중은 양단 핀이면 $P_{cr} = \pi^2 EI/L^2 = 3454$ N, 외팔보 조건($K = 2$)이면 864 N 이다. 위 해석의 10 N 은 그보다 두 자릿수 아래라서 선형 가정이 안전하다. 하중을 864 N 근처까지 올려 놓고 선형 정적 해석 결과만 보고 “변위 작네요”라고 보고하면 그게 바로 사고다.

6. 전단 변형은 언제부터 문제가 되는가

오일러-베르누이 보는 전단 변형을 아예 무시한다. 1921년 티모셴코가 보탠 전단 보정항을 넣으면 자유단 처짐은

$$\delta = \frac{PL^3}{3EI} + \frac{PL}{\kappa G A}$$

가 된다($\kappa = 5/6$: 직사각 단면 전단 보정 계수, $G = E/[2(1+\nu)] = 80.77$ GPa). 위 단면($h = 10$ mm)에서 길이만 바꿔 전단항의 비중을 재면 이렇다.

세장비 L/h	전단항이 전체 처짐에서 차지하는 비율
100	0.008%
20	0.195%
10	0.774%
5	3.03%

즉 세장비 20 이상이면 전단은 소수점 뒤로 사라지고, 5쯤 되면 3%를 먹는다. 실무 감각으로 “짧고 두꺼운 보에는 티모셴코 요소”라는 국룰이 여기서 나온다. 다만 전단 보정 계수 자체가 관례가 갈리는 값이라는 점은 알고 있어야 한다 — 고전적인 $5/6$ 대신 쿠퍼가 제안한 $\kappa = 10(1+\nu)/(12+11\nu) = 0.850$ 을 쓰면 위 비율이 약 2% 작아진다. 표는 $\kappa = 5/6$ 기준이다. 코드마다 기본값이 다르니 전단이 실제로 중요한 영역에서는 설정값을 확인하고 쓸 것. 거꾸로 세장비가 큰데 티모셴코(1차 전단 변형) 요소를 쓰면 전단 잠김(shear locking)으로 구조가 터무니없이 딱딱해져, 축소적분이나 가정 전단 변형장 같은 땜질이 필요해진다.

7. 외연적 적분의 시간 간격이라는 벽

충돌·낙하 해석의 외연적 적분은 요소를 음향파가 건너가는 시간보다 짧은 시간 간격만 허용한다. 강재의 봉파속은

$$c = \sqrt{E/\rho} = \sqrt{210\times10^9 / 7850} = 5172\ \text{m/s}$$

이므로 최소 요소 크기 5 mm 면 $\Delta t \approx L_e/c = 967$ ns 다. 100 ms 짜리 충돌 이벤트를 풀려면 약 10만 3천 스텝이고, 요소를 1 mm 로 줄이면 193 ns · 약 51만 7천 스텝이 된다. 격자를 5배 촘촘하게 하면 스텝 수도 정확히 5배 늘고, 스텝당 비용까지 늘어나므로 총 비용은 그보다 더 뛴다. 외연적 해석에서 단 하나의 찌그러진 작은 요소가 전체 해석 시간을 잡아먹는 이유가 이 식이다. 질량 스케일링(mass scaling)으로 $\Delta t$ 를 억지로 늘리는 기법이 괜히 표준으로 쓰이는 게 아니다.²

8. 주요 소프트웨어

• Abaqus — 비선형과 접촉 해석의 강자. 학계와 첨단 산업에서 선호도가 높다.
• ANSYS Mechanical — 범용성과 워크벤치 통합 환경이 무기인 업계 표준급 코드.
• NASTRAN — NASA가 1960년대에 개발한 원조 격 코드로, 항공우주 분야에서는 여전히 인증의 기준점이다.³
• 그 외 충돌 해석의 LS-DYNA, 오픈소스 진영의 CalculiX와 Code_Aster 등이 있다.

9. CFD와의 연성

유체력이 구조를 변형시키고 그 변형이 다시 유동을 바꾸는 문제는 유체-구조 연성(FSI, Fluid-Structure Interaction)으로 다룬다. 유동 해석 결과를 하중으로 한 번만 넘기는 단방향(one-way) 연성과, 매 시간 단계마다 두 솔버가 값을 주고받는 양방향(two-way) 연성이 있다. 항공기 날개의 플러터(항공탄성학의 간판 문제), 혈관 내 혈류, 풍력 블레이드 해석이 대표적인 응용이다. 두 분야의 격자와 시간 스케일이 서로 달라서, CFD 엔지니어와 구조 엔지니어가 서로의 발산을 탓하는 진풍경이 벌어지기도 한다.⁴

10. 검증과 확인

구조해석은 결과가 설계 승인 문서로 직행하는 분야라서, 검증 및 확인(V&V) 체계가 비교적 이르게 표준화됐다. 계산 고체역학 쪽 기준 문서가 ASME V&V 10 이고, 실무에서 반복해 쓰는 점검은 셋이다.

• 패치 테스트 — 요소 몇 개짜리 패치에 선형 변위장을 경계조건으로 주고 상수 응력이 정확히 재현되는지 본다. 1986년 테일러·시모·지엔키위츠·찬이 이 테스트의 만족이 수렴의 필요조건(일관성과 동치)임을 정리했고, 안정성 점검까지 묶으면 충분조건이 된다. 새 요소를 만들었으면 이것부터 통과해야 한다.
• 메시 독립성 연구 — 격자를 체계적으로 줄이며 관심 응력이 수렴하는지 확인한다. 폰 미제스 응력은 특이점(날카로운 내부 코너)에서 발산하므로, 그 지점의 값은 격자를 줄이면 끝없이 커진다. 수렴하지 않는 양을 수렴시키려 격자를 더 쪼개는 것은 노동이 아니라 오해다.
• 강체 모드 점검 — 구속을 풀고 모드 해석을 돌려 0 Hz 모드가 정확히 6개 나오는지 본다. 더 나오면 요소가 끊겼거나 자유도가 빠진 것이고, 덜 나오면 어딘가 몰래 구속이 걸린 것이다.

11. 관련 문서

• 유한요소법 · 강성행렬 · 응력과 변형률
• 모드 해석 · 좌굴 · 피로 해석 · 접촉 해석
• 뉴마크법 · 뉴턴-랩슨법 · 감쇠
• 유체-구조 연성 · 플러터 · 항공탄성학
• 검증 및 확인 · 메시 독립성 연구 · 폰 미제스 응력
• 전산유체역학 · 격자 · 경계 조건 · CFL 조건 · 전자기해석

12. Footnotes

1. 구속이 없는 자유 상태 모드 해석에서는 고유진동수 0Hz짜리 강체 모드 6개가 먼저 나온다. 이걸 보고 “구조물이 불안정하다”고 보고하는 것이 신입의 전통적인 첫 실수다. ↩

2. 질량 스케일링은 작은 요소의 밀도를 몰래 올려 $\Delta t$ 를 벌어오는 기법이다. 관성이 바뀌므로 공짜가 아니고, 추가된 질량이 전체의 몇 %인지 리포트에 찍히게 되어 있다. 그 숫자를 안 보고 넘기는 순간 해석은 “빨리 끝난 다른 문제”의 답이 된다. ↩

3. NASTRAN의 원본 소스는 NASA의 공공 자산으로 공개되어 있다. 1960년대 포트란을 읽는 고고학적 소양이 필요할 뿐이다. ↩

4. 양방향 FSI가 발산하면 원인은 대개 인터페이스의 부가질량(added mass) 효과지만, 회의실에서는 일단 상대편 솔버 설정부터 의심받는다. ↩