*heading
.inp:  hemispherical shell section under free thermal expansion [STRI65]
*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,1
*nset,nset=xzsymm,generate
1,17,1
*nset,nset=yzsymm,generate
1,1601,100
*element,type=stri65,elset=shell
1,1,3 ,201,    2, 102, 101
2,201,3,203,   102,103,202
*elgen,elset=shell
1,8,2,2,8, 200,200
2,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
*controls,param=field,field=displacement
2.6e-1,
*controls,param=field,field=rotation
2.6e-2,
*temperature
all,430.
*elprint
sinv,sp
*node print,f=0
*node file
u, 
*OUTPUT,FIELD
*NODE OUTPUT
u,
*end step