도입 · 멈추거나, 커지거나

그네는 왜 가만두면 멈추고, 박자에 맞춰 밀면 왜 커질까?

내버려 둔 진자는 느려지다 멈춘다. 같은 진자를 정확히 맞는 박자로 밀면, 점점 더 크게 흔들리다 막대 위까지 휘감는다. 같은 하드웨어. 정반대 결과. 둘 모두를 — 그리고 자동차 서스펜션, 레이저 공동, 무너진 다리까지 — 지배하는 한 식에는 세 항이 있다 — 관성, 복원력, 그리고 속도에 의존하는 저항. 저항이 이기면 운동이 죽고, 외부 밀기가 자연 박자에 도착하면 저항이 못 따라잡고 진폭이 자란다.

감쇠가 어떻게 멈출지를 정하고, 공명이 언제 밀기가 이길지를 정한다.

위젯 — 감쇠 진동

γ (감쇠)0.10

ω₀ (자연)1.00

ω (구동)1.00

상태부족 감쇠

위치 x(t)

진폭 반응 A(ω) — 정상상태

감쇠 γ0.10구동 ω1.00

ω₀ = 1 (fixed)

식: ẍ + 2γ·ẋ + ω₀²·x = F·cos(ω·t). 구동 끔이면 시스템이 변위된 상태에서 자유롭게 — 부족 감쇠 (γ < ω₀) 면 진동하면서 지수적으로 감쇠, 임계/과 감쇠 (γ ≥ ω₀) 면 진동 없이 곧장 정지로 복귀. 구동 켬이면 아래 패널이 정상상태 진폭을 구동 진동수의 함수로 보여준다. ω를 자연 진동수 ω₀ = 1를 가로질러 끌어보자 — 주황 점이 반응 곡선을 오른다. 약한 감쇠에서 피크가 날카롭고 높음. 강한 감쇠에서는 피크가 평평. *같은 식, 세 가지 행동* — 어떤 진동자든 이 곡선 어디에 살게 할지가 공학의 선택이다.

흐름

모든 진동자를 돌리는 식

스프링 위 질량은 뉴턴의 법칙을 따른다 — $m · \ddot{x} = F_net$ . 알짜힘에는 세 조각. 스프링은 변위에 비례해 당긴다 — $−k · x$ . 감쇠 는 속도에 비례해 운동을 거스른다 — $−c · \dot{x}$ . 외부 밀기 (그네에 손, 주기적 구동력) 가 $F(t)$ 를 더한다. 합치고 질량으로 나누면 —

$\ddot{x} + 2γ · \dot{x} + ω₀² · x = F(t)/m$

— 여기서 $γ = c / (2m)$ , $ω₀ = \sqrt{k/m}$ . 같은 식이 작은 각 진자 ( $ω₀ = \sqrt{g/L}$ ), LC 회로 ( $ω₀ = 1/\sqrt{LC}$ ), 바람 속에서 흔들리는 건물, 기타 줄, 빛에 반응하는 원자의 전자 구름 — 모두를 묘사한다. 다른 물리, 같은 수학 — 모든 물리에서 가장 많이 재사용되는 식 중 하나.

자유 감쇠 — 멈추는 세 가지 방식

$F(t) = 0$ (외부 밀기 없음), 질량을 변위시키고 놓아준다. 좌변의 미분들이 우변의 복원항과 경쟁하고, 답은 $γ$ 가 $ω₀$ 에 비해 얼마나 큰지에 달려 있다 —

부족 감쇠 ( $γ < ω₀$ ): 질량이 줄어들면서 진동 — 지수 포락선 $e^(−γt)$ 곱하기 사인. 위젯에서 약한 감쇠가 이 모양 — 자취가 0으로 진동하며 수렴.
임계 감쇠 ( $γ = ω₀$ ): 진동 없이 가장 빠른 정지 복귀. 자동차 충격 흡수기와 도어 클로저가 정확히 여기 맞춰져 있다 — 덜 감쇠하면 흔들리고, 더 감쇠하면 느려진다.
과 감쇠 ( $γ > ω₀$ ): 질량이 임계보다 느리게 정지로 돌아간다. 진동을 절대 원치 않을 때 유용 — 정밀 기기, 일부 제어 시스템.

