*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, 21, 23, 3,   2001, 2021, 2023, 2003,   11, 22, 13, 2,   2011, 2022, 
      2013, 2002,   1001, 1021, 1023, 1003
*ELGEN,ELSET=Eall
 1,   4, 2, 1,   8, 20, 10 
*ELSET,GENERATE,ELSET=Ecircum
  1,  4
 21, 24
 41, 44
 61, 64
*ELSET,GENERATE,ELSET=Eaxial
 11, 14
 31, 34
 51, 54
 71, 74
*ELSET,GENERATE,ELSET=Inside
  1,  4
*ELSET,GENERATE,ELSET=Eoutp
 1, 71, 10  
*SOLID SECTION,ELSET=Ecircum,COMPOSITE,STACKDIR=1
 1.0, 3, Mcircum, Ocircum
*SOLID SECTION,ELSET=Eaxial, COMPOSITE,STACKDIR=1,ORIENTATION=Ocircum
 1.0, 3, Maxial
*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 , 
*MATERIAL,NAME=Maxial
*ELASTIC,TYPE=ENGINEERING CONSTANTS
 10.0E3, 10.0E3, 250.0E3,  0.25, 0.01, 0.01, 2.0E3, 5.0E3,
 5.0E3, 
*ORIENTATION,NAME=Ocircum,SYSTEM=CYLINDRICAL
 0.0, 0.0, 0.0,    0.0, 0.0, 30.0
 1, 0.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
 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