Is my method here for solving this 1D 3-body collision problem correct?

I was attempting to solve a three body problem which I made up that goes like this:

The diagram shows their velocities and masses just before collision.

The coefficient of restitution between each of the objects is 0.5.

The task is to work out their velocities after collision.I got answers, although I'm not sure if my method was correct.

I started off by calculating the total momentum, which was:

$10*20 + 5*15 + 4*(-25) = 175$ kgm/s

So obviously the total momentum after collision must be 175kgm/s.

From this observation, it can only mean that object 3 will move to the right, otherwise momentum couldn't possibly be conserved. This therefore implies that object 3 will have the greatest velocity (not necessarily in magnitude, although if we consider it vectorially with positive and negative numbers), of the three objects.

So we'd have something like this:

Whatever the velocities of objects 1 and 2 are, we know that:

• The velocity of object two will be less than the velocity of object 3
• The velocity of object 1 will be less than the velocity of object 2
• So $v_{3}>v_{2}>v_{1}$

We know that since momentum will be conserved: $10v_{1} + 5v_{2} + 4v_{3} = 175$

We know from Newton's Law of Restitution that:

$\frac{Speed\space of\space separation} {Speed\space of\space Approach}= e$

Now this is the bit I'm not sure about:

If it were just two objects, fine, I'd work out the speed of approach, then I would show the difference between the final speeds (which I could correctly format since I know which velocity is larger)

But there was three objects, so my technique was:

• Use the coefficient of resitution formula for each pair of objects, even for object 1 and object 3 which are not touching, to then form a triple simultaneous equation
• I'd then use the formula I have for the conservation of momentum to then solve for $v_{1}$ , $v_{2}$ and $v_{3}$

The bit I'm not sure about is whether it's valid to directly compare objects 1 and 2 since they are not touching.

This is what I mean:

For objects 1 and 2:

$\frac{v_{2}-v_{1}}{5}= 0.5$

Rearranging we then get

$v_{2}-v_{1} = 2.5$

I then did the same for objects 2 and 3 which seems fine enough:

$\frac{v_{3}-v_{2}} {40} = 0.5$

Rearranging to then get

$v_{3}-v_{2} = 20$

