Hertz Contact Model: Complex Loading Paths on a Single Contact

Problem Statement

The Hertz Model provides a nonlinear elastic force-displacement law, with viscous dashpots. The purpose of this example is to analyze the response of a single contact under complex loading paths. Results for a ball-ball contact and a ball-wall interaction are discussed and compared.

List of files (3D):

• cmhertz-setup.dat (3D)
• cmhertzbb.dat (3D)
• cmhertzbw.dat (3D)

PFC3D Models

Two simple systems are set up. The first system is composed of two identical balls of unit diameter; in the second system, the bottom ball is replaced with a circular wall.

All degrees of freedom are fixed for both bodies in contact, and the top ball is displaced in order to load the contact, with the following loading path (in the global coordinate system):

1. normal loading for \(t=0.1~\text{[time-unit]}\) (loading displacement along \(\mathbf{Y}\))
2. rolling for \(t=0.05~\text{[time-unit]}\) (rotation around \(\mathbf{Z}\))
3. rolling for \(t=0.05~\text{[time-unit]}\) (rotation around \(\mathbf{X}\))
4. twisting for \(t=0.05~\text{[time-unit]}\) (rotation around \(\mathbf{Y}\))
5. normal unloading for \(t=0.1~\text{[time-unit]}\) (unloading displacement along \(\mathbf{Y}\))

For both systems, the same loading path is repeated three times, with different set of properties as follows (refer to the Hertz contact model description for a detailed description and listing of the properties):

1. \(G=1.0 \times 10^9\), \(\nu=0.3\) and \(\mu=0.5\);
2. same as 1) with \(M_s=1\) (activate scaling down of the shear force upon normal unload);
3. same as 2) with \(\beta_n=\beta_s=0.5\) and \(M_d=1\) (activate viscous damping)

Model Results

Ball-Ball Contact

Figure 1: Ball-ball contact model at the end of loading stage 1a.

Two identical balls with unit diameter are submitted to the loading paths detailed above. Fig. Fig. #p3d-verif-hzbb-sys shows the state of the system at the end of loading stage 1a (normal loading). In this figure, the contact is represented twice: as a line perpendicular to the contact plane, and as a disk coplanar with the contact plane.

The histories of the magnitude of the normal and shear forces (both in the springs and the dashpots) monitored during loading 1, 2 and 3 are shown in Figs. Figs. #p3d-verif-hzbb1 to #p3d-verif-hzbb3, respectively.

The evolution of the normal force in the spring is identical in all three cases: it increases non-linearly during normal loading stage a), maintains a constant value during rolling and twisting loading stages b), c) and d), then vanishes non-linearly during normal unloading stage e). Its maximal value, which corresponds to an overlap of \(1.0 \times 10^{-3}~\text{[length-unit]}\), agrees with the analytical solution of \(3.0117 \times 10^4~\text{[force-unit]}\) given by equation (3) and (4) in the Hertz contact model description.

Figure 2: Evolution of the spring and dashpot forces (normal and shear) during loading path 1. Ball-ball contact.

Figure 3: Evolution of the spring and dashpot forces (normal and shear) during loading path 2. Ball-ball contact.

Figure 4: Evolution of the spring and dashpot forces (normal and shear) during loading path 3. Ball-ball contact.

The evolution of the shear force in the spring is identical for cases 2) and 3), but differs for case 1), due to the different value of \(M_s\). In case 1), the shear force remains null during normal loading stage a), then increases linearly during rolling loading stage b). During the latter, the slope is given by the initial tangent shear stiffness, which is constant at constant overlap (equation (11) in the Hertz contact model description). During rolling stage c), the magnitude of the shear force increases non-linearly, although the overlap does not change (thus the initial tangent shear stiffness). This is because the loading direction changed, and the shear force is now made up of two components (along the \(\mathbf{X}\) and \(\mathbf{Z}\) global axes). The shear force magnitude is not affected by twisting (loading stage d)), since no shear displacement is built up during this stage. However, the direction of the shear force vector accommodates the relative twist between the balls - as shown in Fig. Fig. #p3d-verif-hzbb1-shear, where the absolute value of the \(x\)-component of the shear force decreases linearly, while the absolute value of the \(z\)-component of the shear force increases linearly during the twisting loading stage (mechanical time from \(0.2\) to \(0.25~\text{[time-unit]}\)). Finally, since no shear displacement is built up during the normal unloading stage 1e), the shear force remains constant at the beginning of this stage (Fig. Fig. #p3d-verif-hzbb1 - mechanical time from \(0.25\) to \(\approx 0.28~\text{[time-unit]}\)). However, as the unloading proceeds the Coulomb slip criterion is eventually met, and the shear force vanishes with the normal force (mechanical time from \(\approx 0.28\) to \(0.35~\text{[time-unit]}\)).

