- #1
Wavefunction
- 99
- 4
Homework Statement
A yo-yo of mass [itex]M[/itex] is composed of two disks of radius [itex]R[/itex] separated by a shaft of radius [itex]r[/itex]. A
massless string is would on the shaft, and the loose end is held in the hand. Upon release the
yo-yo descends until the string is unwound. The string the begins to rewind, and the yo-yo
climbs. Find the string tension and acceleration of the yo-yo in descent and ascent. Neglect the
mass of the shaft and assume the shaft radius is sufficiently small so that the string is essentially
vertical.
Homework Equations
[itex]J_{ij}=\int_{V}[\rho(r)\delta_{ij}\sum_{k}x_{k}^2-x_ix_j]dV [/itex]
[itex]I_{ij}=J_{ij}-M[\delta_{ij}\sum_{k}a_{k}^2-a_ia_j] [/itex]
[itex]L_i=\sum_{j}I_{ij}\omega_{j}[/itex]
The Attempt at a Solution
Setup: First I'll define the [itex] x_3 [/itex] axis to be the axis of rotation of the yo-yo. Then [itex]\vec{\omega}[/itex] will have the components [itex]\begin{pmatrix}0\\0\\\omega_{3}\end{pmatrix}[/itex]. Now I need to find [itex]\mathbb{I}[/itex] in order to find [itex]\vec{L}[/itex]. Finally, because I can neglect the mass of the small shaft I will be calculating the inertia tensor of a solid cylinder.
Calculation: Using [itex]J_{ij}=\int_{V}[\rho(r)\delta_{ij}\sum_{k}x_{k}^2-x_ix_j]dV [/itex] I'll find the inertia tensor at the bottom of the cylinder. So I'll need the volume of the region in question:
[itex] 0\leq z\leq a, 0\leq r\leq R, 0\leq \theta\leq 2\pi [/itex]
Case 1) [itex] \delta_{ij}=1 \Rightarrow J_{ii}=\int_{0}^{a}\int_{0}^{2\pi}\int_{0}^{R}[x_1^2+x_2^2+x_3^2-x_i^2]rdrd\theta dz[/itex]
where [itex] x_1 = rcos(\theta), x_2 = rsin(\theta), x_3 =z [/itex]
[itex] i=1; \rightarrow J_{11}=\rho\int_{0}^{a}\int_{0}^{2\pi}\int_{0}^{R}[(rsin(\theta))^2+(z)^2]rdrd\theta dz = \rho\int_{0}^{a}\int_{0}^{2\pi}[\frac{R^4sin^2(\theta)}{4}+\frac{R^2z^2}{2}]d\theta dz = \rho\int_{0}^{a}[\frac{\pi R^4}{4}+2\pi \frac{R^2z^2}{2}]dz [/itex]
[itex]= \rho[\frac{a\pi R^4}{4}+\frac{\pi R^2a^3}{3}] \rightarrow
M\frac{\frac{a\pi R^4}{4}+\frac{\pi R^2a^3}{3}}{a\pi R^2} = \frac{MR^2}{4}+\frac{Ma^2}{3}[/itex]
[itex] i=2; [/itex] the only that changes in this calculation is that [itex] sin(\theta) [/itex] becomes [itex] cos(\theta) [/itex]. Since [itex] \int_{0}^{2\pi}sin^2(\theta)d\theta = \int_{0}^{2\pi}cos^2(\theta)d\theta \Rightarrow J_{11}=J_{22} [/itex]
[itex] i=3; \rightarrow J_{11}=\rho\int_{0}^{a}\int_{0}^{2\pi}\int_{0}^{R}[(r)^2]rdrd\theta dz = \rho\int_{0}^{a}\int_{0}^{2\pi}[\frac{R^4}{4}]d\theta dz = \rho\int_{0}^{a}[2\pi\frac{R^4}{4}]dz [/itex]
[itex] = \rho[\frac{a\pi R^4}{2}] \rightarrow M\frac{\frac{a\pi R^4}{2}}{a\pi R^2} = \frac{MR^2}{2} [/itex]
Case 2)[itex] \delta_{ij}=0 \Rightarrow J_{ij}=\int_{0}^{a}\int_{0}^{2\pi}\int_{0}^{R}[-x_ix_j]rdrd\theta dz[/itex]
[itex] i=1,j=2 \rightarrow \int_{0}^{a}\int_{0}^{2\pi}\int_{0}^{R}[-r^2cos(\theta)sin(\theta)]rdrd\theta dz = \frac{1}{2}\int_{0}^{2\pi}sin(2\theta)d\theta[\int_{0}^{a}\int_{0}^{R}[-r^3]dr dz] = 0 [/itex]
[itex] i=1,j=3 \rightarrow \int_{0}^{a}\int_{0}^{2\pi}\int_{0}^{R}[-rcos(\theta)z]rdrd\theta dz = \int_{0}^{2\pi}cos(\theta)d\theta[\int_{0}^{a}\int_{0}^{R}[-zr^2]dr dz] = 0 [/itex]
[itex]i=2,j=3 \rightarrow \int_{0}^{a}\int_{0}^{2\pi}\int_{0}^{R}[-rsin(\theta)z]rdrd\theta dz = \int_{0}^{2\pi}sin(\theta)d\theta[\int_{0}^{a}\int_{0}^{R}[-zr^2]dr dz] = 0
[/itex]
So [itex] \mathbb{J}=\begin{pmatrix}\frac{MR^2}{4}+\frac{Ma^2}{3}&0&0\\0&\frac{MR^2}{4}+\frac{Ma^2}{3}&0\\0&0&\frac{MR^2}{2}\end{pmatrix}[/itex]
Now I want to find the inertia tensor about the center of mass [itex] \mathbb{I} [/itex]. Using [itex] I_{ij}=J_{ij}-M\delta_{ij}\sum_{k}a_{k}^2-a_ia_j [/itex] where [itex] \vec{a} = \begin{pmatrix}0\\0\\-\frac{a}{2}\end{pmatrix} [/itex]
Case 1) [itex] \delta_{ij}=1 \rightarrow I_{ij}=J_{ij}-M[a_1^2+a_2^2+a_3^2-a_i^2] [/itex]
[itex] i=1; I_{11}=J_{11}-M[a_2^2+a_3^2] \rightarrow I_{11} = \frac{MR^2}{4}+\frac{Ma^2}{3}-M\frac{a^2}{4} = \frac{MR^2}{4}+\frac{Ma^2}{12} = M\frac{3R^2+a^2}{12} [/itex]
[itex] i=2; I_{22} =J_{22}-M[a_1^2+a_3^2] \rightarrow I_{22} = \frac{MR^2}{4}+\frac{Ma^2}{3}-M\frac{a^2}{4} = \frac{MR^2}{4}+\frac{Ma^2}{12} = M\frac{3R^2+a^2}{12} [/itex]
[itex] i=3; I_{33} = J_{33}-M[a_1^2+a_2^2] \rightarrow I_{33}=J_{33} [/itex]
Case 2) [itex] \delta_{ij}=0 \rightarrow I_{ij}=J_{ij}+M[a_ia_j] [/itex] since [itex] j≠i \Rightarrow I_{ij} = J_{ij}[/itex] Since the only non-zero component of [itex]\vec{a}[/itex] is [itex] a_3 [/itex].
Therefore: [itex] \mathbb{I} = \begin{pmatrix}M\frac{3R^2+a^2}{12}&0&0\\0&M\frac{3R^2+a^2}{12}&0\\0&0&\frac{MR^2}{2}\end{pmatrix}[/itex]
Then: [itex] \vec{L}=\mathbb{I}\vec{\omega} \rightarrow \begin{pmatrix}L_1\\L_2\\L_3\end{pmatrix} = \begin{pmatrix}M\frac{3R^2+a^2}{12}&0&0\\0&M\frac{3R^2+a^2}{12}&0\\0&0&\frac{MR^2}{2}\end{pmatrix} \begin{pmatrix}0\\0\\\omega_3\end{pmatrix}[/itex]
Clearly the only non-zero component of [itex] \vec{L} [/itex] is [itex] L_3 = \frac{MR^2}{2}\omega_3 [/itex] So what is [itex] \omega_3 [/itex]? its the angular velocity [itex] \dot{\theta} [/itex] about the [itex] x_3 [/itex] axis so:
[itex] \vec{L} = \begin{pmatrix}0\\0\\\frac{MR^2}{2}\dot{\theta}\end{pmatrix} [/itex]. Also I know that [itex] \vec{\omega} [/itex] is parallel to [itex] \vec{L} [/itex] so [itex] \frac{\partial \vec{L}}{\partial t}_{inertial} =\frac{\partial \vec{L}}{\partial t}_{body}=\vec{\tau} [/itex] And [itex]\dot{\vec{L}} = \begin{pmatrix}0\\0\\\frac{MR^2}{2}\ddot{\theta}\end{pmatrix} = \begin{pmatrix}\tau_1\\\tau_2\\\tau_3\end{pmatrix} [/itex]
At this point I think I should look at the sum of the forces and torques on the yo-yo but I'm not quite sure, and I'd like to make sure the above work is correct before moving forward. Thanks for your help in advance.