桩加载fix z range z 0.05 0.15 group pile什么意思

hanbing503

new
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 (0.3,0,0) p9 (0,0.3,0) p10 (0.3,0,-5) p11 (0,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 (0.3,0,-5) p9 (0,0.3,-5) p10 (0.3,0,-8) p11 (0,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 cylin end1 (0,0,0) end2 (0,0,-5.1) radius 0.31 &
                       cylin end1 (0,0,0) end2 (0,0,-5.1) radius 0.29 not
interface 2 face range cylin end1 (0,0,-4.9) end2 (0,0,-5.1) radius 0.31
;
gen zone cyl p0 (0,0,6) p1 (0.3,0,6) p2 (0,0,1) p3 (0,0.3,6) &
             p4 (0.3,0,1) p5 (0,0.3,1) &
             size 3 10 6
gen zone cyl p0 (0,0,6.1) p1 (0.3,0,6.1) p2 (0,0,6) p3 (0,0.3,6.1) &
             p4 (0.3,0,6) p5 (0,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
ini z add -6.0 range group pile
;
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 1 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 -0.1 0.1
fix y range y 7.9 8.1
set grav 0 0 -10
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 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
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
def find_add
   top_head  = null
   gp_pnt = gp_head
   loop while gp_pnt # null
      if gp_zpos(gp_pnt) > 0.05 then
         new        = get_mem(2)
         mem(new)   = top_head
         mem(new+1) = gp_pnt
         top_head   = new
      endif
      gp_pnt = gp_next(gp_pnt)
   endloop
end
find_add
solve rat1.e-6
save pile1.sav
;
ini state 0
ini xdis 0 ydis 0 zdis 0
;监测桩承台的竖向荷载
def zs_top
ad=top_head
zftot=0.0
loop while ad # null
  gp_pnt=mem(ad+1)
  zf=gp_zfunbal(gp_pnt)
  zftot=zftot+zf
  ad=mem(add)
endloop
zs_top=zftot/0.1414
end
fix z range z 0.05 0.15 group pile  固定桩顶速度，用速度来确定位移
这个红色部分的是啥意思啊

sdustqh

1# hanbing503

就是固定Z=0.05到Z=0.15范围内，桩体的z方向速度，默认值为0。

Geolieno

有点疑问，在这个桩基的模拟时，桩受竖向荷载后是向下移动的，桩头部分也是有向下的位移，桩头部分给固定了，好像不行吧？

linzhp

使边界的速率恒定,并非固定边界,只有当速率为0时,才是固定~~fix是固定速率的命令,而非固定位移!!