《ansys断裂例子.docx》由会员分享，可在线阅读，更多相关《ansys断裂例子.docx（8页珍藏版）》请在三一文库上搜索。
1、):第一版（没有加断裂判断准则，强行逐个杀死界面接触单元fini /clear /filn,crack1 /PREP7 !*ET,1,PLANE182 !*KEYOPT,1,1,2 KEYOPT,1,3,1 KEYOPT,1,4,0 KEYOPT,1,6,0 KEYOPT,1,10,0 !* rect,0,100,0,100 rect,0,100,100,110 lesi,1,10 lesi,2,10 esha,2 !*MPTEMP, MPTEMP,1,0 MPDATA,EX,1,210e3 MPDATA,PRXY,1,0.3 MPTEMP, MPTEMP,1,0 MPDATA,EX,2,70
2、 MPDATA,PRXY,2,0.33 amesh,1 lesi,5,10 lesi,6,2 mat,2 amesh,2 lsel,s,3 nsll,s,1 cm,c1,node lsel,s,5 nsll,s,1 cm,t1,node nsel,s,loc,x d,all,ux nsel,s,loc,y d,all,uyd,all,ux mp,mu,3,0/COM, CONTACT PAIR CREATION - STARTCM,_NODECM,NODECM,_ELEMCM,ELEMCM,_LINECM,LINECM,_AREACM,AREA /GSAV,cwz,gsav,tempMP,MU
3、,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! Generate the target surfaceNSEL,S,T1CM,_TARGET,NODETYPE,2ESLN,S
4、,0ESURFALLCMSEL,S,_ELEMCM ! Generate the contact surfaceNSEL,S,C1CM,_CONTACT,NODETYPE,3ESLN,S,0ESURFALLALLSELESEL,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
5、,S,_LINECMCMDEL,_LINECMCMSEL,S,_AREACMCMDEL,_AREACM/GRES,cwz,gsavCMDEL,_TARGETCMDEL,_CONTACT/COM, CONTACT PAIR CREATION - ENDlsel,s,7nsll,s,1cm,s1,node!Gradient surface loadSFGRAD,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.1eki
6、ll,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=100/filn,crack1s1fail=100!定义断裂
7、准则：断裂应力eknum=0flag=0/PREP7!*ET,1,PLANE182i*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!设置单元形状 4边形!*MPTEMP,!定义材料 1MPTEMP,1,0MPDATA,EX,1,210e3MPDATA,PRXY,1,0.3MPTEMP,!定义材料 2MPTEMP,1,0MPDATA,EX,2,70MPDATA,PRXY,2,0.33amesh,1l
8、esi,5,10lesi,6,2mat,2amesh,2lsel,s,3!选择线 3nsll,s,1!选择线3上的所有节点cm,c1,node!定义已选择的点为 c1组lsel,s,5!选择线 5nsll,s,1!选择线5上的所有节点cm,t1,node!定义已选择的点为 t1组nsel,s,loc,x!选择X=0的所有节点d,all,ux!nsel,s,loc,y!选择丫=0的所有节点d,all,uy!d,all,uxmp,mu,3,0!定义材料3/COM, CONTACT PAIR CREATION - START ! 定义接触对，采用接触对向导CM,_NODECM,NODECM,_ELE
9、MCM,ELEMCM,_LINECM,LINECM,_AREACM,AREA/GSAV,cwz,gsav,tempMP,MU,3,0MAT,3R,3REAL,3ET,2,169ET,3,172R,3,100,0.1,0, !Normal penalty stiffness=100RMORE,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,
10、0KEYOPT,3,11,0KEYOPT,3,12,5! Generate the target surfaceNSEL,S,T1!目标面T1组节点CM,_TARGET,NODETYPE,2ESLN,S,0ESURFALLCMSEL,S,_ELEMCM! Generate the contact surfaceNSEL,S,C1!接触面C1组节点CM,_CONTACT,NODETYPE,3ESLN,S,0ESURFALLALLSELESEL,ALLESEL,S,TYPE,2ESEL,A,TYPE,3ESEL,R,REAL,3/PSYMB,ESYS,1/PNUM,TYPE,1/NUM,1EPLO
11、TESEL,ALLESEL,S,TYPE,2ESEL,A,TYPE,3ESEL,R,REAL,3CMSEL,A,_NODECMCMDEL,_NODECM CMSEL,A,_ELEMCM CMDEL,_ELEMCM CMSEL,S,_LINECM CMDEL,_LINECM CMSEL,S,_AREACM CMDEL,_AREACM /GRES,cwz,gsav CMDEL,_TARGETCMDEL,_CONTACT /COM, CONTACT PAIR CREATION - ENDlsel,s,7!选择线 7nsll,s,1!选择线7上的节点cm,s1,node!定义节点组S1!Gradien
12、t surface load定义面载荷梯度!加载!选择所有节点!选择所有元素!读取载荷步结果!选择元素类型3!选择元素类型3,y=100的节点!定义节点组t3!选择所有元素!变量ndnum=节点总数（t3组）!变量nd=节点的最小编号!变量s1nd=节点nd的应力s1!断裂判据SFGRAD,PRES,0,X,0,-0.1, sf,all,pres,-10 nsel,all esel,all !save /solu time,1 deltim,1,1,1 solve esel,s,type,3 *get,ektnum,elem,0,count *dim,ekel,ektnum esel,all
13、loadnum=8 !loadnum=1,10 *do,j,1,loadnum /post1 set,last /solu esel,s,type,3 nsle nsel,r ,loc,y,100 cm,t3,node esel,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,then!第一次求解!选择元素类型3!变量ektnum=元素总素!定义一 21数组ekel!选择所有元素flag=1nsel,s,ndesln!选择节点nd!
14、选择所有被选择的单元esel,r,type,3*get,el,elem,0,num,mineknum=eknum+1!选取元素类型3!选取元素类型3中的最小编号，并存入变量elekel(eknum)=el*endif esel,all cmsel,s,t3 nd=ndnext(nd)*enddo PARSAVall,par1, anty,restPARRES,NEWpar1, *if,flag,eq,0,then cmsel,s,s1!将要杀死的单元存入数组ekel()!选择所有元素!选才I t3组节点!返回节点编号大于nd的节点编号，并存入变量nd!将参数与入文件par1(包括标量和矢量)!重启求解!从文件par1读取参数!如果 flag=0!选择节点s1SFGRAD,PRES,0,X,0,-0.1*(1+j), !加载 sf,all,pres,-10*elseflag=0*endif*if,eknum,ne,0,then*do,i,1,eknum ekill,ekel(i)*enddo!如果 eknum w 0!杀死下点应力 sl人丁 s1fail的单兀*endifesel,all nsel,all time,1+0.1*j deltim,.01,.001,.1 auto,onsolve*enddo