ANSYS裂縫擴(kuò)展模擬程序.doc

！裂紋擴(kuò)展程序/prep7 et,1,plane82mptemp, mptemp,1,0 mpdata,ex,1,2.07e5mpdata,prxy,1,0.3 k,101,0,0, k,102,200,0, k,103,200,50, k,104,0,50,k,105,100,50,k,106,80,0,k,107,80,0, k,108,80,10, l,106,108 l,103,104 lextnd,1,108,100,1 lsbl,2,3lplot a,101,106,108,2,104 a,107,102,103,105,2,108 kscon,108,0.2,1,08,1, mshape,1,2d !自由劃分網(wǎng)格mshkey,0amesh,1amesh,2 dk,101,all dk,102,all fk,105,fy,-200finish /solsolve finish nplot local,11,0,80,10,0,-90, , ,1,1, ！建立以裂尖為原點(diǎn)的局部坐標(biāo)系，裂紋面平行于x軸nsel,s,loc,y,0 !在cs坐標(biāo)系(當(dāng)前坐標(biāo)系)下選擇y=0的結(jié)點(diǎn)nsel,r,loc,x,0,0.2 !根據(jù)kscon命令可知沿裂紋面路徑的x坐標(biāo)nlist !共選擇5個(gè)點(diǎn)*get,nnub,node,count !獲取節(jié)點(diǎn)總數(shù)nnub *dim,aa,array,nnub*get,aa(1),node,num,min*do,i,2,nnub*get,aa(i),node,aa(i-1),nxth*enddoflst,2,nnub,1 *do,i,1,nnubfitem,2,aa(i)*enddopath,101,nnub,30,20, ！定義路徑ppath,p51x,1path,stat kcalc,1,1,3,0 ！計(jì)算裂紋強(qiáng)度因子finish *get,ki1,kcalc,k,1 !分別讀取應(yīng)力強(qiáng)度因子ki1,ki2*get,ki2,kcalc,k,2*afun,deg !改變角度的單位為度theta1=2*atan(0.25*(ki1/ki2-sqrt(ki1/ki2)*2+8) !計(jì)算裂紋擴(kuò)展角 坐標(biāo)the=abs(theta1)*dim,kt,array,20,3kt(1,1)=ki1kt(1,2)=ki2kt(1,3)=theta1*do,i,2,10 !循環(huán)的步數(shù)aclear,1,2 !清除網(wǎng)格劃分adele,1,2,0 !僅刪除面1,2，保留線（如果直接刪除線，關(guān)鍵點(diǎn)也將隨之刪除）ldele,all,0 !刪除所有的線kdele,2 !刪除關(guān)鍵點(diǎn)2*get,kmax,kp,num,max !獲取最大關(guān)鍵點(diǎn)編號(hào)k,kmax+2(i-1)-1,0,0 !在原局部坐標(biāo)系下增加一個(gè)關(guān)鍵點(diǎn)，與之前裂尖關(guān)鍵點(diǎn)坐標(biāo)相同k,kmax+2(i-1),-3*cos(the),3*sin(the) !在原局部坐標(biāo)系下增加另一個(gè)關(guān)鍵點(diǎn)，即為下一個(gè)裂尖關(guān)鍵點(diǎn)kdele,kmax !因?yàn)橐粋€(gè)面上不能定義兩個(gè)KSCON，故需要?jiǎng)h除前面在最大號(hào)關(guān)鍵點(diǎn)定義的KSCON,自動(dòng)刪除局部坐標(biāo)*get,kx,kp,kmax+2(i-1)-1,loc,x*get,ky,kp,kmax+2(i-1)-1,loc,yk,kmax,kx,kycsdele，11*get,knub,kp,count !獲取關(guān)鍵點(diǎn)的數(shù)目*dim,kk,array,knub !建立一個(gè)一維矩陣*get,kk(1),kp,num,max !獲取最大關(guān)鍵點(diǎn)編號(hào)*do,j,2,knub*get,kk(j),kp,kk(j-1),nxtl*enddol,103,104 !連接關(guān)鍵點(diǎn)，自動(dòng)編號(hào)為1，用于截取線段（如果角度過(guò)大，需連接l,102,103,并進(jìn)行判定與哪條直線相交）*do,ii,106,kk(1)l,ii,ii+2 !連接關(guān)鍵點(diǎn)成線ii=ii+1*enddo*get,lmax,line,num,max !獲取線的最大編號(hào)lextnd,lmax,kk(1),200,1 !從裂尖方向延伸最大編號(hào)的線，自動(dòng)生成關(guān)鍵點(diǎn)1*get,lmax,line,num,max !獲取線的最大編號(hào)lsbl,1,lmax !用線1截取生成的延伸線，在1上自動(dòng)生成關(guān)鍵點(diǎn)2，同時(shí)刪除關(guān)鍵點(diǎn)1!通過(guò)判定關(guān)鍵點(diǎn)奇偶數(shù)來(lái)確定合成面需要的點(diǎn)*do,j,2,knub *if,kk(j),gt,108,and,kk(j),lt,kmax,then *if,abs(nint(kk(j)/2)-kk(j)/2),le,0.1,then !判定奇偶數(shù)命令，如果為偶數(shù) *get,kpos,kp,count *dim,os,array,kpos *get,os(1),kp,num,max *do,jj,2,kpos *get,os(jj),kp,os(jj-1),nxtl *enddo flst,2,kpos+7,3 fitem,2,108 fitem,2,106 fitem,2,101 fitem,2,104 fitem,2,105 !必須選擇關(guān)鍵點(diǎn)105，否則荷載無(wú)法施加 fitem,2,2 fitem,2,kk(1) *do,jjj,1,kpos fitem,2,os(jjj) jjj=jjj+1 *enddo a,p51x !利用關(guān)鍵點(diǎn)連接成面 *else *get,kpjs,kp,count *dim,js,array,kpjs *get,js(1),kp,num,max *do,jj,2,kpos *get,js(jj),kp,js(jj-1),nxtl*enddo flst,2,kpjs+5,3 fitem,2,107 fitem,2,102 fitem,2,103 fitem,2,2 fitem,2,kk(1) !必須選擇關(guān)鍵點(diǎn)105，否則荷載無(wú)法施加 *do,jjj,1,kpjs fitem,2,js(jjj) jjj=jjj+1 *enddo a,p51x !利用關(guān)鍵點(diǎn)連接成面 *endif *endif*enddokscon,kk(1),0.2,1,08,1, mshape,1,2d !自由劃分網(wǎng)格mshkey,0amesh,1amesh,2 dk,101,all dk,102,all fk,105,fy,-200FINISH /sol/status,solusolve finish *get,kx,kp,loc,kk(1),x !讀取裂尖關(guān)鍵點(diǎn)的x坐標(biāo)*get,ky,kp,loc,kk(1),y !讀取裂尖關(guān)鍵點(diǎn)的y坐標(biāo)local,11,0,kx,ky,0,-90+theta(i-1), , ,1,1, !旋轉(zhuǎn)角度待修改nsel,s,loc,y,0 !在CS坐標(biāo)系(當(dāng)前坐標(biāo)系)下選擇y=0的結(jié)點(diǎn)nsel,r,loc,x,0,0.2 !根據(jù)kscon命令可知沿裂紋面路徑的x坐標(biāo)nlist !共選擇5個(gè)點(diǎn)*get,nnub,node,count !獲取節(jié)點(diǎn)總數(shù)nnub*dim,aa,array,nnub*get,aa(1),node,num,min*do,i,2,nnub*get,aa(i),node,aa(i-1),nxth*enddoflst,2,nnub,1 *do,i,1,nnubfilem,2,aa(i)*enddopath,101,nnub,30,20, ppath,p51x,1path,stat kcalc,1,1,3,0 finish *get,ki1,kcalc,k,1 !分別讀取應(yīng)力強(qiáng)度因子ki1,ki2*get,k