In this paper, we provide an analytical frame work for the following multi-valued parabolic variational inequality in a cylindrical domain \({Q = \Omega \times (0, \tau)}\) : Find \({{u \in K}}\) and an \({{\eta \in L^{p'}(Q)}}\) such that $$\eta \in f(\cdot,\cdot,u), \quad \langle u_t + Au, v - u\rangle + \int_Q \eta (v - u)\,{\rm d}x{\rm d}t \ge 0, \quad \forall \, v \in K,$$ where \({{K \subset X_0 = L^p(0,\tau;W_0^{1,p}(\Omega))}}\) is some closed and convex subset, A is a time-dependent quasilinear elliptic operator, and the multi-valued function \({{s \mapsto f(\cdot,\cdot,s)}}\) is assumed to be upper semicontinuous only, so that Clarke’s generalized gradient is included as a special case. Thus, parabolic variational–hemivariational inequalities are special cases of the problem considered here. The extension of parabolic variational–hemivariational inequalities to the general class of multi-valued problems considered in this paper is not only of disciplinary interest, but is motivated by the need in applications. The main goals are as follows. First, we provide an existence theory for the above-stated problem under coercivity assumptions. Second, in the noncoercive case, we establish an appropriate sub-supersolution method that allows us to get existence, comparison, and enclosure results. Third, the order structure of the solution set enclosed by sub-supersolutions is revealed. In particular, it is shown that the solution set within the sector of sub-supersolutions is a directed set. As an application, a multi-valued parabolic obstacle problem is treated.