《隧道及地下工程FLAC解析方法》全部命令流汇总.docx
《《隧道及地下工程FLAC解析方法》全部命令流汇总.docx》由会员分享，可在线阅读，更多相关《《隧道及地下工程FLAC解析方法》全部命令流汇总.docx（60页珍藏版）》请在冰点文库上搜索。
《隧道及地下工程FLAC解析方法》全部命令流汇总
第二章FLAC3D原理及入门指南
2-1定义一个FISH函数
new
defabc
abc=25*3+5
End
printabc
2-2使用一个变量
new
defabc
hh=25
abc=hh*3+5
End
Printhh
Printabc
2-3对变量和函数的理解
new
defabc
hh=25
abc=hh*3+5
End
setabc=0hh=0
printhh
printabc
printhh
new
defabc
abc=hh*3+5
end
sethh=25
printabc
setabc=0hh=0
printhh
printabc
printhh
2-4获取变量的历史记录
new
genzonebricksize121
modelmohr
propshear=1e8bulk=2e8cohes=1e5tens=1e10
fixxyzrangey-0.10.1
applyyvel-1e-5rangey1.92.1
plotsetrotation0045
plotblockgroup
defget_ad
ad1=gp_near（0,2,0）
ad2=gp_near（1,2,0）
ad3=gp_near（0,2,1）
ad4=gp_near（1,2,1）
end
get_ad
defload
load=gp_yfunbal（ad1）+gp_yfunbal（ad2）+gp_yfunbal（ad3）+gp_yfunbal（ad4）
end
histload
histgpydis0,2,0
step1000
plothis1vs-2
2-5用FISH函数计算体积模量和剪砌模量
new
defderive
s_mod=y_mod/（2.0*（1.0+p_ratio））
b_mod=y_mod/（3.0*（1.0-2.0*p_ratio））
end
sety_mod=5e8p_ratio=0.25
derive
printb_mod
prints_mod
2-6在FLAC输入中使用符号变量
New
defderive
s_mod=y_mod/（2.0*（1.0+p_ratio））
b_mod=y_mod/（3.0*（1.0-2.0*p_ratio））
end
sety_mod=5e8p_ratio=0.25
derive
genzonebricksize2,2,2
modelelastic
propbulk=b_modshear=s_mod
printzonepropbulk
printzonepropshear
2-7控制循环
New
defxxx
sum=0
prod=1
loopn（1,10）
sum=sum+n
prod=prod*n
end_loop
end
xxx
printsum,prod
new
genzonebrickp0（0,0,0）p1（-10,0,0）p2（0,10,0）p3（0,0,-10）
modelelas
plotsetrotation0045
plotblockgroup
definstall
pnt=zone_head
loopwhilepnt#null
z_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_loop
end
setp_ratio=0.25y_zero=1e7cc=1e8
install
2-8拆分命令行
new;exampleofasumofmanythings
deflong_sum
temp=v1+v2+v3+v4+v5+v6+v7+v8+v9+v10
long_sum=temp+v11+v12+v13+v14+v15
end
2-9变量类型
new
defhaveone
aa=2
bb=3.4
cc=’Haveaniceday’
dd=aa*bb
ee=cc+’,oldchap’
end
haveone
printfish
2-10IF条件语句
new
defabc
ifxx>0then
abc=1000
else
abc=-1000
end_if
end
setxx=10
printabc
setxx=0
printabc
2-11索单元自动生成
new
genzonebricksize1035
plotsetrotation0045
plotblockgroup
defplace_cables
loopn（1,5）
z_d=float（n）-0.5
command
selcablebeg0.0,1.5,z_dend7.0,1.5,z_dnseg7
end_command
end_loop
end
place_cables
plotgridselgeomred
new
genzonebricksize1035
plotsetrotation15060
plotblockgroup
modmohr
propbulk1e8shear.3e8fric35
propcoh1e3tens1e3
inidens1000
setgrav0,0,-10
fixxyzrangez-.1.1
fixyrangey-.1.1
fixyrangey2.93.1
fixxrangex-.1.1
fixxrangex9.910.1
setlarge
histunbal
solve
savecab_str.sav
inixdis0ydis0zdis0
histgpxdisp0,1,5
defplace_cables
loopn（1,5）
z_d=5.5-float（n）
z_t=z_d+0.5
z_b=z_d-0.5
command
freexrangex-.1,.1zz_bz_t
solve
selcablebeg0.0,0.5,z_dend7.0,0.5,z_dnseg7
selcablebeg0.0,1.5,z_dend7.0,1.5,z_dnseg7
selcablebeg0.0,2.5,z_dend7.0,2.5,z_dnseg7
selcablepropemod2e10ytension1e8xcarea1.0&
gr_k2e10gr_coh1e10gr_per1.0
end_command
end_loop
end
place_cables
savecab_end.sav
plotsketchselcableforcered
2-12圆形隧道开挖模拟计算
;建立模型
genzonradcylp0000p1600p2010p3006&
size4284dim3333rat1111.2groupoutsiderock
genzonecshellp0000p1300p2010p3003&
size1284dim2.72.72.72.7rat1111groupconcretlinerfillgroupinsiderock
genzonreflectdip90dd90orig000
genzonreflectdip0dd0ori000
genzonbrickp0006p1606p2016p30013size426groupoutsiderock1
genzonbrickp000-12p160-12p201-12p300-6size425groupoutsiderock2
genzonbrickp0600p12100p2610p3606size1024groupoutsiderock3
genzonreflectdip0dd0orig000rangegroupoutsiderock3
genzonbrickp0606p12106p2616p36013size1026groupoutsiderock4
genzonbrickp060-12p1210-12p261-12p360-6size1025groupoutsiderock5
genzonreflectdip90dd90orig000rangex-0.16.1z6.113.1
genzonreflectdip90dd90orig000rangex-0.16.1z-6.1-12.1
genzonreflectdip90dd90orig000rangex6.121.1z-12.113.1
;绘制模型图
plotblockgroup
plotaddaxesred
;plotsetrotation0045用于显示三维模型
;设置重力
setgravity00-10
;给定边界条件
fixzrangez-12.01,-11.99
fixxrangex-21.01,-20.99
fixxrangex20.99,21.01
fixyrangey-0.010.01
fixyrangey0.99,1.01
;求解自重应力场
modelmohr
inidensity1800;围岩的密度
propbulk=1.47e8shear=5.6e7fric=20coh=5.0e4tension=1.0e4;体积、剪切、摩擦角、凝聚力、抗拉强度
setmechratio=1e-4
solve
saveGravsol.sav
plotcontzdispoutlon
plotcontszz
;毛洞开挖计算
initialxdisp=0ydisp=0zdisp=0
modelnullrangegroupinsiderockanygroupconcretlinerany
plotblockgroup
plotaddaxesred
setmechratio=5e-4
solve
saveKaiwsol.sav
plotcontzdisp
plotcontsdisp
plotcontszz
plotcontxzz
;模筑衬砌计算
modelelasrangegroupconcretlinerany
plotblockgroup
plotaddaxesred
inidensity2500rangegroupconcretlinerany;衬砌混凝土的密度
propbulk=16.67e9,shear=12.5e9rangegroupconcretlinerany;衬砌混凝土的体积弹模、剪切弹模
setmechratio=1e-4
solve
savezhihusol.sav
plotcontzdisp
plotcontsdisp
plotcontszz
plotcontxzz
;完成计算分析
第四章FISH语言及建模技术
4-1数组
new
defafill;fillmatrixwithrandomnumbers
arrayvar（4,3）
loopm（1,4）
loopn（1,3）
var（m,n）=urand
endloop
endloop
end
defashow;displaycontentsofmatrix
loopm（1,4）
hed=’’
msg=’’+string（m）
loopn（1,3）
hed=hed+’’+string（n）
msg=msg+’’+string（var（m,n））
endloop
ifm=1
dum=out（hed）
endif
dum=out（msg）
endloop
end
afill
Ashow
4-2函数操作
new
defxxx
aa=2*3
xxx=aa+bb
end
new
defstress_sum
stress_sum=0.0
pnt=zonehead
loopwhilepnt#null
Stress_sum=stress_sum+z_sxx（pnt）
pnt=z_next（pnt）
end_loop
end
new
defstress_sum
sum=0.0
pnt=zonehead
loopwhilepnt#null
sum=sum+z_sxx（pnt）
pnt=znext（pnt）
endloop
stress_sum=sum
end
new
defh_var_1
ipz=z_near（1,2,3）
H_var_1=sxx（ipz）+fac
H_var_2=syy（ipz）+fac
H_var_3=szz（ipz）+fac
H_var_4=sxy（ipz）+fac
H_var_5=syz（ipz）+fac
H_var_6=sxz（ipz）+fac
end
4-3函数删除与重定义
new
defjoe
ii=out（‘ThisisAFunction’）
end
deffred
joe
end
fred
defjoe
ii=（‘ThisisANewFunction’）
end
fred
4-4字符串
new
defin_def
xx=in（msg+’（’+’default:
’+string（default）+’）:
’）
iftype（xx）=3
in_def=default
else
in_def=xx
endif
end
;
defmoduli_data
default=1.0e9
msg=’InputYoung‘smodulus’
Y_mod=in_def
default=0.25
msg=’InputPoisson‘sratio’
p_ratio=in_def
ifp_ratio=0.5then
ii=out（’BulkmodisundefinedatPoisson‘sratio=0.5’）
ii=out（’Selectadifferentvalue--’）
p_ratio=in_def
endif
s_mod=y_mod/（2.0*（1.0+p_ratio））
b_mod=y_mod/（3.0*（1.0-2.0*p_ratio））
end
;
moduli_data
genzonebricksize2,2,2
modelelastic
propbulk=b_modshear=s_mod
printp_ratioy_modb_mods_mod
pause
printzonepropbulk
pause
printzonepropshear
4-5马蹄形隧道网格
new
genzoneradcylsize510612rat1111.2&
p00,0,0p1100,0,0p20,200,0p30,0,100
genzoneradtunsize510512rat1111.2&
p00,0,0p10,0,-100p20,200,0p3100,0,0
genzonereflectdip90dd270origin0,0,0
plotsetrotation30030
plotblockgroup
4-6复杂形状网格生成
new
;主隧道网格
genzonradcylp01500p12300p215500p31508&
size41064dim4444rat1111fill
genzonreflectdip90dd90orig1500
genzonreflectdip0ori000
;辅助隧道网格
genpointid1（2.969848,0.0,-0.575736）
genpointid2（2.969848,50.0,-0.575736）
genzonradcylp000-1p1700p2050-1p3008p47500&
p50508p6708p77508p8point1p10point2&
size31064dim3333rat1111
genzonradcylp000-1p100-8p2050-1p3700p4050-8&
p57500p670-8p7750-8p9point1p11point2&
size31064dim3333rat1111
;衬砌网格
selshellrangecylend100-1end2050-1rad3
;隧道外围边界网格
genzoneradtunp0700p15000p27500p315050p450500&
p5155050p650050p7505050&
p82300p9708p1023500p117508&
size610310rat1111.1
genzonebrickp0008p1708p20508p30050&
p47508p505050p615050p7155050&
size31010rat111.1
;对称得出1/2模型
genzonreflectdip0ori000rangex023y050z850
genzonreflectdip0ori000rangex2350y050z050
;建立主隧道和辅助隧道块名称
groupservicerangecylend100-1end2050-1rad3
groupmainrangecylend11500end215500rad4
;对称得出完整模型
genzonreflectdip90dd270ori000
;网格显示
plotsetrotation30030
plotblockgroup
4-7网格连接
new
genzonebricksize444p00,0,0p14,0,0p20,4,0p30,0,2
genzonebricksize884p00,0,2p14,0,2p20,4,2p30,0,4
attachfacerangez1.92.1
modelelas
propbulk8e9shear5e9
fixzrangez-.1.1
fixxrangex-.1.1
fixxrangex3.94.1
fixyrangey-.1.1
fixyrangey3.94.1
applyszz-1e6rangez3.94.1x0,2y0,2
histunbal
;网格显示
plotsetrotation30030
plotblockgroup
solve
saveatt.sav
plotcontzdispoutlon
单一网格（小）程序
new
genzonebricksize888p00,0,0p14,0,0p20,4,0p30,0,4
modelelas
propbulk8e9shear5e9
fixzrangez-.1.1
fixxrangex-.1.1
fixxrangex3.94.1
fixyrangey-.1.1
fixyrangey3.94.1
applyszz-1e6rangez3.94.1x0,2y0,2
histunbal
;网格显示
plotsetrotation30030
plotblockgroup
solve
savenoattx.sav
plotcontzdispoutlon
单一网格（大）程序
new
genzonebricksize444p00,0
升级会员 