裂隙岩体的三维非线性有限元分析.docx

裂隙岩体的三维非线性有限元分析 摘要：探讨自然裂隙岩体的机械力学行为及力学问题，建立了相应的有限元格式，介绍了裂隙岩体刚度方程的建立及弹塑性状态下的应力应变关系是非线性的表达式，阐述了非线性问题的分类及解法。 关键词：裂隙岩体；接触单元；非线性问题 工程岩体经受过困难的地质作用后，其内部存在不同规模、性质各异的地质结构面，使得岩体的结构特征与力学性质错综困难。一方面，不连续面的存在使得岩体介质类型不同于其他任何一种力学材料。在进行裂隙岩体力学问题分析时，精确选择岩体力学模型成为必需解决的首要问题。另一方面，岩体中的不连续面使得岩体性质表现出不连续、非匀称、各向异性及尺寸效应等特征。这就使得传统弹性有限元法不能更好地模拟真实状况。 近年来，节理单元被广泛应用于模拟岩体中的断层、剪切带和不连续区域。在节理单元中，由于引入了岩石燕子节理的法向和切向的刚度，使得节理单元的应用受到了肯定的限制。这是因为刚度系数通常都是从野外实测得到的，测试存在困难和不确定性。Katona于11013年提出了一种不用刚度系数的简洁接触摩擦单元，它能模拟两物体之间的滑动摩擦、张开和闭合过程，适用于不计厚度影响的节理、断层和位移不连续场的有限元分析。由于这种单元采纳的是常接触力二节点简洁单元，因而很难适应接触面困难的问题。同时，由于该单元选取节点接触力作为基本未知量，接触应力是由接触力平均得到的，这不仅降低了计算接触应力的精度，而且使得这种简洁求应力的方法很难推广至三维问题。雷晓燕于11014年在此基础上提出了一种新的接触摩擦单元。在这种单元中，干脆选取节点接触应力作为基本未知量，同时采纳六节点的等参单元模拟接触面困难的几何形态。接触单元的几何和静力约束是作为补充方程包含在刚度方程中的。运用虚位移原理导出接触摩擦问题的等效单元刚度约束矩阵和等效载荷向量。 1、非线性问题的分类 引起结构产生非线性变形的缘由许多，大致可以分为三类：材料非线性问题、几何非线性问题、状态非线性问题。 1）材料的非线性问题：指材料的物理定律即其应力应变关系是非线性的，但只考虑小位移和小应变的状况，也就是指结构整个几何形态的改变及结构材料内部的应变与结构尺寸相比，是无限小的，这样可以忽视微元体的局部应变。例如，计算应力时可以采纳原来的、未变形的微元风光积。至于应变位移关系则采纳线性的应变位移关系式。各种小应变、小位移问题的结构弹塑性分析及岩土工程中的问题一般都属于这一类。 2）几何非线性问题：虽然假设线性的应力应变关系，但非线性的应变位移关系和几何形态的有限改变却引起几何非线性问题。大多数的几何非线性问题是小应变和大位移，当然，也包括大位移和大应变的状况。 3）状态非线性问题：例如，一根只能拉伸的电缆可能是松弛的，也可能是张紧的；轴承套可能是接触的，也可能是非接触的；冻土可能是冻结的，也可能是非冻结的；岩土体中的断层结构面间相对运动。由于系统状态的变更，从而引起系统的刚度在不同的值之间突然改变。状态变更或许和载荷干脆有关，也可能由某种外部缘由引起。 岩土体材料的一个重要特征是其应力应变关系具有明显的非线性性质。因此，对于岩土体通常都按材料非线性问题来考虑。材料非线性问题又可以分为两种状况：一种是非线性弹性问题；另一种就是非线性弹塑性问题，系材料超过屈服极限以后就呈现出非线性性质，各种结构和岩土介质若采纳弹塑性的本构模型进行分析就是这类问题。在加载过程中，这两种非线性问题在本质上相同，但是卸载过程就会出现不同的现象，非线性弹性问题是可逆过程，卸载后结构或介质会复原到加载前的位置，非线性弹塑性问题是不行逆的，它将会出现残余变形。 对于材料非线性问题，由于应力应变关系是非线性的，对于这类问题的微分方程的求解在数学上有肯定的困难。然而有限元法处理非线性问题却是非常有效的。用有限元法处理非线性问题的基本思想是用一系列线性问题的解来逐步靠近非线性问题的解。因而，非线性问题可以理解为一系列线性解进行迭代过程的结果。 2材料非线性问题的解法 用有限元解材料非线性问题，通常采纳三种基本的方法来解析，即增量法、迭代法以及增量迭代法。在说明问题时，可以只考虑单个单元的非线性平衡方程： k e=R 这是是非线性出现在单元刚度矩阵ke 中，k e是非线性材料性质D 的函数。可以把ke 写成ke =k，R的形式，因而ke 中的材料参数不再是常数。图1中的载荷R与节点位移相对的非线性应力应变关系。依据这个应力应变关系或本构关系，进而确定非线性分析时改变的矩阵D 。本文主要介绍迭代法。 图1 非线性的载荷位移曲线图 2.1常刚度迭代法 对于方程：k e=R，假设R= R0已知，则有限元解法的实质就是在R-曲线与R= R0直线的交点的横坐标。为此可用摸索法，假如已经找到一个摸索值n，但不满意式，如图1所示，它在曲线上的纵坐标值Rn低于R0值，用Pn表示。现在可以通过点Pn绘一条直线p，其斜率为kn ，则 kn = R0 - Rn 这样，直线p与R0 - Rn的交点的横坐标n+1就是一个改进了的摸索解。假如按n+1 由式计算得的Rn+1仍小于R0，就可以从 Pn+1动身，重复前面过程，双可以得到进一步的改进摸索解R n+2。这样，循环迭代，直到前后两次的位移值接近为止。 所谓常刚度，是指取kn=k0为某个常数，即在迭代时直线p的斜率不变，故称为常刚度法。 图2 常刚度迭代法图 2.2割线变刚度迭代法 用一系列的割线，即直线p与R= R0直线的交点P1S，P2S，P3S，取逐次靠近R-曲线与R= R0直线的交点PS时，这些直线p的斜率并不肯定相同。如对第n条直线，其斜率为kn，则迭代的方程为 R0= Rn+Kn 因而，可得第n+1次近似解为 n+1=n+ kn 这里的割线刚度kn=kns= Rn /n。其中，kns为从原点动身到所迭代时与曲线交点割母的斜率。因而这个过程称为割线迭代法。 图3 割线变刚度迭代法图 2.3 切线变刚度迭代法 变切线刚度的表达式可表示为 kn=kn= 这里kn在物理上代表切线刚度，因而这个过程称为切线变刚度迭代法。 图4 切线变刚度迭代法图 2.4 New-Raphson 法 对应于非线性方程式，可以考虑一个非线性的初始摸索解0，它与真实解*之差为一个小量，任何一个具有一阶导数的连续函数，在0旁边光滑，则可按一阶泰勒级数绽开为 = + += 0 如仅取其线性项，则 = + = 0 所以，=0- 这个迭代过程，直至在某个允许精度内 <。这种迭代方法就称为New-Raphson 法，简称N-R法。 3结语 节理、裂隙、断层及懦弱夹层是地质中常见的现象，并且它们都对岩体的强度和稳定性有极大的影响。通常，对于工程岩体而言，有以下两种处理方法：把节理岩体视作不连续介质，由岩石结构体和结构面组成，因此分别探讨岩石和结构面的力学性质及岩石和结构面共同作用时的耦合原理，这种方法的代表有离散无法、块体理论和刚弹性法；把节理岩体视作宏观上的连续体，建立岩体的等效本构关系，这种方法的代表有当量体法、断裂力学法和损伤力学法。本文依据岩体中含有节理、层理、断层及懦弱夹层等不连续结构面，探讨了塑性有限元分析中常采纳的迭代法的具体过程。 第7页 共7页第 7 页 共 7 页第 7 页 共 7 页第 7 页 共 7 页第 7 页 共 7 页第 7 页 共 7 页第 7 页 共 7 页第 7 页 共 7 页第 7 页 共 7 页第 7 页 共 7 页第 7 页 共 7 页