Numerical modeling of anisotropic granular media