I wrote a batch file for Crack growth model for ANSYS but when the soution is running,ANSYS terminate itself...Why?

Hello everybody,

I am searching precision of crack propogation in a bonded wafer double.I wrote a code but the problem is that when I want to solve it in ANSYS, the program is shutdown. I need help about it. Thanks for your helps...

 Greetings

 Volkan DALMISLI

P.S:The code is given below. 

finish

/title,sphere wafer

/filname,sphere wafer

/clear,start

*CFOPEN,volkan,dat           ! deletes the data file

*CFWRITE,releaserate     ! table heading of the data file

*CFCLOS

/PNUM, volu, 1 

/prep7

!***********************************Modell generation**************************

k,9,0,-50e2,-50   ! Creating keypoints

k,19,0,-50e2,0

k,29,0,-50e2,50

k,11,0,50e2,-50

k,21,0,50e2,0

k,31,0,50e2,50

k,10,-4999,0,-70

k,20,-4999,0,-20

k,22,-4999,0,20

k,30,-4999,0,70

k,55,4999,0,-50

k,56,4999,0,0

k,57,4999,0,50

k,6,100,-60e2,-100

!k,16,100,-50e2,0

k,26,100,-60e2,100

k,5,100,60e2,-100

!k,15,100,50e2,0

k,25,100,60e2,100

!******Creating Volumes at delamination region******

larc,31,29,30   !Drawing arc lines

larc,21,19,22

askin,1,2   !Creating Areas by skining

l,29,31    !Drawing lines

l,19,21

l,9,11

al,3,5,4,6   !Creating Areas by lines

al,5,1    

al,6,2

va,1,2,3,4   !Creating top volume

   

larc,21,19,20    

larc,9,11,10

askin,8,9

al,6,10,7,11

al,9,7

al,8,6

VA,5,6,7,8   !Creating bottom volume

!*********************************

!****Creating minor flat******

k,100,-4918,10e2,100  !Creating area for dividing minor flat

k,101,-4918,-10e2,100

k,102,-4918,10e2,-100

k,103,-4918,-10e2,-100

a,101,103,102,100

vsba,1,9,,,keep

vsba,2,9,,,

vdele,1,,,1

vdele,3,,,1

!*********************************

!********Creating crack front refinement area******

k,100,-100,60e2,100

k,101,-100,60e2,-100

k,102,-100,-60e2,-100

k,103,-100,-60e2,100

a,100,101,102,103 !Creating area for refine mesh around crack front (at delamiation region)

vsba,4,1,,,keep

vsba,5,1,,,

!*********************************

!******Creating rest area of the wafer*****

blc4,0,-50e2,50e2,100e2,50 !I need to create a block volume for brick meshing

blc4,0,-50e2,50e2,100e2,-50 !because Ansys gives error for other geometries

!*********************************

!********Creating crack front refinement area******   

k,100,100,60e2,100 !Creating area for refine mesh around crack front (at bonded region)

k,101,100,60e2,-100

k,102,100,-60e2,-100

k,103,100,-60e2,100

a,100,101,102,103

vsba,5,35,,,keep

vsba,6,35,,,

!*********************************

vglue,5,7,8,9  !Gluing bondend volumes

!******************************************************************************

!***********************************element type*******************************

et,1,95              ! Elementtype: Solid95

!***********************************material properties************************

!***********************************Silicone***********************************

mp,ex,1,168919              ! Young's Modulus in [110] direction

                 ! 168,919 GPa, isotropic

mp,nuxy,1,0.064189             ! possions ratio

                  ! 0,064189, isotropic

 

!***********************************meshing************************************

vsel,all

vsel,s,,,2

vsel,a,,,4

vsel,a,,,10

vsel,a,,,11

ESIZE,500

mshape,0,3d

MSHKEY,1 ! Mapped Mesh

VMESH,ALL ! Meshing des Volumens

allsel

vsel,s,,,1

vsel,a,,,3

vsel,a,,,5

vsel,a,,,6

ESIZE,50

mshape,0,3d

MSHKEY,1 ! Mapped Mesh

VMESH,ALL ! Meshing des Volumens

finish

/solu

antype,static

da,28,all !Fixing back side of the wafer to a wall etc...

da,52,all !Fixing back side of the wafer to a wall etc...

nsel,all

nsel,s,loc,x,-4750,-4700 !At the top volume: Selecting a node line to present touching points of the razor blade

nsel,r,loc,y,50e2,-50e2

nsel,r,loc,z,0,22

f,all,fZ,1e-5   !Applying force to selecting nodes to present blade insertion

nsel,s,loc,x,-4750,-4700 !At the bottom volume: Selecting a node line to present touching points of the razor blade

nsel,r,loc,y,50e2,-50e2

nsel,r,loc,z,-22,0

f,all,fZ,-1e-5   !Applying force to selecting nodes to present blade insertion

solve