This study considers the buckling of orthotropic cylindrical thin shells with material non-homogeneity in the thickness direction, under torsion, which is a power function of time. The dynamic stability and compatibility equations are obtained first. Applying Galerkin's method then applying Ritz type variational method to these equations and taking the large values of loading parameters into consideration, analytic solutions are obtained for critical parameter values. Using those results, the effects of the periodic and power variations of Young's moduli and density, ratio of Young's moduli variations, loading parameters variations and the power of time in the torsional load expression variations are studied via pertinent computations. It is concluded that all these factors contribute to appreciable effects on the critical parameters of the problem in question.