/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
nl=40		!number of line segs
 

/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.15
r,1,tw

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

k,1 
k,2,r 
l,1,2,nl
lmesh,1
dk,2,uy
dk,2,ux,istrn*r
fini
/solu
anty,stat
nlgeom,on
solve

sfl,all,pres,p
solve 
fini
/eof
/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, 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,S,X,0,1.0