*Heading
:  Axisymmetric Elements With Twist [cgax6mh] 
*Restart,Write
*Node,Nset=Nall
 1,
 7,0.,1.5
 801,8.,
 807,8.,1.5
*Ngen,Nset=Nl
 1,7,1
*Ngen,Nset=Nr
 801,807,1
*Nfill,Nset=Nall
 Nl,Nr,8,100
*Nset,Nset=Nbot,Gen
 1,801,100
*Nset,Nset=Ntop,Gen
 7,807,100
*Element,Type=cgax6mh,Elset=Eall
 1,1,201,203,101,202,102
 13,203,3,1,103,2,102
*Elgen,Elset=Eall
 1,3,2,1,4,200,3
 13,3,2,1,4,200,3
*Solid Section,Material=A1, Elset=Eall
 1., 
*Material,Name=A1
*Elastic,Type=Isotropic
 30.0e6,0.3
*Elset,Elset=Load
 10,11,12
*Elset,Elset=El10
 10,22
*Nset,Nset=N1
 201,202,203,204,205,301,302,303,304,305,401,402,403,404,405
*Nset,Nset=N2
 403,404,405,406,407,503,504,505,506,507,603,604,605,606,607
*Boundary
 Nl,1,1
 1,2
 Nall,5,5
*Step
*Static
 .3,1.
*Dload,Op=New
 Load, P2, -1000000.0
*Node File, Nset=N1, Freq=2
 U,
 Coord,
*Node File, Nset=N2, Freq=2
 U,
*Output, Field
*Node Output
 U,
 Coord,
*Element Output
 S,
 E,
*El Print,Elset=El10,Freq=2
 S, 
*Node Print,Nset=N1,Freq=2
 U, 
*Node Print,Nset=N2,Freq=2
 U, 
*End Step
*Step
*Static
 0.3,1.
*Dload,Op=New
*Boundary,Op=New
 Nbot,1,2
 NBot,5,5
 Ntop,5,5,0.01
*Endstep