Contact model for elastically anisotropic bodies and efficient implementation into the discrete element method

We introduce a contact law for the normal force generated between two contacting, elastically anisotropic bodies of arbitrary geometry. The only requirement is that their surfaces be smooth and frictionless. This obtained from simplification exact solution to continuum elasticity problem takes familiar form Hertz’ law, with difference being orientation-dependence material-specific modulus. remarkably accurate when compared solution, wide range materials surface geometries. describe computationally efficient implementation into discrete element method code, taking advantage precomputation modulus over all possible orientations. Finally, we showcase application examples based on real where elastic anisotropy particles induces noticeable effects macroscopic behavior. Notably, second example demonstrates ability engineer tunable vibrational band gaps in one-dimensional granular crystal by mere rotation constituent spherical particles.