카피차 진자

편집 역사 토론
계산물리 물리 수치해석 마지막 수정: 2026-08-23 04:23:07

1. 개요[편집]

카피차 진자
Kapitza pendulum
구성지지점을 연직으로 고속 진동시킨 강체 진자
현상거꾸로 선 평형점이 안정해짐
이론A. Stephenson (1908) · P. Kapitza (1951)
안정 조건(Aω)² > 2gl
상한A/l ≲ 0.45 (주기배가로 재불안정)
도구빠른/느린 분리 · 유효 퍼텐셜 · 플로케 이론
같은 원리파울 트랩 · 폰데로모티브 힘

카피차 진자(Kapitza pendulum)는 지지점을 연직 방향으로 충분히 빠르고 작게 흔들면 원래 불안정하던 거꾸로 선 평형점이 안정해지는 강체 진자로, 고주파 진동이 느린 운동에 남기는 평균 힘 — 즉 동역학적 안정화(dynamic stabilization)의 교과서 표본이다. 표트르 카피차가 1951년에 유효 퍼텐셜로 설명하고 실물 장치를 만들어 보이면서 이름이 붙었지만, 안정화 자체는 1908년 스티븐슨이 이미 수학적으로 지적했다.1

직관을 정면으로 거스른다는 것이 이 계의 매력이다. 흔들기의 시간 평균 힘은 정확히 0이다. 위로 밀린 만큼 아래로 밀린다. 그런데도 진자는 선다. 답은 평균이 0인 것과 효과가 0인 것은 다르다는 데 있다 — 흔들리는 강제력이 자기가 만든 진동 응답과 상관되면서, 평균이 0이 아닌 잔여력을 남긴다. 이 잔여력이 만드는 것이 유효 퍼텐셜(effective potential)이고, 거기서 위쪽이 골짜기가 된다.

이 문서는 현상 자체와 그 정량적 조건, 그리고 근사가 어디서 깨지는지를 다룬다. 평균화 정리의 일반론과 2차 평균의 형식적 기계는 이중 평균법이, 선형화 후의 안정선도·플로케 판정은 마티외 방정식이 소스이므로 그쪽으로 넘긴다.

2. 운동방정식[편집]

길이 ll, 질량 mm의 진자를 생각한다. 지지점의 연직 위치를 yp(t)=Acos⁡ωty_p(t) = A\cos\omega t로 흔든다. 지지점에 붙은 (비관성) 좌표계에서 보면 관성력이 중력에 그대로 더해지므로, 유효 중력가속도가

geff(t)=g+y¨p=g−Aω2cos⁡ωtg_{\text{eff}}(t) = g + \ddot y_p = g - A\omega^2\cos\omega t

로 시간에 따라 요동친다. 아래쪽 연직에서 잰 각을 θ\theta라 하면

θ¨=−1l(g−Aω2cos⁡ωt)sin⁡θ\ddot\theta = -\frac{1}{l}\bigl(g - A\omega^2\cos\omega t\bigr)\sin\theta

가 전부다. 여기서 눈여겨볼 것이 강제 진폭에 ω2\omega^2이 붙는다는 점이다. AA가 아무리 작아도 ω\omega를 올리면 순간 가속도는 얼마든지 커진다. “작게, 그러나 빠르게”라는 처방이 모순이 아닌 이유다.

지배적인 무차원수는 두 개다.

ε=Al≪1,ωω0≫1,ω0=g/l\varepsilon = \frac{A}{l} \ll 1, \qquad \frac{\omega}{\omega_0} \gg 1, \quad \omega_0=\sqrt{g/l}

ε\varepsilon은 흔들기 진폭이 진자보다 훨씬 작다는 것, 두 번째는 흔들기가 진자의 고유진동보다 훨씬 빠르다는 것이다. 뒤에 보겠지만 두 조건이 다 필요하다 — 하나만 만족하면 안정화가 안 되거나, 됐다가 다시 깨진다.

3. 빠른 변수와 느린 변수[편집]

각을 두 조각으로 쪼갠다. 느린 평균 운동 Θ\Theta와, 흔들기 주기에 맞춰 떠는 작은 진동 ξ\xi.

