Calculating the Inertia Tensor of cone with uniform density

[Wavefunction]

Homework Statement

Calculate the moments of inertia I_1, I_2, andI_3 for a homogeneous cone of mass M whose height
is h and whose base has a radius R. Choose the x_3 axis along the axis of symmetry of the cone.
Choose the origin at the apex of the cone, and calculate the elements of the inertia tensor. Then
make a transformation such that the center of mass of the cone becomes the origin, and find the
principal moments of inertia. Hint: Steiner's parallel-axis theorem might be useful for the
second part of the problem.

Homework Equations

J_{ij}=\int_{M}\rho\delta_{ij}\sum_{k}x_{k}^2-x_{i}x_{j}dM

I_{ij}=J_{ij}-M\sum_{k}a_{k}^2-a_{i}a_{j}

The Attempt at a Solution

The region I'll be integrating over will be: 0\leq\theta\leq 2\pi, 0\leq r\leq R, and 0\leq z \leq rcot(\alpha) such that x_1=rcos(\theta), x_2 = rsin(\theta), and x_3 = z. Also \alpha is the angle between the z-axis and the boundary of the cone.

Case 1) \delta_{ij} = 1\rightarrow J_{ii}=\rho\int_V\sum_{k}x_{k}^2-x_{i}^2dV

J_{ii}=\rho\int_{0}^{2\pi}\int_{0}^{R}\int_{0}^{rcot(\alpha)}[\sum_{k}x_{k}^2-x_{i}^2]dzrdrd\theta

i=1 \rightarrow J_{11} = \rho\int_{0}^{2\pi}\int_{0}^{R}\int_{0}^{rcot(\alpha)}[x_{2}^2+x_{3}^2]dzrdrd\theta

J_{11} = \rho\int_{0}^{2\pi}\int_{0}^{R}\int_{0}^{rcot(\alpha)}[(rsin(\theta))^2+(z)^2]dzrdrd\theta = MR^2\frac{2cot^2(\alpha)+3}{10cot(\alpha)}

i=2 \rightarrow J_{22} = \rho\int_{0}^{2\pi}\int_{0}^{R}\int_{0}^{rcot(\alpha)}[x_{1}^2+x_{3}^2]dzrdrd\theta

J_{22} = \rho\int_{0}^{2\pi}\int_{0}^{R}\int_{0}^{rcot(\alpha)}[(rcos(\theta))^2+(z)^2]dzrdrd\theta = MR^2\frac{2cot^2(\alpha)+3}{10cot(\alpha)}

i=3 \rightarrow J_{33} = \rho\int_{0}^{2\pi}\int_{0}^{R}\int_{0}^{rcot(\alpha)}[x_{1}^2+x_{2}^2]dzrdrd\theta

J_{33} = \rho\int_{0}^{2\pi}\int_{0}^{R}\int_{0}^{rcot(\alpha)}[(rcos(\theta))^2+(rsin(\theta))^2]dzrdrd\theta = \rho\int_{0}^{2\pi}\int_{0}^{R}\int_{0}^{rcot(\alpha)}[(r)^2]dzrdrd\theta = \frac{3MR^2}{5}

Case 2 \delta_{ij}=0 \rightarrow J_{ij}=\rho\int_V[-x_{i}x_{j}]dV = \rho\int_{0}^{2\pi}\int_{0}^{R}\int_{0}^{rcot(\alpha)}[-x_{i}x_{j}]dzrdrd\theta

