*HEADING
 CONTOUR INTEGRALS 
 SINGLE EDGE NOTCH SPECIMEN: a/W=0.5, H/W=2
 File name: cintegral2d_1edge_cpe8.inp
*NODE
1, 
7, 10.
87,10., 10.
247, -10., 10.
327, -10.
93, 10., 40.
253, -10.,  40.
321, 
*NGEN,NSET=END1
87,247,20
*NGEN,NSET=END2
93,253,20
*NFILL
END1,END2,6,1
*NGEN,NSET=TIP
1,321,20
*NGEN,NSET=OUTER
7,87,20
87,247,20
247,327,20
*NFILL,SINGULAR=1
TIP,OUTER,6,1
*ELEMENT,TYPE=CPS8
1,  1,3,43,41,  2,23,42,21
*ELGEN,ELSET=EALL
1,3,2,1,8,40,10
23,4,2,1,4,40,10
*NSET,NSET=XAXIS,GENERATE
1,7
*NSET,NSET=YAXIS
7, 
*ELSET,ELSET=PULL
26,36,46,56
*NSET,NSET=CRACK,GENERATE
1,321,20
*SOLID SECTION,MATERIAL=A1,ELSET=EALL
1., 
*MATERIAL,NAME=A1
*ELASTIC
30.E6,.3
*BOUNDARY
XAXIS,2
YAXIS,1
*STEP
*STATIC
*DLOAD
PULL,P2,-1.
*CONTOUR INTEGRAL,FREQUENCY=1,CONTOURS=3,SYMM,OUTPUT=BOTH,TYPE=J
CRACK, 1.0, 0.0
*CONTOUR INTEGRAL,FREQUENCY=1,CONTOURS=3,SYMM,OUTPUT=BOTH,TYPE=K FACTORS
CRACK, 1.0, 0.0
*CONTOUR INTEGRAL,FREQUENCY=1,CONTOURS=3,SYMM,OUTPUT=BOTH,TYPE=T-STRESS
CRACK, 1.0, 0.0
*NODE PRINT,freq=0
U,RF
*EL PRINT,freq=0
S,E
*NODE FILE,freq=0
U,
EL FILE,freq=0
S,
*END STEP