FLAC命令流

ini pp 0 j=0fix pp j=0这样设置后是透水边界了底边界不透水,你可以不用管它了,默认为不透水边界我现在在三维中模拟基坑开挖，x(0~42) y(0~30) z(-40~0),初始水位在－3米处，我的初始水位命令这样写对不对啊？ini pp -3e4 grad 0 0 -1e4 range x 0 42 y 0 30 z -40 -3range -3.1 -2.9;;在开挖中需要降水，降至－7米，需开挖的区域为x(0~8) y(0~30)fix pp range x 8 42 y o 30 z -3.1 -2.9fix pp 0 range x 0 8 y 0 30 z -7.1 -6.9 pl contour pp可以显示孔隙水压力场，一般就认为就是渗流场吧！set flow on 即选择是否要进行渗流分析，其前提是要set fluid on方可以！off就是相反了。

但是，set fluid on以后不一定进行set flow on的计算，这点我的理解（有待证实）是，仅有set fluid on的话，就是把里面的pp当作是静水压力来处理的，而当set flow on 打开以后，就可以进行真正意义上的渗流分析了！再一个，set flow on的一个前提必须设定perm值（渗透系数），否则程序无法继续进行！而set fluid on就没有那么严格，只要是config fluid 以及设定模型以及fmod以后即可以开始计算了。

set mech on就表示把力学计算这部分打开，在默认的情况下，是打开的。

在流固耦合有一节中讲到，当只要使系统产生水压力的分布，而不需要由于机械变形（mechanical deformation）而产生的pp的时候，就必须要set mech off，也就是说这个时候只有fluid的计算，而没有mechanical 的计算了。

