/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=.015		!ring spacing, drift
sc=.0075		!ring spacing, cathode buffer
Rk=.07		!knuckle radius
fk=.85		!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

nd=nint(ld/sd)
nc=nint(((l+fk*rk)-ld)/sc)

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
k,12,Ri,ld+(sc+rca)

/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
CIRCLE,12,rca
nums,area,101
al,101,102,103,104
agen,nc,101,,,-.002,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

/eof

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
lsel,s,line,,4*i+97,4*i+100
dl,all,1000,volt,vc-i*dvc
*enddo

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