陈育民《FLAC3D基础与工程实例》全部命令流.doc

《FLAC/FLAC 3D基础与工程实例》全部命令流

1
gen 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_1
gen 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_2
gen zone bri p0 40 0 0 edge 10 size 10 10 10 rat 1.0 0.9 1.1 group brick_3
plot sur

2-1
new
gen zon bri size 3 3 3
model elas
prop bulk 3e8 shear 1e8
ini dens 2000
fix z ran z -.1 .1
fix x ran x -.1 .1
fix x ran x 2.9 3.1
fix y ran y -.1 .1
fix y ran y 2.9 3.1
set grav 0 0 -10
solve
app nstress -10e4 ran z 3 x 1 2 y 1 2
hist gp vel 0 0 3
hist gp vel 0 3 3
plo hist 1 red
plo add hist 2 blue
solve

3-1
;-------------------------------------工程信息
;Project Record Tree export
;Title:Simple test
;---------------------------------计算第一步
;... STATE: STATE1 ....
config
grid 10,10
model elastic
group 'User:Soil' notnull
model elastic notnull group 'User:Soil'
prop density=1500.0 bulk=3E6 shear=1E6 notnull group 'User:Soil'
fix x y j 1
fix x i 1
fix x i 11
set gravity=9.81
history 999 unbalanced
solve
save state1.sav
;----------------------------------计算第二步
;... STATE: STATE2 ....
initial xdisp 0 ydisp 0
initial xvel 0 yvel 0
model null i 4 7 j 8 10
group 'null' i 4 7 j 8 10
group delete 'null'
history 1 xdisp i=4, j=11
solve
save state2.sav
;--------------------------------绘图命令
;*** plot commands ****
;plot name: syy
plot hold grid syy fill
;plot name: Unbalanced force
plot hold history 999
;plot name: grid
plot hold grid magnify 20.0 lred grid displacement
;plot name: Xdis-A
plot hold history 1 line

5-1
new
; ===============================
; 定义球体半径和半径方向上单元网格数
; ===============================
def parm
rad=10.0
rad_size=5
end
parm
; ===============================
; 建立八分之一球体外接立方体网格
; ===============================
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 1
gen 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 2
gen 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_sphere
p_gp=gp_head
loop 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*px
gp_ypos(p_gp)=k*py
gp_zpos(p_gp)=k*pz
end_if
p_gp=gp_next(p_gp)
end_loop
end
make_sphere
; ===============================
; 利用镜像生成完整球体网格
; ===============================
gen zone ref
gen zone ref dip 90
gen zone ref dip 90 dd 90
; ===============================
; 显示球体网格
; ===============================
plot surf
pl set back wh
pl bl gr

5-3
n
gen zon bri size 1 1 2
group soil ran z 1 1
group rock ran z 0 1
expgrid 1.flac3d

6-1
gen zon bri size 3 3 3
model mohr
prop bu 3e6 sh 1e6 coh 10e3 fric 15
fix z ran z -.1 .1
fix x ran x -.1 .1
fix x ran x 2.9 3.1
fix y ran y -.1 .1
fix y ran y 2.9 3.1
ini dens 2000
hist unbal
set grav 10
solve elastic
save 6-1.sav

6-2
rest 6-1.sav
ini 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
solve
save 6-2.sav

6-3
rest 6-1.sav
ini 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
hist id=2 gp zdis 1 1 3
hist id=3 gp zdis 1 1 2
hist id=4 gp xdis 1 1 3
hist id=5 gp xdis 1 1 3
hist id=6 zone szz 1 1 3
hist id=7 zone szz 1.5 1.5 2.5
hist id=8 zone sxz 1.5 1.5 2.5
solve
save 6-3.sav

6-4
rest 6-3.sav
set log on
set logfile 6-2.log
print zone stress
print gp dis
set log off

6-5
rest 6-1.sav
ini 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
plot set rot 20 0 30
plot con szz ou on magf 10
plot add hist 1
set movie avi step 1 file 6-5.avi
movie start
solve
movie finish

7-1
n
gen zon bri size 1 1 2
model elas
prop bulk 3e7 shear 1e7
fix z ran z 0
fix x ran x 0
fix x ran x 1
fix y ran y 0
fix y ran y 1
ini dens 2000
set grav 0 0 -10
solve
plo con sz

7-2
n
gen zon bri size 1 1 2
model mohr
prop bulk 3e7 shear 1e7 c 1e10 f 15 tension 1e10
fix z ran z 0
fix x ran x 0
fix x ran x 1
fix y ran y 0
fix y ran y 1
ini dens 2000
set grav 0 0 -10
solve
prop bulk 3e7 shear 1e7 c 10e3 f 15 ten 0
solve
plo con sz

