고체역학

고체역학 구조해석 마지막 수정: 2026-07-10 05:41:09

고체역학 Solid Mechanics
분야	고전역학 × 연속체역학
핵심 개념	응력, 변형률, 구성방정식
대표 수치기법	유한요소법
응용	구조해석, 재료공학, 기계설계

1. 개요[편집]

고체역학(Solid Mechanics)은 외력이나 온도 변화, 구속 조건 아래에서 변형 가능한 고체가 어떻게 변형하고 어떤 내부 응력을 갖는지를 다루는 연속체역학의 한 분야이다. 유체가 전단응력에 대해 끝없이 흐른다면, 고체는 전단응력을 받아도 유한하게 변형한 뒤 평형을 이루고 버틴다. 이 “버티는 능력”을 정량화하는 것이 고체역학의 핵심 임무다.

건물이 지진에 무너지지 않을지, 자동차 프레임이 충돌에서 승객을 지킬지, 인공관절이 수백만 번의 하중을 견딜지 — 하중을 받는 고체가 있는 곳이라면 어디든 고체역학이 있다. 그리고 실제 형상은 손으로 못 풀기 때문에, 이 방정식들을 컴퓨터로 푸는 구조해석과 유한요소법이 태어났다.

2. 응력과 변형률[편집]

고체역학의 두 주인공은 응력(stress)과 변형률(strain)이다. 응력은 단위면적당 내부 힘으로, 3차원에서는 방향과 작용면을 모두 지정해야 하므로 2계 텐서인 코시 응력 텐서 $\boldsymbol{\sigma}$ 로 표현된다.

\boldsymbol{\sigma} = \begin{bmatrix} \sigma_{xx} & \tau_{xy} & \tau_{xz} \\ \tau_{yx} & \sigma_{yy} & \tau_{yz} \\ \tau_{zx} & \tau_{zy} & \sigma_{zz} \end{bmatrix}

각 성분은 앞 첨자가 작용면의 법선 방향, 뒤 첨자가 힘의 방향을 뜻한다. 대각 성분은 수직응력, 비대각 성분은 전단응력이다. 각운동량 보존에서 이 텐서는 대칭( $\tau_{xy}=\tau_{yx}$ )이 된다.

변형률은 물체가 얼마나 늘어나거나 뒤틀렸는지를 나타내는 무차원량이다. 미소변형 가정 아래서는 변위장 $\mathbf{u}$ 의 대칭 기울기로 정의한다.

\varepsilon_{ij} = \frac{1}{2}\left(\frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i}\right)

자세한 물리적 해석은 응력과 변형률 문서에서 다룬다.

3. 구성방정식 — 재료의 성격[편집]

응력과 변형률을 이어주는 관계식이 구성방정식(constitutive equation)이며, 이것이 곧 재료의 성격서다. 가장 단순하고 널리 쓰이는 것이 선형 탄성을 기술하는 후크의 법칙이다. 1차원으로는 그 유명한 $\sigma = E\varepsilon$ , 3차원 등방성 재료로는 다음과 같이 쓴다.

\sigma_{ij} = \lambda\, \varepsilon_{kk}\,\delta_{ij} + 2\mu\, \varepsilon_{ij}

여기서 $\lambda, \mu$ 는 라메 상수로, 실무에서 익숙한 영률 $E$ 와 푸아송비 $\nu$ 로부터 변환된다. 이 선형·가역 영역이 탄성(elasticity)이다. 하중을 풀면 원래 모양으로 돌아온다.

4. 탄성을 넘어서 — 소성과 비선형[편집]

하중을 계속 키우면 재료는 어느 순간 항복(yield)해 영구 변형이 남는 소성(plasticity) 영역으로 들어간다. 금속이 항복했는지 판정하는 대표적 기준이 폰 미제스 응력이다. 여러 방향 응력을 하나의 등가 스칼라로 뭉쳐, 이 값이 재료의 항복강도를 넘으면 소성이 시작됐다고 본다.

\sigma_{v} = \sqrt{\tfrac{1}{2}\big[(\sigma_1-\sigma_2)^2 + (\sigma_2-\sigma_3)^2 + (\sigma_3-\sigma_1)^2\big]}

