非对称2D裂纹应力强度因子计算——报告“裂纹面与当前坐标系的X轴不平行”
计算过程中，其他步骤均无问题，只是在后处理过程中，提示“裂纹面与当前坐标系的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 怎么没有高手来？
页:
[1]