next up previous
Next: Moment of inertia Up: Rotational motion Previous: The vector product


Centre of mass

The centre of mass--or centre of gravity--of an extended object is defined in much the same manner as we earlier defined the centre of mass of a set of mutually interacting point mass objects--see Sect. 6.3. To be more exact, the coordinates of the centre of mass of an extended object are the mass weighted averages of the coordinates of the elements which make up that object. Thus, if the object has net mass $M$, and is composed of $N$ elements, such that the $i$th element has mass $m_i$ and position vector ${\bf r}_i$, then the position vector of the centre of mass is given by
\begin{displaymath}
{\bf r}_{cm} = \frac{1}{M} \sum_{i=1,N} m_i {\bf r}_i.
\end{displaymath} (325)

If the object under consideration is continuous, then
\begin{displaymath}
m_i = \rho({\bf r}_i) V_i,
\end{displaymath} (326)

where $\rho({\bf r})$ is the mass density of the object, and $V_i$ is the volume occupied by the $i$th element. Here, it is assumed that this volume is small compared to the total volume of the object. Taking the limit that the number of elements goes to infinity, and the volume of each element goes to zero, Eqs. (325) and (326) yield the following integral formula for the position vector of the centre of mass:
\begin{displaymath}
{\bf r}_{cm} = \frac{1}{M} \int\!\int\!\int \rho {\bf r} dV.
\end{displaymath} (327)

Here, the integral is taken over the whole volume of the object, and $dV = dx dy dz$ is an element of that volume. Incidentally, the triple integral sign indicates a volume integral: i.e., a simultaneous integral over three independent Cartesian coordinates. Finally, for an object whose mass density is constant--which is the only type of object that we shall be considering in this course--the above expression reduces to
\begin{displaymath}
{\bf r}_{cm} = \frac{1}{V} \int\!\int\!\int {\bf r} dV,
\end{displaymath} (328)

where $V$ is the volume of the object. According to Eq. (328), the centre of mass of a body of uniform density is located at the geometric centre of that body.

Figure 72: Locating the geometric centre of a cube.
\begin{figure}
\epsfysize =1.5in
\centerline{\epsffile{cube.eps}}
\end{figure}

For many solid objects, the location of the geometric centre follows from symmetry. For instance, the geometric centre of a cube is the point of intersection of the cube's diagonals. See Fig. 72. Likewise, the geometric centre of a right cylinder is located on the axis, half-way up the cylinder. See Fig. 73.

Figure 73: Locating the geometric centre of a right cylinder.
\begin{figure}
\epsfysize =2in
\centerline{\epsffile{cyl.eps}}
\end{figure}

As an illustration of the use of formula (328), let us calculate the geometric centre of a regular square-sided pyramid. Figure 74 shows such a pyramid. Let $a$ be the length of each side. It follows, from simple trigonometry, that the height of the pyramid is $h=a/\sqrt{2}$. Suppose that the base of the pyramid lies on the $x$-$y$ plane, and the apex is aligned with the $z$-axis, as shown in the diagram. It follows, from symmetry, that the geometric centre of the pyramid lies on the $z$-axis. It only remains to calculate the perpendicular distance, $z_{cm}$, between the geometric centre and the base of the pyramid. This quantity is obtained from the $z$-component of Eq. (328):

\begin{displaymath}
z_{cm} = \frac{\int\!\int\!\int z dx dy dz}{\int\!\int\!\int dx dy dz},
\end{displaymath} (329)

where the integral is taken over the volume of the pyramid.

Figure 74: Locating the geometric centre of a regular square-sided pyramid.
\begin{figure}
\epsfysize =2.5in
\centerline{\epsffile{pyramid.eps}}
\end{figure}

In the above integral, the limits of integration for $z$ are $z=0$ to $z=h$, respectively (i.e., from the base to the apex of the pyramid). The corresponding limits of integration for $x$ and $y$ are $x,y=-a (1-z/h)/2$ to $x,y=+a (1-z/h)/2$, respectively (i.e., the limits are $x,y=\pm a/2$ at the base of the pyramid, and $x,y=\pm 0$ at the apex). Hence, Eq. (329) can be written more explicitly as

\begin{displaymath}
z_{cm} = \frac{\int_0^h z dz\!\int_{-a (1-z/h)/2}^{+a (1...
...a (1-z/h)/2} dy\!
\int_{-a (1-z/h)/2}^{+a (1-z/h)/2} dx}.
\end{displaymath} (330)

As indicated above, it makes sense to perform the $x$- and $y$- integrals before the $z$-integrals, since the limits of integration for the $x$- and $y$- integrals are $z$-dependent. Performing the $x$-integrals, we obtain
\begin{displaymath}
z_{cm} = \frac{\int_0^h z dz\!\int_{-a (1-z/h)/2}^{+a (1...
...h dz\!\int_{-a (1-z/h)/2}^{+a (1-z/h)/2} 
a (1-z/h) dy}.
\end{displaymath} (331)

Performing the $y$-integrals, we obtain
\begin{displaymath}
z_{cm} = \frac{\int_0^h a^2 z (1-z/h)^2 dz}{\int_0^h a^2 (1-z/h)^2 dz}.
\end{displaymath} (332)

Finally, performing the $z$-integrals, we obtain
\begin{displaymath}
z_{cm} = \frac{ a^2\left[z^2/2 - 2 z^3/(3 h) + z^4/(4 h^2...
... h)\right]_0^h} = \frac{a^2 h^2/12}{a^2 h/3} = \frac{h}{4}.
\end{displaymath} (333)

Thus, the geometric centre of a regular square-sided pyramid is located on the symmetry axis, one quarter of the way from the base to the apex.


next up previous
Next: Moment of inertia Up: Rotational motion Previous: The vector product
Richard Fitzpatrick 2006-02-02