Over the past decade, Building Isolation Systems (BIS) have gain significant relevance due to their ability to reduce horizontal acceleration and interstory drifts in structures. Since the 1950s, researchers have focused on developing numerical models to simulate the dissipative behavior of High Damping Rubber Bearings (HDRB) in parallel major earthquake events have highlighted the need for BIS devices in medium and large-scale infrastructure, accentuating the need for further research into accurate models and adding the pressing interest in variability of mass-produced HDRB parameters. This study presents initial results from an identification process using two numerical models, validated using experimental tests at the SISMALAB laboratory. The experimental data involved eight samples subjected to compression forces and horizontal displacement. Optimal values were obtained through a Genetic Algorithm optimization process, minimizing discrepancies between experimental and numerical response. Preliminary variability analysis was conducted on data from 20 independent iterations over the eight samples.