위젯의 자유 감쇠 자취 (구동 끔) 가 세 영역을 즉시 보여준다 — 같은 식, 세 가지 다른 “정지” 모양. 수학적으로는 모두 특성 방정식 $r² + 2γr + ω₀² = 0$ 이 지배 — 판별식 $γ² − ω₀²$ 이 양수, 0, 음수인지가 그 차이.

import numpy as np

# Equation: ẍ + 2γ·ẋ + ω₀²·x = F·cos(ω·t)
# Three free-decay regimes (F = 0), determined by γ vs ω₀.
omega0 = 1.0

def simulate(gamma, omega_drive=0.0, F=0.0, T=30, dt=0.05, x0=1, v0=0):
    """RK4 integration of the damped driven oscillator."""
    n = int(T / dt) + 1
    t = np.zeros(n); x = np.zeros(n); v = np.zeros(n)
    x[0], v[0] = x0, v0
    accel = lambda x_, v_, t_: (-2*gamma*v_ - omega0**2*x_
                                 + F*np.cos(omega_drive*t_))
    for i in range(1, n):
        t[i] = i * dt
        k1x, k1v = v[i-1], accel(x[i-1], v[i-1], t[i-1])
        k2x = v[i-1] + dt/2 * k1v
        k2v = accel(x[i-1] + dt/2 * k1x, v[i-1] + dt/2 * k1v, t[i-1] + dt/2)
        k3x = v[i-1] + dt/2 * k2v
        k3v = accel(x[i-1] + dt/2 * k2x, v[i-1] + dt/2 * k2v, t[i-1] + dt/2)
        k4x = v[i-1] + dt * k3v
        k4v = accel(x[i-1] + dt * k3x, v[i-1] + dt * k3v, t[i])
        x[i] = x[i-1] + dt/6 * (k1x + 2*k2x + 2*k3x + k4x)
        v[i] = v[i-1] + dt/6 * (k1v + 2*k2v + 2*k3v + k4v)
    return t, x

# Three free-decay regimes from displaced start (x=1, v=0):
[(name, simulate(gamma=g)[1][-1]) for name, g in
 (("undamped", 0.0), ("under-damped", 0.1),
  ("critical", 1.0), ("over-damped", 2.0))]
# undamped:    final x oscillates near ±1   (no decay)
# under:       final x near 0 after many cycles, but visibly oscillates
# critical:    final x essentially 0, no oscillation (fastest return)
# over:        final x slightly above 0, slow exponential return

이 식이 작동하는 이유 — 작은 각 선형화

작은 각 진자는 거짓말 위에 산다 — 진자시계 페이지의 선형화 트릭. 진짜 식은 $\sin θ$ 를 가지는데 비선형. 작은 흔들림에서 $\sin θ ≈ θ$ 로 바꾸면 위 선형 진동자가 된다. 일상 공학에서 작동하는 실제 감쇠 진동자들은 거의 전부 이 영역에 산다 — 자동차 서스펜션은 몇 센티미터 변위, 건물 기둥은 천분의 일, 마이크 진동판은 마이크론. 거짓말이 통하는 영역에서 살기 때문에 선형식을 따른다.

거짓말이 깨지면 — 큰 각 진자, 비후크 스프링, 트랜지스터 포화 근처의 큰 전기 신호, 탄성 한계를 넘은 건물 — 선형 공식이 틀려진다. 1940년 타코마 협곡 다리 붕괴가 정확히 이것으로 유명 — 작은 진동은 예측대로였지만 진폭이 커지자 공기역학 구동력이 선형 분석이 못 본 비선형 임계를 넘었다. § 1 식은 영역 안에서만 정확. 밖에서는 다른 문제를 풀고 있는 셈.

강제 — 진폭 = 구동 진동수의 함수

이제 구동을 켜자 — $F(t) = F · \cos(ω · t)$ . 과도 응답 (자유 감쇠 부분) 이 사라지면 질량은 구동 진동수에서 정상상태로 안정 — $x(t) = A(ω) · \cos(ωt − φ)$ . 진폭 $A(ω)$ 는 닫힌 형식 — 정상상태 가정을 식에 넣고 cos / sin 모아 풀면 —

