《岩土力学数值方法》结课论文.docVIP

岩土力学数值方法的应用及发展前景 班级：硕研14级4班 姓名：吴纪达 学号：140431055 专业方向：结构工程 1 引 言 随着计算机技术进入各个学科领域，近年来岩土力学数值分析方法有了迅速的发展，越来越多的生产、设计和科研部门采用数值方法对岩土工程进行分析，出现了一系列的数值方法，这些数值方法可大致分为连续变形分析方法和非连续变形分析方法两大类。本文简要讨论了目前已较为成熟应用于岩土工程中的多种数值方法，介绍了最新尝试应用于岩土工程中几种新的数值方法，并在此基础上对今后岩土工程数值方法的发展作了一些展望。分析了采用连续介质力学分析岩土工程问题的关键，并讨论分析了岩土本构理论发展现状,提出对岩土本构理论发展方向的思考,最后对数值分析在岩土工程分析中的地位作了分析。 分析表明，岩土工程数值分析结果是岩土工程师在岩土工程分析过程中进行综合判断的重要依据之一；采用连续介质力学模型求解岩土工程问题的关键是如何建立岩土的工程实用本构方程；建立多个工程实用本构方程结合积累大量工程经验才能促使数值方法在岩土工程中由用于定性分析转变到定量分析。 2 岩土工程中当前较为常用的几种数值方法 2. 1 有限单元法 有限元法自20世纪70年代提出发展至今，已经相当成熟，是目前最广泛使用的数值方法之一。有限元分析可以处理复杂的边界条件以及材料的非均匀性和各向异性，还可以有效地模拟材料的非线性应力应变关系，得到围岩的整体应力场、位移场，而得出围岩的可能破坏部位，可对围岩的力学特性有全面了解。但如何从有限元分析的结果得出合适的衡量围岩稳定程度的定量指标还需进一步研究。有限单元法也遇到了本身不可克服的困难，如在前处理困难，计算应力应变解答不连续，难以进行任意路径开裂计算等。 2. 2 界面单元有限元法 为了推广有限单元法,使其能够处理简单的非连续问题,人们提出了各种能够反映非连续变形性质的简单力学模型和特殊界面单元,并将其应用于单一性非连续界面力学行为的模拟中。界面单元有限元法将岩土介质视为准连续介质,仍以连续分析为主,但可以对个别具有控制作用的宏观非连续面的变形与破坏力学行为进行重点分析,在工程中得到了广泛应用。同时,该方法也存在诸多的局限性,如只能对原生的非连续界面进行计算,对次生的非连续界面无法进行处理、界面单元数目不能设置太多、界面弹簧刚度的选取较为困难等。因此,该方法在处理复杂的非连续变形问题时仍显得无能为力。 2. 3 　离散元法（DEM）、块体理论及刚体有限元法 这是分析多节理岩体的另一类重要方法 。在这类方法中通常把岩体当成是由节理切割的岩块组成的结构体 。离散元法 、块体理论 、刚体有限元法是三种有代表性的不连续块体方法 。 离散单元法是Cundall于1971年提出来的,该方法与在时域中进行的其他显示计算相似例如与解抛物线型偏微分方程的显示差分格式相似。离散单元法也像有限单元法那样将区域划分成单元。但是单元因受节理等不连续面控制在以后的运动过程中单元节点可以分离即一个单元与其邻近单元可以接触也可以分开适合节理岩体的应力分析。单元之间相互作用的力可以根据力和位移的关系求出而个别单元的运动则完全根据单元所受的不平衡力和不平衡力矩的大小按牛顿运动定律确定。 其基本思想是岩块之间的相互作用同时受表征位移—力的物理方程和反应力—加速度（位移、速度）的运动方程支配通过迭代求解显示岩体的动态破坏过程。 块体理论和刚体有限元法则属于静态分析的方法。仅考虑在恒定荷载作用下系统的平衡及失稳条件。块体理论实质上是一种几何学的方法，根据岩体中实际存在的节理倾角及其方位，利用块体之间的相互作用条件找出具有移动可能性的岩体及其位置，故也称关键块理论。该理论是一种通过判别和描述洞室围岩最危险岩石块体运动来确定岩体稳定性的分析方法，是由石根华提出并在美与 Goodman 合作完善起来的。认为岩体完全被交错节理所切割，形成各种形状的空间块体，而在这种岩体中开挖的洞室，从方向和位置不利且以开挖临空面为一个或若干面的块体的运动开始的，这个块体称为关键块体。对这种理论的一个新发展是边边接触块体理论，它除了与关键块体理论相同的假定之外还假定块体之间接触可以沿节理全长或部分长度分布。 刚体有限元法则把岩块视为刚体块单元，块体之间以具有相应刚度的弹簧（阻尼器）连结形成总的系统，建立相应的基本方程。用刚体有限元法 ( RBFEM) 给出的应力精度不会低于甚至高于位移精度，优于传统的以位移法为基础的有限元数值解 。此外 RBFEM 大大地减小了总体的半带宽及体积，大幅度降低了计算量。这些优点有利于结构分析，特别是非线性问题的数值解。这类模型如有刚体-弹簧模型，刚性有限元，分块刚体位移-界面应力元法,块体-弹簧模型。在这类模型中，界面特性均假设服从 Coulomb摩擦定律，考虑