请教高手—彭文斌书上的边坡实例
;---边坡稳定性分析---new
;-------------------------------
;露天采矿场边坡曲率稳定性分析
;--------------------------------
;创建采矿场（组合矩形网格和环绕柱状网格）
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 68 0 25 p7 add 68 40 25 size 10 12 10 rat 1.1 1 1
;创建地下水位
ini pp 0 grad 0 0 -1e4 range z -15.1 0
def parm
nptabl=7
nprof=10
end
parm
table 1 12 0 18.25 12.5 25 17 35 21 60 24 75 25 150 25
def aux1
loop n(1,nptabl)
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+15,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,nptabl-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 xx2,yy2,zz2 xx4,yy4,zz4
end_loop
end_loop
end
aux2
ret
;设置 MOHR COULOMB 模型及材料参数
mod mo
pro bulk 2e8 she 1e8 fric 45 coh 1e5 ten 1e5
;边界条件
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
;初始化密度；地下水位之上为2500，之下为2600
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
;初始化重应力
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
;监测
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.525
hist gp zdisp 0 24.525
;自动安全因子解算
solve fos
ret
貌似执行到water.fis就有错误，我把water.fis直接合并在一个文本里，但是执行后出现surfaceparallel tocurrentgravity vector请教高手指导！怎么正确运行！ 在线等高手！ 没人理啊？ 看来这世道真没人理了！！！！ :):):):):):)
页:
[1]