二维J积分求解实例命令流

本人做的二维裂缝断裂力学求解实例
BLC4,0,0,0.05,0.075
wpoff,0.01
wprot,,,90
ASBW, 1
wpoff,,0.01
wprot,,90
asbw,all
wpoff,,0.01
wprot,,90
asbw,all
et,1,2
MP,EX,1,2.1e11
MP,PRXY,1,0.3
*CREATE,JIN1,MAC ！创建宏文件，求解J积分
STINFC ！定义数据块名称
SEXP,W,SENE,VOLU,1,-1 ！计算应变能密度
PATH,JINT,4,50,48 ！定义路径
PPATH,1,ARG1 ！定义路径点
PPATH,2,ARG2
PPATH,3,ARG3
PPATH,4,ARG4
PDEF,W,ETAB,W ！将应变能密度映射到路径上
PCALC,INTG,J,W,YG ！应变能密度计算
*GET,JA,PATH,,LAST,J
PDEF,CLEAR ！删除路径变量
PVECT,NORM,NX,NY,NZ ！定义路径单位向量
PDEF,INTR,SX,SX ！将X轴应力映射到路径上
PDEF,INTR,SY,SY ！将Y轴应力映射到路径上
PDEF,INTR,SXY,SXY ！将Z轴应力映射到路径上
PCALC,MULT,TX,SX,NX ！路径相乘操作
PCALC,MULT,C1,SXY,NX
PCALC,ADD,TX,TX,C1 ！路径相加操作
PCALC,MULT,TY,SXY,NX
PCALC,MULT,C1,SY,NY
PCALC,ADD,TY,TY,C1
*GET,DX,PATH,,LAST,S
DX=DX/100
PCALC,ADD,XG,XG,,,,-DX/2
PDEF,INTR,UX1,UX
PDEF,INTR,UY1,UY
PCALC,ADD,XG,XG,,,,DX
PDEF,INTR,UX2,UX
PDEF,INTR,UY2,UY
PCALC,ADD,XG,XG,,,,-DX/2
C=(1/DX)
PCALC,ADD,C1,UX2,UX1,C,-C
PCALC,ADD,C2,UY2,UY1,C,-C
PCALC,MULT,C1,TX,C1
PCALC,MULT,C2,TY,C2
PCALC,ADD,C1,C1,C2
PCALC,INTG,J,C1,S
*GET,JB,PATH,,LAST,J
JINT=2*(JA-JB)
PDEF,CLEAR ！删除路径变量
*END ！结束宏文件
FLST,5,6,5,ORDE,4
FITEM,5,1
FITEM,5,-3
FITEM,5,5
FITEM,5,-7
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT, 1, , 1, 0,
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
KSCON,5,0.001,1,12,0.5,
MSHAPE,0,2D
MSHKEY,0
!*
FLST,5,2,4,ORDE,2
FITEM,5,14
FITEM,5,17
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,10, , , , ,1
!*
FLST,5,4,4,ORDE,3
FITEM,5,10
FITEM,5,-12
FITEM,5,15
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,25, , , , ,1
!*
FLST,5,3,4,ORDE,3
FITEM,5,2
FITEM,5,4
FITEM,5,-5
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,10, , , , ,1
!*
FLST,5,3,4,ORDE,3
FITEM,5,1
FITEM,5,4
FITEM,5,13
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,10, , , , ,1
!*
FLST,5,3,4,ORDE,3
FITEM,5,3
FITEM,5,8
FITEM,5,16
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,10,0.2, , , ,1
!*
FLST,5,2,5,ORDE,2
FITEM,5,3
FITEM,5,5
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
FLST,5,4,5,ORDE,4
FITEM,5,1
FITEM,5,-2
FITEM,5,6
FITEM,5,-7
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
WPCSYS,-1,0
nsel,s,loc,x,0
nplot
dsym,symm,x
nsel,s,loc,x,0.01,0.05
nsel,r,loc,y,0
dsym,symm,y
allsel,all
nsel,s,loc,y,0.075
sf,all,pres,-140e6
allsel,all
wpoff,0.01
CSWPLA,11,0,1,1,
/solu
antype,static
outpr,all,all
slove
finish
/post1
CSYS,11,
FLST,2,3,1
FITEM,2,1
FITEM,2,589
FITEM,2,588
!*
PATH,1,3,30,20,
PPATH,P51X,1
PATH,STAT
KCALC,1,1,0,0
plnsol,u,x
plnsol,u,y
plnsol,u,sum
plnsol,s,eqv
etable,sene,sene ！存储应变能密度到单元列表中
etable,volu,volu ！存储体积到单元列表中
*SET,node20,25
*SET,node21,550
*SET,node22,1228
*SET,node23,573
*ulib,jin1,mac
*use,stinfc,node20,node21,node22,node23 ！调用宏文件
*SET,con1,2.1e11
*SET,ki2,sqrt(con1*jint))
*get,jint
*status,ki1 ！显示KI1
*status,ki2 ！显示KI2