隧道及地下工程FLAC解析方法 命令流Word格式文档下载.docx

隧道及地下工程FLAC解析方法 命令流Word格式文档下载.docx

《隧道及地下工程FLAC解析方法 命令流Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《隧道及地下工程FLAC解析方法 命令流Word格式文档下载.docx（51页珍藏版）》请在冰豆网上搜索。

defload

load=gp_yfunbal（ad1）+gp_yfunbal（ad2）+gp_yfunbal（ad3）+gp_yfunbal（ad4）

histload

histgpydis0,2,0

step1000

plothis1vs-2

2-5用FISH函数计算体积模量和剪砌模量

defderive

s_mod=y_mod/（2.0*（1.0+p_ratio））

b_mod=y_mod/（3.0*（1.0-2.0*p_ratio））

sety_mod=5e8p_ratio=0.25

derive

printb_mod

prints_mod

2-6在FLAC输入中使用符号变量

New

genzonebricksize2,2,2

modelelastic

propbulk=b_modshear=s_mod

printzonepropbulk

printzonepropshear

2-7控制循环

defxxx

sum=0

prod=1

loopn（1,10）

sum=sum+n

prod=prod*n

end_loop

xxx

printsum,prod

genzonebrickp0（0,0,0）p1（-10,0,0）p2（0,10,0）p3（0,0,-10）

modelelas

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）

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

2-9变量类型

defhaveone

aa=2

bb=3.4

cc=’Haveaniceday’

dd=aa*bb

ee=cc+’,oldchap’

haveone

printfish

2-10IF条件语句

ifxx>

0then

abc=1000

else

abc=-1000

end_if

setxx=10

setxx=0

2-11索单元自动生成

genzonebricksize1035

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

place_cables

plotgridselgeomred

plotsetrotation15060

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

z_d=5.5-float（n）

z_t=z_d+0.5

z_b=z_d-0.5

freexrangex-.1,.1zz_bz_t

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

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

绘制模型图

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

saveGravsol.sav

plotcontzdispoutlon

plotcontszz

毛洞开挖计算

initialxdisp=0ydisp=0zdisp=0

modelnullrangegroupinsiderockanygroupconcretlinerany

setmechratio=5e-4

saveKaiwsol.sav

plotcontzdisp

plotcontsdisp

plotcontxzz

模筑衬砌计算

modelelasrangegroupconcretlinerany

inidensity2500rangegroupconcretlinerany;

衬砌混凝土的密度

propbulk=16.67e9,shear=12.5e9rangegroupconcretlinerany;

衬砌混凝土的体积弹模、剪切弹模

savezhihusol.sav

完成计算分析

第四章命令流按照顺序进行

4-1数组

defafill;

fillmatrixwithrandomnumbers

arrayvar（4,3）

loopm（1,4）

loopn（1,3）

var（m,n）=urand

endloop

defashow;

displaycontentsofmatrix

hed=’’

msg=’’+string（m）

hed=hed+’’+string（n）

msg=msg+’’+string（var（m,n））

ifm=1

dum=out（hed）

endif

dum=out（msg）

afill

Ashow

4-2函数操作

aa=2*3

xxx=aa+bb

defstress_sum

stress_sum=0.0

pnt=zonehead

loopwhilepnt#null

Stress_sum=stress_sum+z_sxx（pnt）

pnt=z_next（pnt）

sum=0.0

sum=sum+z_sxx（pnt）

pnt=znext（pnt）

endloop

stress_sum=sum

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

4-3函数删除与重定义

defjoe

ii=out（‘ThisisAFunction’）

deffred

joe

fred

ii=（‘ThisisANewFunction’）

4-4字符串

defin_def

xx=in（msg+’（’+’default:

’+string（default）+’）:

’）

iftype（xx）=3

in_def=default

in_def=xx

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--’）

moduli_data

propbulk=b_modshear=s_mod

printp_ratioy_modb_mods_mod

pause

4-5马蹄形隧道网格

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

4-6复杂形状网格生成

主隧道网格

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&

衬砌网格

selshellrangecylend100-1end2050-1rad3

隧道外围边界网格

genzoneradtunp0700p15000p27500p315050p450500&

p5155050p650050p7505050&

p82300p9708p1023500p117508&

size610310rat1111.1

genzonebrickp0008p1708p20508p30050&

p47508p505050p615050p7155050&

size31010rat111.1

对称得出1/2模型

genzonreflectdip0ori000rangex023y050z850

genzonreflectdip0ori000rangex2350y050z050

建立主隧道和辅助隧道块名称

groupservicerangecylend100-1end2050-1rad3

groupmainrangecylend11500end215500rad4

对称得出完整模型

genzonreflectdip90dd270ori000

网格显示

4-7网格连接

genzonebricksize444p00,0,0p14,0,0p20,4,0p30,0,2

genzonebricksize884p00,0,2p14,0,2p20,4,2p30,0,4

attachfacerangez1.92.1

propbulk8e9shear5e9

fixzrangez-.1.1

fixxrangex3.94.1

fixyrangey3.94.1

applyszz-1e6rangez3.94.1x0,2y0,2

plotblockgroup

saveatt.sav

单一网格（小）程序

genzonebricksize888p00,0,0p14,0,0p20,4,0p30,0,4

savenoattx.sav

单一网格（大）程序

genzonebricksize444p00,0,0p14,0,0p20,4,0p30,0