*HEADING
THERMAL LOADING CONTOUR INTEGRAL CALCULATION, MESH 1, PLANE STRESS
TEMPERATURE DEPENDENT THERMAL EXPANSION
*NODE
*NGEN,NSET=TIP
*NGEN,NSET=NSET1
*NGEN,NSET=NSET2,LINE=C
*NSET,NSET=NSET3
NSET1,NSET2
*NGEN,NSET=TOP
*NFILL,SINGULAR=1
TIP,NSET3,32,2
*NFILL,BIAS=.9
NSET2,TOP,16,4
*NSET,NSET=X0,GENERATE
*NSET,NSET=CLINE,GENERATE
*NSET,NSET=X2,GENERATE
*ELEMENT,TYPE=CPS8R
*ELGEN,ELSET=REGION1
*ELGEN,ELSET=REGION2
*ELSET,ELSET=ALL
REGION1,REGION2
*ELSET,ELSET=CRACK,GENERATE
*ELSET,ELSET=TOP,GENERATE
*NSET,NSET=INDEP1,GENERATE
*NSET,NSET=INDEP2,GENERATE
*NSET,NSET=INDEP3,GENERATE
*NSET,NSET=DEP1,GENERATE
*NSET,NSET=DEP2,GENERATE
*MPC
QUADRATIC,DEP1,INDEP1,INDEP2,INDEP3
QUADRATIC,DEP2,INDEP1,INDEP2,INDEP3
*SOLID SECTION,MATERIAL=MAT1,ELSET=ALL
1.,
*MATERIAL,NAME=MAT1
*ELASTIC
30.E6,.3
*EXPANSION,ZERO=0.
7.5E-6, 0.
12.E-6, 1000.
**
** STEP INFORMATION
**
*STEP
*STATIC
*TEMPERATURE,FILE=jintegraltherm_heatmesh1
*BOUNDARY
CLINE,2,2
TOP,2,2
1729,1
*CONTOUR INTEGRAL,CONTOURS=6,SYMM,OUTPUT=BOTH
TIP,-1.,0.
*CONTOUR INTEGRAL,CONTOURS=6,SYMM,OUTPUT=BOTH,TYPE=K FACTORS
TIP,-1.,0.
*CONTOUR INTEGRAL,CONTOURS=6,SYMM,OUTPUT=BOTH,TYPE=T-STRESS
TIP,-1.,0.
*EL PRINT,ELSET=CRACK
S,MISES,PRESS
E,
*ENERGY PRINT
*END STEP