*HEADING
   SURFACE CRACK IN A RECTANGULAR PLATE UNDER UNIFORM TENSION
   H/t=16, W/t=8, a/t=0.6, a/c=0.24
   a:  MINOR CRACK DEPTH
   c:  MAJOR CRACK WIDTH 
   t:  PLATE THICKNESS
   H:  PLATE HEIGHT
   W:  PLATE WIDTH 
   File name : jktintegral3d.inp
**
**  INCLUDE NODE DATA
*INCLUDE,INPUT=jktintegral3d_node.inp
**
**  INCLUDE ELEMENT DATA
*INCLUDE,INPUT=jktintegral3d_element.inp
*SOLID SECTION,ELSET=EALL,MATERIAL=M
*MATERIAL,NAME=M
*ELASTIC
30.E6,.3
*BOUNDARY
BOT,2
SYM,1
XFIX,3
*ELSET,ELSET=RING,GEN
1,3200
*STEP
*STATIC
*DLOAD
TOP,P2,-1.0
*CONTOUR INTEGRAL,CONTOUR=5,TYPE=J,       SYMM,OUTPUT=BOTH,NORMAL,SURFACE NORMAL=FREE
0.0,-1.0,0.0
T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11,
T12,T13,T14,T15,T16,T17,T18,T19,T20,T21
*CONTOUR INTEGRAL,CONTOUR=5,TYPE=K FACTOR,SYMM,OUTPUT=BOTH,NORMAL,SURFACE NORMAL=FREE
0.0,-1.0,0.0
T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11,
T12,T13,T14,T15,T16,T17,T18,T19,T20,T21
*CONTOUR INTEGRAL,CONTOUR=5,TYPE=T-STRESS,SYMM,OUTPUT=BOTH,NORMAL,SURFACE NORMAL=FREE
0.0,-1.0,0.0
T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11,
T12,T13,T14,T15,T16,T17,T18,T19,T20,T21
*EL PRINT,ELSET=RING,freq=0
S,
E,
*NODE PRINT,freq=0
rf3,
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=TOP
E,S,
*ENDSTEP