Numerical and experimental investigation of an elastoplastic contact model for spherical discrete elements

A contact model for the normal interaction between elastoplastic spherical discrete elements has been investigated in the present paper. The Walton–Braun model with linear loading and unloading has been revisited. The main objectives of the research have been to validate the applicability of the linear loading and unloading models and estimate the loading and unloading stiffness parameters. The investigation has combined experimental tests and finite element simulations. Both experimental and numerical results have proved that the interaction between the spheres subjected to a contact pressure inducing a plastic deformation can be approximated by a linear relationship in quite a large range of elastoplastic deformation. Similarly, the linear model has been shown to be suitable for the unloading. It has been demonstrated that the Storåkers model provides a good evaluation of the loading stiffness for the elastoplastic contact and the unloading stiffness can be assumed as varying linearly with the deformation of the contacting spheres. The unloading stiffness can be expressed in a convenient way as a function of the Young’s modulus and certain scaling factor dependent on the dimensionless parameter defining the level of the sphere deformation.