/com	ANSYS 11.0
! Next100_mesh plane deflection
p=1.54e2	!pressure from electric field
e_ss=19.0e9	!ti modulus (CP)
tw=.0002286	!wire thickness
r=.5		!radius
istrn=.002	!initial strain, axial
h=.002		!height of dome
n=20		!number of line segs
a=h/r**2	!slope for y=ax^2
 

/PNUM,KP,0  
/PNUM,LINE,1
/PNUM,AREA,0
/PNUM,VOLU,0
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0 
/psym,cs,1
/PSF,PRES,NORM,2,0,1
/SHRINK,0   
/ESHAPE,1   
/EFACET,1   
/RATIO,1,1,1

/prep7
et,1,51	!  2d   axi shell


mp,ex,1,e_ss		
MP,PRXY,1,0.3
r,1,tw

! SECTYPE, SECID, Type, Subtype, Name, REFINEKEY
*do,i,0,n
k,i+1,r*(i/n),a*(r*(i/n))**2
*enddo
*do,i,1,n
l,i,i+1,2
*enddo 

lmesh,all

dk,n+1,ux,istrn*r
dk,n+1,uy
dk,n+1,rotz
/com	dl,n,ux,istrn*r

fini
/solu
anty,stat
nlgeom,on
autots,on
solve
sfl,all,pres,p
solve

/com	SFA, AREA, LKEY, Lab, VALUE, VALUE2
/com	SFL, LINE, Lab, VALI, VALJ, VAL2I, VAL2J

solve 
fini

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

/DSCALE,1,100
pldi
! PLNSOL,S,X,0,1.0
! /titl, NEXT100 press vess 15 bar
! /CONT,1,9,-0.11e10, ,.037e9	!use to check quartz limits
! /CONT,1,9,AUTO  
! /titl, 0.371 mm 6H-SiCl stress
PLNSOL,u,sum,0,1.0