Physics Library
 An open source physics library
Encyclopedia | Forums | Docs | Random  

The Inertia Tensor and Products of Inertia

(Topic)

The Inertia Tensor and Products of Inertia

For rotation about one fixed axis, a scalar moment of inertia is sufficient:

    ∫
I =    r2⊥ dm.
(1)

In general three-dimensional rotation, however, one scalar is not enough.

The angular velocity

ω
(2)

can point in any direction, and the angular momentum

L
(3)

need not be parallel to it.

The mass distribution must therefore be described by a second-order tensor:

|--------|
L--=-Iω.--
(4)

The matrix representation of the inertia tensor in Cartesian coordinates is

|---⌊-------------⌋--|
|     Ixx  Ixy Ixz   |
I = ⌈ Iyx  Iyy Iyz⌉ .|
|     I    I    I    |
-------zx---zy---zz----
(5)

For a mass distribution,

|-----∫--------------------|
|        ( 2         )     |
|Iij =     r δij − xixj  dm. |
----------------------------
(6)

This article derives that expression, explains the products of inertia, and connects the tensor back to the familiar scalar moment of inertia.

PIC

Figure 1. A scalar moment of inertia describes rotation about one specified axis. The inertia tensor describes the rotational mass distribution for arbitrary angular-velocity directions.

1 Why a tensor is needed

For a point mass,

L  = r × mv.
(7)

For rigid body rotation about an origin,

v =  ω × r.
(8)

Therefore

L =  r × m (ω × r) .
(9)

Use the vector triple-product identity

a × (b × c) = b(a ⋅ c) − c(a ⋅ b).
(10)

Then

                2
r × (ω ×  r) = r ω − r(r ⋅ ω).
(11)

Thus a particle contributes

-------------------------
|       [ 2           ]  |
-L-=-m---r-ω-−-r(r-⋅ ω-)-.|
(12)

This dependence on the direction of ω is more complicated than multiplication by one scalar.

PIC

Figure 2. For a particle in rigid body rotation, the velocity is ω × r. The resulting angular momentum contains both the magnitude and direction of the particle position relative to the rotation origin.

2 Derivation for a system of particles

For particles indexed by a,

    ∑      [ 2              ]
L =     ma  raω −  ra(ra ⋅ ω ) .
     a
(13)

Write the Cartesian components:

     ∑      [ 2            ]
Li =     ma  raδij − xa,ixa,j ωj.
      a
(14)

Using the summation convention,

|----------|
Li-=--Iijωj,-
(15)

where

|-----∑------[--------------]--|
|Iij =    ma  r2aδij − xa,ixa,j  .|
|       a                      |
-------------------------------
(16)

For a continuous body,

|-----∫--[-----------]-----|
|Iij =     r2δij − xixj  dm. |
---------------------------|
(17)

3 Matrix form

In Cartesian coordinates,

|--------------------------------------------------------|
|    ⌊∫    2    2          ∫               ∫          ⌋  |
|    |   (y + z ) dm     −   xy dm       −    xz dm   |  |
|    |    ∫            ∫                   ∫          |  |
|I = ||  −    xy dm       (x2 + z2)dm     −    yz dm   || .|
|    |⌈    ∫                ∫           ∫              |⌉  |
|       −    xz dm       −    yzdm       (x2 + y2) dm    |
---------------------------------------------------------|
(18)

The diagonal entries are familiar scalar moments:

Ixx = ∫ (y2 + z2) dm, (19)
Iyy = ∫ (x2 + z2) dm, (20)
Izz = ∫ (x2 + y2) dm. (21)

PIC

Figure 3. The diagonal inertia-tensor entries are ordinary moments of inertia about the coordinate axes. The off-diagonal entries encode products of coordinates and couple angular-velocity components to angular-momentum components.

4 Products of inertia and sign convention

The off-diagonal integrals involve coordinate products.

Define the geometric products of inertia

Pxy = ∫ xy dm, (22)
Pxz = ∫ xz dm, (23)
Pyz = ∫ yz dm. (24)

With this convention, the inertia tensor is

|----------------------------|
|    ⌊                    ⌋  |
|       Ixx   − Pxy  − Pxz   |
|I = ⌈ − Pxy   Iyy   − Pyz ⌉ .
|      − Pxz  − Pyz   Izz    |
-----------------------------
(25)

