岩土工程数值法分析实例-有限元原理.docVIP

PAGE 1 岩土工程数值法 班 级 ： 63 专 业 ： 隧道与地下工程 姓 名 ： 学 号 ： 630 吉林大学建设工程学院 年 月 日 目录 一、问题提出………………………………………………………3 二、围岩离散化……………………………………………………4 三、数据准备………………………………………………………5 四、计算过程………………………………………………………6 五、结果初步分析…………………………………………………7 六、图形……………………………………………………………7 七、隧道开挖对围岩的影响………………………………………10 八、结语……………………………………………………………12 九、参考文献………………………………………………………12 随着我国经济快速发展，各种隧道、公路、铁路、房屋以及其它基础设施进入了一个高速建设的阶段，随之而来的是土木工程的跨越式发展。土木工程的设计和研究手段也有了很大的提高，从以前的基于经验的设计理论逐渐过渡到定量与定性相结合的反分析计算理论。目前为止，土木工程的研究方法主要有以下五类：类比法；解析法；模型模拟（物理模型方法）；现场监控量测；数值法（数值模拟）。 通过岩土工程数值法这门课程，我们系统的学习了数值法中的有限单元法的原理以及它在岩土工程中的应用。随着计算机的普及和运算速度的提升，为弹性力学的数值解法开辟了广阔的领域，尤其在隧道工程中，采用有限单元法分析都得到了满意的结果。目前，有限单元法已经是解决不同岩体结构、围岩与支护相互作用、隧道围岩压力、围岩应力和变形、围岩破坏过程与破坏机制的主要方法。本文将针对一个隧道开挖实例，应用有限单元法进行位移、应力等相关参数的分析。 问题提出 在岩土体中修建隧道是一件十分复杂的工程。因为岩土体是地壳内外力长期作用下形成的一种复杂的地质体，具有天然应力、非均质、不连续、各向异性等特点，从而表现在力学性质上具有非线性、剪胀性、蠕变性等。而有限单元法可以将岩土体复杂多变的力学性质，基本地质因素、复杂和混合的边界条件、岩土体与工程结构物的组合作用等问题统筹考虑，以得到接近实际的数值解答。 目前，隧道施工和设计都是基于“新奥法”，新奥法的核心是充分发挥围岩的承载能力，将围岩视为承载的主体。随之而来的是如何确定围岩收敛的极限位移，如何确定衬砌的支护时间，如何判定围岩应力的集中程度等问题。本文基于以下条件进行隧道开挖后的围岩进行分析。 本隧道断面为曲墙式，拱部半径为5m，下部为10m×6m的矩形，跨度为10m，隧道埋深200m。已知岩体参数为：岩体弹性模量，泊松比，岩体初始粘结力，初始摩擦角，岩体残余粘聚力，残余摩擦角，岩体容重，岩体单轴抗压强度。具体见图1。 围岩离散化 有限单元法分析的基础是将弹性体分划为有限大小、彼此只在有限个点相联结、由有限个单元的组合体来研究的，也就是把一个连续的物体或结构分割成有限个有限大小的多边形或者多面体。 因此有限单元类型的选择以及单元的划分数对求解的精度起着至关重要的作用。一般来说，划分的单元愈多，愈细密，精度也就愈高，也就愈能反映结构的实际情况。隧道开挖是一个平面应变问题，可以简化为平面问题进行分析，同时由于隧道断面是一个轴对称图形，选取右半部分进行分析。离散化拟采用四边形单元，因为四边形单元是具有较高次位移模式的单元，能够降低离散化误差，更好地反映岩土体的位移和应力状态。具体的离散过程可以分为以下七个步骤： （1）研究区域大小的确定：对于圆形隧洞，由岩体力学的知识可知，处于6倍洞径以外的岩体不受硐室开挖的影响。本隧道开挖形状不规则，因而考虑取开挖轮廓线外7倍硐径（35m）围岩为研究区域。具体见图2。 （2）单元数量的确定：由于本次使用的有限元计算程序规定的计算结点数为400到500之间，需要对结点数进行控制。对于硐室开挖轮廓线采取1m划分一次，因而硐室底板划为5段，侧边墙划为段，拱部圆弧划为8段。在径向划分为7段。所以总共划分420个结点，380个单元。 （3）单元的疏密变化：尽量注意疏密变化的均匀性，注意单元大小的过渡。沿着隧道径向单元分布由密到疏，单元变长比例为1.1：1.2：1.3：1.4····2.0（逐次增加0.1），这样能更好的反映硐室开挖的影响。 （4）单元编号：单元编号采取逆时针依顺序，共380个单元。 （5）结点编号：结点编号同样采用逆时针，共420个结点。 （6）确定结点坐标：本次只需确定图形轮