We have that $$\tag{1}\int\nabla u_\epsilon\nabla\varphi+\int\beta(u_\epsilon)\varphi=f\varphi,\ \forall\ \varphi\in H_0^1$$
By taking $\varphi=u_\epsilon$, we conclude $$\tag{2}\int|\nabla u_\epsilon|^2\leq\int fu_\epsilon$$
which implies by Holder inequality that $\|\nabla u_\epsilon\|_2$ is bounded, therefore, we can suppose that $u_\epsilon \rightharpoonup u$ in $H_0^1$. Moreover, we can suppose that $u_\epsilon \to u$ in $L^2$. We combine these convergences with $(2)$ to get $$\tag{3}\int |\nabla u|^2\leq \int fu$$
On the other hand, define $\mathcal{K}=\{w\in H_0^1:\ w\geq 0\}$. As you can easily verify, $\int\beta(u_\epsilon)w\leq 0$ in $\mathcal{K}$, hence, from $(1)$, we conclude that $$\tag{4}\int\nabla u_\epsilon\nabla\varphi\geq \int f\varphi,\ \forall\ \varphi\in \mathcal{K}$$
Again, we use the convergences we have, to conclude from $(4)$ $$\tag{5}\int\nabla u\nabla\varphi\geq\int f\varphi,\ \forall\ \varphi\in\mathcal{K}$$
We conclude from $(3)$ and $(5)$ the desired inequality. To conclude that $u\in\mathcal{K}$ and $u$ is unique, you can argue like this:
I - $u$ is the minimizer of some functional $I:\mathcal{K}\to\mathbb{R}$, what is $I$?
II - In this particular setting, for $u$ to minimizes $I$ is equivalently for $u$ to satisfies the variational inequality.
III - The solution of the optimization problems is unique.
As @RayYang suggested this book is a good one to understand it better, however I would like to point out that the main argument here can be found in any good book of convex analysis.
It is worth to note that what we are proving here is that $-I'(u)\in \mathcal{N}_{\mathcal{K}}(u)$, where $\mathcal{N}_{\mathcal{K}}(u)$ is the normal cone of $u$ with respect to $\mathcal{K}$ and when $I$ is a convex differentialbe function, II is valid, i.e. $-I'(u)\in \mathcal{N}_{\mathcal{K}}(u)$ if and only if $u$ minimizes $I$.
arrayenvironment is not usually used here,alignandcasesare more widely used in MSE, also I added the problem number. – Shuhao Cao Jul 04 '13 at 19:33