Couple stress-based unsymmetric 8-node planar membrane elements with good tolerances to mesh distortion

Purpose The paper aims to propose two new 8-node quadrilateral membrane elements with good distortion tolerance for the modified couple stress elasticity based on the unsymmetric finite element method (FEM). Design/methodology/approach The nodal rotation degrees of freedom (DOFs) are introduced into the virtual work principle and constrained by the penalty function for approximating the test functions of the physical rotation and curvature. Therefore, only the C0 continuity instead of C1 continuity is required for the displacement during the element construction. The first unsymmetric element assumes the test functions of the displacement and strain using the standard 8-node isoparametric interpolations, while these test functions in the second model are further enhanced by the nodal rotation DOFs. Besides, the trial functions in these two elements are constructed based on the stress functions that can a priori satisfy related governing equations. Findings The benchmark tests show that both the two elements can efficiently simulate the size-dependent plane problems, exhibiting good numerical accuracies and high distortion tolerances. In particular, they can still exactly reproduce the constant couple stress state when the element shape deteriorates severely into the degenerated triangle. Moreover, it can also be observed that the second element model, in which the linked interpolation technique is used, has better performance than the first one, especially in capturing the steep gradients of the physical rotations. Originality/value As the proposed new elements use only three DOFs per node, they can be readily incorporated into the existing finite element (FE) programs. Thus, they are of great benefit to analysis of size-dependent membrane behaviors of micro/nano structures.