关于一个单向固结例子的问题

lyxl1005

我在论坛上看到一个例子，命令流做到下面这一部分，然后求解，但是为什么最后孔压不能消散到加压前的状态呢，请大家指教。
new
con fl
gen zone brick size 1 1 10 p0 0 0 -10 p1 1 0 -10 p2 0 1 -10 p3 0 0 0
model fl_iso
prop perm 1e-12 poro 0.3
ini fdensity 9.81e2
ini fmod 2e9 ftens -1e-3
set fl off
set mech on
model elas
prop bulk 2512562.814  shear 1152959.262 dens 1335.7
fix x ran x -.1 .1
fix x ran x 0.9 1.1
fix y ran y -.1 .1
fix y ran y 0.9 1.1
fix z ran z -10.1 -9.9
ini szz 0 grad 0 0 15990.3 range z -10 0
ini syy 0 grad 0 0 12604.69 range z -10 0

ini sxx 0 grad 0 0 12604.69 range z -10 0
ini pp 0 grad 0 0 -9810 ran z -10 0
set grav 0 0 -9.81
pl con pp
solve

;耦合求解
set fl on mech on;
model fl_iso
prop perm  1e-12 poro 0.3
ini fdensity 981
ini fmod 2e9 ftens -1e-3
fix pp 0 range z 0.1 -0.1
;fix pp ran z -9.9 -10.1
app nstr -1e4 range z 0 x 0 1 y 0 1
;set mech force 50
;set fluid substep 100
;set mech  substep 100 auto ;slave
solve

lyxl1005

大家还真是冷淡啊，不过我已经解决了，哈哈。因为耦合计算太慢，刚开始由于曼德尔效应，我没有计算完成，所以才会显得不收敛。
        祝大家以后多多回帖。

科技大虾

楼主还真是厉害哦~我想请问你下：
高填方工程自开始施工填筑荷载开始直至最后稳定的过程中，原土地基和填土层始终伴随着固结和流变沉降，请问如何具体模拟填方的填筑加载过程

sdustqh

3# 科技大虾
分层填土堆载的过程，我是这样摸你的，先建好堆载土的网格但是不赋予其性质，到后来再一层一层的激活单元，模拟分层堆载的情况。  我也不知道这样合理不？但是确实发生了变化。。。