Moment of inertia of a cube along the diagonal.

[carllacan]

Homework Statement

Calculate the moment of inertia of a cube which rotates along an axis along its diagonal.

Homework Equations

Moment of inertia definition: I = \int \rho (\vec{r}) \vec{r} ^2 dV
Angular velocity vector; \vec{\omega}=\omega (\frac{1}{\sqrt{3}}, \frac{1}{\sqrt{3}}, \frac{1}{\sqrt{3}})

The Attempt at a Solution

My biggest problem is on finding the position vector for the volume element. My try at it goes like this:
<br /> \vec{r} = \vec{p}- \vec{p}·\vec{\omega} \frac{\vec{\omega}}{\omega}<br />
Where \vec{p} is the point of interest. With this we calculate the difference between the point's vector and a vector with the direction of \vec{\omega} and magnitude equal to that of the component of \vec{p} along \vec{\omega}. This should be the shortest vector from the axis of rotation to the point. Am I right?

I think, however, that there is an easier way to solve this. Isn't there any shortcut similar to the Steiner Theorem, only for rotates axes?

Thank you.

 

[SteamKing]

The moment of inertia values for a body form a tensor. Instead of trying to calculate the moment of inertia directly using the diagonal as an axis, calculate the moments using an axis aligned with the edges of the cube passing through the centroid and rotate them to the diagonal.

This discussion is relevant:

http://www.eng.auburn.edu/~marghitu/MECH2110/staticsC4.pdf

See Eq. 4.8 on page 5

 

[carllacan]

The moment of inertia values for a body form a tensor. Instead of trying to calculate the moment of inertia directly using the diagonal as an axis, calculate the moments using an axis aligned with the edges of the cube passing through the centroid and rotate them to the diagonal.

This discussion is relevant:

http://www.eng.auburn.edu/~marghitu/MECH2110/staticsC4.pdf

See Eq. 4.8 on page 5

The section you refer to talks about translation of axis. Isn't this problem more about a rotation of axes? Anyway, thanks for the link. I used the results on section 4.2 for a general axis and I obtained this:
I = \frac{\rho}{\sqrt{3}}(-\frac{1}{2} ab^2c^2+\frac{1}{2} a^2bc^2-\frac{1}{2} a^2b^2c + \frac{2}{3}b^3)
Does it make any sense? It doesn't for me. Why is the horizontal dimension given appear on its own? Shouldn't it be symmetrical?

 

[SteamKing]

Eq. 4.8 (and Sec. 4.2) come before Sec. 4.3, Translation of Coordinate Axes.

As to the validity of your calculations, I don't know. Why don't you show the full work up, starting with your calculation of the inertia tensor before rotation?

 

[carllacan]

Ah, right, I misread you answer.

My calculations look like this:
<br /> I = \int dV| \vec{u} x \vec{r}|^2 = \int dV \left | \begin{array}{ccc} <br /> \vec{i} &amp; \vec{j} &amp; \vec{k} \\<br /> \frac{1}{\sqrt{3}} &amp; \frac{1}{\sqrt{3}} &amp; \frac{1}{\sqrt{3}} \\<br /> x &amp; y &amp; z \\ \end{array} \right | = \int dV \frac{1}{\sqrt{3}} (z - y, -z + x, y - x)^2 = <br />
<br /> \int dV \frac{1}{\sqrt{3}} (z^2 -2zy + y^2 -z^2 + 2zx -x^2 + y^2 - 2yx + x^2) = \frac{1}{\sqrt{3}} \int dV ( -2zy + y^2 + 2zx + y^2 - 2yx) =

\frac{1}{\sqrt{3}} ( -\frac{ab^2c^2}{2} + \frac{a^2bc^2}{2} - -\frac{a^2b^2c}{2}+ \frac{2}{3}b^3 )<br />

 

[SteamKing]

It's not clear what the length of the edge of the cube is.

You also have 2 i's in your integral determinant.

If the edge length of the cube is say 'a', then the inertia tensor of the cube of uniform density about its centroidal axes will have entries only on the main diagonal (the products of inertia are identically zero because of symmetry). These values on the main diagonal will be

I_{xx} = I_{yy} = I_{zz} = \frac{m}{12}*(a^{2}+a^{2}) = \frac{m}{6}*a^{2}

m = ρ a^{3}, where ρ is the density

and off diagonal

I_{xy} = I_{yz} = I_{xz} = 0

You can apply Eq. 4.8 to this tensor without evaluating the integral.

 

[AlephZero]

Something went wrong in post #5, because the last term in your answer has the different dimensions (length^3) compared with the others (length^5).

And if this is a cube, why are there three different dimensions a b and c?

The really neat way to do this is not Eq 4.8, but look further down at sections 4.4 and 4.5. For a cube, the "ellipsoid of inertia" is not a general ellipsoid. The three axes are the sane length, so it is a sphere. That means the inertia is the same about any line through the centroid. The inertia tensor about any set of axes through the centroid is a diagonal matrix with the three equal terms on the diagonal.