Inertia tensor of triangle

The inertia tensor $\mathbf{𝐉}$ of a triangle (like the inertia tensor of any body) can be expressed in terms of covariance $\mathbf{𝐂}$ of the body:

$\mathbf{𝐉}=\mathrm{t}\mathrm{r}(\mathbf{𝐂})\mathbf{𝐈}-\mathbf{𝐂}$

where covariance is defined as area integral over the triangle:

$\mathbf{𝐂}\triangleq \underset{\mathrm{\Delta }}{\int }\rho \mathbf{𝐱}{\mathbf{𝐱}}^{\mathrm{T}}dA$

Covariance for a triangle in three-dimensional space, assuming that mass is equally distributed over the surface with unit density, is

$\mathbf{𝐂}=a{\mathbf{𝐕}}^{\mathrm{T}}\mathbf{𝐒}\mathbf{𝐕}$

where

• $\mathbf{𝐕}$ represents a 3 × 3 matrix containing triangle vertex coordinates $({\mathbf{𝐯}}_0,{\mathbf{𝐯}}_1,{\mathbf{𝐯}}_2)$ in the rows,
• $a=|({\mathbf{𝐯}}_1-{\mathbf{𝐯}}_0)\times ({\mathbf{𝐯}}_2-{\mathbf{𝐯}}_0)|$ is twice the area of the triangle,
• $\mathbf{𝐒}=\frac{1}{24}\left[\begin{array}{ccc}2& 1& 1\\ 1& 2& 1\\ 1& 1& 2\\ \end{array}\right]$

Substitution of triangle covariance in definition of inertia tensor gives eventually

$\mathbf{𝐉}=\frac{a}{24}({\mathbf{𝐯}}_0^2+{\mathbf{𝐯}}_1^2+{\mathbf{𝐯}}_2^2+({\mathbf{𝐯}}_0+{\mathbf{𝐯}}_1+{\mathbf{𝐯}}_2{)}^2)\mathbf{𝐈}-a{\mathbf{𝐕}}^{\mathrm{T}}\mathbf{𝐒}\mathbf{𝐕}$

A proof of the formula

The proof given here follows the steps from the article.

Covariance of a canonical triangle

Let's compute covariance of the right triangle with the vertices (0,0,0), (1,0,0), (0,1,0).

Following the definition of covariance we receive

${\mathbf{𝐂}}_{xx}^0=\underset{\mathrm{\Delta }}{\int }{x}^{2}dA=\underset{x=0}{\overset{1}{\int }}{x}^2\underset{y=0}{\overset{1-x}{\int }}dydx=\underset{0}{\overset{1}{\int }}{x}^2(1-x)dx=\frac{1}{12}$

${\mathbf{𝐂}}_{xy}^0=\underset{\mathrm{\Delta }}{\int }xydA=\underset{x=0}{\overset{1}{\int }}x\underset{y=0}{\overset{1-x}{\int }}ydydx=\underset{0}{\overset{1}{\int }}x\frac{(1-x{)}^2}{2}dx=\frac{1}{24}$

${\mathbf{𝐂}}_{yy}^0={\mathbf{𝐂}}_{xx}^0$

The rest components of $C$ are zero because the triangle is in $z=0$.

As a result,

${\mathbf{𝐂}}^0=\frac{1}{24}\left[\begin{array}{ccc}2& 1& 0\\ 1& 2& 0\\ 0& 0& 0\\ \end{array}\right]=\frac{1}{48}\left[\begin{array}{c}1\\ -1\\ 0\end{array}\right]{\left[\begin{array}{ccc}1& -1& 0\end{array}\right]}^{\mathrm{T}}+\frac{1}{16}\left[\begin{array}{c}1\\ 1\\ 0\end{array}\right]{\left[\begin{array}{ccc}1& 1& 0\end{array}\right]}^{\mathrm{T}}$

Covariance of the triangle with a vertex in the origin

Consider a linear operator

${\mathbf{𝐱}}^{\prime }=\mathbf{𝐀}{\mathbf{𝐱}}^0$

