“…The mixed finite element proposed by Matsushima et al, 76 referred to here as Q8F4L1, has been tested and successfully used in the literature; 1,4,52,54 thus it is employed to solve bi-dimensional problems in this contribution. It has eight nodes for displacements and geometry interpolation, with quadratic shape functions, a bi-linear interpolation for the relaxed deformation gradient, with four nodes, and a constant Lagrange multiplier field in the element domain, resulting in a total of 36 DOFs (Figure 4A).…”