求助，基坑开挖中碰到的接触面问题？？
40×40的连续墙基坑，连续墙深度20m，基坑开挖10m，连续墙采用衬砌，加在坑外土体表面；在坑内土体外表面加了接触面。可计算结果为什么如此呢？？麻烦各位抽空帮忙分析一下。。谢谢！！！关键程序如下：
; --- 创建连续墙 ---
group wall range x _x1 _x2 y 0 _y2 z _z1 _z4
group wall range x 0 _x1 y _y1 _y2 z _z1 _z4
range name diaphragm1 x=(19.9,20.1) y=(-0.1,20.1) z=(19.9,40.1)
range name diaphragm2 x=(-0.1,20.1) y=(19.9,20.1) z=(19.9,40.1)
range name diaphragm diaphragm1 any diaphragm2 any
sel liner id=100 crossdiag group wallrange diaphragm
sel liner prop iso=(24e9,0.2) thick=0.8 density =2.5e3 cs_ncut=1e3 cs_nk=4e8 &
cs_scoh=20e3 cs_scohres=0.0 cs_sfric=0.0 cs_sk=4e8 slide=on &
range diaphragm
sel group diaphragm1 range diaphragm1
sel group diaphragm2 range diaphragm2
sel group diaphragm range diaphragm
sel link attach zd=rigid xd=free yd=free xrdir=free yrdir=free zrdir=free ran sel liner
;-----------------------施加接触面----------------------
group inner range x 0 20 y 0 20 z 20 40
group outer range group inner not
gen separate inner
int 1 face range group inner x=(19.9,20.1)
int 1 face range group inner y=(19.9,20.1)
interface 1 prop kn=4e9 ks=4e9 tens=1e10 coh=15e3 fric=14.0 bslip=on sbratio=20 ;接触面参数
interface 1 maxedge=4
;interface 1 prop kn=4e8 ks=4e8 tens=1e10 sbratio=100
plot set ba wh
pl ske interface red blue attach cyan green
;-----------------------重新定义墙底约束条件----------------------
def redef_wall_end_link
node_pnt= nd_head
link_id=10000
loop while node_pnt # null
node_id = nd_id(node_pnt)
xx = nd_pos(node_pnt,2,1)
yy = nd_pos(node_pnt,2,2)
zz = nd_pos(node_pnt,2,3)
link_pnt= nd_link(node_pnt)
dist_x = sqrt((xx-20.0)^2+(zz-20.0)^2)
dist_y = sqrt((yy-20.0)^2+(zz-20.0)^2)
dist_tol = 1e-3
if dist_x < dist_tol then
if link_pnt # null then
temp1 = lk_delete(link_pnt)\
link_id = link_id+1
command
sel set link node_tol = dist_tol
sel link id=link_id node_id target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=free yrdir=free zrdir=free range id=link_id
endcommand
endif
endif
if dist_y < dist_tol then
if link_pnt # null then
temp2 = lk_delete(link_pnt)
link_id = link_id+1
command
sel set link node_tol = dist_tol
sel link id=link_id node_id target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=free yrdir=free zrdir=free range id=link_id
endcommand
endif
endif
node_pnt = nd_next(node_pnt)
endloop
end
redef_wall_end_link
附上一些图：
http://home.simwe.com/attachment/200906/1/396989_1243842141JzJJ.gif
z方向位移图
http://home.simwe.com/attachment/200906/1/396989_124384213935b3.gif
接触面法向应力
http://home.simwe.com/attachment/200906/1/396989_1243842142PiaR.gif
衬砌法向应力 接触面参数中tens=1e10大的有点离谱，把这个参数调小了再试试，好像还不止这一点毛病。先试试 调小试了，，不过结果一样！！！:'( 学习一下，一般地下连续墙很少考虑接触问题！ 应该不考虑接触问题吧 学习一下，一般地下连续墙很少考虑接触问题！ 向大家学习 好像这中间有部分程序引用的我的吧，但结果怎么还是不合理？ 兄弟怎么联系，有问题请教！:lol 我认为地连墙还是要考虑接触问题，要不地连墙弯矩是不合理的 我现在纠结于实体和壳单元之间设置接触面 我也做基坑开挖，一般用实体单元做地连墙 难道不需要建立接触面吗 请问高手接触面参数一般怎么取？ 搞基坑的有空找找我的贴子哦，:lol 高手~
学习学习~ 88888888888888888888888 88888888888888888888888 ppppppppppppppp 1111111111111
页:
[1]
2