三维拉格朗日法计算原理..docVIP

1 三维快速拉格朗日法的基本原理 1.1 概述 目前在岩土力学中常用的数值计算方法有差分方法、有限元法、边界元法等几种，特别是后两种方法，随着计算机的发展其应用尤为广泛。但是，这几种方法都是以连续介质为出发点，而且往往囿于小变形的假定。它们虽然也可以用来解决由几种介质所组成的非均质的问题，并且对于个别的断层或弱面，也可以用设置节理单元的办法来解决，但是用以解决富含节理和大变形的岩土力学问题，往往所得的结果与实际的物理图景相差甚远。于是离散单元法和拉格朗日元法就应运而生。 离散单元法是Cundall于上世纪70年代初所提出的。该法将为弱面所切割的岩体视为复杂的块体的集合体，允许各个块体可以平移或转动，甚至相互分离。拉格朗日元法则是由Cundall所加盟的美国ITASCA咨询集团于1986年所开发的。该法将流体力学中跟踪流体运动的拉格朗日方法应用于解决岩体力学的问题获得成功。 三维快速拉格朗日法是一种基于三维显式有限差分法的数值分析方法，它可以模拟岩土或其他材料的三维力学行为。三维快速拉格朗日分析将计算区域划分为若干四面体单元，每个单元在给定的边界条件下遵循指定的线性或非线性本构关系，如果单元应力使得材料屈服或产生塑性流动，则单元网格可以随着材料的变形而变形，这就是所谓的拉格朗日算法，这种算法非常适合于模拟大变形问题。三维快速拉格朗日分析采用了显式有限差分格式来求解场的控制微分方程