中南大学隧道flac2d课程作业.docVIP

FLAC 2D简介 FLAC是一个用于工程力学计算的二维显示有限差分程序。最早由岩土和采矿工程师开发，不过目前应用范围广泛。 适用于模拟土体、岩石或其他弹塑性材料。材料可以由单元以及由网格区域代替，网格区域可形成任意形状以适应模型的要求。 单元的力学行为由预先设定的线性或非线性应力/应变关系确定。FLAC采用的显式拉格朗日差分法和混合离散技术使得材料的塑性变形甚至流动破坏的模拟变得容易且精确。由于计算过程中不产生刚度矩阵，二维计算规模可变得庞大且对内存要求不高。 1.1 FLAC的特点: (1)界面单元能模拟特定面上的滑动和分离；(结构面、断层) (2)平面应变、平面应力和轴对称问题； (3)能处理流固耦合问题； (4)结构单元可模拟多种结构形式(如隧道衬砌、锚杆、桩基等)； (5)动力问题；(地震、冲击、运动的荷载) 粘弹塑性问题； 1.2 FLAC 2D优劣势 : 与有限单元法的对比: 优势： (1)混合离散技术更适应于模拟塑性变形和塑性破坏； (2)计算过程中引入运动方程，使得数值计算的稳定性增强； (3)显式的求解方式使得模拟非线性问题更加高效； (4)对本构模型的适应性强，而有限元方法对不同本构模型须采取 的不同求解方式； (5)网格以行、列的方式排列更有利于数据提取和分析； 劣势： (1)对于线性问题的模拟相比有限元方法稍慢； (2)求解速度与模型的最大和最小周期之比成正比，对于特定问题的求解效率低。 工程概况 依据图1给出的Ⅴ级围岩洞口段隧道衬砌支护结构图，参考图2中示意的数值模型尺寸和边界条件，进行隧道开挖围岩稳定性分析和支护结构强度检算（不考虑系统锚杆）：。 图1 Ⅴ级围岩洞口段隧道衬砌支护结构图 图2 数值模型尺寸和边界条件 2.1计算工况（推荐）： (1)自重应力平衡； (2)全断面开挖隧道，应力释放30%，计算平衡； (3)施做初期支护，应力释放30%，计算平衡； (4)施做二次衬砌，应力释放剩余的40%，计算平衡。 2.2计算参数可参考表1: 表1 数值计算参数 材料名称 计算参数 Ⅴ级围岩 E=1.30GPa γ=22KN/ m3 μ=0.4 c=100KPa φ=25°σt =10KPa C20喷射混凝土 E=2.43e10Pa μ=0.2 C25钢筋混凝土 E=2.95e10Pa μ=0.2 表中：E—弹性模量；γ—容重；μ—泊松比；c —粘结力；φ—内摩擦角；σt —岩体抗拉强度。 说明：隧道采用新奥法施工，其中Ⅴ级围岩段采用全断面开挖施工方法，预留变形量13cm。 初步模型建立与数值计算 3.1 模型建立 根据工程要求在AUTOCAD中根据一定的模型关系画出模型简图（图3-1-1）并记录得到各控制点的坐标记录如下表3-1-1： 表3-1-1 横截面控制点坐标 控制点 X坐标 Y坐标 A 5.6792 0.8995 B -5.6792 0.8995 C -5.6098 -1.7482 D 5.6098 -1.7482 O 0 0 O1 -1.5013 -0.2378 O2 1.5013 -0.2378 O3 0 26.1834 图3-1-1 隧道横截面模型简图 （O3的坐标由于图幅的限制图没显示） 根据以上坐标及洞口附近Ⅴ级围岩段地形，选取其中一个隧道断面建立数值计算模型，本计算建立的数值模型，在FLAC 2D中划分100X80个单元网格，具体模型如图所示： 图3.1.1隧道及其周围岩体的截面图 3.2 计算参数 根据隧道围岩勘察测量情况和支护结构设计材料特性，本计算采用的参数具体见表3-2-1。 表3-2-1 数值计算参数 材料名称 计算参数 Ⅴ级围岩 E=1.30GPa γ=22KN/ m3 μ=0.4 c=100KPa φ=25°σt =10KPa C20喷射混凝土 E=2.43e10Pa μ=0.2 C25钢筋混凝土 E=2.95e10Pa μ=0.2 表中：E—弹性模量；γ—容重；μ—泊松比；c —粘结力；φ—内 摩擦角；σ t —岩体抗拉强度。 初期支护C20喷射混凝土厚度25cm，二次衬砌C25钢筋混凝土厚度45cm。 3.3施工工序 隧道采用新奥法施工中的全断面开挖施工方法，具体施工工序： ①全断面开挖隧道； ②施做初期支护； ③施做二次衬砌； 应力变形模型建立与计算分析 4.1自重应力平衡 利用FLAC 2D模拟计算计算范围内山体在隧道开挖前在自重作用下的应力平衡。数值计算结果如下： 图4-1-1 自重