水平受荷桩加载与接触面单元

fangtaolq

请高手指点：
1、能否用实体单元模拟水平受荷桩？与结构单元的主要差别在什么地方？
2、桩顶水平荷载怎么施加？
3、桩和土之间的接触面怎么模拟？
本就上面三个问题做了些简单的尝试：
用实体单元，wrap命令建立的接触面，水平荷载按照集中荷载的平均值施加在桩顶截面所有节点上；
分析过程中的问题：
发现桩有变形，但是桩周土似乎没有变化，是否是接触面参数设置的问题：
new
rest 30l.sav
model e range group 2
prop bulk 8.333e5 shear 3.846e5 range group 2
model elas range group 1
prop bulk 8.333e6 shear 3.846e6 range group 1
model elas range group 3
prop bulk 8.333e6 shear 3.846e6 range group 3
ini dens 1230 range group 2
ini dens 2230 range group 1
ini dens 2230 range group 3
interface 1 wrap 2 1
interface 1 prop kn 3e7  ks 5e8 fric 30 coh 30000
fix x y z range y -1.69 -1.71
fix x range x 0.59  0.61
fix x range x -0.59  -0.61
fix z range z 0.59  0.61
fix z range z -0.59  -0.61
set grav 0 -9.8 0
hist unbal
solve rat 1.e-5
save 30ph.sav
restore 30ph.sav
ini xdis 0 ydis 0 zdis 0
ini stat 0  
ini xvel 0 yvel 0 zvel 0
interface 1 prop kn 1.0349e7  ks 5.6667e6 fric 16 coh 8500
model mohr range group 2
prop bulk 14.4e6  shear 4.8e6 coh 21348 fric 28.4 range group 2
model elas  range group 1
prop bulk 6.6667e8  shear 3.076e8  range group 1
model elas  range group 3
prop bulk 6.6667e8  shear 3.076e8  range group 3
table 1 name load_displacement
def add_load
p_gp=gp_near(0,0.3,0)
loop n(1,12)
app_load=(421.5+(n-1)*255)/117
file_name='add_step'+string(n)+'.sav'
command
apply xforce app_load range y 0.34 0.36 group 3
solve rat 1.e-5
save file_name
endcommand
xtable(1,n)=app_load
ytable(1,n)=gp_xdisp(p_gp)
endloop
end
add_load
plo set plane dip 90 dd 90 ori 0 0 0
plo con ydisp plane
附上计算力命令流，高手帮忙看看接触面的参数设置有问题吗，竖向受力和横向受力接触面的参数可以设置完全一样的参数吗？

dzd

同问 比较苦恼中

cui841213

跟我遇到的问题很相似，我是在桩顶施加水平动荷载