Physics texts also commonly denote the actual tensor component itself by

        ∫
Ixy = −    xydm  =  − Pxy,
(26)

with analogous definitions for Ixz and Iyz.

Therefore:

|--------------------------------------------|
-Ixy-=-−-Pxy,-----Ixz =-−-Pxz,----Iyz-=-−-Pyz.-
(27)

The sign convention must always be checked when comparing references.

In this article:

  • Pxy means the unsigned coordinate-product integral ∫ xy dm,
  • Ixy means the off-diagonal tensor component −∫ xy dm.

5 The inertia tensor is symmetric

Because

xixj = xjxi,
(28)

|--------|
Iij =-Iji.-
(29)

Thus

|-------|
IT-=-I.--
(30)

Only six independent components exist:

Ixx,  Iyy,  Izz,  Ixy,  Ixz,  Iyz.
(31)

This symmetry will be crucial when principal axes are introduced in M05-10.

6 Angular momentum components

From

L = Iω,
(32)

Lx = Ixxωx + Ixyωy + Ixzωz, (33)
Ly = Ixyωx + Iyyωy + Iyzωz, (34)
Lz = Ixzωx + Iyzωy + Izzωz. (35)

Therefore an angular velocity along one coordinate direction can generate angular momentum components along other directions if the products of inertia are nonzero.

7 When angular momentum is not parallel to angular velocity

Suppose the inertia tensor is diagonal:

    ⌊I    0  0 ⌋
    ⌈ 1        ⌉
I =   0  I2  0   .
      0   0  I3
(36)

Then

L = I1ω1e1 + I2ω2e2 + I3ω3e3.
(37)

If

I1 ⁄= I2
(38)

and both

ω1 ⁄= 0,     ω2 ⁄= 0,
(39)

the direction of L generally differs from the direction of ω.

Thus even zero products of inertia do not guarantee

L ∥ ω
(40)

for every angular velocity.

Parallelism occurs when ω lies along a principal axis or within a degenerate equal-inertia subspace.

PIC

Figure 6. In three-dimensional rotation, angular momentum need not point along angular velocity. The inertia tensor stretches different angular-velocity components by different amounts.

8 Symmetry and vanishing products of inertia

Products of inertia are especially sensitive to symmetry.

Consider

      ∫
P   =    xy dm.
 xy
(41)

If the mass distribution is symmetric under

x →  − x,
(42)

then every mass element at (x,y,z) is paired with an equal mass element at (−x,y,z).

Their contributions to xy cancel.

Therefore

|--------|
-Pxy-=-0-|
(43)

for symmetry about the yz plane.

Similarly:

  • symmetry under y →−y also makes Pxy = 0,
  • symmetry under x →−x or z →−z makes Pxz = 0,
  • symmetry under y →−y or z →−z makes Pyz = 0.

PIC

Figure 4. Products of inertia vanish when symmetry pairs positive and negative coordinate-product contributions. For example, reflection symmetry across the yz plane cancels the integral of xy.

9 Example 1: four point masses in a rectangle

Four equal point masses m lie at

(±a, ±b, 0),
(44)

with all sign combinations present.

The diagonal moments are

Ixx = ∑ my2 = 4mb2, (45)
Iyy = ∑ mx2 = 4ma2, (46)
Izz = ∑ m(x2 + y2) = 4m(a2 + b2). (47)

The product

       ∑
Pxy =     mxy
(48)

contains equal positive and negative terms:

Pxy = 0.
(49)

Also,

Pxz = Pyz = 0
(50)

because

z = 0.
(51)

Thus

|----------------------------------|
|    ⌊4mb2     0         0      ⌋  |
|    ⌈           2              ⌉  |
|I =    0    4ma         02    2   .|
--------0------0----4m-(a--+-b-)---|
(52)

10 Example 2: two asymmetric point masses

Let two point masses have:

m1 = 2 kg, r1 = (1, 2, 0) m, (53)
m2 = 1 kg, r2 = (−1, 1, 0) m. (54)

Then

Ixx = 2(22) + 1(12) = 9 kg m2, (55)
Iyy = 2(12) + 1((−1)2) = 3 kg m2, (56)
Izz = Ixx + Iyy = 12 kg m2. (57)

The geometric product is

