吉林大学岩土工程数值法作业最新版10级.docx

《吉林大学岩土工程数值法作业最新版10级.docx》由会员分享，可在线阅读，更多相关《吉林大学岩土工程数值法作业最新版10级.docx（10页珍藏版）》请在优知文库上搜索。

1、岩土工程数值法班级：专业:姓名：学号：吉林高校建立工程学院书目一、问题提出3二、围岩离散化4三、数据打算6四、计算过程6五、结果初步分析7六、图形7七、隧道开挖对围岩的影响9八、结语12九、参考文献13前言随着我国经济快速开展，各种隧道、马路、铁路、房屋以及其它根底设施进入了一个高速建立的阶段，随之而来的是土木工程的跨越式开展。土木工程的设计和探讨手段也有了很大的提高，从以前的基于经验的设计理论渐渐过渡到定量与定性相结合的反分析计算理论。目前为止，土木工程的探讨方法主要有以下五类：类比法；解析法；模型模拟（物理模型方法）；现场监控量测；数值法（数值模拟）。通过岩土工程数值法这门课程，我们系统的

2、学习了数值法中的有限单元法的原理以及它在岩土工程中的应用。随着计算机的普及和运算速度的提升，为弹性力学的数值解法开拓了广袤的领域，尤其在隧道工程中，采纳有限单元法分析都得到了满足的结果。目前，有限单元法已经是解决不同岩体构造、围岩与支护相互作用、隧道围岩压力、围岩应力和变形、围岩破坏过程与破坏机制的主要方法。本文将针对一个隧道开挖实例，应用有限单元法进展位移、应力等相关参数的分析。一、问题提出在岩土体中修建隧道是一件特别困难的工程。因为岩土体是地壳内外力长期作用下形成的一种困难的地质体，具有自然应力、非均质、不连续、各向异性等特点，从而表现在力学性质上具有非线性、剪胀性、蠕变性等。而有限单元法

3、可以将岩土体困难多变的力学性质，根本地质因素、困难和混合的边界条件、岩土体与工程构造物的组合作用等问题统筹考虑，以得到接近实际的数值解答。目前，隧道施工和设计都是基于“新奥法，新奥法的核心是充分发挥围岩的承载实力，将围岩视为承载的主体。随之而来的是如何确定围岩收敛的极限位移，如何确定衬砌的支护时间，如何判定围岩应力的集中程度等问题。本文基于以下条件进展隧道开挖后的围岩进展分析。本隧道断面为曲墙式，拱部半径为5.5m,下部为IInlX6m的矩形，跨度为Ilm,隧道埋深320mo岩体参数为：岩体弹性模量E=3.2xK）3峥r,泊松比4=0.29,岩体初始粘结力Co=L3的如初始摩擦角%=40。，岩

4、体剩余粘聚力C=LlA剩余摩擦角e=30。，岩体容重q=0.027MN/m3,岩体单轴抗压强度=0.9MP具体见图1。二、围岩离散化1）探讨区域大小的确定：对于圆形隧洞，由岩体力学的学问可知，处于6倍洞径以外的岩体不受碉室开挖的影响。本隧道开挖形态不规那么，因而考虑取开挖轮廓线外8.1倍碉径（）围岩为探讨区域（即拱墙交线最深处坐标（50,0）；对于胴室下方矩形局部，那么按碉室相像的比例，相应考虑开挖轮廓线外的44m围岩为探讨区域（胴室底中心线最深书坐标（0,-50）o考虑具体见图2。（2）单元数量的确定：由于本次运用的有限元计算程序须要对结点数进展限制。因而胴室底板划为8段，侧边墙划为8段，拱

5、部圆弧划为8段。在径向划分为10段。所以总共划分275个结点，240个单TIjo图2三、数据打算工程有限元分析的参数如下表所示:参数弹性模量泊松比内摩擦角剩余内摩内聚力剩余内聚容重(MN/抗拉强度(Mpa)()()擦角(MPa)力(Mpa)Irf)(Mpa)数值3200940301.31.17同时由于隧道埋深320m,垂直方向围岩压力pv=h=0.027320=S.6Mpa,水平方向围岩压力Ph=pv=4pv=8.6=3.5MPa,而深部切向应力%=0。-1-0.29a四、计算过程有限单元法分析计算的过程可以简化为5步：1、划分单元，离散构造；2、单元分析，求单元刚度矩阵；3、整体分析，求整体

6、刚度矩阵；4、约束条件引入和力的移置，解构造平衡方程组；5、求单元应力和位移。下面具体分析各个步骤：1、划分单元，离散构造：确定有限元计算范围。2、单元分析，求单元刚度矩阵：首先确定四边形单元的位移模式=cif1+a2x+cc3y+a4xyV=a5-ahx+a1y+asxy,经过整理得到：二网伤心,其中n称为形态矩阵；再由几何方程得到结点应变和结点位移的关系式：=B,其中网为几何矩阵；紧接着由物理方程可以推出应力矩阵=Di=DlB=sMf其中冈为应力矩阵；最终通过虚功原理得出结点位移与结点力的关系K=J闻ISkwy,K为单元刚度矩阵，表示单元反抗变形的实力，它具有对称性、奇异性和分块性的特征。

7、3、整体分析，求整体刚度矩阵：为了求出整个离散化构造的总体方程组，只需对总体构造中的每一个节点建立平衡方程即可。并将结点力用结点位移便是，即可得出离散体含n个结点时的总体方程组为2n2n阶方程，即ZK3=ZHJ,也可以表示成K5=R,其中Kee就是整体刚度矩阵，它是单元刚度矩阵的总和，是由各单元刚度矩阵的元素在相应位置上迭加而得，即所谓“对号入座的方法。4、约束条件引入和力的移置，解构造平衡方程组：由方程组K监=R可知，要求得同须要知道K和R,K为奇异矩阵，通过引入约束条件，划去对应的行和列可以消退它的奇异性。而因可以通过对力的移置转移到结点上。由此得到一个矩阵方程组，通过求解可以得到各个结点

