RD-E: 0901 - Billiards (Pool)

A pool game is modeled to show the transmission of momentum between one impacting ball and 15 impacted balls.

The purpose of this example is to investigate the transmission of momentum between several balls. Contact with the various interfaces using the Penalty and Lagrange Multipliers' method is analyzed.

Options and Keywords Used

  • 16-node thick shell and sphere mesh (/PROP/TYPE20 (TSHELL))
  • TYPE7 interface using the Lagrange Multipliers method and the Penalty method (/INTER/TYPE7)
  • TYPE16 sliding and tied interface, and quadratic surface contact (/INTER/LAGMUL/TYPE16)
  • Elastic shock
  • Momentum transmission and shock wave
  • Initial velocities (/INIVEL)
An initial velocity of 1.5 ms-1 in X direction is applied to all nodes of the white (cue) ball.

rad_ex_fig_9-7
Figure 1. Initial Translational Velocities of the Impacting Ball

All nodes of the lower face of the table are completely fixed (translations and rotations).

Gravity is considered for all the balls nodes. A function defines the gravity acceleration in the Z direction compared with time. Gravity is activated using /GRAV.

ex9_fig8
Figure 2. Gravity Function (-0.00981 mm.ms-2 )

Input Files

The input files used in this example include:
Billiard_game/Interface_16
<install_directory>/hwsolvers/demos/radioss/example/09_Billiards/Billiards_model/BILLARD*
<install_directory>/hwsolvers/demos/radioss/example/09_Billiards/Billiards_model/Supplement_Interface7Lag/BILLARD*

Model Description

Pool is a game consisting of 16 balls, each 50.8 mm in diameter. It is played on a small billiard table measuring 1800 mm x 900 mm. Fifteen (15) balls are placed in a triangle to enable their tight grouping. The initial velocity of the shooting (cue) ball is presumed equal to 1.5 ms-1. Elastic rebounds are observed.

rad_ex_fig_9-1
Figure 3. Pool Game

Units: mm, g, N, MPa.

The material is subjected to a linear elastic law (/MAT/LAW1) with the following properties:
Material Properties: Balls (pheonolic resin)
Initial density
0.00137 g.mm-3
Young's modulus
10500 [ MPa ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWadaqaai Gac2eacaGGqbGaaiyyaaGaay5waiaaw2faaaaa@3BE6@
Poisson ratio
0.3
Material Properties: Frame (polymer)
Initial density
0.001 g.mm-3
Young's modulus
1000 [ MPa ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWadaqaai Gac2eacaGGqbGaaiyyaaGaay5waiaaw2faaaaa@3BE6@
Poisson ratio
0.49
Material Properties: Plate (slate)
Initial density
0.0028 g.mm-3
Young's modulus
62000 [ MPa ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWadaqaai Gac2eacaGGqbGaaiyyaaGaay5waiaaw2faaaaa@3BE6@
Poisson ratio
0

Model Method

The balls are meshed with 16–node solid shells (quadratic elements) in order to improve the conditions of contact by taking into account the curvatures. The frame of the table is made of 16-node solid shells to comply with the interface used. The plate is modeled using only one solid element. The 16-node thick shells are considered as solid elements. They are defined by a thick TYPE20 shell property (number 16 solid formulation for quadratic 16-node thick shells, fully-integrated with 2x2x2 integration points).

rad_ex_fig_9-2
Figure 4. Pool Game Mesh

rad_ex_fig_9-3
Figure 5. Mesh for Balls

rad_ex_fig_9-4
Figure 6. 16-node Thick Shell Element
The TYPE16 interface with the Lagrange Multipliers method is used to model the ball/ball and balls/table contacts. An interface must be defined for each ball (that is: 16 interfaces in total). An additional interface is used to define the contacts between the balls and the table (plate and frame).

rad_ex_fig_9-5
Figure 7. TYPE16 Interface: Secondary SHEL16 for Balls and Main SHEL16 for the Table

rad_ex_fig_9-6
Figure 8. Example of the TYPE16 Interface Defined for the Contact Between Balls

Secondary nodes (red) are extracted from the external surfaces of the parts.

Results

Curves and Animations

Due to the faceting of the ball, contact between the impacting ball and the impacted balls is not perfectly symmetrical and momentum is not homogeneously transmitted among the balls. An apparent physical strike thus results.

rad_ex_fig_9-8B
Figure 9.

rad_ex_fig_9-9
Figure 10. Collision of the Balls
rad_ex_fig_9-10A rad_ex_fig_9-10B
rad_ex_fig_9-10C rad_ex_fig_9-10D
rad_ex_fig_9-10E rad_ex_fig_9-10F
Figure 11. History of the Balls' Motions (contact control: TYPE16 interface)