A method for numerical simulation of train induced track vibration and wave propagation in subgrade has been proposed. The method uses a mass to simulate the bogie of a train and considers the effect of rail roughness. For this method, rail roughness is considered as a randomly generated signal and a filter is used to block the undesired components. The method predicts the particle velocity around the track and can be applied to many kinds of railroad trackbeds including traditional ballast trackbed and modern Hot mix asphalt (HMA) trackbed. Results from ballast and HMA trackbeds are compared and effects of HMA layer on damping track vibration and dissipating wave propagation are presented. To verify the credibility of the method, in-track measurements were also conducted. Site measurements included performing geophysical tests such as spectral analysis of surface wave test and seismic refraction test to determine the subsurface conditions at the test site. Ballast and HMA samples were tested in the laboratory by resonant column test to obtain the material properties. Particle velocities were measured and analyzed in the frequency domain. Results from in-track tests confirm the applicability of the numerical method. The findings and conclusions are summarized and future research topics are suggested.