\(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\)

MathWiki.png
 
MathWiki: The Dynamic Mathematics E-textbook > Calculus > Vector Calculus > Multiple Integrals > Moments and Centers of Mass

Moments and Centers of Mass

This section shows how to calculate the masses and moments of two- and three- dimensional objects in Cartesian \((x,y,z)\) coordinates.

Mass

We saw before that the double integral over a region of the constant function 1 measures the area of the region. If the region has uniform density 1, then the mass is the density times the area which equals the area. What if the density is not constant. Suppose that the density is given by the continuous function 

\[\text{Density}  =  \rho(x,y)\]

In this case we can cut the region into tiny rectangles where the density is approximately constant. The area of mass rectangle is given by 

Mass  =  (Density)(Area)  =  \(\rho (x,y) \Delta{x} \Delta{y}\)

You probably know where this is going.  If we add all to masses together and take the limit as the rectangle size goes to zero, we get a double integral.

 

Definition: Mass

Let \(\rho(x,y)\) be the density of a lamina (flat sheet) \(R\) at the point \((x,y)\). Then the total mass of the lamina is the double integral

\[ \text{Mass} = \iint \rho (x,y)\, dydx\]

Three-Dimensional Solids

We define \(\rho (x,y,z)\) to be the density (mass per unit volume) of an object occupying a region \(D\) in space. The integral over \(D\) gives us the mass of the object. To see why, imagine partitioning the object into n mass elements. And when summing these mass elements up, it is the total mass.

\[M= \lim_{n\rightarrow\infty}\sum_{k=1}^n \Delta m_k =\lim_{n\rightarrow\infty}\sum_{k=1}^n \delta (x_k,y_k,z_k)\Delta V_k=\iiint_{a}^{b}\delta (x, y,z) dV\]

The integral of \(\rho (x,y,z)\) gives us the mass of the object. 
\[M=\iiint_{a}^{b}\delta dV\]

where \(\rho (x,y,z)\) is the density at point \((x,y,z)\).

Two-Dimensional Plates

\[M=\iint_{a}^{b}\delta dA\]

Example 1

A rectangular metal sheet with \(2 < x < 5\)  and \(0 < y < 3\) has density function 

\[\rho(x,y)  =  x + y\]

Set up the double integral that gives the mass of the metal sheet.

SOLUTION

We just have the integral

\[ \int_2^5 \int _0^3 (x+y)\, dydx\]

Moments and Center of Mass

We have seen in first year calculus that the moments about an axis are defined by the product of the mass times the distance from the axis.

        Mx  =  (Mass)(y)        My  =  (Mass)(x)

If we have a region \(R\) with density function \(\rho (x,y)\), then we do the usual thing. We cut the region into small rectangles for which the density is constant and add up the moments of each of these rectangles. Then take the limit as the rectangle size approaches zero. This will give us the total moment.

 

Definition 2a: Moments of Mass and Center of Mass

Suppose that \(\rho (x,y)\) is a continuous density function on a lamina \(R\). Then the moments of mass are 

\[ M_x = \int_0^1\int_0^2 k(x^2+y^2) y\, dy dx\]

and

\[ M_y = \int_0^1\int_0^2 k(x^2+y^2) x\, dy dx\]

and if \(M\) is the mass of the lamina, then the center of mass is 

\[ (\bar{x},\bar{y})=\left ( \dfrac{M_y}{M},\dfrac{M_x}{M}\right )\]

Example

Set up the integrals that give the center of mass of the rectangle with vertices (0,0), (1,0),(1,1), and (0,1) and density function proportional to the square of the distance from the origin. Use a calculator or computer to evaluate these integrals.

SOLUTION

The mass is given by 

\[M = \int_0^1\int_0^2 k(x^2+y^2)\, dy dx = \dfrac{2k}{3}\]

The moments are given by (definition 2a):

\[ M_x = \int_0^1\int_0^2 k(x^2+y^2) y\, dy dx\]

and

\[ M_y = \int_0^1\int_0^2 k(x^2+y^2) x\, dy dx\]

These evaluate to 

\[M_x  =  \dfrac{5k}{12}\]

and      

\[M_y  =  \dfrac{5k}{12}\]

It should not be a surprise that the moments are equal since there is complete symmetry with respect to \(x\) and \(y\). Finally, we divide to get

        (x,y)  =  (5/8,5/8)

This tells us that the metal plate will balance perfectly if we place a pin at (5/8,5/8)

Moments of Inertia

We often call \(M_x\) and \(M_y\) the first moments. They have first powers of \(y\) and \(x\) in their definitions and help find the center of mass. We define the moments of inertia (or second moments) by introducing squares of \(y\) and \(x\) in their definitions. The moments of inertia help us find the kinetic energy in rotational motion. Below is the definition

Definition: Moments of Intertia

Suppose that \(\rho (x,y)\) is a continuous density function on a lamina \(R\). Then the moments of inertia are

\[I_x = \iint_R \rho(x,y) y^2 \, dydx\]

\[I_y = \iint_R \rho(x,y) x^2 \, dydx\]

Exercise 

Find the moments of inertia for the previous square metal plate.

First moment

The first moment of a 3-D solid region \(D\) about a coordinate plane is defined as the triple integral over \(D\) of the distance from a point \((x,y,z)\) in \(D\) to the plane multiplied by the density of the solid at that point. First moments about the coordinate planes: 

\[M(yz)=\iiint_{a}^{b}\delta x\, dV\]
\[M(xz)=\iiint_{a}^{b}\delta y\,dV\]
\[M(xy)=\iiint_{a}^{b}\delta z \,dV\]

The first moment about the y-axis is the double integral over the region \(R\) forming the 2-D plate of the distance from the axis multiplied by the density.

\[M(y)=\iint_{a}^{b}\delta x\; dV\]
\[M(x)=\iint_{a}^{b}\delta y\; dV\]

Center of mass

We defined center of mass located in \(\bar{x}, \bar{y}, \bar{z}\). Then it is found from the first moments:

\[\bar{x} =\dfrac{M(y)}{M}\]
\[\bar{y} =\dfrac{M(x)}{M}\]

 

Example

Find the center of mass of a solid of constant density below by the disk \(R: Xl +y2 :s 4\) in the plane \(z = 0\) and above by the paraboloid \(z=4- x^2- y^2\). 

\[ \text{Total Mass} = \int_0^1 \int_{x^2}^x 1 \ dy \ dx \]

\[ = \int_0^1 \left. y \right |_{y=x^2}^{y=x} \ dx \]

\[ = \int_0^1 (x-x^2) \ dx \]

\[ = \left(\frac{x^2}{2}-\frac{x^3}{3} \right) |_0^1 \]

\[= \frac{1}{6}\]

Contributors

  • Shengqiao Luo (UCD)

You must to post a comment.
Last Modified
13:36, 29 Jan 2014

Page Rating

Was this article helpful?

Tags

This page has no custom tags set.
Module Vet Level:
Module Target Level:

Creative Commons License UC Davis GeoWiki by University of California, Davis is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 3.0 United States License. Permissions beyond the scope of this license may be available at copyright@ucdavis.edu. Terms of Use