*HEADING
ECA8GYP5: PATCH TEST FOR AXISYMMETRIC ELEMENTS WITH TWIST  [CGAX8RH]
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE
       1,          0.,          0.
       2,          1.,          0.
       3,          1.,          1.
       4,          0.,          1.
       5,         0.3,         0.2
       6,         0.6,         0.3
       7,         0.8,         0.8
       8,         0.2,         0.6
       9,         0.5,          0.
      10,         0.8,        0.15
      11,        0.45,        0.25
      12,        0.15,         0.1
      13,          1.,         0.5
      14,         0.9,         0.9
      15,         0.7,        0.55
      16,         0.5,          1.
      17,         0.1,         0.8
      18,         0.5,         0.7
      19,        0.25,         0.4
      20,          0.,         0.5
*ELEMENT, TYPE=CGAX8RH, ELSET=PATCH
       1,       1,       2,       6,       5,       9,      10,
      11,      12
       2,       2,       3,       7,       6,      13,      14,
      15,      10
       3,       7,       3,       4,       8,      14,      16,
      17,      18
       4,       1,       5,       8,       4,      12,      19,
      17,      20
       5,       5,       6,       7,       8,      11,      15,
      18,      19
*NSET, NSET=BOT
1,9,2
*NSET, NSET=TOP
4,16,3
*NSET, NSET=RIGHT
2,13,3
*NSET, NSET=LEFT
1,20,4
*SOLID SECTION,MATERIAL=A1, ELSET=PATCH
*MATERIAL,NAME=A1
*ELASTIC,TYPE=ISOTROPIC
1.E6,.25
*surface, NAME=CROSS
1,S4
1,S3
2,S3
2,S2
*surface, NAME=USHAPE
4,S1
4,S3
5,S1
5,S2
5,S3
*surface, NAME=INTER
5,S1
5,S2
5,S3
5,S4
*surface, NAME=EXTER
1,S3
2,S3
3,S4
4,S2
*STEP,PERTURBATION
*STATIC
*BOUNDARY
LEFT,1,1,0.
TOP,2,2,0.
BOT,2,5,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, 
*SECTION PRINT, NAME=CROSS, SURFACE=CROSS, AXES=LOCAL
SOAREA,SOF
*SECTION FILE, NAME=CROSS, SURFACE=CROSS, AXES=LOCAL
SOAREA,SOF
*SECTION PRINT, NAME=USHAPE, SURFACE=USHAPE, AXES=LOCAL
SOAREA,SOF
*SECTION FILE, NAME=USHAPE, SURFACE=USHAPE, AXES=LOCAL
SOAREA,SOF
*SECTION PRINT, NAME=INTER, SURFACE=INTER, AXES=LOCAL
SOAREA,SOF
*SECTION FILE, NAME=INTER, SURFACE=INTER, AXES=LOCAL
SOAREA,SOF
*SECTION PRINT, NAME=EXTER, SURFACE=EXTER, AXES=LOCAL
SOAREA,SOF
*SECTION FILE, NAME=EXTER, SURFACE=EXTER, AXES=LOCAL
SOAREA,SOF
*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
9,1,1,5.E-4
9,2,2,5.E-4
13,1,1,1.E-3
13,2,2,1.5E-3
16,1,1,5.E-4
16,2,2,1.5E-3
20,1,1,0.
20,2,2,5.E-4
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
9,1,1,5.E-3
9,2,2,0.
13,1,1,1.E-2
13,2,2,5.E-3
16,1,1,5.E-3
16,2,2,1.E-2
20,1,1,0.
20,2,2,5.E-3
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
9,1,1,5.05E-4
9,2,2,5.05E-4
13,1,1,1.01E-3
13,2,2,1.515E-3
16,1,1,5.05E-4
16,2,2,1.515E-3
20,1,1,0.
20,2,2,5.05E-4
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