결합 진동의 맥놀이와 에너지 흐름

질량의 운동과 에너지를 따로 관찰한다

약하게 연결한 두 진동자는 운동을 주고받는 것처럼 보인다. 한쪽 변위의 느린 포락선이 커지는 동안 다른 쪽은 작아진다. 맥놀이를 이해하는 유용한 그림이지만, 포락선 자체가 에너지 측정값은 아니다. 속도와 연결 용수철에 저장된 에너지까지 확인해야 한다.

이 단원은 정확히 풀 수 있는 두 질량의 계로 실제 궤적·포락선·국소 에너지·부호 있는 에너지 흐름을 비교한다. 익숙한 약한 결합 근사도 같은 초기조건과 명시한 관측 시간에서 정확한 해와 대조한다.

계와 초기 실험을 고정한다

같은 질량 $m>0$ 두 개에 동일한 고정점 강성 $k\ge0$와 결합 강성 $\kappa\ge0$를 사용한다. 단위는 kg, N/m, m, s, J다. 처음에는 첫 질량만 $A\ge0$만큼 옮기고 두 질량 모두 정지 상태에서 놓는다.

$$q_1(0)=A,\qquad q_2(0)=0,\qquad v_1(0)=v_2(0)=0.$$

무감쇠 운동방정식은

$$m\ddot q_1=-(k+\kappa)q_1+\kappa q_2,\qquad m\ddot q_2=\kappa q_1-(k+\kappa)q_2$$

다. 아래의 비교는 이 초기조건을 유지한다. 초기속도를 바꾸거나 한 정상모드만 가진 상태에서 시작하면 다른 실험이다.

대칭 운동과 반대칭 운동을 찾는다

같은 방향의 모양 $(1,1)$은 결합 용수철을 변형시키지 않으므로 각주파수는 $\omega_-=\sqrt{k/m}$다. 반대 방향의 모양 $(1,-1)$은 상대 변위를 두 배로 바꾸므로 $\omega_+=\sqrt{(k+2\kappa)/m}$다.

평균 주파수와 분할을

$$\bar\omega=\frac{\omega_++\omega_-}{2},\qquad \Delta\omega=\omega_+-\omega_-$$

로 정의한다. 결합이 양수이면 가까운 제곱근을 빼는 대신

$$\Delta\omega=\frac{2\kappa/m}{\omega_++\omega_-}$$

로 계산한다. 결합이0이면 분할은 정확히0이다. 두 강성도 모두0이면 $0/0$을 계산하지 말고 자유 운동의 극한을 직접 사용한다.

실제 변위 두 개를 재구성한다

초기 상태는 두 변위 모양을 같은 비율로 더한 것이므로

$$q_1(t)=\frac A2[\cos(\omega_-t)+\cos(\omega_+t)],\qquad q_2(t)=\frac A2[\cos(\omega_-t)-\cos(\omega_+t)]$$

이다. 더하거나 빼면 각 주파수로 움직이는 정상 좌표 하나씩을 되찾는다. 계수 $A/2$는 초기 변위로 정해지며 시각마다 새로 맞추는 값이 아니다.

둘째 질량은 초기 변위가0이어도 이후 움직일 수 있다. 정상모드는 각 질량에 대한 기여를 설명한다. 모드 하나를 질량 하나에 배정하는 설명은 일반적으로 맞지 않는다.

빠른 반송파와 느린 포락선을 드러낸다

삼각함수의 합차 공식에 따라

$$q_1=A\cos(\bar\omega t)\cos(\Delta\omega t/2),\qquad q_2=A\sin(\bar\omega t)\sin(\Delta\omega t/2)$$

이다. 부호 있는 포락선 인자는 $A\cos(\Delta\omega t/2)$와 $A\sin(\Delta\omega t/2)$다. 절댓값은 눈에 보이는 진폭의 경계를 나타내고, 빠른 반송파는 그 경계 안의 실제 위치를 정한다.

두 주파수가 가까우면 포락선이 천천히 바뀌는 동안 반송파가 여러 번 진동한다. 가깝지 않아도 식은 정확하지만, 한 인자를 천천히 변하는 진폭으로 해석하기는 어려워진다.

