*HEADING
J-INTEGRAL FOR SEMI-ELLIPTICAL SURFACE CRACK
IN A SEMI-INFINITE BODY
  A/C=.4      A/T=.1
     A,  MINOR CRACK DEPTH
     T,  MINOR BODY THICKNESS
     C,  MAJOR CRACK WIDTH
*PREPRINT,HISTORY=NO,MODEL=NO
*RESTART,WRITE
*NODE,INPUT=JINT.NOD
*NSET,NSET=XYPL,GENERATE
*NSET,NSET=XYPL,GENERATE
*NSET,NSET=YZPL,GENERATE
*NSET,NSET=YZPL,GENERATE
*NSET,NSET=RAD0,GENERATE
*NSET,NSET=RAD1,GENERATE
*NSET,NSET=RAD2,GENERATE
*NSET,NSET=RAD3,GENERATE
*NSET,NSET=RAD4,GENERATE
*NSET,NSET=RAD5,GENERATE
*NSET,NSET=RAD6,GENERATE
*NSET,NSET=RAD7,GENERATE
*NSET,NSET=RAD8,GENERATE
*NSET,NSET=RAD9,GENERATE
*NSET,NSET=RADA,GENERATE
*NSET,NSET=RADB,GENERATE
*NSET,NSET=RADC,GENERATE
*NSET,NSET=RADD,GENERATE
*NSET,NSET=RADE,GENERATE
*NSET,NSET=RADF,GENERATE
*NSET,NSET=RADG,GENERATE
*NSET,NSET=LEV0,GENERATE
*NSET,NSET=LEV1,GENERATE
*NSET,NSET=LEV2,GENERATE
*NSET,NSET=LEV3,GENERATE
*NSET,NSET=LEV4,GENERATE
*NSET,NSET=LEV5,GENERATE
*NSET,NSET=LEV6,GENERATE
*NSET,NSET=LEV7,GENERATE
*NSET,NSET=LEV8,GENERATE
*NSET,NSET=J0,GENERATE
*NSET,NSET=J1,GENERATE
*NSET,NSET=J2,GENERATE
*NSET,NSET=J3,GENERATE
*NSET,NSET=J4,GENERATE
*NSET,NSET=J5,GENERATE
*NSET,NSET=J6,GENERATE
*NSET,NSET=J7,GENERATE
*NSET,NSET=J8,GENERATE
*NSET,NSET=M0,GENERATE
*NSET,NSET=M1,GENERATE
*NSET,NSET=M2,GENERATE
*NSET,NSET=M3,GENERATE
*NSET,NSET=M4,GENERATE
*NSET,NSET=M5,GENERATE
*NSET,NSET=M6,GENERATE
*NSET,NSET=M7,GENERATE
*NSET,NSET=M8,GENERATE
*MPC
PIN,M0,10
PIN,M1,11
PIN,M2,12
PIN,M3,13
PIN,M4,14
PIN,M5,15
PIN,M6,16
PIN,M7,17
PIN,M8,18
*NSET,NSET=NALL
XYPL,YZPL,
RAD0,RAD1,RAD2,RAD3,RAD4,RAD5,RAD6,RAD7,RAD8,RAD9,
RADA,RADB,RADC,RADD,RADE,RADF,RADG,
LEV0,LEV1,LEV2,LEV3,LEV4,LEV5,LEV6,LEV7,LEV8,
J0,J1,J2,J3,J4,J5,J6,J7,J8
*ELEMENT,TYPE=C3D20R
*ELGEN,ELSET=MAST
*ELGEN,ELSET=FAN1
*ELEMENT,TYPE=C3D20R,ELSET=FAN1
*ELGEN,ELSET=FAN1
*ELGEN
*ELGEN,ELSET=FAN2
*ELGEN,ELSET=FAN3
*ELEMENT,TYPE=C3D20R,ELSET=FAN3
*ELGEN,ELSET=FAN3
*ELGEN,ELSET=FAN4
*ELEMENT,TYPE=C3D20R,ELSET=FAN4
*ELGEN,ELSET=FAN4
*ELGEN,ELSET=FAN5
*ELEMENT,TYPE=C3D20R,ELSET=FAN5
*ELGEN,ELSET=FAN5
*ELGEN,ELSET=FAN6
*ELEMENT,TYPE=C3D20R,ELSET=FAN6
*ELGEN,ELSET=FAN6
*ELGEN,ELSET=FAN7
*ELGEN,ELSET=FAN8
*ELSET,ELSET=LOAD,GENERATE
*ELSET,ELSET=PRINT
*ELSET,ELSET=PLOTA
LOAD,
*ELSET,ELSET=PLOTA,GENERATE
*ELSET,ELSET=PLOTB,GENERATE
*ELSET,ELSET=EALL,GENERATE
*surface, NAME=SURFACE
 LOAD, S5
*MATERIAL,NAME=MAT
*ELASTIC
30.E6,.3
*SOLID SECTION,MATERIAL=MAT,ELSET=EALL
*BOUNDARY
RAD0,2
LEV0,2
138,3
YZPL,1
**
*STEP
*STATIC
*DSLOAD
SURFACE,P,-100.
*CONTOUR INTEGRAL,CONTOURS=3,SYMM,OUTPUT=BOTH,NORMAL
J0,J1,J2,J3,J4,J5,J6,J7,J8
*CONTOUR INTEGRAL,CONTOURS=3,SYMM,OUTPUT=BOTH,NORMAL,TYPE=K FACTORS
J0,J1,J2,J3,J4,J5,J6,J7,J8
*CONTOUR INTEGRAL,CONTOURS=3,SYMM,OUTPUT=BOTH,NORMAL,TYPE=T-STRESS
J0,J1,J2,J3,J4,J5,J6,J7,J8
*EL PRINT,ELSET=PRINT,POSITION=CENTROID
S,MISES
*ENERGY PRINT
*NODE PRINT
U,
*ENERGY FILE
*OUTPUT,FIELD, OP=NEW
*ELEMENT OUTPUT, ELSET=PRINT
S,
*NODE OUTPUT,NSET=NALL
U,
*OUTPUT,HISTORY,OP=NEW
*ENERGY OUTPUT,VARIABLE=ALL
*ELEMENT OUTPUT, ELSET=PRINT
S,
*NODE OUTPUT,NSET=NALL
U,
*END STEP