Since a lot of engineering problems are along with uncertain parameters, stochastic methods are of great importance for incorporating random nature of a system property or random nature of a system input. In this study, the stochastic dynamic analysis of soil mass is performed by finite element method in the frequency domain. Two methods are used for stochastic analysis of soil media which are spectral decomposition and Monte Carlo methods. Shear modulus of soil is considered as a random field and the seismic excitation is also imposed as a random process. In this research, artificial neural network is proposed and added to Monte Carlo method for sake of reducing computational effort of the random analysis. Then, the effects of the proposed artificial neural network are illustrated on decreasing computational time of Monte Carlo simulations in comparison with standard Monte Carlo and spectral decomposition methods. Numerical verifications are provided to indicate capabilities, accuracy and efficiency of the proposed strategy compared to the other techniques.