自己做的一个裂纹接触模型，向各位大神请教

dlutdz

下面是命令流，计算单轴压缩裂纹接触条件下的强度因子，但是计算结果和理论解有很大误差，希望大神们看一下，求指导！！！
FINISH
/CLEAR,START
/PREP7


w=0.0762
h=0.1524
a=0.00635
angle=45
fa=10000000
theta=angle*3.1415926/180
e1=5960000000
v=0.15
ET,1,PLANE183  
KEYOPT,1,3,2



MP,EX,1,e1
MP,PRXY,1,v   



K,1,0,0
K,2,w,0
K,3,w,h
K,4,0,h
K,5,0,h/2-w/2/tan(theta)
K,6,w,h/2+w/2/tan(theta)
K,7,w/2+a*sin(theta),h/2+a*cos(theta)
K,8,w/2-a*sin(theta),h/2-a*cos(theta)
K,9,w/2,h/2
K,19,w/2,h/2
K,21,w/2,h
K,22,w/2,0



L,1,22
L,22,2
L,2,6
L,6,7
L,7,9
L,9,8
L,8,5
L,5,1

L,5,8
L,8,7
L,7,6
L,6,3
L,3,21
L,21,4
L,4,5



AL,1,2,3,4,5,6,7,8
AL,7,9,4,10,11,12,13


KWPAVE,       9
WPRO,,,-90
ASBW,1,2
ASBW,2,2


lesize,ALL,w/40
KSCON,7,a/8,1,10,0.75
KSCON,8,a/8,1,10,0.75
MSHAPE,1,2D
MSHKEY,0
AMESH,ALL




DK,22,UX,0
DL,1,,UY,
DL,2,,UY,
ALLSEL,ALL


SFL,11,PRES,fa,
SFL,12,PRES,fa,
lplot,all

/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,0.363   
MAT,1   
R,3
REAL,3  
ET,2,169
ET,3,172
KEYOPT,3,9,0
KEYOPT,3,10,2   
R,3,
RMORE,  
RMORE,,0
RMORE,0

LSEL,S,,,5  
LSEL,A,,,6  
CM,_TARGET,LINE
TYPE,2  
NSLL,S,1
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM

LSEL,S,,,15
LSEL,A,,,17
CM,_CONTACT,LINE
TYPE,3  
NSLL,S,1
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,2  
ESEL,A,TYPE,,3  
ESEL,R,REAL,,3  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,2  
ESEL,A,TYPE,,3  
ESEL,R,REAL,,3  
CMSEL,A,_NODECM
CMDEL,_NODECM   
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM   
CMSEL,S,_AREACM
CMDEL,_AREACM   
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
/COM, CONTACT PAIR CREATION - END   








/SOLU


CSYS,0
CNODEA=NODE(w/2+a*sin(theta),h/2+a*cos(theta),0)
CNODEB=NODE(w/2-a*sin(theta),h/2-a*cos(theta),0)
NSEL,S,,,CNODEA
CM,CKA,NODE

NSEL,S,,,CNODEB
CM,CKB,NODE
ALLSEL,ALL

LOCAL,11,0,w/2+a*sin(theta),h/2+a*cos(theta),0,45,,,1,1
LOCAL,13,0,w/2-a*sin(theta),h/2-a*cos(theta),0,-135,,,1,1

CINT,NEW,1
CINT,TYPE,SIFS
CINT,CTNC,CKA
CINT,NCON,5
CINT,SYMM,OFF
CINT,NORM,11,2
CINT,NEW,2
CINT,TYPE,SIFS
CINT,CTNC,CKB
CINT,NCON,5
CINT,SYMM,OFF
CINT,NORM,13,2
SOLVE

/POST1
PRCINT,1,,K1
PRCINT,1,,K2
PRCINT,2,,K1
PRCINT,2,,K2