# Advanced Classical Mechanics/Compound ballistic pendulum with spinning ball

The easiest way to create motion with minimum friction is to roll balls down a ramp. For that reason it is worthwhile to develop theories for such systems. Here we consider the ballistic pendulum. This calculation models the stable configuration with a completely inelastic collision. However, it is not difficult to generalized these results to include partially elastic collisions, or even the knockdown of a block of wood by a ball rolled down a ramp.

The intent of this essay is for advanced physics majors to understand the derivation, and for beginners to utilize the results of these calculations. It should not be difficult to develop **Matlab** or **Python** codes that beginners could use in order to correlate classroom experiments to theory.

**Warning**: It must emphasized that the assumptions made regarding an "inelastic" collision made here might not always hold. See **below** for further discussion.

## Unit vectors and multiple coordinate systems[edit | edit source]

Vector notation becomes more powerful when two or more coordinate systems are used, especially when the orthogonality of the unit vectors permit simple transformations between the coordinate systems. We illustrate this with the impulse delivered to the block by the ball. The impulse *to the block* is proportional to the final minus initial velocity of the ball:

Using relations between the unit vectors easily derived by inspection of the figure, we can derive the **Rotation matrix** between the coordinate systems. Defining as the angle between and , we have:

Take the inner product of with the unit vectors and to obtain:

- and

## Four equations of motion.[edit | edit source]

### The impulse torque on the block[edit | edit source]

Following the concepts developed at Wikipedia's * Impulse* we integrate the equations of motion under the assumption that the collision is so brief that the ball and block change velocity but not position during the collision. The block is assumed to be initially at rest, with a vertical wall facing the incoming ball. The first equation relates the impulse of torque to the angular velocity of the block about the axis, which is situated a distance

*R*away from the point of impact:

where *I* is the block's moment of inertia and the minus sign above is due to the "negative" orientation of the unit vector.

### The impulse torque on the ball[edit | edit source]

For a ball that is a solid sphere, we use Wikipedia's list of moments of inertia relate the change in the ball's angular velocity to the appropriate component of the impulse:

where here the minus sign occurs because the impulse to the ball is opposite that to the block.

### The (questionable) assumption of brief stationary contact between ball and block[edit | edit source]

The third and fourth equation arise from a single (two-dimensional) vector equation that arises from the fact that for an instant in time, the ball and block are in contact. At that moment, the portion of the ball and block that are in contact are moving at the same velocity. It must be emphasized that this equality of both objects at the point of contact might not actually occur.^{[1]} We this caveat in mind, we use the concept of **relative velocity** and use the rotation of the ball's surface about the center of mass to find the velocity at the point of contact with the block.

### On nature of these equations[edit | edit source]

Before solving a set equations, it is important to first classify the type of equations they are. In this case, the four equations are *linear*, or *first order*, because in the four unknowns, no polynomial has **degree** higher than one. Here, the "unknowns" are the variables that describe the angular and linear velocities associated with the ball and ramp immediately *after* the collision, namely . These "unknowns" can be deduced if the initial velocities are known. This *linear* set of equations is also **inhomogeneous** , as discussed in Wikipedia's **System of linear equations,**^{[2]} where it is shown that such equations have either no solution, one unique solution, or an infinite number of solutions. We shall later see that one unique solution exists for this system.

## Finding a solution[edit | edit source]

It would be interesting to see if Matlab or Mathematica can solve this problem symbolically. Meanwhile, we shall outline the steps required to obtain a solution that arrives at the final velocities (angular and linear) given the initial velocities. The first step is to write these equations as an inhomogeneous set of equations:

## Algebra[edit | edit source]

The algebra is tedious. In a physics lab we will break down the steps and have individual groups verify two or three steps. We will also see if Matlab symbolic manipulation is up to the task.

{{cot|click to view the steps}}

Step 0:Collect the equations into one place:

Step 1: Eliminate *J* (after factoring)

If everything above this **Algebra** section is true. Then step 0 and the first of Step 1 is confirmed:^{[3]}

Step 2: Move inhomogeneous terms to the RHS (right hand side) and simplify the vector equation.

Step 3: When working by hand it helps to simplify the notation with temporary variables:

Step 4: Substitute

Step 5: Substitute:

Step 6: Reduce to three equations using

Step 7: Substitute :

Step 8: Simplify

Step 9: Substitute

Step 10: Substitute and

{{cob}}

Step 11: Make the small *S* approximation for a thin board

## References and footnotes[edit | edit source]

- ↑ For example, if friction is responsible for bringing the ball and block into contact, the ball might bounce off the block before this friction brings the ball to rest with respect to the block. This situation does not occur for one dimensional collisions because the two objects remain in contact after the inelastic collision.
- ↑ The Wikipedia article expresses these linear equations in terms of the matrix equation, , where
*A*is a matrix that depends on the mass and dimensions of the ball and block. The four final velocities correspond to the unknown column vector,**x**. The equations are inhomogeneous because the column vector**b ≠ 0**. - ↑ phy2400 Feb 26 2015