This paper presents results from thermal buckling analysis and optimization of stiffened composite panels with variable angle tow (VAT) laminates and curvilinear stiffeners. Considering the meshing difficulties for curvilinearly stiffened VAT laminates due to both spatially dependent fiber path orientation and arbitrarily shaped stiffeners in the traditional finite element analysis, the present work models the plate and stiffeners separately without the need to place finite element nodes along the stiffener/plate and stiffener/stiffener interfaces. The displacement compatibility conditions are enforced at these interfaces in the finite element analysis. Convergence and verification study results show that the present method can accurately predict the thermal buckling behaviors of VAT laminates and curvilinearly stiffened VAT laminates. Parametric studies show that either VAT laminates or curvilinear stiffeners alone can improve the buckling performance as compared with using straight fiber laminates and straight stiffeners. VAT laminates and curvilinear stiffeners are found to redistribute the in-plane stress resultants and tailor the buckling mode shape, respectively, to improve the total buckling response. Optimization studies using curvilinear stiffeners and VAT laminates together for maximizing buckling temperature for a stiffened plate are conducted.