《计算机数值模拟实验论文》由会员分享，可在线阅读，更多相关《计算机数值模拟实验论文（14页珍藏版）》请在金锄头文库上搜索。
1、研究生课程考核试卷科目： 计算机数值模拟实验 教师： 姓名： 学号： 专业： 类别： 上课时间： 2016年4月 至2016年5月 t考生成绩：卷面成绩平时成绩课程综合成绩阅卷评语：阅卷教师 (签名)重庆大学课程论文 基于FLAC3D的矩形隧道开挖设计基于FLAC3D的矩形隧道开挖分析1引言FLAC3D是一个开放的系统，通过其命令驱动模式使得用户可以参与从网格模型的建立、边界条件的设置、参数的调试到计算结果的输出等几乎全部求解过程。本文拟从矩形隧道的开挖模拟全过程来更好地熟悉和掌握FLAC3D的使用。隧道工程的开挖，会使开挖区周围岩土体的应力重新进行分布，开挖区的应力释放及变形也对边坡的稳定性
2、造成影响，因此隧道开挖具有一定实际意义。1.1 问题描述本文讨论“矩形断面隧道施工设计”，矩形断面宽5m，高4m，考虑隧道不同埋深，不同水平应力对隧道围岩稳定性影响。模型采用各项同性岩体，Mohr-Coulomb强度准则，岩体力学参数参考相关文献自己选取。注：考虑埋深与水平应力不少于三种情况。表1.1为本文在编程中选取的力学参数。表1.1 岩层物理力学参数岩性/kgm-3K/Pac/MPa/()Rm/MPaG/Pa泥岩25305e90.828.610.11e92 建模过程2.1模型建立由于矩形模型具有对称性，可以采用1/4模型进行分析。首先建立坐标系，坐标系的远点O设置在模型对称轴的点，水平向
3、右为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
4、 ;镜像生成网格plot add surface yellow sav model.sav2.2 初始应力计算网格建立后，将网格赋值为Mohr模型，对边界条件和力学参数进行设定。本设计中设定的力学参数如表1.1中所示，对x两侧、y两侧、z底面的边界进行速度约束。newrestore model.sav;=设置强度准则=model mohrprop 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.1f
5、ix 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.81apply szz -5.02e6 range z 19.9 20.1ini szz -5.02e6ini sxx -1e6ini syy -1e6 pl bl stset mech ratio 1e-4solve save stress_initial.sav由于要考虑不同埋深和不同水平应力的影响，本文中选取埋深分别为200m、400m、600m，垂直应力和水平应力的取值如表2.1所示。表2.1 垂直应力和水平应
6、力取值编号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#初始应力参数设置为:apply szz -5.02e6 range z 19.9 20.1ini szz -5.02e6ini sxx -1e5ini syy -1e5编号2#初始应力参数设置为:app
7、ly szz -5.02e6 range z 19.9 20.1ini szz -5.02e6ini sxx -1e6ini syy -1e6编号3#初始应力参数设置为:apply szz -5.02e6 range z 19.9 20.1ini szz -5.02e6ini sxx -1e7ini syy -1e7编号4#初始应力参数设置为:apply szz -10.4e6 range z 19.9 20.1ini szz -10.4e6ini sxx -1e6ini syy -1e6编号5#初始应力参数设置为:apply szz -15.07e6 range z 19.9 20.1ini
8、 szz -15.07e6ini sxx -1e6ini syy -1e62.3 开挖过程模拟设置好初始应力后，在进行开挖模拟前要将初始应力计算过程中产生的节点位移的速度进行清零处理，程序如下：newrestore stress_initial.sav;开挖;ini xdis 0 ydis 0 zdis 0ini xvel 0 yvel 0 zvel 0model null range gr 巷道solvesave excavation.sav3 结果分析3.1 初始应力分析首先对模型施加初始应力后计算得到的竖向应力和水平应力进行分析，图3.1、3.2为不同埋深和不同水平应力下的初始竖向应力云
9、图和初始水平应力云图。(a) 1# (b) 2#(c) 3# (d) 4#(e) 5#图3.1 初始竖向应力云图通过对图3.1中(a)、(b)、(c)可以看出在同一埋深200m下，1e5、1e6、1e7Pa水平应力下，最大竖向应力先增加后减少，但差别不大，分布在模型底端，最大值大约为6e6Pa，但竖向应力云图的分布稍有差别，说明同一埋深下不同水平应力会导致竖向应力发生变化。通过对图3.1中(b)、(d)、(e)进行比较发现，在同一水平应力1e6Pa下，200m埋深下最大竖向应力为6.07e6Pa，400米埋深下最大竖向应力为1.09e7Pa，600米埋深下最大竖向应力为1.59e7Pa，即随着
10、埋深成倍增加，最大竖向应力也大致成倍增加，这也与事实情况相符。同时我们发现，随着埋深增加，竖向应力开始向模型XZ平面竖直对称面集中，说明竖直应力的增加对模型中部的影响高于对边界处的影响。表3.1 初始竖向应力、初始水平应力最值编号最大初始竖向应力（Pa）最小初始竖向应力（Pa）最大初始水平应力（Pa）最小初始水平应力（Pa）1#6.02 e65.10 e61.95 e61.25 e62#6.06 e65.02 e61.72 e69.93 e53#5.82 e65.01 e61.05 e71.00 e74#1.09 e71.00 e74.24 e63.54 e65#1.59 e71.50 e77
11、.04 e66.38 e6(a) 1# (b) 2#(c) 3# (d) 4#(e) 5#图3.2 初始水平应力云图在图3.2中可以看到随着深度的增加，水平应力增大，在图3.2(a)、(b)、(c)中，同一埋深下，随着水平应力的增加，最大水平应力和最小水平应力均先减小后增大，初始水平应力的变化并不随施加的水平应力的变化而线性变化。图3.2中(b)、(d)、(e)显示同一水平应力下，随着埋深线性增加，最大水平应力和最小水平应力均大致呈线性增加，且水平应力云图分布发生变化，模型水平应力分布呈扩散趋势。3.2 开挖分析3.2.1 开挖后塑性区分析在FLAC3D中采用摩尔库伦模型进行模拟时，可以通过b
12、lock state显示哪些应力符合屈服准则的区域（或称塑性区），以观察潜在破坏区域的范围。塑性区标识以不同的颜色显示两种类型的破坏机制即剪切破坏（shear failure）和拉伸破坏（tensile failure）：某一区域的应力正位于屈服面上，或者说正处于破坏状态时，以shear-n 或tension-n标识；某一区域在计算过程中曾进入过屈服状态，但现在已经退出，以shear-p 或tension-p标识。 (a) 1# (b) 2#(c) 3# (d) 4#(e) 5#图3.3开挖后塑性区分布图从图3.3中可以看到本文中的5个模型都存在塑性区，表明塑性流动正在发生。由于本文中没有对隧
13、道进行监护，隧道破坏也是正常的。但5个模型中塑性区的分布存在差别，由图3.3(a)、(b)、(c)我们可以看到，随着水平应力的增加，塑性区范围缩小，剪切破坏区变小，未破坏区域变大。比较3.2中(b)、(d)、(e)发现同一水平应力下，随着埋深增加，红色和灰色区域增多，蓝色区域和黄色区域变小，塑性流动发生的区域增多，同时可以发现隧道上下底面先进入屈服状态，即上下底面先发生塑性破坏。3.2.2开挖后竖向应力云图分析开挖后竖向应力在隧道左侧和右侧岩壁中心线附近达到最大值，在隧道上底面和下底面中心线附近有最小值。由表3.2可以清楚地看到，随着水平应力的的增加，最大竖向应力减小；随着埋深的增加，最大竖向应力增大。对模型底部来说，四角的竖向应力比底部中部受的竖向应力要大，由隧道左壁向模型左底角、隧道右壁向模型右底角区域有竖向应力集中现象。隧道上底面向模型顶部部分和隧道底部向模型底部部分受竖向应力较小。(a) 1#
