*HEADING
 CONTOUR INTEGRALS FOR AN INTERFACE CRACK 
 SINGLE EDGE NOTCH SPECIMEN: a/W=0.5, H/W=2 with full plane meshes
 File name: cintegral2d_1edge_intf_cpe8.inp 
*NODE
1, 
7, 10.
87,10., 10.
247, -10., 10.
327, -10.
93, 10., 40.
253, -10.,  40.
641, 
407, -10., -10.
567,  10., -10.
413, -10., -40.
573,  10., -40.
647,  10.
*NGEN,NSET=END1
87,247,20
*NGEN,NSET=BND1
407,567,20
*NGEN,NSET=END2
93,253,20
*NGEN,NSET=BND2
413,573,20
*NFILL
END1,END2,6,1
BND1,BND2,6,1
*NGEN,NSET=TIP
1,641,20
*NGEN,NSET=OUTER
7,87,20
87,247,20
247,327,20
327,407,20
407,567,20
567,647,20
*NFILL,SINGULAR=1
TIP,OUTER,6,1
*ELEMENT,TYPE=CPE8
**ELEMENT,TYPE=CPS8
1,  1,3,43,41,  2,23,42,21
*ELGEN,ELSET=EALL
1,3,2,1,16,40,10
23,4,2,1,4,40,10
103,4,2,1,4,40,10
*NSET,NSET=YAXIS
327, 
*NSET,NSET=XAXIS
327, 
*NSET,NSET=CRACK,GENERATE
1,641,20
*ELSET,ELSET=UP,GEN
1,71,10
2,72,10
3,73,10
24,54,10
25,55,10
26,56,10
*ELSET,ELSET=LOW,GEN
81,151,10
82,152,10
83,153,10
104,134,10
105,135,10
106,136,10
*ELSET,ELSET=UPLOAD,GEN
26,56,10
*ELSET,ELSET=LOWLOAD,GEN
106,136,10
*SOLID SECTION,MATERIAL=A1,ELSET=UP
1., 
*SOLID SECTION,MATERIAL=A2,ELSET=LOW
*MATERIAL,NAME=A1
*ELASTIC
1.0E6,0.3
*MATERIAL,NAME=A2
*ELASTIC
2.0E6,0.2
*BOUNDARY
XAXIS,1,2
*STEP
*STATIC
*DLOAD
UPLOAD, P2,-1.0
LOWLOAD,P2,-1.0
**BOUNDARY
**END2,2,, 0.1
**BND2,2,,-0.1
*CONTOUR INTEGRAL,CONTOURS=3
CRACK, -1.0, 0.0
*CONTOUR INTEGRAL,CONTOURS=3,TYPE=K FACTORS
CRACK, -1.0, 0.0
*EL PRINT
S,E
**EL FILE
**S,E
*OUTPUT,FIELD
*ELEMENT OUTPUT
E,S
*END STEP