接触面不起作用，自由侵入

wildhh

模型为上下对称，顶部的接触起作用，为什么底部的接触面不起约束作用，接触面可以自由侵入，C:\Users\ROY\Desktop\QQ截图20140808000235接触分析命令流如下：
fini
/clear
/prep7
et,1,182
et,2,169
et,3,171
keyopt,3,10,2
mp,mu,2,0.2
mp,ex,1,2e8
mp,prxy,1,0.46

k,1,0,1.45
k,2,2,1.45
k,3,2,1.74
k,4,4.39,1.74
k,5,4.39,1
k,6,0.4,1
k,7,0.4,-1
k,8,4.39,-1
k,9,4.39,-1.74
k,10,2,-1.74
k,11,2,-1.45
k,12,0,-1.45
l,1,2
l,2,3
l,3,4
l,4,5
l,5,6
l,6,7
l,7,8
l,8,9
l,9,10
l,10,11
l,11,12
l,12,1
al,all

k,13,0,1.74
k,14,5,1.74
lstr,13,14
k,15,0.4,-1.74
k,16,5,-1.74
lstr,15,16

esize,0.06
type,1
mat,1
smrtsize,5
amesh,all

nsel,s,loc,x,0
d,all,all,0
allsel


/com **** Top and bottom rigid targets ****
r,2

type,2
real,2
mat,2

lesize,13,,,1
lesize,14,,,1
lmesh,13,14

/com **** Contact on the planar seal ****
type,3
lsel,s,line,,9,
lsel,a,line,,3,
nsll,s,1
esurf
allsel,all
finish


/solu
allsel
esel,s,ename,,169
nsle
nsel,r,loc,y,1.74
d,all,uy,-0.240

allsel
esel,s,ename,,169
nsle
nsel,r,loc,y,-1.74
d,all,uy,0.240

nlgeom,on
time,1
nsubst,25,2000,5
outres,all,all
allsel
nropt,unsym ! Unsymmetric due to friction
solve

htbbzzg

从命令流看，你只定义了一个接触啊，线 9 与线 3 ？

htbbzzg

    试了一下你的命令流，确实不收敛。
    分析你的接触对，发现下面的接触单元与目标单元的法线方向不是相对的，而是同一方向。这样，计算时不会形成接触，也就不会收敛了。
    感觉你定义接触对的方式不太好。上下两个接触对最好分开定义。你定义为一个接触对，很容易出现上面说的单元方向问题。需要自己修改下面部分的单元法线方向。
    对你的命令流略做了一些修改，修改位置见以 !**** 作为标记的部分，结果如下，可以收敛了：

fini
/clear

/prep7
et,1,182
et,2,169
et,3,171
keyopt,3,10,2
mp,mu,2,0.2
mp,ex,1,2e8
mp,prxy,1,0.46

k,1,0,1.45
k,2,2,1.45
k,3,2,1.74
k,4,4.39,1.74
k,5,4.39,1
k,6,0.4,1
k,7,0.4,-1
k,8,4.39,-1
k,9,4.39,-1.74
k,10,2,-1.74
k,11,2,-1.45
k,12,0,-1.45
l,1,2
l,2,3
l,3,4
l,4,5
l,5,6
l,6,7
l,7,8
l,8,9
l,9,10
l,10,11
l,11,12
l,12,1
al,all

k,13,0,1.74
k,14,5,1.74
lstr,13,14
k,15,0.4,-1.74
k,16,5,-1.74
lstr,15,16

esize,0.06
type,1
mat,1
smrtsize,5
amesh,all

nsel,s,loc,x,0
d,all,all,0
allsel


/com **** Top and bottom rigid targets ****
r,2

type,2
real,2
mat,2

lesize,13,,,1
lesize,14,,,1
lmesh,13,14

/com **** Contact on the planar seal ****
type,3
lsel,s,line,,9,
lsel,a,line,,3,
nsll,s,1
esurf
allsel,all

!****  颠倒下面部分的接触单元的法线方向：
FLST,5,39,2,ORDE,3  
FITEM,5,1622
FITEM,5,-1659   
FITEM,5,1661
ESEL,S, , ,P51X
ESURF,,REVERSE  

finish


/solu
allsel
esel,s,ename,,169
nsle
nsel,r,loc,y,1.74
d,all,uy,-0.240

allsel
esel,s,ename,,169
nsle
nsel,r,loc,y,-1.74
d,all,uy,0.240

nlgeom,on
time,1
nsubst,25,2000,5
outres,all,all
allsel
nropt,unsym ! Unsymmetric due to friction
solve

wildhh

htbbzzg 发表于 2014-8-8 11:50
试了一下你的命令流，确实不收敛。
    分析你的接触对，发现下面的接触单元与目标单元的法线方向不是 ...

问题得以解决，谢谢帮忙