A quasi-static homogeneous drained triaxial compression test on cohesionless sand under constant lateral pressure was simulated using a three-dimensional DEM model. Grain roughness was modelled by means of symmetric clusters composed of rigid spheres imitating irregular particle shapes. The effect of grain roughness on shear strength, dilatancy, kinetic, elastic and dissipated energies was numerically analyzed. Some numerical results were compared with available experimental results.