i=1,j=2 \rightarrow J_{12} = \rho\int_{0}^{2\pi}\int_{0}^{R}\int_{0}^{rcot(\alpha)}[-r^2cos(\theta)sin(\theta]dzrdrd\theta = 0

i=1,j=3 \rightarrow J_{13} = \rho\int_{0}^{2\pi}\int_{0}^{R}\int_{0}^{rcot(\alpha)}[-rcos(\theta)z]dzrdrd\theta = 0

i=2,j=3 \rightarrow J_{23} = \rho\int_{0}^{2\pi}\int_{0}^{R}\int_{0}^{rcot(\alpha)}[-rsin(\theta)z]dzrdrd\theta = 0

so the Inertia tensor is then:

<br /> J = MR^2\begin{pmatrix}<br /> \frac{2cot^2(\alpha)+3}{10cot(\alpha)} &amp; 0 &amp; 0 \\<br /> 0 &amp; \frac{2cot^2(\alpha)+3}{10cot(\alpha)} &amp; 0 \\<br /> 0 &amp; 0 &amp; \frac{3}{5} <br /> \end{pmatrix}<br />

Now using Steiner's parallel axis thereom:

I_{ij}=J_{ij}-M\delta_{ij}\sum_{k}a_{k}^2-a_{i}a_{j} I'll need the components of \vec{a} which is the vector that points from the center of mass to the origin. So I'll need to find the center of mass of the cone using: R_{i}=\frac{1}{M}\int_{M}[x_{i}]dM. By the symmetry of the cone I already know the center of mass will have have coordinates (0,0,b). So I just have to find the value for x_3. First the total mass of the cone M is:
M=\rho\int_{0}^{2\pi}\int_{0}^{R}\int_{0}^{rcot(\alpha)}[1]dzrdrd\theta = \frac{2R^3\pi}{3}cot(\alpha)

\int_{M}[x_{i}]dM = \rho\int_{0}^{2\pi}\int_{0}^{R}\int_{0}^{rcot(\alpha)}[z]dzrdrd\theta = \frac{R^4\pi}{4}cot^2(\alpha) so the center of mass is located at (0,0,\frac{3}{8}Rcot(\alpha))

Case 1) \delta_{ij}=1 \rightarrow I_{ii} = J_{ii} -M\sum_{k}a_{k}^2-a_{i}^2

i=1 \rightarrow I_{11}=J_{11} -M[a_{2}^2+a_{3}^2] \rightarrow I_{11}=J_{11}-M\frac{9}{64}R^2cot^2(\alpha)

i=2 \rightarrow I{22}=J_{22} -M[a_{1}^2+a_{3}^2] \rightarrow I_{22}=J_{22}-M\frac{9}{64}R^2cot^2(\alpha)

i=3 \rightarrow I{33}=J_{33} -M[a_{1}^2+a_{2}^2] \rightarrow I_{33}=J_{33}-M[0]^2

Case 2) \delta_{ij}=0 \rightarrow I_{ij} = J_{ij} -M[-a_{i}a_{j}] . Since the only non-zero component is a_3 M[a_{i}a_{j}]=0 so I_{ij}=J_{ij} forj≠i. Then the inertia tensor for the center of mass is:

<br /> I = MR^2\begin{pmatrix}<br /> \frac{2cot^2(\alpha)+3}{10cot(\alpha)}-\frac{9}{64}cot^2(\alpha) &amp; 0 &amp; 0 \\<br /> 0 &amp; \frac{2cot^2(\alpha)+3}{10cot(\alpha)}-\frac{9}{64}cot^2(\alpha) &amp; 0 \\<br /> 0 &amp; 0 &amp; \frac{3}{5} <br /> \end{pmatrix}<br />

Finally the inertia tensor has two degenerate eigenvalues (principal moments of inertia) and one distinct eigenvalue:

I_1=I_2 = [\frac{2cot^2(\alpha)+3}{10cot(\alpha)}-\frac{9}{64}cot^2(\alpha)]MR^2

I_3 = \frac{3}{5}MR^2

My result checks out with everything taught in lecture as far as getting the two degenerate eigenvalues and one distinct one due to the symmetry of the cone, but this is the first type of problem like this that I have done so I would appreciate any help ya'll can give me (: thanks.

 

[SteamKing]

Boy, that looks like a lot of work!

There are some things you could have done to make the calculations easier. Writing out and evaluating volume integrals at each step is a bit much. An easier approach, I think, would be to stack circular disks of varying radius along the z-axis. Once you calculate the MOI of a disk of radius r and thickness dz, that eliminates much of your calculation. You can use the parallel axis theorem to refer the moments of each disk to the apex of the cone.

I think instead of expressing the elements of the inertia tensor in terms of the angle α, you should use the given geometric properties of the cone, like the radius of the base, r, and the height, h. That way, you can check your results against tables of moments of inertia:

http://en.wikipedia.org/wiki/List_of_moments_of_inertia#List_of_3D_inertia_tensors

 

[Wavefunction]

Boy, that looks like a lot of work!

There are some things you could have done to make the calculations easier. Writing out and evaluating volume integrals at each step is a bit much. An easier approach, I think, would be to stack circular disks of varying radius along the z-axis. Once you calculate the MOI of a disk of radius r and thickness dz, that eliminates much of your calculation. You can use the parallel axis theorem to refer the moments of each disk to the apex of the cone.

I think instead of expressing the elements of the inertia tensor in terms of the angle α, you should use the given geometric properties of the cone, like the radius of the base, r, and the height, h. That way, you can check your results against tables of moments of inertia:

http://en.wikipedia.org/wiki/List_of_moments_of_inertia#List_of_3D_inertia_tensors

I understand the direction you're going with this and I like it, but I can't reconcile the following:
how can I integrate \delta_{ij}\sum_{k}x_{k}^2-x_ix_j only with respect to z. I understand how the volume will look: V=\int_{0}^{h}[\pi(\frac{Rz}{h})^2]dz, but as I said I can't reconcile that with the definition of the inertia tensor.

 

[SteamKing]

I understand the direction you're going with this and I like it, but I can't reconcile the following:
how can I integrate \delta_{ij}\sum_{k}x_{k}^2-x_ix_j only with respect to z. I understand how the volume will look: V=\int_{0}^{h}[\pi(\frac{Rz}{h})^2]dz, but as I said I can't reconcile that with the definition of the inertia tensor.

While the notation defining the components of the inertia tensor is remarkably compact, some degree of clarity is lost when trying to compute the inertia tensor for certain geometric figures, like those having symmetry, which qualities can be exploited to reduce the overall amount of calculation.

Unless you are some sort of mathematical savant, it takes a certain amount of 'decrypting' the notation \delta_{ij}\sum_{k}x_{k}^2-x_ix_j in order to find simplifications. Sometimes, it pays to do a little analysis of the problem before charging headlong into evaluating various integral expressions.

 

[Wavefunction]

While the notation defining the components of the inertia tensor is remarkably compact, some degree of clarity is lost when trying to compute the inertia tensor for certain geometric figures, like those having symmetry, which qualities can be exploited to reduce the overall amount of calculation.

Unless you are some sort of mathematical savant, it takes a certain amount of 'decrypting' the notation \delta_{ij}\sum_{k}x_{k}^2-x_ix_j in order to find simplifications. Sometimes, it pays to do a little analysis of the problem before charging headlong into evaluating various integral expressions.

Well I guess that's why he gives us the problems he gives us right? which have been helping a lot along with all the help that I get here, but with this problem I think we are suppose to crank through the integrals because we might be given something that does not have nice symmetry like a cone, but I do see I shouldn't have defined \alpha I just redid J_{11} and I got J_{11} = \frac{3MR^3}{10}[\frac{1}{R}+\frac{2h^2}{R^3}]. So I think that when I calculated everything out the first time I was a bit careless when plugging into the inertia tensor.

 

[SteamKing]

Well, usually you are doing these problems for some known geometric shape, which is, after all, what most things are composed of. If you get a truly random object, which cannot be decomposed into simpler pieces, that would be remarkable itself.

 

[Wavefunction]

Well, usually you are doing these problems for some known geometric shape, which is, after all, what most things are composed of. If you get a truly random object, which cannot be decomposed into simpler pieces, that would be remarkable itself.

Okay well let's go back to what you were saying about the disk method. As I said, I understand how to get the volume with such a method and I did some thinking about it and this is what I came up with:

J_{11} = \rho\int_{V}\sum_{k}x_k^2-x_i^2dV \rightarrow J_{11} = \rho\int_{V}[x_1^2+x_2^2+x_3^2-x_1^2]dV the volume of one thin disk is: dV=\pi[r(z)]^2dz. Then I get J_{11}=\rho\int_{0}^{h}[x_2^2+x_3^2]\pi[r(z)]^2dz \rightarrow J_{11}=\rho\int_{0}^{h}[x_2^2+z^2]\pi[r(z)]^2dz So I'm just kind of stuck on how to parameterize x_2.