###### Abstract  
This article presents a numerical study that aims to explore the impacts of the stiffness of elastomeric bearings on the dynamic behavior of railway bridges under train-induced vibrations. For this purpose, a finite element code that considers vehicle-bridge interaction using a coupled approach was developed. The software was validated by comparing the numerical response to the analytical solution. The numerical analysis of single- and multi-span bridges with varying bearing stiffness values under passenger trains showed the interplay between bearing stiffness, its impact on the natural frequency of the bridge and the loading frequency. It is demonstrated that the amplitude of the maximum acceleration on the bridge depends heavily on the stiffness of the bearings. Furthermore, the bearing stiffness significantly impacts the location of the maximum acceleration on the bridge. The results of the extensive numerical analyses improve the understanding of the impact of the bearing stiffness on the dynamic behavior of bridges and highlight the importance of quantifying the boundary conditions correctly for reliable estimation of dynamic response of railway bridges under train-induced vibrations. a 202  
also utilized different modeling approaches to carry out a parametric analysis. In their study, vehicle damping was found to have a negligible influence on the bridge response. [PERSON] et al. [17] analyzed the dynamic behavior of a railway bridge and compared the results to the field measurements. In another aspect, [PERSON]. [18; 19] investigated the vibration characteristics of the bridge-train system using a complex numerical model and field tests. Arguably, the most detailed parametric studies on this topic can be found in [20; 21]. In both studies, the train-bridge interaction was modeled using both simple moving load models and more detailed interaction models. Impact of several parameters on the reduction in acceleration demands using detailed train bridge interaction models were evaluated and discussed. More recently, [PERSON] et al. [22] compared the suitability of uncoupled and coupled methods for solving vehicle-bridge interaction problems. Through numerical analysis of different vehicle and bridge parameters such as mass ratio and frequency ratio, [PERSON]. [22] demonstrated the error introduced by the uncoupled method and how it increases with an increase in the vehicular dynamic force. These studies, similar to most of the others found in literature [9; 15; 17; 18; 19; 23], are based on simply-supported beams and ignores the effect of the bearing stiffness on the vehicle-bridge interaction problem. If the dynamic response of the bridge under moving train loads is sensitive to the bearing stiffness, deviations of these characteristics from the assumed values during the design phase may significantly alter the dynamic response of the bridge during its lifetime. Despite this potentially significant impact, few studies have focused on the interaction between the bearing stiffness and dynamic response. Arguably, [PERSON] and [PERSON] [24] and [PERSON] et. al. [25] present the most detailed work on the subject. In this context, [PERSON] et al. [25] examine the mechanisms of dynamic behaviour regarding resonance and cancellation when using elastic bearings. By implementing a moving load model with equal distance between each axle, a numerical study was conducted with verification from a field test. Although the use of elastic bearings for railway bridges have several benefits, [PERSON] et al. [25] also bring attention to the prevention of transmitting the vehicle-induced vibrations to the ground. This may lead to the accumulation of induced vibrations which promotes greater fatigue of the bridge. Similar studies were also reported for single span beams in [26; 27; 28] while only a handful studies [29; 30] focus on the impact of elastomeric bearings on the behavior of multi-span bridges.