Eigenvectors of Inertia tensor

In summary, the three Eigenvectors of the inertia tensor are mutually orthogonal, and for a system with N>=3, the dot product of the i-th and j-th components of the particle's position vectors is zero.
  • #1
Derivator
149
0
Hi,

I've written a little fortran code that computes the three Eigenvectors [itex]\vec{v}_1[/itex], [itex]\vec{v}_2[/itex], [itex]\vec{v}_3[/itex] of the inertia tensor of a N-Particle system.
Now I observed something that I cannot explain analytically:
Assume the position vector [itex]\vec{r}_i[/itex] of each particle to be given with respect to the center of mass of the system.
Then define three new vectors [itex]\vec{\omega}_j := (\vec{v}_j\times\vec{r}_1,\dots,\vec{v}_j\times\vec{r}_N)[/itex] where j=1,...,3. These new vectors are of length 3*N.
Now, for a non-linear configuration of the [itex]\vec{r}_i[/itex] and N>=3, the [itex]\vec{\omega}_j [/itex] seem to be mutually orthogonal, that is [itex]\vec{\omega}_j \cdot \vec{\omega}_i = 0 [/itex] for [itex]i \neq j[/itex] (At least, I obtain this numerically up to machine precision)

I have no analytical explanation for this...

The most promising ansatz I tried so far is:
[itex]\vec{\omega}_i \cdot \vec{\omega}_j = \sum_{l=1}^N (\vec{v}_i\times \vec{r}_l)(\vec{v}_j \times \vec{r}_l) = \sum_{l=1}^N (\vec{v}_i\cdot\vec{v}_j)(\vec{r}_l\cdot\vec{r}_l)-(\vec{v}_i\cdot\vec{r}_l)(\vec{r}_l\cdot\vec{v}_j)= -\sum_{l=1}^N (\vec{v}_i\cdot\vec{r}_l)(\vec{r}_l\cdot\vec{v}_j)[/itex]

Where I have used the relation
[itex](\mathbf{a \times b})\mathbf {\cdot}(\mathbf{c}\times \mathbf{d}) = (\mathbf{a \cdot c})(\mathbf{b \cdot d}) - (\mathbf{a \cdot d})(\mathbf{b \cdot c})[/itex]
(see: http://en.wikipedia.org/wiki/Quadruple_product)

and the fact that the eigenvectors of the inertia tensor are mutually orthogonal:

[itex]\vec{v}_i\cdot\vec{v}_j = 0[/itex]

Unfortunately, this is the point where I'm stuck. I don't see why [itex]-\sum_{l=1}^N (\vec{v}_i\cdot\vec{r}_l)(\vec{r}_l\cdot\vec{v}_j)[/itex] should vanish.derivator
 
Last edited:
Physics news on Phys.org
  • #2
I'm a bit confused on how you're defining your ##\vec{w}_i##. Can you explain? I'm not quite seeing how you're defining it (maybe it's that I don't understand the physical significance).

Also, the dot product won't vanish in your ansatz. The simplest case I can think of is considering a space in ##\mathbb{R}^3## when ##N=1##. Dotting a random vector in the x-y plane with both the x-axis and y-axis. You'll get a vectors along the positive and/or negative z-axis. The dot product of those two vectors isn't zero.
 
Last edited:
  • #3
Hi,
the [itex]\vec{\omega}_i[/itex] are just the hypervectors build from the cross products of the eigenvectors of the inertia tensor and the particle positions (length 3*N). Frankly, I don't know, if they have any physical significance. They just happen to be an intermediate step in my calculation and I noticed their orthogonality. while playing around with my code.

Of corse, you are right, that for N=1 the dot product will not vanish. However, for N >=3 , as long as the particles are not in one line, it does.

Cheers,
derivator
 
  • #4
Taking the eigenvectors to be normalized your last sum is a sum on the product of the i-th and j-th components (with the eigenvectors as the basis) of the position vectors of the particles. Since you have defined the position vectors relative to the center of mass you know the following holds for each component: $$ \sum_l r_l^x=0. $$ I don't know the solution to the problem, but maybe this will help cast it in a new light.
 
  • #5
Hi Haborix,

thanks a lot for your idea. Unfortunately, I also don't see how to use this relation...

Still hoping, someone might see why the sum is vanishing.

Cheers,
derivator
 
  • #6
Maybe, I just had an idea:
Am I right, that the [itex]-\sum_{l=1}^N (\vec{v}_i\cdot\vec{r}_l)(\vec{r}_l\cdot\vec{v}_j)[/itex] are nothing else than the off-diagonal components of the inertia tensor when being expressed in the coordinate system defined by the eigenvectors of the inertia tensor? (for the components of the inertia tensor see: http://farside.ph.utexas.edu/teaching/336k/Newtonhtml/node64.html)
Now, if I did not overlook anything, then its trivial why the
[itex]-\sum_{l=1}^N (\vec{v}_i\cdot\vec{r}_l)(\vec{r}_l\cdot\vec{v}_j)[/itex] vanish: It's simply because the inertia tensor is (trivially) diagonal in its own `eigenframe'.

Let me know what you think.

Cheers,
derivator
 

Related to Eigenvectors of Inertia tensor

1. What are eigenvectors of the inertia tensor?

Eigenvectors of the inertia tensor are a set of orthogonal vectors that represent the principal axes of an object's mass distribution. They describe the direction and orientation of an object's rotational inertia.

2. How are eigenvectors of the inertia tensor calculated?

Eigenvectors of the inertia tensor can be calculated by solving the eigenvalue problem, where the tensor is multiplied by a vector and the resulting vector is equal to a scalar multiple of the original vector. This scalar multiple is known as the eigenvalue.

3. What is the significance of eigenvectors of the inertia tensor?

Eigenvectors of the inertia tensor are important in physics and engineering because they provide information about the rotational dynamics of an object. They can be used to calculate the moments of inertia and the angular momentum of an object.

4. Can eigenvectors of the inertia tensor change?

The eigenvectors of the inertia tensor are fixed for a given object and do not change unless there is a change in the object's mass distribution. However, the magnitude of the eigenvectors can change if the object's mass distribution changes.

5. How are eigenvectors of the inertia tensor used in real-world applications?

Eigenvectors of the inertia tensor are used in various applications, such as in spacecraft control, robotics, and computer graphics. They are also used in the analysis of rigid body motion and in the design of structures and machines to ensure proper stability and balance.

Similar threads

Replies
1
Views
708
Replies
2
Views
1K
  • Introductory Physics Homework Help
Replies
3
Views
367
Replies
3
Views
925
Replies
3
Views
748
Replies
2
Views
816
Replies
4
Views
839
Replies
42
Views
5K
Back
Top