一下是我自己写的一段小程序,大致就是一个梯形体,高50m,左侧30m处是河水面,假设地下水和水面平齐,一个很简单的模型,但是初始渗流没有一点收敛的意思,另外,如果把左边的边界改成垂直的,就立马可以收敛,实在是困扰我很久了,向大家求助!
new
config fluid
gen zo br p0 0 0 0 p1 50 0 0 p2 0 1 0 p3 20 0 50 p4 50 1 0 p5 20 1 50 p6 50 0 50 p7 50 1 50 size 50 1 50
m m
prop bu 1e9 sh 1e9 c 1e6 fr 42 density 2000 ten 1e10
fix x y z range z -0.1 0.1
fix x range x -0.1 0.1
fix x range x 49.9 50.1
fix y
m fl_iso
prop perm 1e-12 poro 0.3
ini fden 1e3
ini ften -1e10
ini fmod 2e9
set grav 10.0
ini pp 30e4 grad 0 0 -1e4 range z 0 30
ini sat 0
ini sat 1 range z 0 30
fix pp range plane norm -5 0 2 dist 0.01 ori 0 0 0
fix pp range x 50
set fl on mech of
so