*HEADING
FULLY PLASTIC J-INTEGRAL - PLANE STRESS - COARSE MESH 
*SYSTEM
5,0
*NODE
1, 
7,10
87,10,10
247,-10,10
327,-10
91,10,30
251,-10,30
321, 
1891,10,100
2051,-10,100
*NGEN,NSET=OUTER
7,87,20
87,247,20
247,327,20
*NGEN,NSET=CRACK
1,321,20
*NFILL,SINGULAR=1, NSET=NALL
CRACK,OUTER,6,1
*NGEN,NSET=BOT
87,247,20
*NGEN,NSET=MID
91,251,20
*NGEN,NSET=TOP
1891,2051,20
*NFILL, NSET=NALL
BOT,MID,4,1
*NFILL, NSET=NALL
MID,TOP,6,300
*ELEMENT,TYPE=CPS8R
1,1,3,43,41,2,23,42,21
125,91,691,731,131,391,711,431,111
*ELGEN,ELSET=ALL
1,3,2,1,8,40,10
23,3,2,1,4,40,10
125,3,600,100,4,40,10
*NSET,NSET=XAXIS,GENERATE
321,327
*NSET,NSET=YAXIS,GENERATE
247,251
247,327,20
251,2051,300
*ELSET,ELSET=PRINT
1,41
*ELSET,ELSET=PULL
325,335,345,355
*ELSET,ELSET=DRAW,GENERATE
1,71,10
2,72,10
3,73,10
*ELSET,ELSET=MONITOR
DRAW,24,34,44,54
*SOLID SECTION,MATERIAL=A1,ELSET=ALL
*MATERIAL,NAME=A1
*DEFORMATION PLASTICITY
30E6,0.3,40E3,5,.5
*BOUNDARY
XAXIS,2
YAXIS,1
**
*STEP,INC=20
*STATIC,FULLY PLASTIC=MONITOR
.05,1
*DLOAD
PULL,P2,-100.E4
*CONTOUR INTEGRAL,FREQUENCY=1,CONTOURS=3,SYMM,OUTPUT=BOTH
CRACK, -1.0, 0.0
*ENERGY PRINT
*ENERGY FILE
*EL PRINT,ELSET=PRINT
S, 
E, 
*EL FILE,ELSET=PRINT
S, 
E, 
*NODE PRINT
U, 
*NODE FILE
U, 
*OUTPUT,HISTORY
*ENERGY OUTPUT,VARIABLE=ALL
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=PRINT
S, 
E, 
*OUTPUT,HISTORY
*ELEMENT OUTPUT,ELSET=PRINT
S, 
E, 
*OUTPUT,FIELD
*NODE OUTPUT
U, 
*OUTPUT,HISTORY
*NODE OUTPUT, NSET=NALL
U, 
*END STEP