求助：地层中空腔内流固耦合程序

scl5945

建立盐岩地层中一腔体，为圆柱形，如果里边充满流体，计算渗流、同时考虑蠕变效应，这模型该怎么建立呢，荷载如何施加？下边是几何参数部分写的程序。new
title 开挖试算
config creep
def kkk
rr=10
hh=35
rrt1=rr-0.01
rrt2=rr+0.01
rrf=rr+1.75
hhf=hh+1.75
hht1=hh-0.02
hht2=hh+0.02
rrr=50
uhh=hh+30
dhh=-30
lxbrrr=-rrr-0.01;x向模型左边界左
lrxbrrr=-rrr+0.01;x向模型右边界右
rlxbrrr=rrr-0.01
rxbrrr=rrr+0.01
dzbhh=dhh-0.01
uzbhh=dhh+0.01
upzb=uhh-1.75
pmh_0=0
pmh_1=hh/5
pmh_2=2*pmh_1
pmh_3=3*pmh_1
end
kkk
;圆柱体外环绕放射状网格溶腔简化模型
generate zone radcylinder p0 0 0 hh p1 50 0 hh p2 0 0 0 p3 0 50 hh  size 5 5 5 10  rat 1 1 1 1.2 &
dim 10 10 10 10 ;fill group Maincave;
generate zone radcylinder p0 0 0 uhh p1 50 0 uhh p2 0 0 hh p3 0 rrr uhh  size 5 5 5 10 &
rat 1 1 1 1.2 dim 10 10 10 10 fill group Upcov
generate zone radcylinder p0 0 0 0 p1 rrr 0 0 p2 0 0 dhh p3 0 rrr 0  size 5 5 5 10 rat 1 1 1 1.2 &
dim 10 10 10 10 fill group Downcov
gen zon reflect dip 90 dd 90 orig 0 0 0 ;以圆点为准、竖向为对称轴进行对称操作
gen zon reflect dip 90 dd 0 ori 0 0 0 ;以圆点为准、水平方向为对称轴进行对称操作
interface 1 face range cyl end1 0 0 0.1 end2 0 0 hh r rr
interface 2 face range cyl end1 0 0 0.1 end2 0 0 -0.1 r rr
interface 3 face range cyl end1 0 0 hht1 end2 0 0 hht2 r rr
plo inter yel ske

scl5945

1# scl5945 另外，接触面加的是否妥当也请大家帮忙看看好么