To investigate thermal stress within nozzle thermal protection structures,an axisymmetric anisotropic thermal stress solution procedure was developed using the finite volume method.The program capitalizes on the formal similarity between the solid and fluid control equations,using discretization methods from Computational Fluid Dynamics(CFD)for discretizing the solid ther-moelastic equations.The volume source method was utilized to handle axisymmetric aspects,integrating the elastic matrix and area vector to account for material anisotropy.The mechanical parameters at material interfaces were interpolated using the central differ-ence method.An iterative approach was adopted to tackle the nonlinearity arising from temperature-dependent material mechanical properties.To validate the accuriency and efficacy of the developed procedure,four benchmark cases were examined,and computa-tional results were compared against that by using the Finite Element Method(FEM).The results show that the axisymmetric ther-mal stress solution procedure for thermal protection materials,within a cohesive framework,can accurately resolve the thermo-me-chanical coupling issues in anisotropic composite materials,the displacement magnitude at the nozzle's expansion section tail is 2.65 times greater deformation under thermo-mechanical loads compared to pressure loads alone.