θ(t)=Θ(t)+ξ(t),∣ξ∣≪1,⟨ξ⟩=0\theta(t) = \Theta(t) + \xi(t), \qquad |\xi|\ll1, \quad \langle\xi\rangle = 0

⟨⋅⟩\langle\cdot\rangle은 빠른 주기 2π/ω2\pi/\omega에 대한 평균이다.

빠른 조각. ξ\xi의 방정식에서 ξ¨\ddot\xi와 균형을 이룰 만큼 큰 항은 ω2\omega^2이 붙은 강제항뿐이고, 그 계수 안의 Θ\Theta는 이 시간척도에서 상수로 굳어 있다.

ξ¨≃Aω2lcos⁡(ωt)sin⁡Θ⟹ξ(t)=−Alcos⁡(ωt) sin⁡Θ\ddot\xi \simeq \frac{A\omega^2}{l}\cos(\omega t)\sin\Theta \qquad\Longrightarrow\qquad \xi(t) = -\frac{A}{l}\cos(\omega t)\,\sin\Theta

두 번 적분하면서 ω2\omega^2이 정확히 상쇄된다. 즉 떨림의 크기는 진동수와 무관하게 εsin⁡Θ\varepsilon\sin\Theta 정도다. Θ=0\Theta=0이나 π\pi, 즉 연직에 정확히 서 있으면 떨림이 아예 없다는 것도 읽힌다.

느린 조각. 원래 식의 sin⁡θ\sin\theta를 sin⁡Θ+ξcos⁡Θ\sin\Theta + \xi\cos\Theta로 전개하고 한 주기 평균한다. ⟨cos⁡ωt⟩=0\langle\cos\omega t\rangle = 0이라 강제항의 1차 성분은 통째로 죽고, 강제와 떨림의 상관항만 살아남는다.

⟨cos⁡(ωt) ξ(t)⟩=−Alsin⁡Θ ⟨cos⁡2ωt⟩=−A2lsin⁡Θ\bigl\langle \cos(\omega t)\,\xi(t)\bigr\rangle = -\frac{A}{l}\sin\Theta\,\bigl\langle\cos^2\omega t\bigr\rangle = -\frac{A}{2l}\sin\Theta

⟨cos⁡2⟩=1/2\langle\cos^2\rangle = 1/2라는 초등적인 사실 하나가 이 현상 전체를 떠받친다. 결과는

Θ¨=−glsin⁡Θ−A2ω22l2sin⁡Θcos⁡Θ\ddot\Theta = -\frac{g}{l}\sin\Theta - \frac{A^2\omega^2}{2l^2}\sin\Theta\cos\Theta

새로 생긴 둘째 항은 A2ω2A^2\omega^2에 비례한다. 강제력의 평균이 아니라 제곱의 평균이 들어왔다는 것이 요점이며, 평균이 0인데 효과가 0이 아닌 이유가 정확히 이것이다. 형식적으로 이것은 2차 평균화가 남기는 항이고, 왜 1차에서는 아무것도 안 나오는지에 대한 정직한 정리가 이중 평균법에 있다.

4. 유효 퍼텐셜과 안정 조건[편집]

위 방정식은 다음 퍼텐셜 안에서의 보통 운동이다.

Ueff(Θ)=−mglcos⁡Θ+mA2ω24sin⁡2ΘU_{\text{eff}}(\Theta) = -mgl\cos\Theta + \frac{mA^2\omega^2}{4}\sin^2\Theta

Θ¨=− Ueff′(Θ)/(ml2)\ddot\Theta = -\,U_{\text{eff}}'(\Theta)/(ml^2)을 확인하면 바로 맞는다.2 이 꼴은 하전입자가 비균질 고주파 전자기장에서 받는 폰데로모티브 퍼텐셜 ⟨f2⟩/(2mω2)\langle f^2\rangle/(2m\omega^2)와 정확히 같은 구조다 — 진동하는 힘이 센 쪽이 언덕이 된다.

