Exercício Resolvido de Oscilações Harmônicas
publicidade   



Determine a equação de movimento e o período de oscilações para um pêndulo físico, formado por um corpo de forma qualquer e massa m onde é conhecida a posição do centro de massa, no regime de pequenas oscilações.


Dado do problema:
  • Massa do corpo:    m.
Esquema do problema:

No pêndulo atua a seguintes força (Figura 1-A):
  • P: força peso.
Adotamos um sistema de referência em coordenada cilindricas (Figura 1-B), onde er, eθ e ez são os vetores unitários nas direções r, θ e z.

Figura 1

O vetor r localiza o Centro de Massa do corpo a uma distância d em relação ao ponto de fixação.

Solução

Aplicando a 2.ª Lei de Newton para o movimento de rotação
\[ \begin{gather} \bbox[#99CCFF,10px] {\mathbf{N}=I\mathbf{\alpha}} \tag{I} \end{gather} \]
  • N é o torque da força que atua no corpo;
  • I é o momento de inércia do corpo;
  • α é a aceleração angular do corpo.
O torque é dado por
\[ \begin{gather} \bbox[#99CCFF,10px] {\mathbf{N}=\mathbf{r}\times{\mathbf{F}}} \tag{II} \end{gather} \]
a única força que atua no corpo é a força peso P e |r| = d é a distância do ponto de fixação até o Centro de Massa (C.M.) do pêndulo.
Escrevendo a aceleração angular como
\[ \begin{gather} \bbox[#99CCFF,10px] {\mathbf{\alpha}=\frac{d^{2}\mathbf{\theta}}{dt^{2}}} \tag{III} \end{gather} \]
substituindo as expressões (II) e (III) na expressão (I)
\[ \begin{gather} \mathbf{r}\times{\mathbf{P}}=I\frac{d^{2}\mathbf{\theta}}{dt^{2}} \tag{IV} \end{gather} \]
onde

\( {\mathbf{P}}_{T}=-P\operatorname{sen}\theta\;{\mathbf{e}}_{\theta} \)
\( {\mathbf{P}}_{N}=P\cos \theta\;{\mathbf{e}}_{r} \)
\( \mathbf{r}=d\;{\mathbf{e}}_{r} \)
\( \mathbf{\theta}=\theta\;{\mathbf{e}}_{z} \)

Observação: Algumas pessoas acham dífícil entender que o vetor deslocamento angular aponte na direção ez perpendicularmente ao plano de rotação. Quando um corpo se desloca de uma posição r para uma posição r' temos um deslocamento s ao longo da trajetória, o deslocamento angular θ está contido no plano, mas o vetor deslocamento angular θ aponta perpendicularmente ao plano, isto preserva o produtor vetorial (Figura 2). O vetor deslocamento angular indica que o corpo está em rotação, o seu módulo indica o deslocamento angular (escalar) e o sentido do vetor indica o sentido de rotação do corpo, se o vetor θ é positivo o produto vetorial indica que o corpo está se deslocando no sentido horário, se for negativo o corpo está se deslocando no sentido anti-horário.

Figura 2

Este é o mesmo argumento usado para a velocidade angular, o vetor velocidade v é tangente à trajetória, mas o vetor velocidade angular ω é perpendicular à trajetória.

\[ \begin{gather} \mathbf{P}={\mathbf{P}}_{T}+{\mathbf{P}}_{N}\\ \mathbf{P}=-P\operatorname{sen}\theta\;{\mathbf{e}}_{\theta}+P\cos \theta\;{\mathbf{e}}_{r} \end{gather} \]
O produto vetorial   \( \mathbf{r}\times {\mathbf{P}} \)  será
\[ \begin{gather} \mathbf{r}\times{\mathbf{P}}= \left| \begin{matrix} {\mathbf{e}}_{r}& {\mathbf{e}}_{\theta} & {\mathbf{e}}_{z}\\ d & 0 & 0\\ P\cos \theta &-P\operatorname{sen}\theta & 0 \end{matrix} \right|=\\[5pt] =[0.0-0.(-P\operatorname{sen}\theta)]{\mathbf{e}}_{r}-[d.0-0.P\cos\theta]{\mathbf{e}}_{\theta}+[-dP\operatorname{sen}\theta-0.P\cos \theta]{\mathbf{e}}_{z}\\[5pt] \mathbf{r}\times{\mathbf{P}}=-dP\operatorname{sen}\theta\;{\mathbf{e}}_{z} \tag{V} \end{gather} \]
a força peso é dada por
\[ \begin{gather} \bbox[#99CCFF,10px] {P=mg} \tag{VI} \end{gather} \]
substituindo o valor de θ e as expressões (V) e (VI) na expressão (IV)
\[ \begin{gather} -dmg\operatorname{sen}\theta\;{\mathbf{e}}_{z}=I\frac{d^{2}\theta}{dt^{2}}\;{\mathbf{e}}_{z}\\[5pt] \frac{-{dmg}}{I}\operatorname{sen}\theta=\frac{d^{2}\theta}{dt^{2}} \end{gather} \]
a equação só possui componentes na direção ez, e fazendo a seguinte definição   \( \omega ^{2}=\frac{mgd}{I} \)   e escrevendo   \( \frac{d^{2}\theta}{dt^{2}}=\ddot{\theta} \)
\[ \begin{gather} \ddot{\theta}=-\omega ^{2}\operatorname{sen}\theta\\ \ddot{\theta}+\omega ^{2}\operatorname{sen}\theta =0 \end{gather} \]
como estamos trabalhando em um regime de pequenas oscilações podemos expandir a função sen θ em uma série de Taylor.

Expansão em série de Taylor de sen θ
\[ \bbox[#99CCFF,10px] {f(x)=\sum _{n=0}^{\infty}{\frac{f^{n}(a)}{n!}(x-a)^{n}}} \]
fazendo a expansão em torno do ponto de equilíbrio com a = 0, para os 6 primeiros termos da série, temos

\( \displaystyle \frac{f^{0}(0)}{0!}\theta ^{0}=\frac{\operatorname{sen}0}{1}.1=0 \)

Observação: \( f^{0} \) NÃO significa a função f elevada a potência zero, significa a derivada de ordem zero da função f, ou seja, a própria função calculada no pondo a.

\( \displaystyle \frac{f^{\text{I}}(0)}{1!}\theta ^{1}=\frac{\cos 0}{1}\theta =\theta \)

\( \displaystyle \frac{f^{\text{II}}(0)}{2!}\theta^{2}=\frac{-\operatorname{sen}0}{2.1}\theta ^{2}=0 \)

\( \displaystyle \frac{f^{\text{III}}(0)}{3!}\theta ^{3}=\frac{-\cos 0}{3.2.1}\theta^{3}=-{\frac{\theta ^{3}}{6}} \)

\[ \displaystyle \frac{f^{\text{III}}(0)}{3!}\theta ^{3}=\frac{-\cos 0}{3.2.1}\theta^{3}=-{\frac{\theta ^{3}}{6}} \]

\( \displaystyle \frac{f^{\text{IV}}(0)}{4!}\theta^{4}=\frac{-(-\operatorname{sen}0)}{4.3.2.1}\theta ^{4}=0 \)

\[ \displaystyle \frac{f^{\text{IV}}(0)}{4!}\theta^{4}=\frac{-(-\operatorname{sen}0)}{4.3.2.1}\theta ^{4}=0 \]

\( \displaystyle \frac{f^{\text{V}}(0)}{5!}\theta ^{5}=\frac{\cos 0}{5.4.3.2.1}\theta^{5}=\frac{\theta ^{5}}{120} \)

\[ \displaystyle \frac{f^{\text{V}}(0)}{5!}\theta ^{5}=\frac{\cos 0}{5.4.3.2.1}\theta^{5}=\frac{\theta ^{5}}{120} \]

A função seno pode ser representada pela seguinte série de potências
\[ \operatorname{sen}\theta =\theta -\frac{\theta ^{3}}{6}+\frac{\theta^{5}}{120}-... \]
Como estamos considerando θ um ângulo pequeno podemos fazer a aproximação
\[ \operatorname{sen}\theta \approx \theta \]
e desprezamos os termos de potências maiores.
Para um ângulo de \( 10°=\frac{\pi}{18}=0,1745 \), temos \( \operatorname{sen}\frac{\pi}{18}=0,1736 \), a aproximação representa um erro de 0,5%.
\[ \ddot{\theta}+\omega ^{2}\theta =0 \]
Solução da equação diferencial    \( \displaystyle \ddot{\theta}+\omega^{2}\theta =0 \)

A solução é do tipo exponencial, calculando suas derivadas
\[ \begin{array}{l} \theta =\operatorname{e}^{\lambda t} \\ \dot{\theta}=\lambda \operatorname{e}^{\lambda t} \\ \ddot{\theta}=\lambda^{2}\operatorname{e}^{\lambda t} \end{array} \]
substituindo na equação
\[ \begin{gather} \lambda ^{2}\operatorname{e}^{\lambda t}+\omega^{2}\operatorname{e}^{\lambda t}=0\\[5pt] \lambda^{2}+\omega^{2}=0\\[5pt] \lambda ^{2}=-\omega^{2}\\[5pt] \lambda =\pm i\sqrt{\omega^{2}}\\[5pt] \lambda =\pm i\omega \end{gather} \]
a solução é da seguinte forma, onde C1 e C2 são constantes
\[ \theta (t)=C_{1}\operatorname{e}^{i\omega t}+C_{2}\operatorname{e}^{-i\omega t} \]
usando a fórmula de Euler   \( \operatorname{e}^{ix}=\cos x+i\operatorname{sen}x \)
\[ \begin{gather} \theta (t)=C_{1}\left(\cos \omega t+i\operatorname{sen}\omega t\right)+C_{2}\left(\cos \omega t-i\operatorname{sen}\omega t\right)\\[5pt] \theta(t)=\left(C_{1}+C_{2}\right)\cos \omega t-i\left(C_{2}-C_{1}\right)\operatorname{sen}\omega t \end{gather} \]
definindo as seguintes constantes
\[ A=C_{1}+C_{2}\quad ,\quad B=i(C_{2}-C_{1}) \]
\[ \theta (t)=A\cos \omega t+B\operatorname{sen}\omega t \]
fazendo as seguintes definições
\[ \begin{array}{l} \cos \phi=\dfrac{A}{\sqrt{A^{2}+B^{2}}}\\[5pt] \operatorname{sen}\phi=\dfrac{B}{\sqrt{A^{2}+B^{2}}}\\[5pt] \theta_{0}=\sqrt{A^{2}+B^{2}} \end{array} \]
substituindo na equação
\[ \begin{gather} \theta (t)=\left(A\cos \omega t-B\operatorname{sen}\omega t\right)\frac{\sqrt{A^{2}+B^{2}}}{\sqrt{A^{2}+B^{2}}}\\[5pt] \theta(t)=\sqrt{A^{2}+B^{2}}\left(\frac{A}{\sqrt{A^{2}+B^{2}}}\cos \omega t-\frac{B}{\sqrt{A^{2}+B^{2}}}\operatorname{sen}\omega t\right)\\[5pt] \theta (t)=\theta_{0}\left(\cos \phi \cos \omega t-\operatorname{sen}\phi \operatorname{sen}\omega t\right) \end{gather} \]

A equação de movimento será
\[ \bbox[#FFCCCC,10px] {\theta (t)=\theta_{0}\cos \left(\omega t+\phi \right)} \]
O período de oscilações é dado por
\[ \bbox[#99CCFF,10px] {T=\frac{2\pi }{\omega}} \]
substituindo a definição de ω0 feita acima
\[ T=\frac{2\pi }{\sqrt{\dfrac{mgd}{I}}} \]
\[ \bbox[#FFCCCC,10px] {T=2\pi \sqrt{\frac{I}{mgd}}} \]

Observação: O momento de inérica I que aparece neste problema não é o momento calculado em relação ao centro de massa do corpo, mas o momento de inércia em relação ao ponto de fixação do corpo em torno do qual ele oscila.
publicidade   

Licença Creative Commons
Fisicaexe - Exercícios Resolvidos de Física de Elcio Brandani Mondadori está licenciado com uma Licença Creative Commons - Atribuição-NãoComercial-Compartilha Igual 4.0 Internacional .