6.04d - Centres of mass by integration

6.04d - Centres of mass by integration

Build a centre-of-mass calculation from a thin strip or disc. Find both coordinates of a uniform lamina, including a region between curves, and the axial centre of a uniform solid of revolution. Explain why the mass element changes when the same profile is turned into a solid.

A lamina as a sum of strips

A sheet has the profile y=x2y=x^2, between x=0x=0 and x=2x=2, above the horizontal axis. Turn that profile into a solid by rotating it about the axis. Will the sheet and the solid have the same horizontal centre of mass?

Near the wide end, the solid gains mass faster than the sheet does: a strip's area is proportional to its height, but a disc's area is proportional to the square of its radius. We will calculate both centres and measure that difference.

You already know

For components, Mxˉ=∑mixiM\bar x=\sum m_i x_i and Myˉ=∑miyiM\bar y=\sum m_i y_i. Uniform sheet mass is proportional to area. Definite integration adds the contributions of thin elements; integration by parts handles a product such as xexxe^x.

A vertical strip under y=f(x)≥0y=f(x)\ge0, of width Δx\Delta x, has approximate area f(x)Δxf(x)\Delta x. Its own centre is approximately (x,f(x)/2)(x,f(x)/2). To get a first moment, multiply the strip's mass by the coordinate of this centre, then add the moments.

