Second order sliding mode control synthesis is developed for underactuated mechanical systems, operating under uncertainty conditions. In order to locally stabilize an underactuated system around an unstable equilibrium, an output is specified in such a way that the corresponding zero dynamics is locally asymptotically stable. Then, the desired stability property of the closed-loop system is provided by applying a quasihomogeneous second order sliding mode controller, driving the system to the zero dynamics manifold in finite time. Although the present synthesis exhibits an infinite number of switches on a finite time interval, it does not rely on the generation of first order sliding modes, while providing robustness features similar to those possessed by their standard sliding mode counterparts. A second order sliding mode appears on the zero dynamics manifold which is of co-dimension greater than the control space dimension. Performance issues of the proposed synthesis are illustrated in numerical and experimental studies of a cart-Pendulum system.