离散单元法讲稿.pptVIP

由于本章中所讲的单元是刚性的，因此单元之间的相对位移增量等完全由单元的几何尺寸、质心平移和单元绕其质心转角的大小来决定。如图2—5所示为在一个循环中所历经的计算过程。 离散元方法一变形体模型 4‘1变形体动力学 在对变形体单元进行数学推导之前，让我们先来熟悉要研究的问题。 考察受有常力F作用的单个变形体单元(如图4．1),给定单元的几何、密度和弹性模量，需要确定它的运动特性。我们知道用牛顿定律可以确定质心的运动和绕质心的转动。但不能确定单元相对于质心的变形。 在变形体单元中，单元质心是在运动的，固定在单元质心处的一组坐标既有平动又有转动。为了得到相对于这样一组坐标系的运动变形方程．有必要将通常惯性坐标系中的动力平衡方程转换成用单元主轴表示的形式，单元主轴既作平动又作转动，如图4．2。 建立于这种坐标系中的单元具有零相对线动量和角动量，其变形可精确地用单元固有弹性模态的展开式来表示。模态法完全可用于大转动和大应变问题。对后一种情况,我们必须严格选择适当的应力应变度量方法，并注意变形前后的旋转路径。为避免复杂，我们在这里仅限于小应变问题。 动力平衡方程向非惯性坐标系的转换 以下将讨论在惯性系中作任意大转动的弹性体所须遵循的方程。变形体单元的运动可分解为平均刚体运动和相对运动，刚体运动定义了一个非惯性动力参考系，相对运动指相对于该非惯性参考系的变形。 选择满足零相对线动量和角动员的Tisserand条件的动力参考系。对相对位移的均方值而不是对相对动能取极小值，将得到线性形式的零角动量方程． 反之将得到非线性方程。这就使得相对位移能精确地用物体固有弹性模态的展开式来表示。 运用哈米顿原理及广义坐标的概念，可以推导出所有的刚体运动方程和变形方程，在每一时步中采用小应变的假设，则对简化这些方程。这些方程不仅适用于小变形，同时也能正确地处理大变形的问题。 * * 离散单元法 曹 平 岩体是一种地质材料，它经受长期的地质构造作用，在一定的地质环境中形成一定的结构，显现出宽广和多变的材料响应范围。岩体与一殷工程材黔相比，其最大特点是一般具有结构上的不连续比多为层面和节理面等弱面所切割。实际k，岩体结构的连续性和不连续性那是相对的。如图l—1所示，当取样范图较小时，可能碰到完整岩石，逐渐扩大取样范围，则是节理岩体。 一般从宏观意义上说，岩体可以视为连续介质，从而可以用弹性力学或塑性力学的方法来进行分析和计算。但在某些情况下，岩体却不能视做连续介质，如地下节理岩体中的巷道(见图1—2)，这时，就不宜用处理连续介质的力学方法来进行计算。于尼离散单元法作为一种处理节理岩体的数值方法就应运而生。 离散元技术的威力可从下图模拟理想轴向压缩试验中看到。左边的承台推动试件撞击另一个固定的承台．当试件中的应力超过强度极限时，试件出现裂纹。进一次压缩引起大变形，试件断裂，伴随着高速度出现一些飞离的碎片、应当注意断裂过程不用用户干预，断裂方向、新单元的产生及新的相互作用力都足由程序自动计算的。如用常规的有限元或有限差分法进行这种比较直观的离散元分析则是非常困难的。 工程中所见到的岩体其形态常呈非连续结构，所形成的岩石块体运动和受力情况多是几何或材料非线性问题，离散单元法充分考虑到岩体结构的不连续性，适用于解决节理岩石力学问题。 离散单元法除了用于边坡、采场和巷道的稳定性研究以及颗粒介质微观结构的分析外，已扩展到用于研究地震、爆炸等动力过程和地下水渗流、热传导等物理过程。 离散单元法的基本思想 离散单元法也将区域划分成单元(见图1—3(a))。但是，单元因受节理等不连续面控制，在以后的运动过程中，单元结点可以分离，即一个单元与其邻近单元可以接触，也可以分开(见图1—3(b))。 单元之间相互作用的力可以根据力和位移的关系求出，而个别单元的运动则完全根据该单元所受的不平衡力和不平衡力矩的大小按牛顿运动定律确定。 离散单元法是一种显式求解的数值方法。该方法与在时域中进行的共他显式计算相似，例如与解抛物线型偏微分方程的显式差分格式相似。 “显式”是针对一个物理系统进行数值计算时所用的代数方程式的性质而言。在用显式法计算时，所有方程式一侧的量都是已知的，而另一测的量只要用简单的代入法就可求根。这与隐式法不同，隐式法必须求解联立方程组。 2 刚性块体模型 离散学元法是将所研究的区域划分成一个个分立的多边形块体单元，单元之间可以看成是角—角接触、角—边接触或边—边接触，而且随着单元的平移和转孔允许调整各个单元之间的接触关系。最终，块体单元可能达到平衡状态，也可能一直运动下去。 离题单元法的理论基