*HEADING
 INFINITELY LONG CYLINDRICAL TUBE SUBJECTED TO INTERNAL PRESSURE     [C3D20R] 
 UNITS: N, mm
*RESTART,WRITE
*NODE
   1,      60.0,   0.0
   9,       0.0,  60.0
 161,     140.0,   0.0
 169,       0.0, 140.0
*NGEN,LINE=C,NSET=Di
   1,   9, 1, , 0.0, 0.0, 0.0 
*NGEN,LINE=C,NSET=Do
 161, 169, 1, , 0.0, 0.0, 0.0 
*NFILL,NSET=Face1
 Di, Do, 16, 10 
*NCOPY,CHANGE NUMBER=2000,OLD SET=Face1,SHIFT,NEW SET=Face2
 0.0, 0.0, 300.0
 0.0,
*NFILL,NSET=Nall
 Face1, Face2, 2, 1000
*NSET,GENERATE,NSET=Xsymm
    1,  161, 10
 1001, 1161, 10
 2001, 2161, 10
*NSET,GENERATE,NSET=Ysymm
    9,  169, 10
 1009, 1169, 10
 2009, 2169, 10
*NSET,GENERATE,NSET=Noutp
 1, 161, 10
**
*ELEMENT,TYPE=C3D20R
  1,  1, 41, 43, 3,   2001, 2041, 2043, 2003,   21, 42, 23, 2,   2021, 2042, 
      2023, 2002,   1001, 1041, 1043, 1003
*ELGEN,ELSET=Eall
  1,   4, 2, 1,   4, 40, 10 
*ELSET,GENERATE,ELSET=Inside
  1, 4
*ELSET,GENERATE,ELSET=Eoutp
 1, 41, 10
*SOLID SECTION,ELSET=Eall,COMPOSITE,STACKDIR=1
 1.0, 3, Mcircum,  Ocircum
 1.0, 3, Mcircum,  Oaxial
*MATERIAL,NAME=Mcircum
*ELASTIC,TYPE=ENGINEERING CONSTANTS
 10.0E3, 250.0E3, 10.0E3,  0.01, 0.25, 0.25, 5.0E3, 2.0E3, 
 5.0E3 , 
*ORIENTATION,NAME=Ocircum,SYSTEM=CYLINDRICAL
 0.0, 0.0, 0.0,    0.0, 0.0, 30.0
 1, 0.0
*ORIENTATION,NAME=Oaxial,SYSTEM=CYLINDRICAL
 0.0, 0.0, 0.0,    0.0, 0.0, 30.0
 1, 90.0
*BOUNDARY
 Xsymm, 2
 Ysymm, 1
 Face1, 3
 Face2, 3
*******************
*STEP,PERTURBATION
*STATIC
*DLOAD
 Inside, P6, 50.0
*EL PRINT,ELSET=Eoutp
 1,2,3,4,5,6
 COORD, S11, S22, S33
*EL FILE
 S, E
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,E
*NODE PRINT,NSET=Noutp
 U1, 
*NODE FILE
 U,RF
*OUTPUT,FIELD
*NODE OUTPUT
U,RF
*END STEP