*HEADING
.INP:  HEMISPHERICAL SHELL SECTION UNDER FREE THERMAL EXPANSION [S4]
*NODE
999,0.,0.,0.
***  COORD               NORMAL
1, 0.,0., 100.,           0,0,1.
17, 95.106, 0., 30.902,   .95106, 0., .30902
*NGEN,LINE=C,NSET=XZPLANE
1,17,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
*NSET,NSET=XZSYMM,GENERATE
1,17,2
*NSET,NSET=YZSYMM,GENERATE
1,1601,200
*ELEMENT,TYPE=S4,ELSET=SHELL
1,1,3,203,201,
*ELGEN,ELSET=SHELL
1,8,2,2,8, 200,200
*SHELL SECTION,ELSET=SHELL,MAT=ELA
.4,3
*MATERIAL,NAME=ELA
*ELASTIC
6.825E7,.3
*EXPANSION
1.E-6,
*INITIAL CONDITIONS,TYPE=TEMP
ALL,0.0
*BOUNDARY
XYSYMM,3,5,
XZSYMM,2, ,
XZSYMM,4, ,
XZSYMM,6, ,
YZSYMM,1, ,
YZSYMM,5,6,
*STEP,PERTURBATION
*STATIC
*TEMPERATURE
ALL,430.
*ELPRINT
SINV,SP
*NODE PRINT,F=0
*NODE FILE
U, 
*OUTPUT,FIELD
*NODE OUTPUT
U,
*END STEP
*STEP,NLGEOM
***STATIC,DIRECT=NOSTOP
*STATIC
*TEMPERATURE
ALL,430.
*CONTROLS, PARAMETERS=FIELD
, , , 0.01
***CONTROLS, PARAMETERS=TIME INCREMENTATION
**,,,,,,,4
*ELPRINT
SINV,SP
*NODE PRINT,F=0
*NODE FILE
U, 
*OUTPUT,FIELD
*NODE OUTPUT
U,
*ENDSTEP