조화 진동자

'''조화진동자(Harmonic oscillator, 調和振動子)'''는 평형점을 기준으로 물체의 변위에 비례한 복원력이 작용하게 되어 일정한 주기 운동을 하는 계를 말한다. 주로 용수철 진동자가 예로 많이 들어지며, RLC회로 등 조화 진동자와 동일한 양상을 보이는 고전역학을 벗어난 계도 존재한다.

2. 표현법 및 참고사항

이 문서 및 하위문서에서는 시간에 대한 물리량의 미분을 나타낼 때

$$ \displaystyle \dot{A} \equiv \frac{\mathrm{d}A}{\mathrm{d}t},\,\,\ddot{A} \equiv \frac{\mathrm{d}^{2}A}{\mathrm{d}t^{2}} $$[1]

[1] 좌변은 뉴턴식, 우변은 라이프니츠식 미분계수 표기다.

로 나타내져 있다. 여기서 $$A $$는 물리량, $$t $$는 시간이다.
또한, 특별한 말이 없는 이상 모든 마찰과 물체의 크기, 실이나 용수철의 질량은 무시할 수 있다고 가정한다.

3. 고전역학

3.1. 용수철 진동자

3.1.1. 단순 조화 진동

[image]
위 그림과 같이 평형점[2] $$\text{O} $$로 부터 변위 $$+A $$만큼 늘이고[3], 물체를 놓았을 때 물체의 운동이 어떻게 기술되는 지 논의해보자.
훅의 법칙에 따라 용수철이 평형점을 기준으로 $$x $$만큼의 변위로 변화되었을 때의 탄성력은

$$ \displaystyle F=-kx $$

[2] 용수철이 늘어나지 않은 지점을 말한다.[3] 마찬가지 조건에서 압축했을 때도 똑같이 기술된다.

이다. 참고로 마이너스는 운동 방향과 반대 방향의 힘을 나타내며, 복원력이다. 따라서 운동 방정식은

$$ \displaystyle m\ddot{x}=-kx\,\rightarrow \, \ddot{x}=-\frac{k}{m}x $$

이다. 여기서 $$k/m \equiv \omega^{2} $$으로 각진동수의 제곱으로 정의한다. 이렇게 하면, 미분 방정식은 쉽게 풀리며, 그 해는

$$ \displaystyle x(t)=C_{1}\sin(\omega t)+C_{2}\cos(\omega t) $$

이다. 이때, 초기 변위는 $$x(0)=A $$이고, 물체를 놓은 시점이므로 물체의 초기 속도은 $$\dot{x}(0)=0 $$을 만족해야 한다. 따라서 이 조건에 의해 sine 항은 해가 될 수 없고, $$\dot{x}(0)=0 $$, $$C_{2}=A $$를 얻는다. 따라서 용수철의 운동을 기술하는 변위 함수는

$$ \displaystyle x(t)=A\cos(\omega t) $$

로, $$-A \leq x \leq A $$ 사이를 진동하는 것을 알 수 있다.
그런데 위의 상황은 극히 특수한 상황이며, 실제론 관측하는 초기 시점을 뭘로 잡느냐에 따라 운동을 기술하는 변위 함수는 달라지게 된다. 위에서 밝혔듯 이 문제에선 sine 항과 cosine 항의 선형 결합으로 운동을 기술할 수 있고, 두 함수는 평행 이동 관계에 있기 때문에 일반적인 상황에서는 위상차 $$\phi $$를 도입해서,

$$ \displaystyle x(t)=A\sin(\omega t+\phi) $$

로 나타나게 되며, 여기서 $$A $$가 진폭이 되게 된다.
용수철 진자의 속력은 변위 함수를 한번 시간에 대해 미분함으로써 얻는다.

$$ \displaystyle \dot{x}(t)=A\omega\cos(\omega t+\phi) $$

가속도는 두 번 미분함으로써 얻는다.

$$ \displaystyle \ddot{x}(t)=-A\omega^{2}\sin(\omega t+\phi) $$

다음으로는 역학적 에너지가 보존되는지 살펴보자. 맨 위에 같은 상황에서 초기 역학적 에너지 $$E_{0} $$는 탄성에 의한 에너지밖에 없으므로

$$ \displaystyle E_{0}=\frac{1}{2}kA^{2} $$

$$\textrm{O} $$점으로 부터 변위가 $$x $$일 때, 물체는 탄성에 의한 역학적 에너지와 운동 에너지 두 개 모두 갖는다.

