In this work, it is studied on surface wrinkles on buckling surfaces of a sandwich PZT/Metal/PZT rectangular thick plate with edge cracks within the scope of the 3D linearized theory of stability loss in the framework of the piecewise homogeneous body model. It is supposed that between the face and core layers of the sandwich rectangular plate there is an interface edge crack and around these edge cracks an ideal contact conditions are satisfied. Also assume that this plate is subjected to bi-axial uniformly-distributed compressive forces acting on two edge surfaces only and neither mechanical nor electrical load are act on the face surfaces and, also cracks surfaces of the plate. In addition, three of the lateral surfaces of the plate is simply supported mechanically and all of the lateral surfaces of the plate grounded for the PZT layers' surface only. In the analyzing procedure, before the plate is loaded (i.e. in the natural state), it is assume that surfaces of the considered edge cracks have insignificant initial imperfections. Due to action of the aforementioned compressive forces the evolution of these initial imperfections is investigated and as a result of this investigation the values of the critical buckling forces for the considered sandwich plate are found from the initial imperfections criteria. Under buckling delamination of the plate some problem parameters can lead to different surface wrinkles unlike initial imperfection of the cracks’ face surfaces. Our aim is to determine the effect of geometrical and material parameters on the surface wrinkles for the buckling delamination of PZT sandwich rectangular plate by employing the 3D finite element method (3D-FEM).