In this thesis, explicit finite element code is used to simulate the nonlinear behavior of the bridge bearing system. The nonlinear behavior of rubber bearing was investigated in detail. Three types of material models in the LS- DYNA coder were evaluated by comparing with the experimental results of the compression and shearing tests of rubber bearing. It turns out that the Hyperviscoelastic Rubber model can simulate the performance better than the Kinematic Hardening Plastic Material model and the Mooney-Rivlin Rubber model. The material parameters obtained from those validation simulations will be used to conduct the numerical analysis for other types of design of rubber bearing device. The performance curves of the rubber bearing device in the local level can be linked with the global structural analysis of bridge structure under seismic loading.