含水层排水产生固结压缩的讨论

colorfish

用命令采用两种模式可以定义流体的压缩性：（1）定义比奥系数 ，比奥模量；（2）定义流体模量和孔隙率。第一种情况考虑到土颗粒的压缩性，当比奥系数等于1时，颗粒不具有压缩性，第二种情况假定土颗粒不可压缩
    那么用prop biot_c<1，ini biot mod 模拟产生压缩，这时候是不是就不需要定义流体模量和空隙率？
    欢迎于各位交流

hh199422

Iin the second case, Fluid modulus and porosity must be given. it is right!

colorfish

那我想定义一个有含水层和隔水层的模型，怎么编程，另外底部含水层的空隙水压力怎么确定，我编的源程序如下，各位仁兄帮我看看对不对，希望提出更好的建议：
new
set logfile renlian
set log on
config fluid
gen zone brick &
                p0 0,0,0 p1 12,0,0 p2 0,0,-4 p3 0,12,0 &
                size 6 2 6  
group z1 range z -4 0
;
gen zone brick &
                p0 0,0,-4 p1 12,0,-4 p2 0,0,-6 p3 0,12,-4 &
                size 6 1 6  
group z2 range z -6 -4
;
gen zone brick &
                p0 0,0,-6 p1 12,0,-6 p2 0,0,-10 p3 0,12,-6 &
                size 6 2 6  
group z3 range z -10 -6
;
gen zone brick &
                p0 0,0,-10 p1 12,0,-10 p2 0,0,-12 p3 0,12,-10 &
                size 6 1 6  
group z4 range z -12 -10
;
plot surf
model mohr range group z1
prop  bulk 6.173e5 shear 4.065e5  coh 0  fric 30   range group z1
;
model mohr range group z2   ;隔水层
prop  bulk 6.667e6 shear 3.077e6 coh 18000 fric 23.0 range group z2
;
model mohr range group z3  ;含水层
prop bulk 32.051e6 shear 20.161e6 coh 200 fric 25 range group z3
;
model drucker range group z4   ;基岩
prop  bulk 1.344e10 kshear 3.0e7 qdil 1.0 qvol 0.7 shear 1.050e10 tension 4.25e8 &
range group z4
;下面的定义对不对
model fl_iso range group z3
prop perm 1.0e-8  poro 0.3
model fl_iso range group z4
prop perm 1.0e-15  poro 0.1
;
;下面的语句是不是要有个范围才对
ini fmod 2.e9   
ini fdensity 1000  

ini pp  2.0e5  grad 0 0 -10000 range group z3     ;对吗
;fix pp         range z -32.1 -31.9
;fix pp         range z -38.1 -37.9
;
ini dens 1030 range group z1
ini dens 980  range group z2
ini dens 910  range group z3
ini dens 2650 range group z4
;
;--------------------------------边界条件
fix z range z -12.1 -11.9
fix x range x -0.1 0.1
fix x range x 11.9 12.1
fix y range y -0.1 0.1
fix y range y 11.9 12.1
;
set grav 0 0 -10
;一含
ini szz 0 grad 0 0 10300 range z -4 0
ini syy 0 grad 0 0 3090 range z -4 0
ini sxx 0 grad 0 0 3090 range z -4 0
;---------------------------------------------------------
;一隔
ini szz -2000 grad 0 0 9800 range z -6 -4
ini syy 4480 grad 0 0 4210 range z -6 -4
ini sxx 4480 grad 0 0 4210 range z -6 -4
;
;底含
ini szz -6200 grad 0 0 9100 range z -10 -6
ini syy -3320 grad 0 0 2910 range z -10 -6
ini sxx -3320 grad 0 0 2910 range z -10 -6
;
;基岩
ini szz 167800 grad 0 0 26500 range z -12 -10
ini syy 31180 grad 0 0 6360 range z -12 -10
ini sxx 31180 grad 0 0 6360 range z -12 -10
;
set fluid  off
;
hist unbal
; Settings
set large
; Cycling
solve rat 1.e-6
save chu.sav
;
plot sketch
plot add cont zdisp
plot show
;
ini xdis 0 ydis 0 zdis 0
ini state 0
;
set mech force 1.5e2
set mech substep 100 slave
;

hist fltime
hist gp zd 2 6 -8
hist gp zd 6 6 -8
hist gp zd 10 6 -8

set fluid on

solve age 0.5e5
plot cont pp
save 2.sav
;
solve age 1e6
plot cont pp
save 3.sav
plot hist 3 4 5 vs 2