In strongly nonlocal nonlinear media, the evolution of the optical beam is governed by Snyder-Mitchell model. By applying the variational method, we obtain an approximate analytical solution of Snyder-Mitchell model, i.e., the rotating azimuthon, which has an inherent rotation. Under certain conditions, the rotating azimuthon can be reduced to the radially symmetric optical vortex soliton.