DDA中细观参数对宏观特性的影响研究.docVIP

DDA中细观参数对宏观特性的影响研究摘 要：采用单轴拉伸试验计算模型，系统研究了DDA细观参数对宏观特性的影响。研究结果表明：（1）计算模型宏观弹性模量只与接触弹簧刚度和材料弹性模量有关，且与两者均呈对数函数关系；（2）抗拉强度与接触弹簧刚度和材料弹性模量呈现对数函数关系，和节理面粘聚力呈线性关系。并给出了具有一定参考意义的相关关系式。关键词：DDA；细观参数；宏观特性；数值模拟1 引言非连续变形分析方法DDA[1，2] （Discontinuous Deformation Analysis）是石根华博士在极限平衡，有限元法以及离散元法基础上发展的一种新的分析块体系统运动和变形的数值方法，自1986年DDA方法问世以来，受到了国内外岩土工程界学者的广泛关注。国内外学者利用非连续变形分析方法在许多工程方面取得了很好的研究成果[3-7]。然而，在模拟工程问题时，首先必须确定与宏观材料特性相匹配的细观参数，这一调试过程将耗费大量的时间。这一点与离散元法类似，国内学者对于离散元中细观参数与宏观特性的关系已经有了很多研究[8-12]。而DDA方法中细观参数与宏观特性的关系，国内还极少有研究。针对这种情况，本文采用单轴拉伸试验计算模型，系统研究了DDA细观参数对宏观特性的影响，并建立了经验表达式，为快速、合理地确定DDA模型细观参数提供一条很好的途径。2 DDA理论简介DDA的理论要点有以下几个方面的内容[13-14]：2.1 块体的位移和变形DDA方法以位移为未知量，每个块体的位移和变形由6个变量表示：其中（）为块体重心（）的刚体位移； 为块体绕点（）的转动角（弧度）；为该块体的正应变，为切应变。块体内任意一点的位移可以表示为：（1）2.2 块体系统的嵌入判别系统中块体间的接触类型有 3 种：凸形角点与棱边之间的接触；凸形角点与凹形角点之间的接触；凸形角点与凸形角点之间的接触。2.3 方程的迭代与求解假设块体系统中有n个块体，根据势能最小建立的平衡方程：（2）式中，是的子矩阵，和是的子矩阵，代表块体的变形变量。3 单个细观参数对宏观特性的影响研究3.1 单个细观参数对弹性模量的影响3.1.1 块体材料弹性模量对弹性模量的影响采用控制变量法，数值模拟过程中保证其他参数不变，只改变块体材料弹性模量。所得结果如图1所示。由图1可知，当块体材料弹性模量小于600Gpa的时候，弹性模量随着块体材料弹性模量的增加而快速增加，当块体材料弹性模量大于600Gpa，小于1500Gpa的时候，弹性模量随着块体弹性模量的增加缓慢增加，当块体材料弹性模量大于1500Gpa的时候，弹性模量几乎保持120Gpa不再变化。两者之间的相关关系可以表述为（3）3.1.2 块体材料抗拉强度对弹性模量的影响块体材料抗拉强度与弹性模量的关系如图2所示，从图中可以看出，随着块体材料抗拉强度的增加，弹性模量报纸不变，表明宏观弹性模量与块体材料抗拉强度无关。3.1.3 块体材料摩擦角对弹性模量的影响块体材料摩擦角与弹性模量的关系如图3所示，从图中可以看出，当块体材料摩擦角在从15度增加到45度的过程中，宏观弹性模量保持不变，表明宏观弹性模量与块体材料摩擦角无关。3.1.4 块体材料粘聚力对弹性模量的影响块体材料粘聚力与弹性模量的关系如图4所示，从图中可以看出，随着块体材料粘聚力的增加，宏观弹性模量保持不变，表明宏观弹性模量与块体材料粘聚力无关。3.1.5 接触弹簧刚度对弹性模量的影响接触弹簧刚度与弹性模量的关系如图5所示，从图中可以看出，当接触弹簧刚度小于200Gpa的时候，弹性模量随着接触弹簧刚度的增加而迅速增加，当接触弹簧刚度大于200Gpa小于400Gpa的时候，弹性模量随着接触弹簧刚度增加而缓慢增加，当接触弹簧刚度大于400Gpa之后，弹性模量将保持120Gpa几乎不再变化。两者的相关关系式可以表示为：（4）3.1.6 节理面抗拉强度对弹性模量的影响节理面抗拉强度与弹性模量的关系如图6所示，可以看出，宏观弹性模量与节理面抗拉强度没有关系。3.1.7 节理面摩擦角对弹性模量的影响节理面摩擦角与弹性模量的关系如图7所示，可以看出，宏观弹性模量与节理面摩擦角没有关系。3.1.8 节理面粘聚力对弹性模量的影响节理面粘聚力与弹性模量的关系如图8所示，可以看出，宏观弹性模量与节理面粘聚力没有关系。3.2 单个细观参数对抗拉强度的影响3.2.1 块体材料弹性模量对抗拉强度的影响块体材料弹性模量与抗拉强度的关系如图9所示，可以看出，当块体材料弹性模量小于500Gpa时，抗拉强度随着块体材料弹性模量的增加而增加，当块体材料弹性模量大于500Gpa时，抗拉强度几乎不再变化，两者的相关关系式可以表示为（5）3.2.2 块体材料抗拉强度对抗拉强度的影响块体材料抗拉强度与宏观抗拉强度的关系如图10所