Every moment of every day, at least hundreds of thousands of tires roll across a surface throughout the world. Tires are indisputably important in our daily life. The tire's primary component is rubber, which consumes energy when it rotates on a substrate due to the viscoelastic material's internal friction: a phenomenon referred to as rolling resistance. The interaction between the tire and the road surface is one of the most intricate and crucial phenomena in an automobile, because it is responsible for creating forces, moments, and deformation in the tire. Additionally, the road's roughness interacts with the tire and contributes significantly to its performance. This dissertation aims to develop a comprehensive physics-based model for predicting the rolling resistance of a viscoelastic material due to dynamic deformations caused by tire rotation using an analytical approach. The model was developed by proposing a Gaussian wave function propagating across a tire circumference's viscoelastic medium. The wave function was selected to describe the displacement field produced by tire-road interaction. Additionally, by adopting a multi-scale modeling technique, the model was upgraded to estimate rolling resistance while taking into account surface roughness at all length scales, from macroscopic to microscopic. Additionally, another mathematical model was developed using the Fourier series approach to evaluate the steady-state stress response and energy dissipation for any harmonic and non-harmonic periodic strain signals. Additionally, the dissertation strove to build a continuum damage mathematical model using a combined testing/modeling methodology to predict the aging of Styrene-Butadiene Rubber (SBR) after continuous exposure to the atmosphere. The obtained model was developed through the implementation of optimization techniques while formulating a mathematical model, which was then combined with a physics-based model to predict rolling resistance while taking into account rubber aging. Calibration of hyperelastic and viscoelastic material models with testing data was performed using an optimization technique that yielded sufficient results. The results of all mathematical models obtained in this dissertation are reported subsequently. The stress response of a viscoelastic material under harmonic and non-harmonic strain input yielded good agreement with the FEA model obtained using ABAQUS. The rolling resistance behavior under various operating conditions, including texture and aging effects, was reported, and the results aligned with the experimental results found in the literature.