The preferred orientation of ice crystals in polar ice sheets developsas a result of intracrystalline slip. Polycrystalline plasticity theory has been successfully used to simulate development of fabrics in rocks. In this study, we present a simpletheory for plastic deformation and fabric evolution. Each crystal within the aggregateis assumed to deform only by basal glide, and recrystallization processes are not directly taken into account. We have adopted the uniform stress bound, that is, stress is supposed to be homogeneous in the polycrystal. Simulated fabrics in uniaxial compression and tension are similar to those observed in ice sheets. In simple shear,simulations show that c axes rotate toward the direction of the maximal compression,in contradiction to the intuitive conclusion that stable orientations correspond to the easy glide position. But, these results are in agreement with predictions of Taylor theory and self-consistent approaches. We suggest that effects of vertical compression and/or rotation recrystallization are of great importance for fabric evolution in polarice sheets. The proposed model can reproduce viscosities of anisotropic ice samples tested in compression or tension but it is not able to reproduce the low viscosity of ices with a single-maximum fabric when tested in simple shear