Constituting the simplest generalization of spherical Coulomb crystals, assemblies of N equicharged particles confined by radial potentials proportional to the λth power of distance are amenable to rigorous analysis within the recently introduced shell model. Thanks to the power scaling of the confining potential and the resulting pruning property of the shell configurations (i.e., the lists of shell occupancies), the shell-model estimates of the energies and the mean radii of such assemblies at equilibrium geometries follow simple recursive formulas. The formulas greatly facilitate derivations of the first two leading terms in the large-N asymptotics of these estimates, which are given by power series in ξ(4/3) N(-2/3), where -(ξ/2) n(3/2) is the leading angular-correlation correction to the minimum energy of n electrons on the surface of a sphere with a unit radius (the solution of the Thomson problem). Although the scaled occupancies of the outermost shells conform to a universal scaling law, the actual filling of the shells tends to follow rather irregular patterns that vary strongly with λ. However, the number of shells K(N) for a given N decreases in general upon an increase in the power-law exponent, which is due to the (λ + 1)(2) ξ(2) dependence of shell capacities that roughly measure the maximum numbers of particles sustainable within individual shells. Several types of configuration transitions (i.e., the changes in the number of shells upon addition of one particle) are observed in the crystals with up to 10,000 particles and integer values of λ between 1 and 10, but the rule |K(N + 1)-K(N)| ≤ 1 is found to be strictly obeyed.