资源描述：
《计算三维问题的应力强度因子的程序》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、计算三维问题的应力强度因子的程序/COM,ANSYSMEDIAREL.60(090601)REF.VERIF.MANUAL:REL.60/VERIFY,VM143*CREATE,FRACT,MAC/NOPRNSEL,ALL*GET,N,NODE,,NUM,MAX!CURRENTMAXIMUMNODENUMBERCMSEL,S,CRACKTIP!SELECTTHETIPNODESESLN!ANYELEMENTSATTACHED*GET,ELMAX,ELEM,,NUM,MAX!CURRENTMAXIMUMELEM
2、ENTNUMBER*DO,IEL,1,ELMAX!LOOPONMAXELEMENTELMI=IEL*IF,ELMI,LE,0,EXIT!NOMORESELECTED*GET,ELTYPE,ELEM,ELMI,ATTR,TYPE!GETELEMENTTYPE*IF,ELTYPE,NE,ARG1,CYCLE!CHECKFORSELECTEDELEMENTN3=NELEM(ELMI,3)!GETNODE3(K)*IF,NSEL(N3),LE,0,CYCLE!ITMUSTBESELECTEDN7=NELEM(ELMI
3、,7)!GETNODE7(L)*IF,NSEL(N7),LE,0,CYCLE!ITMUSTALSOBESELECTEDN1=NELEM(ELMI,1)!GETNODE1(I)N2=NELEM(ELMI,2)!GETNODE2(J)N5=NELEM(ELMI,5)!GETNODE5(M)N6=NELEM(ELMI,6)!GETNODE6(N)X3=0.75*NX(N3)!WEIGHTEDPOSITIONOFN3Y3=0.75*NY(N3)Z3=0.75*NZ(N3)X=0.25*NX(N2)+X3!QUARTE
4、RPOINTLOCATION(NODE(R))Y=0.25*NY(N2)+Y3Z=0.25*NZ(N2)+Z3N=N+1!NEXTNODEN10=NN,N10,X,Y,Z!MIDSIDENODELOCATIONX=0.25*NX(N1)+X3Y=0.25*NY(N1)+Y3Z=0.25*NZ(N1)+Z3N=N+1N12=NN,N12,X,Y,ZX7=0.75*NX(N7)Y7=0.75*NY(N7)Z7=0.75*NZ(N7)X=0.25*NX(N6)+X7Y=0.25*NY(N6)+Y7Z=0.25*NZ
5、(N6)+Z7N=N+1N14=NN,N14,X,Y,ZX=0.25*NX(N5)+X7Y=0.25*NY(N5)+Y7Z=0.25*NZ(N5)+Z7N=N+1N16=NN,N16,X,Y,ZN4=N3N8=N7NSEL,ALLTYPE,3EN,ELMI,N1,N2,N3,N4,N5,N6,N7,N8!REDEFINETHEELEMENTEMORE,0,N10,0,N12,0,N14,0,N16EMORE,*ENDDOCMSEL,U,CRACKTIP!UNSELECTTHETIPNODESNUMMRG,NO
6、DE!MERGEMIDSIDENODESNSEL,ALL!SELECTALLELEMENTSESEL,ALL!SELECTALLELEMENTS/GOPR*END/PREP7*afun,degInnerRadius=0.1!InnerRadius为裂纹半径OuterRadius=1!OuterRadius为圆柱半径Scaler=0.025!Scaler为裂纹前沿单元范围，已经是最佳BaseHeight=0.51!BaseHeight为基层高度,itmaybethebestLayerHeight=0.18!La
7、yerHeight为扩展层高度,差别越大可能越好LayerAmount=16!LayerAmount为层数,与精度关系不大RotationAngle=6!RotationAngle为单元旋转的角度Rotationtimes=90/RotationAngle!Rotationtimes为旋转的次数Height=LayerAmount*LayerHeight+BaseHeight!Height为总高度SMRT,OFF/TITLE,VM143,FRACTUREMECHANICSSTRESSINTENSITY-CRA
8、CKINAFINITEWIDTHPLATEC***BROWNANDSRAWLEY,ASTMSPECIALTECHNICALPUBLICATIONNO.410./COM,******CRACKIN3-DIMENSIONSUSINGSOLID45ANDSOLID95ANTYPE,STATIC!STATICANALYSISET,1,SOLID45ET,2,SOLID45!ELEMENTSAROUNDTHE