7-3
n
gen zone brick size 1 1 2
model mohr
prop bulk 3e7 shear 1e7 coh 10e3 fri 15 ten 0
fix z ran z 0
fix x ran x 0
fix x ran x 1
fix y ran y 0
fix y ran y 1
ini dens 2000
set grav 0 0 -10
solve elas
plo con sz

7-4
new
gen zone brick size 1 1 2
model mohr
prop bulk 3e7 shear 1e7 coh 10e3 fri 15 ten 0
fix z ran z 0
fix x ran x 0
fix x ran x 1
fix y ran y 0
fix y ran y 1
ini dens 2000
ini szz -40e3 grad 0 0 20e3 ran z 0 2
ini syy -20e3 grad 0 0 10e3 ran z 0 2
ini sxx -20e3 grad 0 0 10e3 ran z 0 2
set grav 0 0 -10
solve
plo con sz

7-5
n
gen zon bri size 1 1 2
model m
prop bulk 3e7 shear 1e7 c 10e10 f 15 ten 1e10
fix z ran z 0
fix x ran x 0
fix x ran x 1
fix y ran y 0
fix y ran y 1
ini dens 2000 ran z 0 1
ini dens 1500 ran z 1 2
ini szz -35e3 grad 0 0 20e3 ran z 0 1
ini syy -17.5e3 grad 0 0 10e3 ran z 0 1
ini sxx -17.5e3 grad 0 0 10e3 ran z 0 1
ini szz -15e3 grad 0 0 15e3 ran z 1 2
ini syy -7.5e3 grad 0 0 7.5e3 ran z 1 2
ini sxx -7.5e3 grad 0 0 7.5e3 ran z 1 2
ini pp 10e3 grad 0 0 -10e3 ran z 0 1
set grav 0 0 -10
solve
plo con sz

7-6
n
gen zon bri size 1 1 2
model m
prop bulk 3e7 shear 1e7 c 10e10 f 15 ten 1e10
fix z ran z 0
fix x ran x 0
fix x ran x 1
fix y ran y 0
fix y ran y 1
ini dens 2000 ran z 0 2
ini szz -50e3 grad 0 0 20e3 ran z 0 1
ini syy -30e3 grad 0 0 10e3 ran z 0 1
ini sxx -30e3 grad 0 0 10e3 ran z 0 1
ini pp 30e3 grad 0 0 -10e3 ran z 0 2
app nstress -10e3 ran z 2
set grav 0 0 -10
solve
plo con sz

7-7
new
gen 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 10
model elas
pro bulk 10e10 she 10e10
ini den 2500
apply 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 120
set grav 0 0 -10
step 30000
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
plo cont szz

7-8
new
gen 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 10
model elas
pro bulk 10e10 she 10e10
ini den 2500
ini sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1
ini sxx -1e9 grad 0 0 6.6666666e6 range x 59.9 60.1
ini syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1
ini syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1
ini szz -1e8 ran z -.1 .1
fix x y z ran z -.1 .1
set grav 0 0 -10
solve
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
plo cont szz

7-9
new
gen 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 clay
gen 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 soil
gen 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 soil
gen 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 dam

fix x y z ran z -49.9 -50.1
fix x ran x -.1 .1
fix x ran x 99.9 100.1
fix y

model mohr ran z -50 0
model null ran z 0 5
prop bulk 7.8e6 shear 3.0e6 coh 10e10 tension 1e10 ran group soil
ini dens 1500 ran group soil
prop bulk 3.91e6 shear 1.5e6 coh 10e10 tension 1e10 ran group clay
ini dens 1800 ran group clay

set grav 0 0 -9.8
hist id=1 unbal
solve
prop bulk 7.8e6 shear 3.0e6 coh 10e3 fric 15 ran group soil
prop bulk 3.91e6 shear 1.5e6 coh 20e3 fric 20 ran group clay
solve
save elastic.sav

ini 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 elastic ran z 0 1 ; ;激活0 m ~ 1 m的单元
prop bulk 7.8e6 shear 3.0e6 ran z 0 1
ini dens 1500 ran z 0 1
solve ;按软件默认精度求解
save fill-1.sav

model elastic ran z 1 2
prop bulk 7.8e6 shear 3.0e6 ran z 1 2
ini dens 1500 ran z 1 2
solve
save fill-2.sav

model elastic ran z 2 3
prop bulk 7.8e6 shear 3.0e6 ran z 2 3
ini dens 1500 ran z 2 3
solve
save fill-3.sav

model elastic ran z 3 4
prop bulk 7.8e6 shear 3.0e6 ran z 3 4
ini dens 1500 ran z 3 4
solve
save fill-4.sav

model elastic ran z 4 5
prop bulk 7.8e6 shear 3.0e6 ran z 4 5
ini dens 1500 ran z 4 5
solve
save fill-5.sav
pau
;plo bl gr
;gen zone brick p0 0 0 0 p1 100 0 0 p2 0 5 0 p3 0 0 5 size
gen 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 clay
gen 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 soil
gen 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 soil
gen 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 dam

