FM2 3.1 - Centres of mass by integration

FM2 3.1 - Centres of mass by integration

Integration extends mass-weighted averages to a continuous body. You will choose a mass element for a wire, lamina, solid or thin shell, integrate its mass moment, and use symmetry to remove unnecessary coordinates. You will also see why the same outline can give different centres when its density or its physical construction changes.

Add small mass moments by integration

Do equal-height bands near the top of a spherical shell contain less material than bands near its equator? Their circles are smaller, but their curved widths are larger. For a uniform thin spherical shell, equal axial heights actually give equal masses. You will build the integral that explains this at the end.

You already know

For discrete masses, Mxˉ=∑mixiM\bar x=\sum m_i x_i. A definite integral adds limiting contributions. Choose a reference coordinate before integrating, and retain the same reference when dividing mass moment by mass. The Core Pure disc method gives volume element πy2 dx\pi y^2\,\mathrm{d}x for rotation about the x-axis.

Divide a body into small pieces. If a piece has mass Δm\Delta m near coordinate xx, its mass moment is approximately xΔmx\Delta m. Summing and refining the pieces produces the continuous moment formula.

Continuous centre of mass

For a body with positive finite total mass,

M=∫dm,Mxˉ=∫x dm,xˉ=∫x dm∫dmM=\int\mathrm{d}m,\qquad M\bar x=\int x\,\mathrm{d}m,\qquad \bar x=\frac{\int x\,\mathrm{d}m}{\int\mathrm{d}m}

Use corresponding moment integrals for other coordinates. Choose dm\mathrm{d}m from the actual mass distribution, not just from the outline. Symmetry of that distribution may give coordinates directly.

Not in the formula booklet: learn it.

The mass element is the modelling decision. A density is mass per unit of the appropriate geometric measure.

ModelDensityMass element
Thin wireLinear density λ\lambdadm=λ ds\mathrm{d}m=\lambda\,\mathrm{d}s
Thin lamina or shellAreal density σ\sigmadm=σ dA\mathrm{d}m=\sigma\,\mathrm{d}A
SolidVolume density ρ\rhodm=ρ dV\mathrm{d}m=\rho\,\mathrm{d}V

These density relations are not printed in the formula booklet. On a straight rod along the x-axis, ds=dx\mathrm{d}s=\mathrm{d}x; on a curved wire its arc length must be used instead. A common constant density cancels from the centre calculation, but a position-dependent density stays inside both integrals.

Worked example 1: a rod denser at one end

Question. A thin straight rod occupies 0≤x≤20\le x\le2, with x measured in metres. Its linear density is k(1+x) kg m−1k(1+x)\,\mathrm{kg\,m^{-1}}, where k>0k>0 and x denotes the numerical coordinate in metres. Find its mass and centre of mass.

Plan. Integrate the linear density for mass and x times that density for mass moment.

M=k∫02(1+x) dx=k[x+x22]02=4k kgMxˉ=k∫02x(1+x) dx=k[x22+x33]02=14k3xˉ=14k/34k=76 m\begin{aligned} M&=k\int_0^2(1+x)\,\mathrm{d}x =k\left[x+\frac{x^2}{2}\right]_0^2=4k\,\mathrm{kg} \\ M\bar x&=k\int_0^2x(1+x)\,\mathrm{d}x =k\left[\frac{x^2}{2}+\frac{x^3}{3}\right]_0^2=\frac{14k}{3} \\ \bar x&=\frac{14k/3}{4k}=\frac76\,\mathrm{m} \end{aligned}

Why this step? Dividing by the rod's length would average a coordinate-density product incorrectly. Divide by the actual total mass.

Answer. Mass is 4k kg4k\,\mathrm{kg}, with centre at x=7/6 mx=7/6\,\mathrm{m}.

Check. The centre lies right of the midpoint x=1x=1, towards the denser end, and remains inside the rod. The scale factor kk cancels.

Change only the density gradient: replace 1+x1+x by 1+2x1+2x, retaining the same rod.

Worked example 2: a stronger density gradient

Question. The same rod has linear density k(1+2x) kg m−1k(1+2x)\,\mathrm{kg\,m^{-1}}. Find its mass and centre.

Plan. Recalculate both integrals using the changed density.

