离散元方法在混合料的研究中应用.pdfVIP

离散元方法在混合料研究中的应用 成晟，倪富健 （东南大学交通学院，江苏 南京 210096) 摘 要：为考察级配、胶结料等微观属性对混合料宏观性质的影响，引入了离散元处理方法。首先 介绍离散元的基本原理和实现方法，并以二维球体为单元，颗粒接触简化为三部分：法向变形、剪 切变形、滑动或开裂，假设系统始终处于准静态，基于平衡方程开发了离散元计算程序。然后利用 离散元程序进行试件的成型和无侧限抗压强度模拟试验，对不同级配和粘结力下试件的宏观表现进 行比较,并记录试件的破坏过程。研究发现：两种级配的混合料空隙率没有显著区别，水泥稳定材料 的模拟破坏过程表现出明显的非线性和脆性特征，和实际情况较为吻合。 关键词：路面材料，离散元，模拟试验 中图分类号：U416.214 文献标识码：A 1. 引言 如今对路面材料研究也已进入微观领域，众所周知级配、胶结料性质等对混合料的影响很大， 以水泥稳定碎石为例，为抑制裂缝，半刚性基层多采用粗级配和低水泥剂量，通常认为强度与裂缝 密度有关系。但现在对这些因素的研究主要集中在对试验数据的分析和整理，缺少力学方面的解答， 而连续介质力学面对这些问题时又无能为力。离散元善于处理颗粒的接触问题，用来模拟混合料不 失为一种途径。通过调整颗粒和接触性质，进行模拟试验，就能在微观和宏观之间建立联系，为进 一步理解级配、胶结料等对混合料性能的影响提供了工具基础。 2. 离散元方法 2.1 基本原理 离散元方法和有限元方法类似，它每一个单元代表一个颗粒。离散元与普通有限元方法的主要 区别是：变形的同时，伴随着颗粒之间的接触的消失和新的接触的产生，从而导致结构刚度发生变 化，因此，结构的刚度矩阵在计算过程中需要不断更新。刚度矩阵需要更新的另一个原因是：无论 法向还是切向，颗粒间的接触对刚度都是非线性的，因此需要不断地重新计算。 2.2 接触性质 混合料中涉及两种接触：颗粒直接接触和颗粒与粘结料的接触，需要分别处理。 2.2.1 颗粒直接接触 [1] 球体间法向接触的力－位移关系采用Hertz 方案（1881） ，微观模型及参数如图1： 图1 微观模型 - 1 - 法向力Fn 和法向位移n 的关系为： 3 F M (n )2 （2-1 ） n 2 2r G 2r r 其中： M µ ；n r =+r −d ；r i j i j 3(1−ν ) r +r µ i j 接触对刚度 M 依赖于颗粒剪切模量G 和泊松比ν 以及两个颗粒的平均直径r ，出现非线性关 µ µ 系是因为接触面的面积和