由颗粒组成的物质在现实世界中广泛存在,这种物质通常包含数量巨大且相互作用的离散颗粒,单个颗粒的材料行为可能比较简单,但是它们组成的物质由于个体之间的接触作用和不连续性使得整体材料行为非常复杂,既具有固体的结构特征又表现出流体的流动能力,基于连续介质的经典力学理论难以适用。离散元法是研究颗粒体物质的有效数值手段,而离散元法的核心在于颗粒之间的接触算法。在一个由众多颗粒构成的大体积离散物质中,接触搜索和对接触作用力的更新耗费大量的计算时间,离散元法的计算效率和计算精度取决于接触算法,对于复杂的三维几何形状颗粒,颗粒间的接触搜索过程通常占整个计算时间的80%以上。因此,三维离散元接触算法在颗粒体物质的计算仿真中占据极其重要的地位。本文首先提出了一种离散颗粒三维几何模型的生成算法,能够产生任意多面体颗粒单元,通过引入了包围球的概念,在有限空间域内能够产生无初始穿透的颗粒体物质模型。随后,提出了基于剖分技术的多面体颗粒单元体积、质心位置和转动惯量的计算方法。随后,本文根据多面体颗粒顶点、边、面之间的几何拓扑关系,以及球形单元的几何特征,研究建立了三维球形单元颗粒间、多面体单元颗粒间、以及颗粒与周围环境物体之间的接触搜索算法,给出接触穿透量、接触点相对位置、接触点相对速度、接触力的计算公式。最后,对所提出的颗粒生成和接触算法进行了编程实现,计算了球形颗粒材料、多面体颗粒材料以及混合形状颗粒材料在有限边界的方形盒内的堆积问题,验证算法的正确性。