반송파뿐 아니라 포락선도 미분한다

$C=\cos(\bar\omega t)$, $S=\sin(\bar\omega t)$, $c=\cos(\Delta\omega t/2)$, $s=\sin(\Delta\omega t/2)$로 쓰면

$$v_1=-A[\bar\omega Sc+(\Delta\omega/2)Cs],\qquad v_2=A[\bar\omega Cs+(\Delta\omega/2)Sc]$$

이다. 포락선을 미분하는 항을 빼면 위치 그림은 그럴듯해도 운동에너지가 틀릴 수 있다. 포락선이0이 되는 순간에 특히 분명하다. 변위가0이어도 포락선이 변하는 속도는0이 아닐 수 있기 때문이다.

포락선의 시각과 전체 상태의 주기를 구분한다

$\Delta\omega>0$, $A>0$이면 첫 포락선 크기 교환은 $t_{\rm swap}=\pi/\Delta\omega$에서 일어난다. 크기는 $2\pi/\Delta\omega$마다 반복되지만 부호 있는 포락선은 $4\pi/\Delta\omega$마다 반복된다.

이 값이 전체 위치–속도 상태의 주기는 아니다. 전체 상태가 같은 시각에 돌아오려면 두 모드 주파수 사이에 적절한 유리수 비가 필요하다. 포락선 교환 때의 실제 반송파 위상을 따로 평가해야 한다.

결합이0이면 유한한 교환 주기가 없다. $A=0$이면 계의 고유 스펙트럼은 남아도 관측되는 교환 운동은 없다. 계산은 스펙트럼에서 정한 포락선 시각과, 0이 아닌 진폭으로 관측하는 교환을 구분한다.

연결 용수철도 에너지 장부에 넣는다

전체 역학적 에너지는

$$E=E_1+E_2+E_c,\quad E_j=\frac12mv_j^2+\frac12kq_j^2,\quad E_c=\frac12\kappa(q_2-q_1)^2$$

이다. 처음에는 $E=(k+\kappa)A^2/2$다. 첫 질량만 옮겼어도 에너지 일부는 연결 용수철에 저장되어 있다. 자체 에너지 $E_1,E_2$만 더하면 보존이 깨지는 것처럼 보일 수 있다.

정확한 모드별 에너지는 일정하다. 실제 질량의 에너지와 결합 에너지가 변하는 것은 같은 보존 총량을 다른 방식으로 나눈 결과다. 이 둘을 같은 에너지 이름표로 섞지 말자.

국소 에너지 관례를 정하고 흐름을 유도한다

총에너지를 두 국소 에너지로 나누기 위해 결합 에너지의 절반씩을 양쪽에 배분하자.

$$E_1^*=E_1+\frac12E_c,\qquad E_2^*=E_2+\frac12E_c.$$

실제 힘을 대입해 미분하면 부호 있는 흐름은

$$J_{1\to2}=\frac\kappa2(q_1-q_2)(v_1+v_2),\qquad \frac{dE_1^*}{dt}=-J_{1\to2},\quad \frac{dE_2^*}{dt}=J_{1\to2}$$

가 된다. 양수는 이 에너지 배분에서 둘째 쪽이 에너지를 얻는다는 뜻이다. 순간 흐름은 음수일 수도 있으며 첫째 쪽을 옮겨 놓은 실험의 초기에도 그렇다. 국소 배분은 명시한 관례이고, 물리 총에너지는 그 관례와 무관하다.

흐름을 끝점 에너지와 독립적으로 적분한다

이 초기조건에서 $q_1-q_2=A\cos(\omega_+t)$, $v_1+v_2=-A\omega_-\sin(\omega_-t)$이므로

$$J_{1\to2}=-\frac{\kappa A^2\omega_-}{2}\cos(\omega_+t)\sin(\omega_-t)$$

이다. 주파수0에서도 정칙인 적분

$$I(\nu,t)=\int_0^t\sin(\nu u)\,du =\frac{\nu t^2}{2}\operatorname{sinc}^2(\nu t/2)$$

을 정의하면 부호 있는 순전달량은

$$W_{1\to2}(t)=\frac{\kappa A^2\omega_-}{4} [I(\Delta\omega,t)-I(\omega_-+\omega_+,t)]$$

