Message Boards Message Boards

0
|
1369 Views
|
1 Reply
|
0 Total Likes
View groups...
Share
Share this post:

How to find the moment of inertia of a symbolic octahedron at an angle?

I want to find the moment of inertia of the following symbolic octahedron:

Octahedron[{\[Theta], \[CurlyPhi]}, \[ScriptL]]

I can find the moment of inertia of a symbolic edge length:

MomentOfInertia[Octahedron[\[ScriptL]]] // MatrixForm

Output:

{{\[ScriptL]^5/(15 Sqrt[2]), 0, 0}, {0, \[ScriptL]^5/(15 Sqrt[2]), 
  0}, {0, 0, \[ScriptL]^5/(15 Sqrt[2])}}

I can't get MomentOfInertia to evaluate with a symbolic edge and an angle:

MomentOfInertia[Octahedron[{1, 2}, l], 
 Assumptions -> {\[ScriptL] > 0}]

I would like to be able to evaluate the following:

MomentOfInertia[Octahedron[{\[Theta], \[CurlyPhi]}, \[ScriptL]]]

Is there a way to find a rotated octahedron's symbolic moment of inertia? I realize this might not be possible. As an example of a concrete case, if the first angle is 1 and the second angle is 2 and the edge length is 3 we have:

MomentOfInertia[Octahedron[{1, 2}, 3]]

The output is a tensor:

{{81/20 (Cos[1]^2 + Cos[2]^2 + Sin[1]^2 + Cos[1]^2 Sin[2]^2 + 
     Sin[1]^2 Sin[2]^2) (Sqrt[2] Cos[1]^2 Cos[2]^2 + 
     Sqrt[2] Cos[2]^2 Sin[1]^2 + Sqrt[2] Cos[1]^2 Sin[2]^2 + 
     Sqrt[2] Sin[1]^2 Sin[2]^2), 0, 
  81/20 (-Cos[2] Sin[2] + Cos[1]^2 Cos[2] Sin[2] + 
     Cos[2] Sin[1]^2 Sin[2]) (Sqrt[2] Cos[1]^2 Cos[2]^2 + 
     Sqrt[2] Cos[2]^2 Sin[1]^2 + Sqrt[2] Cos[1]^2 Sin[2]^2 + 
     Sqrt[2] Sin[1]^2 Sin[2]^2)}, {0, 
  81/20 (1 + Cos[1]^2 + Sin[1]^2) (Sqrt[2] Cos[1]^2 + 
     Sqrt[2] Sin[1]^2) (Cos[2]^2 + Sin[2]^2)^2, 
  0}, {81/20 (-Cos[2] Sin[2] + Cos[1]^2 Cos[2] Sin[2] + 
     Cos[2] Sin[1]^2 Sin[2]) (Sqrt[2] Cos[1]^2 Cos[2]^2 + 
     Sqrt[2] Cos[2]^2 Sin[1]^2 + Sqrt[2] Cos[1]^2 Sin[2]^2 + 
     Sqrt[2] Sin[1]^2 Sin[2]^2), 0, 
  81/20 (Cos[1]^2 + Cos[1]^2 Cos[2]^2 + Sin[1]^2 + Cos[2]^2 Sin[1]^2 +
      Sin[2]^2) (Sqrt[2] Cos[1]^2 Cos[2]^2 + 
     Sqrt[2] Cos[2]^2 Sin[1]^2 + Sqrt[2] Cos[1]^2 Sin[2]^2 + 
     Sqrt[2] Sin[1]^2 Sin[2]^2)}}
POSTED BY: Peter Burbery

Hi Peter,

the result of the expression above is % //TrigToExp//TableForm

{{81 / (5 Sqrt[2]), 0, 0}, {0, 81/(5 Sqrt[2]), 0}, {0, 0, 81/( 5 Sqrt[2])}}

As you may notice this is a multiple of the unit matrix, so as all regular polytopes it has the inertia tensor of a sphere, invariant to rotations.

In general with respect to a set of three orthogonal axis the moment of inertial is a symmetric positive 3x3-matrix, that can be rotated by an orthogonal matirix O ( n, phi) for a rotation of the body around a unit vector axis n and an angle phi

I' = O(n,phi) . I . Transpose[O(n,phi)]

A translation of the fixed point with respect to the body-fixed system is given by Steiners theorem, that is in essence adding the inertia tensor of the point mass fixed at end of the translation vector.

Regards Roland

POSTED BY: Roland Franzius
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard

Group Abstract Group Abstract