This investigation aims to study the random stresses in a functionally graded (FG) thick hollow cylinder with uncertain material properties subjected to mechanical shock loading using a hybrid numerical method. The mechanical properties are considered to vary across thickness of FG cylinder as a nonlinear power function of radius. The stresses are obtained by solving Navier equation and using Galerkin finite element and Newmark finite difference methods. The Monte Carlo simulation is used to generate the random mechanical properties for the problem. The failure probabilities and time history analysis of stresses are determined for various coefficient of variation considering various grading patterns of mechanical properties. The presented hybrid numerical method is effective, with high capability for stochastic analysis of dynamic and transient analysis of FG structures with various boundary conditions.