hi nvmhumg,

ofc, here's the project code file that i've made.

The plot between axial stress vs axial displacement indicates that the maximum capacity is 4e5 N/m2 (i.e. 314kN, considering the cross section area of the pile is 0,785m2), while my calculation gives a pull out capacity of 550kN (considering self weight, friction resistance, and point resistance). Observing the contour of displacement for the soil, my suspicion is that the point resistance is not well generated on the tip of the pile.

I think I have built a right model since I wrote it based on the example from FLAC3D manual. I would like to hear any suggestion to improve / correct it.

Thank you.

<PRE>

new

gen zone radcyl p0 (0,0,0) p1 (8,0,0) p2 (0,0,-3) p3 (0,8,0) p4 (8,0,-3) p5(0,8,-3) &

p6 (8,8,0) p7 (8,8,-3) p8 (0.5,0,0) p9 (0,0.5,0) p10 (0.5,0,-3) &

p11(0,0.5,-3) size 3 3 6 8 rat 1 1 1 1.15

gen zone radcyl p0 (0,0,-3) p1 (8,0,-3) p2 (0,0,-9) p3 (0,8,-3) p4 (8,0,-9) p5(0,8,-9) &

p6 (8,8,-3) p7 (8,8,-9) p8 (0.5,0,-3) p9 (0,0.5,-3) p10 (0.5,0,-9) &

p11(0,0.5,-9) size 3 6 6 8 rat 1 1 1 1.15 fill

gen zone reflect dip 90 dd 270

group zone clay

interface 1 face range cyl end1 (0,0,0) end2 (0,0,-3.1) rad 0.51 &

cyl end1 (0,0,0) end2 (0,0,-3.1) rad 0.49 not

interface 2 face range cyl end1 (0,0,-2.9) end2 (0,0,-3.1)rad 0.51

gen zone cyl p0 (0,0,4.1) p1 (0.5,0,4.1) p2 (0,0,4) p3 (0,0.5,4.1) size 3 1 6

gen zone cyl p0 (0,0,4) p1 (0.5,0,4) p2 (0,0,1) p3 (0,0.5,4) size 3 3 6

gen zone reflect dip 90 dd 270 range z 1 4.1

group zone pile range z 1 4.1

save geo0

ini zpos add -4.0 range group pile

save geo1

model mech mohr range group clay

prop bulk 8.3e7 shear 3.8e7 coh 30000 fric 0 range group clay

model mech elas range group pile

prop bulk 8.3e7 shear 3.8e7 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

fix z range z -9.1 -8.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 -9 0

ini sxx 0. grad 0 0 5271.4 range z -9 0

ini syy 0. grad 0 0 5271.4 range z -9 0

solve rat 1e-6

save ini_cond

model mech elas range group pile

prop bulk 13.9e9 she 10.4e9 range group pile

ini dens 2500 range group pile

;find grid points at pile cap

def find_add

top_head = null

gp_pnt = gp_head

loop while gp_pnt # null

if gp_zpos(gp_pnt) > 0.05 then

mem_head = get_mem(2)

mem(mem_head) = top_head

mem(mem_head+1) = gp_pnt

top_head = mem_head

endif

gp_pnt = gp_next(gp_pnt)

endloop

end

@find_add

solve rat 1e-6

save pile_set

ini state 0

ini xdis 0 ydis 0 zdis 0

;monitor vertical loading (pull out) at pile cap

def pull_out

ad = top_head

zftot = 0.0

loop while ad # null

gp_pnt = mem(ad+1)

local zf = gp_zfunbal(gp_pnt)

zftot = zftot + zf

ad = mem(ad)

endloop

pull_out = zftot/0.3927

end

@pull_out

;apply velocity of 1e-8 using a ramp function to the pile top gridpoints

def ramp

nvel = 1e-7

while_stepping

if step < 30000

udapp = (float(step)/30000)*nvel

else

udapp = nvel

endif

ad = top_head

loop while ad # null

gp_pnt = mem(ad+1)

gp_zvel(gp_pnt) = udapp

ad = mem(ad)

endloop

end

@ramp

hist add gp zdisp 0,0,0

hist add fish @pull_out

set mech damp comb

step 300000

save pile1

plot create view capacity

plot add hist -2 vs 1

plot create view displacement

plot add contour displacement range group clay

plot add displacement colorbymag

</PRE>