*HEADING
:  AXISYMMETRIC ELEMENTS WITH TWIST [CGAX4r]
*RESTART,WRITE
*NODE,NSET=NL
*NODE,NSET=NR
*NFILL,NSET=NALL
NL,NR,4,100
*NSET,NSET=NBOT,GEN
*NSET,NSET=NTOP,GEN
*ELEMENT,TYPE=CGAX4R, ELSET=EALL
*ELGEN,ELSET=EALL
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
 1.0,
*MATERIAL,NAME=A1
*ELASTIC,TYPE=ISOTROPIC
 30.0E6,0.3
*ELSET,ELSET=LOAD
*ELSET,ELSET=EL10
*NSET,NSET=MIDDLE
*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,ELSET=EL10
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