In the first part of this paper, the eddy-viscosity representation for Reynolds stress including the effects of curvature and rotation is derived by using the two-scale Direct Interaction Approximation. This representation is cubic functions of the mean velocity gradients, Coriolis and curvature. In the second part, the anisotropic k-epsilon turbulence model affected by curvature and rotation is proposed and applied to some turbulent flows with curvature and rotation. In the hypothetical case of equilibrium folw, the models are simplified and compared with the previous models.