FLAC3D V4.0 help

Rajeev

我有一个问题，是关于用剑桥模型模拟桩在垂直和水平的荷载同时作用的的横向变形
我在运行例子里的代码后 算出的结果和给的参考结果不一样
new
set fish autocreate off
gen zone brick size 15 15 15 edge=11
model mech cam-clay
prop shear 9e6 bulk_bound 8.7e7

prop mm 0.77 lambda 0.13 kappa 0.05

prop mpc 1.5e5 mp1 1e3 mv_l 5.3

fix z range z=(-0.01, 0.01)

fix x range x=(-0.01, 0.01)

fix x range x=(10.99, 11.01)

fix y range y=(-0.01, 0.01)

fix y range y=(10.99, 11.01)

init szz -1.93e5 grad=(0, 0, 1.75e4)

init syy -1.93e5 grad=(0, 0, 1.75e4)

init sxx -1.93e5 grad=(0, 0, 1.75e4)


def camclay_ini_p
   local p_z = zone_head
   loop while p_z # null
   local mean_p = -(z_sxx(p_z) + z_syy(p_z) + z_szz(p_z))/3.0 - z_pp(p_z)
    z_prop(p_z,'cam_cp') = mean_p
   p_z = z_next(p_z)
   endloop
   end
@camclay_ini_p
def path
   local p_z = zone_head
   global sp  = z_prop(p_z,'cam_cp')
   global sq  = z_prop(p_z,'cq')
   local sqcr= sp*z_prop(p_z,'mm')
   if sp = 0.0 then
    sp = 1.
   endif
   global lnp   = ln(sp)
   global svol  = z_prop(p_z,'cv')
   global mk    = z_prop(p_z,'bulk')
   global mg    = z_prop(p_z,'shear')
   global cpc   = z_prop(p_z,'mpc')
   end

set gravity 0 0 -10
ini density 1750
solve

sel pile id=1 begin=(5.5, 5.5, 12.0) end=(5.5, 5.5, 4.0) nseg=32

sel pile id=1 prop Emod=8.0e10 Nu=0.30 XCArea=0.0982 &
                   XCJ=3.83e-4 XCIy=1.92e-4 XCIz=1.92e-4 &
                   Per=3.14 &
                   CS_sK=1.3e11 CS_sCoh=0.0 CS_sFric=10.0 &
                   CS_nK=1.3e11 CS_nCoh=0.0 CS_nFric=0.0 &
                   CS_nGap=off

hist add id=1 unbal
hist add id=2 sel node zdisp id=1
hist add id=3 sel node ydisp id=1
hist add id=4 sel pilesel force fx cid=1
hist add id=5 sel pilesel force fy cid=1

sel node fix x y range id 1

sel node init zvel=-0.5e-8  yvel=0.5e-8 range id=1


sel set damp combined
save pile-ini


cycle 5000
sav pile-lax1

restore pile-ini
sel delete link range id 29
sel link id=100 2 target zone
sel link attach ydir=free  zdir=free              range id 100
sel link attach xrdir=free yrdir=free zrdir=free  range id 100
sel link attach xdir=nydeform                     range id 100
sel link constit nydeform 1 area=1.0 k=5.4e11 ycomp=2.22e5 range id 100
cycle 5000
sav pile-lax2

return