计算机数值模拟实验论文探究.docx

研究生课程考核试卷科目：计算机数值模拟实验教师：姓名：学号：专业：类别：上课时间：2016年4月至2016年5月t考生成绩：卷面成绩平时成绩课程综合成绩阅卷评语：阅卷教师 (签名)基于FLAC3D的矩形隧道开挖分析1引言FLAC3D是一个开放的系统，通过其命令驱动模式使得用户可以参与从网格模型的建立、边界条件的设置、参数的调试到计算结果的输出等几乎全部求解过程。本文拟从矩形隧道的开挖模拟全过程来更好地熟悉和掌握FLAC3D的使用。隧道工程的开挖，会使开挖区周围岩土体的应力重新进行分布，开挖区的应力释放及变形也对边坡的稳定性造成影响，因此隧道开挖具有一定实际意义。1.1问题描述本文讨论“矩形断面隧道施工设计”，矩形断面宽5m，高4m，考虑隧道不同埋深，不同水平应力对隧道围岩稳定性影响。模型采用各项同性岩体，Mohr-Coulomb强度准则，岩体力学参数参考相关文献自己选取。注：考虑埋深与水平应力不少于三种情况。表1.1为本文在编程中选取的力学参数。表1.1 岩层物理力学参数岩性ρ/kg·m-3K/Pac/MPaφ/(o)Rm/MPaG/Pa泥岩25305e90.828.610.11e92 建模过程2.1模型建立由于矩形模型具有对称性，可以采用1/4模型进行分析。首先建立坐标系，坐标系的远点O设置在模型对称轴的点，水平向右为x方向，竖直向上为z方向，垂直于分析平面的方向为y方向。根据设计要求，建模使用平行六面体隧道外围渐变放射网格单元，如下图2.1，建立网格单元后利用镜像生成完整网格。图2.1 平行六面体隧道外围渐变放射网格单元网格建立的命令如下：new;======建立模型========generate zone radtunnel p0 0 0 0 p1 25 0 0 p2 0 10 0 p3 0 0 20 dimension 2.5 2 2.5 2 size 5 20 5 25 gr 围岩 fill gr 巷道gen zone refgen zone ref dip 90gen zone ref dip 90 dd 90 ;;;;;;镜像生成网格plot add surface yellow savmodel.sav2.2 初始应力计算网格建立后，将网格赋值为Mohr模型，对边界条件和力学参数进行设定。本设计中设定的力学参数如表1.1中所示，对x两侧、y两侧、z底面的边界进行速度约束。newrestoremodel.sav;=========设置强度准则=======modelmohrprop bulk=5e9 shear=1e9 coh 8e5 fri 28.61 ten 1e5 initial density 2.53e3;=========初始边界条件===========fix x range x -24.9 -25.1fix x range x 24.9 25.1fix y range y -10.1 -9.9fix y range y 9.9 10.1fix z range z -19.9 -20.1set gravity 0 0 -9.81applyszz -5.02e6 range z 19.9 20.1iniszz -5.02e6inisxx -1e6inisyy -1e6 plblstsetmech ratio 1e-4solvesavestress_initial.sav由于要考虑不同埋深和不同水平应力的影响，本文中选取埋深分别为200m、400m、600m，垂直应力和水平应力的取值如表2.1所示。表2.1 垂直应力和水平应力取值编号Szz（Pa）sxx（Pa）syy（Pa）1#-5.02e6-1e5-1e52#-5.02e6-1e6-1e63#-5.02e6-1e7-1e74#-10.04e6-1e6-1e65#-15.07e6-1e6-1e6在初始应力设置过程中需要对程序中相应参数进行修改，编号1、2、3三组为同一埋深下三种不同水平应力设置，编号2、4、5为同一水平应力下三组不同埋深（即垂直应力）设置。编号1#初始应力参数设置为:applyszz -5.02e6 range z 19.9 20.1iniszz -5.02e6inisxx -1e5inisyy -1e5编号2#初始应力参数设置为:applyszz -5.02e6 range z 19.9 20.1iniszz -5.02e6inisxx -1e6inisyy -1e6编号3#初始应力参数设置为:applyszz -5.02e6 range z 19.9 20.1iniszz -5.02e6inis