*HEADING
ECA8HHP5: PATCH TEST FOR AXISYMMETRIC ELEMENTS WITH TWIST  [CGAX8HT]
*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=CGAX8HT, 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
*CONDUCTIVITY
4.85E-4, 
*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, 
*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
DGP,
NFORC,
*EL PRINT,POSITION=NODES
S,MISES
E,
DGP,
*EL PRINT,POSITION=AVERAGED AT NODES
S,MISES
E,
DGP,
*NODE PRINT
U,
RF,
CF,
COORD, 
*EL FILE
S, 
E,
ENER,
COORD, IVOL
DGP,
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
*STEP
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
*BOUNDARY,OP=NEW
1,11, ,   0.0
2,11, , 200.0
3,11, , 300.0
4,11, , 100.0
9,11, , 100.0
13,11, ,250.0
16,11, ,200.0
20,11, , 50.0
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.
*DFLUX
PATCH,BFNU
*EL PRINT
S,MISES
E,
ENER,
COORD, IVOL
DGP,
NFORC,
*EL PRINT,POSITION=NODES
S,MISES
E,
DGP,
*EL PRINT,POSITION=AVERAGED AT NODES
S,MISES
E,
DGP,
*NODE PRINT
U,
RF,
CF,
COORD, 
*EL FILE
S, 
E,
ENER,
COORD, IVOL
DGP,
NFORC,
*NODE FILE
U,
RF,
CF,
COORD, 
*NODE PRINT, FREQUENCY=0
*EL PRINT
 HFL,
*EL PRINT,POSITION=NODES
 HFL,
*EL PRINT,POSITION=AVERAGED AT NODES
 HFL,
*EL PRINT,POSITION=CENTROID
 HFL,
*EL FILE
 HFL,
*EL FILE,POSITION=NODES
 HFL,
*EL FILE,POSITION=AVERAGED AT NODES
 HFL,
*EL FILE,POSITION=CENTROID
 HFL,
*END STEP