“…(39), P ni ,P τ i ,θ i , dλ j (i = 1, 2, ..NC, j = 1, 2, ..NumGausP) are taken as unknowns, and the yield function f j and the incremental displacement du j , which should satisfy the equilibrium Eqs. (12) or (13), are the functions of these unknowns. The dimension of H is NumGaus P + 3 × N C. Note that Eqs.…”