higher-order theory;  P-type finite element;  laminates;  Dimensionality reduction method

For thicker laminates, the stress calculation results based on general higher-order theory often cannot meet the engineering requirements. The stress distribution law of laminate was analyzed by the P-type finite element method based on high-order theory. The displacement field is interpolated by the integral of the Legendre polynomial, the three-dimensional laminate problem is reduced into a two-dimensional problem by dimension reduction method, and the internal degrees of freedom of the element are condensed by static condensation, which can effectively calculate the stress distribution of the anisotropic material laminate. The results show that P-type finite elements are used to solve the problem of quadrilateral composite laminates laid in any thickness and direction.