分享自己做的一个二维疲劳裂纹SIF的例子，另外找高手帮我看下三维的分析哪里出错了

ethanxu

RT，以前没怎么接触过ANSYS，自己通过看论坛以及学习VM143，好不容易计算出了一个二维试样疲劳裂纹应力强度因子的算例，结果与公式计算值非常之吻合我很高兴，但是做到三维的时候KI与我二维的计算出来的误差很大（举例二维算出来的255，三维的165），APDL如下，求高手帮帮忙，顺便帮新手们找个例子。PS:我二维的是建立的1/2模型，三维的是1/4模型，加载都是孔上的面载荷。

二维：
finish
/clear
/prep7  
et,1,plane82
KEYOPT,1,3,3
KEYOPT,1,5,0
KEYOPT,1,6,0
R,1,3.8,
MP,EX,1,73100  
MP,PRXY,1,0.33
*SET,cl,24.5                !设置裂纹长度
*SET,cl1,0.01*cl
*SET,cl2,2*cl1
k,1,0,0,
k,2,24.15,0
k,3,cl,0  
k,4,52.5,0  
k,5,10.5,13.85  
k,6,0,18
l,1,3   
l,3,4   
l,1,6   
WPSTYLE,,,,,,,,1
KWPAVE,       2
circle,2,28.35,,,180
circle,5,6.25,,,360
LPTN,3,5
/PNUM,KP,1  
/PNUM,LINE,1
LDELE,11, , ,1  
ldele,13,,,1
numcmp,all  
lplot   
allsel,all  
lglue,all   
al,all  
KSCON,3,cl1,1,16,1,  
WPCSYS,-1,0
wpoff,10.5,13.85
wprota,,90  
asbw,all
wprota,,,90
asbw,all
circle,3,3*cl2,,,180
lcomb,1,9   
asbl,6,1
ESIZE,1,0,  
MSHAPE,0,2D
MSHKEY,0
FLST,5,4,5,ORDE,3   
FITEM,5,1   
FITEM,5,4   
FITEM,5,-6  
CM,_Y,AREA  
ASEL, , , ,all
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
AMESH,_Y1   
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
DL,       9,       ,SYMM   
DL,       21,       ,SYMM
*SET,n4,node(10.5,20.1,0)
D,n4, , , , , ,UX, , , , ,
SFL,4,PRES,21,   
SFL,5,PRES,21，
FINISH  
/SOL
ANTYPE,0
/STATUS,SOLU
SOLVE   
FINISH
/POST1  
*SET,n1,node(cl,0,0)  
*SET,n2,node(cl-cl1,0,0)
*SET,n3,node(cl-cl2,0,0)
PATH,K1,3,30,20,
PPATH,1,n1
PPATH,2,n2
PPATH,3,n3
PATH,STAT
WPCSYS,-1,0
wpoffs,cl
CSWPLA,11,0,1,1,
CSYS,11,
RSYS,11
KCALC,1,1,0,0   

三维：
FINISH
/CLEAR
/PREP7  
ET,1,SOLID45
ET,2,PLANE82
MP,EX,1,73100  
MP,PRXY,1,0.33
*SET,CL,24.5
*SET,CL1,0.01*CL
*SET,CL2,2*CL1
K,1,0,0,
K,2,24.15,0
K,3,CL,0  
K,4,52.5,0  
K,5,10.5,13.85  
K,6,0,18
L,1,3   
L,3,4   
L,1,6   
WPSTYLE,,,,,,,,1
KWPAVE,       2
CIRCLE,2,28.35,,,180
CIRCLE,5,6.25,,,360
LPTN,3,5
LDELE,11, , ,1  
LDELE,13,,,1
NUMCMP,ALL  
LPLOT   
ALLSEL,ALL  
LGLUE,ALL   
AL,ALL  
WPCSYS,-1,0
WPOFF,10.5,13.85
WPROTA,,90  
ASBW,ALL
WPROTA,,,90
ASBW,ALL
CIRCLE,3,3*CL2,,,180
LCOMB,1,9
ASBL,6,1
VEXT,ALL, , ,0,0,-1.9,,,,   
KSCON,3,CL1,1,16,1,
WPCSYS,-1,0
ALLSEL
ASEL,S,LOC,Z,0
APLOT
!基准面上画网格
VATT,2,1,,0
SECNUM,
ESIZE,1,0,  
MSHAPE,0,2D
MSHKEY,0
FLST,5,5,5,ORDE,2   
FITEM,5,1   
FITEM,5,-5  
CM,_Y,AREA  
ASEL, , , ,ALL
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
AMESH,_Y1   
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
！扫略体网格
VATT,1,1,,0
SECNUM,
ESIZE,1,0,  
FLST,5,5,6,ORDE,2   
FITEM,5,1   
FITEM,5,-5  
CM,_Y,VOLU  
VSEL, , , ,ALL
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
VSWEEP,_Y1  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
ALLSEL
APLOT
！约束
DA,14,SYMM
DA,18,SYMM
DL,26, ,UX,
ASEL,S,LOC,Z,0
DA,ALL,SYMM
ALLSEL
APLOT
！加载
SFA,7,1,PRES,21,   
SFA,23,1,PRES,21,
FINISH  
/SOL
ANTYPE,0
/STATUS,SOLU
SOLVE   
FINISH
/POST1  
PRRSOL,FY   
*SET,N1,NODE(CL,0,0)  
*SET,N2,NODE(CL-CL1,0,0)
*SET,N3,NODE(CL-CL2,0,0)
PATH,K1,3,30,20,
PPATH,1,N1
PPATH,2,N2
PPATH,3,N3
PATH,STAT
WPCSYS,-1,0
WPOFFS,CL
CSWPLA,11,0,1,1,
CSYS,11,
RSYS,11
KCALC,0,1,0,0

ethanxu

没有人么？