求助，基坑开挖中碰到的接触面问题？？

绍兴人

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 wall  range 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
附上一些图：

z方向位移图

接触面法向应力

衬砌法向应力

wangweipeng

接触面参数中tens=1e10大的有点离谱，把这个参数调小了再试试，好像还不止这一点毛病。先试试

绍兴人

调小试了，，不过结果一样！！！

无痕20

学习一下，一般地下连续墙很少考虑接触问题！

大浪淘沙420

应该不考虑接触问题吧

try126

学习一下，一般地下连续墙很少考虑接触问题！

try126

向大家学习

ycding

好像这中间有部分程序引用的我的吧，但结果怎么还是不合理？

hjyq520

兄弟怎么联系，有问题请教！

94117239

我认为地连墙还是要考虑接触问题，要不地连墙弯矩是不合理的

94117239

我现在纠结于实体和壳单元之间设置接触面

gaoyuchun

我也做基坑开挖，一般用实体单元做地连墙

shangyu1018

难道不需要建立接触面吗

shangyu1018

请问高手接触面参数一般怎么取？

ycding

搞基坑的有空找找我的贴子哦，

windyxsw

高手~
学习学习~

yangsuxin

88888888888888888888888

yangsuxin

88888888888888888888888

yangsuxin

ppppppppppppppp

yangsuxin

1111111111111