Abstract: Generalized mixed finite element method of pure elastomer was introduced into the static analysis of piezoelectric materials. Because the 8-node hexahedral noncompatible solid element was used to solve the whole structure discretely, many artificial assumptions in the theory of plate and shell were abandoned. The addition of noncompatible terms makes this method exhibit better numerical performance than the same kind of compatible elements. In this method, the stress boundary conditions and displacement boundary conditions were considered simultaneously, and the interlaminar stress and in-plane stress were treated separately in the solution process, and the in-plane stress was calculated according to the constitutive relation of each layer, so that the obtained interlaminar stress and in-plane stress are closer to the exact solution. The accuracy of the proposed method was illustrated by several representative numerical examples of laminates. Compared with the traditional analytical and numerical methods, the theoretical method presented in this paper has advantages in applicability and effectiveness.