Lead-rubber seismic isolation bearings (LRB) have been installed in a number of essential and critical structures, like hospitals, universities and bridges, in order to provide them with period lengthening and the capacity of dissipating a considerable amount of energy to mitigate the effects of strong ground motions. Therefore, studying the damage mechanics of this kind of devices is fundamental to understand and accurately describe their thermo-mechanical behavior, so that seismically isolated structures can be designed more safely. Hitherto, the hysteretic behavior of LRB has been modeled using 1) Newtonian mechanics and empirical curve fitting degradation functions, or 2) heat conduction theories and idealized bilinear curves which include degradation effects. The reason for using models that are essentially phenomenological or that contain some adjusted parameters is the fact that Newton’s universal laws of motion lack the term to account for degradation and energy loss of a system. In this paper, the Unified Mechanics Theory – which integrates laws of Thermodynamics and Newtonian mechanics – is used to model the force-displacement response of LRB. Indeed, there is no need for curve fitting techniques to describe their damage behavior because degradation is calculated at every point using entropy generation along the Thermodynamics State Index (TSI) axis. A finite element model of a lead-rubber bearing was constructed in ABAQUS, where a user material subroutine UMAT was implemented to define the Unified Mechanics Theory equations and the viscoplastic constitutive model for lead. Finite element analysis results were compared with experimental test data.