지지점을 y_s = a cos Ωt 로 흔드는 강체 진자를 RK4(구동주기당 200스텝 이상)로 적분하면서, 같은 화면에 시간평균 유효 퍼텐셜 U_eff/(mgl) = −cosθ + (a²Ω²/4gl)sin²θ 를 함께 그린다. Ω/ω₀ = 20 에서 a/l 을 0.065 → 0.075 로 올리면 (aΩ)² 가 2gl 을 넘어서면서 θ=π 가 U_eff 의 극대에서 극소로 바뀌고, 12초 뒤 max|θ−π| 도 3.140 rad(전복)에서 0.027 rad(직립 유지)로 갈린다. a/l ≈ 0.458 을 넘기면 U_eff 는 골을 더 깊게 그리는데 진자는 다시 쓰러진다 — 유효 퍼텐셜은 a/l ≪ 1 근사다.

첫째 항은 익숙한 중력 퍼텐셜로 Θ=0\Theta=0이 골, Θ=π\Theta=\pi가 마루다. 둘째 항은 sin⁡2\sin^2이라 Θ=π/2\Theta = \pi/2에서 최대이고 연직 방향 양쪽(0,π0,\pi) 모두에서 0이다. 즉 흔들기는 수평 방향에 언덕을 쌓아 올린다. 그 언덕이 충분히 높으면 원래 마루였던 Θ=π\Theta=\pi가 양옆이 더 높아지면서 상대적으로 골짜기가 된다.

정량적으로는 Θ=π\Theta=\pi에서의 2계 도함수 부호다.

Ueff′′(π)=mgl(−1)+mA2ω22>0⟺  (Aω)2>2gl  U_{\text{eff}}''(\pi) = mgl\bigl(-1\bigr) + \frac{mA^2\omega^2}{2} > 0 \qquad\Longleftrightarrow\qquad \boxed{\;(A\omega)^2 > 2gl\;}

AωA\omega는 지지점의 속도 진폭이다. 조건이 진폭이나 진동수 하나가 아니라 그 곱으로 나온다는 것이 실험적으로 중요하다 — 진폭을 반으로 줄이면 진동수를 두 배로 올려야 본전이다. 흔히 인용되는 무차원 형태는 ε2(ω/ω0)2>2\varepsilon^2 (\omega/\omega_0)^2 > 2, 즉

ωω0>2ε\frac{\omega}{\omega_0} > \frac{\sqrt2}{\varepsilon}

ε=0.1\varepsilon = 0.1이면 고유진동수의 14배 이상으로 흔들어야 한다는 뜻이다. 길이 20 cm 진자(ω0≈7 \omega_0\approx7\,rad/s)에 진폭 2 cm면 대략 100 rad/s, 16 Hz 근방. 탁상에서 재현 가능한 숫자라는 점이 이 데모가 100년 넘게 살아남은 이유다.

조건을 넘으면 거꾸로 선 평형 주위의 미소 진동 각진동수도 바로 나온다.

ωeff2=A2ω22l2−gl\omega_{\text{eff}}^2 = \frac{A^2\omega^2}{2l^2} - \frac{g}{l}

영상에서 거꾸로 선 진자를 손가락으로 툭 밀면 느릿하게 되돌아오는데, 그 느릿한 진동수가 이것이다. 흔들기 진동수와는 완전히 다른 시간척도라 눈으로도 구분된다.

5. 평형점 지형이 통째로 바뀐다[편집]

퍼텐셜이 바뀌었으니 평형점 목록 자체가 달라진다. Ueff′=0U_{\text{eff}}'=0은

sin⁡Θ(mgl+mA2ω22cos⁡Θ)=0\sin\Theta\Bigl(mgl + \frac{mA^2\omega^2}{2}\cos\Theta\Bigr) = 0

이고, Θ=0,π\Theta = 0,\pi 외에

cos⁡Θ∗=−2glA2ω2\cos\Theta_* = -\frac{2gl}{A^2\omega^2}

라는 새 평형점이 등장한다. 우변의 절댓값이 1보다 작아야 하므로 존재 조건이 정확히 (Aω)2>2gl(A\omega)^2>2gl — 안정화 조건과 같은 순간이다. 이 Θ∗\Theta_*는 불안정(안장)이고, Θ=π\Theta=\pi의 흡인 유역 경계 역할을 한다.

