*HEADING
 Hemispherical shell section under free thermal expansion [SC8R]
*NODE
  999, 0.,        0.0,   0.0
    1, 0.,      0.0, 100.2          
   17, 95.2959, 0.0,  30.9635
10001, 0.,      0.0,  99.8      
10017, 94.9154, 0.0,  30.8399
*NGEN,LINE=C,NSET=XZPLANE
     1,17,1,999
 10001,10017,1,999
*NCOPY,OLDSET=XZPLANE,NEWSET=ALL,SHIFT,MULT=16,CHANGENUM=100
 0.,0.,0.
 0.,0.,0.,     -1.,0.,0., 5.625
*NSET,NSET=ALL
 ALL,XZPLANE
*NSET,NSET=XYSYMM,GENERATE
  1601,1617,2
 11601,11617,2
*NSET,NSET=XZSYMM,GENERATE
     1,   17,2
 10001,10017,2
*NSET,NSET=YZSYMM,GENERATE
     1, 1601,200
 10001,11601,200
*ELEMENT,TYPE=SC6R,ELSET=SHELL
 1,10001,10003,10203,1,3,203
 10000,10203,10201,10001,203,201,1
*ELGEN,ELSET=SHELL
 1,8,2,2,8, 200,200
 10000,8,2,2,8, 200,200
*SHELL SECTION,ELSET=SHELL,MAT=ELA, POISSON=ELASTIC
.4,3
*MATERIAL,NAME=ELA
*ELASTIC
 6.825e7,.3
*EXPANSION
 1.e-6,
*INITIAL CONDITIONS,TYPE=TEMP
 all,0.0
*BOUNDARY
 xysymm,3,3,
 xzsymm,2, ,
 yzsymm,1, ,
**
**  STEP 1: PERTUBATION STEP
**
*STEP,PERTURBATION
*STATIC
*TEMPERATURE
 ALL,430.
*OUTPUT,FIELD
*NODE OUTPUT
 u,
*OUTPUT,FIELD
*ELEMENT OUTPUT
 E,
 STH,
*END STEP
**
**  STEP 2: NLGEOM STEP
**
*STEP, NLGEOM
 0.1,1.0
*STATIC
*TEMPERATURE
 ALL,430.
*controls, parameters=field
 , , , 1.0
*OUTPUT,FIELD,FREQ=1000
*NODE OUTPUT
 u,
*OUTPUT,FIELD
*ELEMENT OUTPUT
 E,
 STH,
*END STEP