在考虑堆石料真实颗粒形状的基础上，将每个颗粒用一个多边形离散单元表示，基于一种线性搜索算法对多边形-多边形之间的接触详情进行检索计算，引入基于势能原理的法向多边形-多边形接触模型及切向库仑摩擦模型，形成了一种多边形离散元计算方法，并开发了相应的PDEM程序。可以从颗粒尺度层面展示颗粒之间的相互作用以及每个颗粒的位移和转动，帮助进一步揭示颗粒的细观特性（形状、大小、材料特性等）对堆石料的宏观强度和变形的影响。最后采用PDEM程序对某粗粒料的二维模型试验进行了数值模拟，得到了与室内试验一致的应力变形规律，展现了其方法和程序用于研究堆石料的有效性。

Every single particle is simulated by a polygon discrete element to capture the realistic shape of rockfill materials. A polygon discrete element method (PDEM) is developed by adopting a simple contact detection program and a polygon/polygon contact model. A linear program is adopted to detect the contact details between polygons. Then the normal contact force is calculated by a potential energy based polygon/polygon normal contact model, and a polygon discrete element calculation method is formed. Based on this method, a program called PDEM is developed to study the interaction between particles and both the translational and rotational motion of every particle from the microscopic view. The effect of micro-properties (e.g. particle shape, size, material properties et al.) on the macro-strength and deformation is enabled. A two-dimensional model test of a coarse aggregate was carried out by PDEM program. The stress and deformation laws consistent with the lab experiment were obtained, and the method and procedure were used to study the effectiveness of the rockfill.