/com	ANSYS 11.0
! Next100_press vessel design study 

p=1.5e6		!pressure 
eTi=105.0e9	!ti modulus (CP)
tp=.13		!service plate thickness
tv=.01		!vessel wall thickness
tfh=.1		!flange thick,  vessel
tfl=.07 	!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=20		!number ofline segments at contact surfaces

rb=0.5*db	!bolt radius
rv=rvi+tv	!nom radius of vessel
pb=2*(rv/db)**2	!bolt pressure
/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,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
secoff,mid	! set section offset to bottom 

k,1,ri
k,11,rbl-rb
k,12,rbl+rb
k,13,rv
k,2,ro
k,3,ro,-tp
k,4,rbh+rb,-tp
k,5,rbh-rb,-tp
k,6,ri+wfl,-tp
k,7,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


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

l,7,8
l,8,81
l,81,82
l,82,9
l,9,6

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

l,13,14
l,14,15
l,15,16
l,16,17
l,17,2


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


l,18,14
l,8,21


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

/com	add second lines at contact surfaces ( should be l25,l26)
k,113,rv
k,102,ro
k,107,ri,-tp
k,106,ri+wfl,-tp
l,102,113,ne
l,107,106,ne

al,1,2,3,4,5,6,7,8,9,10
al,11,26,15,14,13,12
al,16,17,18,19,20,25
/eof