/com	ANSYS 11.0
! Next100_field cage with rectangular radiused, thin rings and cathode radius
p=15		!pressure, bar
Vc=95000	!cathode voltage
ed=50000	!gradient, drift, v/m
eopa=2e5	!max E over p annulus buffer region v/m*bar
eam=eopa*p	!emax annulus
eopd=0.5e5	!max E over p between drift rings v/m*bar 
edm=eopd*p	!emax drift (between rings)

Ri=.535		!ring major radius
Ro=.57		!ground radius
rd=.005		!ring minor radius, drift
rw=.002		!ring width, radial direction <<rh for rect radiused rings
rca=.0025	!ring minor radius cathode
sd=.02		!ring spacing, drift
rhf=1-((ed*sd)/ed)	!ring height factor,as fraction of spacing dimension
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)

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

/PNUM,KP,1  
/PNUM,LINE,1
/PNUM,AREA,1
/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


/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

/com	blc5,ri,,rw,rh
/com make keypoints ccw from upper right
/com	k,11,Ri+0.5*rw,rh-0.5*rw
/com	k,12,Ri,rh-0.5*rw
/com	k,13,Ri-0.5*rw,rh-0.5*rw
/com	k,14,Ri-0.5*rw,-(rh-0.5*rw)
/com	k,15,Ri,-(rh-0.5*rw)
/com	k,16,Ri+0.5*rw,-(rh-0.5*rw)
/com	larc,11,13,15,0.5*rw
/com	l,13,14,10
/com	larc,14,16,12,0.5*rw
/com	l,16,11,10
numstr,kp,10
k,11,ri,0.5*(rh-rw)
k,12,ri+rw,0.5*(rh-rw)
k,13,ri,-0.5*(rh-rw)
k,14,ri-rw,-0.5*(rh-rw)
circ,11,0.5*rw,,12,180
circ,13,0.5*rw,,14,180
l,16,17
l,19,10

al,6,7,10,8,9,11
/com	LESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV
lsel,s,line,,10,11
lesize,all,,,15,-10
lsel,s,line,,6,9	!select ring ard sections
lesize,all,,5		!use 5 deg arcs per quadrant
alls

/com	AGEN, ITIME, NA1, NA2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE
asba,1,2,,delete,keep
lesize,12,,,20

lesize,13,,,nr,.33
lesize,14,,,8,3
lesize,15,,,8,3

nums,line,30
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*sd,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,-.00025+0.5*(sd+rw),1
k,102,-.0005+0.5*(sd+rw),1
k,103,-.00025+0.5*(sd+rw),170
k,104,0.5*(sd+rw),170
circ,101,.00025,,102,180
circ,103,.00025,,104,180

l,107,108,60,-10
l,110,105,60,-10
lsel,s,line,,101,104
lesize,all,,15
alls

add a second shell at half potential with more ellipticity
CSKP,12,cylin,98,99,100,1.2
k,111,-.00025+0.5*(sd+rw+wa),345
k,112,-.0005+0.5*(sd+rw+wa),345
k,113,-.00025+0.5*(sd+rw+wa),105
k,114,0.5*(sd+rw+wa),105
circ,111,.00025,,112,180
circ,113,.00025,,114,180

l,117,118,120,-10
l,120,115,120,-10
lsel,s,line,,107,110
lesize,all,,15
alls


csys,0
nums,area,101

al,101,102,103,104,105,106
al,107,108,109,110,111,112



/com	build mesh plane
nums,line,301
nums,kp,301
/com	RECTNG,.002,ri-(0.5*rw+.002),(ld+0.5*sd)-.001,ld+0.5*sd
k,301,.005,(ld+0.5*sd)-.001
k,302,.005,ld+0.5*sd
k,303,ri-(0.5*rw+.002),(ld+0.5*sd)-.001
k,304,ri-(0.5*rw+.002),ld+0.5*sd
l,301,302,3
l,302,304,30,-10
l,304,303,3
l,303,301,30,-10
al,301,302,303,304

nums,area,1000
asba,3,all


lsel,s,line,,1,3
dl,all,1000,volt
alls
dl,13,1000,volt,vi
lsel,s,line,,6,7
lsel,a,line,,14,15
dl,all,1000,volt,vi
alls
dl,12,1000,symm
dl,4,1000,symm

*do,i,1,nd
lsel,s,line,,6*i+24,6*i+29
dl,all,1000,volt,vi+i*dvd
*enddo



lsel,s,line,,101,106

dl,all,1000,volt,vc

lsel,s,line,,107,112
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,.05
amesh,1000
/eof

fini
/solu
anty,stat
solve



/post1
set,last
PLNSOL, VOLT,, 0
/EFACET,1     
/ZOOM,1,SCRN,1.192050,0.143972,1.348472,0.057092
PLNSOL, EF,SUM, 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