/com  quartz tube external pressure buckling (from STAR/al_150_long.1.txt) 
lth1=.003	!tube thickness
r=.021		!tube radius, nom
l=1.5		!tube length
ndiv_long=100	!number of longitudinal elements
ndiv_theta=60	!number of elements in circumference
nl= 1	!number of layers
/prep7
/view,1,.5,.5,1
/pbc,all,1
/title, [NEXT] quartz  tube ext. press. e_v_buckle, uz=0 100 kPA unit press.
/com    buckling analysis for STAR beampipe, from external pressure loading
/com 	code written for ansys 5.7

/com ELEMENT TYPES (nl layer shell elements)
et,1,91,nl,1,,,,1,,0,1 		!layered shell, 8 layers,print ilshearstress, store top/bot layers, sandwich on
/com          Real Set 1  al/pitch carbon composite
/com             layer thicknesses

mp,ex,1,70e9	!quartz
mp,prxy,1,0.334

r,1,nl,0,,,,0
rmore,,,,,,,
rmore,2,,lth1	!bottom aluminum(2) or beryllium(4) layer
/com	rmore,2,30,lth2	!tridirectional quasi-isotropic faceskins
/com	rmore,2,-30,lth2
/com	rmore,2,0,lth2
/com	rmore,3,0,lth3	!syntactic material
/com	rmore,2,0,lth2
/com	rmore,2,-30,lth2
/com	rmore,2,30,lth2
/com               COORDINATE SYSTEM DEFINITIONS
local,20,1,0,0,0
csys,20
/com               CYLINDER BODY DEFINITION
type,1
real,1
k,1
k,2,r
k,3,r,,l
k,4,,,l
l,2,3,ndiv_long
arotate,1,,,,,,1,4,360,4
lsel,s,line,,5,12
lesize,all,,3
amesh,1,4
nrotate,all
nsel,s,loc,z,l,l+.000001
d,all,ux,,,,,uy
cp,1,uz,all
nsel,s,loc,z
d,all,ux,,,,,uy,uz
nall
save
sfa,all,2,pres,100000	!unit pressure of 100 kPa
fini
/solu
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
/eof
/prep7
upgeom,.0001,,,al_150_long1,rst
fini
/solu
ANTYPE,0
nlgeom,on
pstress,on
sstif,on
solcontrol,on
nropt,unsym,,off
AUTOTS,on
LNSRCH,0
NEQIT,5   
PRED,0  
TIME,10
DELTIME,0.5,0.05,0.5,on 
/GO
kbc,0
SFA,all,2,PRES,1000000 
OUTPR,all,none   
OUTRES,ALL,ALL 
OUTPR,STAT  
OUTRES,STAT 
/STATUS,SOLU
SOLVE   
fini
/eof
/post26
STORE,NEW,all,
FILE,star_ripple_pipe,rst
TIMERANGE,1,50
NSTORE,1   
NSOL,2,5548,U,X,rad_disp  
VARNAM,uradial,rad_disp 
PLTIME,1,50 
XVAR,0  
SPREAD,0
PLCPLX,0
PLVAR,rad_disp