감쇠

구조해석 고체역학 마지막 수정: 2026-07-13 04:29:44

1. 개요[편집]

감쇠(damping)는 진동하는 시스템에서 운동 에너지가 열이나 다른 형태로 흩어져 진동 진폭이 시간에 따라 줄어드는 현상, 또는 그 소산 효과를 정량화한 계수를 말한다. 감쇠가 없다면 한 번 때린 종은 영원히 울고, 지진에 흔들린 건물은 결코 멈추지 않는다. 현실의 모든 구조물은 재료 내부 마찰, 접합부 미끄러짐, 주변 유체 저항 등으로 에너지를 새어 흘리며, 이를 수치 모델에 어떻게 집어넣느냐가 동적 구조해석의 정확도를 좌우한다. 문제는 감쇠가 강성이나 질량과 달리 물리적 근원이 제각각이라, “측정은 되는데 이론적으로 딱 떨어지지 않는” 구조해석계의 미운 오리 새끼라는 점이다.¹

가장 다루기 쉬운 형태는 속도에 비례하는 점성 감쇠(viscous damping)다. 1자유도 감쇠 진동계의 운동방정식은 질량 $m$ , 감쇠 계수 $c$ , 강성 $k$ 로 이렇게 쓴다.

m \ddot{x} + c \dot{x} + k x = f(t)

이 방정식 하나에 진동공학의 거의 모든 직관이 담겨 있다.

2. 감쇠비와 임계 감쇠[편집]

감쇠의 세기를 무차원으로 나타낸 것이 감쇠비(damping ratio) $\zeta$ 다. 감쇠 계수를 그것의 임계값으로 나눈 값으로 정의된다.

\zeta = \frac{c}{c_{cr}}, \qquad c_{cr} = 2\sqrt{km} = 2 m \omega_n

여기서 $\omega_n = \sqrt{k/m}$ 은 비감쇠 고유진동수다. $\zeta$ 값에 따라 시스템의 성격이 완전히 갈린다.

$\zeta < 1$ (부족 감쇠, underdamped): 진폭이 줄어드는 진동을 한다. 대부분의 구조물이 여기 속하며, 강철 구조는 $\zeta$ 가 대략 0.01~0.05로 매우 작다.
$\zeta = 1$ (임계 감쇠, critically damped): 진동 없이 가장 빠르게 평형으로 되돌아온다. 자동차 서스펜션이나 계측기 바늘이 노리는 지점.
$\zeta > 1$ (과도 감쇠, overdamped): 진동 없이 느릿하게 되돌아온다. 문이 쾅 닫히지 않게 하는 도어 클로저가 대표적.

부족 감쇠계는 매 진동마다 진폭이 지수적으로 감소하며, 감쇠된 고유진동수는 $\omega_d = \omega_n \sqrt{1 - \zeta^2}$ 로 원래보다 살짝 낮아진다. $\zeta$ 를 실험에서 뽑아내는 고전적 수법이 연속 진폭의 비를 로그로 재는 대수 감소율(logarithmic decrement)이다.²

3. 감쇠의 종류[편집]

물리적 근원에 따라 감쇠는 여러 얼굴을 가진다. 수치 모델링의 난이도도 제각각이다.

점성 감쇠: 속도에 비례( $f_d = c\dot{x}$ ). 유체 저항이 원인이며 수학적으로 가장 순하다. 그래서 실제 근원이 뭐든 일단 점성으로 근사하는 게 국룰이다.
구조 감쇠(이력 감쇠, hysteretic damping): 재료가 하중-변형 곡선을 그릴 때 매 주기 이력 루프(hysteresis loop)만큼 에너지를 먹는다. 소산 에너지가 진동수와 무관하고 변형 진폭에만 비례하는 것이 점성과 다른 결정적 특징. 주파수 영역 해석에서 복소 강성 $k(1 + i\eta)$ 로 표현하며 $\eta$ 를 손실 계수(loss factor)라 부른다.
쿨롱 감쇠(Coulomb damping): 두 면이 마른 채로 미끄러질 때의 건마찰. 힘의 크기가 일정하고 속도의 부호만 따라간다( $f_d = -\mu N \,\text{sgn}(\dot{x})$ ). 비선형이라 다루기 까다롭고, 진폭이 선형이 아니라 계단식으로 감소한다.

현실의 구조물은 이 셋이 뒤섞여 있다. 그래서 엔지니어는 대개 계측된 총 감쇠를 “등가 점성 감쇠비” 하나로 뭉뚱그려 처리한다. 정확도보다 실용성이 이기는 대표적 사례.³

4. 레일리 감쇠[편집]

유한요소법으로 이산화한 다자유도 구조의 운동방정식은 행렬 형태로 이렇게 된다.

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

