你看一下,只在边界的孔压改变了,坡内的还是没有改变,计算时间也很长。
def waterdown
fltime1=fltime-259200
pp_2=7.6e3-fltime1*7.6e3/10800
pp_2_2=pp_2-2e3
p3=-1*pp_2
p4=-1*pp_2_2
line=pp_2/1e4
if line>0.2
command
apply remove nstress
free pp ran z 0 x 3.46 4.02
apply pp pp_1 ran z 0 x 3.46 4.02
apply nstress p1 ran z 0 x 3.46 4.02
free pp ran plane norm 1 0 1.5 ori 0 0 3.46 z 0 0.2
apply pp pp_2 grad 0 0 -1e4 ran plane norm 1 0 1.5 ori 0 0 3.46 z 0 0.2
free pp ran x 3.1 3.16 z 0.1999 0.2001
apply pp pp_2_2 ran x 3.1 3.16 z 0.1999 0.2001
free pp ran plane norm 1 0 2 ori 3.1 0 0.2 z 0.2 0.76
apply pp pp_2 grad 0 0 -1e4 ran plane norm 1 0 2 ori 3.1 0 0.2 z 0.2 line
apply nstress p3 grad 0 0 1e4 ran plane norm 1 0 1.5 ori 0 0 3.46 z 0 0.2
apply nstress p4 ran x 3.1 3.16 z 0.1999 0.2001
apply nstress p3 grad 0 0 1e4 ran plane norm 1 0 2 ori 3.1 0 0.2 z 0.2 line
endcommand
else
command
apply remove nstress
free pp ran x 3.1 3.16 z 0.1999 0.2001
apply pp 0 ran x 3.1 3.16 z 0.1999 0.2001
fix pp ran x 3.1 3.16 z 0.1999 0.2001
free pp ran plane norm 1 0 2 ori 3.1 0 0.2 z 0.2 0.76
apply pp 0 ran plane norm 1 0 2 ori 3.1 0 0.2 z 0.2 0.76
fix pp ran plane norm 1 0 2 ori 3.1 0 0.2 z 0.2 0.76
free pp ran plane norm 1 0 1.5 ori 0 0 3.46 z line 0.2
apply pp 0 ran plane norm 1 0 1.5 ori 0 0 3.46 z line 0.2
fix pp ran plane norm 1 0 1.5 ori 0 0 3.46 z line 0.2
free pp ran plane norm 1 0 1.5 ori 0 0 3.46 z line 0.2
apply pp pp_2 grad 0 0 -1e4 ran plane norm 1 0 1.5 ori 0 0 3.46 z 0 line
fix pp ran plane norm 1 0 1.5 ori 0 0 3.46 z 0 line
apply nstress p3 grad 0 0 1e4 ran plane norm 1 0 1.5 ori 0 0 3.46 z 0 line
free pp ran z 0 x 3.46 4.02
apply pp pp_1 ran z 0 x 3.46 4.02
apply nstress p1 ran z 0 x 3.46 4.02
endcommand
endif
end
set fishcall 10 waterdown