In this paper, a novel semi-energy finite strip method (FSM) is developed based on the concept of first order shear deformation theory (FSDT) in order to attempt the post-buckling solution for thin and relatively thick functionally graded (FG) plates under uniform end-shortening. In order to study the effects of through-the-thickness shear stresses on the post-buckling behavior of FG plates, two previously developed finite strip methods, i.e., semi-energy FSM based on the concept of classical laminated plate theory (CLPT) and a CLPT full-energy FSM, are also implemented. Moreover, the effects of aspect ratio on initial post-buckling stiffness of FG rectangular plates are studied. It has been shown that the variation of the ratio of initial post-buckling stiffness to pre-buckling stiffness ($S^*/S$) with respect to aspects ratios is quite independent of volume fractions of constituents in thin FG plates. It has also been seen that the universal curve representing the variation of ($S^*/S$) with aspect ratio of a FG plate demonstrate a saw shape curve. Moreover, it is revealed that for the thin FG plates in contrast to relatively thick plates, the variations of non-dimensional load versus end-shortening is independent of ceramic-metal volume fraction index. This means that the post-buckling behavior of thin FG plates and the thin pure isotropic plates is similar. The results are discussed in detail and compared with those obtained from finite element method (FEM) of analysis. The study of the results may have a great influence in design of FG plates encountering post-buckling behavior.