Moments and center of mass

Videos

Review Videos

Videos, Math Dr. Bob:

Link to original

03 Theory

Theory 1

Moment

The moment of a region to an axis is the total (integral) of mass times distance to that axis:

Moment to x-axis:

Mx=ρydA(general formula)Mx=cdρy(g2(y)g1(y))dy(region between functions g2 and g1)

Moment to y-axis:

My=ρxdA(general formula)My=abρx(f2(x)f1(x))dx(region between functions f2 and f1)

Notice the swap in letters

  • My integrand has x factor
  • Mx integrand has y factor

Notice the total mass

If you remove x or y factors from the integrands, the integrals give total mass M.

These formulas are obtained by slicing the region into rectangular strips that are parallel to the axis in question.

The area per strip is then:

  • f(x)dx — region under y=f(x)
  • (f2f1)dx — region between f1 and f2
  • g(y)dy — region ‘under’ x=g(y)
  • (g2g1)dy — region between g1 and g2

center

center


The idea of moment is related to:

  • Torque balance and angular inertia
  • Center of mass

The center of mass (CoM) of a solid body is a single point with two important properties:

  1. CoM= “average position” of the body
    • The average position determines an effective center of dynamics. For example, gravity acting on every bit of mass of a rigid body acts the same as a force on the CoM alone.
  2. CoM= “balance point” of the body
    • The net torque (rotational force) about the CoM, generated by a force distributed over the body’s mass, equivalently a force on the CoM, is zero.

Centroid

When the body has uniform density, then the CoM is also called the centroid.


Center of mass from moments

Coordinates of the CoM:

x=Mym,y=Mxm

Here M is the total mass of the body.

Link to original

04 Illustration

Example - CoM of a parabolic plate

CoM of a parabolic plate

Find the CoM of the region depicted:

center

Solution

(1) Compute the total mass:

Area under the curve with density factor ρ:

m=02ρx2dxρx33|028ρ3

(2) Compute My:

Formula:

My=abρxdA

Interpret and calculate:

My=02ρxf(x)dxρ02x3dx4ρ=My

(3) Compute Mx:

Formula:

Mx=cdρydA

Width of horizontal strips between the curves:

w(y)=2y

Interpret dA:

dA=(2y)dy

Calculate integral:

Mx=cdρydA04ρy(2y)dy04ρy(2y)dy04ρ2ydy04ρy3/2dy16ρ5=Mx

(4) Compute CoM coordinates from moments:

CoM formulas:

x=Mym,y=Mxm

Insert data:

x=4ρ8ρ/332,y=16ρ/58ρ/365CoM=(x,y)=(32,65) Link to original

05 Theory

Theory 2

A downside of the technique above is that to find Mx we needed to convert the region into functions in y. This would be hard to do if the region was given as the area under a curve y=f(x) but f1(y) cannot be found analytically. An alternative formula can help in this situation.

Midpoint of strips for opposite variables

When the region lies between f1(x) and f2(x), we can find Mx with an x-integral:

Mx=abρ12(f22f12)dx(region between f1 and f2)

When the region lies between g1(y) and g2(y), we can find My with a y-integral:

My=cdρ12(g22g12)dy(region between g1 and g2)

Region under a curve

For the region “under the curve” y=f(x), just set:

f1=0,f2=f

For the region “under the curve” x=g(x), set:

g1=0,g2=g

The idea for these formulas is to treat each vertical strip as a point concentrated at the CoM of the vertical strip itself.

center

The height to this midpoint is 12f(x), and the area of the strip is f(x)dx, so the integral becomes ρ12f(x)2dx.

Link to original

06 Illustration

Example - Computing CoM using only vertical strips

Computing CoM using only vertical strips

Find the CoM of the region:

center

Solution

(1) Compute the total mass m:

Area under the curve times density ρ:

0π/2ρcosxdx=ρsinx|0π/2=ρ

(2) Compute Mx, using vertical strips and the ‘midpoints’ method:

Plug f2(x)=f(x)=cosx and f1(x)=0 into formula:

Mx=0π/2ρ12f22dx0π/2ρ12cos2xdx

Integration by ‘power to frequency conversion’:

0π/2ρ12cos2xdxρ40π/2(1+cos2x)dxρ4x|0π/2+ρsin2x8|0π/2=πρ8

(3) Compute My using vertical strips and the regular method:

Plug f(x)=cosx into formula:

My=abρxf(x)dx0π/2ρxcosxdx

Integration by parts. Set u=x, v=cosx and so u=1, v=sinx:

0π/2ρxcosxdxρxsinx|0π/2ρ0π/2sinxdxπρ21ρ(cosπ2cos0)=ρ(π21)

(4) Compute CoM:

CoM formulas:

x=Mym,y=Mxm

Plug in data:

x=ρ(π/21)ρπ21y=πρ/8ρπ8CoM=(x,y)=(π21,π8) Link to original

07 Theory

Theory 4

CoM of a triangle: “ h/3 trick”

Given a triangle with edge b and altitude h (measured along a perpendicular to b), the CoM of the triangle lies on the line parallel to b and shifted upwards by h/3.

For example, if the triangle is positioned with base b on the x-axis, then its CoM will lie on the line y=h/3.

Link to original

07 Theory

Theory 3

Two useful techniques for calculating moments and (thereby) CoMs:

  • Additivity principle
  • Symmetry

Additivity says that you can add moments of parts of a region to get the total moment of the region (to a given axis).

A symmetry principle is that if a region is mirror symmetric across some line, then the CoM must lie on that line.

Link to original

08 Illustration

Example - Center of mass using moments and symmetry

Center of mass using moments and symmetry

Find the center of mass of the two-part region:

center

Solution

(1) Symmetry: CoM on y-axis

Because the region is symmetric in the y-axis, the CoM must lie on that axis. Therefore x=0.


(2) Additivity of moments:

Write Mx for the total x-moment (distance measured to the x-axis from above).

Write Mxtri and Mxcirc for the x-moments of the triangle and circle.

Additivity of moments equation:

Mx=Mxtri+Mxcirc

(3) Find moment of the circle Mxcirc using Symmetry:

By symmetry we know xcirc=0.

By symmetry we know ycirc=5.

Area of circle with r=2 is A=4π, so total mass is M=4πρ.

Centroid-from-moments equation:

ycirc=Mxcircm

Solve the equation for Mxcirc:

ycirc=Mxcircm5=Mxcirc4πρMxcirc=20πρ

(4) Find moment of the triangle Mxtri using Symmetry and the “h/3 trick”:

By Symmetry, the CoM of the triangle must lie on the vertical line x=0.

By the “h/3 trick,” the CoM of the triangle must lie on the horizontal line at y=h/3. In this case h=3, so it lies on the line y=1.

Therefore the triangle CoM is:

(xtri,ytri)=(0,1)

To get moments, multiply by mass m=ρ12(4)(3)=6ρ:

Mytri=0,Mxtri=(6ρ)(1)=6ρ

(5) Apply additivity:

Mx=Mxtri+Mxcircρ(20π+6)

(6) Total mass of region:

Area of circle is 4π. Area of triangle is 1243=6. Thus m=ρA=ρ(4π+6).


(7) Compute center of mass y from total Mx and total M:

We have Mx=ρ(20π+6) and M=ρ(4π+6). Plug into formula:

y=Mxmρ(20π+6)ρ(4π+6)3.71CoM=(x,y)=(0,3.71) Link to original

Example - Center of mass - two part region

Center of mass - two part region

Find the center of mass of the region which combines a rectangle and triangle (as in the figure) by computing separate moments. What are those separate moments? Assume the mass density is ρ=1.

center

Solution

(1) Apply symmetry to rectangle:

By symmetry, the center of mass of the rectangle is located at (1,2).

Thus xrect=1 and yrect=2.


(2) Find moments of the rectangle:

Total mass of rectangle =mrect=ρ×area=18=8. Thus:

xrect=MyrectmrectMyrect=8yrect=MxrectmrectMxrect=16

(3) Find moments of the triangle:

CoM of triangle using “h/3 trick:”

From the base on the x-axis, we have h=4 so the CoM lies on the line y=4/3.

From the base on the y-axis, we have h=4 so the CoM lies on the line x=4/3.

Multiply by the triangle’s mass to get its moment:

mtri=ρ12(4)(4)8

Therefore:

Mxtri=(8ρ)(4/3)323Mytri=(8ρ)(4/3)323

(4) Add up total moments:

General formulas: Mx=Mxtri+Mxrect and My=Mytri+Myrect

Plug in data: Mx=323+16=803 and My=3238=83


(5) Find center of mass from moments:

Total mass of triangle =Mtri=ρ×area=11244=8.

Total combined mass =m=mtri+mrect=8+8=16.

Apply moment relation:

x=Mym8/31616y=Mxm80/31653CoM=(x,y)=(16,53)Link to original

Improper integrals

Videos

Review Videos

Videos, Math Dr. Bob:

Link to original

03 Theory

Theory 1

Improper integrals are those for which either a bound or the integrand itself become infinite somewhere on the interval of integration.

Examples:

(a)11x2dx,(b)021xdx,(c)1+11x2dx

(a) the upper bound is (b) the integrand goes to as x0+ (c) the integrand is at the point 0[1,1]

The limit interpretation of (a) is this:

11x2dx=limR1R1x2dx

The limit interpretation of (b) is this:

021xdx=limR0+R21xdx