$$ \displaystyle E=\frac{1}{2}kx^{2}+\frac{1}{2}m{\dot{x}}^{2} $$

따라서

$$ \displaystyle E=\frac{1}{2}k\left[ A\sin(\omega t+\phi) \right]^{2}+\frac{1}{2}m\left[ A \omega \cos(\omega t+\phi) \right]^{2} $$

이때, 위에서 한 정의에 따라 $$m \omega^{2} =k $$가 되므로

$$ \displaystyle E=\frac{1}{2}k\left[ A\sin(\omega t+\phi) \right]^{2}+\frac{1}{2}k\left[ A\cos(\omega t+\phi) \right]^{2}=\frac{1}{2}kA^{2}=E_{0} $$

로 역학적 에너지는 보존됨을 알 수 있다.
위 식

$$ \displaystyle E=\frac{1}{2}kx^{2}+\frac{1}{2}m{\dot{x}}^{2} $$

에서 위상공간 즉, $$x \text{-} \dot{x} $$ 공간에서 개형은 적절히 처리하면, 타원임을 알 수 있다. 즉, 단순 조화 진동자의 위상도는 아래와 같이 나오게 된다.[4]
[image]

3.1.1.1. 원운동과의 관계

[image]
위의 그림에서 (a)와 같이 등속 원운동하고 있는 물체를 고려해보자. 이때, 원운동 평면 상에 평행하게 들어오는 평행광선을 고려하고, 이것을 스크린에 비춘다고 생각해보자. 그렇다면, 스크린에는 스크린 평면 상에 투사된 원운동의 자취가 나올 것인데, 이 자취의 위치를 (b)와 같이 $$\theta(t=0)=0$$이라 두고, 시간 $$t$$에 따른 위치를 나타내면, 사인곡선의 그래프가 그려지는데, 이것은 곧, 단순 조화 진동자의 시간에 따른 위치 그래프와 같다.
이상에서 단순 조화 진동자의 운동은 원운동을 위와 같이 투사시킨 자취의 운동이라 볼 수 있는 것이다.
이렇게 생각할 경우, 구심력은 곧 조화 진동자의 최고 지점에서의 복원력과 같게 되므로

$$\displaystyle mr\omega^{2}=kr $$

[4] $$x(0)=+A$$, $$\dot{x}(0)=0$$일 때.

따라서 조화 진동자의 각진동수가 다음과 같고, 위의 미분 방정식을 이용한 해석과 같음을 얻을 수 있다:

$$\displaystyle \omega=\sqrt{\frac{k}{m}} $$

3.1.2. 감쇠 조화 진동

[image]
이번엔 위의 경우에서 속도 $$\dot{x} $$에 비례하는 저항력 $$-b\dot{x} $$가 물체에 작용할 때, 기술되는 운동을 논의해보자. 이때, 운동 방정식은

$$ \displaystyle m\ddot{x}=-kx-b\dot{x}\,\rightarrow \, \ddot{x}+\frac{b}{m} \dot{x}+\frac{k}{m}x=0 $$

이고, 이 방정식은 2계 선형 방정식임에 따라 쉽게 풀린다. 이 방정식의 특성 방정식은

$$ \displaystyle D^2+\frac{b}{m}D+\frac{k}{m}=0 $$

이므로

$$ \displaystyle D=-\frac{b}{2m} \pm \sqrt{ \left( \frac{b}{2m} \right)^2 - \frac{k}{m } } $$

이때, $$b/2m \equiv \beta$$의 감쇠 계수로 정의하면

$$ \displaystyle D=-\beta \pm \sqrt{ \beta^2 - \omega^{2}} $$

이 됨에 따라 이상에서 이 방정식은 판별식 $$\displaystyle \beta^2 - \omega^{2} $$의 부호에 따라 다음과 같이 세 가지 형태의 해를 갖는다는 것을 알 수 있다.

'''형태'''	'''비고'''	'''경우'''
$$ \displaystyle \beta^2 > \omega^{2} $$	두 실근	ⓐ
$$ \displaystyle \beta^2 = \omega^{2} $$	중근	ⓑ
$$ \displaystyle \beta^2 < \omega^{2} $$	두 허근	ⓒ

우선 경우 ⓐ를 따져보도록 하자. 이 경우에 기술되는 변위 함수는

$$ \displaystyle x(t)=e^{-\beta t}[ C_{1} e^{ \sqrt{ \beta^2 - \omega^{2}} \, t } +C_{2} e^{ -\sqrt{ \beta^2 - \omega^{2}} \, t } ] $$