이다. 이 식을 끝점 에너지와 따로 계산한 뒤 $E_2^*(t)-E_2^*(0)-W_{1\to2}(t)$를 검사한다. 전달량을 에너지 차로 정의해 버리면 흐름 계산을 시험하지 않고 잔차를0으로 만들게 된다.

계산 예제: 포락선이 교환되어도 남는 에너지

SI 단위에서 $m=k=A=1$, $\kappa=1/16$ N/m로 두자. $\omega_-=1$, $\omega_+=\sqrt{1.125}\approx1.060660$ rad/s이고 첫 포락선 교환 시각은 약51.790039초다.

그때의 정확한 상태를 반올림하면 다음과 같다.

물리량	값
$q_1$	0 m
$q_2$	0.046223 m
$v_1$	0.030298 m/s
$v_2$	−1.029229 m/s
자체 $E_1$	0.000458974 J
결합 $E_c$	0.000066769 J
총 $E$	0.531250 J

둘째 포락선의 크기는1 m이지만 실제 변위는 반송파의 영점 근처에 있다. 첫 변위가0인 동안에도 속도는0이 아니다. 에너지 대부분이 둘째 쪽에 갔지만 포락선만으로 정확한 전체 전달을 결론 낼 수는 없다.

계산 예제: 정확한 변위 교환이 가능한 특수한 경우

$m=k=A=1$, $\kappa=1.5$로 두면 $\omega_-=1$, $\omega_+=2$ rad/s다. $t=\pi$초에서 $q=(0,-1)$ m이고 두 속도는0이다. 첫째 자체 에너지는 정확히0, 둘째 자체 에너지는0.5 J다.

결합 용수철에는 여전히0.75 J가 있어 총1.25 J다. 절반 배분 관례의 국소 에너지는 처음 $(0.875,0.375)$ J에서 $(0.375,0.875)$ J로 바뀌므로 순전달량은0.5 J다. 특정한 주파수 관계에서는 정확한 자체 에너지 교환이 가능하지만, 모든 포락선 교환의 일반 법칙은 아니다.

약한 결합 근사를 시간에 따라 비교한다

$k>0$, $\kappa/k$가 작을 때

$$\omega_+\approx\omega_0+\frac{\kappa}{m\omega_0},\qquad \omega_-=\omega_0=\sqrt{k/m}$$

로 전개할 수 있다. 이 근사 주파수를 사용해 같은 초기조건의 위치와 속도를 모두 재구성한다. 버린 주파수 보정은 시간이 지나면서 위상 오차로 누적된다. 결합비가 작다는 조건만으로 긴 시간의 위치 정확성까지 보장되지는 않는다.

비교에는 $\kappa/k$와 명시한 관측 창의 오차를 함께 보고한다. $k=0$이면 이 전개는 미정이지만 정확한 두 좌표의 해는 여전히 계산할 수 있다.

비결합·영모드·영진폭을 명시적으로 처리한다

$\kappa=0$이면 둘째 질량은 초기 평형에 남고, $k>0$에서 첫 질량만 자유 진동한다. 순전달량은0이고 유한한 교환 주기는 없다. $k=0$, $\kappa>0$이면 중심 좌표는 영모드이고 상대 좌표가 진동한다.

뒤의 경우 이 초기조건에서 절반 배분한 국소 에너지는 각각 $\kappa A^2/4$로 일정하고 흐름은 항등적으로0이다. 결합 에너지가 줄어드는 동안 두 질량의 운동에너지가 함께 늘 수 있다. 변위 포락선이 달라진다고 모든 국소 배분에서 에너지가 이동하는 것은 아니다.

$A=0$이면 모든 물리 에너지가0이다. 에너지 비율은 분모가0이므로 미정이다. 계의 주파수 정보가 남는다는 사실이 관측된 0이 아닌 운동을 만들어 내지는 않는다.

해설이 있는 연습문제 여섯 개

문제 1: 초기 상태 복원하기

두 코사인 식에 $t=0$을 대입하고 한 번 미분해 보라.

