DDA块体接触开r闭迭代研究.PDFVIP

　 第３１卷第４期 长　江　科　学　院　院　报 Ｖｏｌ．３１　Ｎｏ．４　　 　　２０１４年４月 ＪｏｕｒｎａｌｏｆＹａｎｇｔｚｅＲｉｖｅｒＳｃｉｅｎｔｉｆｉｃＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅ Ａｐｒ．２０１４　 　 　　ＤＯＩ：１０．３９６９／ｊ．ｉｓｓｎ．１００１－５４８５．２０１４．０４．０１２ ２０１４，３１（０４）：５２－５８ ＤＤＡ块体接触开 －闭迭代研究 冯细霞，邬爱清，卢　波 （长江科学院水利部岩土力学与工程重点实验室，武汉　４３００１０） 摘　要：非连续变形方法是一种隐式求解的动力学计算方法，且采用在块体界面加减刚硬弹簧的方式来满足块体 界面无张拉和无嵌入的接触准则。基于建立的ＤＤＡ块体接触分析模型，推导了静力条件下块体的接触位移理论 公式，并且通过接触力的计算，验证了该公式的正确性。运用此公式，定量研究了块体变形模量、弹簧刚度、时间步 长和接触位移之间的相互关系，并绘出了当块体变形模量和弹簧刚度一定时，接触位移随时间步长变化的关系曲 线。通过建立简单的块体系统力学分析模型，剖析ＤＤＡ中开－闭迭代完整的实现过程，证明了石根华博士在ＤＤＡ 计算中采用时间步长来调整和加快开－闭迭代收敛的策略是合理、高效的。该研究也揭示了ＤＤＡ各参数取值应 依据的基本规律，可以推动ＤＤＡ的实际工程应用。 关键词：ＤＤＡ；接触力；时间步；弹簧刚度；变形模量；接触位移 中图分类号：ＴＵ４５　　　文献标志码：Ａ　　　文章编号：１００１－５４８５（２０１４）０４－００５２－０７ 本文块体接触分析模型中，仅考虑２个规则的矩形 １　研究背景 块体接触，下部块体固定，上部块体受均匀分布的铅 直向载荷作用。因此，在该模型中，可只考虑法向嵌 石根华博士提出的非连续变形方法（ＤＤＡ），是 入这种情况，这样仅涉及到法向刚度弹簧，整个计算 继关键块体理论后研究块体系统的一种新的计算方 过程可以手工完成，而且所得到的结果是解析的。 ［１－２］ 法 。该方法通过满足块体系统力系的平衡、 建立了位移解析计算式后，便可以从量化的角度去 开－闭迭代及动力求解收敛，来获得块体系统力与变 研究时间步长、块体弹性模量以及弹簧刚度等控制 形真解。其自提出以来，一直作为岩土工程数值模拟 性参数之间的关系及其对计算结果的影响规律。根 方法研究中的前沿问题，为国内外所重视［３－１１］。 据本文的解析，证明了石根华博士在ＤＤＡ计算中采 非连续变形分析方法是一种隐式求解的动力学 用时间步长来调整和加快开－闭迭代收敛的策略是 计算方法，其中块体接触理论是 ＤＤＡ的核心问题。 合理、高效的。该研究也揭示ＤＤＡ各参数取值应依 由于块体系统的接触是一个高度非线性问题，涉及参 据的基本规律，可以推动ＤＤＡ的实际工程应用。 数众多，如接触面上的弹簧刚度、块体的弹性模量、时 间步长等。关于块体系统接触以及合适参数的取值 等问题的研究，长期以来，主要借助ＤＤＡ程序运行和 ２　ＤＤＡ块体接触开－闭迭代分析模型 具体问题分析，通过参数组合和试算等方式来获取一 为定量研究ＤＤＡ块体接触开－闭迭代的本质 些经验性认识。在某种条件下，建立一个解析的理论 问题，这里建立了如图１所示的２块体接触开－闭 模型，来反映ＤＤＡ计算过程中块体接触位移与上述