즉 문턱을 넘는 순간 벌어지는 일은 이렇다. 불안정하던 Θ=π\Theta=\pi가 안정으로 바뀌면서, 그 대가로 불안정 평형 한 쌍이 π\pi에서 갈라져 나온다. 대칭성을 가진 분기 이론의 갈래질(pitchfork) 그 자체다. 구동을 더 세게 하면 cos⁡Θ∗→0\cos\Theta_*\to0, 즉 Θ∗→π/2\Theta_*\to\pi/2로 벌어지므로 거꾸로 선 상태의 유역은 문턱 바로 위에서 0이었다가 상반구 전체로 자란다. “안정해졌다”와 “웬만큼 밀어도 버틴다”는 별개 질문이고, 후자는 문턱보다 한참 위에서야 참이 된다.

6. 근사는 어디서 깨지는가[편집]

평균화는 공짜가 아니다. (Aω)2>2gl(A\omega)^2>2gl은 최저차 결과이고, 세 방향에서 교정을 받는다.

위쪽 진폭 한계. 거꾸로 선 위치 근방에서 선형화하면 (y=θ−πy = \theta-\pi)

y¨−gly+Aω2lcos⁡(ωt) y=0\ddot y - \frac{g}{l}y + \frac{A\omega^2}{l}\cos(\omega t)\,y = 0

인데, τ=ωt/2\tau=\omega t/2로 재척도하면 마티외 방정식 표준형의

a=−4glω2,∣q∣=2Al=2εa = -\frac{4g}{l\omega^2}, \qquad |q| = \frac{2A}{l} = 2\varepsilon

이 된다. 안정 영역은 인스-스트럿 선도의 첫 번째 흰 띠, 즉 a0(q)a_0(q)와 b1(q)b_1(q) 사이다. 아래쪽 경계 a>a0(q)≈−q2/2a > a_0(q) \approx -q^2/2가 방금의 (Aω)2>2gl(A\omega)^2>2gl이고 — 위쪽 경계 a<b1(q)a < b_1(q)가 평균화 이론이 아예 모르는 조건이다. 빠른 구동 극한(a→0−a\to0^-)에서 그 경계는 q≲0.908q\lesssim0.908이므로

ε=Al≲0.45\varepsilon = \frac{A}{l} \lesssim 0.45

가 된다. 너무 세게 흔들면 다시 넘어진다. 이때의 불안정은 플로케 승수가 −1-1을 통과하는 주기배가라, 진자가 그냥 쓰러지는 것이 아니라 구동 진동수의 절반으로 좌우로 펄럭이기 시작한다(flutter mode). 실험에서 확인된 현상이며, 이 상한이 존재한다는 사실 자체가 유효 퍼텐셜 그림의 한계를 보여 준다.

아래쪽 경계의 고차 보정. 같은 선도에서 a0(q)=−q2/2+7q4/128−⋯a_0(q) = -q^2/2 + 7q^4/128 - \cdots까지 쓰면 조건이

(Aω)2(1−716ε2)>2gl(A\omega)^2\Bigl(1 - \frac{7}{16}\varepsilon^2\Bigr) > 2gl

로 조금 빡빡해진다. ε=0.2\varepsilon=0.2면 보정이 1.8 %, ε=0.4\varepsilon=0.4면 7 %다. 정밀 실험에서는 이 차이가 보인다.

