Based on Hamilton's principle, a new accurate solution methodology is developed to study the torsional bifurcation buckling of functionally graded cylindrical shells in a thermal environment. The effective properties of functionally graded materials (FGMs) are assumed to be functions of the ambient temperature as well as the thickness coordinate of the shell. By applying Donnell's shell theory, the lower-order Hamiltonian canonical equations are established, from which the eigenvalues and eigenvectors are solved as the critical loads and buckling modes of the shell of concern, respectively. The effects of various aspects, including the combined in-plane and transverse boundary conditions, dimensionless geometric parameters, FGM parameters and changing thermal surroundings, are discussed in detail. The results reveal that the in-plane axial edge supports do have a certain influence on the buckling loads. On the other hand, the transverse boundary conditions only affect extremely short shells. With increasing thermal loads, the material volume fraction has a different influence on the critical stresses. It is concluded that the optimized FGM mixtures to withstand thermal torsional buckling are Si 3 N 4/SUS304 and Al 2 O 3/SUS304 among the materials studied in this paper.