$A(ω) = F / \sqrt{(ω₀² − ω²)² + (2γω)²}$

두 극한이 직접 읽힌다. 매우 느린 구동 ( $ω → 0$ ) 에서 $A → F/ω₀²$ — 정적 변위 — 스프링이 천천히 가해지는 힘과 균형. 매우 빠른 구동 ( $ω → ∞$ ) 에서 $A → 0$ — 질량이 따라가지 못함, 관성이 이긴다. 그 사이에 곡선은 피크를 가진다 — 그것이 공명.

위젯 아래 패널이 현재 $γ$ 의 $A(ω)$ 를 보여준다. 감쇠를 0.05까지 끌면 — 피크가 솟구침, 좁고 높음. 감쇠를 1.0까지 끌면 — 피크가 사라짐, 곡선이 단조. 감쇠가 가벼울수록 피크가 더 날카롭고 높음. 진폭은 구동이 펌프하는 에너지 대 감쇠가 빼는 에너지의 균형 — 그 평형이 정상상태 진폭.

# Steady-state amplitude as a function of driving frequency:
# A(ω) = F / sqrt((ω₀² − ω²)² + (2γω)²)
# — a closed-form 'frequency response' that emerges by plugging
# x(t) = A·cos(ωt − φ) into the forced equation and collecting cos / sin.
def amplitude(omega, gamma, omega0=1.0, F=1.0):
    a = omega0**2 - omega**2
    b = 2 * gamma * omega
    return F / np.sqrt(a*a + b*b)

# Peak location (where the response is largest):
# d/dω [A] = 0  →  ω_peak = sqrt(ω₀² − 2γ²)   (for γ < ω₀/√2)
# For very small γ, ω_peak ≈ ω₀ — the natural frequency.
def peak_omega(gamma, omega0=1.0):
    if gamma >= omega0 / np.sqrt(2):
        return 0.0  # no peak — overdamped frequency response
    return np.sqrt(omega0**2 - 2*gamma**2)

[(g, peak_omega(g), amplitude(peak_omega(g) or 0.001, g))
 for g in (0.05, 0.1, 0.3, 0.7, 1.0)]
# γ=0.05  → peak at 0.997, A ≈ 10        sharp resonance
# γ=0.1   → peak at 0.990, A ≈ 5.0
# γ=0.3   → peak at 0.906, A ≈ 1.7
# γ=0.7   → peak at ~0,    A ≈ 1.4       no real resonance peak
# Lighter damping → sharper peak. The Q-factor 1/(2γ/ω₀) measures this
# 'sharpness' directly; a high-Q oscillator (laser cavity, atomic clock)
# is the same equation with γ pushed near zero.

공명 — 박자에 맞춰 민다

왜 피크가 $ω ≈ ω₀$ 에 있을까? 공명 은 매 번의 밀기가 진동자가 이미 같은 방향으로 움직이고 있을 때 도착하는 일 — 매 사이클마다 에너지가 동기적으로 들어간다. 공명에서 벗어나면, 밀기가 번갈아 에너지를 더하고 빼서 순 입력이 평균적으로 0. 아래 수학 — 평균 입력 일률은 $⟨F · \dot{x}⟩$ , 한 사이클에 대한 적분 — 이 평균은 구동력과 속도가 위상이 맞을 때 가장 크고, 그 일이 $ω$ 가 $ω₀$ 에 닿을 때 일어난다.

같은 항등식이 — 오페라 가수가 유리잔을 깨뜨리는 것 (성대가 유리의 자연 진동수 에 구동), 타코마 협곡 다리 붕괴 (바람 구동 공기탄성 플러터, 다리의 자연 진동수 중 하나), MRI 핵 뒤집힘 (라모르 진동수 RF 펄스), 하늘 가득한 신호에서 라디오 다이얼이 한 방송을 골라내기 (한 ω₀에 맞춰진 LC 회로) — 모두를 지배한다. 공명은 특별한 힘이 아니라, 시스템이 이미 가진 박자에 에너지를 주입하는 일.

