The accurate evaluation of interlaminar stresses is of great significance in the analysis and design of laminated and piezoelectric laminated structures because complex behaviours of these stresses near free edges initiate edge delamination that raises concerns about the structural integrity and reliability. This thesis presented 3D hybrid analyses on the interlaminar stresses to investigate the electromechanical coupling and free edge effects of piezoelectric laminated plates with an emphasis on the realistic distributions of the 3D stress and electric fields near free edges. In this research, the state space equations for simply-supported and free-edge piezoelectric laminates under transverse loads and infinite long free-edge piezoelectric laminates under uniaxial extension were obtained in the framework of 3D piezoelasticity by considering all the independent elastic and piezoelectric constants. The equations satisfy the traction-free and open-circuit boundary conditions at free edges and the continuity conditions across all interfaces. On the basis of the transfer matrix and recursive solution approaches, 3D exact solutions were sought by a novel non-uniform layer refinement technique to evaluate the accuracy of the finite element method (FEM), and realistic gradients of interlaminar stresses and electric fields were captured. The FEM results were in good agreement with those from the present solutions except for the regions near free edges. For simply-supported and free-edge laminates, stress variations with material properties, geometries and stacking sequences were obtained. The interlaminar stress τxz was dominant at corners and τyz also tended to contribute to delamination. In the infinite long free-edge laminates, sigmaz, τyz, Ey and Ez exhibited significant gradients near free edges. Furthermore, the considerable influence of the electromechanical coupling effect on interlaminar stresses revealed that piezoelectric laminates were more susceptible to edge delamination and the application of closed-circuited surface conditions might prevent such edge delamination. The present analytical solution demonstrated an improvement in precision over other 2D analytical and numerical solutions and could serve as a benchmark for the determination of interlaminar stresses and electric fields near the free edges of the piezoelectric laminates.