Pxy = 2(1)(2) + 1(−1)(1) (58)
= 3 kg m2. (59)

Therefore

Ixy = − 3kg m2.
(60)

The tensor is

|--------------------------|
|    ⌊ 9   − 3   0⌋        |
|    ⌈            ⌉      2 |
|I =  − 3   3    0  kg m  .|
-------0----0---12---------|
(61)

11 Planar lamina

For a thin lamina in the xy plane,

z = 0.
(62)

Therefore

Pxz =  Pyz = 0.
(63)

The tensor becomes

|-------------------------|
|   ⌊  I     − P     0⌋   |
|   ⌈   xx      xy    ⌉   |
I =   − Pxy   Iyy    0  . |
--------0------0----Izz----
(64)

Also,

      ∫
           2   2
Izz =   (x  + y ) dm =  Ixx + Iyy.
(65)

Thus the perpendicular-axis theorem appears naturally as a component identity of the inertia tensor:

|--------------|
|Izz = Ixx + Iyy|
----------------
(66)

for a planar lamina.

12 Uniform rectangular plate

A rectangular plate of mass M has side lengths a along x and b along y, centered at the origin.

By reflection symmetry,

Pxy = 0.
(67)

The moments are

Ixx =  1
---
12Mb2, (68)
Iyy = -1-
12Ma2, (69)
Izz = -1-
12M(a2 + b2). (70)

Therefore

|---------⌊---------------⌋--|
|      M    b2  0     0      |
ICM  = ---⌈ 0   a2    0   ⌉ .|
|      12   0   0   a2 + b2  |
------------------------------
(71)

13 Uniform rectangular box

For a uniform rectangular box centered at the origin with side lengths:

  • a along x,
  • b along y,
  • c along z,

symmetry makes all products of inertia vanish.

The inertia tensor is

|----------⌊-2----2-----------------⌋--|
|       M   b  + c     0        0      |
|ICM =  ---⌈   0     a2 + c2    0   ⌉ .|
|       12     0       0     a2 + b2   |
---------------------------------------|
(72)

These coordinate directions are principal axes.

The full meaning of that statement is developed in M05-10.

14 Rotated slender rod and nonzero product of inertia

Consider a slender uniform rod of length L and mass M, centered at the origin in the xy plane.

Let the rod make angle ϕ with the x axis.

Define

      1     2
I⊥ = ---M L  .
     12
(73)

A point on the rod can be written

r = s(cos ϕex + sin ϕey) .
(74)

The tensor components become

Ixx = I⊥ sin 2ϕ, (75)
Iyy = I⊥ cos 2ϕ, (76)
Ixy = −I⊥ sin ϕ cos ϕ, (77)
Izz = I⊥. (78)

Thus

|-------⌊-----------------------------⌋--|
|           sin2ϕ      − sinϕ cosϕ  0    |
|       ⌈− sin ϕcos ϕ     cos2ϕ     0 ⌉  |
|I = I⊥                                 .|
---------------0------------0-------1----|
(79)

Unless

               π
ϕ = 0,     ϕ = --,
               2
(80)

or another symmetry-aligned orientation is chosen, the product of inertia is nonzero.

PIC

Figure 5. A slender rod aligned with a rotated direction generally has a nonzero xy product of inertia when described in the original x,y coordinate axes.

15 Example 3: rod at forty-five degrees

Let

ϕ = 45∘.
(81)

Then

   2       2     1-
sin ϕ = cos  ϕ = 2 ,
(82)

and

sinϕ cosϕ =  1.
             2
(83)

Thus

|--------------------------|
|      ⌊                ⌋  |
|         1∕2   − 1∕2  0   |
|I = I⊥⌈ − 1∕2   1∕2   0⌉ .|
|          0      0    1   |
----------------------------
(84)

The off-diagonal terms are not a defect.

They simply indicate that the chosen coordinate axes are not aligned with the rod’s principal directions.

16 Moment of inertia about an arbitrary axis

Let

ˆn
(85)

be a unit vector defining an axis through the tensor origin.

If

ω  = ω ˆn,
(86)

then the rotational kinetic energy is

K  = 1-ω ⋅ Iω.
     2
(87)

Therefore

     1- 2  T
K  = 2 ω ˆn  Iˆn.
(88)

Define the scalar moment about the axis:

