We define, for a density function $\rho(x,y)$, the following to be the center of mass for a wire defined by the curve $C$.

$$\bar{x}= \frac{1}{M_\text{Total}}\int_Cx\rho(x,y)ds$$

$$\bar{y}= \frac{1}{M_\text{Total}}\int_Cy\rho(x,y)ds$$

Where does this formula come from? I'm not clear why we're integrating over the function $x\rho(x,y)$ or $y\rho(x,y)$ and why multiplying the result by the reciprocal of total mass would give us the coordinate for the center of mass.

  • $\begingroup$ What is $C{{}}$? $\endgroup$ – Angina Seng Oct 12 '17 at 5:18
  • $\begingroup$ @LordSharktheUnknown sorry, I should've clarified that. I edited my post to reflect that $C$ is the curve for which a given piece of wire is defined. $\endgroup$ – rb612 Oct 12 '17 at 5:19

Consider a lever...

Place two masses $m_1$ and $m_2$ on opposite ends of a lever of length $L$---where should the fulcrum be placed for this lever to balance?

Letting $d_1$, $d_2$ be the distances from the point masses to the fulcrum, Archimedes tells us that

$$m_1 d_1 = m_2 d_2.$$

In modern terms, we would say that the masses are applying equal and opposite torque to the system under the influence of gravity.

Introduce a coordinate system and let $x_1, x_2$ be the positions of the point masses and $\bar{x}$ be the position of the fulcrum. Our equation becomes \begin{align*}0 &= m_1 (x_1 - \bar{x}) + m_2 (x_2 - \bar{x}) \\ (m_1 + m_2) \bar{x} &= m_1 x_1 + m_2 x_2 \\ \bar{x} &= \frac{1}{m_1 + m_2} \left( m_1 x_1 + m_2 x_2 \right). \end{align*} We can repeat this argument for $n$ point masses $m_1, \dots, m_n$ located at $x_1, \dots, x_n$ to find that $$\bar{x} = \frac{1}{\sum_{i=1}^n m_i} \sum_{i=1}^n m_i x_i.$$

One-dimensional case

Now consider a continuous density distribution $\rho(x)$ along an interval $I=[a,b]$. Take a partition of $I$ by $a=x_0<x_1< \dots<x_{n-1}< x_n=b$ and approximate with a Riemann sum, noting that the mass on each sub-interval is approximately the density at a point times the length of each interval. In effect, we are approximating by collapsing the mass on each interval into a point mass.

$$\bar{x} \approx \frac{1}{\sum_{i=1}^n \rho(x_i)\cdot(x_i - x_{i-1})} \sum_{i=1}^n x_i \cdot\rho(x_i) \cdot(x_i - x_{i-1})$$

which becomes exact upon taking the limit $n \to \infty$:

$$\bar{x} = \frac{1}{\int_I \rho(x) ~\mathrm{d}x} \int_I x \rho(x) ~\mathrm{d}x = \frac{1}{M_{total}} \int_I x \rho(x) ~\mathrm{d}x .$$

Notice that we are integrating the $x$-coordinate against density and then dividing by the total mass to find the center of gravity.

On to your case

Instead of $\bar{x}$ being a point on which our one-dimensional mass rests, we view $\bar{x}$ as being a line on which our bent wire balances. In effect, the torque about this line from a point mass (under the effect of gravity) depends on the $x$-coordinate and the mass located there.

Parametrize the curve $C$ with $\mathbf{r}(t)=( x(t), y(t))$ with $t \in [a,b]$, letting $\mathbf{r}(t)$ be suitably nice, i.e., one-to-one and smooth. Again, take a partition $a=t_0 < t_1 < \dots < t_{n-1} < t_n = b$. The mass of wire associated with each sub-interval $[t_{i-1}, t_i]$ will be approximately the length $\Vert \mathbf{r}(t_i) - \mathbf{r}(t_{i-1}) \Vert$ times the density $\rho ( \mathbf{r}(t_i) )$, with an $x$-coordinate of approximately $x(t_i)$. This leads to the Riemann sum

\begin{align*}\bar{x} &\approx \frac{1}{\sum_{i=1}^n \rho(\mathbf{r}(t_i))\cdot \Vert \mathbf{r}(t_i) - \mathbf{r}(t_{i-1}) \Vert } \sum_{i=1}^n x(t_i)\cdot\rho(\mathbf{r}(t_i)) \cdot \Vert \mathbf{r}(t_i) - \mathbf{r}(t_{i-1})\Vert \\ &= \frac{1}{\sum_{i=1}^n \rho(\mathbf{r}(t_i))\cdot \frac{\Vert \mathbf{r}(t_i) - \mathbf{r}(t_{i-1}) \Vert}{t_i - t_{i-1}} } \sum_{i=1}^n x(t_i)\cdot\rho(\mathbf{r}(t_i)) \cdot \frac{\Vert \mathbf{r}(t_i) - \mathbf{r}(t_{i-1})\Vert}{t_i - t_{i-1}}. \end{align*}

Again, this becomes exact upon taking the limit $n \to \infty$:

$$\bar{x} = \frac{1}{\int_a^b \rho(\mathbf{r}(t)) \Vert \mathbf{r}'(t) \Vert ~\mathrm{d}t} \int_a^b x(t) \cdot \rho(\mathbf{r}(t))\cdot \Vert \mathbf{r}'(t) \Vert ~\mathrm{d}t.$$

Letting $\mathrm{d}s$ be the "arc-length element" $\Vert \mathbf{r}'(t) \Vert ~\mathrm{d}t$, we write this integral in the more familiar form

$$\bar{x} = \frac{1}{\int_C \rho(x,y) ~\mathrm{d}s} \int_C x \rho(x,y) ~\mathrm{d}s = \frac{1}{M_\text{total}} \int_C x \rho(x,y) ~\mathrm{d}s.$$

Repeating this argument for $\bar{y}$ gives

$$\bar{y} = \frac{1}{\int_C \rho(x,y) ~\mathrm{d}s} \int_C y \rho(x,y) ~\mathrm{d}s = \frac{1}{M_\text{total}} \int_C y \rho(x,y) ~\mathrm{d}s.$$

| cite | improve this answer | |
  • $\begingroup$ Nicely explained! (I took the liberty to change $*$ into $\cdot$.) $\endgroup$ – Mårten W Oct 12 '17 at 11:08

The intuition is that the first moment, in any direction, about the centre of mass $(\bar{x},\bar{y})$ should be zero. In the $x$-direction, this means $$ \int_C (x-\bar{x})\rho(x,y)\,ds = 0 \Leftrightarrow \bar{x} \! \underbrace{\int_C \rho(x,y)\,ds}_{=M_\text{Total}} = \int_C x\rho(x,y)\,ds \Leftrightarrow \bar{x} = \frac{1}{M_\text{Total}}\int_C x\rho(x,y)\,ds. $$ Similarly for the $y$-direction.

| cite | improve this answer | |

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.