An important consideration in design involving high temperature variation is the determination of the thermal stresses developed. The numerical solution for thermoelastic transient response of orthotropic cylinder subjected to a constant temperature at the surface is presented. The thermoelastic equations with one relaxation time developed by Lord and Shulman with uncoupled thermoelasticity assumption are used in the present work. The hyperbolic heat conduction model is used for the prediction of the temperature history. Thermally induced displacement and stresses are determined. A numerical method based on implicit finite difference scheme is used to calculate the temperature, displacement, and stress distributions within the cylinder Numerical examples for orthotropic, transverse isotropic, and isotropic cylinders were carried out for the stresses. Furthermore, the results of the numerical solution and the exact solution at the steady state condition are compared.