The mechanical behavior of sand is very complex, and depends on factors including confining pressure, density, and drainage condition. A soil mass can be contractive or dilative when subjected to shear loading, and eventually reaches an ultimate state, referred to as the critical state in soil mechanics. Conventional approach to explore the mechanical behavior of sand mainly relies on the experimental tests in laboratory. This paper gives an alternative view to this subject using discrete element method (DEM), which has attracted much attention in recent years. The implementation of the DEM is carried out by a series of numerical tests on granular assemblies with varying initial densities and confining pressures, under different test configurations. The results demonstrate that such numerical simulations can produce correct responses of the sand behavior in general, including the critical state response, as compared to experimental observations. In addition, the DEM can further provide details of the microstructure evolutions during shearing processes, and the resulting induced anisotropy can be fully captured and quantified in the particle scale. The mechanical behavior of sand is very complex, and depends on factors including confining pressure, density, and drainage condition. A soil mass can be contractive or dilative when subjected to shear loading shear, and ultimately reaches an ultimate state, referred to as the critical state in soil mechanics. Conventional approach to explore the mechanical behavior of sand mainly relies on the experimental tests in laboratory. This paper gives an alternative view to this subject using discrete element method (DEM), which has attracted much attention in recent years. The implementation of the DEM is carried out by a series of numerical tests on granular assemblies with varying initial densities and confining pressures, under different test configurations. The results demonstrate that such numerical simulations can produce correct responses of the sand behavior in general, including the critical state response, as compared to experimental observations details of the microstructure evolutions during shearing processes, and the resulting induced anisotropy can be fully captured and quantified in the particle scale.