7-10
set log on ;打开log记录
set logfile 1.log ;设置记录文件名为：1.log
restore fill-1.sav ;调用保存的文件
print gp dis range id 517 any id 533 any ;输出两个节点的变形值
restore fill-2.sav
print gp dis range id 517 any id 533 any
restore fill-3.sav
print gp dis range id 517 any id 533 any
restore fill-4.sav
print gp dis range id 517 any id 533 any
restore fill-5.sav
print gp dis range id 517 any id 533 any
set log off ;关闭log记录

8-1
def abc
abc = 1 + 2 * 3
abcd = 1.0 / 2.0
end
abc
print fish

8-2
def abc
if aa < 0 then
abc = 0.0
else
abc = 2.0 * aa
endif
end
abc

8-3
def abc
loop aa (1, 2.5)
command
print aa
endcommand
endloop
end

8-4
new
gen zon bri size 3 3 3
model elastic
prop bu 3e7 sh 1e7
ini dens 2000
fix x y z ran z -.1 .1
fix x ran x -.1 .1
fix x ran x 2.9 3.1
fix y ran y -.1 .1
fix y ran y 2.9 3.1
set grav 10
solve
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
save 8-4.sav
8-5
rest 8-4.sav
def E_modify
p_z = zone_head
d_k = 704
d_n = 0.38
d_pa = 101325.0 ;//标准大气压
loop while p_z # null
sigma_3 = -1.0 * z_sig1(p_z)
E_new = d_k * d_pa * (sigma_3 / d_pa) ^ d_n
z_prop(p_z,'young') = E_new
p_z = z_next(p_z)
endloop
end
E_modify

8-6
rest 8-5.sav
table 1 name load_settlement
def add_load
p_gp = gp_near(2,1,3)
loop n (1,5)
app_load = n * (-1000e3)
file_name = '7-6_add_step' + string(n) + '.sav'
command
app nstress app_load ran z 2.9 3.1 x 1 2 y 1 2
solve
save file_name
endcommand
xtable(1,n) = -1.0 * app_load
ytable(1,n) = gp_zdisp(p_gp)
endloop
end
add_load
save 8-6.sav

8-7
rest 8-6.sav
def find_max_disp
p_gp = gp_head
maxdisp_value = 0.0
maxdisp_gpid = 0
loop while p_gp # null
disp_gp = sqrt(gp_xdisp(p_gp) ^ 2 + gp_ydisp(p_gp) ^ 2 + gp_zdisp(p_gp) ^ 2)
if disp_gp > maxdisp_value
maxdisp_value = disp_gp
maxdisp_gpid = gp_id(p_gp)
endif
p_gp = gp_next(p_gp)
endloop
end
find_max_disp
print maxdisp_value maxdisp_gpid

rest 8-6.sav
config zextra 1
def get_sigma_dif
p_z = zone_head
loop while p_z # null
sigma_dif = z_sig3(p_z) - z_sig1(p_z)
z_extra(p_z,1) = sigma_dif
p_z = z_next(p_z)
endloop
end
get_sigma_dif
plot con zextra 1

9-1
;---------------------------------------------------
;
; 移来移去法接触面的建立
;---------------------------------------------------
n
gen zone radcyl 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.15
gen zone radcyl 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 fill
gen zone reflect dd 270 dip 90
group clay
;
interface 1 face range cylinder end1 (0,0,0) end2 (0,0,-5.1) radius .31 &
cylinder end1 (0,0,0) end2 (0,0,-5.1) radius .29 not
interface 2 face range cylinder end1 (0,0,-4.9) end2 (0,0,-5.1) radius .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 6
gen 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 6
gen zone reflect dd 270 dip 90 range z 1 6.1
group pile range z 1 6.1
pause
ini z add -6.0 range group pile
save pile_geom.sav

9-2
;导来导去法
;---------------------------------------------
n
gen zone brick size 3 3 3
group 2 range x 1 2 y 1 2 z 1 2
group 1 range gr 2 not
save 1.sav
del ran group 2 not

interface 1 face range x 1 y 1 2 z 1 2
interface 1 face range x 2 y 1 2 z 1 2
interface 1 face range x 1 2 y 1 z 1 2
interface 1 face range x 1 2 y 1 2 z 1
interface 1 face range x 1 2 y 1 2 z 2

rest 1.sav
del ran group 2
expgrid 1.fac3d

impgrid 1.flac3d

model ela
fix x y z ran z 0
ini den 2000
set grav 0 0 -10
interface 1 prop kn 20e6 ks 20e6 coh 10e3 fri 15
app nstr -200e3 ran x 0 1 y 1 2 z 3
solve

9-3
;切割模型法
;------------------------------
n
gen 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 separate 1
int 1 wrap 1 2
int 1 maxedge 0.5
plo int red

9-4
; Create Material 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.0
group Material
; Create 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 zon