*HEADING
  Center Slant Cracked Plate under Tension
  in isotropic elasticity
  theta=45 degrees, a=1.0, w=2.5, h/w=2, a/w=0.4
  File name: .inp
**
**  INCLUDE NODE DATA
*NODE
*INCLUDE,INPUT=psptskf2d_node.inp
**
**  INCLUDE ELEMENT DATA
*ELEMENT,TYPE=CPE8,ELSET=EALL    
*INCLUDE,INPUT=psptskf2d_element.inp
*SOLID SECTION,ELSET=EALL,MATERIAL=M
*MATERIAL,NAME=M
*ELASTIC
200.E9,.3
*ELSET,ELSET=EL_LOAD,gen
849,853
891,896
867,878
*BOUNDARY
**875,2
**899,2
863,1
887,1
*NSET,NSET=TIPR
1,
*NSET,NSET=TIPL
332,
*NSET,NSET=LOW,GEN
863,869
905,910
2622,2632,2
2706,2716,2
*NSET,NSET=UP,GEN
881,893
2658,2680,2
**
** STEP 1: Calculate J-integrals, stress intensity factors, and T-stresss;
**         and
**         MTS criterion is used to calculate crack propogation direction.
**
*STEP
*STATIC
*BOUNDARY
LOW,2,,-0.01
UP, 2,, 0.01
**DLOAD
**EL_LOAD,P3,-1.0
*CONTOUR INTEGRAL,CONTOUR=5,OUTPUT=BOTH,TYPE=J
TIPR, 0.70710678, 0.70710678
*CONTOUR INTEGRAL,CONTOUR=5,OUTPUT=BOTH,TYPE=K FACTORS
TIPR, 0.70710678, 0.70710678
*CONTOUR INTEGRAL,CONTOUR=5,OUTPUT=BOTH,TYPE=T-STRESS
TIPR, 0.70710678, 0.70710678
*CONTOUR INTEGRAL,CONTOUR=5,OUTPUT=BOTH,TYPE=J
TIPL,-0.70710678,-0.70710678
*CONTOUR INTEGRAL,CONTOUR=5,OUTPUT=BOTH,TYPE=K FACTORS,DIRECTION=MTS
TIPL,-0.70710678,-0.70710678
*CONTOUR INTEGRAL,CONTOUR=5,OUTPUT=BOTH,TYPE=T-STRESS
TIPL,-0.70710678,-0.70710678
*EL PRINT,freq=0
*NODE PRINT,NSET=LOW
RF2, 
*NODE PRINT,NSET=UP
RF2, 
*ENDSTEP
**
** STEP 2: Calculate stress intensity factors;
**         and 
**         MERR criterion is used to calculate crack propogation direction.
**
*STEP
*STATIC
*BOUNDARY
LOW,2,,-0.01
UP, 2,, 0.01
**DLOAD
**EL_LOAD,P3,-1.0
*CONTOUR INTEGRAL,CONTOUR=5,OUTPUT=BOTH,TYPE=K FACTORS,DIRECTION=MERR
TIPR, 0.70710678, 0.70710678
*CONTOUR INTEGRAL,CONTOUR=5,OUTPUT=BOTH,TYPE=K FACTORS,DIRECTION=MERR
TIPL,-0.70710678,-0.70710678
*EL PRINT,freq=0
*NODE PRINT,NSET=LOW,freq=0
RF2, 
*NODE PRINT,NSET=UP,freq=0
RF2, 
*ENDSTEP
**
** STEP 3: Calculate stress intensity factors; 
**         and
** STEP 1: KII0 criterion is used to calculate crack propogation direction.
**
*STEP
*STATIC
*BOUNDARY
LOW,2,,-0.01
UP, 2,, 0.01
**DLOAD
**EL_LOAD,P3,-1.0
*CONTOUR INTEGRAL,CONTOUR=5,OUTPUT=BOTH,TYPE=K FACTORS,DIRECTION=KII0
TIPR, 0.70710678, 0.70710678
*CONTOUR INTEGRAL,CONTOUR=5,OUTPUT=BOTH,TYPE=K FACTORS,DIRECTION=KII0
TIPL,-0.70710678,-0.70710678
*EL PRINT,freq=0
*NODE PRINT,NSET=LOW,freq=0
RF2, 
*NODE PRINT,NSET=UP,freq=0
RF2, 
*ENDSTEP