Q-팩터 $Q = ω₀ / (2γ)$ 가 공명 피크의 좁기를 측정한다. 높은 Q (원자 시계: Q ≈ 10¹⁰; 진자: Q ≈ 100) 면 시스템이 정확히 한 진동수에서 깔끔하게 울린다. 낮은 Q면 넓고 쉽게 구동되는 반응. Q 선택이 진동자 설계자의 첫 결정이고, 나머지는 그 따름결과.

# Resonance with energy bookkeeping: average the power F·v over a few
# cycles and watch energy build up coherently when ω = ω₀.
def energy_input(gamma, omega_drive, F=1.0, T=50, dt=0.01):
    omega0 = 1.0
    t = np.arange(0, T, dt)
    A = amplitude(omega_drive, gamma)
    a = omega0**2 - omega_drive**2
    b = 2 * gamma * omega_drive
    phi = np.arctan2(b, a)
    v = -A * omega_drive * np.sin(omega_drive * t - phi)
    forcing = F * np.cos(omega_drive * t)
    return float(np.mean(forcing * v))

[(omega_drive, energy_input(0.1, omega_drive))
 for omega_drive in (0.5, 0.9, 1.0, 1.1, 1.5)]
# ω=0.5 → avg power ≈ 0.05
# ω=0.9 → avg power ≈ 0.45
# ω=1.0 → avg power ≈ 2.5     ← resonance: 50× more energy in / cycle
# ω=1.1 → avg power ≈ 0.45
# ω=1.5 → avg power ≈ 0.05
# 'Pushing on the beat' is literally an integral identity — F·v averages
# to a positive number only when forcing is in phase with velocity.

이제 깨봐

공명이 항상 재앙은 아니다. “다리 위에서 군인은 발걸음을 풀어야 한다”는 잘해야 반쪽 진실 — 동기화된 발자국이 낮은 Q 다리를 자연 진동수 근처에서 구동할 수도 있지만, 대부분의 다리는 이게 무관할 만큼 감쇠가 충분하다. 공명은 또 라디오를 작동하게 만들고, MRI가 너를 영상화하고, 전자레인지가 물을 데우고 (물의 회전 공명 중 하나에서의 전자기파), 레이저 공동의 원리 전체다. 위험한 공명과 유용한 공명은 같은 식이다. 공학적 질문은 식의 어느 쪽에 살 것인가. 높은 Q 진동자는 원할 때는 도구, 원치 않을 때는 폭탄.

감쇠가 어떻게 멈출지를 정하고, 공명이 언제 밀기가 이길지를 정한다. 한 식 — $\ddot{x} + 2γ·\dot{x} + ω₀²·x = F(t)$ — 이 에너지를 잃으며 맞는 박자를 듣는 모든 진동자를 돌린다. 세 가지 자유 감쇠 영역 (부족 / 임계 / 과) 과 한 정상상태 진폭 곡선 $A(ω) = F / \sqrt{(ω₀² − ω²)² + (2γω)²}$ . 나머지는 공학 — 어떤 γ, 어떤 ω₀, 어떤 Q.

exercises · 손으로 풀기

1감쇠 없는 주기 읽기

위젯에서 감쇠 $γ = 0$ (감쇠 없음 프리셋). 자취가 완전한 코사인이 된다. 주기를 시간축에서 읽어라. $ω₀ = 1$ 로 식 $T = 2π / ω₀$ 와 일치하는지 확인.

2임계 감쇠가 가장 빠르고, 진동 없음계산기 없이

위젯의 세 감쇠 프리셋: 약함 (γ = 0.1), 임계 (γ = 1.0), 과 감쇠 (γ = 2.0). 각각에 대해 자취가 진동할지, 0 근처로 안정되는 데 대략 얼마나 걸릴지 예측. 위젯에서 확인.

3공명 — 더 작은 γ에서 더 날카로움

$γ = 0.1$ 과 $γ = 0.3$ 에 대해 식 $A(ω₀) = F / (2γω₀)$ 로 $A(ω₀)$ 를 계산 ( $F = 1, ω₀ = 1$ ). 감쇠를 반으로 줄이면 왜 피크 진폭이 두 배 이상 되는가?

4사악한 것 · 공명은 항상 위험한가?

