The solution to a multivalued stochastic differential equation is first strongly approximated by its Yosida approximation. Then this approximation is itself approximatively simulated by a semi-implicit Euler scheme and the strong order of convergence is evaluated in terms of the time step.