[DIAGRAM: asset_name: Strip mass and strip centre; asset_slug: h245_lamina_vertical_strip; caption: The strip's centre is halfway up its height, while its mass is proportional to its full height; description: A coordinate diagram shades the region under y equals x squared from zero to two. A narrow vertical strip at x equals 1.3 has height f(x), width dx and midpoint at height f(x) divided by two. The strip midpoint and the complete lamina centre at (1.5,1.2) are separately labelled.]
The strip's centre is halfway up its height, while its mass is proportional to its full height

If the constant mass per unit area is σ\sigma, the strip contributes approximately σf(x)Δx\sigma f(x)\Delta x to the mass and σxf(x)Δx\sigma x f(x)\Delta x to the horizontal first moment. For the vertical first moment it contributes σ[f(x)/2]f(x)Δx\sigma[f(x)/2]f(x)\Delta x. Taking the thin-strip limit and cancelling σ\sigma gives the following result.

A uniform lamina above the horizontal axis

For the region 0≤y≤f(x)0\le y\le f(x), a≤x≤ba\le x\le b, with f(x)≥0f(x)\ge0, finite integrals and positive area AA,

A=∫abf(x) dxAxˉ=∫abxf(x) dxAyˉ=12∫ab[f(x)]2 dx\begin{aligned} A&=\int_a^b f(x)\,dx\\ A\bar x&=\int_a^b x f(x)\,dx\\ A\bar y&=\frac12\int_a^b[f(x)]^2\,dx \end{aligned}

Not in the formula booklet: learn it.

The factor 1/21/2 belongs to the strip's centroid height; it does not belong in the area or horizontal first moment.

Worked example 1: a parabolic lamina

Question. A uniform lamina occupies the region bounded by y=x2y=x^2, the xx-axis and x=2x=2. Determine its centre of mass.

Plan. Find the area, then the horizontal and vertical first moments using vertical strips.

A=∫02x2 dx=[x33]02=83strip areasAxˉ=∫02x3 dx=[x44]02=4multiply by strip xAyˉ=12∫02x4 dx=12[x55]02=165multiply by strip y/2xˉ=48/3=32,yˉ=16/58/3=65divide moments by area\begin{aligned} A&=\int_0^2x^2\,dx=\left[\frac{x^3}{3}\right]_0^2=\frac83 &&\text{strip areas}\\ A\bar x&=\int_0^2x^3\,dx=\left[\frac{x^4}{4}\right]_0^2=4 &&\text{multiply by strip }x\\ A\bar y&=\frac12\int_0^2x^4\,dx=\frac12\left[\frac{x^5}{5}\right]_0^2=\frac{16}{5} &&\text{multiply by strip }y/2\\ \bar x&=\frac4{8/3}=\frac32,\qquad \bar y=\frac{16/5}{8/3}=\frac65 &&\text{divide moments by area} \end{aligned}

Why this step? The lamina's height varies, so averaging the heights of the endpoints would not average the mass. Each strip's contribution must carry its own area.

Answer. G=(32,65)G=\left(\dfrac32,\dfrac65\right).

Check. The centre is to the right of x=1x=1, where the taller strips put more mass, and lies inside the region since 0<6/5<(3/2)20<6/5<(3/2)^2.

Change only the profile's power from x2x^2 to x3x^3, keeping the same interval. The right-hand strips become relatively more important.

Worked example 2: a cubic profile

Question. A uniform lamina occupies 0≤x≤20\le x\le2, 0≤y≤x30\le y\le x^3. Determine its centre of mass.

Plan. Use the new height in the area and both first moments.

A=∫02x3 dx=4new strip heightAxˉ=∫02x4 dx=325horizontal momentAyˉ=12∫02x6 dx=647vertical moment(xˉ,yˉ)=(85,167)divide each by 4\begin{aligned} A&=\int_0^2x^3\,dx=4 &&\text{new strip height}\\ A\bar x&=\int_0^2x^4\,dx=\frac{32}{5} &&\text{horizontal moment}\\ A\bar y&=\frac12\int_0^2x^6\,dx=\frac{64}{7} &&\text{vertical moment}\\ (\bar x,\bar y)&=\left(\frac85,\frac{16}{7}\right) &&\text{divide each by }4 \end{aligned}

Why this step? The vertical moment squares the complete height x3x^3, producing x6x^6; it does not merely increase its power by one.

Answer. G=(85,167)G=\left(\dfrac85,\dfrac{16}{7}\right).

Check. The horizontal centre moves from 3/23/2 to 8/58/5, towards the end containing relatively more mass. Also 0<16/7<(8/5)30<16/7<(8/5)^3.

For an exponential profile, the setup is unchanged; only the antiderivative changes. Keep the area visible as a separate quantity so you can divide both moments by the same value.

Watch out

Wrong: Ayˉ=∫y dxA\bar y=\int y\,dx. That integral is the area, not its first moment. Right: multiply the strip area y dxy\,dx by its centre height y/2y/2, giving Ayˉ=12∫y2 dxA\bar y=\frac12\int y^2\,dx.

Try another power, with a shorter interval.

Exponential profiles and integration by parts

A curved sign can have an exponential upper edge rather than a polynomial one. Its centre still comes from area and first moments; the horizontal moment now calls for integration by parts.

For the profile y=exy=e^x, a vertical strip has area ex dxe^x\,dx. Its horizontal first moment is xex dxxe^x\,dx. Choose the factor xx to differentiate, because one differentiation removes it:

∫xex dx=xex−∫ex dxintegration by parts=(x−1)ex+Cintegrate the remaining exponential\begin{aligned} \int xe^x\,dx&=xe^x-\int e^x\,dx &&\text{integration by parts}\\ &=(x-1)e^x+C &&\text{integrate the remaining exponential} \end{aligned}

Integration by parts is in the formula booklet; this particular antiderivative is not listed separately. The centre-of-mass strip formulas from the previous part are not in the booklet. Writing the limits in the moment integral makes the coordinate system and the piece of material unambiguous.

Worked example 3: an exponential sign

Question. A uniform lamina occupies 0≤x≤ln⁡30\le x\le\ln3, 0≤y≤ex0\le y\le e^x. Determine its centre of mass, giving exact coordinates.

Plan. Integrate for the area, use parts for the horizontal moment, and square the height for the vertical moment.

A=∫0ln⁡3ex dx=[ex]0ln⁡3=2total areaAxˉ=∫0ln⁡3xex dx=[(x−1)ex]0ln⁡3integration by parts=3(ln⁡3−1)−(−1)=3ln⁡3−2substitute both limitsAyˉ=12∫0ln⁡3e2x dx=14[e2x]0ln⁡3=2vertical first moment\begin{aligned} A&=\int_0^{\ln3}e^x\,dx=[e^x]_0^{\ln3}=2 &&\text{total area}\\ A\bar x&=\int_0^{\ln3}xe^x\,dx =[(x-1)e^x]_0^{\ln3} &&\text{integration by parts}\\ &=3(\ln3-1)-(-1)=3\ln3-2 &&\text{substitute both limits}\\ A\bar y&=\frac12\int_0^{\ln3}e^{2x}\,dx =\frac14[e^{2x}]_0^{\ln3}=2 &&\text{vertical first moment} \end{aligned}

Why this step? At the lower limit x=0x=0, (x−1)ex=−1(x-1)e^x=-1, so subtracting the lower value adds 11; dropping it moves the predicted centre.

Answer. G=(3ln⁡3−22,1)G=\left(\dfrac{3\ln3-2}{2},1\right).

Check. The horizontal coordinate is approximately 0.6480.648, between 00 and ln⁡3\ln3, and to the right of the interval midpoint, as expected for increasing strip heights.

Change only the right-hand boundary from x=ln⁡3x=\ln3 to x=ln⁡5x=\ln5. The antiderivatives stay the same, but the new end contributes taller strips.

Worked example 4: extending the same sign

Question. A uniform lamina occupies 0≤x≤ln⁡50\le x\le\ln5, 0≤y≤ex0\le y\le e^x. Determine its centre of mass, giving exact coordinates.

Plan. Evaluate the same area and moment integrals at the new upper limit.

A=[ex]0ln⁡5=4new areaAxˉ=[(x−1)ex]0ln⁡5=5ln⁡5−4new horizontal momentAyˉ=14[e2x]0ln⁡5=6new vertical moment(xˉ,yˉ)=(5ln⁡5−44,32)divide by 4\begin{aligned} A&=[e^x]_0^{\ln5}=4 &&\text{new area}\\ A\bar x&=[(x-1)e^x]_0^{\ln5} =5\ln5-4 &&\text{new horizontal moment}\\ A\bar y&=\frac14[e^{2x}]_0^{\ln5}=6 &&\text{new vertical moment}\\ (\bar x,\bar y)&=\left(\frac{5\ln5-4}{4},\frac32\right) &&\text{divide by }4 \end{aligned}

Why this step? Extending the lamina changes the denominator as well as the moment; reusing the old area would not describe the new body.

Answer. G=(5ln⁡5−44,32)G=\left(\dfrac{5\ln5-4}{4},\dfrac32\right).

Check. Its horizontal centre is approximately 1.011.01, farther right than before, and its height has increased from 11 to 3/23/2 as material was added towards the tall end.

Even when a calculator can evaluate the integral, “Show that” requires the link between the mass model, the integral and the exact answer.

Exam tip

Keep the area, moment integral and substituted endpoint values visible. In a given-answer calculation, writing only the required expression for xˉ\bar x does not establish it.

Extend the same profile once more, using the lower endpoint as carefully as the upper one.

A lamina between two curves

When a strip begins above the axis, its centre is no longer halfway up its height measured from the axis. Locate its two endpoints before forming either moment.

Suppose the upper curve is y=u(x)y=u(x) and the lower is y=l(x)y=l(x). The strip height is u−lu-l, while its midpoint height is (u+l)/2(u+l)/2. Their product is (u2−l2)/2(u^2-l^2)/2, so the vertical first moment subtracts squares of heights. It does not square their difference.

Continue this lesson.

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

Create free account