Granular system commonly encountered in industry or nature is comprised of non-spherical grains. Comparing with spherical particles, high discretization and interlocking among non-spherical particles can effectively dissipate the system energy and improve the buffer capacity. The superquadric element based on continuous function envelop can form the geometric shape of irregular particles accurately, and then contact collision action between particles can be calculated easily. In this paper, we provide a comprehensive introduction to particle-particle and particle-boundary contact collision. In addition, considering different shapes and surface curvatures under various contact patterns between super-quadric particles, the linear contact force model cannot be applied to the accurate calculation of the contact force, and a corresponding non-linear viscoelastic force model is developed. In this model, the equivalent radius of curvature at a local contact point is adopted to calculate the normal contact force, and the tangential contact force is simplified based on the contact model of spherical elements. To examine the validity of the algorithm and this model, we compare the discrete element analytical results with the analytical results for a single cylinder impacting a flat wall and the previous experimental results for spherical granular material under impact load, and this method is verified by good agreement between the simulated results and the previous experimental results. According to the aforementioned method, we study the buffer capacity of non-spherical particles under impact load by the discrete element method, and the influences of granular thickness and particle shapes on the buffer capacity are discussed. The results show that a critical thickness Hc is obtained for different particle shapes. The buffer capacity is improved with increasing the granular thickness when H Hc, but is independent of the granular thickness and particle shapes when H Hc. Moreover, the impact peak and initial packing fraction increase significantly with increasing the blockiness. Rectangular particles account for the highest packing fraction, and the packing fraction of cylindrical particles is higher than the packing fraction of spherical particles. Therefore, Rectangular particles are more likely to form dense face-face contacts and ordered packing structures with high packing fraction. These denser packings prevent the particles from their relatively moving, and thus reducing the buffering capacity of the particles. Furthermore, the impact peak and initial packing fraction decrease with increasing or reducing the aspect ratio of cylindrical particles and the aspect ratio of rectangular particles. The aspect ratio of particle can be used to adjust the dense packing structure and reduce the stability of the system. It means that the particles have more effective buffer capacity for the non-spherical particle system.