ANSYS有限元--冲击载荷作用下的动态应力强度因子计算

/PREP7L1=1.2L2=0.5H1=0.3H2=0.1A=0.04SEITA=55*AFUN,DEGET,1,PLANE183KEYOPT,1,3,2MP,EX,1,2.12E11MP,PRXY,1,0.3MP,DENS,1,7800K,1,0,0,0K,2,L2-H2/TAN(SEITA),0,0K,3,L1,0,0K,4,L1,H1,0K,5,L2+(H1-H2)/TAN(SEITA),H1,0K,6,0,H1,0K,7,L2-A*COS(SEITA),H2-A*SIN(SEITA),0K,8,L2+A*COS(SEITA),H2+A*SIN(SEITA),0K,9,L2,H2,0A,1,2,7,8,5,6A,2,7,9,8,5,4,3wpoff,0,H2,0wpro,,90.000000,ASBW,ALLKSCON,7,A/10,1,12,1,KSCON,8,A/10,1,12,1,ESIZE,0.03LESIZE,7,,,15LESIZE,8,,,15LESIEZE,13,,,15LESIEZE,2,,,15LESIEZE,4,,,25LESIZE,18,,,15LESIZE,14,,,40,0.2,,,,1LESIZE,16,,,40,5,,,,1LESIZE,5,,,40,5,,,,1LESIZE,9,,,40,5,,,,1LESIZE,1,,,40,0.2,,,,1LESIZE,11,,,40,0.2,,,,1AMESH,ALL/SOLANTYPE,4NSEL,S,LOC,X,0D,ALL,,,,,,ALL,,,,,ALLSEL,ALLSFL,19,PRES,-10E6,!施加载荷SFL,15,PRES,-10E6,OUTRES,ALL,ALL,CSYS,0CNODEA=NODE(L2-A*COS(SEITA),H2-A*SIN(SEITA),0)NSEL,S,,,CNODEACM,CKA,NODEALLSEL,ALLCNODEB=NODE(L2+A*COS(SEITA),H2+A*SIN(SEITA),0)NSEL,S,,,CNODEBCM,CKB,NODEALLSEL,ALLLOCAL,11,1,L2-A*COS(SEITA),H2-A*SIN(SEITA),0,180+SEITA,,,1,1,LOCAL,12,1,L2+A*COS(SEITA),H2+A*SIN(SEITA),0,SEITA,,,1,1,CINT,NEW,1CINT,TYPE,SIFSCINT,CTNC,CKACINT,NCON,5CINT,SYMM,OFFCINT,NORM,11,2CINT,LISTCINT,NEW,2CINT,TYPE,SIFSCINT,CTNC,CKBCINT,NCON,5CINT,SYMM,OFFCINT,NORM,12,2CINT,LISTALLSEL,ALL*DIM,KAI,ARRAY,400,1,1,,,*DIM,KAII,ARRAY,400,1,1,,,*DIM,KBI,ARRAY,400,1,1,,,*DIM,KBII,ARRAY,400,1,1,,,I=1*DO,T,0.2E-4,1E-3,0.2E-4TIME,TAUTOTS,1NSUBST,1,,,1KBC,1SOLVESAVE!提取裂纹节点处的I型和II型应力强度因子*GET,KA1,CINT,1,CTIP,CNODEA,,2,DTYPE,K1*GET,KA2,CINT,1,CTIP,CNODEA,,2,DTYPE,K2*GET,KB1,CINT,2,CTIP,CNODEB,,2,DTYPE,K1*GET,KB2,CINT,2,CTIP,CNODEB,,2,DTYPE,K2KAI(I,1)=KA1/10E6/SQRT(3.1415*A)KAII(I,1)=ABS(KA2)/10E6/SQRT(3.1415*A)KBI(I,1)=KB1/10E6/SQRT(3.1415*A)KBII(I,1)=ABS(KB2)/10E6/SQRT(3.1415*A)I=I+1*ENDDOSFL,19,PRES,0,SFL,15,PRES,0,*DO,T,0.0012,0.02,0.0002TIME,TAUTOTS,1NSUBST,1,,,1KBC,1SOLVESAVE*GET,KA1,CINT,1,CTIP,CNODEA,,2,DTYPE,K1*GET,KA2,CINT,1,CTIP,CNODEA,,2,DTYPE,K2*GET,KB1,CINT,2,CTIP,CNODEB,,2,DTYPE,K1*GET,KB2,CINT,2,CTIP,CNODEB,,2,DTYPE,K2KAI(I,1)=KA1/10E6/SQRT(3.1415*A)KAII(I,1)=ABS(KA2)/10E6/SQRT(3.1415*A)KBI(I,1)=KB1/10E6/SQRT(3.1415*A)KBII(I,1)=ABS(KB2)/10E6/SQRT(3.1415*A)I=I+1*ENDDO!将计算的动态断裂因子写入都当前工作目录里的KAI,KAII,KBI,KBII文本里*CREATE,TEMPLE,MAC*CFOPEN,KAI,TXT*DO,I,1,1*VWRITE,KAI(I,1)(F12.6)*ENDDO*CFOPEN,KAII,TXT*DO,I,1,1*VWRITE,KAII(I,1)(F12.6)*ENDDO*CFOPEN,KBI,TXT*DO,I,1,1*VWRITE,KBI(I,1)(F12.6)*ENDDO*CFOPEN,KBII,TXT*DO,I,1,1*VWRITE,KBII(I,1)(F12.6)*ENDDO*ENDTEMPLE