이렇게 응력-변형률 관계가 직선을 벗어나거나, 변형이 커져 형상 자체가 크게 바뀌거나, 접촉이 개입하면 문제는 비선형 구조해석이 된다. 압축 부재가 갑자기 옆으로 휘어버리는 좌굴, 균열이 자라나는 파괴역학, 반복 하중으로 서서히 죽는 피로 해석도 모두 고체역학의 하위 세계다.

5. 지배 방정식 — 평형과 적합성[편집]

고체역학의 뼈대는 세 부류의 방정식이다.¹ 첫째는 평형방정식으로, 물체 내부 어느 점에서나 힘이 균형을 이뤄야 한다는 요구다.

\nabla \cdot \boldsymbol{\sigma} + \mathbf{b} = \mathbf{0}

$\mathbf{b}$ 는 중력 같은 체적력이다. 둘째는 변위와 변형률을 잇는 변형률-변위 관계(운동학), 셋째는 앞서 본 구성방정식이다. 이 셋에 경계조건을 더하면 문제가 닫힌다. 정적 문제에서는 보통 변위장을 미지수로 두고 이 연립계를 푼다.

6. 유한요소법과의 관계[편집]

문제는 실제 부품 형상에서는 이 편미분방정식들이 손으로 안 풀린다는 것. 그래서 등장한 결정적 무기가 유한요소법(FEM)이다. FEM은 물체를 유한 개의 작은 요소로 쪼개고, 평형방정식의 약형식(가상일 원리)을 세워 다음과 같은 거대한 선형 연립방정식으로 바꾼다.

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

여기서 $\mathbf{K}$ 가 재료·형상 정보를 담은 강성행렬, $\mathbf{u}$ 가 절점 변위, $\mathbf{f}$ 가 하중 벡터다. 이 방정식을 가우스 소거법이나 반복법으로 풀면 변위장이 나오고, 그로부터 변형률과 응력을 후처리로 뽑는다. Abaqus, NASTRAN 같은 상용 코드가 바로 이 일을 산업 규모로 수행한다.² 즉 고체역학이 방정식을 세우면, FEM이 그걸 컴퓨터로 부려먹는다.

정적 해석 말고도 갈래는 다양하다. 구조물의 고유진동수와 진동 형상을 뽑는 모드 해석, 짧은 시간에 큰 하중이 걸리는 충돌·낙하를 다루는 명시적 동해석, 부재끼리 맞닿는 비선형성을 처리하는 접촉 해석이 모두 고체역학 방정식 위에 서 있다. 무엇을 미지수로 두고 어떤 시간 적분을 쓰느냐가 다를 뿐, 바닥의 평형·구성·운동학 삼위일체는 공유한다. 국소적으로 응력이 튀는 응력집중 문제도 결국 이 틀에서 정밀 격자로 잡아낸다.

7. 유체역학과의 대칭성[편집]

고체역학과 유체역학은 사실 연속체역학이라는 같은 뿌리에서 갈라진 형제다.³ 둘 다 코시 응력 텐서, 질량·운동량 보존을 공유한다. 차이는 오직 구성방정식 하나 — 고체는 응력이 변형률에 반응하고( $\sigma \sim \varepsilon$ ), 유체는 응력이 변형률’속도’에 반응한다( $\sigma \sim \dot\varepsilon$ ). 이 한 끗 차이가 FEM과 CFD라는 서로 다른 수치해석 문화를 낳았다.

8. 관련 문서[편집]

이 세 부류 — 평형, 운동학, 구성 — 는 고체역학의 삼위일체로 불린다. 학부 재료역학 첫 학기에 이걸 못 외우면 남은 3년이 고달파진다. 반대로 이 셋의 구조만 확실히 잡으면 어떤 응용 문제든 결국 이 틀에 끼워 맞추게 된다. ↩
FEM의 개념적 뿌리는 1940~50년대 항공 구조 해석으로 거슬러 올라간다. 보잉과 UC 버클리의 클러프(R. Clough)가 “finite element”라는 이름을 붙였다. 초창기엔 손으로 강성행렬을 채웠다니, 요즘 격자를 클릭 한 번에 뽑는 세대는 감사한 줄 알아야 한다. ↩
그래서 유체-구조 상호작용(FSI)처럼 둘을 동시에 풀어야 하는 유체-구조 연성 문제가 존재한다. 두 문화의 코드를 한 방에 밀어 넣고 서로 데이터를 주고받게 하는 일인데, 대체로 엔지니어의 수명을 단축시킨다. ↩