*HEADING
 D.E.N  J - INTEGRAL, SUBMODEL ANALYSIS FROM jintegral2d_cpe8 
**RESTART,WRITE,FREQUENCY=1
*SYSTEM
5.0, 
*NODE
1,
21,
801,-5.0,
807,-5.0,5.0
821,5.0,
815,5.0,5.0
*NGEN,NSET=CRACK
1,21,1
*NGEN,NSET=TOP
801,807
807,815
815,821
*NFILL,SINGULAR=1
CRACK,TOP,8,100
*NSET,NSET=XAXIS,GENERATE
1,701,100
*NSET,NSET=DRIVEN,GENERATE
801,821
*ELEMENT,TYPE=CPE8
1,1,3,203,201,2,103,202,101
*ELGEN,ELSET=ALL
1,10,2,1,4,200,11
*SOLID SECTION,MATERIAL=A1,ELSET=ALL
1., 
*MATERIAL,NAME=A1
*ELASTIC
30.E6,.3
*BOUNDARY
XAXIS,2
*SUBMODEL,EXTERIOR TOLERANCE=0.025
DRIVEN,
*STEP
*STATIC
*BOUNDARY,SUBMODEL,STEP=1
DRIVEN,1,2
*CONTOUR INTEGRAL,FREQUENCY=1,CONTOURS=4,SYMM,OUTPUT=BOTH
CRACK, -1.0, 0.0
*CONTOUR INTEGRAL,FREQUENCY=1,CONTOURS=4,SYMM,OUTPUT=BOTH,TYPE=K FACTORS
CRACK, -1.0, 0.0
*CONTOUR INTEGRAL,FREQUENCY=1,CONTOURS=4,SYMM,OUTPUT=BOTH,TYPE=T-STRESS
CRACK, -1.0, 0.0
*ENERGY PRINT
*NODE PRINT
U,RF
*EL PRINT
S,E
*NODE FILE,NSET=DRIVEN
U, 
*OUTPUT,FIELD
*NODE OUTPUT,NSET=DRIVEN
U,
*END STEP