一个基坑流固藕合初始化命令流，请大家帮看看提出的问题

52heming

我从z  0到 -40   0到-3无水  -3到-31有水 下面也没水  
思考很久 可就是不对 欢迎大家 指正  下面是一个错图片 请大家看看  错在哪里
new
set log on
config fluid
def hello
   d=1      ;;;;隧道长度
   h1=-12       ;;;;开挖深度
   h22=-28      ;;;;基底土厚度
   h11=(-1)*h1    ;;;;网格划分
   h2=h1+h22
end
hello
gen zone brick p0 0 0 0 p1 8.2 0 0 p2 0 0 -12 p3 0 d 0  size 2 12 d &
rat 1 1 1    group excavation
gen zone brick p0 0 0 -22 p1 8.2 0 -22 p2 0 0 -40 p3 0 d -22  size 2 18 d &
rat 1 1 1    group  soil_1
gen zone brick p0 0 0 -12 p1 8.2 0 -12 p2 0 0 -22 p3 0 d -12  size 2 10 d &
rat 1 1 1    group  base
gen zone brick p0 8.2 0 0 p1 9.2 0 0 p2 8.2 0 -22 p3 8.2 d 0  size 2 22 d &
rat 1 1 1    group  wall
gen zone brick p0 8.2 0 -22 p1 9.2 0 -22 p2 8.2 0 -40 p3 8.2 d -22 size 2 18 d &
rat 1 1 1 group soil_2
gen zone brick p0 9.2 0 0 p1 40 0 0 p2 9.2 0 -40 p3 9.2 d 0  size 10 40 d &
rat 1.05 1 1 group  soil_3
gen zone ref ori 0 0 0 normal -1 0 0
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;力学模型;;;;;;;;;;;;;;;;;;;;;;;;;
model elas;model mohr
set grav 0,0,-10
prop bulk 2.907e6 shear 1.77e6 range z 0 -1  ;;第1层
prop bulk 2.907e6 shear 1.77e6 range z -1 -5 ;;第2层
prop bulk 3.88e6 shear 2.36e6 range z -5 -7  ;;第3层
prop bulk 3.12e6 shear 1.89e6 range z -7 -8 ;;第4层
prop bulk 5.95e6 shear 3.62e6 range z -8 -12 ;;第5层
prop bulk 3.67e6 shear 2.23e6 range z -12 -14 ;;第6层
prop bulk 5.4e6 shear 3.28e6 range z -14 -15 ;;第6_1层
prop bulk 3.67e6 shear 2.23e6 range z -15 -18 ;;第6层
prop bulk 12.46e6 shear 7.57e6 range z -18 -32 ;;第7层
prop bulk 13.85e6 shear 8.415e6 range z -32 -40 ;;第8层
;;;;;;;密度
ini density 1500 range z 0 -1
ini density 1500 range z -1 -3
ini density 2000 range z -5 -3
ini density 2000 range z -5 -7
ini density 2000 range z -7 -8
ini density 2000 range z -8 -12
ini density 2000 range z -12 -14
ini density 2000 range z -14 -15
ini density 2000 range z -15 -18
ini density 2000 range z -31 -18
ini density 1500 range z -31 -40
ini szz -0.454e5 grad 0 0 15000 ran z 0 -3
ini szz -7.85e5 grad 0 0 20000 ran z -3 -31
ini szz -8.9e5 grad 0 0 15000 ran z -31 -40
ini sxx -0.225e5 grad 0 0 7500 ran z 0 -3
ini sxx -4.425e5 grad 0 0 15000 ran z -3 -31
ini sxx -5.1e5 grad 0 0 7500 ran z -31 -40
ini syy -0.225e5 grad 0 0 7500 ran z 0 -3
ini syy -4.425e5 grad 0 0 15000 ran z -3 -31
ini syy -5.1e5 grad 0 0 7500 ran z -31 -40
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;流力模型;;;;;;;;;;;;;;;;;;;;;;;;;
model fl_iso
;model fl_null range group wall
model fl_null range z -31 -40
prop perm 5.918e-10 poro 0.5 range z -3 -18
prop perm 7.089e-9 poro 0.5  range z -18 -31
;prop perm 5.918e-10 poro 0.5 range z -31 -40
ini sat 0.0 range z 0 -3
ini sat 1 range z -3 -40
ini pp 2.8e5 grad 0 0 -1e4 ran z -3 -31
;fix pp 2.85e5 range z -31
ini fdens 1e3
ini fmod 2e9 ftens -1e-3
;;;;;;;;;;;;;;;;;;;;;;;;;边界条件;;;;;;;;;;;;;;;;
fix x ran x -39.9 -40.1 x 39.9 40.1
fix y ran y 0 y d
fix x y z ran z -40
hist unbal
plot con pp ou on
plot add hist 1
set fluid off
solve

52heming

szzzdpp

zhangyao999

我也正好要学习哈

zhangyao999

我怎么算不收敛啊 算了一晚上还是没结果啊