模型400×300×180，陷落柱为正方体高140，,长30宽30，陷落柱从底部发育。模型顶部加载3兆帕的压力，底部是承压含水层，水压2兆帕。怎么渗流耦合分析?】
config fluid
gen zon brick p0 0 0 0 p1 400 0 0 p2 0 300 0 p3 0 0 18 size 40 30 18 group 石灰岩
gen zon brick p0 0 0 18 p1 400 0 18 p2 0 300 18 p3 0 0 29 size 40 30 11 group 砂岩
gen zon brick p0 0 0 29 p1 400 0 29 p2 0 300 29 p3 0 0 66 size 40 30 37 group 泥岩
gen zon brick p0 0 0 66 p1 400 0 66 p2 0 300 66 p3 0 0 88 size 40 30 22 group 砂岩
gen zon brick p0 0 0 88 p1 400 0 88 p2 0 300 88 p3 0 0 92 size 40 30 4 group 泥岩
gen zon brick p0 0 0 92 p1 400 0 92 p2 0 300 92 p3 0 0 98 size 40 30 6 group 煤层
gen zon brick p0 0 0 98 p1 400 0 98 p2 0 300 98 p3 0 0 105 size 40 30 7 group 泥岩
gen zon brick p0 0 0 105 p1 400 0 105 p2 0 300 105 p3 0 0 141 size 40 30 36 group 砂岩
gen zon brick p0 0 0 141 p1 400 0 141 p2 0 300 141 p3 0 0 159 size 40 30 18 group 泥岩
gen zon brick p0 0 0 159 p1 400 0 159 p2 0 300 159 p3 0 0 180 size 40 30 21 group 砂岩
group 陷落柱 ran x 180 220 y 130 170 z 0 159
m e
prop bu 3e10 sh 1.8e10 ran group 砂岩 ;co 3.8e6 fric 40 ten 3.1e6
prop bu 2.7e10 sh 1.5e10 ran group 泥岩 ;co 3.98e6 fric 24 ten 2.4e6
prop bu 7e9 sh 2e9 ran group 煤层 ;co 0.7e6 fric 27 ten 1.5e6
prop bu 3e10 sh 2e10 ran group 石灰岩 ;co 2.6e6 fric 35 ten 3.2e6
prop bu 1.8e9 sh 0.94e9 ran group 陷落柱 ;co 0.8e6 fric 28 ten 0.28e6
model fl_iso
prop perm 1.5e-15 poro 0.2 ran group 砂岩
prop perm 2.55e-15 poro 0.18 ran group 泥岩
prop perm 1.02e-15 poro 0.26 ran group 煤层
prop perm 3.672e-15 poro 0.25 ran group 石灰岩
prop perm 4.386e-14 poro 0.5 ran group 陷落柱
ini fden 1000 fmod 2e9 ften -1e10
ini pp 3e6 grad 0 0 -1.7e4 ran z 0 ;承压水压？怎么设置？
fix pp 0 ran z 0 ;底部透水边界？怎么设置？
ini density 2700 ran group 砂岩
ini density 2300 ran group 泥岩
ini density 1400 ran group 煤层
ini density 2000 ran group 陷落柱
ini density 2600 ran group 石灰岩
fix x ran x -0.1 0.1
fix x ran x -399.9 400.1
fix y ran y -0.1 0.1
fix y ran y 299.9 300.1
fix z ran z 179.9.1 180.1
app szz -3.2e6 ran z -0.1 180.1
app sxx -3.84e6 ran x -.1 .1
app sxx -3.84e6 ran x 399.9 400.1
app syy -3.84e6 ran y -.1 .1
app syy -3.84e6 ran y 299.9 300.1
apply nsterss -3.25e6 ran z 179.9 180.1 ;施加顶部载荷
set grav 0 0 -10
solve
save 39.sav
想要得到的结果如下图、
，
内容来自[短消息]
