1. The problem statement, all variables and given/known data
a) Write down the Lagrangian and show that the equations of motion are:
[itex] m\ddot{θ_1}+k[2θ_1-θ_2-θ_3] = 0 [/itex]
[itex] m\ddot{θ_2}+k[2θ_2-θ_1-θ_3] = 0 [/itex]
[itex] m\ddot{θ_3}+k[2θ_3-θ_1-θ_2] = 0 [/itex]
b) Show that the mode in which [itex] θ_1 = θ_2 = θ_3 [/itex] corresponds to constant total angular momentum [itex] \vec{L} [/itex].
c) Assume that the total angular momentum is zero and that [itex] θ_1 = θ_2 = θ_3 = 0 [/itex]. Find two degenerate oscillatory modes and their frequency. Hint: you might want to express the three coupled equations as a single equation for the vector [itex]\vec{θ}[/itex], then assume that the eigenvectors of the tensor appearing in this equation have a time dependence [itex]e^{iωt}[/itex] and find the values of [itex]ω[/itex] for each eigenvector.
2. Relevant equations
eq(1) [itex] \frac{∂L}{∂q_j}-\frac{d}{dt}\frac{∂L}{∂\dot{q_j}}=0 [/itex]
3. The attempt at a solution
Part A
Step 1) The Lagrangian:
[itex] T = \frac{mR^2}{2}[\dot{θ_1}^2+\dot{θ_2}^2+\dot{θ_2}^2] [/itex]
[itex] U = \frac{kR^2}{2}[θ_1^2+θ_2^2+θ_1^2-θ_1θ_2-θ_1θ_3-θ_2θ_3] [/itex]
[itex] L = T-U = \frac{mR^2}{2}[\dot{θ_1}^2+\dot{θ_2}^2+\dot{θ_2}^2] - \frac{kR^2}{2}[θ_1^2+θ_2^2+θ_1^2-θ_1θ_2-θ_1θ_3-θ_2θ_3] [/itex]
[itex] \frac{∂L}{∂θ_1} = -kR^2[2θ_1-θ_2-θ_3] [/itex]
[itex] \frac{d}{dt}\frac{∂L}{∂\dot{θ_1}} = \frac{d}{dt}[mR^2\dot{θ_1}] = mR^2\ddot{θ_1} [/itex]
[itex] -kR^2[2θ_1-θ_2-θ_3]-mR^2\ddot{θ_1}=0 [/itex]
eq(2) [itex] m\ddot{θ_1}+k[2θ_1-θ_2-θ_3] [/itex]
Similarly for [itex] θ_2 [/itex] and [itex] θ_3 [/itex]
eq(3) [itex] m\ddot{θ_2}+k[2θ_2-θ_1-θ_3] [/itex]
eq(4) [itex] m\ddot{θ_3}+k[2θ_3-θ_1-θ_2] [/itex]
Part B
Combining eqs(2,3,and 4) into a vector equation:
[itex] \ddot{\vec{θ}}+ω_0^2\begin{pmatrix} 2&-1&-1\\-1&2&-1\\-1&-1&2\end{pmatrix}\vec{θ}=\vec{0} [/itex]
Let the above matrix be called [itex] \mathbf A [/itex]
Now when [itex] θ_1= θ_2= θ_3 [/itex] the matrix vector product [itex] \mathbf A [/itex][itex]\vec{θ}=\vec{0}[/itex]
This in turn yields [itex] \ddot{\vec{θ}} = \vec{0} →\dot{\vec{θ}}=\vec{C} [/itex] where [itex] \vec{C} [/itex] is a constant vector.
Part C
First I'll need the eigenvalues and the resulting eigenvectors of [itex] \mathbf A [/itex] which without too much trouble are:
[itex] λ_1=0 [/itex] and a degenerate eigenvalue [itex] λ_2 = λ_3 = 3 [/itex]
Which corresponds to the eigenvectors which have the time dependence [itex] e^{iωt} [/itex]:
[itex] \vec{e_1} = (1,1,1)e^{iωt}, \vec{e_2} = (-1,0,1)e^{iωt}[/itex], and [itex] \vec{e_3} = (-1,1,0)e^{iωt} [/itex]
When [itex] \mathbf A [/itex] acts on [itex]\vec{θ} [/itex] it produces two equations:
eq(5) [itex]\ddot{\vec{θ}}+(0)ω_0^2\vec{θ}=0[/itex]
eq(6) [itex] \ddot{\vec{θ}}+(3)ω_0^2\vec{θ}=0[/itex]
When I plug in eigenvectors 2 and 3 into eq(6) I get that [itex] ω^2 = 3ω_0^2 [/itex].
Physically, I think this scenario would correspond to one particle being held still while the other two oscillate back and forth with the same angular speed but in the opposite angular direction. Its been forever since I have dealt with systems of DE's and I seem to remember needing to construct another vector when a matrix had a repeated eigenvalue (like this one does) in order to form a complete solution space for the system. Also the scenario I give satifies that total angular momentum is zero, but it doesn't satisfy the second condition that [itex] θ_1=θ_2=θ_3=0 [/itex] which is what really kind of has me worried. Thank you for any help in advance!(:
a) Write down the Lagrangian and show that the equations of motion are:
[itex] m\ddot{θ_1}+k[2θ_1-θ_2-θ_3] = 0 [/itex]
[itex] m\ddot{θ_2}+k[2θ_2-θ_1-θ_3] = 0 [/itex]
[itex] m\ddot{θ_3}+k[2θ_3-θ_1-θ_2] = 0 [/itex]
b) Show that the mode in which [itex] θ_1 = θ_2 = θ_3 [/itex] corresponds to constant total angular momentum [itex] \vec{L} [/itex].
c) Assume that the total angular momentum is zero and that [itex] θ_1 = θ_2 = θ_3 = 0 [/itex]. Find two degenerate oscillatory modes and their frequency. Hint: you might want to express the three coupled equations as a single equation for the vector [itex]\vec{θ}[/itex], then assume that the eigenvectors of the tensor appearing in this equation have a time dependence [itex]e^{iωt}[/itex] and find the values of [itex]ω[/itex] for each eigenvector.
2. Relevant equations
eq(1) [itex] \frac{∂L}{∂q_j}-\frac{d}{dt}\frac{∂L}{∂\dot{q_j}}=0 [/itex]
3. The attempt at a solution
Part A
Step 1) The Lagrangian:
[itex] T = \frac{mR^2}{2}[\dot{θ_1}^2+\dot{θ_2}^2+\dot{θ_2}^2] [/itex]
[itex] U = \frac{kR^2}{2}[θ_1^2+θ_2^2+θ_1^2-θ_1θ_2-θ_1θ_3-θ_2θ_3] [/itex]
[itex] L = T-U = \frac{mR^2}{2}[\dot{θ_1}^2+\dot{θ_2}^2+\dot{θ_2}^2] - \frac{kR^2}{2}[θ_1^2+θ_2^2+θ_1^2-θ_1θ_2-θ_1θ_3-θ_2θ_3] [/itex]
[itex] \frac{∂L}{∂θ_1} = -kR^2[2θ_1-θ_2-θ_3] [/itex]
[itex] \frac{d}{dt}\frac{∂L}{∂\dot{θ_1}} = \frac{d}{dt}[mR^2\dot{θ_1}] = mR^2\ddot{θ_1} [/itex]
[itex] -kR^2[2θ_1-θ_2-θ_3]-mR^2\ddot{θ_1}=0 [/itex]
eq(2) [itex] m\ddot{θ_1}+k[2θ_1-θ_2-θ_3] [/itex]
Similarly for [itex] θ_2 [/itex] and [itex] θ_3 [/itex]
eq(3) [itex] m\ddot{θ_2}+k[2θ_2-θ_1-θ_3] [/itex]
eq(4) [itex] m\ddot{θ_3}+k[2θ_3-θ_1-θ_2] [/itex]
Part B
Combining eqs(2,3,and 4) into a vector equation:
[itex] \ddot{\vec{θ}}+ω_0^2\begin{pmatrix} 2&-1&-1\\-1&2&-1\\-1&-1&2\end{pmatrix}\vec{θ}=\vec{0} [/itex]
Let the above matrix be called [itex] \mathbf A [/itex]
Now when [itex] θ_1= θ_2= θ_3 [/itex] the matrix vector product [itex] \mathbf A [/itex][itex]\vec{θ}=\vec{0}[/itex]
This in turn yields [itex] \ddot{\vec{θ}} = \vec{0} →\dot{\vec{θ}}=\vec{C} [/itex] where [itex] \vec{C} [/itex] is a constant vector.
Part C
First I'll need the eigenvalues and the resulting eigenvectors of [itex] \mathbf A [/itex] which without too much trouble are:
[itex] λ_1=0 [/itex] and a degenerate eigenvalue [itex] λ_2 = λ_3 = 3 [/itex]
Which corresponds to the eigenvectors which have the time dependence [itex] e^{iωt} [/itex]:
[itex] \vec{e_1} = (1,1,1)e^{iωt}, \vec{e_2} = (-1,0,1)e^{iωt}[/itex], and [itex] \vec{e_3} = (-1,1,0)e^{iωt} [/itex]
When [itex] \mathbf A [/itex] acts on [itex]\vec{θ} [/itex] it produces two equations:
eq(5) [itex]\ddot{\vec{θ}}+(0)ω_0^2\vec{θ}=0[/itex]
eq(6) [itex] \ddot{\vec{θ}}+(3)ω_0^2\vec{θ}=0[/itex]
When I plug in eigenvectors 2 and 3 into eq(6) I get that [itex] ω^2 = 3ω_0^2 [/itex].
Physically, I think this scenario would correspond to one particle being held still while the other two oscillate back and forth with the same angular speed but in the opposite angular direction. Its been forever since I have dealt with systems of DE's and I seem to remember needing to construct another vector when a matrix had a repeated eigenvalue (like this one does) in order to form a complete solution space for the system. Also the scenario I give satifies that total angular momentum is zero, but it doesn't satisfy the second condition that [itex] θ_1=θ_2=θ_3=0 [/itex] which is what really kind of has me worried. Thank you for any help in advance!(:
0 commentaires:
Enregistrer un commentaire