감쇠. 실물에는 마찰이 있다. 감쇠 2γθ˙2\gamma\dot\theta가 붙으면 θ=e−γtu\theta = e^{-\gamma t}u 치환으로 무감쇠 마티외로 정확히 되돌아가고(마티외 방정식의 그 항등변환), 안정 영역이 aa축으로 γ2\gamma^2만큼 이동한다. 그런데 두 경계가 서로 다르게 반응한다는 것이 이 계의 재미있는 대목이다.

  • 아래쪽 경계 a0a_0는 안 움직인다. γ2\gamma^2만큼 밀려 들어간 만큼, θ\theta가 유계이려면 uu가 eγte^{\gamma t}보다 느리게만 자라면 되므로 허용치도 딱 그만큼 늘어난다. 두 효과가 최저차에서 정확히 상쇄된다. 확인은 q=0q=0에서 한 줄로 된다 — θ¨+2γθ˙+(g ⁣∗/l)θ=0\ddot\theta + 2\gamma\dot\theta + (g_{\!\ast}/l)\theta=0의 특성근이 −γ±γ2−g ⁣∗/l-\gamma\pm\sqrt{\gamma^2 - g_{\!\ast}/l}인데, 뒤집힌 평형(g∗<0g_\ast<0)이면 γ\gamma가 아무리 커도 한쪽 근이 항상 양수다. 마찰로는 넘어지는 것을 막지 못한다는, 생각해 보면 당연한 사실. 따라서 (Aω)2>2gl(A\omega)^2>2gl은 감쇠가 있어도 그대로 살아남는다.
  • 위쪽 경계 b1b_1은 올라간다. 이쪽은 파라메트릭 공진형(주기배가) 불안정이라 성장률이 감쇠를 이겨야 켜진다. 즉 마찰이 있으면 A/lA/l 상한이 0.45보다 높아지고, 안정 띠가 위로 넓어진다.

덤으로 감쇠는 진짜 점근 안정을 만들어 준다. 무감쇠 마티외의 “안정”은 유계 진동일 뿐 수렴이 아니라서, 이상적인 무마찰 카피차 진자는 한 번 밀면 영원히 흔들린다. 실험 영상에서 진자가 결국 똑바로 멎는 것은 마찰 덕이다.

7. 왜 중요한가 — 동역학적 안정화[편집]

카피차 진자가 물리학 교과서에 남은 이유는 진자 자체보다 “고주파 흔들기로 불안정한 평형을 안정화한다”는 처방이 도처에서 재사용되기 때문이다.

  • 파울 트랩과 이온 트랩. 라플라스 방정식 때문에 정전기장만으로는 하전입자를 3차원으로 가둘 수 없다 — 어느 방향으로든 안장점이 생긴다. 그 안장을 RF로 흔들어 가둔 것이 파울 트랩이고, 그 유효 퍼텐셜이 이 문서의 UeffU_{\text{eff}}와 문자 그대로 같은 폰데로모티브 퍼텐셜이다. 볼프강 파울이 1989년 노벨상 강연에서 카피차 진자 데모 장치를 들고 나온 것도 우연이 아니다.
  • 플라스마와 레이저 트랩. 비균질 고주파장 속 하전입자에 작용하는 폰데로모티브 힘, 광쌍극자 덫(optical dipole trap)의 갇힘, 자기 미러의 고주파 안정화가 전부 같은 계산이다.
  • 유체 계면. 무거운 유체가 가벼운 유체 위에 놓인 레일리-테일러 불안정 배치도 용기를 연직으로 진동시키면 안정화된다. 뒤집힌 진자의 연속체 판이며, 조건도 유사한 꼴이다. 반대로 진폭을 더 올리면 패러데이 파가 터진다 — 안정선도의 위쪽 경계를 넘은 것이다.
  • 진동 제어(vibrational control). 되먹임 없이 고주파 입력만으로 계를 안정화하는 제어 기법의 원형이다. 센서가 필요 없다는 것이 장점이고, 원하는 상태를 정확히 겨냥할 수 없다는 것이 단점이다. 되먹임으로 역진자를 세우는 것과는 철학이 정반대다.
  • 다관절 확장. 애치슨과 멀린이 1993년에 보인 것이 압권이다. 관절이 여러 개인 사슬 진자도 지지점을 흔들면 통째로 곧게 선다. 인도의 로프 마술처럼 보인다고 해서 별명이 붙었고, 조건은 각 관절의 유효 퍼텐셜이 동시에 골이 되는 것이다.3

8. 수치적으로 다루기[편집]

