The discrete element method (DEM) simulation of discs leads to typically low internal friction angles compared with real sand. The paper introduces a developed discrete element program NS2D

of which the subject is two-dimensional assemblages of elliptical particles. The basic contact models of the developed NS2D program are described; and the force-displacement formulations of the contacts of ellipse-ellipse and ellipse-wall have been deduced. Assemblages of elliptical particles with aspect ratio1.1:1 and 1.4:1 were generated with undercompaction method (UCM); then consolidated isotropically

and then sheared with biaxial compression at constant confining pressure

respectively. Their mechanical responses are both qualitatively and quantitatively in good agreement with the reported literature; and the simulated internal friction angles fall within the typical range of real sand. Therefore

the developed DEM program NS2D can be used to study the mechanical properties of real sand.