ANSYS分析巴西圆盘接触问题

sunkang0221

D:\Program Files\ANSYS Inc\SUN\BD(OPEN CRACK PROPAGATION )\OPEN CRACK IMAGE\OPEN CRACK GIF.gif
发现在计算过程中圆盘发生了转动或者平移，但是计算时收敛的
请问是网格划分有问题吗？下面是命令流
FINISH
/CLEAR
/TITLE,The Brazilian Disc Uniaxial Compression Test with Contact
/PREP7
/RGB,INDEX,100,100,100, 0                   !屏幕背景色设置为反相（白色）
/RGB,INDEX, 80, 80, 80,13
/RGB,INDEX, 60, 60, 60,14
/RGB,INDEX, 0, 0, 0,15
/REPLOT

D=40                        !单位：mm
L=8
BB=0.3                      ! 单元宽度
H=10
PI=ACOS(-1)
BETA=45
ANGLE=BETA*PI/180                       

!设置材料属性
YOUNG = 1.4E3               !MPa
NU = 0.3                    ! POISSION'S RATIO
COF=0.5
MP,EX,1,YOUNG
MP,PRXY,1,NU
MP,MU,1,COF

MP,EX,2,100*YOUNG
MP,PRXY,2,NU
MP,MU,2,COF

MP,EX,3,YOUNG
MP,PRXY,3,NU
MP,MU,0

!设置单元属性以及实常数
ET,1,PLANE183
KEYOPT,1,3,0
ET,2,TARGE169
ET,3,CONTA172
KEYOPT,3,2,1
KEYOPT,3,9,0
KEYOPT,3,10,1
R,2
REAL,2
R,2,0,0,10,0,0,0
R,3,0,0,10,0,0,0


!建立模型
CSYS,0
K,1,0,0
K,2,0,0
K,3,L/2*COS(ANGLE),L/2*SIN(ANGLE)
K,4,-L/2*COS(ANGLE),-L/2*SIN(ANGLE)
K,5,D/2*COS(ANGLE),D/2*SIN(ANGLE)
K,6,-D/2*COS(ANGLE),D/2*SIN(ANGLE)
K,7,-D/2*COS(ANGLE),-D/2*SIN(ANGLE)
K,8,D/2*COS(ANGLE),-D/2*SIN(ANGLE)

LARC,5,6,1,D/2
LARC,6,7,1,D/2
LARC,7,8,1,D/2
LARC,8,5,1,D/2
L,1,3$L,2,3$L,3,5$L,6,1
L,1,4$L,2,4$L,4,7$L,8,2

AL,5,7,1,8
AL,8,2,11,9
AL,10,11,3,12
AL,6,7,4,12
BLC4,-D/2,D/2,D,H,
BLC4,-D/2,-(D/2+H),D,H,

!进行圆盘网格划分
ASEL,S,,,1,4,1
TYPE,1         
MAT,1
ESIZE,BB
KSCON,3,0.4,1,30,0.5         !裂尖半径可以设稍微大一点，0.5左右
KSCON,4,0.4,1,30,0.5         !定义裂纹尖端
MSHKEY,0
AMESH,ALL

!进行加载平台网格划分
ALLSEL,ALL
ASEL,S,,,5,6,1
TYPE,1
MAT,2
ESIZE,,60
MSHKEY,0
AMESH,ALL


ALLSEL,ALL
!进行裂纹面接触设置
TYPE,3
MAT,3
REAL,2
LSEL,S,,,5,9,4
NSLL,S,1
ESLN,S,0
ESURF
ALLSEL,ALL

TYPE,2
MAT,3
LSEL,S,,,6,10,4
NSLL,S,1
ESLN,S,0
ESURF
ALLSEL,ALL


!进行圆盘与加载平台接触设置
TYPE,3
MAT,1
REAL,3
LSEL,S,,,1
LSEL,A,,,3
NSLL,S,1
ESLN,S,0
ESURF
ALLSEL,ALL

TYPE,2
MAT,2   
LSEL,S,,,13
LSEL,A,,,19
NSLL,S,1
ESLN,S,0
ESURF
ALLSEL,ALL


/SOLU
CSYS,0
LSEL,S,LOC,Y,-D/2-H
DL,ALL,,ALL,0
ALLSEL,ALL
LSEL,S,LOC,Y,-D/2-H/2
DL,ALL,,UX,0
ALLSEL,ALL
LSEL,S,LOC,Y,D/2+H/2
DL,ALL,,UX,0
ALLSEL,ALL

LSEL,S,LOC,Y,D/2+H
DL,ALL,,UY,-0.5
ALLSEL,ALL


SOLCONTROL,1,1,,          !非线性求解设置
OUTRES,ALL,ALL,
NLGEOM,ON                !考虑大变形
NSUBST,20
TIME,1
AUTOTS,ON
CNVTOL,F,2000,0.005,2,
SOLVE
FINISH