Plate-like structures had been thoroughly studied in literature over years to reduce the computational space from 3D to 2D. Many of these theories suffer either from satisfying the free traction condition or thickness extensibility in addition to the consistency of transverse shear strain energy. This work presents a higher order shear deformation thickness-extensible plate theory (eHSDT) for the analysis of plates. The proposed eHSDT satisfies the condition of free traction as other theories do but it also satisfies the condition of consistency of transverse shear strain energy which is neglected by many theories in the area of plates and shells. The implementation of the proposed theory in displacement-based finite element procedure requires continuity of derivatives across elements. This necessary condition was achieved using the penalty enforcement method for derivative-based nodal degrees of freedom across the standard 9-nodes Lagrange element. The theory was tested for elastic bending deformation of Polyether-ether-ketone (PEEK) which is one of the basic materials for medical implants. The theory showed good accuracy compared to experimental data of the three-points bending test. The present eHSDT was also tested for different conditions with a wide range of aspects ratios (thin to thick plates) and different boundary conditions. The accuracy of the proposed eHSDT was verified against exact solutions for these conditions which showed the advantage over other approaches and commercial finite element packages.