帮助的第ftd138.pdfbulk modulus, K 390 MPa 体积模量（土质比较硬）shear modulus, G 280 MPa 剪切模量soil dry density, ρd 1200 kg/m3 土的密度water density, ρw 1000 kg/m3 水密度wall density, ρwal 1500 kg/m3 档土墙的密度permeability, k 10−12 m2/Pa-s 渗透系数porosity, n 0.3 孔隙率fluid bulk modulus, Kf 2.0 GPa 流体的体积模量;-----以下是命令流config fluid ;设置流体算法gen zone brick p1 12 0 0 p2 0 12 0 p3 0 0 12 size 12 12 12 rat 1 1 1group soil ;设置土组group excavate range x 0 4 y 0 4 z 0 5 ;设置开挖范围组group wal1 range x 4 5 y 0 5 z 0 7 ; 设置档土墙1 组group wal2 range x 0 4 y 4 5 z 0 7 ;设置档土墙2 组group wall range group wal1 any group wal2 any ;档土墙1 和2 合并为一个档土墙组; --- fluid flow model --- 流体模型model fl_iso ;各项同性流体prop perm 1e-12 poro 0.3 ;设置渗透系数和？？比or ？？率ini fdensity 1e3 ; 初始化水密度ini fmod 2.0e9 ftens -1e-3 ; 初始化水的体积模量及？？model fl_null range group excavate ; 开挖部分没有水流model fl_null range group wall ;档土墙不渗水jini pp 0 grad 0 0 1e4 ;初始化孔隙水压力fix pp range z -0.1 0.1 ;约束应该是整个模型的底部吧？fix pp range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1 ;开挖部分也给约束住了，而其他的地方不管，为什么要约束开挖部分？model elas ; 弹塑性prop bul 3.9e6 shea 2.8e6 ;设置体积模量和剪切模量，这两个必须model null range group excavate ;力学模型的开挖ini density 1.2e3 ;初始化密度，应该是土的密度ini density 1.5e3 range group wall ;初始化墙（档土墙）的密度fix x range x -.1 .1 ;施加约束，据说差分法约束0 的时候，必须取一个-0.1<0<0.1,范围fix x range x 11.9 12.1 ;施加约束，或者叫做边界条件fix y range y -.1 .1 ;施加约束fix y range y 11.9 12.1fix z range z 11.9 12.1ini szz 0 grad 0 0 -1.5e4 ;z方向的应力大一些ini sxx 0 grad 0 0 -1.2e4 ;水平方向的，小一些ini syy 0 grad 0 0 -1.2e4apply nstress 0 grad 0 0 -1.2e4 range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0apply nstress 0 grad 0 0 -1.2e4 range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0apply nstress -7.5e4 range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1; --- setting --- 设置set gravity 0 0 10 ;设置重力加速度solve force 1 ; check initial equilibrium 求解平衡; --- excavation --- 开挖set fluid off ;关闭水apply nstress 0 grad 0 0 -1.e4 range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0apply nstress 0 grad 0 0 -1.e4 range x 3.9 4.1 y 0.0 4.0 z 0.0 5.０apply nstress -5.e4 range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1solve ;force 1 ;求解save exc1.sav ; 保存; --- drainage --- 排水apply remove nstress ;撤掉刚才的压力def relaxsetup ;定义一个释放函数的参数step0 = stependrelaxsetup ;调用def relax ;定义一个释放函数rstep = step - step0if rstep < ncyc thenrelax=1.0-(float(rstep)/float(ncyc))elserelax=0.0endifendset ncyc=1000；先赋值，然后调用apply nstress 0 grad 0 0 -1e4 his relax range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0 apply nstress 0 grad 0 0 -1e4 his relax range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0 apply nstress -5e4 his relax range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1cyc ncycsolve; --- percolation --- 渗透set fluid on ;设置水算法Fluentfix pp 0 range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1 ;施加孔隙水压力cyc 9000 ;循环save exc3.sav ;保存你是想做渗流分析还是静力分析渗流分析的话用config fluid做流固耦合分析把手册中几个例子搞明白应该收获很大提醒生成初始应力场时水平应力＝k0×竖向有效应力＋水应力如果仅仅是静力分析只要设定稳定水位手动干湿密度就行与Interface相关的命令总结（望续）最近用到这个了，因为初学尝尽苦头，以前整理的贴出来，望有别的还可以跟帖，simwe9j2g2K0`$^.p;j9@仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM7v6Q8E*O1M%R4~ 1）interface后面可以使用range关键字，用来限定其范围；range id用来设定界面元或节点的开始和最后的标识号码；2l5M:Y!T$asimwe2）interface命令用来定义两个或多个子网格之间的界面，这些界面可以滑动、分离。

《FLA C/FLA C 3D基础与‎工程实例》全部命令流‎1gen zone bri p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 &p4 15 15 0 p5 0 15 15 p6 15 0 10 p7 20 20 20 &size 10 10 10 rat 1.0 0.9 1.1 group‎brick‎_1gen zone bri p0 20 0 0 p1 add 10 0 0 p2 add 0 20 0 p3 add 0 0 15 & size 10 10 10 rat 1.0 0.9 1.1 group‎brick‎_2gen zone bri p0 40 0 0 edge 10 size 10 10 10 rat 1.0 0.9 1.1 group‎brick‎_3 plot surnewgen zon bri size 3 3 3model‎elasprop bulk 3e8 shear‎1e8ini dens 2000fix z ran z -.1 .1fix x ran x -.1 .1fix x ran x 2.9 3.1fix y ran y -.1 .1fix y ran y 2.9 3.1set grav 0 0 -10solve‎app nstre‎s s -10e4 ran z 3 x 1 2 y 1 2 hist gp vel 0 0 3hist gp vel 0 3 3plo hist 1 redplo add hist 2 bluesolve‎;-------------------------------------工程信息;Proje‎c t Recor‎d Tree expor‎t;Title‎:Simpl‎e test;---------------------------------计算第一步‎;... STATE‎: STATE‎1 ....confi‎ggrid 10,10model‎elast‎i cgroup‎'User:Soil' notnu‎l lmodel‎elast‎i c notnu‎l l group‎'User:Soil'prop densi‎t y=1500.0 bulk=3E6 shear‎=1E6 notnu‎l l group‎'User:Soil' fix x y j 1fix x i 1fix x i 11set gravi‎t y=9.81histo‎r y 999 unbal‎a nced‎solve‎save state‎1.sav;----------------------------------计算第二步‎;... STATE‎: STATE‎2 ....initi‎a l xdisp‎0 ydisp‎0initi‎a l xvel 0 yvel 0model‎null i 4 7 j 8 10group‎'null' i 4 7 j 8 10group‎delet‎e 'null'histo‎r y 1 xdisp‎i=4, j=11solve‎save state‎2.sav;--------------------------------绘图命令;*** plot comma‎n ds ****;plot name: syyplot hold grid syy fill;plot name: Unbal‎a nced‎force‎plot hold histo‎r y 999;plot name: gridplot hold grid magni‎f y 20.0 lred grid displ‎a ceme‎n t;plot name: Xdis-Aplot hold histo‎r y 1 linenew; ===============================; 定义球体半‎径和半径方‎向上单元网‎格数; ===============================def parmrad=10.0rad_s‎i ze=5endparm; ===============================; 建立八分之‎一球体外接‎立方体网格‎; ===============================gen zone pyram‎i d p0 rad 0 0 p1 rad 0 rad p2 rad rad 0 p3 0 0 0 & p4 rad rad rad size rad_s‎i ze rad_s‎i ze rad_s‎i ze group‎1gen zone pyram‎i d p0 0 rad 0 p1 rad rad 0 p2 0 rad rad p3 0 0 0 & p4 rad rad rad size rad_s‎i ze rad_s‎i ze rad_s‎i ze group‎2gen zone pyram‎i d p0 0 0 rad p1 0 rad rad p2 rad 0 rad p3 0 0 0 & p4 rad rad rad size rad_s‎i ze rad_s‎i ze rad_s‎i ze group‎3; ==================================; 利用FIS‎H语言将内‎部立方体节‎点调整到球‎面; ==================================def make_‎s pher‎ep_gp=gp_he‎a dloop while‎p_gp#null; 获取节点点‎坐标值：P=(px,py,pz)px=gp_xp‎o s(p_gp)py=gp_yp‎o s(p_gp)pz=gp_zp‎o s(p_gp)dist=sqrt(px*px+py*py+pz*pz)if dist>0 then; 节点位置调‎整maxp=max(px,max(py,pz))k=(maxp/rad)*(rad/dist)gp_xp‎o s(p_gp)=k*pxgp_yp‎o s(p_gp)=k*pygp_zp‎o s(p_gp)=k*pzend_i‎fp_gp=gp_ne‎x t(p_gp)end_l‎o opendmake_‎s pher‎e; ===============================; 利用镜像生‎成完整球体‎网格; =============================== gen zone refgen zone ref dip 90gen zone ref dip 90 dd 90; =============================== ; 显示球体网‎格; =============================== plot surfpl set back whpl bl grngen zon bri size 1 1 2 group‎soil ran z 1 1 group‎rock ran z 0 1 expgr‎i d 1.flac3‎dgen zon bri size 3 3 3model‎mohrprop bu 3e6 sh 1e6 coh 10e3 fric 15 fix z ran z -.1 .1fix x ran x -.1 .1fix x ran x 2.9 3.1fix y ran y -.1 .1fix y ran y 2.9 3.1ini dens 2000hist unbal‎set grav 10solve‎elast‎i csave 6-1.savrest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstre‎s s -100e3‎ran z 2.9 3.1 x 1 2 y 1 2 solve‎save 6-2.savrest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstre‎s s -100e3‎ran z 2.9 3.1 x 1 2 y 1 2 hist id=2 gp zdis 1 1 3hist id=3 gp zdis 1 1 2hist id=4 gp xdis 1 1 3hist id=5 gp xdis 1 1 3hist id=6 zone szz 1 1 3hist id=7 zone szz 1.5 1.5 2.5hist id=8 zone sxz 1.5 1.5 2.5solve‎save 6-3.savrest 6-3.savset log onset logfi‎l e 6-2.log print‎zone stres‎s print‎gp disset log offrest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstre‎s s -100e3‎ran z 2.9 3.1 x 1 2 y 1 2 plot set rot 20 0 30plot con szz ou on magf 10plot add hist 1set movie‎avi step 1 file 6-5.avimovie‎start‎solve‎movie‎finis‎hngen zon bri size 1 1 2 model‎elasprop bulk 3e7 shear‎1e7 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solve‎plo con szngen zon bri size 1 1 2model‎mohrprop bulk 3e7 shear‎1e7 c 1e10 f 15 tensi‎o n 1e10 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solve‎prop bulk 3e7 shear‎1e7 c 10e3 f 15 ten 0solve‎plo con szngen zone brick‎size 1 1 2model‎mohrprop bulk 3e7 shear‎1e7 coh 10e3 fri 15 ten 0 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solve‎elasplo con sznewgen zone brick‎size 1 1 2model‎mohrprop bulk 3e7 shear‎1e7 coh 10e3 fri 15 ten 0 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000ini szz -40e3 grad 0 0 20e3 ran z 0 2ini syy -20e3 grad 0 0 10e3 ran z 0 2ini sxx -20e3 grad 0 0 10e3 ran z 0 2set grav 0 0 -10solve‎plo con szngen zon bri size 1 1 2model‎mprop bulk 3e7 shear‎1e7 c 10e10‎f 15 ten 1e10 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000 ran z 0 1ini dens 1500 ran z 1 2ini szz -35e3 grad 0 0 20e3 ran z 0 1ini syy -17.5e3 grad 0 0 10e3 ran z 0 1ini sxx -17.5e3 grad 0 0 10e3 ran z 0 1ini szz -15e3 grad 0 0 15e3 ran z 1 2ini syy -7.5e3 grad 0 0 7.5e3 ran z 1 2ini sxx -7.5e3 grad 0 0 7.5e3 ran z 1 2ini pp 10e3 grad 0 0 -10e3 ran z 0 1set grav 0 0 -10solve‎plo con szngen zon bri size 1 1 2model‎mprop bulk 3e7 shear‎1e7 c 10e10‎f 15 ten 1e10 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000 ran z 0 2ini szz -50e3 grad 0 0 20e3 ran z 0 1ini syy -30e3 grad 0 0 10e3 ran z 0 1ini sxx -30e3 grad 0 0 10e3 ran z 0 1ini pp 30e3 grad 0 0 -10e3 ran z 0 2app nstre‎s s -10e3 ran z 2set grav 0 0 -10solve‎plo con sznewgen zone brick‎p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 & p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & size 6 6 10model‎elaspro bulk 10e10‎she 10e10‎ini den 2500apply‎sxx -1e9 grad 0 0 1.11111‎11e7 range‎x -.1 .1 apply‎sxx -1e9 grad 0 0 6.66666‎66e6 range‎x 59.9 60.1 apply‎syy -1e9 grad 0 0 8.33333‎33e6 range‎y -.1 .1 apply‎syy -1e9 grad 0 0 8.33333‎33e6 range‎y 59.9 60.1 apply‎szz -1e8 grad 0 0 8.33333‎33e5 ran z 0 120set grav 0 0 -10step 30000‎ini xdisp‎0 ydisp‎0 zdisp‎0ini xvel 0 yvel 0 zvel 0plo cont szznewgen zone brick‎p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 & p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & size 6 6 10model‎elaspro bulk 10e10‎she 10e10‎ini den 2500ini sxx -1e9 grad 0 0 1.11111‎11e7 range‎x -.1 .1ini sxx -1e9 grad 0 0 6.66666‎66e6 range‎x 59.9 60.1ini syy -1e9 grad 0 0 8.33333‎33e6 range‎y -.1 .1ini syy -1e9 grad 0 0 8.33333‎33e6 range‎y 59.9 60.1ini szz -1e8 ran z -.1 .1fix x y z ran z -.1 .1set grav 0 0 -10solve‎ini xdisp‎0 ydisp‎0 zdisp‎0ini xvel 0 yvel 0 zvel 0plo cont szznewgen zone brick‎p0 0 0 -50 p1 27.5 0 -50 p2 0 5 -50 p3 0 0 -10 size 8 1 10 group‎claygen zone brick‎p0 27.5 0 -50 p1 100 0 -50 p2 27.5 5 -50 p3 27.5 0 -10 ratio‎1.1 1 1 size 12 1 10 group‎claygen zone brick‎p0 0 0 -10 p1 27.5 0 -10 p2 0 5 -10 p3 0 0 0 ratio‎1 1 0.8 size 8 1 4 group‎soilgen zone brick‎p0 27.5 0 -10 p1 100 0 -10 p2 27.5 5 -10 p3 27.5 0 0 ratio‎1.1 1 0.8 size 12 1 4 group‎soilgen zone brick‎p0 0 0 0 p1 27.5 0 0 p2 0 5 0 p3 0 0 5 p4 27.5 5 0 &p5 0 5 5 p6 20 0 5 p7 20 5 5 size 8 1 5 group‎damfix x y z ran z -49.9 -50.1fix x ran x -.1 .1fix x ran x 99.9 100.1fix ymodel‎mohr ran z -50 0model‎null ran z 0 5prop bulk 7.8e6 shear‎3.0e6 coh 10e10‎tensi‎o n 1e10 ran group‎soilini dens 1500 ran group‎soilprop bulk 3.91e6 shear‎1.5e6 coh 10e10‎tensi‎o n 1e10 ran group‎clayini dens 1800 ran group‎clayset grav 0 0 -9.8hist id=1 unbal‎solve‎prop bulk 7.8e6 shear‎3.0e6 coh 10e3 fric 15 ran group‎soilprop bulk 3.91e6 shear‎1.5e6 coh 20e3 fric 20 ran group‎claysolve‎save elast‎i c.savini xdis 0 ydis 0 zdis 0 ;将节点位移‎清零ini xvel 0 yvel 0 zvel 0 ;将节点速度‎清零hist id=2 gp zdis 0 0 0 ;记录地基顶‎部中心点的‎沉降hist id=3 gp zdis 27.5 0 0 ;记录路基坡‎脚处的沉降‎hist id=4 gp xdis 27.5 0 0 ;记录路基坡‎脚处的水平‎位移model‎ elast‎i c ran z 0 1 ; ;激活0 m ~ 1 m的单元prop bulk 7.8e6 shear‎3.0e6 ran z 0 1ini dens 1500 ran z 0 1solve‎;按软件默认‎精度求解save fill-1.savmodel‎elast‎i c ran z 1 2prop bulk 7.8e6 shear‎3.0e6 ran z 1 2ini dens 1500 ran z 1 2solve‎save fill-2.savmodel‎elast‎i c ran z 2 3prop bulk 7.8e6 shear‎3.0e6 ran z 2 3ini dens 1500 ran z 2 3solve‎save fill-3.savmodel‎elast‎i c ran z 3 4prop bulk 7.8e6 shear‎3.0e6 ran z 3 4ini dens 1500 ran z 3 4solve‎save fill-4.savmodel‎elast‎i c ran z 4 5prop bulk 7.8e6 shear‎3.0e6 ran z 4 5ini dens 1500 ran z 4 5solve‎save fill-5.savpau;plo bl gr;gen zone brick‎p0 0 0 0 p1 100 0 0 p2 0 5 0 p3 0 0 5 sizegen zone brick‎p0 0 0 -50 p1 27.5 0 -50 p2 0 5 -50 p3 0 0 -10 size 8 1 10 group‎clay gen zone brick‎p0 27.5 0 -50 p1 100 0 -50 p2 27.5 5 -50 p3 27.5 0 -10 ratio‎1.1 1 1 size 12 1 10 group‎claygen zone brick‎p0 0 0 -10 p1 27.5 0 -10 p2 0 5 -10 p3 0 0 0 ratio‎1 1 0.8 size 8 1 4 group‎soilgen zone brick‎p0 27.5 0 -10 p1 100 0 -10 p2 27.5 5 -10 p3 27.5 0 0 ratio‎1.1 1 0.8 size 12 1 4 group‎soilgen zone brick‎p0 0 0 0 p1 27.5 0 0 p2 0 5 0 p3 0 0 5 p4 27.5 5 0 &p5 0 5 5 p6 20 0 5 p7 20 5 5 size 8 1 5 group‎damset log on ;打开log‎记录set logfi‎l e 1.log ;设置记录文‎件名为：1.logresto‎re fill-1.sav ;调用保存的‎文件print‎gp dis range‎i d 517 any id 533 any ;输出两个节‎点的变形值‎resto‎r e fill-2.savprint‎gp dis range‎id 517 any id 533 anyresto‎r e fill-3.savprint‎gp dis range‎id 517 any id 533 anyresto‎r e fill-4.savprint‎gp dis range‎id 517 any id 533 anyresto‎r e fill-5.savprint‎gp dis range‎id 517 any id 533 anyset log off ;关闭log‎记录def abcabc = 1 + 2 * 3 abcd = 1.0 / 2.0 endabcprint‎fishdef abcif aa < 0 thenabc = 0.0 elseabc = 2.0 * aa endif‎endabcdef abcloop aa (1, 2.5) comma‎n dprint‎aaendco‎m mand‎ endlo‎o pendnewgen zon bri size 3 3 3model‎elast‎i cprop bu 3e7 sh 1e7ini dens 2000fix x y z ran z -.1 .1fix x ran x -.1 .1fix x ran x 2.9 3.1fix y ran y -.1 .1fix y ran y 2.9 3.1set grav 10solve‎ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0save 8-4.sav8-5rest 8-4.savdef E_mod‎i fyp_z = zone_‎h eadd_k = 704d_n = 0.38d_pa = 10132‎5.0 ;//标准大气压‎loop while‎p_z # nullsigma‎_3 = -1.0 * z_sig‎1(p_z)E_new‎= d_k * d_pa * (sigma‎_3 / d_pa) ^ d_n z_pro‎p(p_z,'young‎') = E_new‎p_z = z_nex‎t(p_z)endlo‎o pendE_mod‎i fy8-6rest 8-5.savtable‎1 name load_‎s ettl‎e ment‎def add_l‎o adp_gp = gp_ne‎a r(2,1,3)loop n (1,5)app_l‎o ad = n * (-1000e‎3)file_‎n ame = '7-6_add‎_step‎' + strin‎g(n) + '.sav' comma‎n dapp nstre‎s s app_l‎o ad ran z 2.9 3.1 x 1 2 y 1 2 s olve‎s ave file_‎n ameendco‎m mand‎xtabl‎e(1,n) = -1.0 * app_l‎o adytabl‎e(1,n) = gp_zd‎i sp(p_gp)endlo‎o pendadd_l‎o adsave 8-6.sav8-7rest 8-6.savdef find_‎m ax_d‎i spp_gp = gp_he‎a dmaxdi‎s p_va‎l ue = 0.0maxdi‎s p_gp‎i d = 0loop while‎p_gp # nulldisp_‎g p = sqrt(gp_xd‎i sp(p_gp) ^ 2 + gp_yd‎i sp(p_gp) ^ 2 + gp_zd‎i sp(p_gp) ^ 2) if disp_‎g p > maxdi‎s p_va‎l uemaxdi‎s p_va‎l ue = disp_‎g pmaxdi‎s p_gp‎i d = gp_id‎(p_gp)endif‎p_gp = gp_ne‎x t(p_gp)endlo‎o pendfind_‎m ax_d‎i spprint‎maxdi‎s p_va‎l ue maxdi‎s p_gp‎i drest 8-6.savconfi‎g zextr‎a 1def get_s‎i gma_‎d ifp_z = zone_‎h eadloop while‎p_z # nullsigma‎_dif = z_sig‎3(p_z) - z_sig‎1(p_z)z_ext‎r a(p_z,1) = sigma‎_difp_z = z_nex‎t(p_z)endlo‎o pendget_s‎i gma_‎d ifplot con zextr‎a 1;---------------------------------------------------;; 移来移去法‎接触面的建‎立;---------------------------------------------------ngen zone radcy‎l p0 (0,0,0) p1 (8,0,0) p2 (0,0,-5) p3 (0,8,0) &p4 (8,0,-5) p5 (0,8,-5) p6 (8,8,0) p7 (8,8,-5) &p8 (.3,0,0) p9 (0,.3,0) p10 (.3,0,-5) p11 (0,.3,-5) & size 3 10 6 15 ratio‎1 1 1 1.15gen zone radcy‎l p0 (0,0,-5) p1 (8,0,-5) p2 (0,0,-8) p3 (0,8,-5) &p4 (8,0,-8) p5 (0,8,-8) p6 (8,8,-5) p7 (8,8,-8) &p8 (.3,0,-5) p9 (0,.3,-5) p10 (.3,0,-8) p11 (0,.3,-8) & size 3 6 6 15 ratio‎1 1 1 1.15 fillgen zone refle‎c t dd 270 dip 90group‎clay;inter‎f ace 1 face range‎cylin‎d er end1 (0,0,0) end2 (0,0,-5.1) radiu‎s .31 &cylin‎d er end1 (0,0,0) end2 (0,0,-5.1) radiu‎s .29 notinter‎f ace 2 face range‎cylin‎d er end1 (0,0,-4.9) end2 (0,0,-5.1) radiu‎s .31;pause‎gen zone cyl p0 (0,0,6) p1 (.3,0,6) p2 (0,0,1) p3 (0,.3,6) &p4 (.3,0,1) p5 (0,.3,1) &size 3 10 6gen zone cyl p0 (0,0,6.1) p1 (.3,0,6.1) p2 (0,0,6) p3 (0,.3,6.1) &p4 (.3,0,6) p5 (0,.3,6) &size 3 1 6gen zone refle‎c t dd 270 dip 90 range‎z 1 6.1group‎pile range‎z 1 6.1pause‎ini z add -6.0 range‎group‎pilesave pile_‎g eom.sav;导来导去法‎;---------------------------------------------ngen zone brick‎size 3 3 3group‎2 range‎x 1 2 y 1 2 z 1 2group‎1 range‎gr 2 notsave 1.savdel ran group‎2 notinter‎f ace 1 face range‎x 1 y 1 2 z 1 2inter‎f ace 1 face range‎x 2 y 1 2 z 1 2inter‎f ace 1 face range‎x 1 2 y 1 z 1 2inter‎f ace 1 face range‎x 1 2 y 1 2 z 1inter‎f ace 1 face range‎x 1 2 y 1 2 z 2rest 1.savdel ran group‎2expgr‎i d 1.fac3d‎impgr‎i d 1.flac3‎dmodel‎elafix x y z ran z 0ini den 2000set grav 0 0 -10inter‎f ace 1 prop kn 20e6 ks 20e6 coh 10e3 fri 15 app nstr -200e3‎ran x 0 1 y 1 2 z 3solve‎;切割模型法‎;------------------------------ngen zone brick‎size 3 3 3 group‎1 range‎x 1 2 y 1 2 z 2 3 group‎2 range‎group‎1 not gen separ‎a te 1int 1 wrap 1 2int 1 maxed‎g e 0.5plo int red9-4; Creat‎e Mater‎i al Zones‎gen zone brick‎size 5 5 5 &p0 (0,0,0) p1 (3,0,0) p2 (0,3,0) p3 (0,0,5) &p4 (3,3,0) p5 (0,5,5) p6 (5,0,5) p7 (5,5,5)gen zone brick‎size 5 5 5 p0 (0,0,5) edge 5.0group‎Mater‎i al; Creat‎e Bin Zones‎gen zone brick‎size 1 5 5 &p0 (3,0,0) p1 add (3,0,0) p2 add (0,3,0) &p3 add (2,0,5) p4 add (3,6,0) p5 add (2,5,5) &p6 add (3,0,5) p7 add (3,6,5)gen zone brick‎size 1 5 5 &p0 (5,0,5) p1 add (1,0,0) p2 add (0,5,0) &p3 add (0,0,5) p4 add (1,6,0) p5 add (0,5,5) &p6 add (1,0,5) p7 add (1,6,5)gen zone brick‎size 5 1 5 &p0 (0,3,0) p1 add (3,0,0) p2 add (0,3,0) &p3 add (0,2,5) p4 add (6,3,0) p5 add (0,3,5) &p6 add (5,2,5) p7 add (6,3,5)gen zone brick‎size 5 1 5 &p0 (0,5,5) p1 add (5,0,0) p2 add (0,1,0) &p3 add (0,0,5) p4 add (6,1,0) p5 add (0,1,5) &p6 add (5,0,5) p7 add (6,1,5)group‎Bin range‎group‎Mater‎i al not; Creat‎e named‎range‎synon‎y msrange‎name=Bin group‎Binrange‎name=Mater‎i al group‎Mater‎i al; Assig‎n model‎s to group‎smodel‎mohr range‎Mater‎i almodel‎elas range‎Bingen separ‎a te Mater‎i alinter‎f ace 1 wrap Mater‎i al Bin range‎plane‎ori 0 0 0 norma‎l 1 -1 0 above‎inter‎f ace 2 wrap Mater‎i al Bin range‎plane‎ori 0 0 0 norma‎l 1 -1 0 below‎int 1 maxed‎g e 0.55int 2 maxed‎g e 0.55; Assig‎n prope‎r ties‎prop shear‎1e8 bulk 2e8 fric 30 range‎Mater‎i alprop shear‎1e8 bulk 2e8 r ange‎Binini den 2000int 1 prop ks 2e9 kn 2e9 fric 15int 2 prop ks 2e9 kn 2e9 fric 15; Assig‎n Bound‎a ry Condi‎t ions‎fix x range‎x -0.1 0.1 any x 5.9 6.1 any fix y range‎y -0.1 0.1 any y 5.9 6.1 any fix z range‎z -0.1 0.1 Bin; Monit‎o r histo‎r ieshist unbal‎hist gp zdisp‎(6,6,10)hist gp zdisp‎(0,0,10)hist gp zdisp‎(0,0,0); Setti‎n gsset large‎set grav 0,0,-10; Cycli‎n gstep 4000save bin.sav9-5;;------------------------------------------------------------------;rest pile_‎g eom.sav‎group‎claymodel‎mohr rangeprop bulk 8.333e7‎ shear‎3.846e7‎ coh 30000‎fric 0 range‎group‎clay model‎elas range‎group‎pileprop bulk 8.333e7‎ shear‎3.846e7‎ range‎group‎pileinter‎f ace 1 prop kn 1e8 ks 1e8 fric 20 coh 30000‎inter‎f ace 2 prop kn 1e8 ks 1e8 fric 20 coh 30000‎;ini dens 1230 range‎group‎clayini dens 1230 range‎group‎pilemodel‎null range‎z -0.1 0.15;fix z range‎z -8.1 -7.9fix x range‎x -8.1 -7.9fix x range‎x 7.9 8.1fix y range‎y -.1 .1fix y range‎y 7.9 8.1set grav 0 0 -10ini szz 0. grad 0 0 12300‎. range‎z -5.5 0.ini szz 17600‎grad 0 0 15500‎range‎z -8 -5.5ini sxx 0. grad 0 0 5271.4 range‎z -5.5 0.ini sxx 7542.86 grad 0 0 6642.86 range‎z -8 -5.5ini sxx add 31428‎.6 grad 0 0 5714.3 range‎z -8 -5.5ini syy 0. grad 0 0 5271.4 range‎z -5.5 0.ini syy 7542.86 grad 0 0 6642.86 range‎z -8 -5.5ini syy add 31428‎.6 grad 0 0 5714.3 range‎z -8 -5.5;water‎densi‎t y 1000water‎table‎origi‎n 0,0,-5.5 norma‎l 0 0 -1ini dens 1550 range‎z -8 -5.5hist unbal‎;solve‎rat 1.e-6save pile0‎.sav;model‎elas range‎group‎pileprop bulk 13.9e9 shear‎10.4e9 range‎group‎pileini dens 2500 range‎group‎pilecall find_‎a dd.fissolve‎rat 1.e-6save pile1‎.savrest pile1‎.sav ;调用保存文‎件ini state‎0ini xdis 0.0 ydis 0.0 zdis 0.0 ;位移清零apply‎szz -0.4e6 range‎z 0.05 0.15 group‎pile ;桩顶加第一‎级荷载solve‎save app0.4.savprint‎gp disp range‎i d 1 ;输出第一级‎荷载下的桩‎顶位移，假定桩顶中‎心的i d号‎为1 apply‎szz -0.6e6 range‎z 0.05 0.15 group‎pile ;桩顶加第二‎级荷载solve‎save app0.6.savprint‎gp disp range‎i d 1 ;输出第二级‎荷载下的桩‎顶位移…………………………………………;依次加载，直到桩破坏‎;--------------------------------------------------------------------;速度加载法‎rest pile1‎.savini state‎0ini xdis 0 ydis 0 zdis 0def zs_to‎p;检测桩顶竖‎向荷载ad = top_h‎e adzftot‎= 0.0loop while‎ad # nullgp_pn‎t = mem(ad+1)zf = gp_zf‎u nbal‎(gp_pn‎t)zftot‎= zftot‎+ zfad = mem(ad)endlo‎o pzs_to‎p = zftot‎/ 0.1414endfix z range‎z 0.05 .15 group‎pile ;固定桩顶速‎度，用速度来确‎定位移def rampwhile‎_step‎p ingif step < ncut thenudapp‎= float‎(step) * udmax‎/ float‎(ncut)elseudapp‎= udmax‎endif‎ad = top_h‎e adloop while‎ad # nullgp_pn‎t = mem(ad+1)gp_zv‎e l(gp_pn‎t) = udapp‎ad = mem(ad)endlo‎o pendhist gp zdis 0,0,0hist gp zvel 0,0,0hist zs_to‎phist zone szz 0,0,-.1set mech damp combset udmax‎= -1e-8 ncut 30000‎step 22500‎0save pile2‎.sav;;------------------------------------------------------------------------;位移控制法‎def solve‎_step‎sloop n (1,21)save_‎f ile = strin‎g(n) + '-step.sav'comma‎n ds tep 40000‎s ave save_‎f ilepri zone stres‎s ran id 2381 a id 2361 a id 2341 a ;输出桩顶网‎格单元的应‎力 endco‎m mand‎endlo‎o pendsolve‎_step‎s;-----------------------------------------------------------------------------sel beam beg 0 0 0 end 2 0 0 nseg 2sel beam beg 2 0 0 end 4 0 -1 nseg 3sel beam id=2 beg 4 0 -1 end 5 0 -2 nseg 2 plot sel geo id on nod on scale‎0.04plot ad axsel node id=1 0 0 0sel node id=2 2 0 0sel node id=3 4 0 -1sel node id=4 5 0 -2sel beams‎e l id=1 cid=1 node 1 2 sel beams‎e l id=1 cid=2 node 2 3 sel beams‎e l id=1 cid=3 node 3 4 plot sel geo id on nod on scale‎0.04 plot ad axdef set_v‎a lsptA = 25.0 * sin( 40.0*degra‎d )ptB = 25.0 * cos( 40.0*degra‎d )endset_v‎a lsgen zone cylin‎d er p0=( 0.0, 0.0, 0.0 ) &p1=( ptA, 0.0, ptB ) &p2=( 0.0, 25.0, 0.0 ) &p3=( 0.0, 0.0, 25.0 ) &p4=( ptA, 25.0, ptB ) &p5=( 0.0, 25.0, 25.0 ) &size=(1, 2, 2)sel shell‎id=5 range‎cylin‎d er end1=(0.0, 0.0,0.0) &e nd2=(0.0,25.0,0.0) radiu‎s=24.5 not plot blo groplot ad sel geom black‎black‎cid on scale‎=0.03plot ad axpaudelet‎e ; delet‎e all zones‎sel node init zpos add -25.0newgen zone brick‎size 6 8 8model‎mohrprop bulk 1e8 shear‎0.3e8 fric 35prop coh 1e10 tens 1e10set grav 0 0 -9.81ini dens 1000fix x range‎x -0.1 0.1fix x range‎x 5.9 6.1fix y range‎y -0.1 0.1fix y range‎y 7.9 8.1fix z range‎z -0.1 0.1hist n 5hist unbal‎set mech force‎50solve‎save beam-brace‎0.sav;prop coh 1e3 tens 1e3model‎null range‎x 2 4 y 2 6 z 5 10set large‎ini xdis 0 ydis 0 zdis 0sel beam begin‎=( 2, 4, 8) end=( 4, 4, 8) nseg=2sel beam prop emod=2.0e11 nu=0.30sel beam prop XCAre‎a=6e-3 XCIz=200e-6 XCIy=200e-6 XCJ=0.0 hist gp zdisp‎4 4 8solve‎save beam-brace‎1.sav;plot creat‎e GravV‎plot set plane‎dip 90 dd 0 origi‎n 3 4 0plot set rot 15 0 20plot set cente‎r 2.5 4.2 4.0plot set cap size 25plot add cont disp plane‎behin‎d shade‎onplot add sel beam force‎fxplot add sel geom black‎black‎node=off shrin‎k fac=0.0plot add axesplot show;非全长锚固‎、预紧力锚杆‎（锚索）模拟;方法1、通过删除－建立lin‎k链接来模‎拟托盘gen zone radtu‎n p0 0,0,0 p1 25,0,0 p2 0,50, 0 p3 0,0,25 size 4 25 4 10 dim 4 4 4 4 ratio‎1 1 1 1.1 fillgen zone refle‎c t norma‎l 1 0 0 ori 0 0 0gen zone refle‎c t norma‎l0 0 1 ori 0 0 0mo mohrpro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5ini dens 2000fix x range‎x -25.1,-24.9fix x range‎x 24.9 25.1fix y range‎y 49.9 50.1fix z range‎z -25.1 -24.9fix z range‎z 24.9 25.1sel cable‎id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10sel cable‎id=1 beg 0,29,0 end 0,35,0 nseg 6sel cable‎id=1 prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 1 gr_k 1 gr_pe‎r 0.0785 range‎cid 1,10sel cable‎id=1 prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 10e5 gr_k 2e7 range‎cid 11,17sel delet‎e link range‎id 1sel link id=100 1 targe‎t zonesel link attac‎h xdir=rigid‎ydir=rigid‎zdir=rigid‎xrdir‎=rigid‎yrdir‎=rigid‎zrdir‎=rigid‎range‎id 100sel cable‎id=1 prete‎n sion‎60e3 range‎cid 1,10step 2000sav 10-5.sav;非全长锚固‎、预紧力锚杆‎（锚索）模拟;方法2、通过设置极‎大锚固剂参‎数模拟托盘‎gen zone radtu‎n p0 0,0,0 p1 25,0,0 p2 0,50, 0 p3 0,0,25 size 4 25 4 10 dim 4 4 4 4 ratio‎1 1 1 1.1 fillgen zone refle‎c t norma‎l 1 0 0 ori 0 0 0gen zone refle‎c t norma‎l0 0 1 ori 0 0 0mo mohrpro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5ini dens 2000fix x range‎x -25.1,-24.9fix x range‎x 24.9 25.1fix y range‎y 49.9 50.1fix z range‎z -25.1 -24.9fix z range‎z 24.9 25.1sel cable‎id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10sel cable‎id=1 beg 0,29,0 end 0,35,0 nseg 6sel cable‎ prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 1 gr_k 1 gr_pe‎r 0.0785 range‎cid 2,10sel cable‎ prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 10e5 gr_k 2e7 range‎cid 11,17sel cable‎ prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 10e8 gr_k 2e10 range‎cid 1,1sel cable‎id=1 prete‎n sion‎60e3 range‎cid 1,10step 2000sav 10-6.sav;非全长锚固‎、预紧力锚杆‎（锚索）模拟;方法3、借助别的结‎构单元（如line‎r单元）来模拟托盘‎gen zone radtu‎n p0 0,0,0 p1 25,0,0 p2 0,50, 0 p3 0,0,25 size 4 25 4 10 dim 4 4 4 4 ratio‎1 1 1 1.1 fillgen zone refle‎c t norma‎l 1 0 0 ori 0 0 0gen zone refle‎c t norma‎l0 0 1 ori 0 0 0mo mohrpro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5ini dens 2000fix x range‎x -25.1,-24.9fix x range‎x 24.9 25.1fix y range‎y 49.9 50.1fix z range‎z -25.1 -24.9fix z range‎z 24.9 25.1sel cable‎id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10sel cable‎id=1 beg 0,29,0 end 0,35,0 nseg 6sel cable‎id=1 prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 1 gr_k 1 gr_pe‎r 0.0785 range‎cid 1,10sel cable‎id=1 prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 10e5 gr_k 2e7 range‎cid 11,17sel liner‎range‎y=-.1, .1 x=-1,1 z=-1,1sel liner‎PROP iso=( 25e9, 0.15) thick‎=0.1 ; concr‎e tesel liner‎PROP cs_nk‎=8e8 cs_sk‎=8e8 &cs_nc‎u t=0.0 cs_sc‎o h=0.0 cs_sc‎o hres‎=0.0 cs_sf‎r ic=0.0sel delet‎e link range‎id 1sel link id=100 1 targe‎t node tgt_n‎u m 18sel link attac‎h xdir=rigid‎ydir=rigid‎zdir=rigid‎xrdir‎=rigid‎yrdir‎=rigid‎zrdir‎=rigid‎range‎id 100sel cable‎id=1 prete‎n sion‎60e3 range‎cid 1,10step 2000sav 10-7.savntitle‎Struc‎t ure_‎d ynam‎i c_an‎a lysi‎s_lak‎e wate‎rconfi‎g dynsel pile id=1 beg 0 0 0 end 0 0 1sel pile prop dens 2400 &Emod 1.0e10 Nu 0.3 XCAre‎a 0.3 &XCJ 0.16375‎XCIy 0.00625‎XCIz 0.01575‎&Per 2.8 &CS_sK‎1.3e11 CS_sC‎o h 0.0 CS_sF‎r ic 10.0 &CS_nK‎1.3e11 CS_nC‎o h 0.0 CS_nF‎r ic 0.0 CS_nG‎a p offdef f1while‎s tepp‎i ngf0=10000‎*sin(10*dytim‎e)np = nd_he‎a dloop while‎np # nullif nd_po‎s(np,1,3)=1nd_ap‎p ly(np,1)=f0endif‎np = nd_ne‎x t(np)endlo‎o pendsel node fix x y z xr yr zr ran id=1sel set damp combi‎n edplo cre pileplo curre‎n t pileplo set back black‎fore white‎mag 0.8plo add sel geo id on sca .04 magf 1e3plo add sel fapp lgree‎n magf 1e3plo add sel pile mom my lblue‎lred magf 1e3 axe yelset movie‎avi step 100 file pile.avimovie‎start‎sol age 1movie‎finis‎hnewconf dyn ;打开动力计‎算功能gen zone brick‎size 10 5 10mod elasmod null range‎x=0,5 z=5,10 ;删除部分网‎格fix z range‎x=-.1 .1 z=.1 10.1 ;设置静力边‎界条件fix z range‎x=9.9,10.1 z=.1 10.1fix y range‎y=-.1 .1fix y range‎y=4.9 5.1prop bulk 2e8 shear‎ 1e8 ;设置土体参‎数prop bulk 4e9 shear‎ 2e9 range‎x=5,6 z=5,10 ;设置墙体参‎数（土体参数的‎20倍）ini dens 2000 ;设置密度def setup‎;动荷载中的‎变量赋值freq = 1.0omega‎= 2.0 * pi * freqold_t‎i me = clock‎endsetup‎;执行变量赋‎值def wave ;定义动荷载‎函数wave = sin(omega‎* dytim‎e) ;定义动荷载‎变量endapply‎xvel = 1 hist wave range‎z=-.1 .1 ;施加动荷载‎apply‎zvel = 0 r ange‎z=-.1 .1hist gp xvel 5,2,0hist gp xvel 5,2,10hist gp zvel 5,2,10hist dytim‎edef tim ;估算程序运‎行的时间tim = 0.01 * (clock‎- old_t‎i me)endset dyn multi‎ on ;设置动态多‎步solve‎age 1.0print‎ tim ;输出计算时‎间print‎dyn ;输出动力计‎算相关信息‎save mult1‎.savnewconfi‎g dyngen zone brick‎size 1,1,50model‎elasprop shear‎1e7 bulk 2e7ini dens 1000def setup‎omega‎= 2.0 * pi * freqpulse‎= 1.0 / freqendset freq=4.0setup‎def waveif dytim‎e > pulse‎wave = 0.0elsewave = 0.5 * (1.0 - cos(omega‎* dytim‎e))endif‎endrange‎name botto‎m z=-.1 .1fix z range‎z=.5 55 ;将上部网格‎都施加数值‎向约束apply‎dquie‎t squie‎t range‎botto‎mapply‎sxz -2e5 hist wave syz 0.0 szz 0.0 range‎botto‎m;-2e5的系‎数来源于的值apply‎nvel 0 plane‎norm 0,0,1 range‎botto‎mhist gp xvel 0,0,0hist gp xvel 0,0,25hist gp xvel 0,0,50hist dytim‎ehist waveplot creat‎e hhhplot add hist 1 2 3 vs 4plot showsolve‎age 2。

FLAC-命令流及其解释大全1. FLAC3D是有限元程序吗？答：不是！是有限差分法。

2. 最先需要掌握的命令有哪些？答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。

3. 怎样看模型的样子？答：plo blo gro可以看到不同的group的颜色分布4. 怎样看模型的边界情况？答：plo gpfix red sk5. 怎样看模型的体力分布？答：plo fap red sk6. 怎样看模型的云图？答：位移：plo con dis (xdis, ydis, zdis)应力：plo con sz (sy, sx, sxy, syz, sxz)7. 怎样看模型的矢量图？答：plo dis (xdis, ydis, zdis)8. 怎样看模型有多少单元、节点？答：plo info9. 怎样输出模型的后处理图？答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件10. 怎样调用一个文件？答：File/call或者call命令10. 如何施加面力？答：app nstress11. 如何调整视图的大小、角度？答：综合使用x, y, z, m, Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键12. 如何进行边界约束？答：fix x ran （约束的是速度，在初始情况下约束等效于位移约束）13. 如何知道每个单元的ID？答：用鼠标双击单元的表面，可以知道单元的ID和坐标14. 如何进行切片？答：plo set plane ori (点坐标) norm (法向矢量)plo con sz plane (显示z方向应力的切片)15. 如何保存计算结果？答：save +文件名16. 如何调用已保存的结果？答：rest +文件名；或者File / Restor17. 如何暂停计算？答：Esc18. 如何在程序中进行暂停，并可恢复计算？答：在命令中加入pause命令，用continue进行继续19. 如何跳过某个计算步？答：在计算中按空格键跳过本次计算，自动进入下一步20. Fish是什么东西？答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能21. Fish是否一定要学？答：可以不用，需要的时候查Mannual获得需要的变量就可以了22. FLAC3D允许的命令文件格式有哪些？答：无所谓，只要是文本文件，什么后缀都可以23. 如何调用一些可选模块？答：config dyn (fluid, creep, cppudm)后注：这个工作很繁琐，需要的时间很多，希望广大网友能够将自己曾经遇到的常见问题在后续跟贴，也为了将这个FAQ进行很好的充实。

《FLAC3D基础与工程实例》全部命令流《FLAC/FLAC 3D基础与工程实例》全部命令流1gen zone bri p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 &p4 15 15 0 p5 0 15 15 p6 15 0 10 p7 20 20 20 &size 10 10 10 rat 1.0 0.9 1.1 group brick_1gen zone bri p0 20 0 0 p1 add 10 0 0 p2 add 0 20 0 p3 add 0 0 15 &size 10 10 10 rat 1.0 0.9 1.1 group brick_2gen zone bri p0 40 0 0 edge 10 size 10 10 10 rat 1.0 0.9 1.1 group brick_3 plot surnewgen zon bri size 3 3 3model elasprop bulk 3e8 shear 1e8ini dens 2000fix z ran z -.1 .1fix x ran x -.1 .1fix x ran x 2.9 3.1fix y ran y -.1 .1fix y ran y 2.9 3.1set grav 0 0 -10solveapp nstress -10e4 ran z 3 x 1 2 y 1 2 hist gp vel 0 0 3hist gp vel 0 3 3plo hist 1 redplo add hist 2 bluesolve;-------------------------------------工程信息;Project Record Tree export;Title:Simple test;---------------------------------计算第一步;... STATE: STATE1 ....configgrid 10,10model elasticgroup 'User:Soil' notnullmodel elastic notnull group 'User:Soil'prop density=1500.0 bulk=3E6 shear=1E6 notnull group 'User:Soil' fix x y j 1fix x i 1fix x i 11set gravity=9.81history 999 unbalancedsolvesave state1.sav;----------------------------------计算第二步;... STATE: STATE2 ....initial xdisp 0 ydisp 0initial xvel 0 yvel 0model null i 4 7 j 8 10group 'null' i 4 7 j 8 10group delete 'null'history 1 xdisp i=4, j=11solvesave state2.sav;--------------------------------绘图命令;*** plot commands ****;plot name: syyplot hold grid syy fill;plot name: Unbalanced forceplot hold history 999;plot name: gridplot hold grid magnify 20.0 lred grid displacement;plot name: Xdis-Aplot hold history 1 linenew; ===============================; 定义球体半径和半径方向上单元网格数; ===============================def parmrad=10.0rad_size=5endparm; ===============================; 建立八分之一球体外接立方体网格; ===============================gen zone pyramid p0 rad 0 0 p1 rad 0 rad p2 rad rad 0 p3 0 0 0 & p4 rad rad rad size rad_size rad_size rad_size group 1gen zone pyramid p0 0 rad 0 p1 rad rad 0 p2 0 rad rad p3 0 0 0 & p4 rad rad rad size rad_size rad_size rad_size group 2gen zone pyramid p0 0 0 rad p1 0 rad rad p2 rad 0 rad p3 0 0 0 & p4 rad rad rad size rad_size rad_size rad_size group 3; ==================================; 利用FISH语言将内部立方体节点调整到球面; ==================================def make_spherep_gp=gp_headloop while p_gp#null; 获取节点点坐标值：P=(px,py,pz)px=gp_xpos(p_gp)py=gp_ypos(p_gp)pz=gp_zpos(p_gp)dist=sqrt(px*px+py*py+pz*pz)if dist>0 then; 节点位置调整maxp=max(px,max(py,pz))k=(maxp/rad)*(rad/dist)gp_xpos(p_gp)=k*pxgp_ypos(p_gp)=k*pygp_zpos(p_gp)=k*pzend_ifp_gp=gp_next(p_gp)end_loopendmake_sphere; ===============================; 利用镜像生成完整球体网格; =============================== gen zone refgen zone ref dip 90gen zone ref dip 90 dd 90; =============================== ; 显示球体网格; =============================== plot surf pl set back whpl bl grngen zon bri size 1 1 2 group soil ran z 1 1 group rock ran z 01 expgrid 1.flac3dgen zon bri size 3 3 3model mohrprop bu 3e6 sh 1e6 coh 10e3 fric 15 fix z ran z -.1 .1fix x ran x -.1 .1fix x ran x 2.9 3.1fix y ran y -.1 .1fix y ran y 2.9 3.1ini dens 2000hist unbalset grav 10solve elasticsave 6-1.savrest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2 solvesave 6-2.savrest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2 hist id=2 gp zdis 1 1 3hist id=3 gp zdis 1 1 2hist id=4 gp xdis 1 1 3hist id=5 gp xdis 1 1 3hist id=6 zone szz 1 1 3hist id=7 zone szz 1.5 1.5 2.5hist id=8 zone sxz 1.5 1.5 2.5solvesave 6-3.savrest 6-3.savset log onset logfile 6-2.log print zone stress print gp disset log offrest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2 plot set rot 20 0 30plot con szz ou on magf 10plot add hist 1set movie avi step 1 file 6-5.avimovie startsolvemovie finishngen zon bri size 1 1 2 model elasprop bulk 3e7 shear 1e7 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solveplo con szngen zon bri size 1 1 2model mohrprop bulk 3e7 shear 1e7 c 1e10 f 15 tension 1e10 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solveprop bulk 3e7 shear 1e7 c 10e3 f 15 ten 0solveplo con szngen zone brick size 1 1 2model mohrprop bulk 3e7 shear 1e7 coh 10e3 fri 15 ten 0 fix z ran z 0 fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solve elasplo con sznewgen zone brick size 1 1 2model mohrprop bulk 3e7 shear 1e7 coh 10e3 fri 15 ten 0 fix z ran z 0 fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000ini szz -40e3 grad 0 0 20e3 ran z 0 2ini syy -20e3 grad 0 0 10e3 ran z 0 2ini sxx -20e3 grad 0 0 10e3 ran z 0 2set grav 0 0 -10solveplo con szngen zon bri size 1 1 2model mprop bulk 3e7 shear 1e7 c 10e10 f 15 ten 1e10 fix z ran z 0 fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000 ran z 0 1ini dens 1500 ran z 1 2ini szz -35e3 grad 0 0 20e3 ran z 0 1ini syy -17.5e3 grad 0 0 10e3 ran z 0 1ini sxx -17.5e3 grad 0 0 10e3 ran z 0 1ini szz -15e3 grad 0 0 15e3 ran z 1 2ini syy -7.5e3 grad 0 0 7.5e3 ran z 1 2ini sxx -7.5e3 grad 0 0 7.5e3 ran z 1 2ini pp 10e3 grad 0 0 -10e3 ran z 0 1set grav 0 0 -10solveplo con szngen zon bri size 1 1 2model mprop bulk 3e7 shear 1e7 c 10e10 f 15 ten 1e10 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000 ran z 0 2ini szz -50e3 grad 0 0 20e3 ran z 0 1ini syy -30e3 grad 0 0 10e3 ran z 0 1ini sxx -30e3 grad 0 0 10e3 ran z 0 1ini pp 30e3 grad 0 0 -10e3 ran z 0 2app nstress -10e3 ran z 2set grav 0 0 -10solveplo con sznewgen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 & p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & size 6 6 10model elaspro bulk 10e10 she 10e10ini den 2500apply sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1 apply sxx -1e9 grad 0 0 6.6666666e6 range x 59.9 60.1 apply syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1 apply syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1 apply szz -1e8 grad 0 0 8.3333333e5 ran z 0 120set grav 0 0 -10step 30000ini xdisp 0 ydisp 0 zdisp 0ini xvel 0 yvel 0 zvel 0plo cont szznewgen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 & p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & size 6 6 10model elaspro bulk 10e10 she 10e10ini den 2500ini sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1ini sxx -1e9 grad 0 0 6.6666666e6 range x 59.9 60.1ini syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1ini syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1ini szz -1e8 ran z -.1 .1fix x y z ran z -.1 .1set grav 0 0 -10solveini xdisp 0 ydisp 0 zdisp 0ini xvel 0 yvel 0 zvel 0plo cont szz。

FLAC 实例分析教程刘波韩彦辉（美国）编著《FLAC原理实例与应用指南》北京：人民交通出版社，2005.9 Appendix(附录) 版权所有：人民交通出版社，Itasca Consulting Group, USA说明：1．本实例分析教程是为方便读者学习、应用FLAC和FLAC3D而编写的，作为《FLAC原理、实例与应用指南》一书的附录。

2．计算算例参考了Itasca Consulting Group的培训算例，命令流的解析旨在方便读者理解FLAC和FLAC3D建模及求解问题的一般原则与步骤。

3．实例分析的算例中，FLAC算例是基于FLAC 5.00版本、FLAC 3D算例是基于FLAC 3D 3.00版本实现计算分析的。

读者在学习和研究相关算例时，请务必采用Itasca授权的合法版本进行分析计算。

4．本实例分析教程仅供读者参考，读者在参考本教程算例进行工程分析时，编者对可能产生的任何问题概不负责。

编者2005.10.18Project: [tunnel.prj] 隧道分析 -- Example 1-1Record Tree•[new]•[tun1.sav]o configo grid 10,10o；10*10建立网格o model elastico；设定为弹性模型o gen circle 5.0,5.0 2.0o；生成圆, 该圆圆心位置为（5.0 5.0），半径为2.0o group 'Tunnel:strong rock' notnullo；设定非零的区域为group 'Tunnel:strong rock'o model mohr notnull group 'Tunnel:strong rock'o；设定group'Tunnel:strong rock'的非零模型的区域为弹性模型o prop density=2000.0 bulk=1E8 shear=3E7 cohesion=1000000.0 friction=35.0 dilation=0.0 tension=0.0 notnull group 'Tunnel:strong rock'o；对group 'Tunnel:strong rock'的非空区域，设定模型材料参数。

《隧道及地下工程FLAC解析方法》全部命令流汇总隧道及地下工程FLAC解析方法是一种用于模拟隧道和地下工程行为的计算程序。

它可以对不同条件下的隧道和地下工程进行分析和优化设计。

以下是隧道及地下工程FLAC解析方法的全部命令流汇总：1.开始命令：FLAC-这个命令启动了FLAC程序，告诉计算机将要运行隧道及地下工程FLAC解析方法。

2.模型概况命令：MODELDIMENSION-这个命令设置了模型的维度，可以是2D或3D，具体取决于模型的需求。

3.材料定义命令：MAT-这个命令定义了隧道或地下工程中使用的材料的参数，如弹性模量、泊松比、摩擦角等。

4.边界条件命令：BOUNDARY-这个命令设置了模型的边界条件，包括加载、固定或自由应力等。

5.网格设置命令：GRID-这个命令定义了模型的网格，并且可以对网格进行细分或剖分，以适应复杂的地质条件。

6.初始条件命令：INITIAL-这个命令设置了模型的初始条件，如应力、位移或速度。

7.应力平衡命令：STRESS/EQUILIBRIUM-这个命令用于检查模型中的应力平衡情况，并调整模型的参数以满足平衡要求。

8.载荷应用命令：LOAD/APPLY-这个命令指定了要应用到模型中的加载条件，可以是恒定载荷、动态载荷或变化载荷。

9.运行命令：RUN-这个命令启动了模型的运行，并进行了迭代求解以得到模型的响应。

10.取样命令：SAMPLE-这个命令用于对模型中的节点或单元进行采样，以获取特定时间或位置的应力、位移或应变等信息。

11.结果输出命令：PLOT-这个命令用于设置结果的输出方式，可以选择输出为图形、数据文件或报表的形式。

12.模型评估命令：EVALUATE-这个命令用于对模型的结果进行评估，可以比较不同模型或不同条件下的结果。

13.结束命令：END-这个命令结束了FLAC程序的运行。

以上是隧道及地下工程FLAC解析方法的全部命令流汇总，通过这些命令可以对隧道及地下工程进行模拟和分析，得到相关的结果和优化设计。

【常规显示】PLOT con syy ;表示Y方向应力云图PLOT con sxy ;表示XY方向的应力云图PLOT con szz ou on ;ou是outline的缩写，用来显示网格PLOT con szz ef on ;ef是effective的缩写，用来显示z方向有效应力的云图PLOT con szz inter 1e4 ; inter是interval的缩写，以改变云图显示中的增量大小，本例中默认为5.0E3PLOT con szz max -10e3 ;max表示显示的应力数值最大不超过-10e3，本例中默认为0 PLOT con szz min -40e3 ;min表示显示的应力数值最小值为-40e3，本例中默认为-60e3 PLOT bcon prop bu ;显示体积模量bulk的赋值块云图【面荷载作用下的变形】rest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0 初始应变归零app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2 施加面荷载solvesave 6-2.savPLOT con zd 分析面荷载作用下的沉降云图PLOT con zd ou on magf 20 变形后的网格图（20为放大倍数）PLOT con szz ou on magf 20 将变形后的网格与应力云图结合起来plot block state 显示塑性区域plot block state shear ;获得剪切屈服的单元，包括shear-n和shear-pplot block state tension-p ;获得过去拉伸屈服的单元plot block state now ;获得当前处于塑性状态的单元，包括shear-n和tension-nplot block state past ;获得过去处于塑性状态的单元，包括shear-p和tension-p【历史跟踪】rest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2hist id=2 gp zdis 1 1 3 记录此节点的Z方向位移，记录号为2hist id=3 gp zdis 1 1 2hist id=4 gp xdis 1 1 3hist id=5 gp xdis 1 1 2hist id=6 zone szz 1 1 3hist id=7 zone szz 1.5 1.5 2.5hist id=8 zone sxz 1.5 1.5 2.5solvesave 6-3.savPLOT hist 2 查看记录的2号记录PLOT hist 7 v 8 ;v是vs的缩写hist write 7 v 8 file 6-3hist.txt 使用该命令，程序会在默认文件夹中生成一个名为6-3hist.txt 的文本文件，读者可以打开这个文本文件观察所得的数据结果。

;****************************************group命令group soil range z 1 2 ;定义groupprop bulk 7.8e6 shear 3.0e6 coh 10e3 fric 15 ran group soil ;给group赋值model null range group soilmodel elastic range group soilplot block group range group dam ;只显示damrange命令range name trench x 0 1 y 0 4 z 0 2model null range trenchmodel null range x=2,4 y=2,6 z=5,10range name Big_Brick x -3 3 y -2 2 z -1 1model elastic range Big_Brickprop bulk 1e8 shear 1e8 range Big_Brickrange name Layer1 plane dip 0 dd 0 ori 0 0 0 aboverange name Layer2 plane dip 0 dd 0 ori 0 0 0 belowrange cylinder end1 x1 y1 z1 end2 x2 y2 z2 radius rcylindrical range with one end of the cylinder axis (end1) at location (x1, y1, z1), the other end (end2) at location (x2, y2, z2), and with a cylinder radius of r '由(x1, y1, z1) 、(x2, y2, z2)两点确定旋转轴macro命令macro Sand 'bulk 1e8 shear 0.5e8 coh 0 tens 0 fric 35'macro Clay 'bulk 1e7 shear 0.3e7 coh 1e7 tens 0 fric 0'prop sand range Layer1prop clay range layer2macro Pt0 'p0 0 0 0'macro Pt1 'p1 add 10 0 0'macro Pt2 'p2 add 0 10 0'macro Pt3 'p3 add 0 0 10'macro Model_Size 'size 4 5 6'macro Big_Brick 'gen zone brick Pt0 Pt1 Pt2 Pt3 Model_Size'Big_Brickmacro 'Pt0' 'p0 15 15 15'gen Big_Brick ; this will cause an errorgroup range macro 命令Using different object types to do the same jobusing a RANGE object ...range name Big_Brick x = (-3,3) y = (-2,2) z = (-1,1)model null range Big_Brickusing a GROUP object ...group Big_Brick range x = (-3,3) y = (-2,2) z = (-1,1)model null range group Big_Brickusing a macro object ...macro Big_Brick 'x = (-3,3) y = (-2,2) z = (-1,1)'model null range Big_Brickgen zone brick size 10,10,10macro SiltySand 'bulk 1.5e8 shear 0.3e8'macro ClayeyGravel 'bulk 1.5e8 shear 0.6e8 fric 30 coh 5e6 ten 8.66e6'model elas range z 5 10prop SiltySand range z 5 10model mohr range z 0 5prop ClayeyGravel range z 0 5;****************************************视图操作：x y z mctr+r 还原ctr+p 保存图片ctr+z 鼠标选择ctr+g 彩图变为灰色图shift + （x y z m）旋转缩放edit---copy to clipboard----粘贴到word;****************************************ini xdis 0 ydis 0 zdis 0set grav 0 0 -9.81set mech force=50 '最大不平衡力小于50n，停止计算;****************************************hist命令hist reset '清除已有历史信息set hist_rep 1 '每一个时步记录一次hist n=5 '每5个时步记录一次hist write 1 3 vs 2 begin 150 end 375 file xx.txthist gp zdisp 4,4,8hist id=3 gp zdis 1 1 3 'id为监测变量的编号，默认的是从1开始编号。

gen zone brick size 6 8 8 p0 -10 -10 -20 &p1 10 -10 -20 &p2 -10 10 -20 &p3 -10 -10 0plot surfplot his -2 vs 1 为输出一二的关系图，如果应力值都为负数，则加-2是为了应力坐标轴为正plot add cont szz out on shade on 打开szz的等值线图plo con zdisp outline on shade on 打开zdisp的等值线云图plot set magf 1.0视图的放大倍数为1.0plo con szz z方向应力云图gen merge 1e-5 range z 0此命令是接触面单元合并成一个整体，1e-5是容差渐变应力施加：apply nstress -1e6 gradient 0,0,1e5 range z 3.464,0 plane dip 60dd 270 origin .1 0 0；施加法向应力：apply nstress -1e6 range plane dip 60 dd 270 origin .1 0 print gp position range id=14647 输出节点坐标print zone state都可以用来显示模型的塑性区Hist gp zdisp 1；Hist gp szz 2输出2和1的hist关系图用：plot his -2 vs 1 如果应力值都为负数，则加-2是为了应力坐标轴为正！或者hist id=？ Gp zdis x y z 是记录xyz坐标的位移值并顶为ID=？print zone stress 在命令窗口中会输出模型全部单元的6个方向应力数值print gp dis 该命令会显示模型中全部节点三个方向的变形大小hist unbal监测不平衡力，并保留历史记录hist gp zdisp 4,4,8监测网格坐标点（4，4，8）在z方向的位移，并保留历史纪录HIST gpvel<x y z>记录某个节点的速度；PLOT vel命令绘制完整的速度矢量场图；如果都在水下，则孔隙水压力为ini pp 30e3 grad 0 0 -10e3 ran z 0 3渐变的过程删除其它群组只保留群组2 del ran group 2 not桩单元：sel pile id=1 beg 0 0 0 end 0 0 10 nseg 4 ID为1 长度10m，五个结构节点（端面两个，杆件三个）；四段pile单元构件（结构杆件）CID结构构件的编号 Node结构节点的编号SEL pile prop rockbolt on 激活锚杆特性锚索单元：sel cable id=() begin=() end=() nseg=() 类似桩单元怎样查看模型的矢量图？答：pl grid vel (disp) red4.问：怎么能把一个PLOT的图像数据导出来以便用其他软件绘图？回答：用set log on 命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。

第二章FLAC3D原理及入门指南2-1定义一个FISH函数newdef abcabc = 25 * 3 + 5Endprint abc2-2使用一个变量newdef abchh = 25abc = hh * 3 + 5EndPrint hhPrint abc2-3对变量和函数的理解newdef abchh = 25abc = hh * 3 + 5Endset abc=0 hh=0print hhprint abcprint hhnewdef abcabc = hh * 3 + 5endset hh=25print abcset abc=0 hh=0print hhprint abcprint hh2-4获取变量的历史记录newgen zone brick size 1 2 1model mohrprop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10fix x y z range y -0.1 0.1apply yvel -1e-5 range y 1.9 2.1plot set rotation 0 0 45plot block groupdef get_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_addef loadload=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4) endhist loadhist gp ydis 0,2,0step 1000plot his 1 vs -22-5用FISH函数计算体积模量和剪砌模量newdef derives_mod = y_mod / (2.0 * (1.0 + p_ratio))b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio))endset y_mod = 5e8 p_ratio = 0.25deriveprint b_modprint s_mod2-6 在FLAC输入中使用符号变量Newdef derives_mod = y_mod / (2.0 * (1.0 + p_ratio))b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio))endset y_mod = 5e8 p_ratio = 0.25derivegen zone brick size 2,2,2model elasticprop bulk=b_mod shear=s_modprint zone prop bulkprint zone prop shear2-7 控制循环Newdef xxxsum = 0prod = 1loop n (1,10)sum = sum + nprod = prod * nend_loopendxxxprint sum, prodnewgen zone brick p0 (0,0,0) p1 (-10,0,0) p2 (0,10,0) p3 (0,0,-10) model elasplot set rotation 0 0 45plot block groupdef installpnt = zone_headloop while pnt #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_loopendset p_ratio=0.25 y_zero=1e7 cc=1e8install2-8 拆分命令行new ;example of a sum of many thingsdef long_sumtemp = v1 + v2 + v3 + v4 + v5 + v6 + v7 + v8 + v9 + v10long_sum = temp + v11 + v12 + v13 + v14 + v15end2-9 变量类型newdef haveoneaa = 2bb = 3.4cc = ’Have a nice day’dd = aa * bbee = cc + ’, old chap’endhaveoneprint fish2-10 IF条件语句newdef abcif xx > 0 thenabc = 1000elseabc = -1000end_ifendset xx = 10print abcset xx = 0print abc2-11 索单元自动生成newgen zone brick size 10 3 5plot set rotation 0 0 45plot block groupdef place_cablesloop n (1,5)z_d = float(n) - 0.5commandsel cable beg 0.0,1.5,z_d end 7.0,1.5,z_d nseg 7 end_commandend_loopendplace_cablesplot grid sel geom rednewgen zone brick size 10 3 5plot set rotation 15 0 60plot block groupmod mohrprop bulk 1e8 shear .3e8 fric 35prop coh 1e3 tens 1e3ini dens 1000set grav 0,0,-10fix x y z range z -.1 .1fix y range y -.1 .1fix y range y 2.9 3.1fix x range x -.1 .1fix x range x 9.9 10.1set largehist unbalsolvesave cab_str.savini xdis 0 ydis 0 zdis 0hist gp xdisp 0,1,5def place_cablesloop n (1,5)z_d = 5.5 - float(n)z_t = z_d + 0.5z_b = z_d - 0.5commandfree x range x -.1,.1 z z_b z_tsolvesel cable beg 0.0,0.5,z_d end 7.0,0.5,z_d nseg 7sel cable beg 0.0,1.5,z_d end 7.0,1.5,z_d nseg 7sel cable beg 0.0,2.5,z_d end 7.0,2.5,z_d nseg 7sel cable prop emod 2e10 ytension 1e8 xcarea 1.0 & gr_k 2e10 gr_coh 1e10 gr_per 1.0end_commandend_loopendplace_cablessave cab_end.savplot sketch sel cable force red2-12圆形隧道开挖模拟计算;建立模型gen zon radcyl p0 0 0 0 p1 6 0 0 p2 0 1 0 p3 0 0 6 &size 4 2 8 4 dim 3 3 3 3 rat 1 1 1 1.2 group outsiderockgen zone cshell p0 0 0 0 p1 3 0 0 p2 0 1 0 p3 0 0 3 &size 1 2 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group concretliner fill group insiderock gen zon reflect dip 90 dd 90 orig 0 0 0gen zon reflect dip 0 dd 0 ori 0 0 0gen zon brick p0 0 0 6 p1 6 0 6 p2 0 1 6 p3 0 0 13 size 4 2 6 group outsiderock1 gen zon brick p0 0 0 -12 p1 6 0 -12 p2 0 1 -12 p3 0 0 -6 size 4 2 5 group outsiderock2 gen zon brick p0 6 0 0 p1 21 0 0 p2 6 1 0 p3 6 0 6 size 10 2 4 group outsiderock3 gen zon reflect dip 0 dd 0 orig 0 0 0 range group outsiderock3gen zon brick p0 6 0 6 p1 21 0 6 p2 6 1 6 p3 6 0 13 size 10 2 6 group outsiderock4gen zon brick p0 6 0 -12 p1 21 0 -12 p2 6 1 -12 p3 6 0 -6 size 10 2 5 group outsiderock5 gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z 6.1 13.1gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z -6.1 -12.1gen zon reflect dip 90 dd 90 orig 0 0 0 range x 6.1 21.1 z -12.1 13.1;绘制模型图plot block groupplot add axes red;plot set rotation 0 0 45 用于显示三维模型;设置重力set gravity 0 0 -10;给定边界条件fix z range z -12.01,-11.99fix x range x -21.01,-20.99fix x range x 20.99,21.01fix y range y -0.01 0.01fix y range y 0.99,1.01;求解自重应力场model mohrini density 1800 ;围岩的密度prop bulk=1.47e8 shear=5.6e7 fric=20 coh=5.0e4 tension=1.0e4 ;体积、剪切、摩擦角、凝聚力、抗拉强度set mech ratio=1e-4solvesave Gravsol.savplot cont zdisp outl onplot cont szz;毛洞开挖计算initial xdisp=0 ydisp=0 zdisp=0model null range group insiderock any group concretliner anyplot block groupplot add axes redset mech ratio=5e-4solvesave Kaiwsol.savplot cont zdispplot cont sdispplot cont szzplot cont xzz;模筑衬砌计算model elas range group concretliner anyplot block groupplot add axes redini density 2500 range group concretliner any ;衬砌混凝土的密度prop bulk=16.67e9,shear=12.5e9 range group concretliner any ;衬砌混凝土的体积弹模、剪切弹模set mech ratio=1e-4solvesave zhihusol.savplot cont zdispplot cont sdispplot cont szzplot cont xzz;完成计算分析第四章FISH语言及建模技术4-1数组newdef afill ;fill matrix with random numbersarray var(4,3)loop m (1,4)loop n (1,3)var(m,n) = urandendloopendloopenddef ashow ;display contents of matrixloop m (1,4)hed = ’ ’msg = ’ ’+string(m)loop n (1,3)hed = hed + ’ ’+string(n)msg = msg + ’ ’+string(var(m,n))endloopif m = 1dum = out(hed)endifdum = out(msg)endloopendafillAshow4-2函数操作newdef xxxaa = 2 * 3xxx = aa + bbendnewdef stress_sumstress_sum = 0.0pnt = zone headloop while pnt # nullStress_sum = stress_sum + z_sxx(pnt)pnt=z_next(pnt)end_loopendnewdef stress_sumsum = 0.0pnt = zone headloop while pnt # nullsum = sum + z_sxx(pnt)pnt=z next(pnt)end loopstress_sum = sumendnewdef h_var_1ipz = z_near(1,2,3)H_var_1 = sxx(ipz) + facH_var_2 = syy(ipz) + facH_var_3 = szz(ipz) + facH_var_4 = sxy(ipz) + facH_var_5 = syz(ipz) + facH_var_6 = sxz(ipz) + facend4-3函数删除与重定义newdef joeii=out(‘This is A Function’)enddef fredjoeendfreddef joeii=(‘This is A New Function’)endfred4-4字符串newdef in_defxx = in(msg+’(’+’default:’+string(default)+’):’)if type(xx) = 3in_def = defaultelsein_def = xxendifend;def moduli_datadefault = 1.0e9msg=’Input Young‘s modulus ’Y_mod = in_defdefault = 0.25msg=’Input Poisson‘s ratio ’p_ratio = in_defif p_ratio = 0.5 thenii = out(’ Bulk mod is undefined at Poisson‘s ratio = 0.5’)ii = out(’ Select a different val ue --’)p_ratio = in_defendifs_mod = y_mod / (2.0 * (1.0 + p_ratio))b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio))end;moduli_datagen zone brick size 2,2,2model elasticprop bulk = b_mod shear = s_modprint p_ratio y_mod b_mod s_modpauseprint zone prop bulkpauseprint zone prop shear4-5马蹄形隧道网格newgen zone radcyl size 5 10 6 12 rat 1 1 1 1.2 &p0 0,0,0 p1 100,0,0 p2 0,200,0 p3 0,0,100gen zone radtun size 5 10 5 12 rat 1 1 1 1.2 &p0 0,0,0 p1 0,0,-100 p2 0,200,0 p3 100,0,0gen zone reflect dip 90 dd 270 origin 0,0,0plot set rotation 30 0 30plot block group4-6复杂形状网格生成new; 主隧道网格gen zon radcyl p0 15 0 0 p1 23 0 0 p2 15 50 0 p3 15 0 8 &size 4 10 6 4 dim 4 4 4 4 rat 1 1 1 1 fillgen zon reflect dip 90 dd 90 orig 15 0 0gen zon reflect dip 0 ori 0 0 0; 辅助隧道网格gen point id 1 (2.969848,0.0,-0.575736)gen point id 2 (2.969848,50.0,-0.575736)gen zon radcyl p0 0 0 -1 p1 7 0 0 p2 0 50 -1 p3 0 0 8 p4 7 50 0 &p5 0 50 8 p6 7 0 8 p7 7 50 8 p8 point 1 p10 point 2 & size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1gen zon radcyl p0 0 0 -1 p1 0 0 -8 p2 0 50 -1 p3 7 0 0 p4 0 50 -8 & p5 7 50 0 p6 7 0 -8 p7 7 50 -8 p9 point 1 p11 point 2 & size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1;衬砌网格sel shell range cyl end1 0 0 -1 end2 0 50 -1 rad 3; 隧道外围边界网格gen zone radtun p0 7 0 0 p1 50 0 0 p2 7 50 0 p3 15 0 50 p4 50 50 0 & p5 15 50 50 p6 50 0 50 p7 50 50 50 &p8 23 0 0 p9 7 0 8 p10 23 50 0 p11 7 50 8 &size 6 10 3 10 rat 1 1 1 1.1gen zone brick p0 0 0 8 p1 7 0 8 p2 0 50 8 p3 0 0 50 &p4 7 50 8 p5 0 50 50 p6 15 0 50 p7 15 50 50 &size 3 10 10 rat 1 1 1.1;对称得出1/2模型gen zon reflect dip 0 ori 0 0 0 range x 0 23 y 0 50 z 8 50gen zon reflect dip 0 ori 0 0 0 range x 23 50 y 0 50 z 0 50;建立主隧道和辅助隧道块名称group service range cyl end1 0 0 -1 end2 0 50 -1 rad 3group main range cyl end1 15 0 0 end2 15 50 0 rad 4;对称得出完整模型gen zon reflect dip 90 dd 270 ori 0 0 0;网格显示plot set rotation 30 0 30plot block group4-7网格连接newgen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,2gen zone brick size 8 8 4 p0 0,0,2 p1 4,0,2 p2 0,4,2 p3 0,0,4attach face range z 1.9 2.1model elasprop bulk 8e9 shear 5e9fix z range z -.1 .1fix x range x -.1 .1fix x range x 3.9 4.1fix y range y -.1 .1fix y range y 3.9 4.1apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2hist unbal;网格显示plot set rotation 30 0 30plot block groupsolvesave att.savplot cont zdisp outl on单一网格（小）程序newgen zone brick size 8 8 8 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,4 model elasprop bulk 8e9 shear 5e9fix z range z -.1 .1fix x range x -.1 .1fix x range x 3.9 4.1fix y range y -.1 .1fix y range y 3.9 4.1apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2hist unbal;网格显示plot set rotation 30 0 30plot block groupsolvesave noattx.savplot cont zdisp outl on单一网格（大）程序newgen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,4 model elasprop bulk 8e9 shear 5e9fix z range z -.1 .1fix x range x -.1 .1fix x range x 3.9 4.1fix y range y -.1 .1fix y range y 3.9 4.1apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2hist unbal;网格显示plot set rotation 30 0 30plot block groupsolvesave noattx.savplot cont zdisp outl on4-8立方体洞穴网格生成newdef parmrad=4.0len=10.0in_size=6rad_size=10endparmgen zone radbrick edge len size in_size in_size in_size rad_size &rat 1.0 1.0 1.0 1.2 dim rad rad rad;对称得出完整模型gen zon reflect dip 0 dd 0 ori 0 0 0gen zon reflect dip 90 dd 270 ori 0 0 0;网格显示plot set rotation 30 0 30plot block group4-9球体洞穴网格生成newdef parmrad=4.0len=10.0in_size=6rad_size=10endparmdef make_spherep_gp=gp_headloop while p_gp#nullpx=gp_xpos(p_gp)py=gp_ypos(p_gp)pz=gp_zpos(p_gp)dist=sqrt(px*px+py*py+pz*pz)if dist>0 thenk=rad/distax=px*kay=py*kaz=pz*kmaxp=max(px,max(py,pz))k=len/maxpbx=px*kby=py*kbz=pz*ku=(maxp-rad)/(len-rad)gp_xpos(p_gp)=ax+u*(bx-ax)gp_ypos(p_gp)=ay+u*(by-ay)gp_zpos(p_gp)=az+u*(bz-az)end_ifp_gp=gp_next(p_gp)end_loopendgen zone radbrick edge len size in_size in_size in_size rad_size & rat 1.0 1.0 1.0 1.2 dim rad rad radmake_sphere;对称得出完整模型gen zon reflect dip 0 dd 0 ori 0 0 0gen zon reflect dip 90 dd 270 ori 0 0 0;网格显示plot set rotation 15 0 30plot block group4-10应力边界newgen zone brick size (4,4,4) p0 (0,0,0) p1 (4,0,0) &p2 (0,4,0) p3 (2,0,3.464)model elasticprop bulk 1e8 shear .3e8apply nstress -1e6 range plane dip 60 dd 270 origin 0.1,0,0 above step 1;网格显示plot set rotation 30 0 30plot block groupplot grid red fap green4-11改变应力边界条件newgen zone brick size 6 6 6model elasprop bulk 1e8 shear 7e7fix x range x -0.1 0.1def supersteploop ns (1,n_steps)x_stress = stress_inccommandapply sxx add x_stress range x 5.9,6.1 y 0,6 z 0,2step 100end_commandend_loopendset n_steps=100 stress_inc=-1e3hist zone sxx 6,0,0plot create sxx_histplot showSuperstep;网格显示plot set rotation 30 0 30plot block groupplot grid red fap greennewgen zone brick size 6 6 6model elasprop bulk 1e8 shear 7e7fix x range x -0.1 0.1def x_stressx_stress = stress_inc * stependset stress_inc = -1e3apply sxx 1.0 hist x_stress range x 5.9,6.1 y 0,6 z 0,2 hist zone sxx 6,0,0hist x_stressstep 100;网格显示plot set rotation 30 0 30plot grid red fap greennewgen zone brick size 6 6 6model elasprop bulk 1e8 shear 7e7fix x range x -0.1 0.1table 1 0,0 100,-1e5apply sxx 1.0 hist table 1 range x 5.9,6.1 y 0,6 z 0,2 hist zone sxx 6,0,0step 100;网格显示plot set rotation 30 0 30plot grid red fap greennewgen zone brick size 6,6,6 p1 6,0,-1model elasprop bulk 8e9 shear 5e9apply sxx -2e6 range x -0.1 0.1apply sxx -2e6 range x 5.9 6.1step 500plot grid dispnewgen zone brick size 5,5,5model elasprop bulk 8e9 shear 5e9set grav 0 0 -10ini dens 1000fix x range x -.1 .1fix x range x 4.9 5.1fix y range y -.1 .1fix y range y 4.9 5.1ini szz -5e4 grad 0 0 -1e4app szz -5e4 range z -.1 .1solvemodel null range x 1,4 y 1 4 z 3 5step 100plot set plane dip 90 dd 180 origin 0,2.5,0plot add surf plane behind yellplot add vel plane behind blackplot set rotation 30 0 30plot block groupplot show4-12位移边界newgen zone brick size 4,4,4 p3 2,0,3.464model elasprop bulk 1e8 shear .3e8macro left_boun ’plane dip 60 dd 270 origin 0.1,0,0 above’macro right_boun ’plane dip 60 dd 270 origin 3.9,0,0 below’apply nvel 0.1 plane dip 60 dd 270 range left_bounapply nvel 0.1 plane dip 120 dd 90 range right_bounstep 1;网格显示plot set rotation 30 0 30plot block groupplot grid velnewnewgen zone brick size 4 4 4model elasticprop bulk 1e8 shear .3e8apply nstress -1e6 plane dip 0 dd 0 range y 3.9 4.1apply nvel 0.0 plane dip 0 dd 0 range y -.1 .1apply nvel 0.0 plane dip 0 dd 0 range x -.1 .1step 10;网格显示plot set rotation 30 0 30plot block groupplot grid velnewgen zone brick size 10 5 5mod elprop shear 1e8 bulk 2e8fix x y z range x -.1 .1 y 0 5 z 0 5fix x y z range x 0 10 y 0 5 z -.1 .1fix x y z range x 9.9 10.1 y 0 5 z 0 5table 1def find_addhead = nullp_gp = gp_headloop while p_gp # nullx_pos = gp_xpos(p_gp)if x_pos = width thennew = get_mem(2)mem(new) = headmem(new+1) = p_gphead = newendifp_gp = gp_next(p_gp)endloopendset width=10.0find_adddef apply_velwhile_steppingad = headloop while ad # nullp_gp = mem(ad+1)gp_xvel(p_gp) = vel_max * gp_zpos(p_gp) / heightgp_zvel(p_gp) = -vel_max * (gp_xpos(p_gp) - width) / height ad = mem(ad)endloopendset large vel_max=1e-2 height=5.0step 100;网格显示plot set rotation 30 0 30plot block groupplot grid vel4-13不考虑重力影响的均匀应力new ;set grav 0 0 0gen zone brick size 6 6 6model elasini sxx=-5e6 syy=-1e7 szz=-2e7apply sxx=-5e6 range x -0.1 0.1apply sxx=-5e6 range x 5.9 6.1apply syy=-1e7 range y -0.1 0.1apply syy=-1e7 range y 5.9 6.1apply szz=-2e7 range z -0.1 0.1apply szz=-2e7 range z 5.9 6.1prop bulk 8e9 shear 5e9ini dens 2000fix x range z -.1 .1step 10plot set rotation 30 0 30 ;网格三维显示plot cont zdisp outl on4-14考虑应力梯度的均匀材料newgen zone brick size 10 10 10 p1 20,0,0 p2 0,20,0 p3 0,0,20 model mohrprop bulk 5e9 shear 3e9 fric 35ini density 2500set gravity 0,0,-10fix x y z range z -0.1 0.1ini szz -5.0e6 grad 0,0,2.5e4ini syy -2.5e6 grad 0,0,1.25e4ini sxx -2.5e6 grad 0,0,1.25e4apply szz -4.5e6 range z 19.9 20.1apply szz -5.0e6 range z -0.1 0.1apply sxx -2.5e6 grad 0,0,1.25e4 range x -0.1 0.1apply sxx -2.5e6 grad 0,0,1.25e4 range x 19.9 20.1apply syy -2.5e6 grad 0,0,1.25e4 range y -0.1 0.1apply syy -2.5e6 grad 0,0,1.25e4 range y 19.9 20.1step 10;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-15考虑应力梯度的非均匀材料newgen zone brick size 10 10 10 &p0 0,-25,0 p1 20,-25,0 p2 0,0,0 p3 0,-25,20model elasprop bulk 5e9 shear 3e9ini density 1600 range y -10,0ini density 2000 range y -15,-10ini density 2200 range y -25,-15set gravity 0,-10,0fix x range x -.1 .1fix x range x 19.9 20.1fix z range z -.1 .1fix z range z 19.9 20.1fix y range y -25.1 -24.9ini syy 0.0 grad 0,1.6e4,0 range y -10,0ini syy 4e4 grad 0,2.0e4,0 range y -15,-10ini syy 7e4 grad 0,2.2e4,0 range y -15,-25step 10;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-16非均匀网格应力初始化newgen zone radcyl size 3 8 4 5 fill p1 10,0,0 p2 0,10,0 p3 0,0,10 mode elasticprop shear 3e8 bulk 5e8fix x range x -.1 .1fix x range x 9.9 10.1fix y range y -.1 .1fix y range y 9.9 10.1fix z range z -.1 .1ini szz = -2.5e5 grad 0,0,2.5e4ini density 2500set grav 0,0,-10step 10;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-17不规则自由面应力初始化gen zone brick size 15 15 10 p0 0,0,0 edge=100.0model elasticprop shear 3e8 bulk 5e8def mountaingp = gp_headloop while gp # nullzz = sqrt(gp_xpos(gp)ˆ2 + gp_ypos(gp)ˆ2)dz = 0.06 * sin(0.2 * zz + 100.0) ; Sum Fourier terms fordz = dz + 0.06 * sin(0.22 * zz - 20.3) ; quasi-random surfacedz = dz - 0.04 * sin(0.33 * zz + 33.3) ; topology.gp_zpos(gp) = 0.5 * gp_zpos(gp) * (1.0 + dz)gp = gp_next(gp)end_loopendmountainfix x range x -.1 .1fix x range x 99.9 100.1fix y range y -.1 .1fix y range y 99.9 100.1fix z range z -.1 .1set grav 0,0,-10ini density=2000ini szz=-2.0e6 (grad 0,0,2.0e4) sxx=-4.0e6 (grad 0,0,4.0e4) syy=-4.0e6 (grad 0,0,4.0e4) step 100plot set rotation 30 0 30plot cont zdisp outl on;网格三维显示4-18非均网格内部压实newgen zone brick size 8 8 10 ratio 1.2 1 1model mohrini dens 2000prop bulk 2e8 shear 1e8prop fric 30fix x range x -.1 .1fix x range x 7.9 8.1fix y range y -.1 .1fix y range y 7.9 8.1fix z range z -.1 .1;ini szz -2.0e5 grad 0,0,2e4;ini sxx -1.5e5 grad 0,0,1.5e4;ini syy -1.5e5 grad 0,0,1.5e4set grav 10step 1000;pause;prop tens 1e10 coh 1e10;step 750;prop tens 0 coh 0;step 250;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-19模型改变后初始应力变化newgen zone brick size 5 5 5model elasprop sh 2e8 bu 3e8fix x y z range z -.1 .1set grav 0 0 -10ini dens 2000Solve;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on;模型改变model mohr range x 0 2 y 0 5 z 0 2prop sh 2e8 bu 3e8 fric 35 range x 0 2 y 0 5 z 0 2 Step 2000;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-20应力与孔隙压力的初始化newconfig fluidgen zone brick size 8 5 10model elasmodel fl_isoini dens 2000prop bulk 1e9 shear 5e8prop poros 0.5 perm 1e-10ini fmod 2e9ini fdensity 1e3ini sat 0ini sat 1 range z -.1 5.1set grav 0 0 -10fix x range x -.1 .1fix x range x 7.9 8.1fix y range y -.1 .1fix y range y 4.9 5.1fix z range z -.1 .1ini pp 5.e4 grad 0,0,-1.e4 range z 0.0 5.ini szz -20e4 grad 0,0,20e3ini szz add -1.5e4 grad 0,0,.25e4 range z 5,6 ini szz add -2.5e4 grad 0,0,.5e4 range z 0,5 ini szz add -.25e4 range z 0,5solve;网格三维显示plot set rotation 30 0 30plot cont zdisp outl onnewgen zone brick size 1 1 10model elasini dens 2500 range z 0 5ini dens 2250 range z 5 6ini dens 2000 range z 6 10prop bulk 1e9 shear 5e8set grav 10water dens 1000water table ori 0 0 5 normal 0 0 1fix x range x -.1 .1fix x range x 7.9 8.1fix y range y -.1 .1fix y range y 4.9 5.1fix z range z -.1 .1ini szz -20e4 grad 0,0,20e3ini szz add -1.5e4 grad 0,0,.25e4 range z 5,6 ini szz add -2.5e4 grad 0,0,.5e4 range z 0,5 ini szz add -.25e4 range z 0,5solve;网格三维显示plot set rotation 30 0 30plot cont zdisp outl on4-21加载顺序new; 主隧道建模gen zon radcyl p0 15 0 0 p1 23 0 0 p2 15 50 0 p3 15 0 8 &size 4 10 6 4 dim 4 4 4 4 rat 1 1 1 1 fillgen zon reflect dip 90 dd 90 orig 15 0 0gen zon reflect dip 0 ori 0 0 0;辅助隧道建模gen point id 1 (2.969848,0.0,-0.575736)gen point id 2 (2.969848,50.0,-0.575736)gen zon radcyl p0 0 0 -1 p1 7 0 0 p2 0 50 -1 p3 0 0 8 p4 7 50 0 &p5 0 50 8 p6 7 0 8 p7 7 50 8 p8 point 1 p10 point 2 &size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1 fillgen zon radcyl p0 0 0 -1 p1 0 0 -8 p2 0 50 -1 p3 7 0 0 p4 0 50 -8 &p5 7 50 0 p6 7 0 -8 p7 7 50 -8 p9 point 1 p11 point 2 &size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1 fill;隧道边界网格gen zone radtun p0 7 0 0 p1 50 0 0 p2 7 50 0 p3 15 0 50 p4 50 50 0 & p5 15 50 50 p6 50 0 50 p7 50 50 50 &p8 23 0 0 p9 7 0 8 p10 23 50 0 p11 7 50 8 &size 6 10 3 10 rat 1 1 1 1.1gen zone brick p0 0 0 8 p1 7 0 8 p2 0 50 8 p3 0 0 50 &p4 7 50 8 p5 0 50 50 p6 15 0 50 p7 15 50 50 &size 3 10 10 rat 1 1 1.1gen zon reflect dip 0 ori 0 0 0 range x 0 23 y 0 50 z 8 50gen zon reflect dip 0 ori 0 0 0 range x 23 50 y 0 50 z 0 50group main1 range cyl end1 15 0 0 end2 15 25 0 rad 4group main2 range cyl end1 15 25 0 end2 15 50 0 rad 4group service range cyl end1 0.0,0.0,-0.575736 &end2 0.0,50.0,-0.575736 rad 3.0save tun0.savrest tun0.sav; 初始应力场模量model mohrprop shear 0.36e9 bulk 0.6e9 coh 1e5 fric 20 tens 1e5apply szz -1.4e6 range z 49.9 50.1fix z range z -50.1 -49.1fix x range x -.1 .1fix x range x 49.9 50.1fix y range y -.1 .1fix y range y 49.9 50.1ini sxx -1.4e6 syy -1.4e6 szz -1.4e6hist unbalhist gp xdis 3,0,-1hist gp zdis 0,0,2hist gp xdis 3,25,-1hist gp zdis 0,25,2step 1000save tun1.sav; 辅助隧道开挖25 m计算rest tun1.savini xdis 0.0 ydis 0.0 zdis 0.0model null range group service y 0,25step 1000save tun2.sav; 在辅助隧道上施加衬砌restore tun2.savsel shell id=1 range cyl end1 0 0 -1 end2 0 25 -1 rad 3sel shell prop iso=(25.3e9, 0.266) thick = 0.5sel node fix y xr zr range y -0.1 0.1 ; symmetry cond.sel node fix x yr zr range x -0.1 0.1 ; symmetry cond.model mohr range group service y 25 50prop shear 0.36e9 bulk 0.6e9 coh 1e5 fric 20 tens 1e5; 主隧道开挖25mmodel null range group main1step 1000save tun3.savplot set rotation 0 0 60plot cont zdisp outl on第六章双线铁路隧道施工过程分析6-2 V级围岩施工过程模拟（1）建立隧道1/4圆周模型gen zon radcyl p0 0 0 0 p1 10.55 0 0 p2 0 1 0 p3 0 0 10.55 &size 5 2 10 4 dim 5.55 5.55 5.55 5.55 rat 1 1 1 1.2 group outsiderocksyplot block groupplot add axes redgen zone cshell p0 0 0 0 p1 5.55 0 0 p2 0 1 0 p3 0 0 5.55 &size 1 2 10 4 dim 5.05 5.05 5.05 5.05 rat 1 1 1 1 group concretlinersy fill group insiderocksy gen zon radcyl p0 0 0 0 p1 0 0 -9 p2 0 1 0 p3 10.55 0 0 &size 5 2 10 4 dim 4 5.55 4 5.55 rat 1 1 1 1.2 group outsiderockxygen zone cshell p0 0 0 0 p1 0 0 -4 p2 0 1 0 p3 5.55 0 0 &size 1 2 10 4 dim 3.5 5.05 3.5 5.05 rat 1 1 1 1 group concretlinerxy fill group insiderockxy plot block groupplot add axes redgen zon brick p0 0 0 10.55 p1 10.55 0 10.55 p2 0 1 10.55 p3 0 0 25.55 &size 5 2 8 group outsiderock1gen zon brick p0 0 0 -34.55 p1 10.55 0 -34.55 p2 0 1 -34.55 p3 0 0 -9 &size 5 2 12 group outsiderock2gen zon brick p0 10.55 0 0 p1 50 0 0 p2 10.55 1 0 p3 10.55 0 10.55 &size 20 2 5 group outsiderock3gen zon brick p0 10.55 0 -9 p1 50 0 -9 p2 10.55 1 -9 p3 10.55 0 0 &size 20 2 5 group outsiderock3gen zon brick p0 10.55 0 10.55 p1 50 0 10.55 p2 10.55 1 10.55 p3 10.55 0 25.55 & size 20 2 8 group outsiderock4gen zon brick p0 10.55 0 -34.55 p1 50 0 -34.55 p2 10.55 1 -34.55 p3 10.55 0 -9 & size 20 2 12 group outsiderock5gen zon reflect dip 90 dd 90 orig 0 0 0plot set rotation 30 0 30 ;显示三维图，绕x轴负方向转30度set gravity 0 0 -10 ;设置重力加速度为z方向-10fix z range z -34.56,-34.54 ;设置底边界fix x range x -50.01,-49.99 ;设置左边界fix x range x 49.99,50.01 ;设置右边界fix y range y -0.01 0.01 ;设置前边界fix y range y 0.99,1.01 ;设置后边界model mohr ;莫尔～库仑模型ini density 2000 ;围岩的密度prop bulk=7.14e8 shear=3.333e9 fric=25 coh=2e5 tension=1.0e5Setp 6000 ;求解6000次Solve ;求解计算save Gsol.sav ;计算结果保存在Gsol.sav文件中plot cont zdisp ;绘制竖向位移场，如图6-14所示plot cont xdisp ;绘制竖向位移场，如图6-15所示plot cont szz ;绘制竖向应力场，如图6-16所示plot cont sxx ;绘制竖向应力场，如图6-17所示;1-左上半断面开挖rest Gsol.savplot block groupplot add axes redini xdis 0.0 ydis 0.0 zdis 0.0ini density 2200 range group outsiderocksy any group outsiderocksz anyprop bulk=1.923e9, shear=4.167e9 fric=35, coh=0.5e6 &tension=5e5, range group outsiderocksy any group outsiderocksz anymodel null range group insiderocksz any group concretlinersz anystep 2000save stepp1.savplot cont zdisprest stepp1.savsel shell id=1 range x -0.01 0.01 z -0.01 5.56 cyl end1 0 0 0 end2 0 1 0 rad 5.55sel shell id=1 range x -0.01 0.01 z -0.01 5.56sel shell id=1 prop iso=(21.0e9,0.2) thick=0.25step 2000plot cont zdispsave stepp2.sav;3-左下半断面开挖rest stepp2.savplot block groupplot add axes redmodel null range group insiderockxz any group concretlinerxz anystep 2000save stepp3.savplot cont zdisprest stepp3.savsel shell id=1 range x -0.01 0.01 z -4.01 0.01sel shell id=1 prop iso=(21.0e9,0.2) thick=0.25step 2000plot cont zdispsave stepp4.sav;5-右上半断面开挖rest stepp4.savplot block groupplot add axes redmodel null range group insiderocksy any group concretlinersy anystep 2000save stepp5.savplot cont zdisprest stepp5.savsel shell id=1 range x -0.01 5.56 z -0.01 5.56 cyl end1 0 0 0 end2 0 1 0 rad 5.55 sel shell id=1 prop iso=(21.0e9,0.2) thick=0.25step 2000plot cont zdispsave stepp6.sav;7-右下半断面开挖rest stepp6.savplot block groupplot add axes redmodel null range group insiderockxy any group concretlinerxy anystep 2000save stepp7.savplot cont zdisprest stepp7.savsel delete shell id=1 range x -0.01 0.01sel shell id=1 prop iso=(21.0e9,0.2) thick=0.25step 2000save stepp8.savrest stepp8.savmodel elas range group concretlinerxymodel elas range group concretlinerxzmodel elas range group concretlinersymodel elas range group concretlinerszini density 2600 range group concretlinerxyini density 2600 range group concretlinerxzini density 2600 range group concretlinersyini density 2600 range group concretlinerszprop bulk=12.5e9, shear=16.667e9 range group concretlinerxyprop bulk=12.5e9, shear=16.667e9 range group concretlinerxzprop bulk=12.5e9, shear=16.667e9 range group concretlinersyprop bulk=12.5e9, shear=16.667e9 range group concretlinerszstep 2000plot cont zdispsave stepp9.savplot block groupplot add axes red6.3 IV级围岩施工过程模拟new;1/4 上半部分gen zon radcyl p0 0 0 0 p1 10.55 0 0 p2 0 1 0 p3 0 0 10.55 size 5 2 10 4&dim 5.55 5.55 5.55 5.55 rat 1 1 1 1.2 group outsiderocksygen zone cshell p0 0 0 0 p1 5.55 0 0 p2 0 1 0 p3 0 0 5.55 size 1 2 10 4&dim 5.15 5.15 5.15 5.15 rat 1 1 1 1 group concretlinersy fill group insiderocksy;1/4 下半部分gen zon radcyl p0 0 0 0 p1 0 0 -9 p2 0 1 0 p3 10.55 0 0 &size 5 2 10 4 dim 4 5.55 4 5.55 rat 1 1 1 1.2 group outsiderockxygen zone cshell p0 0 0 0 p1 0 0 -4 p2 0 1 0 p3 5.55 0 0 size 1 2 10 4&dim 3.6 5.15 3.6 5.15 rat 1 1 1 1 group concretlinerxy fill group insiderockxy;1/2绘制上下地层网格gen zon brick p0 0 0 10.55 p1 10.55 0 10.55 p2 0 1 10.55 p3 0 0 25.55 size 5 2 8 group outsiderock1gen zon brick p0 0 0 -34.55 p1 10.55 0 -34.55 p2 0 1 -34.55 p3 0 0 -9 size 5 2 12 group outsiderock2;1/2右侧中间部分土体网格gen zon brick p0 10.55 0 0 p1 50 0 0 p2 10.55 1 0 p3 10.55 0 10.55 size 20 2 5 group outsiderock3gen zon brick p0 10.55 0 -9 p1 50 0 -9 p2 10.55 1 -9 p3 10.55 0 0 size 20 2 5 group outsiderock3。

地铁车站FLAC命令流newtitle桩锚支护基坑初始应力模拟;调入模型重新指定分组impgrid 10mmodel.flac3Dgenerate merge 0.02group 粉质粘土 ran group 3group 杂填土 ran group 1group 粉土 ran group 2group 中砂 ran group 4group 粉粘 ran group 5group 开挖1 ran group 11group 开挖2 ran group 10group 开挖3 ran group 9group 开挖4 ran group 8group 开挖5 ran group 7group 开挖6 ran group 6;指定材料模型model elasdefine mqb1;杂填土me1=7.8e7mv1=0.35;mbulk1=me1/(3.0*(1-2.0*mv1)) mshear1=me1/(2.0*(1+mv1));粉土me2=6.36e7mv2=0.4;mbulk2=me2/(3.0*(1-2.0*mv2))mshear2=me2/(2.0*(1+mv2));粉粘上me3=8e7mv3=0.35;mbulk3=me3/(3.0*(1-2.0*mv3))mshear3=me3/(2.0*(1+mv3));中砂me4=20e7mv4=0.4;mbulk4=me4/(3.0*(1-2.0*mv4))mshear4=me4/(2.0*(1+mv4));粉粘下me5=9.74e7mv5=0.3;mbulk5=me5/(3.0*(1-2.0*mv5))mshear5=me5/(2.0*(1+mv5))endmqb1;定义材料参数property density 1850 bulk mbulk1 shear mshear1 ran group 杂填土 any group 开挖1 anyproperty density 1950 bulk mbulk2 shear mshear2 ran group 粉土 any group 开挖2 any粉质粘土 any group 开挖3 anyproperty density 1980 bulk mbulk3 shear mshear3 ran group 开挖4 any group 开挖5 anyproperty density 1980 bulk mbulk3 shear mshear3 ran group 开挖6property density 2000 bulk mbulk4 shear mshear4 ran group 中砂property density 1970 bulk mbulk5 shear mshear5 ran group 粉粘;定义边界条件fix x y z ran z 0fix x ran x 0fix x ran x 140fix y;初始应力状态set gravity 0 0 -10history unbalsolve;塑性计算model mohrproperty density 1850 bulk mbulk1 shear mshear1 friction 15 cohesion 1e4 tension 0 ran group 杂填土 any group 开挖1 any property density 1950 bulk mbulk2 shear mshear2 friction 28 cohesion 2.46e4 tension 0 ran group 粉土 any group 开挖2 any property density 2000 bulk mbulk4 shear mshear4 friction 35 cohesion 0 tension 0 ran group 中砂property density 1970 bulk mbulk5 shear mshear5 friction 19.2 cohesion 2.41e4 tension 0 ran group 粉粘16.3 cohesion 2.31e4 tension 0 ran group 粉质粘土 any group 开挖3 anyproperty density 1980 bulk mbulk3 shear mshear3 friction 16.3 cohesion 2.31e4 tension 0 ran group 开挖4 any group 开挖5 anyproperty density 1980 bulk mbulk3 shear mshear3 friction 16.3 cohesion 2.31e4 tension 0 ran group 开挖6solvesave initial.sav;调入初始应力计算结果rest initial.savtitle基坑分步开挖支护效果set log onset logfile jisuanguocheng.loginitial xdisp 0 ydisp 0 zdisp 0initial xvel 0 yvel 0 zvel 0initial sta 0;;;;;;;;;开挖;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;设置监控变量;竖向位移hist id=2 gp zdisp 93.1 5 50hist id=3 gp zdisp 89.1 5 50hist id=4 gp zdisp 84.1 5 50hist id=5 gp zdisp 79.1 5 50histid=6 gp zdisp 74.1 5 50hist id=7 gp zdisp 69.1 5 50hist id=8 gp zdisp 59.1 5 50hist id=9 gp zdisp 49.1 5 50hist id=10 gp zdisp 39.1 5 50hist id=11 gp zdisp 29.1 5 50;桩的水平位移hist id=12 gp xdisp 93.1 5 50hist id=13 gp xdisp 100 0 43hist id=14 gp xdisp 100 0 41hist id=15 gp xdisp 100 0 39hist id=16 gp xdisp 100 0 37hist id=17 gp xdisp 100 0 35hist id=18 gp xdisp 100 0 33hist id=19 gp xdisp 100 0 31hist id=20 gp xdisp 100 0 29hist id=21 gp xdisp 100 0 27hist id=22 gp xdisp 100 0 25hist id=23 gp xdisp 100 0 23hist id=24 gp xdisp 100 0 21 ;;;;;;;;第一步开挖;;;;;;;;;;;model null ran group 开挖1;安装第一道土钉;定义土钉相关参数def canshu;---可以修改的量----_mLength=9.0 ;锚杆1,2长度_mLength3=6.0 ;锚杆3长度_mAngle=5.0 ;锚杆与水平方向的夹角_msInterval=2.0 ;锚杆水平方向的间隔_mzInterval=1.5 ;锚杆纵向的间隔_mnseg=9 ;12锚杆划分单元数_mnseg3=6 ;3锚杆划分单元数_zm1t=48 ;第一排锚杆顶点高度，每排锚杆是在同一高度的。

隧道及地下工程FLAC解析方法命令流1.确认工程文件：首先，确认需要解析的FLAC文件和相关的工程文件是否齐全，并保存在同一个文件夹中。

2. 安装解析软件：在电脑上安装合适的FLAC解析软件，比如Free Lossless Audio Codec（FLAC）或者Exact Audio Copy（EAC）等。

3.打开软件：运行安装好的FLAC解析软件，点击“打开文件”或者“导入文件”按钮，选择需要解析的FLAC文件。

4.设置解析参数：在打开的软件界面中，可以设置一些解析参数，比如编码格式、采样率、音频通道等。

根据不同的需求来进行设置。

推荐采用默认参数进行解析，以保持最高的音质。

5.开始解析：点击“开始解析”按钮，软件将开始解析FLAC文件。

这个过程可能需要几分钟甚至更长的时间，具体的解析速度取决于电脑的配置和FLAC文件的大小。

6. 验证解析结果：解析完成后，软件将生成一个解析后的音频文件，一般是以WAV格式保存。

使用音频播放软件，比如Windows MediaPlayer或者VLC media player等，打开解析后的音频文件，验证解析结果是否符合预期。

7.转换为其他格式（可选）：如果需要将解析后的音频文件转换为其他格式，比如MP3或者AAC等，可以使用FLAC解析软件内置的转换功能，选择对应的音频输出格式，然后再点击“开始转换”按钮进行转换。

8.导出解析结果：最后，将解析后的音频文件导出到指定的文件夹中。

可以选择保存为原文件名，或者自定义新的文件名。

以上是一个FLAC解析方法的命令流示例。

根据不同的软件和需求，具体的步骤可能会有所不同。

在进行FLAC解析时，建议先了解所使用软件的具体功能和操作方法，并根据实际情况进行相应的设置和调整。