搜索
利用ANSYS进行断裂分析初次试做断裂分析，希望有这方面经验的高手能发表些经验之谈！这个模型由两种材料组成：表面镀层为铝,基底为钢。目的是对表面镀层的剥离过程进行分析。目前这个模型是个假想的简化模型，初步目标是实现剥离过程的模拟。裂纹扩展是通过接触单元生死功能实现的。基层和镀层由接触单元连接，然后通过断裂判断准则确定要杀死的失效的接触单元。第一版(没有加断裂判断准则，强行逐个杀死界面接触单元)：fini/clear/filn,crack1/PREP7!*ET,1,PLANE182!*KEYOPT,1,1,2KEYOPT,1,3,1KEYOPT,1,4,0KEYOPT,1,6,0KEYOPT,1,10,0!*rect,0,100,0,100rect,0,100,100,110lesi,1,,,10lesi,2,,,10esha,2!*MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,210e3MPDATA,PRXY,1,,0.3MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,2,,70MPDATA,PRXY,2,,0.33amesh,1lesi,5,,,10lesi,6,,,2mat,2amesh,2lsel,s,,,3nsll,s,1cm,c1,nodelsel,s,,,5nsll,s,1cm,t1,nodensel,s,loc,xd,all,uxnsel,s,loc,yd,all,uyd,all,uxmp,mu,3,0/COM,CONTACTPAIRCREATION-STARTCM,_NODECM,NODECM,_ELEMCM,ELEMCM,_LINECM,LINECM,_AREACM,AREA/GSAV,cwz,gsav,,tempMP,MU,3,0MAT,3R,3REAL,3ET,2,169ET,3,172R,3,,,100,0.1,0,RMORE,,,1.0E20,0.0,1.0,RMORE,0.0,0,1.0,,1.0,0.5RMORE,0,0.5,1.0,0.0,KEYOPT,3,2,0KEYOPT,3,3,0KEYOPT,3,4,0KEYOPT,3,5,0KEYOPT,3,7,0KEYOPT,3,8,0KEYOPT,3,9,0KEYOPT,3,10,0KEYOPT,3,11,0KEYOPT,3,12,5!GeneratethetargetsurfaceNSEL,S,,,T1CM,_TARGET,NODETYPE,2ESLN,S,0ESURF,ALLCMSEL,S,_ELEMCM!GeneratethecontactsurfaceNSEL,S,,,C1CM,_CONTACT,NODETYPE,3ESLN,S,0ESURF,ALLALLSELESEL,ALLESEL,S,TYPE,,2ESEL,A,TYPE,,3ESEL,R,REAL,,3/PSYMB,ESYS,1/PNUM,TYPE,1/NUM,1EPLOTESEL,ALLESEL,S,TYPE,,2ESEL,A,TYPE,,3ESEL,R,REAL,,3CMSEL,A,_NODECMCMDEL,_NODECMCMSEL,A,_ELEMCMCMDEL,_ELEMCMCMSEL,S,_LINECMCMDEL,_LINECMCMSEL,S,_AREACMCMDEL,_AREACM/GRES,cwz,gsavCMDEL,_TARGETCMDEL,_CONTACT/COM,CONTACTPAIRCREATION-ENDlsel,s,,,7nsll,s,1cm,s1,node!GradientsurfaceloadSFGRAD,PRES,0,X,0,-0.1,sf,all,pres,-0.1nsel,allesel,all!save/solutime,1deltim,1,1,1solve/post1plns,s,1/soluanty,,resttime,1.1ekill,140solve/post1plns,s,1/soluanty,,resttime,1.2ekill,140ekill,139solve/post1plns,s,1/soluanty,,resttime,1.3ekill,140ekill,139ekill,138solve/post1plns,s,1/soluanty,,resttime,1.4ekill,140ekill,139ekill,138ekill,137solve/post1plns,s,1第二版(加了断裂自动判断准则)。整整花了半天时间才调试明白。fini/clear/filn,crack1s1fail=100eknum=0flag=0/PREP7!*ET,1,PLANE182!*KEYOPT,1,1,2KEYOPT,1,3,1KEYOPT,1,4,0KEYOPT,1,6,0KEYOPT,1,10,0!*rect,0,100,0,100rect,0,100,100,110lesi,1,,,10lesi,2,,,10esha,2!*MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,210e3MPDATA,PRXY,1,,0.3MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,2,,70MPDATA,PRXY,2,,0.33amesh,1lesi,5,,,10lesi,6,,,2mat,2amesh,2lsel,s,,,3nsll,s,1cm,c1,nodelsel,s,,,5nsll,s,1cm,t1,nodensel,s,loc,xd,all,uxnsel,s,loc,yd,all,uyd,all,uxmp,mu,3,0/COM,CONTACTPAIRCREATION-STARTCM,_NODECM,NODECM,_ELEMCM,ELEMCM,_LINECM,LINECM,_AREACM,AREA/GSAV,cwz,gsav,,tempMP,MU,3,0MAT,3R,3REAL,3ET,2,169ET,3,172R,3,,,100,0.1,0,RMORE,,,1.0E20,0.0,1.0,RMORE,0.0,0,1.0,,1.0,0.5RMORE,0,0.5,1.0,0.0,KEYOPT,3,2,1KEYOPT,3,3,0KEYOPT,3,4,0KEYOPT,3,5,0KEYOPT,3,7,2KEYOPT,3,8,0KEYOPT,3,9,0KEYOPT,3,10,0KEYOPT,3,11,0KEYOPT,3,12,5!GeneratethetargetsurfaceNSEL,S,,,T1CM,_TARGET,NODETYPE,2ESLN,S,0ESURF,ALLCMSEL,S,_ELEMCM!GeneratethecontactsurfaceNSEL,S,,,C1CM,_CONTACT,NODETYPE,3ESLN,S,0ESURF,ALLALLSELESEL,ALLESEL,S,TYPE,,2ESEL,A,TYPE,,3ESEL,R,REAL,,3/PSYMB,ESYS,1/PNUM,TYPE,1/NUM,1EPLOTESEL,ALLESEL,S,TYPE,,2ESEL,A,TYPE,,3ESEL,R,REAL,,3CMSEL,A,_NODECMCMDEL,_NODECMCMSEL,A,_ELEMCMCMDEL,_ELEMCMCMSEL,S,_LINECMCMDEL,_LINECMCMSEL,S,_AREACMCMDEL,_AREACM/GRES,cwz,gsavCMDEL,_TARGETCMDEL,_CONTACT/COM,CONTACTPAIRCREATION-ENDlsel,s,,,7nsll,s,1cm,s1,node!GradientsurfaceloadSFGRAD,PRES,0,X,0,-0.1,sf,all,pres,-10nsel,allesel,all!save/solutime,1deltim,1,1,1solveesel,s,type,,3*get,ektnum,elem,0,count*dim,ekel,,ektnumesel,allloadnum=8!loadnum=[1,10]*do,j,1,loadnum/post1set,last/soluesel,s,type,,3nslensel,r,loc,y,100cm,t3,nodeesel,all*get,ndnum,node,0,count*get,nd,node,0,num,min*do,i,1,ndnum*GET,s1nd,NODE,nd,S,1*if,s1nd,ge,s1fail,thenflag=1nsel,s,,,ndeslnesel,r,type,,3*get,el,elem,0,num,mineknum=eknum+1ekel(eknum)=el*endifesel,allcmsel,s,t3nd=ndnext(nd)*enddoPARSAV,all,'par1','',''anty,,restPARRES,NEW,'par1','',''*if,flag,eq,0,thencmsel,s,s1SFGRAD,PRES,0,X,0,-0.1*(1+j),sf,all,pres,-10*elseflag=0*endif*if,eknum,ne,0,then*do,i,1,eknumekill,ekel(i)*enddo*endifesel,allnsel,alltime,1+0.1*jdeltim,.01,.001,.1auto,onsolve*enddo