Based on the generalized Huygens Fresnel integral, we derive the analytical formula of the cross-spectral density of a twisted partially coherent array beam propagating in non-Kolmogorov anisotropic turbulence, and investigate the changes in orbital angular momentum (OAM). The results show that the anisotropy of the turbulence causes different effects in horizontal and vertical directions. The spectral density distribution of twisted partially coherent array beam in turbulence presents self-splitting and rotation, which combines the interesting effects of the twist phase and coherent structure. Although OAM is conserved, the spatial distribution of OAM flux density can be changed by changing the propagation distance, power and anisotropy of turbulence, and the modulation of the twist phase affects not only the magnitude of OAM but also its distribution. Our work is helpful for exploring new forms of OAM sources, and promote the application of free-space optical communications and optical field modulation.