The goal of this paper is to study an evolution inclusion problem with fractional derivative in the sense of Caputo, and Clarke’s subgradient. Using the temporally semi-discrete method based on the backward Euler difference scheme, we introduce a discrete approximation system of elliptic type corresponding to the fractional evolution inclusion problem. Then, we employ the surjectivity of multivalued pseudomonotone operators and discrete Gronwall’s inequality to prove the existence of solutions and its priori estimates for the discrete approximation system. Furthermore, through a limiting procedure for solutions of the discrete approximation system, an existence theorem for the fractional evolution inclusion problem is established. Finally, as an illustrative application, a complicated quasistatic viscoelastic contact problem with a generalized Kelvin–Voigt constitutive law with fractional relaxation term and friction effect is considered.