This paper deals with second order statistics of post buckling load of shear deformable laminated composite plates resting on linear elastic foundation with random system properties. The formulation is based on higher order shear deformation plate theory in general von Karman sense, which includes foundation effect using two-parameter Pasternak model. The random system equations are derived using the principal of virtual work. A finite element method is used for spatial descretization of the laminate with a reasonable accuracy. A perturbation technique has been the first time successfully combined with direct iterative technique by neglecting the changes in nonlinear stiffness matrix due to random variation of transverse displacements during iteration. The numerical results for the second order statistics of post buckling loads are obtained. A detailed study is carried out to highlight the characteristics of the random response and its sensitivity to different foundation parameters, the plate thickness ratio, the plate aspect ratio, the support condition, the stacking sequence and the lamination angle on the post buckling response of the laminate. The results have been compared with existing results and an independent Monte Carlo simulation.