资源描述：
《计算三维问题的应力强度因子的程序《ansys工程有限元分析范例精选》中第24章 基于ansys的三维应力强度因子》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、计算三维问题的应力强度因子的程序文章来源：hotmold.com文章作者：发布时间：2006-05-28 字体：[大中小] /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!
2、CURRENTMAXIMUMELEMENTNUMBER*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,7)!GETNODE
3、7(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!QUARTERPOINTLOCATION(NODE(R))Y=0.25*NY(N2)Y3Z=0
4、.25*NZ(N2)Z3N=N1!NEXTNODEN10=NN,N10,X,Y,Z!MIDSIDENODELOCATIONX=0.25*NX(N1)X3Y=0.25*NY(N1)Y3Z=0.25*NZ(N1)Z3N=N1N12=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(N6)Z7N=N1N14=NN,N14,X,Y,ZX=0.25*NX(N5)X7Y=0.25*NY(N5)Y7Z=0.25*NZ(N5)Z7N=N1N16=
5、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,NODE!MERGEMIDSIDENODESNSEL,ALL!SELECTALLELEMENTSESEL,ALL!SELECTALLELEMENTS/GOPR*END/PREP7*afun,degInnerRadius=0.1!In
6、nerRadius为裂纹半径OuterRadius=1!OuterRadius为圆柱半径Scaler=0.025!Scaler为裂纹前沿单元范围，已经是最佳BaseHeight=0.51!BaseHeight为基层高度,itmaybethebestLayerHeight=0.18!LayerHeight为扩展层高度,差别越大可能越好LayerAmount=16!LayerAmount为层数,与精度关系不大RotationAngle=6!RotationAngle为单元旋转的角度Rotationtimes=90/RotationAngle!Rotationtimes为
7、旋转的次数Height=LayerAmount*LayerHeightBaseHeight!Height为总高度SMRT,OFF/TITLE,VM143,FRACTUREMECHANICSSTRESSINTENSITY-CRACKINAFINITEWIDTHPLATEC***BROWNANDSRAWLEY,ASTMSPECIALTECHNICALPUBLICATIONNO.410./COM,******CRACKIN3-DIMENSIONSUSINGSOLID45ANDSOLID95ANTYPE,STATIC!STATICANALYSISET,1,SOLID45
