采用离散元等数值模拟方法研究粗颗粒土体材料的细观力学行为,是近年来岩土力学的热点研究课题.提出了一种基于计算接触力学的粗颗粒土体材料细观力学行为模拟分析新方法.该法将土体颗粒剖分成一定数量的单元,通过计算接触力学方法模拟颗粒间的多体接触行为,是一种基于有限元变分原理的隐式数值求解方法,具有严格的理论基础.和离散元法方法相比,该方法在描述颗粒本身的力学特性方面具有优势,可以计算得到颗粒内部的应力分布情况,从而为颗粒破碎等细观行为的模拟计算提供依据.利用所发展的粗颗粒土体多体接触有限元计算程序,进行了不同围压的常规三轴数值试验,模拟计算结果符合粗颗粒土体材料三轴试验的一般规律,初步验证了所提出的计算方法的适用性.