Constraints are all that' needed. Just make sure the generalized coordinates are selected so that the constraints ie in this case use spherical coordinates. So ρ= constant and all the derivatives with respect to ρ are zero. You can calculate the actual constraint force after you solve the equations of motion.
Jim