求助：填充墙框架接触分析时，柱墙分离！

xinrenshao

一个简单的单层单跨填充墙框架，在梁端部施加水平位移荷载，计算过程中出现了几次某节点UX方向位移过大的警告，最后得到的变形图，施加位移荷载一侧，填充墙压缩，另一侧与柱分离，基本没有变形。接触分析中，用Targe170单元模拟砼面作为目标面，contac173单元模拟填充墙作为接触面，单元关键选项，采用默认设置。
请高手看看命令流有什么问题！     命令流和变形图如下所示
finish
/clear
/filname,danceng gangxinglianjie,1
/com,structural
/prep7
et,1,solid65
et,2,targe170
et,3,contac173
!定义混凝土材料------------
mp,ex,1,1.36e10
mp,dens,1,2500
mp,prxy,1,0.2
tb,miso,1,,11
tbpt,,0.0002,2.72E6
tbpt,,0.0004,5.15E6
tbpt,,0.0006,7.29E6
tbpt,,0.0008,9.15E6
tbpt,,0.001,1.07E7
tbpt,,0.0012,1.20E7
tbpt,,0.0014,1.30E7
tbpt,,0.0016,1.37E7
tbpt,,0.0018,1.42E7
tbpt,,0.002,1.43E7
tbpt,,0.0033,1.22E7
!定义砌体材料------------------
mp,ex,2,3.685E9
mp,dens,2,1900
mp,prxy,2,0.15
tb,miso,2,,13
tbpt,,0.0002,7.37E5
tbpt,,0.0004,1.38E6
tbpt,,0.0006,1.94E6
tbpt,,0.0008,2.40E6
tbpt,,0.001,2.77E6
tbpt,,0.0011,2.92E6
tbpt,,0.0012,3.04E6
tbpt,,0.0013,3.15E6
tbpt,,0.0014,3.23E6
tbpt,,0.0015,3.28E6
tbpt,,0.0016,3.32E6
tbpt,,0.0017,3.33E6
tbpt,,0.0027,2.18E6
mp,mu,3,0.7     
r,1$r,2$r,3
!建立有限元模型---------------------
k,1,-0.05,0,-0.075
k,2,0.05,0,-0.075
k,3,0.05,0,0.075
k,4,-0.05,0,0.075
k,5,-0.05,1.09,-0.075
k,6,0.05,1.09,-0.075
k,7,0.05,1.09,0.075
k,8,-0.05,1.09,0.075
v,1,2,3,4,5,6,7,8
vgen,2,1,,,1.5,,,,0
k,,-0.41,0,0.075
k,,-0.41,0,-0.075
k,,-0.41,-0.25,-0.075
k,,-0.41,-0.25,0.075
k,,1.91,0,0.075
k,,1.91,0,-0.075
k,,1.91,-0.25,-0.075
k,,1.91,-0.25,0.075
v,17,18,19,20,21,22,23,24
k,,0.05,1.09,0.05
k,,0.05,1.09,-0.05
k,,0.05,0.91,-0.05
k,,0.05,0.91,0.05
k,,1.45,1.09,0.05
k,,1.45,1.09,-0.05
k,,1.45,0.91,-0.05
k,,1.45,0.91,0.05
v,25,26,27,28,29,30,31,32
vglue,all
k,,0.05,0.91,0.02875
k,,0.05,0.91,-0.02875
k,,0.05,0,-0.02875
k,,0.05,0,0.02875
k,,1.45,0.91,0.02875
k,,1.45,0.91,-0.02875
k,,1.45,0,-0.02875
k,,1.45,0,0.02875
v,33,34,35,36,37,38,39,40
wpoff,,0.91
wprota,,90
vsbw,all
wpcsys,-1
wpoff,,,-0.05
vsbw,all
wpoff,,,0.1
vsbw,all
wpcsys,-1
wprota,,,-90
wpoff,,,0.05
vsbw,all
wpoff,,,-0.1
vsbw,all
wpoff,,,-1.4
vsbw,all
wpoff,,,-0.1
vsbw,all
wpcsys
wpstyle
numcmp,volu
vsel,s,,,2,29
vatt,1,,1
esize,0.05
mshape,0,3d
mshkey,1
vmesh,all
allsel,all
vsel,s,,,1
vatt,2,,1
esize,0.04
mshape,0,3d
mshkey,1
vmesh,all
allsel,all
!有限元模型-----------------------------
asel,s,,,13,22,9    !定义目标单元
asel,a,,,93
nsla,s,1
real,3
mat,3
type,2
esurf,all
allsel,all
asel,s,,,3,6,3     !选择墙体面，定义接触单元
asel,a,,,15
nsla,s,1
real,3
type,3
esurf,all
allsel,all
/solu
asel,s,loc,y,-0.25 !定义约束及荷载
asel,a,loc,y,0
da,all,all
nsel,s,loc,y,0.91,1.09 !面上施加位移荷载
nsel,r,loc,x,1.55
d,all,ux,-0.03
allsel,all
!施加位移荷载时非线性分析选项
!cnvtol,U,,0.05,0 !定义位移收敛条件,使用5%，放松条件
!cnvtol,F,,0.05,1  !定义力收敛条件,使用50%，放松条件
outres,all,all !输出每一子步的结果
autots,1 !打开自动时间步控制
lnsrch,1 !打开线性搜索
ncnv,2 !如果不收敛时结束而不退出
neqit,50 !每一子步中方程的迭代次数限值
pred,on !打开预测器
time,1 !定义时间
nsubst,80 !定义子步数
solve !求解

xinrenshao

xinrenshao

有没有哪位高手看看，这是什么原因啊！

陕西王堃

想问一下楼主现在的问题解决了没有，我也在做填充墙，也是出现了分离的问题，但我用的是实体单元，不知道楼主有什么建议没有。