1
$\begingroup$

Suppose that $\Omega$ is a body of revolution of the function $y=f(x), a\le x \le b$ around the $x$-axis, where $f(x)>0$ is continuous.

How to compute the inertia tensor ${\bf{J}} _{\Omega}$?

After computing ${\bf{J}} _{\Omega}$, how to solve ${\bf{J}} _{\Omega} \dot w={\bf{J}} _{\Omega}w \times w$?

$\endgroup$
  • $\begingroup$ for the different components of $J$ there are definitions, see e.g. en.wikipedia.org/wiki/Moment_of_inertia#The_inertia_tensor, where you probably need to go to the integral form. as for the equation i don't think you have a $J$ on the RHS, but rather $Jw$. $\endgroup$ – Phoenix87 Feb 10 '15 at 16:31
  • 1
    $\begingroup$ Re After computing ${\bf{J}} _{\Omega}$, how to solve ${\bf{J}} _{\Omega} \dot w={\bf{J}} _{\Omega} \times w$? -- You can't. Look at the units. The left hand side has units of mass×length$^2$×time$^{-2}$, but the right hand side has units of mass×length$^2$×time$^{-1}$. The correct expression is ${\bf{J}} _{\Omega} \dot w=\boldsymbol\tau_\text{ext} - w \times ({\bf{J}} _{\Omega} w)$. $\endgroup$ – David Hammen Feb 10 '15 at 16:33
  • $\begingroup$ @DavidHammen Sorry, it is my mistake.:) $\endgroup$ – Sherry Feb 14 '15 at 19:57
1
$\begingroup$

I will answer the mass moment of inertia tension question. For an infinitesimal clump of mass ${\rm d}m$ located at $\vec{r}$ its effect on the inertia tensor is $${\rm d}{\bf J} = -[\vec{r}\times][\vec{r}\times]{\rm d}m$$

where $[\vec{r}\times]$ is the skew symmetric cross product operator $$ \begin{pmatrix}x\\y\\z \end{pmatrix}\times = \begin{bmatrix}0&-z&y\\z&0&-x\\-y&x&0\end{bmatrix}$$

For a axisymmetric part $\vec{r}(x,y,\theta) = (x,y \cos\theta,y \sin \theta)$ where $x=a\ldots b$, $y=0\ldots f(x)$ and $\theta=0\ldots 2\pi$.

The total mass is $m = \iiint \rho {\rm d}V$ where ${\rm d}V = (y{\rm d}\theta)({\rm d}y)({\rm d}x)$ $$m = \int \limits_a^b \int \limits_0^{f(x)} \int \limits_0^{2\pi} \rho\,y{\rm d}\theta\,{\rm d}y\,{\rm d}x$$ $$\rho = \frac{m}{\pi \int_a^b {f(x)}^2\,{\rm d}x}$$

The inertia tensor is ${\bf J}=\rho \iiint -[\vec{r}\times][\vec{r}\times]\,{\rm d}V$. You can do some math to find that the inner integral is

$$\int \limits_0^{2\pi} -[\vec{r}\times][\vec{r}\times]\,{\rm d}\theta = \begin{bmatrix} 2\pi y^2&0&0\\0&\pi(2 x^2+y^2)&0\\0&0&\pi(2 x^2+y^2)\end{bmatrix} $$

In the end, I get $${\bf J} = \rho\, \begin{bmatrix} \frac{\pi}{2}\int_a^b{f(x)}^4\,{\rm d}x&0&0\\ 0&\frac{\pi}{4}\int_a^b\left({f(x)}^4+4 x^2 {f(x)}^2\right)\,{\rm d}x &0\\0&0&\frac{\pi}{4}\int_a^b\left({f(x)}^4+4 x^2 {f(x)}^2\right)\,{\rm d}x \end{bmatrix}$$

$\endgroup$
  • $\begingroup$ Hi, ja72. Could you help me look at this problem? Thanks! physics.stackexchange.com/questions/165094/… $\endgroup$ – Sherry Feb 14 '15 at 21:16
  • $\begingroup$ $$\int \limits_0^{2\pi} -[\vec{r}\times][\vec{r}\times]\,{\rm d}\theta = \begin{bmatrix} 2\pi y^2&0&0\\0&\pi(2 x^2+y^2)&0\\0&0&\pi(2 x^2+y^2)\end{bmatrix} $$ Why this matrix is diagonal? I tried to compute it and it is not diagonal? $\endgroup$ – Sherry Feb 21 '15 at 19:39
  • $\begingroup$ Did you use $\vec{r}(x,y,\theta) = (x,y \cos\theta,y \sin \theta)$ and $[\vec{r}\times] = \begin{bmatrix}0&-y \sin\theta&y \cos \theta \\ y \sin\theta & 0 & -x \\ -y \cos\theta & x & 0 \end{bmatrix}$ the 3×3 skew symmetric matrix? Then perform the integral over $\theta$ and it will become diagonal since the off diagonal terms are of $\cos \theta$ or $\sin \theta$ form. $\endgroup$ – ja72 Feb 21 '15 at 23:35

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.