The coupled vehicle-track system (CVTS) dynamics have been extensively investigated for decades. However, the calculation accuracy of prevailing vehicle-track coupling models needs to be improved in the high frequency range due to the inappropriate model simplification and neglect of material nonlinearity. In this study, we propose a refined numerical model of the CVTS that considers the nonlinear properties of the railpads and primary suspension using the fraction derivative Zener model. Furthermore, we more realistically simulate the wheelset, rail and railpad configuration with the elastic axle, solid finite element and surface-support models, respectively, and improve the computation efficiency by employing the mode superposition method. The results demonstrate that the refined CVTS model is more accurate than the classical model in simulating vehicle-track coupling dynamics above 2 kHz. In particular, there are significant differences in the dynamic response of the elastic wheelset model compared to the rigid model over a broad frequency range, with an 11% difference in the bogie acceleration response at the first dominant frequency. When the railpads are modeled using the surface-support model, the rail acceleration differences exceed 41% near 1 kHz and 44% near 2650 Hz, compared to the point-support model. Additionally, the rail response at various locations across the rail cross section can be calculated using the finite element method in this refined model. Overall. the proposed CVTS model provides high accuracy and efficiency for random vibration analysis, especially in the high frequency domain.