해설: 코사인이 모두1이므로 $(q_1,q_2)=(A,0)$이고 사인으로 된 미분항은 모두0이므로 초기속도도0이다. 모드 계수를 정한 뒤 네 초기조건을 확인하지 않으면 다른 실험을 풀 수 있다.

문제 2: 초기 에너지 세기

$m=k=A=1$, $\kappa=1/16$에서 처음 $E_1,E_2,E_c$는 얼마인가?

해설: 각각0.5, 0, 0.03125 J로 총0.53125 J다. 둘째 자체 에너지는0이지만 절반 배분 관례에서는 결합 에너지의 반인0.015625 J를 처음부터 둘째 쪽에 배정한다.

문제 3: 포락선 영점에서 미분하기

$\Delta\omega t/2=\pi/2$에서 $q_1=0$이면 $v_1=0$인가?

해설: 아니다. 반송파 미분항은 $c=0$과 함께 사라지지만 포락선 미분항 $v_1=-A(\Delta\omega/2)\cos(\bar\omega t)$가 남는다. 반송파 위상이나 진폭에 추가 조건이 있어야 속도까지0이 된다.

문제 4: 전달 방향 해석하기

양의 $k,\kappa,A$와 이 초기조건에서 작은 양의 시각의 $J_{1\to2}$ 부호는 무엇인가?

해설: 정확한 식을 전개하면 $J_{1\to2}\approx-\kappa A^2\omega_-^2t/2<0$이다. 절반 배분 관례에서 초기 흐름은 첫째 쪽을 향한다. 처음 변위를 준 위치만으로 이후 에너지 흐름의 부호를 모두 정할 수 없다.

문제 5: 고정점 용수철 없애기

$k=0$, $\kappa>0$에서 두 운동에너지가 변하는데도 흐름이0일 수 있는 이유는?

해설: 질량중심이 고정되어 속도의 합이0이므로 흐름 식이0이다. 연결 용수철에서 풀리는 에너지를 두 질량이 함께 받는다. 대칭 배분 관례의 두 국소 에너지는 서로 같다.

문제 6: 근사 성능 공정하게 평가하기

약한 결합 근사가2초 동안 정확해 보이면100초까지 검증한 것인가?

해설: 아니다. 주파수 오차는 위상 오차로 누적된다. 같은 시각의 정확한 위치·속도와 근삿값을 비교하고, 결합비와 관측 시간을 명시하며 전체 창의 오차를 읽어야 한다. 초기조건은 동일하게 유지한다.

그림을 읽고 AI 검증에 연결한다

실제 위치와 부호 있는 포락선, 자체·결합 에너지와 배분한 국소 에너지, 순간 흐름과 누적 적분, 정확한 궤적과 약한 결합 근사를 각각 분리해 표시한다. 모든 그림에 같은 선택 시각을 표시한다. 포락선 교환 시각이 관측 창 밖이면 그 사실을 유지하며 마지막 그림 점으로 옮기지 않는다.

교육용 계산은0.25 kg 간격의 질량과1/16 N/m 간격의 강성, 0–4 m 진폭을 사용한다. 관측 창은128초까지, 실제 한 시점 계산은512초까지 지원한다. 그리기 구간당 가장 빠른 위상 진행은0.125 rad 이하다. 포락선 주기를 수치 적분 단계 수와 혼동하지 않는다.

학습 시뮬레이터는 위치 손실이 작아도 속도·결합 에너지·흐름 부호를 잘못 재현할 수 있다. 학습에 없는 결합 강성, 긴 관측 시간, 영모드 경계를 시험하자. 국소 에너지의 물리 배분 관례를 보고해야 한다. 임의 잠재 에너지의 분할 하나에 맞았다는 사실만으로 실제 수송을 검증할 수 없다.

일차 자료와 모델의 범위

정상모드 중첩과 초기값 재구성은 MIT 8.03SC 교재3장 3.3.3–3.4.3절을 참고한다. 힘의 정식화는 MIT 8.03SC 강의4, 1–5쪽에도 나온다. 위의 에너지 배분·흐름 적분·반례는 명시한 두 질량의 초기 실험에서 유도했다. 감쇠·비선형 결합·다체 수송에는 추가 분석이 필요하다. 다음 원래 주제는 파동 방정식이다.