미분방정식 · 감쇠

감쇠

감쇠는 진동에서 에너지를 빼앗습니다. ζ를 조정하여 진동하며 돌아오는 경우와 평형점을 넘지 않고 돌아오는 경우를 비교하세요. 시간 단위 T₀는 감쇠가 없을 때의 주기입니다.

01

이 시각화에서 보는 감쇠

감쇠는 에너지 손실로 진동이 사라지는 현상입니다. 감쇠기가 달린 용수철에 매단 물체는 mx′′+cx′+kx=0m x'' + c x' + k x = 0를 따릅니다. 고유진동수 ω0=k/m\omega_0 = \sqrt{k/m}와 감쇠비 ζ=c/(2mk)\zeta = c/(2\sqrt{mk})을 쓰면 x′′+2ζω0x′+ω02x=0x'' + 2\zeta\omega_0 x' + \omega_0^2 x = 0입니다. 감쇠비가 1보다 작으면 진동하는 부족감쇠, 1이면 임계감쇠, 1보다 크면 천천히 돌아오는 과감쇠입니다. 수직으로 매단 경우 중력은 정지 위치를 mg/k만큼 낮출 뿐이며, x는 그 위치에서 잽니다. 이 페이지에서는 x = 1의 정지 상태에서 물체를 놓고, 감쇠가 강할수록 더 빨리 멈추는지 살펴봅니다.

근이 운동을 결정하며, 중근이면 t·e^(λt) 항도 생깁니다
λ2+2ζω0λ+ω02=0,λ=ω0(−ζ±ζ2−1),σ=−max⁡Re⁡λ\lambda^2 + 2\zeta\omega_0\lambda + \omega_0^2 = 0, \qquad \lambda = \omega_0\left(-\zeta \pm \sqrt{\zeta^2 - 1}\right), \qquad \sigma = -\max \operatorname{Re}\lambda

그림 오른쪽은 근 평면입니다. 두 근의 곱이 ω02\omega_0^2이므로 진동하는 동안에는 반지름 ω₀인 원 위에 있습니다. 실수부 −ζω0-\zeta\omega_0는 감소를, 허수부 ωd=ω01−ζ2\omega_d = \omega_0\sqrt{1-\zeta^2}는 진동수를 정하므로 감쇠는 한 번의 진동도 느리게 만듭니다. ζ = 1에서 두 근이 −ω₀에서 만납니다. 그 이후 한 근은 왼쪽으로 가고 다른 근은 0 쪽으로 천천히 돌아옵니다. ζ가 크면 후자는 약 −ω0/(2ζ)=−k/c-\omega_0/(2\zeta) = -k/c입니다. 운동의 가장 느린 부분은 e−σte^{-\sigma t}처럼 줄어듭니다. 이 감소율 σ는 가장 오른쪽 근과 허수축의 거리이며, 정확히 ζ = 1에서 최대입니다. x(t)의 파란 점선이 이를 보여 줍니다. 부족감쇠에서는 포락선 ±e−ζω0t/1−ζ2\pm e^{-\zeta\omega_0 t}/\sqrt{1-\zeta^2}(곡선이 봉우리 사이에서 접함), 과감쇠에서는 느린 근의 항입니다.

‘멈춘다’에는 기준이 필요합니다. 이 페이지는 |x|가 초기 변위의 1 % 이내에 들어와 이후 계속 머무르기 시작하는 가장 이른 시각 T1%T_{1\%}를 사용합니다. x(t) 아래의 띠는 이 범위를 확대한 것입니다. 극값의 시점 tn=nπ/ωdt_n = n\pi/\omega_d과 값 ∣x(tn)∣=e−ζω0tn|x(t_n)| = e^{-\zeta\omega_0 t_n}로 어느 봉우리가 띠를 벗어나는지 정확히 판정한 뒤, 닫힌형식 해에 이분법을 적용해 마지막 진입 시각을 구하므로 수치는 근삿값입니다. ζ가 조금 작아져 봉우리 하나가 더 띠를 벗어날 때마다 정착 시간이 늦어져, ζ에 따른 T1%T_{1\%} 곡선에는 점프가 생깁니다. 최솟값은 임계감쇠가 아니라 ζ∗=ln⁡100/π2+ln⁡2100≈0.826\zeta^* = \ln 100/\sqrt{\pi^2 + \ln^2 100} \approx 0.826에서 나타납니다. 한 번의 지나침이 정확히 −1 %에 닿는 경우로, 1.057 T₀ 대신 0.667 T₀입니다. 임계감쇠가 가장 빠른 것은 x = 0을 전혀 넘지 않는 운동 중에서, 정지 상태로 출발할 때뿐입니다. 띠의 폭을 바꾸면 ζ*도 달라집니다. 2 % 띠에서는 약 0.780입니다.

02

살펴볼 점

  • ‘가장 빠르다’는 두 가지 뜻이 있습니다. 감소율 σ는 정확히 ζ = 1에서 최대지만, 정착 시간은 그보다 조금 작은 ζ*에서 최소입니다. 직접 입력한 0.826은 ζ*가 아닙니다. 1.0015 %만큼 넘어가므로 정착 시간이 0.896 T₀로 뛰어오릅니다. 이 페이지는 ζ*를 정확히 유지합니다.
  • 이미 큰 감쇠비를 두 배로 늘리면 돌아오는 시간도 대략 두 배가 됩니다. 느린 근의 시간상수는 약 c/k입니다.
  • x(t)는 세 영역 모두에서 정확한 해입니다. 입력한 감쇠비의 영역은 그 수 그대로 판정하므로 0.99999999999999999는 부족감쇠입니다. 원통 안 액체는 개략적인 표현일 뿐입니다.
03

조작 방법

  • ζ 슬라이더를 움직이거나(0.1, 1, 10에 맞춰 붙음), 분수·소수를 입력하거나, ζ = 1 또는 ζ*를 누르세요. 근은 경로를 따라 드래그할 수 있습니다. −ω₀에서 붙으며 그곳에서만 가지를 바꿉니다. 빨간 시간 커서를 드래그하여 물체를 움직이세요.
  • 발표 모드는 전체 화면을 엽니다. ← →로 ζ를 슬라이더 한 단계씩 바꾸고, H로 감소율과 정착 시간을 숨깁니다.
04

관련 개념

위상 초상 페이지의 대각합–행렬식 지도에서 감쇠 진동자는 수평선 Δ = ω₀²입니다.

더 읽어 보기: Wikipedia: Damping (영어); OpenStax, University Physics Volume 1, 15.5 Damped Oscillations (영어).