FLAC3D在边坡稳定上的应用

shanxixiaoxin

这是一个边坡实例，与大家共享一下。命令流如下：
;------------------------------------------------------- -----
; sl 3d.dat influence of slope curvature on stability
; of an open pit
;------------------------------------------------------ -----
; create model by joining radial cylinder and brick primitives
gen zone radcyl &
p0 0 0 25 p1 add 80 0 0 p2 add 0 0 -25 p3 add 0 80 0 &
dim 24.5 24.5 12 12 ratio 1 1 1 1.1
gen zone radcyl &
p0 0 0 0 p1 add 80 0 0 p2 add 0 0 -15 p3 add 0 80 0 &
dim 12 12 12 12 ratio 1 1 1 1.1 fill size 5 7 10 10
gen zone brick &
p0 0 -40 -15 p1 add 12 0 0 p2 add 0 40 0 p3 add 0 0 15 &
size 5 12 7
gen zone brick &
p0 12 -40 -15 p1 add 68 0 0 p2 add 0 40 0 p3 add 0 0 15 &
size 10 12 7 ratio 1.1 1 1
gen zone brick &
p0 12 -40 0 p1 add 68 0 0 p2 add 0 40 0 p3 add 12.5 0 25 &
p4 add 68 40 0 p5 add 12.5 40 25 p6 add 68 0 25 p7 add 68 40 25 &
size 10 12 10 rat 1.1 1 1
; install water table
ini pp 0 grad 0 0 -1e4 range z -15.1 0
call water.fis
; assign Mohr Coulomb model and properties
mod mo
pro bulk 2e8 she 1e8 fric 45 coh 1e6 ten 1e6
; boundary conditions
fix x range x -.1 .1
fix x range x 79.9 80.1
fix y range y -40.1 -39.9
fix y range y 79.9 80.1
fix x range z -15.1 -14.9
fix y range z -15.1 -14.9
fix z range z -15.1 -14.9
; initialize density: saturated and unsaturated
def ini dens
pnt = zone head
loop while pnt # null
if z pp(pnt) # 0.0 then
z density(pnt) = 2600
else
z density(pnt) = 2500
endif
pnt = z next(pnt)
endloop
end
ini dens
set g 0 0 -10
; initialize gravity stresses
ini syy -6.25e5 grad 0 0 2.5e4
ini sxx -6.25e5 grad 0 0 2.5e4
ini szz -6.25e5 grad 0 0 2.5e4
; histories
hist unbal
hist gp xdisp 24.5 -40 25
hist gp zdisp 24.5 -40 25
hist gp xdisp 24.5 0 25
hist gp ydisp 24.5 0 25
hist gp zdisp 24.5 0 25
hist gp ydisp 0 24.5 25
hist gp zdisp 0 24.5 25
; equilibrate
solve
save sl eq.sav
; reduce strength
pro coh 3.7e4
ini xd 0
ini yd 0
ini zd 0
ini xv 0
ini yv 0
ini zv 0
step 4000
save sl fail.sav
;
; create plot views of results
;
; view : sl geom (model grid)
plot create sl geom
set out sl geom.plt
plot set rot 20 0 50
plot set ang 20
plot set cent 40 20 5
plot add surface lcyan
plot add axes brown
plot hard
;
; view : sl water (water table with model grid)
plot create sl water
set out sl water.plt
plot set rot 20 0 50
plot set ang 20
plot set cent 40 20 5
plot add water lgra
plot add sketch black
plot hard
;
; view : sl ppres (pore pressure contours)
plot create sl ppres
set out sl ppres.plt
plot set rot 20 0 50
plot set ang 20
plot set cent 40 20 5
plot add cont pp int 2.5e4 sh on
plot hard
;
; view : sl disp (displacement contours)
plot create sl disp
set out sl disp.plt
plot set rot 20 0 50
plot set ang 20
plot set cent 40 20 5
plot add cont disp sh on
plot hard
;
; view : sl hist (displacement histories)
plot create sl hist
set out sl hist.plt
plot add hist 2 3 7 8
plot hard
ret

;Data File "SL PS.DAT"
config extra 5
tit
Plane Strain Simulation
grid 15 17
mod mo
mo null i 1 5 j 8 17
gen 0 -15 0 0 12 0 12 -15 i 1 6 j 1 8
gen s s 80 0 80 -15 rat 1.1 1 i 6 16 j 1 8
gen s 24.5 25 80 25 s rat 1.1 1 i 6 16 j 8 18
fix x i 1
fix x y j 1
fix x i 16
pro bul 2e8 she 1e8 den 2500 fric 45 coh 1e6 ten 1e6
set grav 10
water den 1000
table 1 0 0 12 0 18.25 12.5 25 17 35 21 50 24 74.5 25 80 25
water table 1
def ini dens
loop i (1,izones)
loop j (1,jzones)
if pp(i,j) # 0 then
density(i,j) = 2600.
end if
end loop
end loop
end
ini dens
ini syy -1e6 var 0 1e6
ini sxx -1e6 var 0 1e6
ini szz -1e6 var 0 1e6
his unbal
hist xd i 6 j 18
hist yd i 6 j 18
solve
save sl ps1.sav
ini xd 0 yd 0
ini xv 0 yv 0
pro coh 3.7e4
step 1250
save sl ps2.sav
set out slps1.plt
plot pen h 2 3
set out slps2.plt
plot pen bou ydisp z fill abs
tit
flac model grid
set out slgria.plt
plot pen grid water iwhite
ret

