*HEADING
 NAFEMS R0031/2 WRAPPED THICK CYLINDER
*RESTART,WRITE
*NODE
 10001, 
 1,0.,23.
 101,0.,24.
 201,0.,25.
 301,0.,26.
 401,0.,27.
 9,23.,0.
 109,24.,0.
 209,25.,0.
 309,26.,0.
 409,27.,0.
*NGEN,LINE=C,NSET=LHS
 1,9,1,10001
 101,109,1,10001
 201,209,1,10001
 301,309,1,10001
 401,409,1,10001
*NCOPY,CHANGE NUMBER=8000,OLD SET=LHS,SHIFT,NEW SET=RHS
 0.,0.,-100.
 1.,0.,0.,0.,1.,0.,0.
*ELEMENT,TYPE=SC8R,ELSET=ELIN
 1,1,2,1002,1001,101,102,1102,1101,
*ELGEN,ELSET=ELIN
 1,8,1,1,8,1000,1000,2,100,100
*ELEMENT,TYPE=SC8R,ELSET=ELOUT
 201,201,202,1202,1201,301,302,1302,1301
*ELGEN,ELSET=ELOUT
 201,8,1,1,8,1000,1000,2,100,100
*ELSET,ELSET=OUTPUT,GEN
7,307,100
*ELSET,ELSET=ELOAD, GENERATE
    1,   8,1
 1001,1008,1
 2001,2008,1
 3001,3008,1
 4001,4008,1
 5001,5008,1
 6001,6008,1
 7001,7008,1
*NFILL,NSET=NALL
 LHS,RHS,8,1000
*NSET,NSET=AB,GEN
   9,8009,1000
 109,8109,1000
 209,8209,1000
 309,8309,1000
 409,8409,1000
*NSET,NSET=CD,GEN
 1,8001,1000
 101,8101,1000
 201,8201,1000
 301,8301,1000
 401,8401,1000
*SHELL SECTION,ELSET=ELIN,MATERIAL=INNER,ORIENTATION=ORIENT
 1,
*MATERIAL,NAME=INNER
*ELASTIC
 2.1E5,.3
*EXPANSION
 2.E-5, 
*DENSITY
 1.E-4, 
*SHELL SECTION,ELSET=ELOUT,MATERIAL=OUTER,ORIENTATION=ORIENT
 1,
*MATERIAL,NAME=OUTER
*ELASTIC,TYPE=ENGINEERING CONSTANTS
 1.3E5,5.E3,5.E3,.25,.25,0.,1.E4,1.E4
 5.E3, 
*EXPANSION,TYPE=ORTHO
 3.E-6,2.E-5,2.E-5
*DENSITY
 1.E-4, 
*ORIENTATION,NAME=ORIENT,SYSTEM=CYLINDRICAL,DEFINITION=COORDINATES
 0.,0.,0.,0.,0.,-100.
 1,0.
*TRANSFORM,NSET=NALL,TYPE=C
 0.,0.,0.,0.,0.,-100.
*BOUNDARY
 LHS,3
 AB,2
 CD,2
*STEP
*STATIC
*DLOAD
 ELOAD,P1,200.
*OUTPUT, FIELD
*NODE OUTPUT
 U,
*ELEMENT OUTPUT
 1,2,3,4,5
 S,
*END STEP
*STEP
*STATIC
*TEMPERATURE
 NALL,130.
*END STEP