边坡计算实例

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 de

n 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
endcommand
end loop
end loop
end
aux2
ret