37 frictionalStressModel,
38 SchaefferFrictionalStress,
76 *
pow(
Foam::max(alpha - alphaMinFriction, scalar(0)), 10.0);
93 *
pow(
Foam::max(alpha - alphaMinFriction, scalar(0)), 9.0);
106 const scalar I2Dsmall = 1.0e-15;
129 if (alpha[celli] > alphaMax.
value() - 5
e-2)
134 sqrt(1.0/6.0*(
sqr(D[celli].
xx() - D[celli].yy())
135 +
sqr(D[celli].yy() - D[celli].zz())
136 +
sqr(D[celli].zz() - D[celli].
xx()))
137 +
sqr(D[celli].xy()) +
sqr(D[celli].xz())
138 +
sqr(D[celli].yz())) + I2Dsmall