/com	ANSYS 11.0
! Next100_field cage 
Vc=95000	!cathode voltage
dvdyd=50000	!gradient, drift
Ri=.535		!ring major radius
Ro=.57		!ground radius
rd=.005	!ring minor radius drift
rca=.0025	!ring minor radius cathode
sd=.02		!ring spacing, drift
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=10		!number of line segs
nr=60
nRk=20
nR=40
wa=ro-(ri+rd)	!annulus width
f=.5	!voltage add back in factor (max) for cathode 1= dvc,.5=0.5dvc, etc

nd=nint(ld/sd)
/com	nc=nint(((l+fk*rk)-ld)/sc)
nc=3	!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
/com	LARC, P1, P2, PC, RAD

larc,2,4,3,Rk,nRk
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

CIRCLE,11,rd
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,ld
k,99,Ri+.01,ld
k,100,ri,l+.01

/com	CSKP, KCN, KCS, PORIG, PXAXS, PXYPL, PAR1, PAR2
CSKP,11,cylin,98,99,100,(0.7*rk)/(ro-ri)
nums,kp,101
*do,i,1,nc-1
k,97+4*i,(i/nc)*wa+rd,0
k,98+4*i,(i/nc)*wa+rd,90
k,99+4*i,(i/nc)*wa+.001+rd,0
k,100+4*i,(i/nc)*wa+.001+rd,90
l,97+4*i,98+4*i,40
l,99+4*i,100+4*i,40
l,97+4*i,99+4*i,5
l,98+4*i,100+4*i,5
*enddo
csys,0

nums,area,101
*do,i,1,nc-1
al,97+4*i,98+4*i,99+4*i,100+4*i
*enddo
/com	agen,nc,101,,,,sc,,,1

/com	BLC4, XCORNER, YCORNER, WIDTH, HEIGHT, DEPTH
/com	RECTNG, X1, X2, Y1, Y2
nums,line,301
RECTNG,.001,ri-.001,ld+rd,(ld+rd)+.0001
nums,area,1000
asba,3,all


esiz,.0003
amesh,1000
lsel,s,line,,1,9
lsel,u,line,,4
lsel,u,line,,6,7
dl,all,1000,volt
dl,6,1000,volt,vi
dl,7,1000,volt,vi
dl,11,1000,volt,vi
dl,10,1000,symm
dl,4,1000,symm

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



*do,i,1,nc-1
lsel,s,line,,4*i+97,4*i+100
/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-(vc*(i/nc))
*enddo

/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
/com	lsel,s,line,,12,999
/com	dl,all,2,volt,v
alls



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