Turbulent channel and Couette flows are studied numerically by using an anisotropic K- epsilon model. A feature of this model lies in an anisotropic expression for the Reynolds stress and the deviation of the Reynolds stress from an isotropic eddy-viscosity representation is incorporated. Only one kind of wall damping function is introduced to impose the no-slip boundary condition on solid walls. The results obtained show that turbulence quantities of channel and Couette flows are in good agreement with experimental data and numerical results from large-eddy simulation. The anisotropy of turbulent intensities, which the usual K- epsilon model cannot predict, is well reproduced.