이다.
다음으로 경우 ⓑ를 따져보면, 그 해는 다음과 같다.

$$ \displaystyle x(t)= e^{-\beta t} [ C_{3}t+C_{4} ] $$

또한, 경우 ⓒ를 따져보면, 다음을 얻는다.

$$ \displaystyle x(t)= e^{-\beta t} [ C_{5}\sin{ ( \sqrt{ \beta^2 - \omega^{2}} \, t ) }+ C_{6}\cos{ ( \sqrt{ \beta^2 - \omega^{2}} \, t ) } ] $$

여기서 $$C_{1} \sim C_{6}$$는 각각 상수이며, 초기조건으로 결정할 수 있다.
경우 ⓐ~ⓒ에 대하여, 마찰 계수 $$b $$만을 변화시키고, $$x(0)=A $$, $$\dot{x}(0)=0 $$경우에 대해서 시뮬레이션 해보면, 아래와 같은 지수적으로 감소하는 개형이 나오게 된다.
[image]
특히 ⓒ에 대하여 변위 함수 $$x(t)$$는 아래와 같이 두 함수 $$e^{-\beta t}$$, $$-e^{-\beta t}$$ 사이에 위치하게 된다.
[image]
이때, ⓑ와 같이 평형점으로 돌아가는 감쇠를 '''임계 감쇠(Critical damping)''', ⓐ와 같이 ⓑ와 유사하지만 ⓑ에 비해 느리게 평형점으로 되돌아가는 감쇠를 '''과대 감쇠(Over damping)''', ⓒ와 같이 진동 형태를 유지하다, 평형점으로 되돌아가는 감쇠를 '''과소 감쇠(Under damping)'''라 한다. 이때, ⓐ, ⓑ는 계가 비진동성으로 응답하는 것을 알 수 있다.
이때, 아래와 같은 조건에 따른 감쇠 형태로 분류된다.

'''경우'''	'''감쇠'''	'''경우'''
$$ \displaystyle \beta>\omega $$	과대 감쇠	ⓐ
$$ \displaystyle \beta=\omega $$	임계 감쇠	ⓑ
$$ \displaystyle 0<\beta<\omega $$	과소 감쇠	ⓒ

감쇠 조화 진동자의 위상공간 즉, $$x \text{-} \dot{x} $$ 공간에서 그래프는 아래와 같은 개형을 가진다.[5]
[image]
위상 공간의 개형을 보면 역학적 에너지는 보존되지 않는다는 사실[6]을 알 수 있는데, 이것은 마찰력이라는 비보존력이 계에 작용하고 있음을 상기하면, 당연한 결과라는 사실을 알 수 있다.

3.1.3. 강제 조화 진동

