/prep7
dis1=0.9
dis2=12.0
n1=1000
n2=1000
n3=10
dl=100
dh=3
a0=30
nel=200
neh=6
toler=0.1e-5
et,1,182 ! 2d 4-node structural solid element
keyopt,1,1,2 ! enhance strain formulation
keyopt,1,3,2 !* plane strain
et,2,182
keyopt,2,1,2
keyopt,2,3,2
et,3,202 ! 2d 4-node cohesive zone element
!keyopt,3,2,2 ! element free option
keyopt,3,3,2 !* plane strain
mp,ex,1,1.353e5 ! e11 = 135.3 gpa
mp,ey,1,9.0e3 ! e22 = 9.0 gpa
mp,ez,1,9.0e3 ! e33 = 9.0 gpa
mp,gxy,1,5.2e3 ! g12 = 5.2 gpa
mp,prxy,1,0.24
mp,prxz,1,0.24
mp,pryz,1,0.46
g1c=0.28 ! critical energy-release rate
g2c=0.80
g3c=0.80
tb,cgcr,1,,3,linear ! linear fracture criterion
tbdata,1,g1c,g2c,g3c
! fe model
rectng,0,dl,dh/2 ! define areas
rectng,0,dl,0,-dh/2
lsel,s,line,,2,8,2 ! define line division
lesize,all,dh/neh
lsel,inve
lesize,all, , ,nel
allsel,all
type,1 ! mesh area 2
mat,1
local,11,0,0,0,0
esys,11
amesh,2
csys,0
type,2 ! mesh area 1
esys,11
amesh,1
csys,0
nsel,s,loc,x,a0-toler,dl
nummrg,nodes
esln
type,3
mat,5
czmesh,,,1,y,0, ! generate interface elements
allsel,all
nsel,s,loc,x,dl ! apply constraints
d,all,all
nsel,all
!
esel,s,ename,,202 ! select interface element to
cm,cpath,elem ! define crack growth path
Release 15.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential informationVCCT-Based Crack Growth Simulation