/com	ANSYS 11.0
! mesh plane
! 1 layer 2d-axi shell208

p=12.3		!electric field pressure 
prestr=10000	!initial stress 
e_mesh=4e9	!mesh modulus


! t_mesh=.000025

r1=.070	!ring radius



/prep7
et,1,208	!2d-axi shell208 (non-harmonic)
! KEYOPT, ITYPE, KNUM, VALUE
keyopt,1,1,1	!membrane stiffnes only
KEYOPT, 1,3,0	!exclude extra internal node


mp,ex,1,e_mesh		
	

MP,PRXY,1,0.3


! SECTYPE, SECID, Type, Subtype, Name, REFINEKEY

SECTYPE,1,shell,,mesh
! secdata,TK, MAT, THETA, NUMPT, LayerName

secdata,t_mesh,1,0,3,mesh



k,1
k,2,r1

l,1,2,40



type,1
secnum,1
lmesh,1
dk,2,ux,0,,,uy

! inis,defi,,,,,prestr 	!apply radial prestress


/eshape,1

nummrg,node,1e-6
numc,all
fini
/solu
anty,stat
nlgeom,on
! SFL, LINE, Lab, VALI, VALJ, VAL2I, VAL2J
SFL,1,pres,p
/com	DK, KPOI, Lab, VALUE, VALUE2, KEXPND, Lab2, Lab3, Lab4, Lab5, Lab6




solve 
fini
/post1
set,last
! ETABLE, Lab, Item, Comp
/PSF,PRES,NORM,2,0,1
/PBF,DEFA, ,1 
/PBC,ALL,,1   


! PLNSOL,S,X,0,1.0
! /titl, 0.3mm 6H-SiC, 0.1mm epoxy, 0.2mm quartz, radial stress
! /CONT,1,9,-0.11e10, ,.037e9	!use to check quartz limits
! /CONT,1,9,AUTO  
! /titl, 0.371 mm 6H-SiCl stress
! PLNSOL,S,X,0,1.0