;Data File "SL AXI.DAT"
config extra 5 axi
tit
Axisymmetry Simulation
grid 15 17
mod mo
mo null i 1 5 j 8 17
gen 0 -15 0 0 12 0 12 -15 i 1 6 j 1 8
gen s s 80 0 80 -15 rat 1.1 1 i 6 16 j 1 8
gen s 24.5 25 80 25 s rat 1.1 1 i 6 16 j 8 18
fix x i 1
fix x y j 1
fix x i 16
pro bul 2e8 she 1e8 den 2500 fric 45 coh 1e6 ten 1e6
set grav 10
water den 1000
table 1 0 0 12 0 18.25 12.5 25 17 35 21 50 24 74.5 25 80 25
water table 1
def ini dens
loop i (1,izones)
loop j (1,jzones)
if pp(i,j) # 0 then
density(i,j) = 2600.
end if
end loop
end loop
end
ini dens
ini syy -1e6 var 0 1e6
ini sxx -1e6 var 0 1e6
ini szz -1e6 var 0 1e6
his unbal
hist xd i 6 j 18
hist yd i 6 j 18
solve
save sl ax1.sav
ini xd 0 yd 0
ini xv 0 yv 0
pro coh 3.7e4
step 1250
save sl ax2.sav
set out slax1.plt
plot pen h 2 3 min -1.7
ret

;Data File "WATER.FIS"
def parm
nptab1=7
nprof=10
end
parm
table 1 12 0 18.25 12.5 25 17 35 21 50 24 75 25 150 25
def aux1
loop n(1,nptab1)
rr=xtable(1,n)
zz=ytable(1,n)
loop k(1,nprof)
case of k
alfa=0.5*(pi+.1)*float(k-2)/float(nprof-2)
xtable(n+10,k)=rr*cos(alfa)
ytable(n+10,k)=rr*sin(alfa)
case 1
xtable(n+10,k)=rr
ytable(n+10,k)=-40.0
case 2
xtable(n+10,k)=rr
ytable(n+10,k)=0.0
end case
end loop
end loop
end
aux1
def aux2
command
set grav 0 0 -10
water dens 1000
end command
loop n(1,nptab1-1)
loop k (1,nprof-1)
xx1=xtable(n+10,k)
yy1=ytable(n+10,k)
zz1=ytable(1,n)
xx2=xtable(n+10,k+1)
yy2=ytable(n+10,k+1)
zz2=zz1
xx3=xtable(n+11,k+1)
yy3=ytable(n+11,k+1)
zz3=ytable(1,n+1)
xx4=xtable(n+11,k)
yy4=ytable(n+11,k)
zz4=zz3
command
WATER TABLE face xx1,yy1,zz1 xx2,yy2,zz2 xx3,yy3,zz3
WATER TABLE face xx1,yy1,zz1 xx3,yy3,zz3 xx4,yy4,zz4
end command
end loop
end loop
end
aux2
ret

phonixs

很好，还有什么可以共享交流的，欢迎阿

haca7841

命令中好多都看不懂,有谁能把这个命令流的详细的解释下啊?(建模的部分可以不解释)
haca7841@126.com
期待ＩＮＧ，谢谢大虾！

guogjf

hp79 wrote:
为什么我运行了几次老是出现"error opening file water.fis"啦?

缺少water的fish语言！

图拉

guogjf wrote:
缺少water的fish语言！

这个fish语言是自己编的,还是软件自己的?

islandgirl

这个是手册里的例子，里面有编好的fish语言的

hja2000415

如果我没有搞错得话，这个好像是手册里得例子吧？我也运算过也出现过这样得问题，你可以把water.fis里得内容拷贝到"SL PS.DAT"里然后运行，应该就没有什么错误提示了！

zjl16

honestzyf

上面的兄弟，我现在用flac做边坡分析的，苦于没有实例，你能否提供哪个手册里东东！
qq：115549297
email:honestzyf@sina.com
谢谢了！

yjchen8

4r 这个是手册中的第一例题呀，不过里面的孔隙水压力这一块，我也不怎么能看懂。

mifan

ms 这就是manual里的第一个例子，不过还是要感谢lz的辛勤劳作。

luhuan

学习，顶一下

winneryzj

关键还得看自己学习。。。。。。。有好东西大家一起学习

li_2001

命令中好多都看不懂,有谁能把这个命令流的详细的解释下啊?(建模的部分可以不解释)