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