*HEADING
 INFINITELY LONG CYLINDRICAL TUBE SUBJECTED TO  
 INTERNAL PRESSURE
 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
 11,  41, 161, 163, 43, 2041, 2161, 2163, 2043, 101, 162, 103, 42,  
      2101, 2162, 2103, 2042, 1041, 1161, 1163, 1043
*ELGEN,ELSET=Ei
  1,   4, 2, 1
*ELGEN,ELSET=Eo
 11,   4, 2, 1
*ELSET,GENERATE,ELSET=Inside
  1, 4
*ELSET,ELSET=Eioutp
 1, 
*ELSET,ELSET=Eooutp
 11 , 
*SOLID SECTION,ELSET=Ei,COMPOSITE,ORIENTATION=Ocircum,STACKDIR=1
 1.0, 3, Mcircum
 1.0, 3, Maxial
*SOLID SECTION,ELSET=Eo,COMPOSITE,ORIENTATION=Ocircum,STACKDIR=1
 1.0, 3, Mcircum
 1.0, 3, Maxial
 1.0, 3, Mcircum
 1.0, 3, Maxial
 1.0, 3, Mcircum
 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=Eioutp
 1,2,3,4,5,6
 COORD, S11, S22, S33
*EL PRINT,ELSET=Eooutp
 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16
 COORD, S11, S22, S33
*EL PRINT,ELSET=Eooutp
 17,18
 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