The topology of magnetic islands with m = 1 dominant poloidal mode number, as observed in the Frascati Tokamak Upgrade, cannot be explained by a single harmonic perturbation of an axisymmetric configuration. We present an analytical, nonlinear model that solves a helical modification of the Grad–Shafranov equation in cylindrical approximation and gives results comparable with experiments.