Stephanie Deboeuf,Olivier Dauchot,Lydie Staron,Anne Mangeney,Jean-Pierre Vilotte
Abstract:Discrete numerical simulations are performed to study the evolution of the micro-structure and the response of a granular packing during successive loading-unloading cycles, consisting of quasi-static rotations in the gravity field between opposite inclination angles. We show that internal variables, e.g., stress and fabric of the pile, exhibit hysteresis during these cycles due to the exploration of different metastable configurations. Interestingly, the hysteretic behaviour of the pile strongly depends on the maximal inclination of the cycles, giving evidence of the irreversible modifications of the pile state occurring close to the unjamming transition. More specifically, we show that for cycles with maximal inclination larger than the repose angle, the weak contact network carries the memory of the unjamming transition. These results demonstrate the relevance of a two-phases description -strong and weak contact networks- for a granular system, as soon as it has approached the unjamming transition.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: **To study the evolution of microstructure and response during the cyclic tilting of particle packing, especially the memory effect and hysteresis exhibited by the system near the unjamming transition**.
Specifically, the author studied the behavior of particle packing during quasi - static rotation in a gravitational field through discrete numerical simulations. These rotations include loading - unloading cycles, that is, the quasi - static rotation of the particle pile between opposite tilt angles. The focus of the study is:
1. **The hysteresis behavior of internal variables (such as stress and structural texture) during the cycle**.
2. **The influence of different maximum tilt angles on the hysteresis behavior**, especially when the maximum tilt angle exceeds the angle of repose, how the weak contact network carries the memory of the unjamming transition.
3. **The contribution of the strong contact network and the weak contact network to the overall response**, verifying the relevance of the two - phase description (strong and weak contact networks) for the unjamming transition.
### Main problem analysis
- **Memory effect of unjamming transition**: Research shows that when the maximum tilt angle of the cycle exceeds the angle of repose, the weak contact network can remember the process of the unjamming transition. This indicates that near the unjamming transition, the microstructure of the system has undergone irreversible changes.
- **Hysteresis phenomenon**: Internal variables (such as stress and structural texture) show obvious hysteresis loops during the loading and unloading processes, and this hysteresis phenomenon depends on the maximum tilt angle of the cycle. In particular, when the maximum tilt angle is close to or exceeds the angle of repose, the hysteresis loop becomes more significant.
- **Validity of the two - phase description**: The study found that the weak contact network plays a key role near the unjamming transition, verifying the importance of the two - phase description (i.e., strong and weak contact networks) for understanding the behavior of particle systems.
### Conclusion
This study reveals the complex behavior of particle systems near the unjamming transition, especially the memory effect and hysteresis phenomenon exhibited during the cyclic loading - unloading process. These findings are of great significance for understanding the behavior of granular materials in geological and industrial processes and provide new insights into the unjamming transition in condensed matter physics.
### Formula representation
Some formulas involved in the paper are as follows:
- **Inertial Number**:
\[
I_{\text{max}}=\frac{\gamma_{\text{max}}d}{\sqrt{\rho / P}}
\]
where $\gamma_{\text{max}}$ is the maximum shear strain rate, $d$ is the particle diameter, $\rho$ is the particle volume density, and $P$ is the typical pressure caused by gravity.
- **Stress tensor calculation**:
\[
\sigma_{ij}=\frac{1}{V}\sum_{\alpha = 1}^{n_c}f_{\alpha}l_il_j
\]
where $i$ and $j$ represent spatial dimensions, $\alpha$ represents the contact point, $\vec{f}$ is the contact force, $\vec{l}$ is the vector in the normal direction of the contact surface, and $n_c$ is the total number of contacts in the representative volume $V$.
- **Normal and tangential stress**:
\[
\sigma_N=\|(\sigma\vec{n}_S)\vec{n}_S\|
\]
\[
\sigma_T=\|\sigma\vec{n}_S-\sigma_N\vec{n}_S\|
\]
where $\vec{n}_S$ is the unit normal vector of the free surface.
These formulas are used to describe the stress state and deformation behavior of particle systems during the cyclic loading - unloading process.