How to Calculate Moment of Inertia for Polygons in a 2D System?

In summary, to calculate the moment of inertia of a polygon in a 2D system, one can use a double integral without the need for triangulation. It is possible to use the same method for complex polygons with crossing edges. The moment of inertia is calculated relative to the polygon's centroid, which is already known. The polygon should be "star-shaped" with respect to the centroid for the formula to be valid. The formula involves the position vectors of the vertices, the density of the polygon, and a summation. The norm of a vector is used in the formula. When using different units for mass and density, the calculated moment of inertia may be off by a factor.
  • #1
Aidman
16
0
How can one calculate the moment of inertia of a polygon? Assuming that one knows the polygon’s vertexes which in turn are connected by straight lines in a 2D system? If the calculation is possible without triangulating the polygon, is it then also possible to use the same method with complex polygons, where some edges cross? Any helpful information is highly appreciated.

EDIT: vertexes = vertices :redface:
 
Last edited:
Mathematics news on Phys.org
  • #2
it's 'just' a double integral, no triangulation nreeded. it might not be a nice integral though
 
  • #3
But the polygon has a continuous area/mass shouldn’t it then be possible to use a method which resolves the inertia by calculating the polygon's total mass and centroid?
 
Last edited:
  • #4
1. What do you mean by "continuous mass"?
a) Continuous mass distribution (i.e the density function is continuous)?
b) Constant density?

2. What axis do you wish to calculate the moment of inertia about, or do you want to find the matrix of inertia?
 
  • #5
By continuous mass I mean that the polygon consists of mass elements/points which are connected to each other with constant density.

And the moment of inertia should be relative to the polygon's centroid, which is already known.
 
  • #6
Are your polygons 3-D structures which are symmetric around each of the 3 independent axes they may rotate about?
 
  • #7
Solid polygons, frames, or vertex points?

It seems like another easy option is to work out the moment of inertia for a general triangle, and then disect the polygon.

Note that you fail to specify the axis of rotation, which is significant.
 
  • #8
Aidman said:
And the moment of inertia should be relative to the polygon's centroid, which is already known.
Just to clarify:
"moment of inerta" is calculated relative to a specified axis (going through a point), not relative to a point, as NateTG comments.

Unless you specifically assume that the (instantaneous) rotation axis has a constant direction (through time, that is), it is insufficient to calculate only the moments of inertae (around the 3 axes); you must calculate the full, inertial matrix, which include the products of inertiae, and not only the moments.
 
  • #9
As I mentioned before the polygons are in a 2D system, which limits their rotation to one axle, therefore calculating a single inertia is sufficient.

And yes dealing with triangles are a lot simpler but the process of triangulating a polygon seems very difficult, especially if it also should handle concave polygons.
 
  • #10
actually triangulation is easy, and can be done using many methods, i believe erdos had a proof of some bounds
 
  • #11
matt grime said:
actually triangulation is easy, and can be done using many methods, i believe erdos had a proof of some bounds

Considering that inverse triangles are acceptable for this, it's extremely easy.
 
  • #12
I am having a hard time believing it’s very easy considering I have found very little information about triangulation.
 
Last edited:
  • #13
Aidman:
If I have understood you correctly, you wish to express the moment of inertia of your polygon in terms of the position vectors of the vertices (given that we place the centroid in the origin).

Let [tex]\rho[/tex] be the (constant) density of the polygon, let N be the number of vertices, and let [tex]\vec{P}_{n}[/tex] be the position vector of the n'th vertex.
(We also define [tex]\vec{P}_{N+1}\equiv\vec{P}_{1}[/tex])

Then we have, that the moment of inertia I about the centroid can be expressed as:
[tex]I=\sum_{n=1}^{N}\frac{\rho}{12}||\vec{P}_{n+1}\times\vec{P}_{n}||(\vec{P}_{n+1}^{2}+\vec{P}_{n+1}\cdot\vec{P}_{n}+\vec{P}_{n}^{2})[/tex]

Note:
I have assumed that the polygon is "star-shaped" with respect to the centroid in the derivation of the formula; i.e., in particular, that the whole triangle lying between the centroid and two adjacent vertices is included in the polygon.
 
Last edited:
  • Like
Likes trytodoit
  • #14
I am sorry, but what exactly does the double lines [tex]||...||[/tex] mean?
 
Last edited:
  • #15
That means the norm of the vector inside the double lines
 
  • #16
So I can ignore the double lines for 2D vectors, as the vectorproduct returns a value and not a vector?
 
  • #17
