*HEADING
:  AXISYMMETRIC ELEMENTS WITH TWIST [CGAX4r]
*RESTART,WRITE
*NODE,NSET=NL
 1,  0.0, 0.0
 2,  0.0, 0.5
 3,  0.0, 1.0
 4,  0.0, 1.5
*NODE,NSET=NR
 401,  8.0, 0.0
 402,  8.0, 0.5
 403,  8.0, 1.0
 404,  8.0, 1.5
*NFILL,NSET=NALL
NL,NR,4,100
*NSET,NSET=NBOT,GEN
1,401,100
*NSET,NSET=NTOP,GEN
4,404,100
*ELEMENT,TYPE=CGAX4R, ELSET=EALL
 1,1,101,102,2
*ELGEN,ELSET=EALL
 1,3,1,1,4,100,3,
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
 1.0, 
*MATERIAL,NAME=A1
*ELASTIC,TYPE=ISOTROPIC
 30.0E6,0.3
*ELSET,ELSET=LOAD
 10,11,12
*ELSET,ELSET=EL10
 10, 
*NSET,NSET=MIDDLE
 401,402,403,404,405,406,407,408
*BOUNDARY
 NL,1,1
 1,2
 NALL,5,5
*STEP
*STATIC
1.,1.
*DLOAD,OP=NEW
 LOAD, P2, -1000000.0
*EL FILE,ELSET=EL10
 S,E
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,E
*NODE FILE
 U, 
*OUTPUT,FIELD
*NODE OUTPUT
U,
*EL PRINT,ELSET=EL10
 E, 
*NODE PRINT
 U, 
*END STEP
*STEP
*STATIC
*DLOAD,OP=NEW
*BOUNDARY,OP=NEW
NBOT,1,5
NTOP,5,5,0.01
*ENDSTEP