“…In parallel, the form of Equation (22) makes the deviatoric component f evolve only upon dilation (due to the Macauley brackets), thereby simulating conclusions from micromechanical studies depicting important changes in the sand particle contact normal orientation distribution only upon the dilative phase of monotonic shearing (e.g. Nemat-Nasser and Tobita [44]). Because of the negative sign in front of N in Equation (22), the f component evolves in the opposite sense to tensor [Cn+f], and in this way, the f : n term in the denominator of Equation (20) becomes non-zero, only upon shear reversal following a dilative shearing phase, and this is due to the change in the direction of n. In such a case, the response is characterized by larger plastic strains (due to a decrease of K p ), which under undrained conditions provides the well-known intense contraction following dilation at small p values, which may eventually lead to liquefaction or cyclic mobility.…”