친구가 말한다 — “공명은 위험하다 — 다리가 무너지고 가수가 유리를 깨뜨리는 이유다.” 유용한 공명 세 곳을 들어 한 문장으로 반박하라.

왜 교과서는 이렇게 안 가르치나

입문 물리는 감쇠 없는 조화 진동자 (후크 법칙 + 뉴턴 II) 를 장식 예제로 가르치고, 그다음 파동 방정식, 광학, 양자로 건너뛴다. 감쇠 강제 케이스 — 세 영역, 닫힌 형식 진폭 반응, 정직한 공명 정의 — 는 2학년 역학 또는 전기공학 강의에 남겨지고, 거기서 라플라스 변환에 빠진다. Lemma는 그것을 무너뜨린다 — 세 항의 한 식, 판별식 부호로 갈리는 세 영역, 한 진폭 공식, 그리고 모든 것을 직접 보여주는 위젯. Q-팩터와 공명 진동수가 한 계산으로 따라 나오고, 다리·라디오·레이저·자동차 서스펜션이 같은 축 위의 다른 선택으로 등장한다. 이 식은 전기공학보다 오래됐고 그 안 모든 곳에서 쓰인다 — ‘특별하고 고급’이라 가르치는 건 커리큘럼의 인공물.

용어집 · 이 페이지에서 쓰임 · 3

damping·감쇠

진동자에서 에너지를 빼내 결국 정지에 이르게 하는 모든 메커니즘. 가장 단순한 모형은 *속도에 비례*하는 것 — 운동 반대 방향으로 힘 `−c · ẋ` — 으로, 식 `m ẍ + c ẋ + k x = 0`을 선형이고 풀기 쉽게 만든다. `c`가 `√(4mk)`와 비교해 어떤지에 따라 세 영역 — _부족 감쇠_ (줄어들면서 진동), _임계 감쇠_ (가장 빠르게 정지, 진동 없음), _과 감쇠_ (느리게 정지, 진동 없음). 감쇠는 *공기 저항*보다 더 일반적이다 — 공기 저항은 한 특정 원천 (유체 저항) 이고, 감쇠는 마찰, 점성 흐름, 전기 저항, 구조적 이력 등 진동자의 에너지를 가져가서 돌려주지 않는 모든 것을 포괄한다.

resonance·공명

외부 구동력이 진동자의 _자연 진동수_ 근처에서 (또는 정확히 그곳에서) 밀 때 일어나는 일. 매 번의 밀기가 진동자가 이미 같은 방향으로 움직이고 있을 때 도착 — 에너지가 동기적으로 들어가고, 진폭이 자란다. 공명에서 벗어나면, 밀기가 번갈아 에너지를 더하고 빼면서 반응이 작게 유지된다. 구동 진동수 `ω`의 함수로서 정상상태 진폭은 `A(ω) = F / √((ω₀² − ω²)² + (2γω)²)`, 가벼운 감쇠에서는 `ω = ω₀` 근처에 피크. _감쇠가 가벼울수록 더 날카롭고 높은 피크._ 같은 항등식이 오페라 가수가 유리잔을 깨뜨리는 것, 타코마 협곡 다리 붕괴, MRI 핵의 뒤집힘, 하늘 가득한 신호에서 라디오 다이얼이 한 방송을 골라내는 것 — 모두를 지배한다. 공명은 *특별한 힘*이 아니라, 시스템이 이미 가지고 있는 박자에 에너지를 주입하는 일.

natural frequency·자연 진동수

외부 구동도 감쇠도 없이 시스템이 _스스로_ 진동하는 진동수. 질량-스프링: `ω₀ = √(k/m)`. 작은 각 진자: `ω₀ = √(g/L)`. LC 회로: `ω₀ = 1/√(LC)`. 다른 물리, 같은 모양 — `ω₀² = 강성 / 관성`. 감쇠는 겉보기 진동수를 약간 아래로 옮긴다. 외부 구동력은 `ω₀` 자체를 바꾸지 _않지만_, 구동 진동수가 `ω₀`에 가까워지면 큰 반응을 일으킨다 (_공명_). 모든 시스템은 적어도 하나의 자연 진동수를 가지고, 복잡한 시스템 (다리, 기타 줄, 건물) 은 여러 개를 가진다.