The cross product of two vectors always return a vector;
if we have V1=(x1,y1,0), V2=(x2,y2,0), then the crossproduct V1(cross)V2=(0,0,x1y2-x2y1).
Just remember that the norm is non-negative, so the norm above is simply abs(x1y2-x2y1), where abs designates the absolute value.
 
  • #18
Ok, many thanks. I have tested the algorithm but for some reason it calculates a moment of inertia that is 10000 times greater then the excepted value. For example when I use it for a square, with the sides 100, it resolves the value 16666666.66... but if I am not mistaken the correct value should be 1666.66...
 
Last edited:
  • #19
Have you used unit mass or unit density in that calculation?
If you have used the algorithm with unit density, while you calculated the known moment of inertia with unit mass, then the expression from the algorithm should be 10000 times bigger for a 100 times 100 square.
 
  • #20
Sorry, I am not following... If I am using [tex][m^2][/tex] as unit area and [tex][kg][/tex] as unit mass how can I regulate the calculated inertia [tex][kgm^2][/tex], for any given density [tex][kg/m^2][/tex]?
 
Last edited:
  • #21
Allright, I'll make what I said a bit clearer:
1) Usual fomulae for moment of inertia:
Usually the moment of inertia is given in the form:
[tex]I=kmr^{2}[/tex], where m is the mass of the object, r is a typical length scale in the object, while k is some number.
(If k is 1, then r is the gyration radius of the object)
Examples:
Sphere:[tex]\frac{2}{5}mR^{2}[/tex], where R is the radius of the sphere.
Rod about center of mass/centroid::[tex]\frac{1}{12}mL^{2}[/tex], where L is the length of the rod.

2) I have expressed my moment of inertia with a density, rather than with mass, that is, I haven't the mass as an explicit parameter in the expression.

3) So, my question was if you had used a mass of 1kg in the calculation from the known formula, while you had used a density of 1kg/m^{2} in the algorithm.
(That would have explained the factor of 10000)

4. If you want to have a formula for the moment of inertia where the mass is explicit, here you have it:

[tex]I=\frac{m}{6}\frac{\sum_{n=1}^{N}||\vec{P}_{n+1}\times\vec{P}_{n}||(\vec{P}^{2}_{n+1}+\vec{P}_{n+1}\cdot\vec{P}_{n}+\vec{P}_{n}^{2})}{\sum_{n=1}^{N}||\vec{P}_{n+1}\times\vec{P}_{n}||}[/tex]
 
  • #22
Ah yes, that was my error. I missed the total mass parameter in the formula for the square inertia. :redface:

Thanks again arildno!
 
  • #23
I would like to emphasize/clarify a few points:

1)
You wanted to find the moment of inertia around a known centroid.
In my derivation, I made no simplifications based on that the origin was the centroid, hence, the moment of inertia formula derived is valid for rotation about an axis going through any interior point of the polygon that the polygon is "star-shaped" with respect to (more of that in 3).

Clearly, if the polygon is described in an arbitrary coordinate system, the vertex vectors are the relative vectors to the point through which the rotation axis passes.

2)
The formula cannot be used directly in calculating the moment of inertia around a point on the boundary of the polygon.
I suggest using the parallell axis theorem in that case.

3)
I must emphasize that the formula is not valid for an arbitrarily constructed polygon; the polygon must be "star-shaped" with regard to the point through which the rotation axis passes.
This means, in particular, that the full triangle drawn between between the rotation point and any 2 adjacent vertices is included in the polygon.

If you want to develop an algorithm for a fully general polygon, the best way is to split the polygon into a set of triangles.
 
  • #24
Sorry to bring the subject back up again, but if I where to write an algorithm for triangulated polygons how would I calculate the polygon’s inertia from the triangles?
 
  • #25
Calculate the moment of inertia around the center of masses of an individual triangle, and then utilize the parallell axis theorem to calculate the triangle's moment of inertia about the polygon's rotation axis. Adding these moments from all triangles gives the polygon's moment of inertia about the rotation axis.
 
  • #26
How do I actually go about and calculate the moment of inertia for an individual triangle, assuming there is a simpler method than using the algorithm for “star-shaped” polygons?
 
  • #27
Aidman said:
How do I actually go about and calculate the moment of inertia for an individual triangle, assuming there is a simpler method than using the algorithm for “star-shaped” polygons?

I'm not inclined to go through and do the work for you, but I expect that, if you assume a constant density, there will be a closed form for calculating the moment of inertia of a triangular prism about a particular axis.
 
  • #28
