离散元模拟方法具有微观特征可控、研究成本低廉、结果可重复等优点,目前已成为从颗粒尺度研究土体宏观力学性质的重要手段。然而由于黏土颗粒形状不规则、粒间相互作用复杂等原因,使得黏土离散元研究相对较弱。本论文基于课题前期开发出的黏土离散元模拟程序框架,围绕黏土颗粒具有粒径变化范围大且易形成平行团聚体的特点,从黏土颗粒单元几何分析、接触检查、粒间力作用模型与由粒间力控制的黏土团聚体生成以及程序合理性检验等方面对黏土离散元方法展开深入研究。基于黏土特殊的微观结构形态、复杂的颗粒间力学作用效果以及颗粒单元的实际尺寸分析,本文建立较符合黏土颗粒实际的板状模型单元,并成功生成宽粒径范围的初始黏土试样。在此基础上,针对单元模型特征从接触形式判断、接触向量求解两个方面,提出新的粗略接触检测与接触信息求解方法,为确保黏土离散元模拟研究的可靠性、准确度及运算效率提供重要基础。通过对黏土SEM图像特征进行二维化处理,统计获得黏土微观试样粒径占比信息,并结合对双电层、范德华非接触作用力以及Born接触作用力研究,得出刚度系数随颗粒尺寸变化规律。综合结论,生成更符合黏土实际微观结构的宽粒径范围团聚体试样,并通过与赫兹理论对比分析,论证求解方法的合理性以及团聚体生成的可靠性。在建立宽粒径范围黏土团聚体初始试样并获取适用于特殊板状几何形态的接触信息检测与求解方法基础上:1)通过对颗粒间合力曲线的引力、斥力作用区间及峰值点等特征进行分析,添加与赫兹模型类似且与Born力连续的接触斥力模型,大幅增加试样峰值承载力,更进一步贴合实际情况下黏土的微观力学作用效果;2)通过对试样不同接触形式碰撞研究,利用二分法获取速度阈值v₀,在不改变单元运动规律的同时对其端点最大速度进行合理限制;3)通过设置多组速度清零模拟对比试验,综合验证Anandarajah与商用离散元程序所用系统能量消减处理方法的有效性,大幅提高模拟加载效率,并避免接触响应模块中颗粒单元交叉穿透等错误异常的发生。参考圆形离散元中的小球弹落模型与斜坡滚动模型合理性标准检验思路,以黏土微观力学势能与单元动能之间的转化为出发点,通过对颗粒与颗粒、颗粒与团聚体、颗粒与边框等检验模型赋值不同初始动能,分析系统从发生碰撞到最终平衡前后的能量变化,对开发程序的合理性和科学性进行检验;并且,本文最终通过控制边框位移的方式,实现对宽粒径范围团聚体试样的一维压缩模拟,验证了程序的有效性。综上,本文研究在确保黏土离散元模拟研究的可靠性、准确度及运算效率的同时,拓展离散元在黏土中应用广度和深度,为建立黏土力学性质的宏微观联系奠定方法基础,具有明确和普遍的学术意义。