가장 간단한 강제 조화 진동의 예인 윗 문단의 상황에서 외부 cosine형 구동력 $$F_{0}\cos{\omega't} $$이 주어진다고 하자. 이렇게 되면, 운동 방정식은 다음과 같이 세워지게 된다.

$$ \displaystyle m\ddot{x}=-kx-b\dot{x}+F_{0}\sin{\omega't}\,\rightarrow \, \ddot{x}+\frac{b}{m} \dot{x}+\frac{k}{m}x=\frac{F_{0}}{m}\cos{\omega't} $$

[5] 아래의 위상 공간에 대한 초기 조건은 위 그래프의 상황과 동일하다.[6] 역학적 에너지가 보존되는 단순 조화 진동자와 비교해보라.

이때, $$\sqrt{k/m} \equiv \omega $$, $$F_{0}/m \equiv G_{0} $$, 윗 문단에서 언급한 감쇠 계수 $$\beta $$를 도입하면, 위의 미분 방정식은

$$ \displaystyle \ddot{x}+2\beta \dot{x}+\omega^{2} x=G_{0}\cos{\omega't} $$

이 되며, 이 방정식은 특이해 $$x_{p}(t) $$와 정상해 $$x_{c}(t) $$의 선형 결합으로 이루어진다. 이것의 의미는 후술하기로 한다.
우선 특이해를 구하자. 오일러 공식을 이용하여,

$$ \displaystyle \ddot{X}+2\beta \dot{X}+\omega^{2} X=G_{0}e^{i \omega ' t} $$

으로 식을 바꾼 뒤, 예상되는 특이해 $$X_{p}(t)=C e^{i \omega ' t} $$를 대입하면,

$$ \displaystyle (-{\omega '}^{2}+2 \beta \omega ' i+\omega^{2})Ce^{i \omega ' t}=G_{0}e^{i \omega ' t} $$

이상에서

$$ \displaystyle C=\frac{G_{0}}{(\omega^{2}-{\omega '}^{2})+2{\beta}i \omega '}=\frac{ [ (\omega^{2}-{\omega '}^{2})-2{\beta}i \omega ' ]G_{0}}{(\omega^{2}-{\omega '}^{2})^{2}+4{\beta}^{2} {\omega '}^{2}} $$

이때, $$C $$를 극형식으로 나타냄으로써, 다음을 얻는다.

$$ \displaystyle C=\frac{G_{0}}{\sqrt{{(\omega^{2}-{\omega '}^{2})}^{2}+4{\beta}^{2} {\omega '}^{2} } }e^{-i \phi} \qquad \qquad \phi=\arctan{\left ( \frac{2{\beta} \omega '}{\omega^{2}-{\omega '}^{2}} \right )} $$

여기서 $$\phi $$는 구동력과 그 응답으로 나타나는 운동에 대한 위상차이다. 따라서 특이해는

$$ \displaystyle X_{p}(t)=\frac{G_{0}}{\sqrt{{(\omega^{2}-{\omega '}^{2})}^{2}+4{\beta}^{2} {\omega '}^{2} } }e^{i(\omega ' t-\phi)} $$

인데, 우리가 구하는 것은 물리현상이므로 실수부만 해로 취급하면, 본래의 미분방정식의 특이해는

$$ \displaystyle x_{p}(t)=\frac{G_{0}}{\sqrt{{(\omega^{2}-{\omega '}^{2})}^{2}+4{\beta}^{2} {\omega '}^{2} } }\cos{(\omega ' t-\phi)} $$

정상해 $$x_{c}(t) $$는 미분 방정식의 우변이 $$0 $$이 되는 해로써, 이것은 감쇠 진동에서 구했던 해와 같다. '''즉, 정상해는 진동계 자체의 진동 효과를 기술한다는 것을 알 수 있다.''' 또한, 정상해의 경우 윗 문단에서 보았듯, 감쇠항 $$\exp{(-\beta t)} $$이 있어 $$t \,\rightarrow \, \infty $$일 때, $$x_{c}(t) \, \rightarrow \, 0 $$이 됨에 따라 시간이 매우 지난 후에 계의 운동을 기술하는 변위 함수는

$$ \displaystyle x(t) \,\rightarrow \, x_{p}(t)=\frac{G_{0}}{\sqrt{{(\omega^{2}-{\omega '}^{2})}^{2}+4{\beta}^{2} {\omega '}^{2} } }\cos{(\omega ' t-\phi)} $$

가 된다. 이것은 곧 초반에는 정상해 항이 존재하기 때문에 계 자체의 진동과 외부 구동력에 의한 진동 두 효과가 동시에 나타나다($$x(t)=x_{p}(t)+x_{c}(t)$$) 시간이 많이 지나면, 결국 계는 외부 구동력의 효과만 남게 된다.($$x(t)\to x_{p}(t)$$)
아래의 그림은 위 과정을 시각화한 것이다. 아래와 같이 강제 조화 진동은 진동계 자체의 진동 효과인 $$x_{c}(t)$$와 외부 구동력에 의한 진동 효과 $$x_{p}(t)$$의 선형결합으로 주어진다는 것을 알 수 있다.
[image]
이번에는 이 운동의 위상차에 대한 것을 고려해보도록 하자. 강제 조화 진동자의 위상차는

$$ \displaystyle \phi=\arctan{\left ( \frac{2{\beta} \omega '}{\omega^{2}-{\omega '}^{2}} \right )} $$

임을 보았다. 따라서 이는 아래와 같이 요약할 수 있다.

$$\boldsymbol{\omega \gg \omega'}$$: $$\phi \to 0$$이므로 외부 구동력과 그 응답으로 나타나는 진동계의 운동의 위상은 거의 비슷하다.
$$\boldsymbol{\omega \approx \omega'}$$: $$\phi \approx \pi/2$$이므로 그 응답으로 나타나는 진동계의 운동의 위상은 외부 구동력의 위상의 $$-\pi/2$$만큼 차이난다.
$$\boldsymbol{\omega \ll \omega'}$$: $$\phi \to \pi$$이므로 그 응답으로 나타나는 진동계의 운동의 위상은 구동력의 위상에 반전되어 나타난다.($$-\pi$$의 위상차를 가진다.)

아래는 강제 조화 진동자의 외부 구동력의 진동수에 따른 위상차를 나타낸 그래프이다.
[image]
다음으로는 이 계의 공명에 대해 생각해보자. 시간이 매우 지난 뒤, 계는 진폭

$$ \displaystyle \frac{G_{0}}{\sqrt{{(\omega^{2}-{\omega '}^{2})}^{2}+4{\beta}^{2} {\omega '}^{2} } } \equiv A $$

로 진동하게 된다. 이것이 어떤 각진동수 $$\omega_{r} $$에서 최대가 될 조건은

$$ \displaystyle \left. \frac{dA}{d \omega '} \right|_{\omega '=\omega_{r}}=0 $$

을 만족해야 하므로

$$ \displaystyle \omega_{r}=\sqrt{{\omega}^{2}-{\beta}^{2}} $$

일 때, 공명이 일어나게 된다.
아래는 외부 구동력의 진동수에 따른 계의 진폭을 그래프로 나타낸 것이다.
[image]
흥미로운 것은 만약 계의 감쇠 계수가 매우 작아지면, 공명 진동수는 $$\omega $$로 수렴하게 되며, 이때, 진폭은 $$\infty $$로 발산한다. 즉, 진폭이 급격히 커져 진동계 자체가 파괴될 수 있다.
이와 관련해서, 물체의 기본 진동수로 접근하는 구동력을 물체에 가하면, 물체는 급격한 진동을 일으켜, 물체가 파괴될 수 있는 것을 간접적으로 알 수 있게 한다.[7]

3.1.3.1. 주기적이고 임의의 외부 구동력이 작용할 때

이 문단에서는 cosine형 구동력이 아닌 임의의 주기적인 구동력 $$F(t)$$가 주어지는 경우를 보고자한다. 이 경우 운동 방정식은

$$ \displaystyle m\ddot{x}=-kx-b\dot{x}+F(t)\,\rightarrow \, \ddot{x}+\frac{b}{m} \dot{x}+\frac{k}{m}x=\frac{F(t)}{m}$$

[7] 이것의 예로는 모 방송이 실험했던 '목소리로 와인잔 깨기' 등이 있다.

와 같이 세워질 것이다. 우리는 앞서 $$F(t)$$를 주기적인 구동력이라 가정했고, 이에 따라

$$ \displaystyle \begin{aligned} F(t)&=\sum_{n=0}^{\infty}a_{n} \sin{(n \omega' t)}+\sum_{n=0}^{\infty}b_{n} \cos{(n \omega' t)} \\ a_{n}&=\dfrac{\omega'}{\pi} \int_{0}^{2\pi/\omega'} F(t')\sin{(n \omega' t')}\,dt' \\ b_{n}&=\dfrac{\omega'}{\pi} \int_{0}^{2\pi/\omega'} F(t')\cos{(n \omega' t')}\,dt' \end{aligned}$$

의 푸리에 급수로 전개할 수 있다. 따라서 우변은 곧 무한한 개수의 sine형 구동력이 중첩되어있는 상황으로 해석할 수 있고, 해는

$$ \displaystyle x(t)=x_{c}(t)+\sum_{n=0}^{\infty} X_{n}(t)+\sum_{n=0}^{\infty}Y_{n}(t) $$

의 형태로 주어짐을 알 수 있다.[8] 여기서 $$X_{n}$$, $$Y_{n}$$은 각각 sine항, cosine항의 특이해이며, $$x_{c}$$는 정상해이다.
그런데 위의 결과를 살펴보면,

$$ \displaystyle \begin{aligned} X_{n}(t)&=A_{n} \cos{(n \omega' t-\phi_{X_{n}})}, \qquad &\begin{cases} A_{n}&= \dfrac{a_{n}}{m}\dfrac{1}{ \sqrt{ (\omega^{2}-n^2 \omega'^{2} )^{2} +4 \beta^2 n^2 \omega'^2 }}\\ \phi_{X_{n}}&=\arctan{\left( \dfrac{2 \beta n \omega'}{\omega^2-n^2 \omega'^2} \right )} \end{cases} \\ Y_{n}(t)&=B_{n} \cos{(n \omega' t-\phi_{Y_{n}})}, \qquad &\begin{cases} B_{n}&= \dfrac{b_{n}}{m}\dfrac{1}{ \sqrt{ (\omega^{2}-n^2 \omega'^{2} )^{2} +4 \beta^2 n^2 \omega'^2 }}\\ \phi_{Y_{n}}&=\arctan{\left( \dfrac{2 \beta n \omega'}{\omega^2-n^2 \omega'^2} \right )} \end{cases} \end{aligned} $$

[8] 이는 운동을 기술하는 방정식이 선형 미분방정식이기 때문이다.

으로 쓸 수 있다.[9] $$\phi$$와 관련된 것은 외부 구동력과 그 응답으로 나타나는 운동에 대한 위상차이다.
이상에서 우리는 이 경우

$$ \displaystyle x(t)=x_{c}(x)+\sum_{n=0}^{\infty} A_{n} \cos{(n \omega' t-\phi_{X_{n}})}+\sum_{n=0}^{\infty}B_{n} \cos{(n \omega' t-\phi_{Y_{n}})} $$

[9] 외부 구동력이 sine형일 때는 이 문서에서 다루지는 않았으나, cosine형과 동일한 형태의 특이해를 얻는다.

임을 알 수 있다.
이 결과를 이용하면, 외부 구동력이 사각파나 삼각파, 톱니파 형태인 경우에도 변위함수를 구할 수 있다.
아래는 $$F(t)$$가 아래와 같은 사각파 형태[10]이고, $$\omega=4\,\mathrm{rad/s}$$, $$\omega'=2^{-1}\pi\,\mathrm{rad/s}$$, $$m=1\,\mathrm{kg}$$, $$\beta=0.5$$, $$x_{c}(0)=0$$, $$\dot{x}_{c}(0)=0.1\,\mathrm{m/s}$$일 때를 시뮬레이션한 것이다.
[image]

3.1.4. 용수철 진자

[image]
용수철 진자는 용수철 진동자를 세운 버전이라 볼 수 있다.
위 그림의 (가)와 같이 용수철 상수가 $$k $$인 가벼운 용수철만을 매달아 놓는다면, 용수철은 자연 길이[11]를 유지한다. 그러나, (나)와 같이 질량 $$m $$의 물체를 조심스레 용수철에 연결하게 되면, 중력 때문에 용수철은 늘어나고, 탄성력과 중력이 평형을 이룰 때까지 늘어난다. 즉, 위 그림에서

$$ \displaystyle mg=k \alpha $$

[10] 푸리에 급수로 전개 시 $$\displaystyle \begin{aligned} F(t)=\frac{1}{2}+\sum_{n=1}^{\infty} \dfrac{2}{n \pi} \sin{\left( \frac{\pi n }{2} \right)} \cos{\left( \frac{\pi n t}{2} \right)} \end{aligned}$$를 얻을 수 있고, 이 시뮬레이션에서는 100항까지의 합만 사용되었다. 더 간단한 형태인 $$(u \circ \cos)(x)$$ (단, $$u(x)$$는 헤비사이드 계단 함수)도 있긴 하지만 이 함수가 불연속함수이므로 여기서는 쓰지 않았다.[11] 본래의 자신의 길이

이며, 실제 실험할 때는 용수철 상수를 모르기 때문에 이렇게 용수철에 물체를 매달아 늘어난 길이 $$\alpha $$를 측정하여,

$$ \displaystyle k=\frac{mg}{\alpha} $$

를 이용하여 용수철 상수를 구한다.
반동이 가해지지 않는 이상 (나) 상태에서 용수철 진자는 힘의 평형이 이루어져, 정지하게 된다. 따라서 용수철 진자의 평형점은 (나)에서 O이다. 그런데 (다)처럼 평형점을 기준으로 $$ \displaystyle x $$만큼 늘어났을 때, 물체의 운동은 어떻게 기술되는 지 알아보자.
물체가 받는 힘은 자신의 중력과 탄성력이다. 중력 방향을 양으로 두면, 따라서 운동 방정식은

$$ \displaystyle m \ddot{x}=-k(x+\alpha)+mg $$

이고, 이것을 다시 쓰면,

$$ \displaystyle \ddot{x}+\frac{k}{m}x=-\frac{k\alpha}{m} +g=0 $$

으로 용수철 진동자와 같게 기술되는 것을 알 수 있다. 따라서 용수철 진자는 주기

$$ \displaystyle T=2\pi \sqrt{\frac{m}{k}} $$

으로 진동하며, 초기에 점 O를 기준으로, $$ \displaystyle x_{0} $$만큼 늘렸다면,

$$ \displaystyle x(t)=x_{0}\sin{\left( \sqrt{\frac{k}{m}} t+\phi \right)} $$

가 된다. 여기서 $$\phi $$는 위상차이다.
이때, 물체의 속도와 가속도는 각각 다음과 같다.

$$ \displaystyle \dot{x}(t)=x_{0} \sqrt{\frac{k}{m}} \cos{\left( \sqrt{\frac{k}{m}} t+\phi \right)} $$

$$ \displaystyle \ddot{x}(t)=- \frac{k}{m}x_{0} \sin{\left( \sqrt{\frac{k}{m}} t+\phi \right)} $$

이번에는 역학적 에너지가 보존되는 지 확인해보자. (나) 상황에서 O의 높이를 $$ \displaystyle 0 $$이라 잡으면, O를 기준으로 $$ \displaystyle x_{0} $$만큼 늘렸을 때, 물체의 역학적 에너지는 용수철의 탄성 에너지와 중력 퍼텐셜 에너지의 합이다.

$$ \displaystyle E_{0}=\frac{1}{2}k(x_{0}+\alpha)^{2}-mgx_{0} $$

이때, (다) 상황에서 물체의 역학적 에너지는

$$ \displaystyle E=\frac{1}{2}k (x+\alpha)^{2}+\frac{1}{2}m {\dot{x}}^{2}-mgx $$

위에서 구한 변위 함수와 속도 함수를 대입해서 전개해보면, $$E=E_{0} $$의 결과를 얻으므로 용수철 진자에서도 역학적 에너지는 보존된다는 사실을 확인할 수 있다.

3.2. 진자

3.3. 관련 논의

3.3.1. 지구 관통 터널

[image]
위 그림과 같이 지구 중심 $$\text{O} $$를 지나가는 무지막지하게 뜨거운 직선형 터널을 뚫었고, 물체를 가만히 놓았을 때, 물체의 운동을 분석해보자. 문제를 간단히하기 위해 모든 마찰은 무시하고, 지구는 구형이며, 밀도 $$\rho $$는 균일하며[12], 터널의 직경 $$L $$은 지구 반경에 비해 무시 가능할 정도로 작다고 생각할 것이다. $$m $$은 물체의 질량이다.
지구 중심 $$\text{O} $$를 기준으로 하여, 윗 방향을 $$+\hat{\mathbf{x}} $$라 놓자. 이때, 물체가 변위 $$x $$에 위치할 때, 터널의 직경이 매우 작으므로 물체는 지구 내부에 위치한다고 생각할 수 있고, 힘은 물체를 지나고, 중심이 $$\text{O} $$인 지구의 부분 질량에 의한 만유인력이다.[13][14]

$$ \displaystyle \mathbf{F}=- G \frac{m}{x^2} \left( \rho \cdot \frac{4}{3} \pi x^{3} \right) \hat{\mathbf{x}}=- \frac{4G \pi \rho m}{3} {\mathbf{x}} $$

[12] 즉 맨틀, 외핵, 내핵 등이 없다고 가정하고[13] 그림에서 점선에 해당하는 지구의 질량에 의한 만유인력을 말한다.[14] 이것의 증명은 뉴턴의 구각 정리로 할 수 있으며, 자세한 것은 고전역학 책을 참조한다.

이 된다. $$G $$는 만유인력 상수이다.
따라서 물체의 운동 방정식은

$$ \displaystyle m \ddot{x}=- \frac{4G \pi \rho m}{3} x \,\rightarrow \, \ddot{x}+ \frac{4G \pi \rho}{3} x=0 $$

따라서 물체는 주기

$$ \displaystyle T=2\pi \sqrt{\frac{3}{4G \pi \rho }}=\sqrt{\frac{3 \pi}{G\rho}} $$

로 진동운동 함을 알 수 있다. 따라서 물체를 초기에 $$x_{0} \leq R $$인 지점에서 떨어뜨렸다면, $$-x_{0} \leq x \leq x_{0} $$에서 진동하게 된다.
따라서 위 결과로 지구의 한 표면에서 대척점의 표면까지 이동하는데 걸리는 시간은 위에서 구한 주기의 반이므로

$$ \displaystyle \frac{T}{2}=\frac{1}{2}\sqrt{\frac{3 \pi}{G\rho}} $$

가 되고, 각종 상수

$$G=6.67384 \times 10^{-11} \, \mathrm{N} \cdot \text{m}^{2} / \mathrm{kg}^{2} $$
$$\rho=5515 \,\mathrm{kg}/\mathrm{m}^{3} $$

를 대입하면,

$$ \displaystyle \frac{T}{2} \approx 42.2\,\mathrm{min} $$

가 됨을 알 수 있다.
더 나아가 터널을 비스듬히 파서, 중심을 지나지 않더라도 주기는 모두 같게 나온다. 증명은 위와 같은 방법으로 하면 되니, 관심있는 위키러들은 해보도록 하자.

3.3.2. 액체관 진동

[image]
그림과 같이 단면적이 $$A $$인 U자관에 밀도가 $$\rho $$인 액체가 길이 $$l $$만큼 차있는 것을 고려하자. 문제를 단순하게 생각하기 위해 모든 마찰과 액체는 점성이 없다[15]고 가정하자. 따라서 관에 든 액체의 질량은 밀도와 부피의 곱으로 쓸 수 있으므로

$$ \displaystyle m = \rho Al $$

[15] 점성을 넣는 순간 식에 나비에-스토크스 방정식이라는 일반해조차 아직 안 나온 골칫거리를 넣어야 하기 때문.

이때, 한쪽이 평형점 O를 기준으로 $$x $$ 만큼 눌러졌을 때, 작용하는 복원력은 $$2x $$만큼의 액체 기둥에 작용하는 중력이다. 따라서 운동 방정식은 다음과 같이 쓸 수 있다.

$$ \displaystyle \rho Al \ddot{x}=-\rho A g \cdot 2x $$

이것을 다시 쓰면,

$$ \displaystyle \ddot{x}+\frac{2g}{l} x=0 $$

이고, 이것은 명백히 주기

$$ \displaystyle T= \pi \sqrt{\frac{2l}{g}} $$

로 진동하는 용수철 진동자의 진동 방정식과 같다.
따라서 액체 기둥은 위와 같은 주기로 진동함을 알 수 있다.

4. 전자기학

4.1. 전자기 진동: ''RLC'' 직렬 회로

[image]
그림과 같이 저항값, 인덕턴스, 전기 용량이 각각 $$R $$, $$L $$, $$C $$인 저항, 인덕터, 축전기로 연결된 ''RLC'' 직렬 회로를 생각해보자.
우선, 저항이 연결되어 있지 않다고 가정하고, 회로에 키르히호프 법칙을 적용하면,

$$ \displaystyle L \frac{di}{dt}+\frac{q}{C}=0 $$

이때, $$i \equiv dq/dt $$임을 이용하고, 식을 다시 쓰면,

$$ \displaystyle \frac{d^2q}{dt^2}+\frac{1}{LC}q=0 $$

이것은 명백히 주기

$$ \displaystyle T=2\pi \sqrt{LC} $$

로 진동하는 조화 진동자의 방정식이다.
따라서 축전기에 저장된 전하는 시간에 따라 진동하는 것을 알 수 있고, 더 나아가 회로에 흐르는 전류 또한 진동함을 쉽게 알 수 있다.[16]
이번에는 저항을 연결했다고 하면,

$$ \displaystyle L \frac{di}{dt}+Ri+\frac{q}{C}=0 $$

[16] $$i \equiv dq/dt $$

가 되고, 위처럼 식을 다시 쓰면,

$$ \displaystyle \frac{d^2q}{dt^2}+\frac{R}{L}\frac{dq}{dt}+\frac{1}{LC}q=0 $$

이 되고, 이것은 용수철 진동자에서 '감쇠 조화 진동'을 기술하는 방정식임을 발견할 수 있고, 감쇠에 영향을 주는 항은 저항임을 쉽게 이해할 수 있다.
또한, 외부 전원을 연결했다면, 방정식은

$$ \displaystyle L \frac{di}{dt}+Ri+\frac{q}{C}=V(t) $$

로 주어지고, 이것이 각진동수 $$\omega $$의 cosine 형 교류라면, $$V(t)=V_{0}\cos{(\omega t)} $$이므로

$$ \displaystyle L \frac{di}{dt}+Ri+\frac{q}{C}=V_{0}\cos{(\omega t)} $$

이것을 다시 쓰면,

$$ \displaystyle \frac{d^2q}{dt^2}+\frac{R}{L}\frac{dq}{dt}+\frac{1}{LC} q=\frac{V_{0}}{L}\cos{(\omega t)} $$

이고, 이것은 명백히 위에서 봤던 용수철 진동자에서 '강제 조화 진동'을 기술하는 방정식임을 알 수 있다.
따라서 전자기 진동 또한, 고전역학적 조화 진동자와 동일한 계로 기술될 수 있음을 얻는다.