*HEADING
THERMAL LOADING CONTOUR INTEGRAL CALC., MESH 2, 3-D STRESS ELEMENTS
*NODE
**
**
*NSET,NSET=RIGHT,GENERATE
*NSET,NSET=YAXIS,GENERATE
*NGEN,NSET=CRFR
*NGEN,NSET=CRMD
*NGEN,NSET=CRBK
*NSET,NSET=CRACK
CRFR,CRMD,CRBK
*NGEN,NSET=OUTER
*NSET,NSET=LOW,GENERATE
*NGEN,NSET=TOP
*NFILL,BIAS=.9
LOW,TOP,14,20
*NFILL,SINGULAR=1
CRACK,OUTER,12,30
*NSET,NSET=FRONT,GENERATE
*NSET,NSET=BACK,GENERATE
*ELEMENT,TYPE=C3D20R
***
*ELGEN,ELSET=ELS
*ELSET,ELSET=PRINT,GENERATE
*NSET,NSET=XAXIS,GENERATE
*SOLID SECTION,MATERIAL=MAT1,ELSET=ELS
*MATERIAL,NAME=MAT1
*ELASTIC
30.E6,.3
*EXPANSION,ZERO=0.
7.5E-6,
*NSET,NSET=FRONT,GENERATE
*NSET,NSET=BACK,GENERATE
*NSET,NSET=XAXIS,GENERATE
**
** STEP INFORMATION
**
*STEP
*STATIC
*TEMPERATURE,FILE=jintegraltherm_3dheat.odb
*BOUNDARY
*BOUNDARY
XAXIS,2
TOP,2
379,1
1379,1
2379,1
FRONT,3
BACK,3
*CONTOUR INTEGRAL,FREQUENCY=1,CONTOURS=6,SYMM,OUTPUT=BOTH,NORMAL
CRFR,CRMD,CRBK,
*CONTOUR INTEGRAL,FREQUENCY=1,CONTOURS=6,SYMM,OUTPUT=BOTH,NORMAL,TYPE=K FACTORS
CRFR,CRMD,CRBK,
*CONTOUR INTEGRAL,FREQUENCY=1,CONTOURS=6,SYMM,OUTPUT=BOTH,NORMAL,TYPE=T-STRESS
CRFR,CRMD,CRBK,
*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