*heading
.inp:  hemispherical shell section under free thermal expansion [S4]
*RESTART,WRITE,FREQ=999
*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,offset=sneg
.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