*HEADING
 : J-INTEGRALS
*RESTART,WRITE
*NODE
1,1.
21,6.
821,5.,3.
1601,1.
1621,2.5,5.
3001,
3021,0.,5.
*NGEN,NSET=CRACKTIP
1,1601,100
*NGEN,NSET=OUTER
21,821,100
821,1621,100
*NFILL,SINGULAR=1
CRACKTIP,OUTER,10,2
*NSET,NSET=N3001
3001,
*NSET,NSET=N3021
3021,
*NFILL,SINGULAR=1,NSET=NX0
N3001,N3021,10,2
*NSET,NSET=NB1,GEN
1601,1621,2
*NFILL,BIAS=.925
NB1,NX0,14,100
*NSET,NSET=NB3,GEN
1701,3001,100
*NSET,NSET=NB4,GEN
1705,3005,100
*NFILL
NB3,NB4,2,2
*NSET,NSET=N1601
1601,
*NSET,NSET=N1801
1801,
*NFILL,SINGULAR=1
N1601,N1801,2,100
*NSET,NSET=NB5,GEN
21,2721,100
*NSET,NSET=NB5B
2821,
*NSET,NSET=NB5A
3021,
*NCOPY,OLD SET=NB5,CHANGE NUMBER=10000,POLE,NEW SET=FAR
1,
*NCOPY,OLD SET=NB5A,CHANGE NUMBER=10000,POLE,NEW SET=FARA
3001,
*NCOPY,OLD SET=NB5B,CHANGE NUMBER=10000,POLE,NEW SET=FARB
2401,
*ELEMENT,TYPE=CPE8R
1,1,5,205,201,3,105,203,101
*ELGEN,ELSET=CONT
1,5,4,4,15,200,200
*ELEMENT,TYPE=CINPE5R
221,221,21,10021,10221,121
*ELGEN,ELSET=INF
221,15,200,200
*ELSET,ELSET=EALL
CONT,INF
*SOLID SECTION,ELSET=EALL,MATERIAL=M1,ORIENTATION=ORN1
*MATERIAL,NAME=M1
*ELASTIC
200.E9,.3
*ELSET,ELSET=LOAD,GEN
1601,2801,200
*NSET,NSET=XAXIS,GEN
1,21,2
*NSET,NSET=XAXIS1
XAXIS,10021
*NSET,NSET=YAXIS
*NSET,NSET=CRTIP,GEN
101,1601,100
*ORIENTATION,NAME=ORN1,SYSTEM=CYLINDRICAL
1.0,0.0,0.0, 1.0,0.0,1.0
*EQUATION
2,
CRTIP,     1,     1.
    1,     1,    -1.
2,
CRTIP,     2,     1.
    1,     2,    -1.
**
** STEP 1 - UNIFORM LOAD
**
*STEP
*STATIC
*DLOAD
LOAD,P4NU
*BOUNDARY
XAXIS1,2
21,1,1,0.
*CONTOUR INTEGRAL,TYPE=J,CONTOUR=5,SYMM,OUTPUT=FILE
CRACKTIP,1.0,
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*ENDSTEP
**
** STEP 2 - LINEAR LOAD
**
*STEP
*STATIC
*DLOAD
LOAD,P4NU
*BOUNDARY
XAXIS1,2
21,1,1,0.
*CONTOUR INTEGRAL,TYPE=J,CONTOUR=5,SYMM,OUTPUT=FILE
CRACKTIP,1.0,
*ENDSTEP
**
** STEP 3 - QUADRATIC LOAD
**
*STEP
*STATIC
*DLOAD
LOAD,P4NU
*BOUNDARY
XAXIS1,2
21,1,1,0.
*CONTOUR INTEGRAL,TYPE=J,CONTOUR=5,SYMM,OUTPUT=FILE
CRACKTIP,1.0
*ENDSTEP
**
** STEP 4 - CUBIC LOAD
**
*STEP
*STATIC
*DLOAD
LOAD,P4NU
*BOUNDARY
XAXIS1,2
21,1,1,0.
*CONTOUR INTEGRAL,TYPE=J,CONTOUR=5,SYMM,OUTPUT=FILE
CRACKTIP,1.0,
*ENDSTEP
**
** STEP 5 - QUARTIC LOAD
**
*STEP
*STATIC
*DLOAD
LOAD,P4NU
*BOUNDARY
XAXIS1,2
21,1,1,0.
*CONTOUR INTEGRAL,TYPE=J,CONTOUR=5,SYMM,OUTPUT=FILE
CRACKTIP,1.0,
*ENDSTEP