/com	ANSYS 11.0
! Next100_press vessel design study 

p=1.5e6		!pressure 
eTi=105.0e9	!ti modulus (CP)
tp=.07		!service plate thickness
tv=.01		!vessel wall thickness
tfh=.05		!flange thick,  vessel
tfl=.04 	!flange thick. lower lid
rvi=.665	!vessel inner radius
ri=.46		!iR
ro=.75		!outer vesselr
wflh=.1		!vessel flange width
wfl=.05		!lower flange width
rbh=.71		!bolt circle radius upper
rbl=.48		!bolt circle lower
db=.015		!bolt diameter
lh=1		!cyl length, higher
ll=.2		!cyl length, lower
ne=30		!number ofline segments at contact surfaces
olh=.01		!overlap at contact plane, flange over plat
oll=.01	
rb=0.5*db	!bolt radius
rv=rvi+tv	!nom radius of vessel
pb=2*(rv/db)**2	!bolt pressure
/PNUM,KP,1 
/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,182	!higher order 2d plane
et,2,208	!axi 2 node finite strain shell element
/com	et,3,12		!contact, point to point, with friction
/com	KEYOPT, 2,3,2	!include extra internal node
/com	KEYOPT, 2,8,2	! store top, mid and bottom of each layer

mp,ex,1,eTi		
MP,PRXY,1,0.3

! SECTYPE, SECID, Type, Subtype, Name, REFINEKEY
SECTYPE,1,shell,,pv
! secdata,TK, MAT, THETA, NUMPT, LayerName
secdata,tv,1
/com	secoff,top	! set section offset to bottom 

k,1,ri
k,11,rbl-rb
k,12,rbl+rb
k,131,rv
k,13,rv
k,2,ro
k,201,ro
k,3,ro,-tp
k,4,rbh+rb,-tp
k,5,rbh-rb,-tp
k,61,ri+wfl,-tp
k,6,ri+wfl,-tp
k,7,ri,-tp
k,71,ri,-tp



l,1,11
l,11,12
l,12,13
l,13,2,ne
l,2,3
l,3,4
l,4,5
l,5,6
l,6,7,ne
l,7,1
al,1,2,3,4,5,6,7,8,9,10

k,8,ri,-(tp+tfl)
k,81,rbl-rb,-(tp+tfl)
k,82,rbl+rb,-(tp+tfl)
k,9,ri+wfl,-(tp+tfl)

l,71,8
l,8,81
l,81,82
l,82,9
l,9,61
l,61,71,ne

al,11,16,15,14,13,12

k,14,rv,tfh
k,15,rbh-rb,tfh
k,16,rbh+rb,tfh
k,17,ro,tfh

l,131,14
l,14,15
l,15,16
l,16,17
l,17,201
l,201,131,ne
al,17,18,19,20,21,22


k,18,rv,lh
k,19,0,lh+0.5*rv
k,20,0,lh
k,21,ri,-(tp+(tfl+ll))
k,22,0,-(tp+(tfl+ll))
k,23,0,-(tp+((tfl+ll)+0.5*ri))


l,18,14
l,8,21


/com	CSKP, KCN, KCS, PORIG, PXAXS, PXYPL, PAR1, PAR2
cskp,11,1,20,18,19,0.5
/com	circ,pcent,rqd,paxis,pzero,arc,nseg
/com	circ,20,rv,,18,90,1
l,19,18
cskp,12,1,22,23,21,2
/com	circ,pcent,rqd,paxis,pzero,arc,nseg
/com	circ,22,ri,,23,90,1
l,21,23
csys,0

/com	asbl,na,nl,--,keepa,keepl



type,1
mat,1
esiz,.002
amesh,1,3
type,2
esiz,.02
lmesh,11
lmesh,23,24
/com	secoff,bot
lmesh,17
lmesh,25,26


/com	esiz,.003
/com	type,3
/com	lmesh,4
/com	lmesh,9
/com	lsel,type,item,comp,vmin,vmax,vinc,kswp

lsel,s,line,,1,3
lsel,a,line,,10,11
lsel,a,line,,17
sfl,all,pres,p


lsel,s,line,,23,26
esll
sfe,all,2,pres,1,p


lsel,s,line,,7
lsel,a,line,,2
lsel,a,,,13
lsel,a,,,19
sfl,all,pres,pb

alls

nummrg,node,1e-6
numc,all
fini
/solu
anty,stat
! SFL, LINE, Lab, VALI, VALJ, VAL2I, VAL2J

dk,2,uy,0

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,S,X,0,1.0