Abstract Delamination and matrix cracking are common damage modes in composite laminates, significantly impacting their structural integrity and performance. Under elevated temperatures, the collision of anisotropic thermal expansion and nonlinear stiffness responses affect the onset of thermal buckling in flying vehicles. This study implements a multi-scale analysis framework, incorporating finite element modeling and analytical methods to obtain the reduction in mechanical properties of the composite laminate and its stiffness degradation based on an effective stiffness model due to pre-existing matrix cracking in plies. Based on non-linear equations, obtained from the principle of virtual work, in combination with first-order shear deformable theory and von Kármán strain-displacement relations, the thermal buckling temperature, post-buckling behavior of composite laminates are investigated. This work paves the way for developing a model to predict the influence of the laminate stiffness, delamination size, transverse crack location, stacking sequence, and plies thickness on the laminate.