新手求解，输出显示错误，请大神指点！！！
config fluidgen zone brick p0 0 0 -25 p1 12.4 0 -25 p2 0 3 -25 p3 0 0 -20 size 8 3 4 group zsand6
plot add surface red
plot add axes yellow
plot set magnification 0.1
plot show
gen zone brick p0 12.4 0 -25 p1 40 0 -25 p2 12.4 3 -25 p3 12.4 0 -20 ratio 1.1 1 1 size 12 3 4 group zsand
gen zone brick p0 0 0 -20 p1 12.4 0 -20 p2 0 3 -20 p3 0 0 -0.5 ratio 1 1 0.98 size 8 3 20 group clay
gen zone brick p0 12.4 0 -20 p1 40 0 -20 p2 12.4 3 -20 p3 12.4 0 -0.5 ratio 1.1 1 0.98 size 12 3 20 group clay
gen zone brick p0 0 0 -0.5 p1 12.4 0 -0.5 p2 0 3 -0.5 p3 0 0 0 size 8 3 1 group stone
gen zone brick p0 12.4 0 -0.5 p1 40 0 -0.5 p2 12.4 3 -0.5 p3 12.4 0 0 ratio 1.1 1 0.95 size 12 3 1 group stone
gen zone brick p0 0 0 0 p1 12.4 0 0 p2 0 3 0 p3 0 0 3.3 p4 12.4 3 0 &
p5 0 3 3.3 p6 7.45 0 3.3 p7 7.45 3 3.3 size 8 3 10 group huitian
gen zone brick p0 0 0 3.3 p1 7.45 0 3.3 p2 0 3 3.3 p3 0 0 4.1 p4 7.45 3 3.3 &
p5 0 3 4.1 p6 6.25 0 4.1 p7 6.25 3 4.1 size 7 3 3 group huitian2
gen zone reflect n (-1 0 0) o = (0 0 0)
plot block group
fix x ran x -39.9 -40.1
fix x ran x 39.9 40.1
fix y ran y -0.1 0.1
fix y ran y 2.9 3.1
model mohr ran z -25 0
model null ran z 0 4.1
prop bulk 17.5e6 shear 13.125e6 coh 1e10 tension 1e10 ran group zsand
ini dens 2100 ran group zsand
prop bulk 7.7e6 shear 1.6e6 coh 1e10 tension 1e10 ran group clay
ini dens 1950 ran group clay
prop bulk 11e6 shear 8.25e6 coh 1e10 tension 1e10 ran group stone
ini dens 2300 ran group stone
model fl_iso
ini pp 0 grad 0 0 -10e3 ran z -6 -25
ini ftens -1e-3
ini fdens 1000
ini fmod 2e8
prop perm 1e-11 poro 0.45 ran z -6 -25
fix pp 0 range z -5.9 -6.1
fix pp range z -25.1 -24.9
set fluid off
set grav 0 0 -9.8
plot con zdis sk
solve
prop bulk 17.5e6 shear 13.125e6 coh 0 fric 36 ran group zsand
prop bulk 7.7e6 shear 1.6e6 coh 10e3 fric 8 ran group clay
prop bulk 11e6 shear 8.25e6 coh 0 fric 30 ran group stone
solve
ini xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0
;土体中心创建桩
sel pile id 1 beg=(0.75 0.2 0.0) end (0.75 0.2 -18) nseg 18
sel pile id 1 beg=(2.25 0.2 0.0) end (2.25 0.2 -18) nseg 18
sel pile id 1 beg=(3.75 0.2 0.0) end (3.75 0.2 -18) nseg 18
sel pile id 1 beg=(5.25 0.2 0.0) end (5.25 0.2 -18) nseg 18
sel pile id 1 beg=(6.75 0.2 0.0) end (6.75 0.2 -18) nseg 18
sel pile id 1 beg=(8.25 0.2 0.0) end (8.25 0.2 -18) nseg 18
sel pile id 1 beg=(9.75 0.2 0.0) end (9.75 0.2 -18) nseg 18
sel pile id 1 beg=(11.25 0.2 0.0) end (11.25 0.2 -18) nseg 18
sel pile id 1 beg=(12.75 0.2 0.0) end (12.75 0.2 -18) nseg 18
sel pile id 1 beg=(0.0 1.5 0.0) end (0.0 1.5 -18) nseg 18
sel pile id 1 beg=(1.5 1.5 0.0) end (1.5 1.5 -18) nseg 18
sel pile id 1 beg=(3.0 1.5 0.0) end (3.0 1.5 -18) nseg 18
sel pile id 1 beg=(4.5 1.5 0.0) end (4.5 1.5 -18) nseg 18
sel pile id 1 beg=(6.0 1.5 0.0) end (6.0 1.5 -18) nseg 18
sel pile id 1 beg=(7.5 1.5 0.0) end (7.5 1.5 -18) nseg 18
sel pile id 1 beg=(9.0 1.5 0.0) end (9.0 1.5 -18) nseg 18
sel pile id 1 beg=(10.5 1.5 0.0) end (10.5 1.5 -18) nseg 18
sel pile id 1 beg=(12 1.5 0.0) end (12 1.5 -18) nseg 18
sel pile id 1 beg=(0.75 2.8 0.0) end (0.75 2.8 -18) nseg 18
sel pile id 1 beg=(2.25 2.8 0.0) end (2.25 2.8 -18) nseg 18
sel pile id 1 beg=(3.75 2.8 0.0) end (3.75 2.8 -18) nseg 18
sel pile id 1 beg=(5.25 2.8 0.0) end (5.25 2.8 -18) nseg 18
sel pile id 1 beg=(6.75 2.8 0.0) end (6.75 2.8 -18) nseg 18
sel pile id 1 beg=(8.25 2.8 0.0) end (8.25 2.8 -18) nseg 18
sel pile id 1 beg=(9.75 2.8 0.0) end (9.75 2.8 -18) nseg 18
sel pile id 1 beg=(11.25 2.8 0.0) end (11.25 2.8 -18) nseg 18
sel pile id 1 beg=(12.75 2.8 0.0) end (12.75 2.8 -18) nseg 18
sel pile id 1 beg=(-0.75 0.2 0.0) end (-0.75 0.2 -18) nseg 18
sel pile id 1 beg=(-2.25 0.2 0.0) end (-2.25 0.2 -18) nseg 18
sel pile id 1 beg=(-3.75 0.2 0.0) end (-3.75 0.2 -18) nseg 18
sel pile id 1 beg=(-5.25 0.2 0.0) end (-5.25 0.2 -18) nseg 18
sel pile id 1 beg=(-6.75 0.2 0.0) end (-6.75 0.2 -18) nseg 18
sel pile id 1 beg=(-8.25 0.2 0.0) end (-8.25 0.2 -18) nseg 18
sel pile id 1 beg=(-9.75 0.2 0.0) end (-9.75 0.2 -18) nseg 18
sel pile id 1 beg=(-11.25 0.2 0.0) end (-11.25 0.2 -18) nseg 18
sel pile id 1 beg=(-12.75 0.2 0.0) end (-12.75 0.2 -18) nseg 18
sel pile id 1 beg=(-1.5 1.5 0.0) end (-1.5 1.5 -18) nseg 18
sel pile id 1 beg=(-3.0 1.5 0.0) end (-3.0 1.5 -18) nseg 18
sel pile id 1 beg=(-4.5 1.5 0.0) end (-4.5 1.5 -18) nseg 18
sel pile id 1 beg=(-6.0 1.5 0.0) end (-6.0 1.5 -18) nseg 18
sel pile id 1 beg=(-7.5 1.5 0.0) end (-7.5 1.5 -18) nseg 18
sel pile id 1 beg=(-9.0 1.5 0.0) end (-9.0 1.5 -18) nseg 18
sel pile id 1 beg=(-10.5 1.5 0.0) end (-10.5 1.5 -18) nseg 18
sel pile id 1 beg=(-12 1.5 0.0) end (-12 1.5 -18) nseg 18
sel pile id 1 beg=(-0.75 2.8 0.0) end (-0.75 2.8 -18) nseg 18
sel pile id 1 beg=(-2.25 2.8 0.0) end (-2.25 2.8 -18) nseg 18
sel pile id 1 beg=(-3.75 2.8 0.0) end (-3.75 2.8 -18) nseg 18
sel pile id 1 beg=(-5.25 2.8 0.0) end (-5.25 2.8 -18) nseg 18
sel pile id 1 beg=(-6.75 2.8 0.0) end (-6.75 2.8 -18) nseg 18
sel pile id 1 beg=(-8.25 2.8 0.0) end (-8.25 2.8 -18) nseg 18
sel pile id 1 beg=(-9.75 2.8 0.0) end (-9.75 2.8 -18) nseg 18
sel pile id 1 beg=(-11.25 2.8 0.0) end (-11.25 2.8 -18) nseg 18
sel pile id 1 beg=(-12.75 2.8 0.0) end (-12.75 2.8 -18) nseg 18
sel pile id=1 prop emod=8e10 nu=0.3 xcarea=1.256e-3 &
xcj=1.57e-4 xciy=7.85e-5 xciz=7.85e-5 &
per=1.256 dens 2500 &
cs_sk=1.3e11 cs_scoh=1e4 cs_sfric=12 &
cs_nk=1.3e11 cs_ncoh=0 cs_nfric=0 &
cs_ngap=off
sel geogrid id 1range x = (-11.9 11.9 ) z = (-0.01 0.01)
sel node init zpos add -0.3 range z = (-0.01 0.01)
sel geogrid id 2range x = (-11.9 11.9 ) z = (-0.01 0.01)
sel node init zpos add -0.15 range z = (-0.01 0.01)
sel geogrid prop iso = (26e9 0.33) thick 5e-3 &
cs_sk = 2.3e6 cs_scoh = 0.0
sel geogrid prop cs_sfric 30
set fluid on
set mech force 1.5e3
set mech subs 10000 auto
set fluid subs 30
def day
day=fltime/24/3600
day
hist id=1 day
model elastic ran z 0 0.8 ;激活0 m ~ 0.8 m的单元
prop bulk 5e6 shear 7.5e6 ran z 0 0.81 j
ini dens 1900 ran z 0 0.8
hist id=2 gp zdis 0 0 0 ;记录底部中心点的沉降
hist id=3 gp zdis 0 0 0.8 ;记录顶部中心点的沉降
hist id=4 gp zdis 12.4 0 0 ;记录坡脚处的沉降
hist id=5 gp xdis 12.4 0 0;记录坡脚处的水平位移
hist id=6 gp zdis 11.2 0 0.8 ;记录路肩处的沉降
hist id=7 gp xdis 11.2 0 0.8 ;记录路肩处的水平位移
solve age 7.776e6 ;按软件默认精度求解
model elastic ran z 0.8 1.6 ;激活0.8m ~ 1.6m的单元
prop bulk 7e6 shear 7.5e6 ran z 0.8 1.67
ini dens 1900 ran z 0.8 1.65
hist id=8 gp zdis 0 0 0 ;记录底部中心点的沉降
hist id=9 gp zdis 0 0 1.6 ;记录顶部中心点的沉降
hist id=10 gp zdis 12.4 0 0 ;记录坡脚处的沉降
hist id=11 gp xdis 12.4 0 0;记录坡脚处的水平位移
hist id=12 gp zdis 10 0 1.6 ;记录路肩处的沉降
hist id=13 gp xdis 10 0 1.6;记录路肩处的水平位移
solve age 7.776e6 ;按软件默认精度求解
model elastic ran z 1.6 2.4 ;激活1.6m ~ 2.4m的单元
prop bulk 7e6 shear 7.5e6 ran z 1.6 2.4
ini dens 1900 ran z 1.6 2.4
hist id=14 gp zdis 0 0 0 ;记录底部中心点的沉降
hist id=15 gp zdis 0 0 3 ;记录顶部中心点的沉降
hist id=16 gp zdis 12.4 0 0 ;记录坡脚处的沉降
hist id=17 gp xdis 12.4 0 0;记录坡脚处的水平位移
hist id=18 gp zdis 8.8 0 2.4 ;记录路肩处的沉降
hist id=19 gp zdis 8.8 0 2.4;记录路肩处的水平位移
solve age 7.776e6 ;按软件默认精度求解
model elastic ran z 2.4 3.3 ;激活2.4m ~ 3.3m的单元
prop bulk 7e6 shear 7.5e6 ran z 2.4 3.3
ini dens 1900 ran z 2.4 3.31
hist id=20 gp zdis 0 0 0 ;记录底部中心点的沉降
hist id=21 gp zdis 0 0 3.3 ;记录顶部中心点的沉降
hist id=22 gp zdis 12.4 0 0 ;记录坡脚处的沉降
hist id=23 gp xdis 12.4 0 0;记录坡脚处的水平位移
hist id=24 gp zdis 7.45 0 3.3 ;记录路肩处的沉降
hist id=25 gp zdis 7.45 0 3.3;记录路肩处的水平位移
solve age 7.776e6 ;按软件默认精度求解
model elastic ran z 3.3 3.95 ;激活3.3m ~ 3.95m的单元
prop bulk 7e6 shear 7.5e6 ran z 3.3 3.95
ini dens 1900 ran z 3.3 3.95
hist id=26 gp zdis 0 0 0 ;记录底部中心点的沉降
hist id=27 gp zdis 0 0 3.95 ;记录顶部中心点的沉降
hist id=28 gp zdis 12.4 0 0 ;记录坡脚处的沉降
hist id=29 gp xdis 12.4 0 0;记录坡脚处的水平位移
hist id=30 gp zdis 6.475 0 3.95 ;记录路肩处的沉降
hist id=31 gp zdis 6.475 0 3.95 ;记录路肩处的水平位移
solve age 7.776e6 ;按软件默认精度求解
也是做流固耦合的，同病相怜 好像问题一下就看出来了, end?
页:
[1]