In plasmas with strongly anisotropic distribution functions T kb =T ?b 1 a Harris-like collective instability may develop if there is sufficient coupling between the transverse and longitudinal degrees of freedom. Such anisotropies develop naturally in accelerators and may lead to a deterioration of beam quality. This paper extends previous numerical studies [E. A. Startsev, R. C. Davidson, and H. Qin, Phys. Plasmas 9, 3138 (2002)] of the stability properties of intense non-neutral charged particle beams with large temperature anisotropy T ?b T kb to allow for nonaxisymmetric perturbations with @=@ Þ 0. The most unstable modes are identified, and their eigenfrequencies, radial mode structure, and nonlinear dynamics are determined. The simulation results clearly show that moderately intense beams with s b ! ! 2 pb =2 2 b ! 2 ? * 0:5 are linearly unstable to short-wavelength perturbations with k 2 z r 2 b * 1, provided the ratio of longitudinal and transverse temperatures is smaller than some threshold value. Here,! ! Both the simulations and the analytical theory predict that the dipole mode (azimuthal mode number m 1) is the most unstable mode. In the nonlinear stage, tails develop in the longitudinal momentum distribution function, and the kinetic instability saturates due to resonant wave-particle interactions.