接触问题

ldhltw

各位朋友好。我最近在算O形密封圈，下面是我的文件。问题：计算老是不收敛；当不施加耦合约束时，发生很大的刚体运动。请大家帮我看下，错在哪里哈。小弟万分感谢。
fini
/cle
/prep7
et,1,plane183, , ,2
KEYOPT,1,6,1                   ! Mixed U-P FORMULATION
et,2,plane42
KEYOPT,2,6,1  
ET,3,TARGE169                  ! 2-D TARGET ELEMENT
ET,4,CONTA172                  ! 2-D CONTACT ELEMENT
KEYOPT,4,5,4
KEYOPT,4,10,1
MP,EX,1,14.04                  ! YOUNG'S MODULUS [MPA]
MP,NUXY,1,0.49967              ! POISSON'S RATIO
C10 = 1.87
C01 = 0.47
NU1 = 0.49967
DD  = 2*(1-2*NU1)/(C10+C01)
TB,HYPER,1,1,2,MOONEY
TBDATA,1,C10,C01,DD
mp,ex,2,280               !滑块--聚四氟乙烯
mp,nuxy,2,0.4
mp,ex,3,2e5               !轴套及轴-钢材料
mp,nuxy,3,0.3
k,1
k,2,19.62,0
k,3,19.62,2.5
k,4,15.6,2.5
k,5,15.6,10.5
k,6,19.62,10.5
k,7,19.62,12.5
k,8,15.6,12.5
k,9,15.6,20.5
k,10,19.62,20.5
k,11,19.62,23
k,12,6.16,23
!-----------o形圈圆心（19.1,6.5）（19.1,16.5）
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,1, 12
lwplan,,12,0
csys,4
k,13,0,0,-4
k,14,0,0,27.81
l,13,14
csys,0
wpcsys,-1
k,15,-5.034832, -3.863822      
k,16,  3.194672, 26.86322  
     