원식은 강성(stiff)이 아니지만 시간척도가 두 개라는 것이 함정을 만든다.

  • 스텝 크기. 빠른 주기 2π/ω2\pi/\omega를 최소 수십 스텝으로 나눠야 한다. 이걸 어기면 재앙이 아주 그럴듯한 얼굴로 온다 — 스텝이 빠른 진동을 언더샘플링하면 적분기가 강제항을 저주파 별칭으로 읽고, 문턱 아래에서도 진자가 얌전히 서 있는 가짜 안정화를 보여 준다. “돌려 보니 되던데요”가 가장 위험한 계 중 하나다. 구동 주기당 스텝 수를 두 배로 늘려도 답이 같은지 반드시 확인할 것.
  • 적분기. 감쇠 없는 경우는 심플렉틱 적분기(예: 속도 벨레) 쪽이 장시간 에너지 표류가 없어 유리하고, 감쇠가 있으면 룽게-쿠타법 4차로 충분하다.
  • 안정성 판정은 적분이 아니라 플로케로. “궤적이 안 넘어지네”는 유한 시간 관찰일 뿐이다. 선형화 방정식을 한 주기 [0,2π/ω][0,2\pi/\omega] 동안 두 초기조건 (1,0),(0,1)(1,0),(0,1)으로 적분해 모노드로미 행렬 MM을 만들고 ∣tr⁡M∣|\operatorname{tr}M|을 재는 것이 정확한 판정이다. 무감쇠면 det⁡M=1\det M=1이라 ∣tr⁡M∣<2|\operatorname{tr}M|<2가 곧 안정이다(플로케 이론).
  • 검증 루틴. 평균화가 맞는지 보려면 완전 비선형 방정식의 궤적에서 2π/ω2\pi/\omega 이동 평균을 뽑아 Θ(t)\Theta(t)를 얻고, 유효 퍼텐셜 방정식의 해와 겹쳐 그린다. 문턱 근처에서 두 곡선이 갈라지기 시작하는 ε\varepsilon이 곧 근사의 유효 한계다. 미소 진동 진동수 ωeff\omega_{\text{eff}}를 FFT로 재서 예측식과 맞춰 보는 것도 좋은 정량 검사다.

9. 관련 문서[편집]

10. Footnotes[편집]

  1. A. Stephenson, On induced stability, Phil. Mag. 15 (1908) 233. 카피차의 1951년 논문(Dynamic stability of a pendulum when its point of suspension vibrates)이 이름을 가져간 것은 그가 유효 퍼텐셜이라는 설명 도구를 세웠기 때문이다. 스티븐슨은 “이 해가 유계다”를 보였고 카피차는 “왜 그런지”를 그림 한 장으로 말했다. 물리학사에서 이름은 대체로 후자에게 간다. ↩

  2. 부호가 헷갈리기 쉬운 자리다. UeffU_{\text{eff}}를 ml2ml^2로 나눈 것이 각가속도의 퍼텐셜이므로 Θ¨=−Ueff′/(ml2)\ddot\Theta = -U_{\text{eff}}'/(ml^2)이고, ddΘ(A2ω24sin⁡2Θ)=A2ω24⋅2sin⁡Θcos⁡Θ\tfrac{d}{d\Theta}\bigl(\tfrac{A^2\omega^2}{4}\sin^2\Theta\bigr) = \tfrac{A^2\omega^2}{4}\cdot2\sin\Theta\cos\Theta를 l2l^2로 나누면 본문의 A2ω22l2sin⁡Θcos⁡Θ\tfrac{A^2\omega^2}{2l^2}\sin\Theta\cos\Theta가 나온다. 계수 2와 4가 어디서 왔는지만 놓치지 않으면 된다. ↩

  3. D. J. Acheson & T. Mullin, Upside-down pendulums, Nature 366 (1993) 215. 애치슨은 임의의 NN-링크 진자가 뒤집힌 상태로 안정화될 수 있음을 보였고, 멀린이 실험으로 확인했다. 영상으로 보면 스파게티가 접시 위에 서 있는 것처럼 보이는데, 물리학과 학부 데모 중 “이게 왜 되는지 설명해 보라”는 질문으로 대학원생을 가장 오래 붙잡아 두는 부류에 속한다. ↩