that maps the canonical triangle in the triangle $\mathbf{𝐯}{\text{'}}_0=\mathbf{𝟎}$, $\mathbf{𝐯}{\text{'}}_1={\mathbf{𝐯}}_1-{\mathbf{𝐯}}_0$, $\mathbf{𝐯}{\text{'}}_2={\mathbf{𝐯}}_2-{\mathbf{𝐯}}_0$. The first two columns of $\mathbf{𝐀}$ contain $\mathbf{𝐯}{\text{'}}_1$ and $\mathbf{𝐯}{\text{'}}_2$ respectively, while the third column is arbitrary. The target triangle is equal to the triangle in question (in particular their areas are equal), but shifted with its zero vertex in the origin.

${\mathbf{𝐂}}^{\prime }=\underset{{\mathrm{\Delta }}^{\prime }}{\int }{\mathbf{𝐱}}^{\prime }\mathbf{𝐱}{\text{'}}^{\mathrm{T}}d{A}^{\prime }=\underset{{\mathrm{\Delta }}^0}{\int }\mathbf{𝐀}{\mathbf{𝐱}}^0{\mathbf{𝐱}}^{0\mathrm{T}}{\mathbf{𝐀}}^{\mathrm{T}}ad{A}^0=a\mathbf{𝐀}{\mathbf{𝐂}}^0{\mathbf{𝐀}}^{\mathrm{T}}$

${\mathbf{𝐂}}^{\prime }=\frac{a}{48}({\mathbf{𝐯}}_1-{\mathbf{𝐯}}_2)({\mathbf{𝐯}}_1-{\mathbf{𝐯}}_2{)}^{\mathrm{T}}+\frac{a}{16}({\mathbf{𝐯}}_1+{\mathbf{𝐯}}_2-2{\mathbf{𝐯}}_0)({\mathbf{𝐯}}_1+{\mathbf{𝐯}}_2-2{\mathbf{𝐯}}_0{)}^{\mathrm{T}}$

Covariance of the triangle in question

The last thing remaining to be done is to conceive how covariance is changed with the translation of all points on vector ${\mathbf{𝐯}}_0$.

$\mathbf{𝐂}=\underset{\mathrm{\Delta }}{\int }({\mathbf{𝐱}}^{\prime }+{\mathbf{𝐯}}_0)({\mathbf{𝐱}}^{\prime }+{\mathbf{𝐯}}_0{)}^{\mathrm{T}}dA={\mathbf{𝐂}}^{\prime }+\frac{a}{2}({\mathbf{𝐯}}_0{\mathbf{𝐯}}_0^{\mathrm{T}}+{\mathbf{𝐯}}_0\overline{\mathbf{𝐱}}{\text{'}}^{\mathrm{T}}+{\overline{\mathbf{𝐱}}}^{\prime }{\mathbf{𝐯}}_0^{\mathrm{T}})$

where

${\overline{\mathbf{𝐱}}}^{\prime }=\underset{{\mathrm{\Delta }}^{\prime }}{\int }{\mathbf{𝐱}}^{\prime }d{A}^{\prime }=\frac{1}{3}(\mathbf{𝐯}{\text{'}}_1+\mathbf{𝐯}{\text{'}}_2)=\frac{1}{3}({\mathbf{𝐯}}_1+{\mathbf{𝐯}}_2-2{\mathbf{𝐯}}_0)$

is the centroid of the dashed triangle.

It's easy to check now that all coefficients in $\mathbf{𝐂}$ before ${\mathbf{𝐯}}_i{\mathbf{𝐯}}_i^{\mathrm{T}}$ is $\frac{a}{12}$ and before ${\mathbf{𝐯}}_i{\mathbf{𝐯}}_j^{\mathrm{T}}(i\ne j)$ is $\frac{a}{24}$. This can be expressed in matrix form with $\mathbf{𝐒}$ as above.

References

This article "Inertia tensor of triangle" is from Wikipedia. The list of its authors can be seen in its historical and/or the page Edithistory:Inertia tensor of triangle. Articles copied from Draft Namespace on Wikipedia could be seen on the Draft Namespace of Wikipedia and not main one.

📰 Article(s) of the same category(ies)

• Parallel force
• Jacobi point
• Amber Nova
• Bibhorr angle
• Trigonometry
• Mechanics of planar particle motion
• Body moment

---