高等岩石力学讲解.ppt

8.2 非连续变形分析方法（DDA） 特点：在DDA问世之前，岩石力学数值方法中处理节理问题最早是Goodman等人的节理单元，随后是Cundall提出的离散单元法。前者只能模拟岩体中的少量节理，后者用于块体化的岩体，但两种方法都采用了假想力来控制滑动、防止块体重叠等。DDA方法以位移作为未知量，并按有限元法中结构矩阵分析的同样方法求解平衡方程。DDA用于块体系统的力和位移的相互作用，允许每个块体有位移、变形和应变，还允许整个块体系统在块体交界面上有滑动、张开与闭合，但不允许块体互相侵入。 与有限单元法的比较：DDA的建模和求解方法与有限元比较接近，如使总势能最小来建立平衡方程，选择位移作为联立方程的未知量，把刚度、质量及荷载子矩阵加到联立方程的系数矩阵中等。但DDA的块体刚度矩阵比有限元的单元刚度矩阵简单，其最大的优点就是块体边界不一定连续，而在有限元方法中处理不连续问题则十分困难。 与离散元的比较：DDA是位移法，平衡方程中的未知量是位移，而离散单元法是力法，它企图用阻尼调节接触力使其成为常数。另外，DDA是一种隐式方法，其快速、精确及具有正确能量消耗的特点使其易于推广应用。 8.3 数值流形方法（NNM） 简介：岩体是一种非常特殊的介质，其既不是完全的连续体，也不是完全的非连续体，而是介于二者之间的一种部分非连续介质。前述的有限元法等适用于完全连续介质，要处理非连续问题则十分困难，且数量有限；而离散元或DDA等则处理完全非连续的问题，走向了另一个极端，虽然可以处理块体系统的问题，但与工程实际岩体之间仍存在一定的差距。数值流形方法正是基于解决部分非连续介质问题而提出来的，采用“有限覆盖”技术，统一解决连续和非连续变形的力学问题。其基本思想是在求解域上构造一组函数，称为覆盖函数，这种覆盖函数有两个基本性质：（1）局部非零性，即在一个局部区域内不为零；（2）单位分解性，即在求解域内这组覆盖函数之和恒为1。数值流形方法中，对于连续变形区域采用分片可微的覆盖函数，对于非连续变形部分，即块体之间的接触面处，则采用非连续的覆盖函数。因此通过采用连续和非连续覆盖函数的办法可以把连续和非连续的力学问题的计算统一起来。连续体的有限元方法（FEM）和块体系统的非连续变形分析方法（DDA）都是数值流形方法的特例。 关键：覆盖函数的构造。目前常用的是Lagrange型分片插值函数，这也是有限元方法中经常采用的一种插值函数。 8.3 数值流形方法（NNM） 由数学网格和物理网格生成的有限覆盖：所谓“流形”就是把许多个別的重叠的区域连接在一起，去覆盖全部材料体。因此，总体形状可以用局部覆盖所定义的函数来计算。此法有分开的且独立的数学覆盖和物理网格：数学覆盖只定义解的精度，无物理网格，作为实际的材料边界，定义其积分区域。 数学覆盖：用户选择，有许多有限重叠覆盖组成，如有限元网格、规则格子等能转换为数学覆盖。 物理网格：材料边界、裂缝、块体和不同材料区域的交界面，代表材料条件，不能认为选择。 物理覆盖系统：＝数学覆盖＋物理网格。如裂缝或块体边界把一个数学覆盖分成两个或更多个的完全不连续的区域，这些区域定义为物理覆盖。物理覆盖实际上是不连续面对学覆盖的再剖分。 单元：两个或更多的物理覆盖的共同部分。在DDA中，每个块体是一个数学覆盖和物理覆盖，任何两个覆盖之间没有重叠。 简单示例：见图。 9 岩石力学与工程中的反分析方法 概述 反分析的目标未知数 反分析的基本方法 9.1 概述 岩石力学参数问题：岩石力学问题的分析是基于岩石力学参数（包括弹性常数、强度参数、地应力等）进行的，只有参数准确，分析结果才能符合实际。长期以来，岩石力学参数“测不准”一直困扰着岩石力学工作者，从理论和试验手段的各方面想法解决这一问题。但不管是室内试验，还是现场试验，所得的岩石力学参数都只能代表局部，而推广到整体就存在误差。经过理论研究和工程实践积累，工程界和学术界逐步认识到通过现场监测结果反演岩石力学参数，也许是一条可行的途径。 反分析概念：所谓反分析，就是利用工程实际已经取得的量测信息，来反推必要的计算参数。 由于工程实际的量测信息一般为位移，所以反分析也称为位移反分析，也有称为反演分析的。实际上就是岩石力学问题的反问题，将待求量（位移）以实测结果代替，作为已知量，而将必须输入的参数（已知量）作为未知量来求解。然后利用反演所得参数，进一步进行工程其他部分的分析、预测。 反分析应用：反分析最早是从地下工程中发展起来的。由于地下工程位于地层介质中，地下结构的荷载与地层介质的力学参数、地应力等密切相关，而实际工程中为了保证围岩稳定，一般要进行相应量测，获得了大量的信息，所以工程界和学术界逐步认识到这些量测信息的重要性，逐步形成了反分析的理论和方法。随着反分析理论的逐步完善及在地下工程实