The limit interpretation of (c) is this:

1+11x2dx=101x2dx+0+11x2dx=limR01R1x2dx+limR0+R+11x2dx

These limits are evaluated using the usual methods.

An improper integral is said to be convergent or divergent according to whether it may be assigned a finite value through the appropriate limit interpretation.

For example, (a) converges while (b) diverges.

Link to original

04 Illustration

Example - Improper integral - infinite bound

Improper integral - infinite bound

Show that the improper integral 2dxx3 converges. What is its value?

Solution

(1) Improper integral definition:

2dxx3limR2Rdxx3

(2) Replace infinity with a new symbol R:

2Rdxx312x2|2R1812R2

(3) Take limit as R:

limR2Rdxx3limR1812R218converges Link to original

Improper integral - infinite integrand

Improper integral - infinite integrand

Show that the improper integral 09dxx converges. What is its value?

Solution

(1) Improper integral definition:

09dxxlimR0+R9dxx

(2) Switch 0 to R and integrate:

R9dxxR9x1/2dx2x+1/2|R962R

(3) Take limit as R0+:

limR0+R9dxxlimR0+62R6converges Link to original

Example - Improper integral - infinity inside the interval

Example - Improper integral - infinity inside the interval

Does the integral 1+11xdx converge or diverge?

Solution

WRONG APPROACH:

It is tempting to compute the integral incorrectly, like this:

1+11xdxln|x||1+1ln|2|ln|2|=0

But this is wrong! There is an infinite integrand at x=0. We must break it into parts!


RIGHT APPROACH:

(1) Break apart at x=0 discontinuity:

1+11xdx101xdx+0+11xdx

(2) Improper integral definition, both terms separately:

limR01R1xdx+limR0+R+11xdx

(3) Integrate:

1R1xdxln|R|ln|1|ln|R|R+11xdxln|1|ln|R|lnR

(4) Limits:

limR0ln|R|=,limR0+lnR=+diverges

Neither limit is finite. For 1+11xdx to exist we’d need both of these limits to be finite. So the original integral diverges.

Link to original

05 Theory

Theory 2

Two tools allow us to determine convergence of a large variety of integrals. They are the comparison test and the p-integral cases.

Comparison test - integrals

The comparison test says:

  • When an improper integral converges, every smaller integral converges.
  • When an improper integral diverges, every bigger integral diverges.

Here, smaller and bigger are comparisons of the integrand at all values (accounting properly for signs), and the bounds are assumed to be the same.

For example, 2dxx3 converges, and x4>x3 implies 1x4<1x3 (when x>1), therefore the comparison test implies that 2dxx4 converges.

p-integral cases

Assume p>0 and a>0. We have:

p>1:adxxpconvergesand0adxxpdivergesp<1:adxxpdivergesand0adxxpconvergesp=1:adxxdivergesand0adxxdiverges

Proving the p-integral cases

It is easy to prove the convergence / divergence of each p-integral case using the limit interpretation and the power rule for integrals. (Or for p=1, using 1xdx=lnx+C.)

Additional improper integral types

The improper integral af(x)dx also has a limit interpretation:

af(x)dx=limRRaf(x)dx

The double improper integral f(x)dx has this limit interpretation:

f(x)dx=limRRaf(x)dx+limRaRf(x)dx

Where a is any finite number. This double integral does not exist if either limit does not exist for any value of a.

Double improper is not simultaneous!

Watch out! This may happen:

f(x)dxlimRRRf(x)dx

This simultaneous limit might exist only because of internal cancellation in a case where the separate individual limits do not exist! We do not say ‘convergent’ in these cases!

Link to original

06 Illustration

Example - Comparison to p-integrals

Comparison to p-integrals

Determine whether the integral converges:

(a) 2x3x41dx

(b) 11x2+x+1dx

Solution

(a) (1) Observe large x tendency:

Consider large values. Notice the integrand tends toward 1/x for large x.

x3x41x3x4forx,andx3x4=1x

(2) Try comparison to 1/x:

x3x41>?1x

Validate. Notice x41>0 and x>0 when x2.

x3x41>?1xx3x>?1(x41)x4>x41

(3) Apply comparison test:

We know:

x3x41>1x,21xdxdiverges

We conclude:

2x3x41dxdiverges

(b) (1) Observe large x tendency:

Consider large values. Notice the integrand tends toward 1/x2 for large x.

1x2+x+11x2forx

(2) Try comparison to 1/x2:

1x2+x+1<?1x2

Validate. Notice x2+x+1>0 and x2>0 when x1.

1x2+x+1<?1x21x2<?1(x2+x+1)x2<x2+x+1

(3) Apply comparison test:

We know:

1x2+x+1<1x2,11x2dxconverges

We conclude:

11x2+x+1dxconvergesLink to original