config fluid 渗流模式打开
set fluid biot on 比奥计算开启，说明考虑了土颗粒的可压缩性
gen zone bri p0 0 0 0 p1 5 0 0 p2 0 2 0 p3 0 0 20 size 3 1 10 模型网格划分
group Face range z 18 20 分组砂层
group Soft range z 10 18 ;soft soil 8m 分组软土
group Clay range z 0 10 分组粘土
model mohr 定义摩尔库伦模型
prop bulk 33e6 shear 7e6 c 0 f 30 range group Face 设置三种材料参数
prop bulk .83e6 shear .17e6 c 7.0e3 f 5 range group Soft
prop bulk 6.66e6 shear 1.42e6 c 8.0e3 f 20 rang group Clay
ini dens 1500 设置土体密度 这里土体密度为什么设置了一个，明明三种嘛？？？
model fl_iso 定义各向同性
ini ftens -5e5 流体抗压强度，抗拉强度怎么得来的？
ini fdens 1000 流体密度 也就是水的密度嘛？
ini biot_mod 4e9
prop perm 10e-13 biot_c 1 定义渗透系数 三种材料同性，所以定义一个，那么具体是哪一个，是砂，软土，粘土？，定义比奥系数，前面set fluid biot on 比奥计算开启，说明考虑了土颗粒的可压缩性，那比奥系数怎么为1？
prop perm 10e-10 ran grou Face
fix x y z range z -0.1 0.1
fix y rang y -.01 0.01
fix y rang y 1.9 2.1

fix x range x -.01 .01
fix x range x 4.9 5.1
;Initial Stress K0=1 这里的ko是说泊松比都定义为0.5 的意思嘛？
set grav 0 0 -10
ini szz -400e3 grad 0 0 20e3
ini syy -400e3 grad 0 0 20e3
ini sxx -400e3 grad 0 0 20e3
ini pp 200e3 grad 0 0 -10e3 这些初始应力和梯度是怎么得来的
set fluid off 渗流计算关闭
step 5000
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
;for fixed pp on the PVD position
def fixPP
loop i(1,10)
pv = -80e3 +(i-1)*10e3
zt = 20.1-(i-1)
zb = 20.1-i
command
fix pp pv rang x -.1 .1 z zb zt
end_command
end_loop
end fixpp
step 5000 （这一部分是真空度的传递规律嘛）
set fluid on
ini biot_mod 4.3e7 rang group Soft not ;modify biot_mod
ini biot_mod 9.0e6 rang group Soft
set mech force 1.5e3
set mech subs 100 auto ;slave
set fluid subs 30
; --- histories ---
def day day=fltime/24/3600
end
day
hist id=1 day
hist id=2 gp zd 2.5,0,20
hist id=3 gp zd 2.5,0,10
sol age 7.776e6

欢迎大神们来讲解

这只是一个例子，为了演示执行的过程，里面的参数是次要的，不是按照工程实例给出的，不必纠结

三种土设置一个密度这不就是为了简化吗？比奥系数为1就是不可压缩，为什么计算考虑可压缩性又设置系数为1这不明摆着只是演示一个执行过程，至于流体的抗拉强度，在书中其他例子里你可以看看流体的抗拉强度给了好几个值，岩土工程中一般遇到的都是水为什么值不一样，比奥系数为1还有一个可能就是方便计算比奥模量，因为只有比奥系数为1时才有公式可以准确计算比奥模量

@如果没有你luc

prop perm 10e-13 biot_1这里的渗透系数是指那个土体的？并没有指出来啊？

请问通常粘聚力和内摩察角是通过什么试验得来的？

@如果没有你luc

如何通过打印操作把沉降数值拷贝出来

三种土体设置一个密度，的确是为了简化，这里的密度是指干密度

这里也没有考虑地下水位的变化