ansys断裂例子.docx
《ansys断裂例子.docx》由会员分享，可在线阅读，更多相关《ansys断裂例子.docx（10页珍藏版）》请在冰点文库上搜索。
ansys断裂例子
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,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,,70MPDATA,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,xd,all,ux
nsel,s,loc,yd,all,uy
d,all,ux
mp,mu,3,0
/COM,CONTACTPAIRCREATION-START
CM,_NODECM,NODECM,_ELEMCM,ELEMCM,_LINECM,LINECM,_AREACM,AREA/GSAV,cwz,gsav,,temp
MP,MU,3,0
MAT,3
R,3
REAL,3
ET,2,169
ET,3,172
R,3,,,100,0.1,0,RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,0.5,1.0,0.0,
KEYOPT,3,2,0
KEYOPT,3,3,0
KEYOPT,3,4,0
KEYOPT,3,5,0
KEYOPT,3,7,0
KEYOPT,3,8,0
KEYOPT,3,9,0
KEYOPT,3,10,0KEYOPT,3,11,0KEYOPT,3,12,5!
Generatethetargetsurface
NSEL,S,,,T1
CM,_TARGET,NODETYPE,2
ESLN,S,0
ESURF,ALL
CMSEL,S,_ELEMCM!
Generatethecontactsurface
NSEL,S,,,C1
CM,_CONTACT,NODETYPE,3
ESLN,S,0
ESURF,ALL
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3ESEL,R,REAL,,3/PSYMB,ESYS,1/PNUM,TYPE,1/NUM,1
EPLOT
ESEL,ALL
ESEL,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-END
lsel,s,,,7
nsll,s,1
cm,s1,node
!
Gradientsurfaceload
SFGRAD,PRES,0,X,0,-0.1,
sf,all,pres,-0.1nsel,all
esel,all
!
save
/solu
time,1
deltim,1,1,1solve
/post1
plns,s,1
/solu
anty,,rest
time,1.1
ekill,140
solve
/post1
plns,s,1
/solu
anty,,resttime,1.2
ekill,140ekill,139solve
/post1
plns,s,1
/solu
anty,,resttime,1.3
ekill,140ekill,139ekill,138solve
/post1
plns,s,1
/solu
anty,,resttime,1.4
ekill,140ekill,139ekill,138ekill,137solve
/post1
plns,s,1
第二版（加了断裂自动判断准则）。
整整花了半天时间才调试明白。
fini
/clear
/filn,crack1s1fail=100～定义断裂准则:
断裂应力=100
eknum=0
flag=0
/PREP7
!
*
ET,1,PLANE182!
*
KEYOPT,1,1,2KEYOPT,1,3,1
KEYOPT,1,4,0
KEYOPT,1,6,0
KEYOPT,1,10,0
!
*
rect,0,100,0,100rect,0,100,100,110lesi,1,,,10
lesi,2,,,10
esha,2～设置单元形状4边形!
*
MPTEMP,,,,,,,,～定义材料1MPTEMP,1,0
MPDATA,EX,1,,210e3MPDATA,PRXY,1,,0.3MPTEMP,,,,,,,,!
定义材料2MPTEMP,1,0
MPDATA,EX,2,,70MPDATA,PRXY,2,,0.33amesh,1
lesi,5,,,10
lesi,6,,,2
mat,2
amesh,2
lsel,s,,,3!
选择线3
nsll,s,1!
选择线3上的所有节点cm,c1,node!
定义已选择的点为c1组lsel,s,,,5!
选择线5
nsll,s,1!
选择线5上的所有节点cm,t1,node!
定义已选择的点为t1组nsel,s,loc,x!
选择X=0的所有节点d,all,ux!
nsel,s,loc,y!
选择Y=0的所有节点d,all,uy!
d,all,ux
mp,mu,3,0!
定义材料3
/COM,CONTACTPAIRCREATION-START!
定义接触对,采用接触对向导
CM,_NODECM,NODECM,_ELEMCM,ELEMCM,_LINECM,LINECM,_AREACM,AREA/GSAV,cwz,gsav,,tempMP,MU,3,0
MAT,3
R,3
REAL,3
ET,2,169
ET,3,172
R,3,,,100,0.1,0,!
Normalpenaltystiffness=100
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,0.5,1.0,0.0,
KEYOPT,3,2,1
KEYOPT,3,3,0
KEYOPT,3,4,0
KEYOPT,3,5,0
KEYOPT,3,7,2
KEYOPT,3,8,0
KEYOPT,3,9,0
KEYOPT,3,10,0KEYOPT,3,11,0KEYOPT,3,12,5!
Generatethetargetsurface
NSEL,S,,,T1!
目标面T1组节点CM,_TARGET,NODETYPE,2
ESLN,S,0
ESURF,ALL
CMSEL,S,_ELEMCM!
Generatethecontactsurface
NSEL,S,,,C1!
接触面C1组节点CM,_CONTACT,NODETYPE,3
ESLN,S,0
ESURF,ALL
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,2ESEL,A,TYPE,,3ESEL,R,REAL,,3/PSYMB,ESYS,1/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,2ESEL,A,TYPE,,3ESEL,R,REAL,,3CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
/COM,CONTACTPAIRCREATION-END
lsel,s,,,7!
选择线7
nsll,s,1!
选择线7上的节点
cm,s1,node!
定义节点组S1
!
GradientsurfaceloadSFGRAD,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!
选择元素类型3*get,ektnum,elem,0,count!
变量ektnum=元素总素*dim,ekel,,ektnum!
定义一维数组ekelesel,all!
选择所有元素
loadnum=8!
loadnum=[1,10]*do,j,1,loadnum
/post1
set,last!
读取载荷步结果/solu
esel,s,type,,3!
选择元素类型3nslensel,r,loc,y,100!
选择元素类型3,y=100的节点cm,t3,node!
定义节点组t3esel,all!
选择所有元素*get,ndnum,node,0,count!
变量ndnum=节点总数（t3组）*get,nd,node,0,num,min!
变量nd=节点的最小编号*do,i,1,ndnum*GET,s1nd,NODE,nd,S,1!
变量s1nd=节点nd的应力s1*if,s1nd,ge,s1fail,then!
断裂判据
flag=1
nsel,s,,,nd!
选择节点ndesln!
选择所有被选择的单元esel,r,type,,3!
选取元素类型3*get,el,elem,0,num,min!
选取元素类型3中的最小编号,并存入变量el
eknum=eknum+1
ekel（eknum）=el!
将要杀死的单元存入数组ekel（）
*endif
esel,all!
选择所有元素cmsel,s,t3!
选择t3组节点nd=ndnext（nd）!
返回节点编号大于nd的节点编号,并存入变量nd
*enddo
PARSAV,all,'par1','',''!
将参数写入文件par1（包括标量和矢量）
anty,,rest!
重启求解
PARRES,NEW,'par1','',''!
从文件par1读取参数*if,flag,eq,0,then!
如果flag=0cmsel,s,s1!
选择节点s1SFGRAD,PRES,0,X,0,-0.1*（1+j）,!
加载
sf,all,pres,-10
*else
flag=0
*endif
*if,eknum,ne,0,then!
如果eknum?
0*do,i,1,eknum
ekill,ekel（i）!
杀死节点应力sl大于s1fail的单元
*enddo
*endif
esel,all
nsel,all
time,1+0.1*j
deltim,.01,.001,.1
auto,on
solve
*enddo
升级会员 