辽工大FLAC考试资料.docVIP

FLAC n;（新建） pl bl gr;（显示组信息） pl sh;（画出显示） ;1—建模 gen zone radcyl p0 0 0 0 p1 add 50 0 0 p2 add 0 50 0 p3 add 0 0 50 ;（gen生成，zone区域，定义起始点和X、Y、Z方向长度，add代表在P0坐标基础上加） dim 1.75 1.75 1.75 1.75 ratio 1.0 1.0 1.0 1.5 size 8 25 8 12 fill;（确定隧洞，前两个1.75代表x=0面上X、Y方向半径是1.75，后两个1.75代表X=50的面上X、Y方向半径是1.75；ratio表示网格长度比例，顺序为X、Y、Z和围岩径向；size表示网格个数，顺序为径向、轴向、周向和围岩径向；fill表示填充） gen zone reflect orig 0 0 0 norm 0 0 1;（沿过(0.0.0)法线为(0.0.1)的面做对称，reflect是对称） gen zone brick p0 0 0 50 p1 add 50 0 0 p2 add 0 50 0 p3 add 0 0 250 ;（brick块体） ratio 1.0 1.0 1.5 size 4 25 4 ;2—指定本构模型及参数mm model moh;（本构关系：模量 摩尔库伦准则） def chushi;（定义 初始） jian=yang/(2*(1+passion));（剪切模量表达式=杨氏模量/(2*(1+泊松比))） tiji=yang/(3*(1-2*passion));（体积模量表达式=杨氏模量/(3*(1-2*泊松比))） end def cailiao1;（定义 材料1） command;（命令开始） set yang=3.2e9 set passion=0.25 chushi;（调用 初始） print jian;（打印 即屏幕输出） print tiji prop shear jian bulk tiji ran z -4,4;（-4到4之间为同种材料，prop是定义，shear是剪切模量，bulk是体积模量，ran是范围） prop coh=2.3e6 fric=45.0 tens=0.5e6 ran z -4,4;（coh是粘结力，fric是内摩擦角，tens是抗拉强度） ini den 0.0025e6 ran z -4,4;（ini是初始，den是密度） endcommand;（命令结束） end cailiao1;（调用 材料1） def cailiao2 command set yang=4.8e9 set passion=0.24 chushi print jian print tiji prop shear jian bulk tiji ran z 4,50 prop shear jian bulk tiji ran z -4,-50 prop coh=2.8e6 fric=48.0 tens=0.7e6 ran z 4,50 prop coh=2.8e6 fric=48.0 tens=0.7e6 ran z -4,-50 ini den 0.0025e6 ran z 4,50 ini den 0.0025e6 ran z -4,-50 endcommand end cailiao2 def cailiao3 command set yang=5.5e9 set passion=0.23 Chushi print jian print tiji prop shear jian bulk tiji ran z 50,30 prop coh=3.2e6 fric=50.0 tens=0.95e6 ran z 50,300 ini den 0.0025e6 ran z 50,300 endcommand end cailiao3 ;3—施加边界条件 set grav 0 0 -10;（施加重力） ini szz=-7.5e6 grad 0,0,2.5e4;（Z方向应力=szz+容重*z，Z方向应力为地表应力得0，容重=密度*g，所以0=szz+0.0025e6*(-10)*300，szz=-7.5e6，grad是梯度，梯度=容重） ini sxx=-11.25e6 grad 0,0,3.75e4;（7.5*1.5=11.25 2.5*1.5=3.75，定义初始应力场 lamdx=1.5 lamdy=1.25 为侧压力系数） ini syy=-9.375e6 grad 0,0,3.125e4;（7.5*1.25=9.375 2.5*1.25=3.125） fix x range x -0.1 0.1;（施加位移、速度场边界条件） f