We study energy minimization of a continuum Landau-de Gennes energy functional for nematic liquid crystals, in three-dimensional axisymmetric domains and in a restricted class of $mathbb{S}^1$-equivariant (i.e., axially symmetric) configurations. We assume smooth and nonvanishing $mathbb{S}^1$-equivariant (e.g. homeotropic) Dirichlet boundary conditions and a physically relevant norm constraint (Lyuksyutov constraint) in the interior. Relying on results in cite{DMP1} in the nonsymmetric setting, we prove partial regularity of minimizers away from a possible finite set of interior singularities lying on the symmetry axis. For a suitable class of domains and boundary data we show that for smooth minimizers (torus solutions) the level sets of the signed biaxiality are generically finite union of tori of revolution. Concerning nonsmooth minimizers (split solutions), we characterize their asymptotic behavior around any singular point in terms of explicit $mathbb{S}^1$-equivariant harmonic maps into $mathbb{S}^4$, whence the generic level sets of the signed biaxiality contains invariant topological spheres. Finally, in the model case of a nematic droplet, we provide existence of torus solutions, at least when the boundary data are suitable uniaxial deformations of the radial anchoring, and existence of split solutions for boundary data which are suitable linearly full harmonic spheres.