Sorry didn't meant for anyone to do the work for me, just point at a formula for calculating the inertia for a triangle around its center of mass. Because I have had no luck finding one.
 
  • #29
The closed form expression of the moment of inertia is:
[tex]I=\frac{M_{P}}{6A_{P}}\sum_{i=1}^{M}A_{T,i}(\vec{T}_{i,1}^{2}+\vec{T}_{i,1}\cdot\vec{T}_{i,2}+\vec{T}_{i,1}\cdot\vec{T}_{i,3}+\vec{T}_{i,2}\cdot\vec{T}_{i,3}+\vec{T}_{i,2}^{2}+\vec{T}_{i,3}^{2})[/tex]

where:
1. The polygon P (of mass [tex]M_{P}[/tex])consists of M triangles [tex]T_{i}[/tex]
2. [tex]\vec{T}_{i,s}[/tex] is the position vector of the s'th vertex in the i'th triangle (measured from the rotation axis)
3. [tex]A_{T,i}=\frac{1}{2}||(\vec{T}_{i,2}-\vec{T}_{i,1})\times(\vec{T}_{i,3}-\vec{T}_{i,1})||[/tex] is the area of triangle [tex]T_{i}[/tex]
4. [tex]A_{P}=\sum_{i=1}^{M}A_{T,i}[/tex]
 
Last edited:
  • #30
arildno said:
4. If you want to have a formula for the moment of inertia where the mass is explicit, here you have it:

[tex]I=\frac{m}{6}\frac{\sum_{n=1}^{N}||\vec{P}_{n+1}\times\vec{P}_{n}||(\vec{P}^{2}_{n+1}+\vec{P}_{n+1}\cdot\vec{P}_{n}+\vec{P}_{n}^{2})}{\sum_{n=1}^{N}||\vec{P}_{n+1}\times\vec{P}_{n}||}[/tex]
Sorry to comment on an old thread but i was just wondering why the ||\vec{P}_{n+1}\times\vec{P}_{n}|| is not canceled out in the equation seeing as it is in the numerator and denominator.

I'm trying to find out the moment of inertia of a regular pentagon with a side length of 265mm and mass of 42.9kg.

EDIT: So using this equation, I calculated the moment of inertia to be 0.835960338, which seems about right :)
 
Last edited:
  • #31
GINGERBEER said:
Sorry to comment on an old thread but i was just wondering why the ||\vec{P}_{n+1}\times\vec{P}_{n}|| is not canceled out in the equation seeing as it is in the numerator and denominator.

You can only cancel common factors in the numerator and denominator.

Just because a given term in the numerator happens to include a factor that is identical to a term in the numerator does not make this into a common factor.
 

1. What is the formula for calculating moment of inertia for polygons in a 2D system?

The formula for calculating moment of inertia for polygons in a 2D system is I = ∑mr², where I is the moment of inertia, m is the mass of the object, and r is the distance from the axis of rotation to the mass element.

2. How do I determine the axis of rotation for a polygon in a 2D system?

The axis of rotation for a polygon in a 2D system is typically located at the centroid of the polygon. This can be calculated by finding the average of the x and y coordinates of all the vertices of the polygon.

3. Can the moment of inertia be negative for polygons in a 2D system?

No, the moment of inertia cannot be negative for polygons in a 2D system. It is always a positive value as it represents the resistance of an object to changes in its rotational motion.

4. How does the shape of a polygon affect its moment of inertia in a 2D system?

The shape of a polygon has a significant impact on its moment of inertia in a 2D system. Generally, a larger moment of inertia indicates a greater resistance to changes in rotational motion. This means that polygons with a greater spread of mass away from the axis of rotation will have a larger moment of inertia.

5. Can I use the same formula to calculate moment of inertia for all polygons in a 2D system?

Yes, the formula for calculating moment of inertia is the same for all polygons in a 2D system. However, the values for mass and distance from the axis of rotation will vary depending on the specific shape and dimensions of the polygon.

Similar threads

Replies
2
Views
2K
Replies
1
Views
2K
Replies
1
Views
1K
Replies
2
Views
14K
  • Linear and Abstract Algebra
Replies
9
Views
1K
Replies
4
Views
4K
  • Introductory Physics Homework Help
2
Replies
40
Views
2K
  • Mechanical Engineering
Replies
11
Views
1K
  • Classical Physics
2
Replies
49
Views
2K
Back
Top