  We investigate the phase diagram of a spin-$1/2$ Ising model on a cubic
lattice, with competing interactions between nearest and next-nearest neighbors
along an axial direction, and fully connected spins on the sites of each
perpendicular layer. The problem is formulated in terms of a set of
noninteracting Ising chains in a position-dependent field. At low temperatures,
as in the standard mean-feild version of the Axial-Next-Nearest-Neighbor Ising
(ANNNI) model, there are many distinct spatially commensurate phases that
spring from a multiphase point of infinitely degenerate ground states. As
temperature increases, we confirm the existence of a branching mechanism
associated with the onset of higher-order commensurate phases. We check that
the ferromagnetic phase undergoes a first-order transition to the modulated
phases. Depending on a parameter of competition, the wave number of the striped
patterns locks in rational values, giving rise to a devil's staircase. We
numerically calculate the Hausdorff dimension $D_{0}$ associated with these
fractal structures, and show that $D_{0}$ increases with temperature but seems
to reach a limiting value smaller than $D_{0}=1$.
