ANSYS crack propagation(growth) under cycle load simulation.

I want to simulate the crack growth under cycle load through ANSYS

I have thought two way to simulate it, first way is cohesive model(debonding)

second way is remesh(rezone at certain load step), this way is important for me, since this method mapsolves the new mesh into the old mesh, so the plastic strain information near crack tip element could transfer to new mesh.

I need a example file(code-command file of ansys) or comment

sanghuyck [at] naver.com (sanghuyck[at]naver[dot]com) please send me e-mail

below text is my simulation model

FINISH 

/CLEAR,START

/COM,ANSYS RELEASE 12.0.1  UP20090415       10:10:38    03/21/2011

/input,start120,ans,'C:\Program Files\ANSYS Inc\v120\ANSYS\apdl\',,,,,,,,,,,,,,,,1 

/PREP7 

ET,1,PLANE182 

ET,2,TARGE169          ! 2-D TARGET ELEMENTS

ET,3,CONTA175          ! 2-D CONTACT ELEMENTS

!7075-T6 property (Mpa)

MP,EX,1,40122.1146 

MP,PRXY,1,.33  

TB,kinh,1,1,10,0

TBTEMP,0

TBPT,,0.01343, 538.843

TBPT,,0.01852, 581.8182

TBPT,,0.02778, 608.2645

TBPT,,0.04444, 647.9339

TBPT,,0.06435, 690.9091

TBPT,,0.07778, 717.3553

TBPT,,0.08935, 737.1901

TBPT,,0.09491, 747.1074

!specimen geometry (mm)

CYL4,158.5,150,5,180

CYL4,158.5,150,5,-180

cyl4,146.5,150,5,180

cyl4,146.5,150,5,-180

rect,0,152.5*2,0,150

rect,0,152.5*2,150,300

ASBA,6,3,,DELETE,KEEP

ASBA,7,1,,DELETE,KEEP

ASBA,5,4,,DELETE,KEEP

ASBA,7,2,,DELETE,KEEP

!kp murge

ksel,s,loc,x,158.5

ksel,r,loc,y,150

nummrg,all

ksel,s,loc,x,146.5

ksel,r,loc,y,150

nummrg,all

ksel,s,loc,x,146.5

ksel,r,loc,y,150

cm,y,kp

KSCON,y,.1,0,30

ksel,s,loc,x,158.5

ksel,r,loc,y,150

cm,y,kp

KSCON,y,.1,0,30

!mesh

amesh,1

amesh,2

amesh,3

amesh,4

esize,10

amesh,6

amesh,5

!murge of node

nsel,s,loc,y,150

nsel,r,loc,x,0,146.5

nummrg,all

nsel,s,loc,y,150

nsel,r,loc,x,158.5,152.5*2

nummrg,all

!contact element

lsel,s,line,,9

lsel,a,line,,23

lsel,a,line,,6

nsll,,1

REAL,1

TYPE,2

ESURF 

lsel,s,line,,11

lsel,a,line,,15

lsel,a,line,,2

nsll,,1

REAL,1

TYPE,3

ESURF

!constrain

nsel,s,loc,y,300

d,all,all

!load control

endtime=4

*dim,loadsign,arrary,endtime

in=1

*Do,i,1,endtime,1

   *if,i,eq,in,THEN

       loadsign(i)=-4

       in=i+2

   *else

      loadsign(i)=-1

   *endif

*enddo

 /solu

   nsel,s,loc,y

    sf,all,pres,0*loadsign(1)

    time,1

    deltim,1,0.1,0.5

    solve

!solve along load cycle

*do,i,2,endtime,1

  /solu

 

     nsel,s,loc,y

    sf,all,pres,80*loadsign(i)

 

 time,i

 nsub,8,7,3

 solve

*enddo

/post1

!failure criteria

FC,1,TEMP,,

FC,1,S,XTEN,750

FC,1,S,XCMP,-750,

FC,1,S,YTEN,750

FC,1,S,YCMP,-750,

FC,1,S,ZTEN,750

FC,1,S,ZCMP,-750,

FC,1,S,XY,750

FC,1,S,YZ,750,

FC,1,S,XZ,750,

FC,1,S,XYCP,

FC,1,S,YZCP,

FC,1,S,XZCP,