非对称2D裂纹应力强度因子计算——报告“裂纹面与当前坐标系的X轴不平行”

qp1012

计算过程中，其他步骤均无问题，只是在后处理过程中，提示“裂纹面与当前坐标系的X轴不平行”，但我所建的局部坐标系与裂纹面绝对是垂直的。为此我特地令裂纹面的的斜率为3.5，目的是防止ansys在计算过程中忽略小数点造成的问题。但计算时仍不能解决问题。特跪求高手解决。
我建的模型的命令流如下：
/filename,2Dcrack
/title,qp
/prep7
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2e10   
MPDATA,PRXY,1,,0.3  
ET,1,PLANE82
!*  
KEYOPT,1,3,2
KEYOPT,1,5,0
KEYOPT,1,6,0
!*  
!裂纹长
a=0.08
!裂纹倾角
B=75/180*3.1415926
!裂纹半角
Q=7.5/180*3.1415926
!裂纹尖端长
L=0.01
rw=0.15
a1=rw*sqrt(3)/2+0.01
b1=rw*1/2+0.01
a6=a1+L*cos(B-Q)
b6=b1+L*sin(B-Q)
a2=a1+L*cos(B+Q)
b2=b1+L*sin(B+Q)

a4=a1+a*cos(B)
b4=b1+3.5*(a4-a1)
a5=a6+(a-2*L*cos(Q))*cos(B)
b5=b6+3.5*(a5-a6)
a3=a2+(a-2*L*cos(Q))*cos(B)
b3=b2+3.5*(a3-a2)
a100=a1-0.01*cos(B)
b100=b1+3.5*(a100-a1)
!裂纹的关键点


k,1,a1,b1
k,2,a2,b2
k,3,a3,b3
k,4,a4,b4
k,5,a5,b5
k,6,a6,b6
k,100,a100,b100

LSTR,      1,      2
LSTR,      2,      3  
LSTR,      3,      4  
LSTR,      4,      5
LSTR,      5,      6  
LSTR,      6,      1
!线连接成面
FLST,2,6,4  
FITEM,2,1   
FITEM,2,2   
FITEM,2,3   
FITEM,2,4   
FITEM,2,5   
FITEM,2,6   
AL,P51X
!area2
RECTNG,-1,1,-1,1,
!area3
CYL4,0,0,rw,

!area subtract
ASBA,2,3
asba,4,1

!compress numbers
NUMCMP,AREA

!define local
CSKP,11,0,1,100,2,1,1,

!set crack tip's singularity
!KSCON,3,0.0025,1,16,0,  
KSCON,1,0.005,1,10,0,

FLST,5,8,4,ORDE,2   
FITEM,5,7   
FITEM,5,-14
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
LESIZE,_Y1, , ,15, , , , ,1
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
AMESH,_Y1   
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
FLST,2,1,4,ORDE,1   
FITEM,2,9   
SFL,P51X,PRES,8,
FLST,2,1,4,ORDE,1   
FITEM,2,8   
SFL,P51X,PRES,12,
FLST,2,1,4,ORDE,1   
FITEM,2,7   
SFL,P51X,PRES,8,
FLST,2,1,4,ORDE,1   
FITEM,2,10  
SFL,P51X,PRES,12,

FLST,2,4,4,ORDE,2   
FITEM,2,11  
FITEM,2,-14
/GO
SFL,P51X,PRES,5,
  
/SOL
ANTYPE,0
/STATUS,SOLU
SOLVE   

/POST1  
FLST,2,5,1  
FITEM,2,242
FITEM,2,246
FITEM,2,305
FITEM,2,250
FITEM,2,258
!*  
PATH,L1,5,30,20,
PPATH,P51X,1
PATH,STAT   
!*  

CSYS,11,
!*  
RSYS,11
AVPRIN,0,0  
AVRES,2,
/EFACET,1   
LAYER,0
FORCE,TOTAL
!*  
KCALC,0,1,3,0

qp1012

怎么没有高手来？