## 1 Introduction
In addition, the method and regulations for measuring the shear modulus of elastomeric bearings for land and offshore are different, so the results of existing studies on bridge bearings cannot be applied to offshore applications. Therefore, the necessity of this research work is to analyze the shear modulus of offshore elastomeric bearings according to the procedure mentioned in EN 1337-3 regulation. Offshore elastomeric bearings must be physically tested with full-size products in accordance with EN 1337-3 regulation. However, the experimental set-up incurs high cost due to large capacity of the test equipment and a dedicated shear test device. On the other hand, the numerical simulation can analyze various conditions at low cost, but in order to accurately analyze elastomers, which are hyperelastic (nonlinear) materials, it is necessary to find suitable analysis conditions. Therefore, finding the optimum analysis conditions for the numerical simulation technique that reproduce corresponding experimental results is necessary. The most crucial step in predicting the nonlinear behavior of a hyperelastic material through numerical simulations is considering the correct nonlinear model. When performing numerical simulations for hyperelastic materials such as rubber, the first-order, two-parameter (2P) Mooney-Rivlin (M-R) model is mainly used, which has a good response in compression within 30% and tension within 200%. On the contrary, the error is substantial in tensions of 200% or more. In this case, the responsiveness can be improved by using the 5P and 9P M-R models for the material behavior [18; 19; 20]. In this study, three types of M-R models, such as 2P (first-order), 5P (second-order), and 9P (third-order), are used to increase the reliability of numerical simulation for offshore elastomeric bearings. The main objective of this research work is to propose a suitable M-R model for simulating the behavior of an elastomeric bearing under compression and shear tests by comparing and analyzing load-displacement graphs obtained with 2P, 5P, and 9 Pmodels with the experiment. Overall aim is to match the simulation results of elastomeric bearing's shear modulus values with the experimental results. The research procedure for this is as follows: First, actual-sized offshore elastomeric bearing specimens are fabricated, and the experiment was conducted according to the procedure specified in EN 1337-3 Annex F: Shear Modulus Test. Second, the numerical simulation of the elastomeric bearing with respect to the behavior model of hyperelastic material was performed for the same conditions as the experiment. Third, after obtaining a shear force-displacement graph in both the experiment and the numerical simulation, the shear modulus is calculated and compared by using the slope of the graph. Last, the shear modulus of the elastomeric bearing with respect to the shore hardness was verified by comparing the numerical simulation values with the values presented in EN 1337-3 Annex D. A 2500-ton press equipped with a shear test device with a capacity of 300-tons is used for the experiment, and ANSYS Mechanical APDL [21], a general-purpose finite element analysis program, is used for numerical simulation.