FLAC3D-实例命令流1

第1部分命令流按照顺序进行2-1定义一个FISH函数newdefabcabc=25*3+5Endprintabc2-2使用一个变量newdefabchh=25abc=hh*3+5EndPrinthhPrintabc2-3对变量和函数的理解newdefabchh=25abc=hh*3+5Endsetabc=0hh=0printhhprintabcprinthhnewdefabcabc=hh*3+5endsethh=25printabcsetabc=0hh=0printhhprintabcprinthh2-4获取变量的历史记录newgenzonebricksize121modelmohrpropshear=1e8bulk=2e8cohes=1e5tens=1e10fixxyzrangey-0.10.1applyyvel-1e-5rangey1.92.1plotsetrotation0045plotblockgroupdefget_adad1=gp_near(0,2,0)ad2=gp_near(1,2,0)ad3=gp_near(0,2,1)ad4=gp_near(1,2,1)endget_addefloadload=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4)endhistloadhistgpydis0,2,0step1000plothis1vs-22-5用FISH函数计算体积模量和剪砌模量newdefderives_mod=y_mod/(2.0*(1.0+p_ratio))b_mod=y_mod/(3.0*(1.0-2.0*p_ratio))endsety_mod=5e8p_ratio=0.25deriveprintb_modprints_mod2-6在FLAC输入中使用符号变量Newdefderives_mod=y_mod/(2.0*(1.0+p_ratio))b_mod=y_mod/(3.0*(1.0-2.0*p_ratio))endsety_mod=5e8p_ratio=0.25derivegenzonebricksize2,2,2modelelasticpropbulk=b_modshear=s_modprintzonepropbulkprintzonepropshear2-7控制循环Newdefxxxsum=0prod=1loopn(1,10)sum=sum+nprod=prod*nend_loopendxxxprintsum,prodnewgenzonebrickp0(0,0,0)p1(-10,0,0)p2(0,10,0)p3(0,0,-10)modelelasplotsetrotation0045plotblockgroupdefinstallpnt=zone_headloopwhilepnt#nullz_depth=-z_zcen(pnt)y_mod=y_zero+cc*sqrt(z_depth)z_prop(pnt,’shear’)=y_mod/(2.0*(1.0+p_ratio))z_prop(pnt,’bulk’)=y_mod/(3.0*(1.0-2.0*p_ratio))pnt=zone_next(pnt)end_loopendsetp_ratio=0.25y_zero=1e7cc=1e8install2-8拆分命令行new;exampleofasumofmanythingsdeflong_sumtemp=v1+v2+v3+v4+v5+v6+v7+v8+v9+v10long_sum=temp+v11+v12+v13+v14+v15end2-9变量类型newdefhaveoneaa=2bb=3.4cc=’Haveaniceday’dd=aa*bbee=cc+’,oldchap’endhaveoneprintfish2-10IF条件语句newdefabcifxx0thenabc=1000elseabc=-1000end_ifendsetxx=10printabcsetxx=0printabc2-11索单元自动生成newgenzonebricksize1035plotsetrotation0045plotblockgroupdefplace_cablesloopn(1,5)z_d=float(n)-0.5commandselcablebeg0.0,1.5,z_dend7.0,1.5,z_dnseg7end_commandend_loopendplace_cablesplotgridselgeomrednewgenzonebricksize1035plotsetrotation15060plotblockgroupmodmohrpropbulk1e8shear.3e8fric35propcoh1e3tens1e3inidens1000setgrav0,0,-10fixxyzrangez-.1.1fixyrangey-.1.1fixyrangey2.93.1fixxrangex-.1.1fixxrangex9.910.1setlargehistunbalsolvesavecab_str.savinixdis0ydis0zdis0histgpxdisp0,1,5defplace_cablesloopn(1,5)z_d=5.5-float(n)z_t=z_d+0.5z_b=z_d-0.5commandfreexrangex-.1,.1zz_bz_tsolveselcablebeg0.0,0.5,z_dend7.0,0.5,z_dnseg7selcablebeg0.0,1.5,z_dend7.0,1.5,z_dnseg7selcablebeg0.0,2.5,z_dend7.0,2.5,z_dnseg7selcablepropemod2e10ytension1e8xcarea1.0&gr_k2e10gr_coh1e10gr_per1.0end_commandend_loopendplace_cablessavecab_end.savplotsketchselcableforcered2-12圆形隧道开挖模拟计算;建立模型genzonradcylp0000p1600p2010p3006&size4284dim3333rat1111.2groupoutsiderockgenzonecshellp0000p1300p2010p3003&size1284dim2.72.72.72.7rat1111groupconcretlinerfillgroupinsiderockgenzonreflectdip90dd90orig000genzonreflectdip0dd0ori000genzonbrickp0006p1606p2016p30013size426groupoutsiderock1genzonbrickp000-12p160-12p201-12p300-6size425groupoutsiderock2genzonbrickp0600p12100p2610p3606size1024groupoutsiderock3genzonreflectdip0dd0orig000rangegroupoutsiderock3genzonbrickp0606p12106p2616p36013size1026groupoutsiderock4genzonbrickp060-12p1210-12p261-12p360-6size1025groupoutsiderock5genzonreflectdip90dd90orig000rangex-0.16.1z6.113.1genzonreflectdip90dd90orig000rangex-0.16.1z-6.1-12.1genzonreflectdip90dd90orig000rangex6.121.1z-12.113.1;绘制模型图plotblockgroupplotaddaxesred;plotsetrotation0045用于显示三维模型;设置重力setgravity00-10;给定边界条件fixzrangez-12.01,-11.99fixxrangex-21.01,-20.99fixxrangex20.99,21.01fixyrangey-0.010.01fixyrangey0.99,1.01;求解自重应力场modelmohrinidensity1800;围岩的密度propbulk=1.47e8shear=5.6e7fric=20coh=5.0e4tension=1.0e4;体积、剪切、摩擦角、凝聚力、抗拉强度setmechratio=1e-4solvesaveGravsol.savplotcontzdispoutlonplotcontszz;毛洞开挖计算initialxdisp=0ydisp=0zdisp=0modelnullrangegroupinsiderockanygroupconcretlineranyplotblockgroupplotaddaxesredsetmechratio=5e-4solvesaveKaiwsol.savplotcontzdispplotcontsdispplotcontszzplotcontxzz;模筑衬砌计算modelelasrangegroupconcretlineranyplotblockgroupplotaddaxesredinidensity2500rangegroupconcretlinerany;衬砌混凝土的密度propbulk=16.67e9,shear=12.5e9rangegroupconcretlinerany;衬砌混凝土的体积弹模、剪切弹模setmechratio=1e-4solvesavezhihusol.savplotcontzdispplotcontsdispplotcontszzplotcontxzz;完成计算分析