運用flac3d模拟沟渠开挖前后岩体的应力应变报告.docVIP

运用flac3d模拟沟渠开挖前后岩体的应力应变报告 一、问题的提出 在天然山体与平原的的结合处开挖一条10m×4m×3m的沟渠，进行应力应变场的分析。 二、问题的分析与处理 本模型的分析共分为三个阶段： 第一阶段——前处理阶段：创建初始几何模型并划分网格，并在flac3d中显示网格体，然后定义材料的模型和材料参数，再定义边界条件和初始条件，这样便用flac3d建立起了计算的模型。 第二阶段——求解：在第一阶段建立好模型的基础上对模型进行变量的监控和求解运算。 第三阶段——后处理阶段：通过建立各种视口显示出各种运算的结果，以此判断模型的求解是否收敛。 对于沟渠的开挖模型及其运算，分三个阶段说明其生成及运算过程： 1、前处理阶段 由图可知，该模型可由两个radtunnel模型合并而成，分别以左侧模型p0（10，0,20）和右侧模型p0（10,10,10）为起始点。 步骤： 1.打开flac3d软件，输入new命令全部清空原来数据，接着输入代码，如下：generate zone radtunnel p0 10 0 20 p1 0 0 20 p2 10 10 20 p3 10 0 0 p4 0 10 20 p5 10 10 0 p6 0 0 0 p7 0 10 0 p8 8 0 20 p9 10 0 7 p10 8 10 20 p11 10 10 7 p12 8 0 7 p13 8 10 7 generate zone radtunnel p0 10 10 10 p1 20 10 10 p2 10 0 10 p3 10 10 0 p4 20 0 10 p5 10 0 0 p6 20 10 0 p7 20 0 0 p8 12 10 10 p9 10 10 7 p10 12 0 10 p11 10 0 7 p12 12 10 7 p13 12 0 7 plot add surface blue show 生成网格模型如图1 图1 几何网格模型 2.定义材料的本构模型，本题为摩尔-库伦模型，然后对材料参数（体积模量，内摩擦角和剪胀角，粘聚力，抗拉强度）进行赋值，如没有输入其中一种参数，则系统默认为0。代码如下: model mohr property bulk=1e8 shear=0.3e8 friction=35 property cohesion=1e10 3.定义边界条件和初始条件，用到set,fix,ini,apply命令，先加载重力再定义材料的密度，然后固定模型的5个面，使其保持稳定（不固定z轴的正面），代码如下： initial density=2000 set gravity 0 0 -9.81 fix x range x -0.1 0.1 fix x range x 19.9 20.1 fix y range y -0.1 0.1 fix y range y 9.9 10.1 fix z range z -0.1 0.1 求解阶段 该题研究三种状态下岩体模型的应力应变，即天然状态，饱和状态，地震状态。分别讨论如下： 1.在天然状态下，以上的代码即可求解，为了记录一些变量或参数的变化，在计算过程中运用History命令记录了最大不平衡力和点（5,5,5）的位移，令系统运算500步，每隔8步记录一次，以此来模拟模型的变化，代码如下： hist n=8 hist unbal hist gp zdis 5,5,5 step 500 再输入plot命令 再输入history 1 命令，显示位移等线值图如图2 图2 最大不平衡力采样记录图 输入history 2 显示点（5,5,5）z轴位移采样记录如图3 图3 点（5,5,5）z轴位移采样记录图 再输入以下命令创建任意剖面上的等值线图 create trench create grav1 set plane dip 90 dd 0 origin 3 5 0 set rot 50 25 130 set center 7 10 5 add bound behind add bcontour szz plane add axes black show 如图4 图4 剖面上的垂直应力图 2.在饱和状态下，更改密度，然后再重复上述运算，得到结果。代码如下： initial density=3000 hist n=10 hist unbal hist gp zdis 5,5,5 step 500 plot history 1 显示最大不平衡力采样记录如图5 图5最大不平衡力采样记录图 再输入history 2 显示点（5,5,5）z轴位移采样记录如图6 图6点（5,5,5）z轴位移采样记录图 再输入以下命令流创建任意剖面上的等值线图