*HEADING
ECA4GYP5: PATCH TEST FOR AXISYMMETRIC ELEMENTS WITH TWIST  [CGAX4RH]
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE,NSET=NALL
1, 0.0, 0.0
2, 1.0, 0.0
3, 1.0, 1.0
4, 0.0, 1.0
5, 0.3, 0.2
6, 0.6, 0.3
7, 0.8, 0.8 
8, 0.2, 0.6
*ELEMENT, TYPE=CGAX4RH, ELSET=PATCH
1,1,2,6,5
2,2,3,7,6
3,7,3,4,8
4,1,5,8,4
5,5,6,7,8
*NSET, NSET=BOT
1,2
*NSET, NSET=TOP
3,4
*NSET, NSET=RIGHT
2,3
*NSET, NSET=LEFT
1,4
*SOLID SECTION,MATERIAL=A1, ELSET=PATCH
*MATERIAL,NAME=A1
*ELASTIC,TYPE=ISOTROPIC
1.E6,.25
*STEP,PERTURBATION
*STATIC
*BOUNDARY
LEFT,1,1,0.
BOT,2,5,0.
TOP,2,2,0.
TOP,5,5,0.01
*EL PRINT
S,MISES
E,
ENER,
COORD, IVOL
NFORC,
*EL PRINT,POSITION=NODES
S,MISES
E,
*EL PRINT,POSITION=AVERAGED AT NODES
S,MISES
E,
*NODE PRINT
U,
RF,
CF,
COORD, 
*EL FILE
S, 
E,
ENER,
COORD, IVOL
NFORC,
*NODE FILE
U,
RF,
CF,
COORD, 
*END STEP
*STEP,PERTURBATION
*STATIC
*BOUNDARY
1,1,1,0.
1,2,2,0.
2,1,1,1.E-3
2,2,2,1.E-3
3,1,1,1.E-3
3,2,2,2.E-3
4,1,1,0.
4,2,2,1.E-3
TOP,5,5,0.
BOT,5,5,0.
*DLOAD
PATCH,BZNU
*EL PRINT
S,MISES
E,
ENER,
COORD, IVOL
NFORC,
*EL PRINT,POSITION=NODES
S,MISES
E,
*EL PRINT,POSITION=AVERAGED AT NODES
S,MISES
E,
*NODE PRINT
U,
RF,
CF,
COORD, 
*EL FILE
S, 
E,
ENER,
COORD, IVOL
NFORC,
*NODE FILE
U,
RF,
CF,
COORD, 
*END STEP
*STEP,NLGEOM
*STATIC
*BOUNDARY
1,1,1,0.
1,2,2,0.
2,1,1,1.E-2
2,2,2,0.
3,1,1,1.E-2
3,2,2,1.E-2
4,1,1,0.
4,2,2,1.E-2
TOP,5,5,0.
BOT,5,5,0.
*EL PRINT
S,MISES
E,
ENER,
COORD, IVOL
NFORC,
*EL PRINT,POSITION=NODES
S,MISES
E,
*EL PRINT,POSITION=AVERAGED AT NODES
S,MISES
E,
*NODE PRINT
U,
RF,
CF,
COORD, 
*EL FILE
S, 
E,
ENER,
COORD, IVOL
NFORC,
*NODE FILE
U,
RF,
CF,
COORD, 
*END STEP
*STEP,PERTURBATION
*STATIC
*BOUNDARY
1,1,1,0.
1,2,2,0.
2,1,1,1.01E-3
2,2,2,1.01E-3
3,1,1,1.01E-3
3,2,2,2.02E-3
4,1,1,0.
4,2,2,1.01E-3
TOP,5,5,0.
BOT,5,5,0.
*DLOAD,OP=NEW
PATCH,BZNU
*EL PRINT
S,MISES
E,
ENER,
COORD, IVOL
NFORC,
*EL PRINT,POSITION=NODES
S,MISES
E,
*EL PRINT,POSITION=AVERAGED AT NODES
S,MISES
E,
*NODE PRINT
U,
RF,
CF,
COORD, 
*EL FILE
S, 
E,
ENER,
COORD, IVOL
NFORC,
*NODE FILE
U,
RF,
CF,
COORD, 
*END STEP