;问题描述：本问题为模拟土体上部加载，土体为软土，在荷载作用下，发生大变形，土体侧向限制变形，只允许竖向变形，土体模型为摩尔库伦模型。拟采用流固耦合模型解决问题，求得土体中各点的孔压、应力及变形量。
new
title
hydrulic filling model
config fluid
; --- mode geometry ---
gen zone brick size 20 1 10
; --- mechnical mode ---
model mohr
set fluid biot on
prop bulk=5e8, shear=1e8, friction=35, cohesion=2e10 &
tension=1e10
fix x y range x -.1 .1 x 19.9 20.1
fix x y z range z -.1 .1
fix y
ini den 2000 ;土颗粒密度2000，手册中很多例题都没有考虑土体密度，是不是因为忽略了重力的固结因素？
apply szz 2.e6 range z 19.9 20.1 ; 土体上部加载
; --- fluid flow model ---
model fl_iso
prop perm 2.3e-6 biot_c 1. ;渗透系数应该是指水在土体中的渗透系数吧？为什么手册例题上的系数取值都那么小呢？
ini biot_mod 4.e9 ;如果没有这一行命令，计算很快，但是没有结果，有了这一行之后，计算又不收敛，这一行命令到底是做什么的呢？
ini pp 0
ini den 1000 ;水密度1000
; --- setting ---
set gravity 0 0 -9.8
set large ;大变形计算
; --- fish functions ---
def point
pnt = gp_near(10.,0.,5.)
zpnt = gp_znear(10.,0.,5.)
end
point
def pp10
pp10 = gp_pp(pnt)
ft = fltime
c_szz = z_szz(zpnt)
end
; --- drained response ---
fix pp 0 range z 19.9 20.1
set fluid on
set mech force 0 ratio 1e-4
set hist_rep 200
set fluid substep 1
set mech substep 1 auto ;slave
solve age 500
; --- histories ---
hist pp10
hist c_szz
hist ft
; --- plot ---
plot create pp
plot set title text
Vertical displacement versus time at column mid-height
plot add his 1 vs 3
plot add his 2 vs 3
plot show
;运行不收敛！哪里出了问题呢？
