FLAC3D模拟实例循环开挖与支护nres ini.savset geometry=0.001ini ydis0ini xdis0ini zdis0ini yvel0ini xvel0ini zvel0m mprop bulk 4.0e9shear 2.5e9fri32coh 2.0e6& range grou diban-shayan;prop bulk 1.8e9shear 1.2e9fri25coh 1.0e6& range grou diban-niyan any grou hangdao any;prop bulk 1.2e9shear0.8e9fri22coh0.8e6& range grou diban-gentuyan;prop bulk 1.9e9shear 1.3e9fri24coh 1.0e6& range grou diban-tniyan;prop bulk0.7e9shear0.8e9fri21coh0.7e6& range grou mc any grou gzm any;prop bulk 3.0e9shear 2.5e9fri30coh 1.8e6& range grou dingban-fenshayan;prop bulk 1.5e9shear 1.2e9fri25coh 1.1e6& range grou dingban-niyan;prop bulk 3.5e9shear 2.5e9fri34coh 1.4e6& range grou dingban-shayan;添加接触面gen separate gzminterface1wrap mc gzminterface2wrap dingban-fenshayan gzm interface1prop kn20e9ks10e9tens1e9 interface2prop kn20e9ks10e9tens1e9set mech ratio=5e-4def excav_mcloop n(excav_p,excav_p_z+cut_liang);每次开挖量cut_0=excav_pcut_1=excav_p+cut_liang;开挖commandm null range grou gzm z cut_0cut_1step100end_commandn=excav_p+cut_Liangexcav_p=excav_p+cut_Liang;条件判断保存文件,这里判断条件必须和cut_liang对应上，否则不能得到想要的文件。
if cut_1=40thencommandsave kwmc40m.savend_commandend_ifend_loopend;每次需要修改set cut_liang=5excav_p=0excav_p_z=80excav_mcsolvepl block groupl add axes redpl set rotation9000save kwmc2.savcall kwhdno.txt;n;res kwmc.savini ydis0ini xdis0ini zdis0ini yvel0ini xvel0ini zvel0set largeset mech ratio=5e-4pl cont ydis shad on outline onhis gp ydis901740;40m处巷道巷道顶部中点his gp xdis77.7514.7540;40m处巷道巷道左帮中点his gp xdis82.2514.7540;40m处巷道巷道右帮中点def excav_hangdaoloop n(excav_p,excav_p_z+cut_liang);每次开挖量cut_0=excav_pcut_1=excav_p+cut_liang;开挖commandm null range grou hangdao z cut_0cut_1step50end_commandn=excav_p+cut_Liangexcav_p=excav_p+cut_Liang;条件判断保存文件,这里判断条件必须和cut_liang对应上，否则不能得到想要的文件。
if cut_1=40thencommandsave kwhdno40m.savend_commandend_ifend_loopendset cut_liang=8excav_p=0excav_p_z=80excav_hangdaosolvesave kwhdno2.savcall kwhdzh.txt;n;res kwmc2.sav;pl set rot901800;从后视图看，零点面对着;pl set rot9000;从前视图看，终点面对着pl set rot907515pl cont syy range grou hangdaopl add sel geom node off fill on black blue;显示支护单元;pl add sel geom cable blue;sel set damp combined;sel liner prop slide on;添加接触面;gen separate2;inter1wrap12;inter2wrap32;inter1face range plane norm010ori011.08150dist 0.1;inter2face range plane norm010ori010.71930dist 0.1;interface1prop kn1e10ks1e10;interface2prop kn1e10ks1e10set largeini ydis0ini xdis0ini zdis0ini yvel0ini xvel0ini zvel0set largeset mech ratio=5e-4his gp ydis901740;40m处巷道巷道顶部中点his gp xdis77.7514.7540;40m处巷道巷道左帮中点his gp xdis82.2514.7540;40m处巷道巷道右帮中点;监控锚杆索命令;his sel cable force x y z;his sel cable stress x y z;his sel cable grout stress x y z;监控锚杆索构件端部水泥浆应力(端部);his sel cable grout slip x y z;滑动状态（端部）（0未屈服1正在屈服2曾屈服）;his sel cable grout disp x y z;监控锚杆索构件端部在水泥浆中的位移;———————————————————————————————————;—————————————————————————————————————————————;锚杆（索）、初撑参数及剖面位置参数;; ;—————————————————————————————————————————————;———————————————————————-————————————def maogan_canshucable_seg=10;构件数emod_e=180e9;弹性模量xcarea_x=3.46185e-4;锚杆（索）横截面积ytens_y=25e4;抗拉强度gr_per_gp=0.11;水泥浆外圈长度gr_k_gk=6e6;单位长度上水泥浆刚度gr_fric_gf=30.0;水泥浆摩擦角gr_coh_gc=1.0e5;单位长度上水泥浆的粘结力;ycomp_y=;抗压强度enddef maosuo_canshus_cable_seg=12;构件数emod__e=200e9;弹性模量xcarea__x=3.46185e-4;锚杆（索）横截面积ytens__y=35e4;抗拉强度gr_per__gp=0.11;水泥浆外圈长度gr_k__gk=6e6;单位长度上水泥浆刚度gr_fric__gf=30;水泥浆摩擦角;gr_coh__gc=;单位长度上水泥浆的粘结力;ycomp__y=;抗压强度endmaogan_canshu maosuo_canshudef maogan(suo)_p ;锚杆位置;maogan1xp1=77.75yp1=14.0xp1_2=72.25yp1_2=14.0;maogan2xp2=77.8686yp2=15.4709xp2_2=75.5004yp2_2=16.2719;maogan3xp3=79.0593yp3=16.7939xp3_2=78.0140yp3_2=19.0649;maogan4xp4=80.9407yp4=16.7939xp4_2=81.9680 yp4_2=19.0649 ;maogan5xp5=82.1314yp5=15.4709xp5_2=84.4996 yp5_2=16.2719 ;maogan6xp6=82.25yp6=14.0xp6_2=84.75yp6_2=14.0;;锚索位置;maosuo1xps1=77.75yps1=13.3536 xps1_2=71.25 yps1_2=13.3536 ;maosuo2xps2=77.75yps2=14.75xps2_2=71.25 yps2_2=14.75;maosuo3xps3=78.3554 yps3=16.2638 xps3_2=73.5264 yps3_2=20.6368 ;maosuo4xps4=80.0yps4=17.0xps4_2=80.0yps4_2=23.5;maosuo5xps5=81.6646 yps5=16.2638 xps5_2=86.4736 yps5_2=20.63668 ;maosuo6xps6=82.25yps6=14.75xps6_2=88.75yps6_2=14.75;maosuo7xps7=82.25yps7=13.3336xps7_2=88.75yps7_2=13.3536endmaogan(suo)_p;;;开挖量、范围、支护排距控制等参数def excav_canshucut_x=80.0cut_y=14.75cut_r=2.25;开挖巷道半圆拱的中心坐标及半径cut_zhiqiang=1.5;直墙壁高度;cut_0=0.0;开挖起始位置cut_Liang=10;一次开挖量cut_1=cut_0+cut_liangexcav_p_z=80.0;开挖终点位置;paiju=1.0;支护断面间的排距first_p=0.0;开挖起始段第一排支护的位置，参数first_p需要根据前一开挖段剩下的排距大小调整endexcav_canshugen separate hangdao;添加初衬liner或shell单元前需要产生一分离单元面;开挖与支护def excavandzhihuexcav_p=cut_0loop n(excav_p,excav_p_z+2*cut_liang);（2是参数,防止循环跳出，未装锚杆索，下个循环中亦是如此）;每次开挖量cut_0=excav_pcut_1=excav_p+cut_liang;;下面的参数是控制添加初衬单元范围，1.1为参数，可更改。