辽工大FLAC考试资料

FLACn;（新建）plblgr;（显示组信息）plsh;（画出显示）;1—建模genzoneradcylp0000p1add5000p2add0500p3add0050&;（gen生成，zone区域，定义起始点和X、Y、Z方向长度，add代 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 在P0坐标基础上加）dim1.751.751.751.75ratio1.01.01.01.5size825812fill;（确定隧洞，前两个1.75代表x=0面上X、Y方向半径是1.75，后两个1.75代表X=50的面上X、Y方向半径是1.75；ratio表示网格长度比例，顺序为X、Y、Z和围岩径向；size表示网格个数，顺序为径向、轴向、周向和围岩径向；fill表示填充）genzonereflectorig000norm001;（沿过(0.0.0)法线为(0.0.1)的面做对称，reflect是对称）genzonebrickp00050p1add5000p2add0500p3add00250&;（brick块体）ratio1.01.01.5size4254;2—指定本构模型及参数mmmodelmoh;（本构关系：模量摩尔库伦准则）defchushi;（定义初始）jian=yang/(2*(1+passion));（剪切模量表达式=杨氏模量/(2*(1+泊松比))）tiji=yang/(3*(1-2*passion));（体积模量表达式=杨氏模量/(3*(1-2*泊松比))）enddefcailiao1;（定义材料1）command;（命令开始）setyang=3.2e9setpassion=0.25chushi;（调用初始）printjian;（打印即屏幕输出）printtijipropshearjianbulktijiranz-4,4;（-4到4之间为同种材料，prop是定义，shear是剪切模量，bulk是体积模量，ran是范围）propcoh=2.3e6fric=45.0tens=0.5e6ranz-4,4;（coh是粘结力，fric是内摩擦角，tens是抗拉强度）iniden0.0025e6ranz-4,4;（ini是初始，den是密度）endcommand;（命令结束）endcailiao1;（调用材料1）defcailiao2commandsetyang=4.8e9setpassion=0.24chushiprintjianprinttijipropshearjianbulktijiranz4,50propshearjianbulktijiranz-4,-50propcoh=2.8e6fric=48.0tens=0.7e6ranz4,50propcoh=2.8e6fric=48.0tens=0.7e6ranz-4,-50iniden0.0025e6ranz4,50iniden0.0025e6ranz-4,-50endcommandendcailiao2defcailiao3commandsetyang=5.5e9setpassion=0.23Chushiprintjianprinttijipropshearjianbulktijiranz50,30propcoh=3.2e6fric=50.0tens=0.95e6ranz50,300iniden0.0025e6ranz50,300endcommandendcailiao3;3—施加边界条件setgrav00-10;（施加重力）iniszz=-7.5e6grad0,0,2.5e4;（Z方向应力=szz+容重*z，Z方向应力为地表应力得0，容重=密度*g，所以0=szz+0.0025e6*(-10)*300，szz=-7.5e6，grad是梯度，梯度=容重）inisxx=-11.25e6grad0,0,3.75e4;（7.5*1.5=11.252.5*1.5=3.75，定义初始应力场lamdx=1.5lamdy=1.25为侧压力系数）inisyy=-9.375e6grad0,0,3.125e4;（7.5*1.25=9.3752.5*1.25=3.125）fixxrangex-0.10.1;（施加位移、速度场边界条件）fixxrangex49.950.1fixyrangey-0.10.1fixyrangey49.950.1fixxyzrangez-50.1-49.9;（固定z=-50面上X、Y、Z方向的速度）;4—生成初始/原岩地应力场sol;（sol是求解）sav05-000.sav;（sav是保存，05-000.sav是文件名，后缀为.sav）;5—开挖支护前模型位移场、速度场清零inixdisp0ydisp0zdisp0xvel0yvel0zvel0;(xdisp表示X方向位移为0，xvel表示X方向速度为0);6—开挖过程的模拟modelnullrangecylend10,0,0end2040radius1.75;（开挖，cyl表示圆柱形轴线起点（0,0,0）挖到（0,4,0），radius表示半径）histgpzdisp021.751;（拱顶竖向位移，监测深2米，高1.751米处的围岩情况，hist表示监测，gp表示节点）histgpzvel021.751;（拱顶竖向速度）histgpxdisp1.75120;（腰部水平位移，监测深2米，宽1.751米处的围岩情况）histgpxvel1.75120;（腰部水平速度）histzonesmax021.751;（拱顶最大主应力）histzonesmax1.75120;（腰部最大主应力）histzonesmin021.751;（拱顶最小主应力）histzonesmin1.75120;（腰部最小主应力）histunba;（unba表示不平衡）solvesave05-tunnel-4m-00.sav;7—支护modelnullrangecylend10,4,0end2080radius1.751;（开挖，4到8米）selshellrangecylend1000end2040radius1.751;（sel表示支护）selshellpropiso=(25.3e9,0.266)thick=0.25;（iso表示各向同性，25.3e9是弹性模量，0.266是泊松比，thick表示厚度）selrecoversurfacesurfx010;（定义局部坐标系）;selcableid1begin021.75end024.25nseg5;（nseg表示份数）;selcableid2begin1.7520end4.2520nseg5;selcableproemod45e9xcarea1.57e-3gr_per1.&;（emod是弹性模量，xcarea是横截面积，gr_per是水泥浆外圈周长）;yten25e5gr_k17.5e6gr_c20e5ranid=1;（yten是抗拉强度，gr_k是单位长度上锚固层剪切刚度，gr_c是锚固层剪切强度，内聚力）;selcableproemod45e9xcarea1.57e-3gr_per1.&;yten25e5gr_k17.5e6gr_c20e5ranid=2defsel_cable;（cable是锚杆）a=pi/4.0;（定义锚杆间夹角）r=1.75;（定义半径）y0=2.0;（锚杆位置）s=2.5;（锚杆长度）loopn(1,5);（循环）x0=r*sin(a*(n-1));（定义锚杆起点坐标）y0=y0z0=r*cos(a*(n-1))x1=(r+s)*sin(a*(n-1));（定义锚杆终点坐标）y1=y0z1=(r+s)*cos(a*(n-1))commandselcableid1beginx0y0z0endx1y1z1nseg5selcableproemod45e9xcarea1.57e-3gr_per1.&;（gr_per要大于xcarea）yten25e5gr_k17.5e6gr_c20e5ranid=1endcommandendloopendsel_cablesolsave05-tunnel-8m-00.savmodelnullrangecylend10,8,0end20120radius1.751selshellrangecylend1040end2080radius1.751selshellpropiso=(25.3e9,0.266)thick=0.25selrecoversurfacesurfx010defsel_cable1a=pi/4.0r=1.75y0=4.0s=2.5loopn(1,5)x0=r*sin(a*(n-1))y0=y0z0=r*cos(a*(n-1))x1=(r+s)*sin(a*(n-1))y1=y0z1=(r+s)*cos(a*(n-1))commandselcableid1beginx0y0z0endx1y1z1nseg5selcableproemod45e9xcarea1.57e-3gr_per1.&yten25e5gr_k17.5e6gr_c20e5ranid=1endcommandendloopendsel_cable1defsel_cable2a=pi/4.0r=1.75y0=6.0s=2.5loopn(1,5)x0=r*sin(a*(n-1))y0=y0z0=r*cos(a*(n-1))x1=(r+s)*sin(a*(n-1))y1=y0z1=(r+s)*cos(a*(n-1))commandselcableid1beginx0y0z0endx1y1z1nseg5selcableproemod45e9xcarea1.57e-3gr_per1.&yten25e5gr_k17.5e6gr_c20e5ranid=1endcommandendloopendsel_cable2solsave05-tunnel-12m-00.savret;（返回）基本操作及命令：*（必考）显示锚杆轴力：plselcableforce*（必考）计算安全系数：solvefos显示坐标轴，激活显示窗口：Plotitems->Add->Axes->X=-1,Y=0,Z=-1；敲击X、Y、Z图像旋转，shift+X、Y、Z反向旋转；敲击M放大、shift+M缩小，或Ctrl+Z（不是按住），鼠标变成+，按住左键圈住相放大区域，松开左键即可放大；调用文件File->Call；调整至合适大小Ctrl+R；注释前要加英文分号“;”；显示分区：点击显示窗口，Plotitems->Add->Block->Property->选择分类->OK->OK显示计算结果：激活显示窗口->Plotitems->Add->Contour->Displacement->Magnitude读取上次运行后保存的文件用res，例如：res05-000.sav显示：激活显示窗口->Plotitems->Add->Block->State显示锚杆受力：激活显示窗口->Plotitems->Add->StructuralElements->Cables->Force显示拉：plblsttension显示剪：plblstshear显示总：plblst显示总位移（outlon表示显示网格）：plcondispoutlon显示X位移：plconxdispoutlon显示X方向正应力：plconsxxoutlon显示最大主应力：plconsmaxoutlon显示壳单元：plselshall