Figure 5: Evolution of the spring shear force components during loading path 1. Ball-ball contact.

The evolution of the shear force in cases 2) and 3) is similar to that of case 1), except during the normal unloading stage e). For those cases, scaling-down of the shear force upon normal unload is activated (\(M_s=1\)). Thus the shear force starts decreasing linearly while the normal force is decreasing, until the Coulomb criterion is met and the shear force vanishes proportionally to the normal force (Fig. Fig. #p3d-verif-hzbb1 - mechanical time from \(0.25\) to \(0.35~\text{[time-unit]}\)). Because the shear force magnitude at the beginning of slip is therefore lower in cases 2) and 3) than in case 1), the energy dissipated during sliding will also be lower. This can be seen by comparing Figs. Figs. #p3d-verif-hzbb1-energy, #p3d-verif-hzbb2-energy and #p3d-verif-hzbb3-energy, which show the evolution of the contact energies (strain energy and slip and dashpot work) during loading cases 1), 2) and 3), respectively. Note that viscous damping is activated only for loading case 3), which is the only case where the dashpot forces (Fig. Fig. #p3d-verif-hzbb3) and dashpot work (Fig. Fig. #p3d-verif-hzbb1-energy) are nonzero.

Figure 6: Evolution of the contact energies (strain energy, slip work and dashpot work) during loading path 1. Ball-ball contact.

Figure 7: Evolution of the contact energies (strain energy, slip work and dashpot work) during loading path 2. Ball-ball contact.

Figure 8: Evolution of the contact energies (strain energy, slip work and dashpot work) during loading path 3. Ball-ball contact.

Ball-Wall Interaction

Figure 9: Ball-wall interaction model at the end of loading stage 1a.

Fig. Fig. #p3d-verif-hzbw-sys shows the state of the system at the end of loading stage 1a (normal loading), for the ball-wall interaction case. In this case, multiple contacts exist between the ball and the facets that make up the wall. However since full contact resolution mode is active (default setting), only one contact is deemed active at a time, although more than one contact may exhibit a positive overlap (refer to the description of the faceted wall logic for further details).

The histories of the magnitude of the normal and shear forces (both in the springs and the dashpots) monitored during loading 1, 2 and 3 are shown in Figs. Figs. #p3d-verif-hzbw1 to #p3d-verif-hzbw3 respectively, and the histories of the contact energies are shown in Figs. Figs. #p3d-verif-hzbw1-energy to #p3d-verif-hzbw3-energy.

All quantities evolve qualitatively similarly to the ball-ball contact case. However they are quantitatively different, although the contact shear modulus \(G\) and Poisson ratio \(\nu\) are identical, since the contact effective radius is different (equation (6) in the Hertz contact model description).

Figure 10: Evolution of the spring and dashpot forces (normal and shear) during loading path 1. Ball-facet contact.

Figure 11: Evolution of the spring and dashpot forces (normal and shear) during loading path 2. Ball-facet contact.

Figure 12: Evolution of the spring and dashpot forces (normal and shear) during loading path 3. Ball-facet contact.

Figure 13: Evolution of the contact energies (strain energy, slip work and dashpot work) during loading path 1. Ball-facet contact.

Figure 14: Evolution of the contact energies (strain energy, slip work and dashpot work) during loading path 2. Ball-facet contact.

Figure 15: Evolution of the contact energies (strain energy, slip work and dashpot work) during loading path 3. Ball-facet contact.

Conclusion

In this example, a ball-ball contact and a ball-wall interaction, with the Hertz contact model, are exercised under complex loading paths. The evolution of the forces and energies are discussed and compared.