We propose a particle simulation method derived from a Hamiltonian with a mechanism of energy dissipation through particle-particle contact friction. The total energy dissipation variable S is paired with a function \(F(\mu )\) which works at each contact among constituent particles. A two dimensional quasi-static system is simulated to make clear the possibilities and the limitations of the proposed method. By changing the macroscopic dissipation coefficients \(\gamma \) systematically, we have shown how the macroscopic energy dissipation \(\gamma S\) and the microscopic contact dissipation \(F(\mu )\) is connected through the state of the granular system. In the compression/decompression cycle, pressure shows hysterisis but have an opposite trend depending whether it is in the direction of compression or not.
Graphical abstract
