Hertzian theory has so far been employed extensively for modeling dynamic contact interactions between rough surfaces. First significant theoretical work on contact dynamics was by Nayak (1972), who developed some of the theoretical groundwork necessary for detailed physical explanations of experimentally observed phenomena in vibratory point contact by modelling a single-degree-of-freedom dynamic system and obtained analytical solutions by using a single term harmonic balancing method (HBM). Hess and Soom (1991a, 1991b) studied nonlinear vibrations at a Hertzian contact excited by the dynamic component of an externally applied normal load by using the method of multiple scales (MMS) (Nayfeh and Mook, 1979). The authors obtained solutions for both the average and instantaneous contact deflections and quantified the amount of friction reduction due to contact vibration. Perret-Liaudet(1997, 1998) investigated a sphere-plane (Hertzian) contact problem for sub-harmonic and super-harmonic resonance of order two using MMS and expanding the nonlinearityby third order Taylor series. The condition for contact loss was taken into consideration for the first time in this work. Sabot et al. (1998) studied the non-linear vibrations of a sphere–plane contact excited by a normal load equal to the sum of a static load and a harmonic load and predicted the contact natural frequency, frequency contents and softening behavior using numerical and analytical methods. Perret-Liaudet and Sabot (1999) analyzed the sphere-plane model of Hertzian contact by numerical methods and confirmed the significance of clearance at the contact on super-harmonic and sub-harmonic resonances, which lead to vibro-impact responses. Rigaud and Perret-Liaudet (2003) and Perret-Liaudet and Rigaud (2003) investigated dynamic behaviour of a preloaded double sphere–plane Hertzian contact under sinusoidal excitation and random normal excitation and solution was obtained using the shooting method in conjunction with parametric continuation technique. The authors also analyzed the dynamic response of an impacting Hertzian contact subjected to an order-2 subharmonic excitation and an order-2 superharmonic excitation [Perret-Liaudet and Rigaud, (2006, 2007)]. Ma et al. (2007) studied the sphere-plane contact model with possible contact loss experimentally and analytically using a multi-term harmonic balancing method. Tian and Xie(2008) analyzed the dynamic contact stiffness at the interface between a vibrating rigid sphere and a semi-infinite transversely isotropic viscoelastic solid with a harmonic force superimposed onto a static compressive force. Xiao et al. (2011) used exact method, MMS and HBM to determine the natural frequency of a mass interacting with a nonlinear contact stiffness. Zili et al. (2013) employed fractal contact model in conjunction with thin layer elements to capture the effects of lightly clamped joint interfaces on the dynamics of assembled structures taking the effects of surface topography, preload, and material properties into account. Xiao et al. (2015) modelled an elastic fractal surface contacting with a rigid flat surface. The authors analyzed its force-displacement relation and studied the free and forced vibration responses of the system.