*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
*NGEN,NSET=END1
*NGEN,NSET=BND1
*NGEN,NSET=END2
*NGEN,NSET=BND2
*NFILL
END1,END2,6,1
BND1,BND2,6,1
*NGEN,NSET=TIP
*NGEN,NSET=OUTER
*NFILL,SINGULAR=1
TIP,OUTER,6,1
*ELEMENT,TYPE=CPE8
**ELEMENT,TYPE=CPS8
*ELGEN,ELSET=EALL
*NSET,NSET=YAXIS
*NSET,NSET=XAXIS
*NSET,NSET=CRACK,GENERATE
*ELSET,ELSET=UP,GEN
*ELSET,ELSET=LOW,GEN
*ELSET,ELSET=UPLOAD,GEN
*ELSET,ELSET=LOWLOAD,GEN
*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,OUTPUT=BOTH
CRACK, -1.0, 0.0
*CONTOUR INTEGRAL,CONTOURS=3,OUTPUT=BOTH,TYPE=K FACTORS
CRACK, -1.0, 0.0
*EL PRINT
S,E
*EL FILE
S,E
*END STEP