In order to investigate the stress concentration phenomenon and analyze the distribution characteristics of the interlaminar stress in the hole edge region of composite laminates. Based on the generalized mixed variational principle, the generalized mixed finite element model for laminated plates with various stacking modes were established. The stress field variables were divided into the interlaminar stress and the in-plane stress, with the introduction of stress boundary conditions to ensure the physical continuity of interlaminar stresses between layers and the discontinuity of in-plane stresses between layers.The interlaminar stresses at the edge of the laminated plate hole were respectively analyzed through the thickness direction and the circumferential direction. Numerical examples demonstrated that the incompatible generalized mixed element could obtain more accurate stress singularity results than the 8-node three-dimensional solid incompatible displacement element results solved by the finite element software Abaqus. Stresses on both upper and lower surfaces of the laminated plate consistently reflected actual situations. The research indicates that compared with the displacement element, the incompatible generalized mixed element can more effectively capture the high stress gradient of the interlaminar stresses at the edge of the laminated plate hole, which provides a new idea for the optimal design of the laminate.