Mirror boundary conditions were imposed on each edge of the mesh, with stress and constraining equal displacements imposed on edges BC and CD, and with edges AB and DA constrained in the orthogonal direction.