-------------
|       T    |
-Iˆn-=-ˆn--Iˆn.-|
(89)

This is the tensor form of

    ∫
I =    r2⊥ dm.
(90)

PIC

Figure 7. The scalar moment of inertia about any unit axis n is obtained from the quadratic form In = nT In.

17 Direct proof of the arbitrary-axis formula

The squared perpendicular distance from r to the axis

ˆn
(91)

is

 2    2         2
r⊥ = r  − (r ⋅ ˆn) .
(92)

Therefore

In = ∫ [r2 − (r ⋅ ˆn)2] dm. (93)

In component form,

(r ⋅ ˆn)2 = nixinjxj.
(94)

Thus

In = ni[∫                  ]
    (r2δij − xixj)dmnj (95)
= niIijnj. (96)

Hence

|------------|
|Iˆn = ˆn TIˆn. |
-------------
(97)

18 Example 4: arbitrary axis for a diagonal tensor

Suppose

    ⌊        ⌋
      2  0  0
I = ⌈ 0  5  0⌉ kg m2.
      0  0  8
(98)

Let the unit axis lie in the xy plane:

ˆn = cos 𝜃ex + sin𝜃 ey.
(99)

Then

In = 2 cos 2𝜃 + 5 sin 2𝜃. (100)

Therefore

|---------2---------2--------2-|
-Iˆn =-2cos--𝜃 +-5sin-𝜃--kg-m--.-
(101)

At

𝜃 = 0,
(102)

the result is

      2
2kg m  .
(103)

At

      ∘
𝜃 = 90 ,
(104)

it is

5kg m2.
(105)

19 Rotational kinetic energy

For a system of particles,

       ∑
K  = 1-    ma |ω  × ra|2.
     2  a
(106)

The tensor form is

|--------------|
|K  = 1-ωT Iω. |
------2--------|
(107)

Equivalently,

|------------|
|     1      |
|K =  -ω ⋅ L.|
------2-------
(108)

This relation remains valid even when

L  ⁄∥ ω.
(109)

20 Positive character of the inertia tensor

For any angular velocity,

2K  =  ωT Iω.
(110)

Since kinetic energy cannot be negative,

|----------|
ωT-Iω--≥-0.-
(111)

Thus the inertia tensor is positive semidefinite.

For a genuinely three-dimensional body with nonzero mass away from every possible axis through the origin, it is positive definite.

For idealized lower-dimensional objects, one eigenvalue can be zero.

For example, an ideal slender rod has zero moment of inertia about its own longitudinal axis if its thickness is neglected.

21 Trace of the inertia tensor

Add the diagonal components:

tr(I) = Ixx + Iyy + Izz (112)
= ∫ [  2   2      2    2      2   2 ]
 (y  + z ) + (x +  z ) + (x + y )dm. (113)

Therefore

|---------∫--------|
|            2     |
|tr(I) = 2   r  dm. |
--------------------
(114)

The trace is unchanged by a rotation of coordinate axes.

It equals the sum of the principal moments and provides a useful consistency check.

22 Changing coordinate axes

Suppose coordinate components transform under an orthogonal matrix Q:

r′ = Qr.
(115)

Then

ω ′ = Q ω,
(116)

and

L ′ = QL.
(117)

Because

L = Iω,
(118)

the tensor components in the rotated coordinates satisfy

|------------|
|I′ = QIQT  .|
-------------
(119)

The physical mass distribution has not changed.

Only its matrix representation has changed.

23 Products of inertia depend on the coordinate axes

A body can have

Pxy = 0
(120)

in one coordinate system and a nonzero product in another rotated system.

Therefore a nonzero product of inertia does not mean anything is physically wrong or unbalanced.

It indicates that the coordinate axes are not aligned with directions that diagonalize the inertia tensor.

M05-10 finds those directions systematically.

24 Tensor parallel-axis theorem

M05-05 introduced the scalar parallel-axis theorem.

The three-dimensional tensor generalization is

|--------------[-----------]-|
|IO = ICM +  M  d21 − ddT   ,|
------------------------------
(121)

where:

  • ICM is the inertia tensor about the center of mass,
  • O is the new origin,
  • d is the vector from the new origin to the center of mass,
  • 1 is the identity matrix.

In components,

|---------------------(-----------)--|
|(IO)ij = (ICM )ij + M  d2δij − didj .|
-------------------------------------
(122)

