Fixed boundary three-dimensional anisotropic pressure equilibria are computed that are consistent with a family of distribution functions that load particles preferentially towards v parallel to similar or equal to v and model the effects of balanced tangential neutral beam injection. The total pressure surfaces remain coupled with the magnetic surfaces even when the parallel beam pressure greatly exceeds its perpendicular counterpart. The anisotropy factors and the Shafranov shift, which increases with the anisotropy, are calculated in a currentless two-field period quasiaxisymmetric stellarator device.