ANSYS有限元--界面裂纹能量释放率的计算

/PREP7
L=1
L1=0.6
H=0.3
H1=0.1
A=0.05
ET,1,PLANE183
KEYOPT,1,3,3
R,1,0.003,
MP,EX,1,2.1E11
MP,PRXY,1,0.3
MP,EX,2,5E11
MP,PRXY,2,0.25
K,1,0,0,0
K,2,L1,0,0
K,3,L,0,0
K,4,L,H,0
K,5,L1,H,0
K,6,0,H,0
K,7,L1,H1-A,0
K,8,L1,H1,0
K,9,L1,H1+A,0
A,1,2,7,9,5,6
A,2,3,4,5,9,8,7
wpoff,0,H1,0
wpro,,90.000000,
ASBW,ALL
KSCON,7,A/15,1,12,1,
KSCON,9,A/15,1,12,1,
ESIZE,0.005
MAT,1
AMESH,5
AMESH,3
MAT,2
AMESH,4
AMESH,6
/SOL
DL,17, ,ALL,
DL,12, ,ALL,
SFL,19,PRES,-2E6,
SFL,15,PRES,-2E6,
CSYS,0
CNODEA=NODE(L1,H1+A,0)
CNODEB=NODE(L1,H1-A,0)
NSEL,S, , ,CNODEA
CM,CKA,NODE
ALLSEL,ALL
NSEL,S, , ,CNODEB
CM,CKB,NODE
ALLSEL,ALL
LOCAL,11,1,L1,H1+A,0,90, , ,1,1,
LOCAL,12,1,L1,H1-A,0,-90, , ,1,1,
CINT,NEW,1
CINT,TYPE,VCCT
CINT,CTNC,CKA
CINT,SYMM,OFF
CINT,NORM,11,2
CINT,LIST
ALLSEL,ALL
CINT,NEW,2
CINT,TYPE,VCCT
CINT,CTNC,CKB
CINT,SYMM,OFF
CINT,NORM,12,2
CINT,LIST
ALLSEL,ALL
SOLVE
/POST1
PRCINT,1,
PRCINT,2,