陈育民书上例题中单桩静载试验模拟的疑问

粉色人1986

这是陈老师书上第9章中9.6节的例子，单桩静载试验模拟
本人的问题在下面有红色标出了，请大家赐教！
谢谢！




“移来移去”法接触面的建立
gen zone radcyl p0
(0,0,0)
p1
(8,0,0)
p2
(0,0,-5)
p3
(0,8,0)
&

p4
(8,0,-5) p5
(0,8,-5) p6
(8,8,0)
p7
(8,8,-5) &

p8
(.3,0,0)
p9
(0,.3,0)
p10 (.3,0,-5)
p11 (0,.3,-5) &

size 3 10 6 15
ratio 1 1 1 1.15
gen zone radcyl p0
(0,0,-5)
p1
(8,0,-5)
p2
(0,0,-8)
p3
(0,8,-5)
&

p4
(8,0,-8)
p5
(0,8,-8)
p6
(8,8,-5)
p7
(8,8,-8)
&

p8
(.3,0,-5) p9
(0,.3,-5) p10 (.3,0,-8) p11 (0,.3,-8) &


size 3 6 6 15
ratio 1 1 1 1.15 fill
gen zone reflect dd 270 dip 90
group clay
;
interface 1 face range cylinder end1 (0,0,0)
end2 (0,0,-5.1) radius .31 &

cylinder end1 (0,0,0)
end2 (0,0,-5.1) radius .29 not
interface 2 face range cylinder end1 (0,0,-4.9)
end2 (0,0,-5.1) radius .31
;
pause
gen zone cyl p0 (0,0,6) p1 (.3,0,6) p2 (0,0,1)
p3 (0,.3,6)
&

p4 (.3,0,1) p5 (0,.3,1) &

size 3 10 6
gen zone cyl p0 (0,0,6.1) p1 (.3,0,6.1) p2 (0,0,6)
p3 (0,.3,6.1)
&

p4 (.3,0,6) p5 (0,.3,6) &

size 3 1 6
gen zone reflect dd 270 dip 90 range z 1 6.1
group pile range z 1 6.1
pause
ini z add -6.0 range group pile
save pile_geom.sav

;------------------------------------------------------------------

rest pile_geom.sav
为什么桩跟土的体积模量跟剪切模量值一样，密度也一样大？
model mohr
range group clay
prop bulk 8.333e7
shear 3.846e7
coh 30000 fric 0
range group clay

model elas
range group pile
prop bulk 8.333e7
shear 3.846e7
range group pile

interface 1 prop kn 1e8
ks 1e8
fric 20 coh 30000
interface 2 prop kn 1e8
ks 1e8
fric 20 coh 30000
;
ini dens 1230 range group clay
ini dens 1230 range group pile
；最大疑问的地方，这为什么要加个空单元呢？
model null range z -0.1 0.15
;
fix z range z -8.1 -7.9
fix x range x -8.1 -7.9
fix x range x 7.9 8.1
fix y range y -.1 .1
fix y range y 7.9 8.1
set grav 0 0 -10
    下面描述的是设置初始应力的初始地应力场生成，为什么不在z=5m处作为分界呢？假设在z=5.5m处分界成立，这里面szz的值是怎么算出来的呢？
ini szz 0. grad 0 0 12300. range z -5.5 0.
ini szz 17600 grad 0 0 15500 range z -8 -5.5

ini sxx 0. grad 0 0 5271.4 range z -5.5 0.
ini sxx 7542.86 grad 0 0 6642.86 range z -8 -5.5
ini sxx add 31428.6 grad 0 0 5714.3 range z -8 -5.5

ini syy 0. grad 0 0 5271.4 range z -5.5 0.
ini syy 7542.86 grad 0 0 6642.86 range z -8 -5.5
ini syy add 31428.6 grad 0 0 5714.3 range z -8 -5.5

;
water density 1000
water table origin 0,0,-5.5
normal 0 0 -1

ini dens 1550 range z -8 -5.5

hist unbal
;
solve rat 1.e-6
save pile0.sav
;

model elas
range group pile
prop bulk 13.9e9
shear 10.4e9
range group pile

ini dens 2500
range group pile
call find_add.fis
solve rat 1.e-6
save pile1.sav

rest pile1.sav
;调用保存文件
ini state 0
ini xdis 0.0 ydis 0.0 zdis 0.0
;位移清零
apply szz -0.4e6 range z 0.05 0.15 group pile
;桩顶加第一级荷载
solve
save app0.4.sav
print gp disp range id 1
;输出第一级荷载下的桩顶位移，假定桩顶中心的id号为1
apply szz -0.6e6 range z 0.05 0.15 group pile
;桩顶加第二级荷载
solve
save app0.6.sav
print gp disp range id 1
;输出第二级荷载下的桩顶位移
…………………………………………
;依次加载，直到桩破坏

mthinking

现开始都把它们作为土体看啊，自重应力生成以后在改桩的参数，就会产生类似与桩-土相互的作用，从而应力场会发生改变
这是我个人的理解。