k,17,22.6,-3.863822
k,18,26.6,-3.863822
k,19,22.6,26.86322
k,20,26.6,26.86322
k,21,22.6,12
!------------------
FLST,2,12,4
FITEM,2,1   
FITEM,2,2   
FITEM,2,3   
FITEM,2,4   
FITEM,2,5   
FITEM,2,6   
FITEM,2,7   
FITEM,2,8   
FITEM,2,9   
FITEM,2,10  
FITEM,2,11  
FITEM,2,12  
AL,P51X
a,13,14,16,15
a,17,21,19,20,18
!-------------------
CYL4,19.1,6.5,3.5
CYL4,19.1,16.5,3.5
!
FLST,3,5,5,ORDE,2   
FITEM,3,1   
FITEM,3,-5  
AGEN, ,P51X, , ,135.38, , , , ,1
!---------------------o形圈圆心（19.1,6.5）（19.1,16.5）
wpoff,135.38+19.1,6.5
asel,s,area,,4
lsla,s
wprota,,90
wprota,,,45
lsbw,all
wprota,,,90
lsbw,all
wpcsys,-1
wpoff,135.38+19.1,16.5
asel,s,area,,5
lsla,s
wprota,,90
wprota,,,45
lsbw,all
wprota,,,90
lsbw,all
!*********************
alls
type,1
mat,1
esize,0.5
amesh,4,5
type,2
mat,2
esize,0.8
amesh,1
type,2
mat,3
esize,0.8
amesh,2,3
!---------建立接触对
r,1
real,1
type,3
lsel,s,line,,3
nsll,,1
esln,s,0
esurf,all
type,4
lsel,s,line,,33,35,2
NSLL,,1
esln,s,0
esurf,all
r,2
ET,5,TARGE169                  ! 2-D TARGET ELEMENT
ET,6,CONTA172                  ! 2-D CONTACT ELEMENT
KEYOPT,6,5,4
KEYOPT,6,10,1
real,2
type,5
lsel,s,line,,4
NSLL,,1
esln,s,0
esurf,all
type,6
lsel,s,line,,32,34,2
NSLL,,1
esln,s,0
esurf,all
r,3
ET,7,TARGE169                  ! 2-D TARGET ELEMENT
ET,8,CONTA172                  ! 2-D CONTACT ELEMENT
KEYOPT,8,5,4
KEYOPT,8,10,1
real,3
type,7
lsel,s,line,,17
nsll,s,1
esln,s,0
esurf,all
type,8
lsel,s,line,,22,30,8
nsll,s,1
esln,s,0
esurf,all
r,4
ET,9,TARGE169                  ! 2-D TARGET ELEMENT
ET,10,CONTA172                  ! 2-D CONTACT ELEMENT
KEYOPT,10,5,4
KEYOPT,10,10,1
real,4
type,9
lsel,s,line,,7
nsll,s,1
esln,s,0
esurf,all
type,10
lsel,s,line,,37,39,2
nsll,s,1
esln,s,0
esurf,all
r,5
ET,11,TARGE169                  ! 2-D TARGET ELEMENT
ET,12,CONTA172                  ! 2-D CONTACT ELEMENT
KEYOPT,12,5,4
KEYOPT,12,10,1
real,5
type,11
lsel,s,line,,8
nsll,s,1
esln,s,0
esurf,all
type,12
lsel,s,line,,36,38,2
nsll,s,1
esln,s,0
esurf,all
r,6
ET,13,TARGE169                  ! 2-D TARGET ELEMENT
ET,14,CONTA172                  ! 2-D CONTACT ELEMENT
KEYOPT,14,5,4
KEYOPT,14,10,1
real,6
type,13
lsel,s,line,,9
nsll,s,1
esln,s,0
esurf,all
type,14
lsel,s,line,,25,28,3
nsll,s,1
esln,s,0
esurf,all
r,7
ET,15,TARGE169                  ! 2-D TARGET ELEMENT
ET,16,CONTA172                  ! 2-D CONTACT ELEMENT
KEYOPT,16,5,4
KEYOPT,16,10,1
real,7
type,15
lsel,s,line,,18
nsll,s,1
esln,s,0
esurf,all
type,16
lsel,s,line,,23,26,3
nsll,s,1
esln,s,0
esurf,all
r,8
ET,17,TARGE169                  ! 2-D TARGET ELEMENT
ET,18,CONTA172                  ! 2-D CONTACT ELEMENT
KEYOPT,18,5,4
KEYOPT,18,10,1
real,8
type,17
lsel,s,line,,5
nsll,s,1
esln,s,0
esurf,all
type,18
lsel,s,line,,24,31,7
nsll,s,1
esln,s,0
esurf,all
r,9
ET,19,TARGE169                  ! 2-D TARGET ELEMENT
ET,20,CONTA171                  ! 2-D CONTACT ELEMENT
KEYOPT,20,5,4
KEYOPT,20,10,1
real,9
type,19
lsel,s,line,,12
nsll,s,1
esln,s,0
esurf,all
type,20
lsel,s,line,,13
nsll,s,1
esln,s,0
esurf,all
alls
cp,,all,,38,1645    !不耦合时，O形圈发生刚体位移。
cp,,all,,823,1670
fini
/solu
lsel,s,,,15
nsll,s
d,all,all
lsel,s,,,1
nsll,s
d,all,all
lsel,s,,,20
nsll,s
d,all,ux,-0.5
lsel,s,,,21
nsll,s
d,all,uy
alls
ANTYPE,STATIC
NLGEOM,ON                     ! INCLUDE LARGE DEFORMATION EFFECTS
NSUBST,6                      ! SPECIFY NUMBER OF SUBSTEPS IN LOAD STEP
OUTRES,,1                     ! WRITE SOLUTION FOR EVERY SUBSTEP
SOLVE
FINISH