PIC

Figure 8. The tensor parallel-axis theorem shifts the inertia tensor from the center of mass to a new origin by adding M[d21 − ddT ].

25 Recovering the scalar parallel-axis theorem

Let the desired axis have unit direction

ˆn.
(123)

Use

Iˆn,O =  ˆnTIO ˆn.
(124)

Then

In,O = In,CM + M[  2         2]
 d  − (ˆn ⋅ d ). (125)

The bracket is the squared perpendicular separation of the two parallel axes:

 2    2         2
d⊥ = d  − (ˆn ⋅ d) .
(126)

Therefore

|------------------2-|
Inˆ,O-=-Iˆn,CM-+-M--d⊥.-
(127)

This exactly reproduces the scalar parallel-axis theorem.

26 Example 5: tensor shift along the x direction

Suppose

       ⌊          ⌋
        Ix  0   0
I   =  ⌈ 0  I   0 ⌉
 CM          y
         0  0   Iz
(128)

and the new origin is displaced from the center of mass by distance d along the x direction.

Then

d =  dex.
(129)

The shift matrix is

                      ⌊ 0  0  0⌋
  [ 2       T ]      2⌈        ⌉
M  d 1 − dd    =  M d   0  1  0  .
                        0  0  1
(130)

Thus

|-----⌊-------------------------⌋--|
|       Ix      0          0       |
|I  = ⌈ 0   I + M  d2      0    ⌉ .|
| O          y                 2   |
--------0-------0------Iz +-M-d----|
(131)

The moment about the x axis is unchanged because shifting the origin along the axis does not change perpendicular distances to that axis.

27 Example 6: tensor shift in two coordinate directions

Let

d  = aex + bey.
(132)

Then

d2 = a2 + b2.
(133)

The shift term is

   ⌊                   ⌋
     b2   − ab     0
M  ⌈− ab   a2      0   ⌉ .
      0     0   a2 + b2
(134)

Notice that shifting the reference origin can create nonzero off-diagonal tensor components even when the center of mass tensor is diagonal.

28 Products of inertia and balance

A nonzero product of inertia should not be confused automatically with static or dynamic imbalance.

Balance depends on both:

  • the center of mass location,
  • the orientation of the inertia tensor relative to the intended rotation axis.

For a rotor intended to spin about a fixed shaft, desirable design often places:

  • the center of mass on the shaft,
  • a principal axis along the shaft.

The second requirement concerns products of inertia and principal directions.

29 Connection to principal axes

Because the inertia tensor is real and symmetric, it has orthogonal eigenvectors.

Those eigenvectors are the principal axes.

In coordinates aligned with them,

|------------------|
|    ⌊I    0  0 ⌋  |
|    ⌈ 1        ⌉  |
|I =   0  I2  0   .|
-------0---0--I3---|
(135)

The eigenvalues

I ,     I ,    I
 1       2      3
(136)

are the principal moments of inertia.

For rotation exactly about a principal axis,

|------|
L  ∥ ω.|
--------
(137)

M05-10 develops this eigenvalue problem in detail.

30 Common mistakes

  1. Treating moment of inertia as one scalar in arbitrary three-dimensional rotation.
  2. Assuming L is always parallel to ω.
  3. Forgetting the minus sign between Pxy = ∫ xy dm and the tensor component Ixy = −Pxy.
  4. Comparing formulas from different references without checking their product-of-inertia sign convention.
  5. Assuming all products of inertia are automatically zero.
  6. Failing to use symmetry to eliminate odd coordinate-product integrals.
  7. Assuming a diagonal tensor guarantees L ∥ ω for every ω.
  8. Confusing a diagonal coordinate representation with a spherical inertia tensor having equal diagonal entries.
  9. Using In without normalizing the axis vector.
  10. Forgetting that inertia-tensor components depend on both the origin and orientation of the coordinate axes.
  11. Applying the scalar parallel-axis theorem independently to matrix entries without using the tensor shift formula.
  12. Forgetting that an origin shift can create off-diagonal terms.
  13. Treating nonzero products of inertia as proof of center of mass offset.
  14. Confusing the inertia tensor with the stress tensor or area moment tensor.
  15. Forgetting that the inertia tensor is symmetric.

