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

岩土工程数值法 班 级 ： 专 业 ： 姓 名 ： 学 号 ： 吉林大学建设工程学院 目录一、问题提出 3二、围岩离散化 4三、数据准备 6四、计算过程 6五、结果初步分析 7六、图形 7七、隧道开挖对围岩的影响 9八、结语 12九、参考文献 13前 言随着我国经济快速发展，各种隧道、公路、铁路、房屋以及其它基础设施进入了一个高速建设的阶段，随之而来的是土木工程的跨越式发展。土木工程的设计和研究手段也有了很大的提高，从以前的基于经验的设计理论逐渐过渡到定量与定性相结合的反分析计算理论。目前为止，土木工程的研究方法主要有以下五类：类比法；解析法；模型模拟（物理模型方法）；现场监控量测；数值法（数值模拟）。通过岩土工程数值法这门课程，我们系统的学习了数值法中的有限单元法的原理以及它在岩土工程中的应用。随着计算机的普及和运算速度的提升，为弹性力学的数值解法开辟了广阔的领域，尤其在隧道工程中，采用有限单元法分析都得到了满意的结果。目前，有限单元法已经是解决不同岩体结构、围岩与支护相互作用、隧道围岩压力、围岩应力和变形、围岩破坏过程与破坏机制的主要方法。本文将针对一个隧道开挖实例，应用有限单元法进行位移、应力等相关参数的分析。1、 问题提出在岩土体中修建隧道是一件十分复杂的工程。因为岩土体是地壳内外力长期作用下形成的一种复杂的地质体，具有天然应力、非均质、不连续、各向异性等特点，从而表现在力学性质上具有非线性、剪胀性、蠕变性等。而有限单元法可以将岩土体复杂多变的力学性质，基本地质因素、复杂和混合的边界条件、岩土体与工程结构物的组合作用等问题统筹考虑，以得到接近实际的数值解答。目前，隧道施工和设计都是基于“新奥法”，新奥法的核心是充分发挥围岩的承载能力，将围岩视为承载的主体。随之而来的是如何确定围岩收敛的极限位移，如何确定衬砌的支护时间，如何判定围岩应力的集中程度等问题。本文基于以下条件进行隧道开挖后的围岩进行分析。本隧道断面为曲墙式，拱部半径为5.5m，下部为11m6m的矩形，跨度为11m，隧道埋深320m。已知岩体参数为：岩体弹性模量，泊松比，岩体初始粘结力，初始摩擦角，岩体残余粘聚力，残余摩擦角，岩体容重，岩体单轴抗压强度。具体见图1。2、 围岩离散化1）研究区域大小的确定：对于圆形隧洞，由岩体力学的知识可知，处于6倍洞径以外的岩体不受硐室开挖的影响。本隧道开挖形状不规则，因而考虑取开挖轮廓线外8.1倍硐径（45.5m）围岩为研究区域（即拱墙交线最深处坐标（50,0）；对于硐室下方矩形部分，则按硐室相似的比例，相应考虑开挖轮廓线外的44m围岩为研究区域（硐室底中心线最深书坐标（0，-50）。考虑具体见图2。 （2）单元数量的确定：由于本次使用的有限元计算程序需要对结点数进行控制。因而硐室底板划为8段，侧边墙划为8段，拱部圆弧划为8段。在径向划分为10段。所以总共划分275个结点，240个单元。图 23、 数据准备 工程有限元分析的参数如下表所示:参数弹性模量（Mpa）泊松比（）内摩擦角（）残余内摩擦角内聚力（Mpa)残余内聚力（Mpa)容重（MN/)抗拉强度（Mpa)数值32000.2940301.31.10.0270.9 同时由于隧道埋深320m，垂直方向围岩压力，水平方向围岩压力，而深部切向应力。4、 计算过程 有限单元法分析计算的过程可以简化为5步：1、划分单元，离散结构；2、单元分析，求单元刚度矩阵；3、整体分析，求整体刚度矩阵；4、约束条件引入和力的移置，解结构平衡方程组；5、求单元应力和位移。下面具体分析各个步骤： 1、划分单元，离散结构：确定有限元计算范围。 2、单元分析，求单元刚度矩阵：首先确定四边形单元的位移模式，经过整理得到：，其中称为形态矩阵；再由几何方程得到结点应变和结点位移的关系式：，其中为几何矩阵；紧接着由物理方程可以推出应力矩阵，其中为应力矩阵；最后通过虚功原理得出结点位移与结点力的关系，为单元刚度矩阵，表示单元抵抗变形的能力，它具有对称性、奇异性和分块性的特征。 3、整体分析，求整体刚度矩阵：为了求出整个离散化结构的总体方程组，只需对总体结构中的每一个节点建立平衡方程即可。并将结点力用结点位移便是，即可得出离散体含n个结点时的总体方程组为2n2n阶方程，即，也可以表示成，其中就是整体刚度矩阵，它是单元刚度矩阵的总和，是由各单元刚度矩阵的元素在相应位置上迭加而得，即所谓“对号入座”的方法。 4、约束条件引入和力的移置，解结构平衡方程组：由方程组可知，要求得需要知道,为奇异矩阵，通过引入约束条件，划去对应的行和列可以消除它的奇异性。而可以通过对力的移置转移到结点上。由此得到一个矩阵方程组，通过求解可以得到各个结点的位移。 5、求单元应力和位移：求得结点位移后，通过物理方程可以得到单元的应力。5、 结果初步分析通过学院自编有限元软件计算，得到位移变化图、应力矢量图以及相关的数据。通过位移变化图形和相关数据分析得出：隧道开挖后，硐周围岩体产生了向下及向临空面方向的位移。拱顶最大下沉量为29.32mm，底部围岩向上最大变形量为20.05mm。而对于应力矢量，最大主应力出现在洞壁出，在硐室周围岩体径向应力很小，环向应力很大，应力集中程度最大，随着远离硐室，围岩集中系数减小，逐渐接近天然应力状态。6、 图形 隧道开挖后，围岩位移变化图和应力矢量图如下所示。（完整成图见附图） 图3 硐室位移变化整体及局部构造图 图4 围岩应力矢量整体及局部构造图7、 隧道开挖对围岩的影响 隧道开挖对围岩的影响主要从围岩的位移变化特征及原因；应力特征及集中程度分析；围岩压力特征和隧道施工建议，这三个方面进行详细分析。 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.635-0.242-0.379-0.510-1.077-2.310-2.932表2 洞底中心线y变化结点号123457911洞底y位移（cm）2.0051.4771.0340.6990.0170.012-0.0820.000表3 拱、墙交点x变化结点号177178179180181182183184拱墙交点x（cm）-0.554-0.249-0.114-0.049-0.0050.0100.0140.011在硐室开挖后，由于围岩压力的作用，硐周围各点均产生向洞内的收缩位移，具体表现为洞顶产生向下的位移，出现沉降，而且拱顶中线处沉降量较大，并向边墙逐渐减小。同时，拱底面上各个结点均向上抬升，产生竖向位移，总体上呈现从边墙底部向底板中心逐渐增大的趋势。边墙产生向洞内收缩的位移。具体表1所示。同时，为了得到围岩收敛位移随位置的变化，取洞顶和洞底以及拱墙交点作为控制点，比较它们的特征位移，如表2和表3所示，得到相关规律：随着结点离围岩开挖轮廓线距离逐渐增大，结点向临空面收缩的位移量逐渐减小。 隧道开挖后，在洞壁上由于径向的约束解除，围岩向着临空面收缩产生位移。当位移达到围岩的极限位移过后，围岩产生松动坍塌。隧道围岩收敛量测是新奥法施工的重要量测项目，通过测试确定围岩变化趋势，确定何时进行支护。2、 围岩应力变化特征及集中程度分析 硐室开挖后，FINIT2.DAT文件可以反映出，在洞底有多个单元（1、12、23、34、45、67、78、89）压应力非常小（小于1MPa）。如表4和所示，应力有由负值（受压）变为正值（受拉）的可能性很大，这说明底部会有很大的变形，有破坏的可能性，因此需要适时的施做仰拱，防止衬砌陷落和抵御底部围岩压力。 同时观察，在硐室开挖后，在洞壁上的单元的最大和最小主应力如下表所示： 表5 洞壁单元最大和最小主应力单元1112131415161718191-2.04641-2.56026-3.43304 -4.61755 -6.11349-7.93552-10.12187-12.76577-15.09847 -14.61616 -0.77115-0.65056 -.57727-.61768 -.79350 -1.12713 -1.68088 -2.62647-3.41000-3.17247单元101111121131141151161171181191 -14.23406 -13.93872-13.73701 -13.64836 -13.69609 -13.89834-14.06139-13.76980-12.75043-11.15828 -2.69969 -2.17435 -1.71644 -1.40095 -1.27363 -1.37311 -1.66539 -1.87459 -1.79001 -1.540752由此，我们发现，在洞壁围岩上，最大主应力基本是呈现沿着硐室的环向，而最小主应力呈径向，且值非常小。这是由于洞壁围岩开挖使得径向应力解除，因此如果通过喷射混凝土提供切向摩擦力，那么会很大程度上加固围岩，这就是新奥法三大支柱之一（喷锚支护）的原理。对于应力集中，由硐室围岩应力图，我们得出在洞壁上围岩的应力集中程度最大，随着向深部围岩延伸，围岩应力集中程度降低，接近围岩天然应力状态。同时观察应力变化情况，在拱墙交错处、拱底和边墙交错处，围岩的应力集中现象最为明显，这些区域由于处于围岩不规则的地方，容易产生应力集中。对于拱底和边墙交汇处，可以通过施做仰拱，封闭断面的同时，形成圆弧形过渡区，减小应力集中程度。3、 围岩压力特征及施工建议 经过以上分析，得出围岩压力有以下特征：硐室开挖使得径向应力释放，环向应力增大，在靠近洞壁的地方围岩应力集中程度突出，尤其是拱墙交汇处和墙、底板交汇处。在研究区域中，最大应力也出现在洞壁处，同时在围岩底板附近压应力变得非常小并极可能转化为拉张应力。在拱顶和拱底向临空面收缩的位移量很大，当超过一定数值会产生松动坍塌，围岩由形变压力转化为松动压力。针对围岩单轴抗压强度低，围岩变形大的特点，在施工时尽量采用分部开挖法，如环形开挖留核心土法和双侧壁导坑法，并及时施做衬砌。爆破时，可以采用控制爆破的方法，减小对围岩的扰动。对拱顶根据位移收敛量测结果，在适当的时候打入锚杆，施做钢拱架，也可以铺设钢筋网。对于洞底，需要及时施做仰拱，尽快形成完整的圆形断面，提高衬砌的整体稳定性。当施工区段达到了二衬的安全距离后，整体施做二衬，作为强度储备。8、 结语通过岩土工程数值法的学习，我们学习了数值分析的过程。尤其是通过对有限单元法的详细讲解，我知道了大型通用有限元软件ansys、abaqus等的计算原理，不再是感