RD-E: 0902 Collision between Two Balls

Two balls are now considered in order to study the behavior of impacting spherical balls.

Input Files

The input files used in this example include:

Collision study: <install_directory>/hwsolvers/demos/radioss/example/09_Billiards/Collision_simulation/COLLISION*

Model Description

Study on Trajectories

The balls' behavior is described using the parameters (angles and velocities) shown in Figure 1. The numerical results are compared with the analytical solution, assuming a perfect elastic rebound (coefficient of restitution is equal to 1).

Initial Values
V₁: 0.7m.s^-1
V₂: 1m.s^-1
$θ$ ₁: 40°
$θ$ ₂: 30
mass_ball: 44.514g

Model Method

The balls and the table have the same properties, previously defined for a pool game. The dimensions of the table are 900 mm x 450 mm x 25 mm and the balls' diameter is 50.8 mm. The balls and the table are meshed with 16-node thick shell elements for using the TYPE16 Lagrangian interface.

The initial translational velocities are applied to the balls in the /INIV Engine option. Velocities are projected on the X and Y axes.

Gravity is considered for the balls (0.00981 mm.ms^-2 ).

The ball/ball and balls/table contact is modeled using the TYPE16 interface (secondary nodes/main 16-node thick shells contact). The interface defining the ball/ball contact is shown in Figure 4.

Analytical Solution

Take two balls, 1 and 2 from masses

m_{1}

and

m_{2}

, moving in the same plane and approaching each other on a collision course using velocities

V_{1}

and

V_{2}

, as shown in Figure 5.

Velocities are projected onto the local axes

n

and

t

. To obtain the velocities and their direction after impact, the momentum conservation law is recorded for the two balls:(1)

- m_{1} V_{1 n} + m_{2} V_{2 n} = m_{1} V_{1 n}^{'} - m_{2} V_{2 n}^{'}

Or(2)

- m_{1} V_{1} \sin θ_{1} + m_{2} V_{2} \sin θ_{2} = m_{1} V_{1}^{'} \sin θ_{1}^{'} - m_{2} V_{2}^{'} \sin θ_{2}^{'}

The shock is presumed elastic and without friction. Maintaining the translational kinetic energy is respected as there is no rotational energy:(3)

\frac{1}{2} m_{1} (V_{1 n}^{' 2} + V_{1 t}^{' 2}) + \frac{1}{2} m_{2} (V_{2 n}^{' 2} + V_{2 t}^{' 2}) = \frac{1}{2} m_{1} (V_{1 n}^{2} + V_{1 t}^{2}) + \frac{1}{2} m_{2} (V_{2 n}^{2} + V_{2 t}^{2})

Such equality implies that the recovering capacity of the two balls corresponds to their tendency to deform.

This condition equals one of the elastic impacts, with no energy loss. Maintaining the system's energy gives:(4)

(V_{2 n}^{'} - V_{1 n}^{'}) = - (V_{2 n} - V_{1 n})

This relation means that the normal component of the relative velocity changes into its opposite during the elastic shock (coefficient of restitution value e is equal to the unit).

The following equations must be checked for normal components:(5)

\begin{array}{l} V_{2 n}^{'} = V_{1 n}^{'} = - (V_{2 n} - V_{1 n}) \\ m_{2} V_{2 n}^{'} + m_{1} V_{1 n}^{'} = m_{2} V_{2 n} + m_{1} V_{1 n} \end{array}

The equations system using V'₁ and V'₂ as unknowns is easily solved:(6)

\begin{array}{l} V_{2 n}^{'} = (\frac{m_{2} - m_{1}}{m_{2} + m_{1}}) V_{2 n} + (\frac{2 m_{1}}{m_{1} + m_{2}}) V_{1 n} \\ V_{1 n}^{'} = (\frac{m_{1} - m_{2}}{m_{1} + m_{2}}) V_{1 n} + (\frac{2 m_{2}}{m_{1} + m_{2}}) V_{2 n} \end{array}

It should be noted that these relations depend upon the masses ratio.

As the balls do not suffer from velocity change in the t-direction, maintaining the tangential component of each sphere's velocity provides:(7)

\begin{array}{l} V_{1 t}^{'} = V_{1 t} = V_{1} \cos θ_{1} \\ V_{2 t}^{'} = V_{2 t} = V_{2} \cos θ_{2} \end{array}

The norms of velocities after shock result from the following relations.(8)

V_{1}^{'} = {({(V_{1 n}^{'})}^{2} + {(V_{1 t}^{'})}^{2})}^{1 / 2}

(9)

V_{2}^{'} = {({(V_{2 n}^{'})}^{2} + {(V_{2 t}^{'})}^{2})}^{1 / 2}

In this example, balls have the same mass: m₁ = m₂.

Therefore, $V_{2}^{'} = V_{1 n}$ and $V_{1 n}^{'} = V_{2 n}$ .

The norms of the velocities are given using the following relations, depending on the initial velocities and angles. Used to determine the analytical solutions (angles and velocities after collision):(10)

V_{1}^{'} = {({(V_{2})}^{2} \sin^{2} θ_{2} + (V_{1}) \cos^{2} θ_{1})}^{1 / 2}

(11)

V_{2}^{'} = {({(V_{1})}^{2} \sin^{2} θ_{1} + (V_{2}) \cos^{2} θ_{2})}^{1 / 2}

By recording the projection of the velocities, directions after shock can be evaluated using relation. Used to determine the analytical solutions (angles and velocities after collision):(12)

θ_{1}^{'} = \arcsin (\frac{V_{2}}{V_{1}} \sin θ_{2})

(13)

θ_{2}^{'} = \arcsin (\frac{V_{1}}{V_{2}} \sin θ_{1})

Results

Numerical Results Comparison with the Analytical Solution

Figure 6 shows the trajectories of the balls' center point obtained using numerical simulation before and after collision.

For given initial values of

V_{1}

V_{2}

θ

₁ and

θ

₂, simulation results are reported in Table 1.

Table 1. Comparison of Results for After Collision
	Numerical Results	Analytical Solution
$θ$ _1'	42.27°	44.72°
$θ$ _2'	26.75°	26.48°
V _1'	0.731 m/s	.731 m/s
V _2'	0.969 m/s	0.977 m/s

Conclusion

The simulation corroborates with the analytical solution. The 16-node thick shells are fully-integrated elements without hourglass energy. This modeling provides a good transmission of momentum. However, the TYPE16 interface does not take into account the quadratic surface on the secondary side (ball 2), due to the node to thick shell contact. Accurate results are obtained for a collision without penetrating the quadratic surface of the secondary side in order to confirm impact between the spherical bodies.

A fine mesh could improve the results.