The shape of the free boundary arising from the solution of a variational inequality is controlled by the shape of the domain where the variational inequality is defined. Shape and topological sensitivity analysis is performed for the obstacle problem and for a regularized version of its primal-dual formulation. The shape derivative for the regularized problem can be defined and converges to the solution of a linear problem. These results are applied to an inverse problem and to the electrochemical machining problem.