*HEADING
THERMAL LOADING CONTOUR INTEGRAL CALCULATION, MESH 1, PLANE STRESS
TEMPERATURE DEPENDENT THERMAL EXPANSION
*NODE
1,10.
6401,10.
65,
1665,,6.666
4865,20.,6.666
6465,20.
9999,10.,-10.
1729,,40.
4929,20.,40.
*NGEN,NSET=TIP
1,6401,200
*NGEN,NSET=NSET1
65,1665,200
4865,6465,200
*NGEN,NSET=NSET2,LINE=C
1665,4865,200,9999
*NSET,NSET=NSET3
NSET1,NSET2
*NGEN,NSET=TOP
1729,4929,200
*NFILL,SINGULAR=1
TIP,NSET3,32,2
*NFILL,BIAS=.9
NSET2,TOP,16,4
*NSET,NSET=X0,GENERATE
1665,1729,4
65,1665,200
*NSET,NSET=CLINE,GENERATE
1,65,2
*NSET,NSET=X2,GENERATE
4865,6465,200
4865,4929,8
*ELEMENT,TYPE=CPS8R
1,1,401,409,9,201,405,209,5
1665,1665,2465,2481,1681,2065,2473,2081,1673
*ELGEN,ELSET=REGION1
1,8,8,8,16,400,400
*ELGEN,ELSET=REGION2
1665,4,16,16,4,800,800
*ELSET,ELSET=ALL
REGION1,REGION2
*ELSET,ELSET=CRACK,GENERATE
1,6401,400
9,6409,400
17,6417,400
*ELSET,ELSET=TOP,GENERATE
1697,4097,800
*NSET,NSET=INDEP1,GENERATE
1665,4065,800
*NSET,NSET=INDEP2,GENERATE
2065,4465,800
*NSET,NSET=INDEP3,GENERATE
2465,4865,800
*NSET,NSET=DEP1,GENERATE
1865,4265,800
*NSET,NSET=DEP2,GENERATE
2265,4665,800
*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