This paper deals with the investigation of nonlinear free vibration behavior of elastically supported carbon nanotube reinforced composite (CNTRC) beam subjected to thermal loading with random system properties. Material properties of each constituent’s material, volume fraction exponent and foundation parameters are considered as uncorrelated Gaussian random input variables. The beam is supported by a Pasternak foundation with Winkler cubic nonlinearity. The higher order shear deformation theory (HSDT) with von-Karman non-linearity is used to formulate the governing equation using Hamilton principle. Convergence and validation study is carried out through the comparison with the available results in the literature for authenticity and accuracy of the present approach used in the analysis. First order perturbation technique (FOPT),Second order perturbation technique (SOPT) and Monte Carlo simulation (MCS) methods are employed to investigate the effect of geometric configuration, volume fraction exponent, foundation parameters, distribution of reinforcement and thermal loading on nonlinear vibration characteristics CNTRC beam.The present work signifies the accurate analysis of vibrational behaviour influences by different random variables. Results are presented in terms of mean, variance (COV) and probability density function (PDF) for various aforementioned parameters.