Finally, this part I wasn't sure about: Do I then proceed to do the same thing for objects 1 and 3 even though they don't directly touch (assuming the method I'm using is correct to begin with)? Or do I then solve the set of equations I already have.

I then proceeded with the first idea to end up with $\frac{v_{3}-v_{1}}{45} = 0.5$

Then rearranging to get $v_{3}-v_{1} = 22.5$

Long story short, through substitution into the momentum equation and rearranging we get:

$v_{1} = \frac{+145}{38}$ m/s

$v_{2} = \frac{+120}{19}$ m/s

$v_{3} = \frac{+500}{19}$ m/s

Which seems to look right (you'd expect object 3 to go veering off to the right)

But I don't know if this method is valid, since I made it up, and it doesn't feel rigorous.

So is my method correct? If not, how would I solve this and what are the answers?

Edit: the collisions are all simultaneous, units don't change. I initially accidentally mixed up the formula for the coefficient of restitution, although I have updated my answer accordingly. Momentum is conserved, and kinetic energy lost (since the coefficient of restitution is less than 1).

• It would help if you would state the initial positions of the three objects. With that information, readers can determine the order of the collisions. It would also help if you showed units on all numbers. Commented Dec 10, 2016 at 16:42
• @David White it's a simultaneous collision, all at the same time. i left out the units later on since is showed them to be m/s earlier on, Commented Dec 10, 2016 at 16:48
• This can be checked with a PhET simulation. PhET's answers for final velocities: v1 = 4.17 m/s to the right; v2 = 12.22 m/s to the right; v3 = 18.06 m/s to the right. Commented Dec 10, 2016 at 17:21
• @David White Explanation? Commented Dec 10, 2016 at 18:13
• My explanation is that your method appears to be invalid. I set up equations in the center-of-mass reference frame for conservation of momentum and loss of kinetic energy while accounting for the coefficient of restitution, and I note that I still have one degree of freedom in the equations. I suspect that PhET iteratively solved this problem for two objects at a time until all velocities were constant, and I have some degree of confidence in the PhET solution. Commented Dec 10, 2016 at 18:38

3 Answers

The order in which the collisions occur is important. The outcome is different if 1&2 collide first or if 2&3 collide first. In your scenario both collisions occur simultaneously. But collisions between rigid bodies are instantaneous, so that one of the 2 collisions will happen first, if we look on a small enough time scale. Therefore I think that a simulataneous collision, in which 3 or more rigid objects are in contact simultaneously, must be ambiguous. Because of this I think your method is not valid.

The fault is the rigid-body model, which implies that collisions are instantaneous. Real bodies are deformable to some extent, and deformation takes a finite time. Unlike in the rigid body model, it is possible for a 2nd collision to start before a 1st has ended, with the result that 3 colliding bodies can be in contact simultaneously. Forces can then be transmitted between bodies 1&3 through body 2, at the speed of sound in body 2. If this transit time is less than the duration of each collision, the 2 collisions will affect each other.

A simultaneous collision between more than 2 rigid objects is indeterminate, even in 1D. There are 3 unknown speeds in the outcome, but only 2 constraining equations (conservation of momentum, the restitution law). This indeterminacy also arises when considering forces in static structures, eg Calculating distribution of force in a N-body system of balls.

The equations governing the collision between 2 rigid bodies are
$m_1(u_1-v_1)+m_2(u_2-v_2)=0$ ... (conservation of momentum)
$v_2-v_1=-e(u_2-u_1)$ ... (the restitution law)
Combining these we get $$v_1=\frac{(m_1-em_2)u_1+(1+e)m_2u_2}{m_1+m_2}, v_2=\frac{(1+e)m_1u_1+(m2-em_1)u_2}{m_1+m_2}$$

As suggested by David White, I used PhET Collision Lab to simulate the collisions between your objects. If 1&2 collide first the final velocities are $(3.75,14.167,16.667)$ whereas if 2&3 collide first they are $(4.167,12.222,18.056)$. In both cases there are 3 collision events.

Using your method the equations are
$10v_1+5v_2+4v_3=175$
$v_2-v_1=\frac52$
$v_3-v_2=20$
with the resulting velocities $(3.816, 6.316, 26.316)$.

Although the outcome of each constituent 2-body collision can be determined, the overall outcome of a 3-body collision is difficult to determine because of the possibility of multiple collisions, especially when the middle body is the lightest. The special case of mass ratios $(1,\alpha,1)$ with $\alpha \lt 1$ is considered in Magic mass ratios of complete energy-momentum transfer in one-dimensional elastic three-body collisions.

• There is a similar problem in Mathematics SE : Infinite sequence of nested, falling, colliding spheres. The answers to this question about Newton's Cradle are also quite useful. Commented Dec 11, 2016 at 5:48
• Can you prove that the order of collisions is important for a 1D problem? Commented Jul 20, 2017 at 21:30
• @ja72 The results mentioned in the 2nd part of my answer were made using PhET Collision Lab in 1D. They differ depending on order of collision. ... Intuitively, if $m_2 \gg m_1, m_3$ then $m_1$ collides only once with $m_2$. Its final speed depends on the speed of $m_2$, which depends on whether or not it has already collided with $m_3$. Commented Jul 20, 2017 at 22:33
• I learned something new today. If the collision is plastic then the order doesn't matter, otherwise it does. Commented Jul 20, 2017 at 23:30
• Plastic collision means $e=0$. So I am agreeing with you. Commented Jul 21, 2017 at 17:24

I just finished a program for colliding balls in 2D and was looking for info relating to simultaneous collisions (as in a cue ball hitting a rack of pool balls) when I saw this question submitted long ago. Just wanted to offer a short answer using a change of reference frames for the rare case of a simultaneous collision. In any Galilean reference frame with no outside forces acting on the objects within it, any change in momentum and kinetic energy in elastic collisions is conserved. So, let's attach a reference frame to object 2 moving at $$v_{2}$$. To get the motion of objects 1 and 3 in this frame, subtract $$v_{2}$$ from each so $$v_{1, new} = v_{1} - v_{2}$$ and $$v_{3, new} = v_{3} - v_{2}$$ and of course $$v_{2, new} = v_{2} - v_{2} = 0$$ and is stationary in this frame.

Now object 3 is coming at motionless object 2 from the right at $$v_{3, new}$$ and object 1 is coming at object 2 from the left at $$v_{1, new}$$. Using the formula for elastic collisions in 1D we get:

$$v_{1} = \frac{ \left( m_{1} - m_{3} \right) v_{1, new}}{m_{1} + m_{3}} + 2 \frac{m_{3} v_{3, new}}{m_{1} + m_{3}}$$

$$v_{3} = \frac{ \left( m_{3} - m_{1} \right) v_{3, new}}{m_{1} + m_{3}} +2 \frac{m_{1} v_{1, new}}{m_{1} + m_{3}}$$

To transfer back into the so-called 'laboratory frame' where object 2 is moving at $$v_{2}$$ we simply add $$v_{2}$$ to the speed of all objects in object 2's frame to get:

$$v_{1} = \frac{ \left( m_{1} - m_{3} \right) v_{1}}{m_{1} + m_{3}} + 2 \frac{m_{3} v_{3}}{m_{1} + m_{3}} = \frac{6 \cdot 20}{14} + 2 \frac{4 \cdot \left( -25 \right)}{14} = \frac{-80}{14}$$

$$v_{2} = v_{2}$$ (is unchanged) $$= 15$$

$$v_{3} = \frac{ \left( m_{3} - m_{1} \right) v_{3}}{ \left( m_{1} + m_{3} \right)} + 2 \frac{m_{1} v_{1}}{m_{1} + m_{3}} = \frac{ -6 \cdot \left( -25 \right) }{14} + 2 \frac{10 \cdot 20}{14} = \frac{550}{14}$$

$$v_{1} \approx -5.714 ~ \frac{m}{s}$$ (to the left)

$$v_{2} = 15 ~ \frac{m}{s}$$

$$v_{3} \approx 39 ~ \frac{m}{s}$$ (to the right)

However, order of collisions is important and you can use the equations above for the colliding pairs. If object 1 and 2 collide first, you actually get 4 collisions total:

1. After first collision between 1 and 2: v1 = 16.67 m/s and v2 is 21.67 m/s
2. Then 2 collides with 3 and you get: v2 = -19.8 m/s and v3=26.85 m/s
3. Then 2 collides again with 1 and you get: v1 = -7.65 m/s and v2=28.27 so final v1 = -7.65
4. But notice that v2 > v3 so 2 finally catches 3 and you get: v2 = 27 m/s and v3 = 29 m/s for final velocities.

If 2 hits 3 first you get three collisions:

1. After first collision between 2 and 3: v2 = -20.56 m/s, v3 = 19.44 m/s
2. Then 1 collides with 2: v1 = -7.03 m/s, v2 = 33.5 m/s so final v1 = -7.03 m/s
3. Finally, v2 catches up with 3: v2 = 21 m/s, v3 = 35 m/s (final velocities).

In all cases kinetic energy and momentum are conserved.

Total linear momentum is conserved, so however much momentum is added to one body due to a collision must also be removed from another body. These exchanges of momentum are called impulses.

With three bodies there are 2 impulses shown below as $$J_{12}$$ and $$J_{23}$$

To find their magnitudes we need to solve a 2×2 system.

1. First, find the two relative velocities

• $$v_{12} = v_1 - v_2$$
• $$v_{23} = v_2 - v_3$$

with the convention of positive for vectors to the right.

2. During collisions, the relative velocity change is similarly defined by the change in velocities for each body.

• $$\Delta v_{12} = \Delta v_1 - \Delta v_2$$
• $$\Delta v_{23} = \Delta v_2 - \Delta v_3$$
3. Also due to the law of contacts the change in relative velocity is proportional to the relative velocities

• $$v_{12} + \Delta v_{12} =-\epsilon v_{12}$$
• $$v_{23} + \Delta v_{23} =-\epsilon v_{23}$$
4. Then we look at the change in velocity for each body due to the two impulses

• $$\Delta v_1 = - \frac{J_{12}}{m_1}$$
• $$\Delta v_2 = \frac{J_{12} - J_{23}}{m_2}$$
• $$\Delta v_3 = \frac{J_{23}}{m_3}$$
5. Combine all of the above to create the system of equations for the solution. First, use steps 1 and 2 in step 3 to get

• $$\Delta v_1 - \Delta v_2 =-(1+ \epsilon) (v_1 - v_2)$$
• $$\Delta v_2 - \Delta v_3 =-(1+ \epsilon) (v_2 - v_3)$$

Then use step 4 above to expressed the changes in velocities in terms of the impulses

• $$- \frac{J_{12}}{m_1} - \frac{J_{12} - J_{23}}{m_2} =-(1+ \epsilon) (v_1 - v_2)$$
• $$\frac{J_{12} - J_{23}}{m_2} - \frac{J_{23}}{m_3} =-(1+ \epsilon) (v_2 - v_3)$$

The above in matrix form is

$$\begin{vmatrix} \frac{1}{m_1} + \frac{1}{m_2} & - \frac{1}{m_2} \\ - \frac{1}{m_2} & \frac{1}{m_2} + \frac{1}{m_3} \end{vmatrix} \begin{vmatrix} J_{12} \\ J_{23} \end{vmatrix} = (1+\epsilon) \begin{vmatrix} v_1-v_2 \\ v_2-v_3 \end{vmatrix}$$

6. Once the above is solved and the impulses $$J_{12}$$ and $$J_{23}$$ are known, then we use them in step 4 to get the changes in velocities

$$\begin{vmatrix} \Delta v_1 \\ \Delta v_2 \\ \Delta v_3 \end{vmatrix} = \begin{vmatrix} -\frac{1}{m_1} & 0 \\ \frac{1}{m_2} & -\frac{1}{m_2} \\ 0 & \frac{1}{m_3} \end{vmatrix} \begin{vmatrix} J_{12} \\ J_{23} \end{vmatrix}$$

7. The final velocities are found from the initial velocities and the changes in velocity

• $$v_1^\text{after} = v_1 + \Delta v_1$$
• $$v_2^\text{after} = v_2 + \Delta v_2$$
• $$v_3^\text{after} = v_3 + \Delta v_3$$

With some linear algebra, you can find the final velocities as

$$\small \begin{vmatrix} v_{1}^\text{after} \\ v_2^\text{after} \\ v_3^\text{after} \end{vmatrix} = \frac{1}{m_1+m_2+m_3} \begin{vmatrix} m_1 - \epsilon (m_2+m_3) & (1+\epsilon) m_2 & (1+\epsilon) m_3 \\ (1+\epsilon) m_1 & m_2 - \epsilon (m_1+m_3) & (1+\epsilon) m_3 \\ (1+\epsilon) m_1 & (1+\epsilon) m_2 & m_3 - \epsilon ( m_1 + m_2) \end{vmatrix} \begin{vmatrix} v_1 \\ v_2 \\ v_3 \end{vmatrix}$$

Even though the above is mathematically correct, it does not necessarily reflect reality.

Read this very informative article about how to simulate simultaneous collisions and you will find how difficult it is to get realistic results. Some of the key points are

• The order of collisions is important.
• Collisions are more realistic if handled sequentially even for cases where they should be simultaneous, like a Newton's cradle.
• Impulse-based collisions, and force-based collisions each have advantages and disadvantages and neither can handle all cases without problems.
• A hybrid approach seems to yield the best results.