裂隙岩体三维非线性有限元研究

《裂隙岩体三维非线性有限元研究》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、裂隙岩体三维非线性有限元研究摘要：讨论天然裂隙岩体的机械力学行为及力学问题，建立了相应的有限元格式，介绍了裂隙岩体刚度方程的建立及弹塑性状态下的应力一一应变关系是非线性的表达式，阐述了非线性问题的分类及解法。关键词：裂隙岩体；接触单元；非线性问题工程岩体经受过复杂的地质作用后，其内部存在不同规模、性质各异的地质结构面（节理、裂隙、断层等），使得岩体的结构特征与力学性质错综复杂。一方面，不连续面的存在使得岩体介质类型不同于其他任何一种力学材料。在进行裂隙岩体力学问题分析时，准确选择岩体力学模型成为必须解决的首要问题。另一方面，岩体中的不连续面使得岩体性质表现

2、出不连续、非均匀、各向异性及尺寸效应等特征。这就使得传统弹性有限元法不能更好地模拟真实情况。近年来，节理单元被广泛应用于模拟岩体中的断层、剪切带和不连续区域。在节理单元中，由于引入了岩石燕子节理的法向和切向的刚度，使得节理单元的应用受到了一定的限制。这是因为刚度系数通常都是从野外实测得到的，测试存在困难和不确定性。Katona于1983年提出了一种不用刚度系数的简单接触摩擦单元，它能模拟两物体之间的滑动摩擦、张开和闭合过程，适用于不计厚度影响的节理、断层和位移不连续场的有限元分析。由于这种单元采用的是常接触力二节点简单单元，因而很难适应接触面复杂的问题。同

3、时,由于该单元选取节点接触力作为基本未知量，接触应力是由接触力平均得到的，这不仅降低了计算接触应力的精度，而且使得这种简单求应力的方法很难推广至三维问题。雷晓燕于1994年在此基础上提出了一种新的接触摩擦单元。在这种单元中，直接选取节点接触应力作为基本未知量，同时采用六节点的等参单元模拟接触面复杂的几何形状。接触单元的几何和静力约束是作为补充方程包含在刚度方程中的。运用虚位移原理导出接触摩擦问题的等效单元刚度约束矩阵和等效载荷向量。1、非线性问题的分类引起结构产生非线性变形的原因很多，大致可以分为三类：材料非线性问题、几何非线性问题、状态非线性问题。1）材

4、料的非线性问题：指材料的物理定律即其应力应变关系是非线性的，但只考虑小位移和小应变的情况，也就是指结构整个几何形状的变化及结构材料内部的应变与结构尺寸相比，是无限小的，这样可以忽略微元体的局部应变。例如，计算应力时可以采用原来的、未变形的微元体面积。至于应变一一位移关系则采用线性的应变一一位移关系式。各种小应变、小位移问题的结构弹塑性分析及岩土工程中的问题一般都属于这一类。2）几何非线性问题：虽然假设线性的应力一一应变关系，但非线性的应变一一位移关系和几何形状的有限变化却引起几何非线性问题。大多数的几何非线性问题是小应变和大位移，当然，也包括大位移和大应变

5、的情况。3）状态非线性问题：例如，一根只能拉伸的电缆可能是松弛的，也可能是张紧的；轴承套可能是接触的，也可能是非接触的；冻土可能是冻结的，也可能是非冻结的；岩土体中的断层结构面间相对运动。由于系统状态的改变，从而引起系统的刚度在不同的值之间突然变化。状态改变也许和载荷直接有关（如电缆情况），也可能由某种外部原因引起（如冻土情况）。岩土体材料的一个重要特征是其应力一一应变关系具有明显的非线性性质。因此，对于岩土体通常都按材料非线性问题来考虑。材料非线性问题又可以分为两种情况：一种是非线性弹性问题；另一种就是非线性弹塑性问题，系材料超过屈服极限以后就呈现出非线

6、性性质，各种结构和岩土介质若采用弹塑性的本构模型进行分析就是这类问题。在加载过程中，这两种非线性问题在本质上相同，但是卸载过程就会出现不同的现象，非线性弹性问题是可逆过程，卸载后结构或介质会恢复到加载前的位置，非线性弹塑性问题是不可逆的，它将会出现残余变形。对于材料非线性问题，由于应力一一应变关系是非线性的，对于这类问题的微分方程的求解在数学上有一定的困难。然而有限元法处理非线性问题却是十分有效的。用有限元法处理非线性问题的基本思想是用一系列线性问题的解来逐步逼近非线性问题的解。因而，非线性问题可以理解为一系列线性解进行迭代过程的结果。2材料非线性问题的解

7、法用有限元解材料非线性问题，通常采用三种基本的方法来解析，即增量法、迭代法以及增量迭代法(或混合法)。在说明问题时，可以只考虑单个单元的非线性平衡方程：[k(6)]e{8}={R}(1)这是是非线性出现在单元刚度矩阵[k]e中，k[8]e是非线性材料性质[D(e)]的函数。可以把[k]e写成[k]e={R}]的形式，因而[k]e中的材料参数不再是常数。图1中的载荷{R}与节点位移{5}相对的非线性应力一一应变关系。根据这个应力一一应变关系或本构关系，进而确定非线性分析时变化的矩阵[D(£)]。本文主要介绍迭代法。图1非线性的载荷一一位移曲线图2.1常刚度迭

8、代法对于方程(1)：：k(8)]e{8}={R},假设R二RO已知