辽工大FLAC考试资料

FLAC

n;（新建）

pl bl gr;（显示组信息）

pl sh;（画出显示）

;1—建模

gen zone radcyl p0 0 0 0 p1 add 50 0 0 p2 add 0 50 0 p3 add 0 0 50 &;（gen生成，zone区域，定义起始点和X、Y、Z方向长度，add代表在P0坐标基础上加）

dim 1.75 1.75 1.75 1.75 ratio 1.0 1.0 1.0 1.5 size 8 25 8 12 fill;（确定隧洞，前两个1.75代表x=0面上X、Y方向半径是1.75，后两个1.75代表X=50的面上X、Y方向半径是1.75；ratio 表示网格长度比例，顺序为X、Y、Z和围岩径向；size表示网格个数，顺序为径向、轴向、周向和围岩径向；fill表示填充）

gen zone reflect orig 0 0 0 norm 0 0 1;（沿过(0.0.0)法线为(0.0.1)的面做对称，reflect是对称）gen zone brick p0 0 0 50 p1 add 50 0 0 p2 add 0 50 0 p3 add 0 0 250 &;（brick块体）

ratio 1.0 1.0 1.5 size 4 25 4

;2—指定本构模型及参数mm

model moh;（本构关系：模量摩尔库伦准则）

def chushi;（定义初始）

jian=yang/(2*(1+passion));（剪切模量表达式=杨氏模量/(2*(1+泊松比))）

tiji=yang/(3*(1-2*passion));（体积模量表达式=杨氏模量/(3*(1-2*泊松比))）

end

def cailiao1;（定义材料1）

command;（命令开始）

set yang=3.2e9

set passion=0.25

chushi;（调用初始）

print jian;（打印即屏幕输出）

print tiji

prop shear jian bulk tiji ran z -4,4;（-4到4之间为同种材料，prop是定义，shear是剪切模量，bulk是体积模量，ran是范围）

prop coh=2.3e6 fric=45.0 tens=0.5e6 ran z -4,4;（coh是粘结力，fric是内摩擦角，tens是抗拉强度）

ini den 0.0025e6 ran z -4,4;（ini是初始，den是密度）

endcommand;（命令结束）

end

cailiao1;（调用材料1）

def cailiao2

command

set yang=4.8e9

set passion=0.24

chushi

print jian

print tiji

prop shear jian bulk tiji ran z 4,50

prop shear jian bulk tiji ran z -4,-50

prop coh=2.8e6 fric=48.0 tens=0.7e6 ran z 4,50

prop coh=2.8e6 fric=48.0 tens=0.7e6 ran z -4,-50

ini den 0.0025e6 ran z 4,50

ini den 0.0025e6 ran z -4,-50

endcommand

end

cailiao2

def cailiao3

command

set yang=5.5e9

set passion=0.23

Chushi

print jian

print tiji

prop shear jian bulk tiji ran z 50,30

prop coh=3.2e6 fric=50.0 tens=0.95e6 ran z 50,300

ini den 0.0025e6 ran z 50,300

endcommand

end

cailiao3

;3—施加边界条件

set grav 0 0 -10;（施加重力）

ini szz=-7.5e6 grad 0,0,2.5e4;（Z方向应力=szz+容重*z，Z方向应力为地表应力得0，容重=密度*g，所以0=szz+0.0025e6*(-10)*300，szz=-7.5e6，grad是梯度，梯度=容重）

ini sxx=-11.25e6 grad 0,0,3.75e4;（7.5*1.5=11.25 2.5*1.5=3.75，定义初始应力场lamdx=1.5 lamdy=1.25 为侧压力系数）

ini syy=-9.375e6 grad 0,0,3.125e4;（7.5*1.25=9.375 2.5*1.25=3.125）

fix x range x -0.1 0.1;（施加位移、速度场边界条件）

fix x range x 49.9 50.1

fix y range y -0.1 0.1

fix y range y 49.9 50.1

fix x y z range z -50.1 -49.9;（固定z=-50面上X、Y、Z方向的速度）

;4—生成初始/原岩地应力场

sol;（sol是求解）

sav 05-000.sav;（sav是保存，05-000.sav是文件名，后缀为.sav）

;5—开挖支护前模型位移场、速度场清零

ini xdisp 0 ydisp 0 zdisp 0 xvel 0 yvel 0 zvel 0;(xdisp表示X方向位移为0，xvel表示X方向速度为0)

;6—开挖过程的模拟

model null range cyl end1 0,0,0 end2 0 4 0 radius 1.75;（开挖，cyl表示圆柱形轴线起点（0,0,0）挖到（0,4,0），radius表示半径）

hist gp zdisp 0 2 1.751;（拱顶竖向位移，监测深2米，高1.751米处的围岩情况，hist表示监测，gp表示节点）

hist gp zvel 0 2 1.751;（拱顶竖向速度）