M=k∫02(1+2x) dx=k[x+x2]02=6k kgMxˉ=k∫02(x+2x2) dx=k[x22+2x33]02=22k3xˉ=22k/36k=119 m\begin{aligned} M&=k\int_0^2(1+2x)\,\mathrm{d}x =k[x+x^2]_0^2=6k\,\mathrm{kg} \\ M\bar x&=k\int_0^2(x+2x^2)\,\mathrm{d}x =k\left[\frac{x^2}{2}+\frac{2x^3}{3}\right]_0^2=\frac{22k}{3} \\ \bar x&=\frac{22k/3}{6k}=\frac{11}{9}\,\mathrm{m} \end{aligned}

Why this step? The density affects the mass denominator as well as the moment numerator. Changing it in only one integral would describe two different bodies.

Answer. M=6k kgM=6k\,\mathrm{kg}, xˉ=11/9 m\bar x=11/9\,\mathrm{m}.

Check. 11/9>7/611/9>7/6, consistent with the stronger bias towards the right end.

A zero density at an isolated endpoint is harmless if the body still has positive finite mass; a negative density over a region is not an admissible mass model.

Watch out

Wrong: cancelling a variable factor such as 1+x1+x from inside two different integrals. Right: integrate it in both. Only common constant factors may be cancelled before integration.

Change only the density gradient again.

A lamina strip has both area and its own centre

For a lamina under y=f(x)≥0y=f(x)\ge0, a vertical strip has area approximately f(x)Δxf(x)\Delta x. If its density depends on x alone, it is constant along that strip. Its centre is halfway up, at height f(x)/2f(x)/2. The x and y mass moments therefore use different factors.

Vertical strips in a lamina

For a≤x≤ba\le x\le b, 0≤y≤f(x)0\le y\le f(x), and areal density σ(x)\sigma(x) constant across each vertical strip,

M=∫abσ(x)f(x) dxM=\int_a^b\sigma(x)f(x)\,\mathrm{d}x Mxˉ=∫abxσ(x)f(x) dx,Myˉ=12∫abσ(x)[f(x)]2 dxM\bar x=\int_a^bx\sigma(x)f(x)\,\mathrm{d}x,\qquad M\bar y=\frac12\int_a^b\sigma(x)[f(x)]^2\,\mathrm{d}x

The one-half in the y moment comes from the strip's own centre. If density varies along the strip as well, its centre is not automatically halfway up; choose elements appropriate to that density.

Not in the formula booklet: learn it.

Sketch the strip before choosing an integrand. It makes the extra x in one moment and the half-height in the other visible.

[DIAGRAM: asset_name: A vertical lamina strip and its centre; asset_slug: 9fm0-fm153-lamina-strip; caption: The strip mass is density times f(x) dx, acting halfway up when density is constant vertically; description: The region under y equals x squared from zero to two has a narrow vertical strip at x. The strip height is f(x), width dx and centre at x f(x) over two. Horizontal and vertical coordinate arms to that centre are labelled. A side annotation distinguishes the x moment x dm and y moment f(x) dm over two.]
The strip mass is density times f(x) dx, acting halfway up when density is constant vertically

Worked example 3: a uniform curved plate

Question. A uniform lamina occupies 0≤x≤20\le x\le2, 0≤y≤x20\le y\le x^2, with coordinates in metres and constant areal density kk. Find its centre of mass.

Plan. Use vertical strips, calculating their mass, x moment and half-height y moment.

M=k∫02x2 dx=8k3Mxˉ=k∫02x3 dx=4kMyˉ=k2∫02x4 dx=16k5xˉ=48/3=32,yˉ=16/58/3=65\begin{aligned} M&=k\int_0^2x^2\,\mathrm{d}x=\frac{8k}{3} \\ M\bar x&=k\int_0^2x^3\,\mathrm{d}x=4k \\ M\bar y&=\frac{k}{2}\int_0^2x^4\,\mathrm{d}x=\frac{16k}{5} \\ \bar x&=\frac{4}{8/3}=\frac32,\qquad \bar y=\frac{16/5}{8/3}=\frac65 \end{aligned}

Why this step? A strip's whole height determines its area, but only half that height locates its centre. Using the top boundary as the strip centre doubles the y moment.

Answer. G=(3/2,6/5) mG=(3/2,6/5)\,\mathrm{m}.

Check. Both coordinates lie within the lamina's coordinate ranges. The region is wider in vertical extent near the right, consistent with xˉ>1\bar x>1.

Change only the density to kxkx, retaining the same boundary and limits. It still stays constant along each vertical strip.

Worked example 4: density varying across the same plate

