Discrete element method investigation on thermally-induced shakedown of granular materials

Temperature change, as a common kind of internal perturbation performed on granular materials, has a significant effect on the bulk properties of granular materials. However, few studies on thermally-induced shakedown under a long-term thermal cycling were reported. In this work, the discrete element method was used to give insight into the thermally-induced shakedown on the fabric and stress states within non-cohesive, frictional granular assemblies. Assemblies were submitted to thermal cycling at a stationary boundary condition after experiencing a one-dimensional compression. Evolution of coordination number, entropy and anisotropy was investigated as well as boundary forces and contact forces. At the same time, effects of the heating rate, the initial vertical load and the magnitude of temperature change were examined. It demonstrates that thermal cycling induces a significant force relaxation within granular materials, while the corresponding granular fabric has a small change. In addition, the entropy and anisotropy decreases with thermal cycling. Moreover, the initial vertical load can constrain the development of thermally-induced fabric change, thereby limiting force relaxation to some degree. Both high heating rate and larger magnitudes of temperature change contribute to more significant force relaxation. However, they cause smaller fabric changes even though they provide larger perturbations.