*HEADING
THERMAL LOADING CONTOUR INTEGRAL CALCULATION, MESH 2 PLANE STRAIN
*NODE
1, 
25, 
361,10.,0.
367,10.,10.
379,-10.,10.
385,-10,0.
647,10.,40.
659,-10.,40.
*NSET,NSET=RIGHT,GENERATE
361,367
367,527,20
527,647,20
*NSET,NSET=YAXIS,GENERATE
379,385
379,539,20
539,659,20
*NGEN,NSET=CRACK
1,25
*NGEN,NSET=OUTER
361,367
367,379
379,385
*NFILL,SINGULAR=1
CRACK,OUTER,12,30
*NGEN,NSET=LOW
367,379
*NGEN,NSET=TOP
647,659
*NFILL,BIAS=.9
LOW,TOP,14,20
*ELEMENT,TYPE=CPE8R
1,1,61,63,3,31,62,33,2
73,367,407,409,369,387,408,389,368
*ELGEN,ELSET=ELS
1,6,60,1,12,2,6
73,6,2,1,7,40,6
*ELSET,ELSET=PRINT,GENERATE
1,6
*NSET,NSET=XAXIS,GENERATE
25,385,30
*SOLID SECTION,ELSET=ELS,MATERIAL=ELAS
*MATERIAL,NAME=ELAS
*ELASTIC
30.E6,.3
*EXPANSION,ZERO=0.
7.5E-6, 
*BOUNDARY
XAXIS,2
TOP,2
658,1
*STEP
*STATIC
*TEMPERATURE,FILE=jintegraltherm_heatmesh2
*CONTOUR INTEGRAL,FREQUENCY=1,CONTOURS=6,SYMM,OUTPUT=BOTH
CRACK, -1.0, 0.0
*CONTOUR INTEGRAL,FREQUENCY=1,CONTOURS=6,SYMM,OUTPUT=BOTH,TYPE=K FACTORS
CRACK, -1.0, 0.0
*CONTOUR INTEGRAL,FREQUENCY=1,CONTOURS=6,SYMM,OUTPUT=BOTH,TYPE=T-STRESS
CRACK, -1.0, 0.0
*EL PRINT,ELSET=PRINT
S,MISES,PRESS
E, 
*EL FILE,ELSET=PRINT
S,SINV
E, 
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=PRINT
S,SINV
E,
*ENERGY PRINT
*ENERGY FILE
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*END STEP