Yang Ming, Taiebat Mahdi, Mutabaruka Patrick, Radjaï Farhang
Department of Civil Engineering, University of British Columbia, Vancouver BC V6T1L7, Canada.
CNRS, University of Montpellier, LMGC, 34090 Montpellier, France.
Phys Rev E. 2021 Mar;103(3-1):032904. doi: 10.1103/PhysRevE.103.032904.
By means of 3D particle dynamics simulations, we analyze the microstructure of granular materials subjected to isochoric (constant volume) cyclic shearing, which drives the system towards a liquefaction state characterized by loops of jamming-unjamming transition with periodic loss of strength and irreversible accumulation of shear strain. We first show that the macroscopic response obtained by these simulations agrees well with the most salient features of the well-known cyclic behavior of granular materials both before and after liquefaction. Then we investigate the evolution of particle connectivity, force transmission, and anisotropies of contact and force networks. The onset of liquefaction is marked by partial collapse of the force-bearing network with rapid drop of the coordination number and nonrattler fraction of particles, and significant broadening of the contact force probability density function, which begins in the preliquefaction period. We find that the jamming transition in each cycle occurs for a critical value of the coordination number that can be interpreted as the percolation threshold of the contact network and appears to be independent of the initial mean stress, void ratio, and cyclic shear amplitude. We show that upon unjamming in each cycle an isotropic loss of contacts occurs and is followed by the development of high contact anisotropy and a large proportion of particles with only two or three contacts. The higher mobility of the particles also involves a lower degree of frustration of particle rotations and thus lower friction mobilization and tangential force anisotropy. These findings are relevant to both undrained cyclic deformations of saturated soils and rheology of dense non-Brownian suspensions where volume change is coupled with pore liquid drainage conditions.
通过三维颗粒动力学模拟,我们分析了等容(恒定体积)循环剪切作用下颗粒材料的微观结构,这种剪切作用会使系统趋向于液化状态,其特征是出现阻塞 - 非阻塞转变的循环,伴随着强度的周期性损失和剪切应变的不可逆积累。我们首先表明,这些模拟得到的宏观响应与颗粒材料在液化前后众所周知的循环行为的最显著特征吻合良好。然后我们研究了颗粒连通性、力传递以及接触和力网络各向异性的演变。液化的开始以承载网络的部分坍塌为标志,伴随着配位数和非游动颗粒分数的快速下降,以及接触力概率密度函数的显著展宽,这一过程在液化前期就已开始。我们发现每个循环中的阻塞转变发生在配位数的一个临界值处,该临界值可解释为接触网络的渗流阈值,并且似乎与初始平均应力、孔隙率和循环剪切幅度无关。我们表明,在每个循环的非阻塞阶段,会出现各向同性的接触损失,随后会发展出高接触各向异性以及很大比例仅具有两三个接触的颗粒。颗粒较高的迁移率还涉及颗粒旋转的较低受阻程度,从而导致较低的摩擦 mobilization 和切向力各向异性。这些发现与饱和土的不排水循环变形以及体积变化与孔隙液体排水条件耦合的致密非布朗悬浮液的流变学都相关。 (注:原文中“lower friction mobilization”里“mobilization”在这里含义不太明确,暂保留英文,根据上下文可能是“较低的摩擦发挥程度”之类的意思)