여기서 골칫거리는 감쇠 행렬 $\mathbf{C}$ 다. 질량 행렬 $\mathbf{M}$ 과 강성행렬 $\mathbf{K}$ 는 재료·형상에서 자연스럽게 조립되지만, $\mathbf{C}$ 를 채울 물리적 레시피가 없다. 가장 널리 쓰는 편법이 레일리 감쇠(Rayleigh damping), 즉 감쇠 행렬을 질량과 강성의 선형 결합으로 가정하는 것이다.

\mathbf{C} = \alpha \mathbf{M} + \beta \mathbf{K}

이 형태가 사랑받는 이유는 순전히 수학적 편의다. 이렇게 두면 $\mathbf{C}$ 가 모드 좌표에서 대각화되어(비례 감쇠, proportional damping) 각 진동 모드가 서로 얽히지 않고 독립적인 1자유도계로 분리되기 때문이다. 각 모드의 감쇠비는 다음처럼 진동수의 함수가 된다.

\zeta_i = \frac{\alpha}{2\omega_i} + \frac{\beta \omega_i}{2}

$\alpha$ 항(질량 비례)은 저주파 모드를, $\beta$ 항(강성 비례)은 고주파 모드를 주로 잡는다. 실무에서는 관심 있는 두 진동수에서 원하는 $\zeta$ 를 지정해 $\alpha, \beta$ 를 역산한다. 다만 그 두 점 사이 진동수에서만 얼추 맞고 나머지는 어긋난다는 근본적 한계가 있다. 더 정밀하게 각 모드마다 감쇠비를 따로 주고 싶으면 모드 해석으로 얻은 모드 형상을 이용해 모드 감쇠(modal damping)를 직접 지정한다.⁴

5. 시간 적분과 주파수 응답에서의 역할[편집]

감쇠는 두 갈래의 동적 해석에서 서로 다른 방식으로 작동한다.

시간 영역에서는 뉴마크법 같은 시간 적분법으로 매 시간 스텝의 응답을 순차적으로 계산한다. 이때 감쇠 행렬 $\mathbf{C}$ 는 매 스텝의 유효 강성에 녹아 들어가 응답을 감쇠시킨다. 흥미로운 점은 시간 적분 기법 자체가 만들어내는 수치적 감쇠(numerical damping)도 있다는 것 — 뉴마크법의 파라미터를 조절하면 고주파 수치 노이즈를 인위적으로 걷어낼 수 있는데, 이건 물리적 감쇠가 아니라 알고리즘의 부산물이다. 이 둘을 헷갈리면 결과를 오독하기 십상이다.

주파수 영역에서는 하모닉 하중에 대한 정상 응답, 즉 주파수 응답 함수(FRF)를 구한다. 감쇠는 여기서 공진(resonance) 봉우리의 높이를 결정하는 핵심 인자다. 감쇠가 없으면 공진에서 응답이 무한대로 발산하지만, 감쇠비 $\zeta$ 가 있으면 공진 진폭이 대략 $1/(2\zeta)$ 배로 유한하게 억제된다. 즉 감쇠가 작을수록 공진 봉우리는 날카롭고 높아진다. 지진·바람·기계 진동에서 구조물이 살아남느냐를 가르는 게 바로 이 공진 근처 감쇠라, 감쇠는 동적 구조해석에서 결코 대충 넘길 수 없는 파라미터다.

6. 관련 문서[편집]

질량과 강성은 재료 물성표에서 바로 따오지만, 감쇠는 “일단 실제로 흔들어 보고” 역산하는 경우가 대부분이다. 그래서 설계 단계에서 감쇠비는 경험값이나 규준값(예: 콘크리트 5%, 강 2%)을 그냥 갖다 쓴다. 진동공학에서 가장 자주 “적당히” 처리되는 값. ↩
대수 감소율 $\delta = \ln(x_n / x_{n+1})$ 로 정의하며, $\zeta$ 가 작을 때 $\zeta \approx \delta / 2\pi$ 로 간단히 환산된다. 종을 한 번 치고 소리가 잦아드는 곡선만 녹음해도 감쇠비가 나온다는 뜻이라, 실험이 우아하다. ↩
등가 점성 감쇠는 “한 주기 동안 소산되는 에너지가 같도록” 실제 감쇠를 가상의 점성으로 치환한 것이다. 진폭이 바뀌면 등가값도 바뀌니 엄밀히는 근사지만, 선형 해석 도구를 그대로 쓸 수 있다는 실용성이 모든 이론적 찜찜함을 이긴다. ↩
레일리 감쇠의 $\alpha \mathbf{M}$ 항은 물리적으로 좀 이상하다. 구조 전체가 마치 끈적한 매질에 잠긴 것처럼 강체 운동조차 감쇠시키기 때문. 그래서 순수 $\beta \mathbf{K}$ 만 쓰거나, 모드별로 감쇠비를 직접 박는 모드 감쇠를 선호하는 해석자도 많다. ↩

감쇠