8、的位移。5、求单元应力和位移：求得结点位移后，通过物理方程b=s界可以得到单元的应力。五、结果初步分析通过学院自编有限元软件计算，得到位移改变图、应力矢量图以及相关的数据。通过位移改变图形和相关数据分析得出：隧道开挖后，碉四周岩体产生了向下及向临空面方向的位移。拱顶最大下沉量为29.32mm,底部围岩向上最大变形量为20.05mmo而对于应力矢量，最大主应力出现在洞壁出，在胴室四周岩体径向应力很小，环向应力很大，应力集中程度最大，随着远离胴室，围岩集中系数减小，渐渐接近自然应力状态。六、图形隧道开挖后，围岩位移改变图和应力矢量图如下所示。（完整成图见附图）图3胴室位移改变整体及局部构造图图4围

9、岩应力矢量整体及局部构造图七、隧道开挖对围岩的影响隧道开挖对围岩的影响主要从围岩的位移改变特征及缘由；应力特征及集中程度分析；围岩压力特征和隧道施工建议，这三个方面进展具体分析。1、围岩的位移改变特征及缘由表1胴室周边特别点位移位置碉室底中线碉室底1/4占八、碉室边墙底部边墙距底部1/4占八、墙底1/2占、边距部点拱墙交占拱部1/4点拱部中线占/、y位移(cm)y位移(cm)X位移(cm)X位移(cm)X位移(cm)y位移(cm)y位移(cm)y位移(cm)数值2.0051.6350.2420.3790.5101.0772.3102.932表2洞底中心线y改变结点号123457911洞底y位移

10、(cm)2.0051.4771.0340.69917128200表3拱、墙交点X改变结点1771781793184结点1771781793184号拱墙交点X（cm）0.5540.2490.114-49-0510140.011在胴室开挖后，由于围岩压力的作用，碉四周各点均产生向洞内的收缩位移，具体表现为洞顶产生向下的位移，出现沉降，而且拱顶中线处沉降量较大，并向边墙渐渐减小。同时，拱底面上各个结点均向上抬升，产生竖向位移，总体上呈现从边墙底部向底板中心渐渐增大的趋势。边墙产生向洞内收缩的位移。具体表1所示。同时，为了得到围岩收敛位移随位置的改变，取洞顶和洞底以及拱墙交点作为限制点，比拟它们的特征

11、位移，如表2和表3所示，得到相关规律：随着结点离围岩开挖轮廓线距离渐渐增大，结点向临空面收缩的位移量渐渐减小。隧道开挖后，在洞壁上由于径向的约束解除，围岩向着临空面收缩产生位移。当位移到达围岩的极限位移过后，围岩产生松动坍塌。隧道围岩收敛量测是新奥法施工的重要量测工程，通过测试确定围岩改变趋势，确定何时进展支护。2、围岩应力改变特征及集中程度分析胴室开挖后，FINIT2.DAT文件可以反映出，在洞底有多个单元（1、12、23、34、45、67、78、89）压应力特别小（小于IMPa）。如表4和所示，应力有由负值（受压）变为正值（受拉）的可能性很大，这说明底部会有很大的变形，有破坏的可能性，因此

12、须要适时的施做仰拱，防止衬砌陷落和抵挡底部围岩压力。同时视察，在胴室开挖后，在洞壁上的单元的最大和最小主应力如下表所示:表5洞壁单元最大和最小主应力单元111281913.4330414.61616%0.771150.65056.617681.12713单元1118119113.696092.174351.400951.373111.665391.79001由此，我们发觉，在洞壁围岩上，最大主应力根本是呈现沿着胴室的环向，而最小主应力呈径向，且值特别小。这是由于洞壁围岩开挖使得径向应力解除，因此假如通过喷射混凝土供应切向摩擦力，那么会很大程度上加固围岩，这就是新奥法三大支柱之一（喷锚支护）的原

13、理。对于应力集中，由胴室围岩应力图，我们得出在洞壁上围岩的应力集中程度最大，随着向深部围岩延长，围岩应力集中程度降低，接近围岩自然应力状态。同时视察应力改变状况，在拱墙交织处、拱底和边墙交织处，围岩的应力集中现象最为明显，这些区域由于处于围岩不规那么的地方，简洁产生应力集中。对于拱底和边墙交汇处，可以通过施做仰拱，封闭断面的同时，形成圆弧形过渡区，减小应力集中程度。3、围岩压力特征及施工建议经过以上分析，得出围岩压力有以下特征：嗣室开挖使得径向应力释放，环向应力增大，在靠近洞壁的地方围岩应力集中程度突出，尤其是拱墙交汇处和墙、底板交汇处。在探讨区域中，最大应力也出现在洞壁处，同时在围岩底板旁边压应力变得特别小并极可能转化为拉张应力。在拱顶和拱底向临空面收缩的位移量很大，当超过肯定数值会产生松动坍塌，围岩由形变压力转化为松动压力。针对围岩单轴抗压强度低，围岩变形大的特点，在施工时尽量采纳分部开挖法，如环形开挖留核心土法和双侧壁导坑法，并刚好施做衬砌。爆破时，可以采纳限制爆破的方法，减小对围岩的扰动。对拱顶依据位移收敛量测结果，在适当的时候打入锚杆，施做钢拱架，也可以铺设钢筋网。对于洞底，须要刚好施做仰拱，尽快形