/com	ANSYS 11.0
! Next100_mesh plane deflection
p=1.54e2	!pressure from electric field
e_ss=24.0e9	!ti modulus (CP)
tw=.0005		!wire thickness
r=.5		!radius
istrn=.002	!initial strain, axial
h=.02		!height of dome
n=20		!number of line segs
a=h/r**2	!slope for y=ax^2
 

/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,63	!  3d  shell


mp,ex,1,e_ss		
MP,PRXY,1,0.3
r,1,tw

! SECTYPE, SECID, Type, Subtype, Name, REFINEKEY
*do,i,0,n
k,i+1,r*(i/n),a*(r*(i/n))**2
*enddo
*do,i,1,n
l,i,i+1,2
*enddo 
k,99,,0.1
/com	LESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV

/com	AROTAT, NL1, NL2, NL3, NL4, NL5, NL6, PAX1, PAX2, ARC, NSEG
arotate,all,,,,,,1,99
lesi,all,,5
amesh,all

csys,1
nrot,all
/com	DL, LINE, AREA, Lab, Value1, Value2
dl,100,20,ux
dl,120,40,ux
dl,140,60,ux
dl,160,80,ux
dl,100,20,uy
dl,120,40,uy
dl,140,60,uy
dl,160,80,uy
dl,100,20,uz
dl,120,40,uz
dl,140,60,uz
dl,160,80,uz
/com	dl,n,ux,istrn*r
csys,0
sfa,all,2,pres,1
fini
/solu
/com	anty,stat
/com	nlgeom,on
/com	autots,on
/com	solve
/com	SFA, AREA, LKEY, Lab, VALUE, VALUE2

anty,0
eqslv,sparse
pstress,on
solve
save
fini
/solu

antype,buckle
bucopt,lanb,1	!block lanczos option extract to mode 1
mxpand,1,,,,yes
outpr/nsol,all
save
solve
fini
/solu
expass,on
mxpand,1
outres,all
solve
save
fini
solve 
fini

/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,u,sum,0,1.0