2离散单元法讲稿的.ppt

离散单元法; 岩体是一种地质材料，它经受长期的地质构造作用，在一定的地质环境中形成一定的结构，显现出宽广和多变的材料响应范围。岩体与一殷工程材黔相比，其最大特点是一般具有结构上的不连续比多为层面和节理面等弱面所切割。实际k，岩体结构的连续性和不连续性那是相对的。如图l—1所示，当取样范图较小时，可能碰到完整岩石，逐渐扩大取样范围，则是节理岩体。; 一般从宏观意义上说，岩体可以视为连续介质，从而可以用弹性力学或塑性力学的方法来进行分析和计算。但在某些情况下，岩体却不能视做连续介质，如地下节理岩体中的巷道(见图1—2)，这时，就不宜用处理连续介质的力学方法来进行计算。于尼离散单元法作为一种处理节理岩体的数值方法就应运而生。; 离散元技术的威力可从下图模拟理想轴向压缩试验中看到。左边的承台推动试件撞击另一个固定的承台．当试件中的应力超过强度极限时，试件出现裂纹。进一次压缩引起大变形，试件断裂，伴随着高速度出现一些飞离的碎片、应当注意断裂过程不用用户干预，断裂方向、新单元的产生及新的相互作用力都足由程序自动计算的。如用常规的有限元或有限差分法进行这种比较直观的离散元分析则是非常困难的。;Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd.;Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd.;Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd.; 工程中所见到的岩体其形态常呈非连续结构，所形成的岩石块体运动和受力情况多是几何或材料非线性问题，离散单元法充分考虑到岩体结构的不连续性，适用于解决节理岩石力学问题。 离散单元法除了用于边坡、采场和巷道的稳定性研究以及颗粒介质微观结构的分析外，已扩展到用于研究地震、爆炸等动力过程和地下水渗流、热传导等物理过程。;离散单元法的基本思想 离散单元法也将区域划分成单元(见图1—3(a))。但是，单元因受节理等不连续面控制，在以后的运动过程中，单元结点可以??离，即一个单元与其邻近单元可以接触，也可以分开(见图1—3(b))。 单元之间相互作用的力可以根据力和位移的关系求出，而个别单元的运动则完全根据该单元所受的不平衡力和不平衡力矩的大小按牛顿运动定律确定。; 离散单元法是一种显式求解的数值方法。该方法与在时域中进行的共他显式计算相似，例如与解抛物线型偏微分方程的显式差分格式相似。 “显式”是针对一个物理系统进行数值计算时所用的代数方程式的性质而言。在用显式法计算时，所有方程式一侧的量都是已知的，而另一测的量只要用简单的代入法就可求根。这与隐式法不同，隐式法必须求解联立方程组。;Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd.;2 刚性块体模型 离散学元法是将所研究的区域划分成一个个分立的多边形块体单元，单元之间可以看成是角—角接触、角—边接触或边—边接触，而且随着单元的平移和转孔允许调整各个单元之间的接触关系。最终，块体单元可能达到平衡状态，也可能一直运动下去。 离题单元法的理论基础是结合不同本构关系的牛顿第二定律，因而可以采用动态松弛法或静态松弛法进行求解。离散单元法的原理比较简单，解决非连续介质大变形问题时非常实用。; 离散单元法的单元，从性质上分可以是刚性的，也可以是非刚性的；从几何形状上分，可以是任意多边形，也可以是圆形。 下面将要介绍的是刚性块体模型，块体可以是任意多边形。刚性假设对于应力水平比较低的问题(如边坡稳定和放矿模拟等)是合理的。;2．2 离散单元法的基本方程 在解决连续介质力学问题时，除了边界条件外，还有3个方程必须满足，即平衡方程、变形协调方程和本构方程。变形协调方程保证介质的变形连续；本构方程即物理方程，它表征介质应力和应变间的物理关系。对于离散单元法而言，由于介质一开始就假定为离散块体的集合(图2—1(a))，故块与块之间没有变形协调的约束，但平衡方程