*HEADING
THERMAL LOADING CONTOUR INTEGRAL CALCULATION, MESH 2 PLANE STRAIN
*NODE
*NSET,NSET=RIGHT,GENERATE
*NSET,NSET=YAXIS,GENERATE
*NGEN,NSET=CRACK
*NGEN,NSET=OUTER
*NFILL,SINGULAR=1
CRACK,OUTER,12,30
*NGEN,NSET=LOW
*NGEN,NSET=TOP
*NFILL,BIAS=.9
LOW,TOP,14,20
*ELEMENT,TYPE=CPE8R
*ELGEN,ELSET=ELS
*ELSET,ELSET=PRINT,GENERATE
*NSET,NSET=XAXIS,GENERATE
*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