/com	ANSYS 11.0
! curved window study
!   2d-axi  
pi=acos(-1)
p=1.54e6		!pressure 
 
e_6h_sic=350.0e9	!moissanite
e_si2o3=73.0e9	!quartz
e_epoxy=3.0e9	!epoxy
e_sapphire=300e9
e_ti=100e9


! t_6h_sic=.000371
t_6h_sic=.0005
t_si2o3 =.00010
t_epoxy =.000050


r1=.040	!quartz or
r2=.04	!seal ir

rc1=.15		!radius of curvature,outer surface
t=.005		!center thickness
rc2=rc1-t	!radius of curvature,inner surface
tc=.003		!can thickness
lc=.04		!can length
lch=.005		!chamfer length
nch=nint(3000*lch)	!number of elements


ang1=asin(r1/rc1)*(180/pi)
ang2=asin(r2/rc2)*(180/pi)

/prep7
et,1,82	!2d solid
! KEYOPT, ITYPE, KNUM, VALUE
KEYOPT,1,3,1	!axisymmetric

mp,ex,1,e_6h_sic		
mp,ex,2,e_epoxy		
mp,ex,3,e_si2o3	
mp,ex,4,e_sapphire
mp,ex,5,e_ti	

MP,PRXY,1,0.3
MP,PRXY,2,0.3
MP,PRXY,3,0.3

k,1,0,rc1
k,2,0,(rc1)
k,3

circle,1,rc1,,3,ang1
circle,2,rc2,,3,ang1

lesi,all,,,100
l,4,6,6
l,5,7,6
al,1,2,3,4
mat,3

amesh,1
*get,x5,kp,5,loc,x
*get,y5,kp,5,loc,y
*get,x7,kp,7,loc,x
*get,y7,kp,7,loc,y

k,15,x5+tc,y5+lch
k,16,x5+tc,lc
k,17,x7,lc


l,5,15,nch
l,15,16,70
l,7,17,80
l,16,17,6
al,4,5,6,8,7
mat,5
esiz,.0005
amesh,2


/com	dl,4,1,ux
/com	dl,4,1,uy
dl,8,2,uy
sfl,1,pres,p
sfl,5,pres,p
sfl,6,pres,p


 

fini
/solu
anty,stat

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, 0.3mm 6H-SiC, 0.1mm epoxy, 0.2mm quartz, radial stress
!/CONT,1,9,-0.11e10, ,.037e9	!use to check quartz limits
! /CONT,1,9,AUTO  
! /titl, 0.371 mm 6H-SiCl stress
PLNSOL,S,1,0,1.0




/eof
da1z_n0Hig2