31 Practice exercises

  1. Starting from L = r × m(ω × r), derive the particle inertia tensor.
  2. Extend the particle result to a discrete system and then to a continuous body.
  3. Write the inertia tensor explicitly in Cartesian coordinates and identify all three products of inertia.
  4. Explain the distinction between Pxy = ∫ xy dm and the tensor component Ixy used in this article.
  5. Prove that the inertia tensor is symmetric.
  6. Show that symmetry under x →−x forces Pxy = Pxz = 0.
  7. Compute the inertia tensor of four equal masses at (±a,±b, 0).
  8. Compute the tensor of the two-mass asymmetric example in this article and verify Izz = Ixx + Iyy.
  9. Derive the inertia tensor of a slender rod lying at angle ϕ in the xy plane.
  10. Evaluate the rotated-rod tensor at ϕ = 45∘.
  11. Prove that the scalar moment about an arbitrary unit axis is In = nT In.
  12. Derive K = 1
2ωT Iω.
  13. Show that tr(I) = 2 ∫ r2 dm.
  14. Derive the tensor parallel-axis theorem and recover the scalar theorem from it.
  15. Explain why principal axes correspond to eigenvectors of the inertia tensor.

32 Summary

The angular momentum of a rigid body rotating about a chosen origin is

|--------|
-L-=-Iω.--
(138)

The inertia tensor is

|-----∫--------------------|
|        ( 2         )     |
|Iij =     r δij − xixj  dm. |
----------------------------
(139)

Its Cartesian form is

|----⌊--------------------⌋--|
|       Ixx   − Pxy  − Pxz   |
|I = ⌈ − Pxy   Iyy   − Pyz ⌉ .
|      − Pxz  − Pyz   Izz    |
-----------------------------|
(140)

The tensor is symmetric:

|T------|
I--=-I.--
(141)

The scalar moment about an arbitrary unit axis is

|-----------|
Iˆn = ˆnT Iˆn. |
-------------
(142)

Rotational kinetic energy is

|--------------|
|     1   T    |
|K  = 2-ω  Iω. |
---------------
(143)

Under a coordinate rotation,

|-′-------T--|
-I-=--QIQ---.|
(144)

Under an origin shift from the center of mass,

|--------------[-2-------T-]-|
-IO-=-ICM-+--M--d-1-−-dd----.-
(145)

M05-10 uses the symmetry of I to identify principal axes and principal moments of inertia.

References

References

[1]   J. R. Taylor, Classical Mechanics, University Science Books, 2005.

[2]   H. Goldstein, C. Poole, and J. Safko, Classical Mechanics, 3rd ed., Addison-Wesley, 2002.

[3]   L. D. Landau and E. M. Lifshitz, Mechanics, 3rd ed., Butterworth-Heinemann, 1976.

[4]   K. R. Symon, Mechanics, 3rd ed., Addison-Wesley, 1971.


"The Inertia Tensor and Products of Inertia" is owned by bloftin.
(view preamble)
View style:
Other names:  M05-09
Also defines:  inertia tensor, product of inertia, tensor component of inertia, inertia dyadic, tensor parallel-axis theorem
Keywords:  inertia tensor, products of inertia, rigid body, angular momentum, angular velocity, rotational kinetic energy, moment of inertia, tensor parallel-axis theorem, symmetry, coordinate transformation, principal axes

Attachments:
GRE Physics Companion: The Inertia Tensor and Products of Inertia (Example) by bloftin

Cross-references: theorem, forces, formulas, spin, static, center of mass, parallel-axis theorem, M05-05, coordinate system, trace, kinetic energy, relation, system, unit, rotational kinetic energy, unit vector, perpendicular-axis theorem, principal axis, M05-10, summation convention, position, magnitude, velocity, particle, identity, vector, rigid body, Cartesian coordinates, representation, matrix, tensor, mass, angular momentum, angular velocity, moment of inertia, scalar
There are 9 references to this object.

This is version 1 of The Inertia Tensor and Products of Inertia, born on 2026-10-17.
Object id is 1480, canonical name is InertiaTensorAndProductsOfInertia.
Accessed 9 times total.

Classification:
Physics Classification: 45.40.-f (Dynamics and kinematics of rigid bodies)
 45.20.Dd (Newtonian mechanics)

Pending Errata and Addenda

None.

Discussion

No messages.

Interact