*HEADING
  PSPTSKF2D: Center Slant Cracked Plate under Tension
  in orthotropic elasticity specified by the elastic stiffness constants 
  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
*ORIENTATION,NAME=RECT
 1.0, 0.0, 0.0,  0.0, 1.0, 0.0
 1, 0.0
*SOLID SECTION,ELSET=EALL,MATERIAL=M,ORIENT=RECT
*MATERIAL,NAME=M
*ELASTIC,TYPE=ORTHOTROPIC
1.23E+11,2.24E+11,4.79E+11,4.74E+10,4.21E+10,1.21E+11,7.69E+10,7.69E+10,
9.00E+9, 
*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
*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=J
TIPL,-0.70710678,-0.70710678
*CONTOUR INTEGRAL,CONTOUR=5,OUTPUT=BOTH,TYPE=K FACTORS
TIPL,-0.70710678,-0.70710678
*EL PRINT,freq=0
*NODE PRINT,NSET=LOW
RF2, 
*NODE PRINT,NSET=UP
RF2, 
*ENDSTEP