Multi-layered elastomeric bearings (EB) have been applied in structural engineering as a mitigation solution for improving the seismic resilience of structural configurations. Such devices have been extensively verified for different loading conditions with laboratory tests and numerical simulations during the last decades. However, their behavior under the mutual effect of axial and shear load is still under assessment. This paper presents the results of 3D numerical simulations of EBs when subjected to large horizontal movements by considering the role of the shear modulus. The paper considers several numerical models of EBs. The main goal is to implement the proposed analytical formulation to numerical platforms that can describe the combined axial -shear macroscopic behavior of the EBs.