我做的接触分析的例子,但是最后的结果是错的?希望高手指点一下

fengerchai

FINISH
/CLEAR
/COM,  Structural
*SET,x,12  ! 接触部位划分的线数
/PREP7  
ET,1,PLANE42
ET,2,TARGE169
ET,3,CONTA175
R,1,,,0.1, 0.1,,
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2.06e11
MPDATA,PRXY,1,,0.3  
save
PCIRC,150,332/2,-85,-95,   ! 画截面
PCIRC,420/2,392/2,-85,-95,
wpoff,0,-181
CSYS,4  
PCIRC,15, ,0,360,   
CSYS,1  
/PNUM,AREA,1
/REPLOT
SAVE
csys,0              ! 生成接触面和目标面
WPAVE,0,0,0
wpoff,0,-332/2   
CSYS,4  
PCIRC,3, ,0,360,
CSYS,0  
WPAVE,0,0,0
CSYS,0  
wpoff,0,-392/2+0.2  
CSYS,4  
PCIRC,3, ,0,360,
save
FLST,2,5,5,ORDE,2           ! 对滚子OVERLAP并用工作台切割生成接触区域最后再粘接
FITEM,2,1   
FITEM,2,-5  
AOVLAP,P51X
FLST,2,4,5,ORDE,4   
FITEM,2,6   
FITEM,2,8   
FITEM,2,10  
FITEM,2,12  
ADELE,P51X, , ,1
wpro,,,-90.000000   
FLST,2,5,5,ORDE,5   
FITEM,2,7   
FITEM,2,9   
FITEM,2,11  
FITEM,2,13  
FITEM,2,16  
ASBW,P51X   
FLST,2,6,5,ORDE,6   
FITEM,2,2   
FITEM,2,4   
FITEM,2,6   
FITEM,2,8   
FITEM,2,12  
FITEM,2,17  
AGLUE,P51X  
save
FLST,5,6,4,ORDE,6                   ! 划分接触区域的线并mesh面
FITEM,5,7   
FITEM,5,10  
FITEM,5,22  
FITEM,5,28  
FITEM,5,34  
FITEM,5,37  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,x,10, , , ,1
!*  
FLST,5,4,4,ORDE,2   
FITEM,5,49  
FITEM,5,-52
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,30, , , , ,1
!*  
MSHAPE,0,2D
MSHKEY,1
!*  
FLST,5,4,5,ORDE,4   
FITEM,5,2   
FITEM,5,4   
FITEM,5,6   
FITEM,5,8   
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
MSHKEY,0
!*  
FLST,5,2,5,ORDE,2   
FITEM,5,12  
FITEM,5,17  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
save
FLST,2,3,5,ORDE,3   ! 对内外圈进行粘接
FITEM,2,1   
FITEM,2,3   
FITEM,2,14  
AGLUE,P51X  
FLST,2,3,5,ORDE,3   
FITEM,2,5   
FITEM,2,10  
FITEM,2,15  
AGLUE,P51X  
save
FLST,5,3,4,ORDE,3   ! 划分内圈接触区域网格
FITEM,5,1   
FITEM,5,21  
FITEM,5,25  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,x,10, , , ,1
!*  
MSHKEY,1
!*  
FLST,5,2,5,ORDE,2   
FITEM,5,1   
FITEM,5,3   
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
save
FLST,5,3,4,ORDE,3   ! 划分外圈接触区域网格
FITEM,5,9   
FITEM,5,33  
FITEM,5,40  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,x,10, , , ,1
!*  
FLST,5,2,5,ORDE,2   
FITEM,5,5   
FITEM,5,10  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
MSHKEY,0
MSHAPE,1,2D
!*  
FLST,5,2,5,ORDE,2   
FITEM,5,14  
FITEM,5,-15
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
SAVE
lsel,s,line,,28           ! 下面是接触对的设置
nsll,,1
cm,cotac1,node
real,1
type,3
esurf
lsel,s,line,,25
nsll,,1
cm,targe1,node
real,1
type,2
esurf
lsel,s,line,,22            
nsll,,1
cm,cotac,node
real,1
type,3
esurf
lsel,s,line,,21
nsll,,1
cm,targe,node
real,1
type,2
esurf
lsel,s,line,,34
nsll,,1
cm,cotac2,node
real,1
type,3
esurf
lsel,s,line,,33
nsll,,1
cm,targe2,node
real,1
type,2
esurf
lsel,s,line,,37
nsll,,1
cm,cotac3,node
real,1
type,3
esurf
lsel,s,line,,40
nsll,,1
cm,targe3,node
real,1
type,2
esurf
ALLSEL,ALL
SAVE
FLST,2,1,4,ORDE,1   !添加约束和载荷
FITEM,2,5   
!*  
/GO
DL,P51X, ,ALL,0
FLST,2,3,5,ORDE,3   
FITEM,2,1   
FITEM,2,3   
FITEM,2,14  
!*  
/GO
DA,P51X,UX,0
FLST,2,1,4,ORDE,1   
FITEM,2,3   
/GO
!*  
SFL,P51X,PRES,1000,
CSYS,1  
FLST,5,6,5,ORDE,6   
FITEM,5,2   
FITEM,5,4   
FITEM,5,6   
FITEM,5,8   
FITEM,5,12  
FITEM,5,17  
ASEL,S, , ,P51X
NSLA,S,1
NPLOT   
FLST,2,4478,1,ORDE,2
FITEM,2,1   
FITEM,2,-4478   
NROTAT,P51X
ALLSEL,ALL
save
FINISH                        ! 求解设置
/SOL
NSUBST,10,50,5  
AUTOTS,1
TIME,1  
/STATUS,SOLU
SOLVE


附件是最后的结果,似乎发生了穿透!