善良的喵 发表于 2016-12-8 08:54:33

渗流情况下为什么最后收敛那里不动了

new

config fluid
set fluid off
gen zone bri p0 0 0 0 p1 3 0 0 p2 0 3 0 p3 0 0 3 p4 3 3 0 p5 0 3 3 p6 3 0 3 p7 3 3 3 size 3 3 3

group 1 range z 0 1
group 2 range z 1 3

plot set plane ori 1 1 1 norm 0 1 0

model elas
prop density 1700 bulk 1e8 shear 4e7 range group 1
prop density 1700 bulk 1e8 shear 4e7 range group 2

model fl_iso
ini fdens 1000

ini ftens -5e5


prop por 0.4 perm 10e-17 range group 1
prop por 0.4 perm 10e-17 range group 2



fix x y z range z -0.1 0.1
fix y range y -0.1 0.1
fix y range y 2.9 3.1
fix x range x -0.1 0.1
fix x range x 2.9 3.1

set grav 0 0 -10

water dens 1000
water table ori 0 0 1300 norm 0 0 1


plocon pp plane

set mech ratio 1e-4
solve

save 初始应力场.sav

ini xdis=0 ydis=0 zdis=0
ini xvel=0 yvel=0 zvel=0

model mohr range group1 any group2 any

prop bulk 1e8 shear 4e7 fric 10 coh 6e5 tens 3e5 range group1 any group2 any


app nstress -3.2e4 range z 2 3 x 1 2 y 1 2

set mech ratio 1e-4
solve
save 未加载

这个程序有什么错误吗,模拟渗流模式
页: [1]
查看完整版本: 渗流情况下为什么最后收敛那里不动了