/com	ANSYS 11.0
! Next100_field cage with rectangular rings and cathode radius
Vc=95000	!cathode voltage
dvdyd=50000	!gradient, drift
Ri=.535		!ring major radius
Ro=.57		!ground radius
rd=.005		!ring minor radius drift
rw=.01		!ring width, radial direction
rhf=.9		!ring height factor, as percentage of spacing dimension
rca=.0025	!ring minor radius cathode
sd=.02		!ring spacing, drift
rh=sd*rhf	!ring height
sc=.01		!ring spacing, cathode buffer
Rk=.07		!knuckle radius
fk=.7		!knuckle factor
Rc=1.14		!crown radius
l=0.1		!straight section length
ld=0.1		!cathode distance up
nl=40		!number of line segs
nr=60
nRk=40
nR=40
wa=ro-(ri+rw)	!annulus width
f=.5	!voltage add back in factor (max) for cathode 1= dvc,.5=0.5dvc, etc
s2=1.7	!factor for second ellipsotoroid
nd=nint(ld/sd)
/com	nc=nint(((l+fk*rk)-ld)/sc)
nc=1	!number of cathode divisions (not graders)

dvdyc=vc/((l+rk)-ld)
dvd=dvdyd*sd
dvc=dvdyc*sc
vi=Vc-ld*dvdyd

/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,121	!2d 8-node  estat solid
keyopt,1,3,1	!axi symmetric
mp,perx,1,1

/com	define pv shell
k,1,Ro 
k,2,Ro,l
k,3,Ro-Rk,l
k,4,Ro-Rk,l+Rk
k,5,,l+Rk
k,6
l,1,2,nl,.33
/com	LARC, P1, P2, PC, RAD

larc,2,4,3,Rk,nRk,3
l,4,5,nR,4
l,5,6,nl
l,6,1,nR,.25

al,1,2,3,4,5
k,11,Ri

/com	CIRCLE, PCENT, RAD, PAXIS, PZERO, ARC, NSEG

/com	CIRCLE,11,rd
/com	define field cage rings, drift region
/com	BLC5, XCENTER, YCENTER, WIDTH, HEIGHT, DEPTH

blc5,ri,,rw,rh
al,6,7,8,9
/com	AGEN, ITIME, NA1, NA2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE
asba,1,2,,delete,keep
/com	LESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV
lesize,10,.001
lesize,11,,,nr,.25
nums,line,12
agen,nd+1,2,,,,sd,,,1
adel,2
nums,line,101

/com	instead of full toroids
/com	CIRCLE,12,rca
/com 	use quarter-ellipto-toroids
/com	k,12,Ri-(rd-rca),ld+(sc+rca)
/com 	x axis forlocal corrd sys
k,98,Ri-0.5*rw,ld+0.5*sd
k,99,Ri+.01,ld+0.5*sd
k,100,ri-0.5*rw,ld+sd

/com	build cathode arc toroid
/com	CSKP, KCN, KCS, PORIG, PXAXS, PXYPL, PAR1, PAR2
CSKP,11,cylin,98,99,100,1.15
nums,kp,101

k,101,rw,0
k,102,rw,170
k,103,-.0005+rw,0
k,104,-.0005+rw,170
l,101,102,60,4
l,103,104,60,4
l,101,103,3
l,102,104,3

/com	add a second shell at half potential with more ellipticity
CSKP,12,cylin,98,99,100,1.35

k,105,rw+0.5*wa,345
k,106,rw+0.5*wa,110
k,107,.0005+(rw+0.5*wa),345
k,108,.0005+(rw+0.5*wa),110
l,105,106,60,4
l,107,108,60,4
l,105,107,3
l,106,108,3

csys,0
nums,area,101

al,101,102,103,104
al,105,106,107,108

/com	agen,nc,101,,,,sc,,,1
/com	draw mesh plane
/com	BLC4, XCORNER, YCORNER, WIDTH, HEIGHT, DEPTH
/com	RECTNG, X1, X2, Y1, Y2

/com	build mesh plane
nums,line,301
RECTNG,.002,ri-(0.5*rw+.002),(ld+rw)-.001,(ld+rw)-.0005
nums,area,1000
asba,3,all


lsel,s,line,,1,3
dl,all,1000,volt
alls
dl,8,1000,volt,vi
lsel,s,line,,11,13
dl,all,1000,volt,vi
alls
dl,10,1000,symm
dl,4,1000,symm

*do,i,1,nd
lsel,s,line,,4*i+10,4*i+13
dl,all,1000,volt,vi+i*dvd
*enddo



lsel,s,line,,101,104
/com	divide voltage linearly then add back in a progressive voltage 
/com	dl,all,1000,volt,(vc-i*dvc)+((nc-i)/nc)*(f*dvc)
dl,all,1000,volt,vc

lsel,s,line,,105,108
dl,all,1000,volt,0.55*vc


/com	use a series of elliptic arc sections
/com	*do,i,1,nc



/com	use a series of elliptic arcs
lsel,s,line,,301,304
dl,all,1000,volt,vc

alls


esiz,.0003
amesh,1000


fini
/solu
anty,stat
solve



/post1
set,last
PLNSOL, VOLT,, 0
/eof
! 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