Question. The same lamina has areal density kxkx, where k>0k>0 and x is the numerical coordinate in metres. Find its centre.

Plan. Include the extra x density factor in all three integrals.

M=k∫02x3 dx=4kMxˉ=k∫02x4 dx=32k5Myˉ=k2∫02x5 dx=16k3xˉ=85,yˉ=43\begin{aligned} M&=k\int_0^2x^3\,\mathrm{d}x=4k \\ M\bar x&=k\int_0^2x^4\,\mathrm{d}x=\frac{32k}{5} \\ M\bar y&=\frac{k}{2}\int_0^2x^5\,\mathrm{d}x=\frac{16k}{3} \\ \bar x&=\frac85,\qquad\bar y=\frac43 \end{aligned}

Why this step? The boundary height and the density are different ingredients. The y moment is still half the height times strip mass; the varying density does not remove the one-half.

Answer. G=(8/5,4/3) mG=(8/5,4/3)\,\mathrm{m}.

Check. Both coordinates increase towards the tall, denser right-hand part of the lamina.

Horizontal strips are equally valid when their lengths and centres are easier to express. Use the corresponding coordinate and density consistently; do not mix a vertical-strip mass with a horizontal-strip centre.

Watch out

Wrong: Myˉ=∫σy2 dxM\bar y=\int\sigma y^2\,\mathrm{d}x for these full-height strips. Right: include the factor 1/21/2. A moment integrand is coordinate of the element's centre multiplied by its mass.

Change only the density again.

Use circular slices for a uniform solid

Rotate a region about the x-axis. A disc at coordinate x has volume πy2 dx\pi y^2\,\mathrm{d}x and its centre lies on that axis at x. If density is uniform, rotational symmetry gives yˉ=zˉ=0\bar y=\bar z=0, leaving one coordinate to calculate.

Uniform solid of revolution

For a uniform solid formed by rotating 0≤y≤f(x)0\le y\le f(x), a≤x≤ba\le x\le b, about the x-axis,

V=π∫ab[f(x)]2 dx,Vxˉ=π∫abx[f(x)]2 dxV=\pi\int_a^b[f(x)]^2\,\mathrm{d}x,\qquad V\bar x=\pi\int_a^bx[f(x)]^2\,\mathrm{d}x

A common volume density cancels. The denominator is volume, not the area under the original curve. For washers, replace f(x)2f(x)^2 by outer-radius squared minus inner-radius squared.

Not in the formula booklet: learn it.

A profile diagram is a section through the solid, not a lamina whose centroid can be reused. Rotation changes the relative masses of the slices.

[DIAGRAM: asset_name: Disc slices for a solid's mass moment; asset_slug: 9fm0-fm153-solid-slices; caption: The disc area is pi times radius squared and its moment arm is the axial coordinate; description: The generating profile y equals square root x between zero and four is shown above and below the x-axis. A thin disc at x has radius square root x and axial thickness dx; the circle's area is pi x. The axial centre eight thirds is marked, and symmetry places the other coordinates on the axis.]
The disc area is pi times radius squared and its moment arm is the axial coordinate

Worked example 5: a uniform curved solid

Question. The region 0≤x≤40\le x\le4, 0≤y≤x0\le y\le\sqrt{x} is rotated about the x-axis to form a uniform solid. Coordinates are in metres. Find its volume and centre of mass.

Plan. Use disc volumes and multiply each by its x-coordinate for the axial moment.

V=π∫04x dx=8π m3Vxˉ=π∫04x2 dx=64π3xˉ=64π/38π=83 m\begin{aligned} V&=\pi\int_0^4x\,\mathrm{d}x=8\pi\,\mathrm{m^3} \\ V\bar x&=\pi\int_0^4x^2\,\mathrm{d}x=\frac{64\pi}{3} \\ \bar x&=\frac{64\pi/3}{8\pi}=\frac83\,\mathrm{m} \end{aligned}

Why this step? Each slice has mass proportional to its circular area πx\pi x, rather than its profile height x\sqrt{x}.

Answer. Volume is 8π m38\pi\,\mathrm{m^3}, with centre (8/3,0,0) m(8/3,0,0)\,\mathrm{m}.

Check. 8/3>28/3>2, towards the broader end; the other coordinates are zero by rotational mass symmetry.

Change only the upper endpoint, extending the same profile from x=4x=4 to x=6x=6.

Continue this lesson.

You’re previewing NovaLearn without an account. Create a free account to read the rest.

Create free account