大家帮忙看看我做的抗拔桩模型为什么出问题
想做一个简单的抗拔桩模型,但是最后加载总是发生单元过度变形的错误而退出。大家帮我看看出了什么问题吧。说明一下1.桩土截面摩擦角为0,是我特意设置的。我想看下摩擦角为0的时候,其他参数的关系。2、应该不是所加荷载过大的原因,因为我只加1牛的力,也出错。
new
set log on
;生成土体
gen zone radbrick p0 0 0 0 &
p1 0.5 0 0 &
p2 0 0 -1 &
p3 0 0.5 0 &
dimension 0.25 0.8 0.25 &
size 5 16 5 5 &
ratio 1 1 1 1.15 &
fill group inner
;映射
generate zone reflect normal -1 0 0 origin 0 0 0 range z -1 0
generate zone reflect normal 0 -1 0 origin 0 0 0 range z -1 0
plot show
plot ske
attach face
group subgrade_block2 range z -1,0
;设置土体材料参数
model cam-clay range group subgrade_block2
prop bulk_bound=352101384poisson=0.36 range group subgrade_block2
prop mm=1.07 lambda=0.0186 kappa=0.0067 range group subgrade_block2
prop mpc=160709.645 mp1=1 mv_l=1.9134 range group subgrade_block2
ini dens=1860.0 range group subgrade_block2
save soilmap.sav
;设置边界条件
fix x y z range z -1.001 -0.999
fix x y range x -0.501 -0.499
fix x y range y -0.501 -0.499
fix x y range x 0.501 0.499
fix x y range y 0.501 0.499
set grav 0 0 -10
set mech ratio 5e-5
;initaial stress
def ini_szz
; --- layer 2 is dry ---
_grad = -18600 ; (a positive value)
_gradfu =18600
_szzl = _grad*1
_szz0 = 0
command
ini szz add _szz0 grad 0 0 _gradfu range z -1 0
end_command
end
ini_szz
def ini_conf
pnt = zone_head
loop while pnt # null
val =0.65*z_szz(pnt)
z_sxx(pnt) = val
z_syy(pnt) = val
pnt=z_next(pnt)
end_loop
end
ini_conf
def camclay_ini_p
p_z = zone_head
loop while p_z # null
mean_p = -(z_sxx(p_z) + z_syy(p_z) + z_szz(p_z))/3.0
z_prop(p_z,'cam_cp') = mean_p
p_z = z_next(p_z)
endloop
end
camclay_ini_p
hist id=1 unbal
solve
save subgrade_ph0.sav
ini state 0
ini xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0
;生成桩,设置参数
sel pile id=1 begin=(0,0,0) end=(0,0,-0.5) nseg=10
sel pile prop emod=3e10 Nu=0.2 XCAREA=0.196349538 XCJ=0.006135923 XCIy=0.003067962 &
XCIz=0.003067962 cs_sk=1.25e11 cs_nK=1.6667e11 cs_sCoh=1.5e5 cs_nCoh=1.5e5 &
cs_sFric=0 cs_nfric=0 &
cs_ngap=off per=0.1570795 dens=2500
set logfile 10.log
;施加荷载
sel node apply force 0 0 10.0 range id=1
set large
solve
save pile_load10.sav
页:
[1]