/batch,list
/verify,demonstration_problem2_281
JPGPRF,500,100,1
/title,Piping Model meshed with SHELL281 elements
/filnam,281
/prep7
fact = 1.5
out_dia = 7.288*fact ! Outer Diameter
wall_thk = 0.241 ! Wall Thickness
out_rad = out_dia/2 ! Outer radius
in_rad = out_dia/2 - wall_thk ! Inner radius
hf_thick = wall_thk/2 ! Half Thickness
midd = in_rad+hf_thick
et,1,shell281 ! SHELL281 elements
sectype,1,shell ! Section definition
secdata,wall_thk,1,0,3
mp,ex,1,24e6
mp,nuxy,1,0.3
mp,dens,1,0.000125
/com, keypoints
/com,*************
k,1,0.0,0.0,0.0
k,2,0.0,54.45,0.0
k,3,0.0,108.9,0.0
k,4,10.632,134.568,0.0
k,5,36.3,145.2,0.0
k,6,54.15,145.2,0.0
k,7,72.0,145.2,0.0
k,8,97.668,145.2,10.632
k,9,108.3,145.2,36.3
k,10,108.3,145.2,56.80
k,11,108.3,145.2,77.3
k,12,2.7631,122.79,0
k,13,22.408,142.44,0
k,14,85.9,145,2.76
k,15,106,145,22.4
k,101,midd,0,0
k,102,0,0,midd
k,103,-midd,0,0
k,104,0,0,-midd
/com, lines
/com,*************
larc,101,102,1,midd
larc,102,103,1,midd
larc,103,104,1,midd
larc,104,101,1,midd
lsel,s,line,,1,4,1
lesize,all,,,12
lsel,all
l,1,3,12
adrag,1,2,3,4,,,5
larc,3,4,12,-36.3
adrag,6,13,11,9,,,14
larc,4,5,13,-36.3,
adrag,18,15,22,20,,,23
l,5,7,12
adrag,31,27,29,24,,,32
larc,7,8,14,-36.3
adrag,33,40,38,36,,,41
larc,8,9,15,-36.3
adrag,42,49,47,45,,,50
l,9,11,12
adrag,51,58,56,54,,,59
lesize,all,,,12
allsel,all
amap,1,16,17,101,102
amap,2,17,18,102,103
amap,3,18,19,103,104
amap,4,16,19,101,104
amap,5,16,17,20,21
amap,6,17,18,21,22
amap,7,18,19,22,23
amap,8,16,19,20,23
amap,9,21,22,24,25
amap,10,22,23,25,26
amap,11,20,23,26,27
amap,12,20,21,24,27
amap,13,24,27,28,29
amap,14,24,25,29,30
amap,15,25,26,30,31
amap,16,26,27,28,31
amap,17,28,29,32,33
amap,18,29,30,33,34
amap,19,30,31,34,35
amap,20,28,31,32,35
amap,21,32,33,36,37
amap,22,33,34,37,38
amap,23,34,35,38,39
amap,24,32,35,36,39
amap,25,36,37,40,41
amap,26,37,38,41,42
amap,27,38,39,42,43
amap,28,36,39,40,43
/com, constraints
/com,*************
nsel,s,loc,y,0
d,all,all,0
alls
lsel,s,line,,60
lsel,a,line,,63
lsel,a,line,,65
lsel,a,line,,67
nsll,s,1
d,all,all,0
alls
finish
/solution
antype,modal
outres,all,all
modopt,lanb,15 ! LANB mode extraction method
mxpand,15,,,yes ! Expand all 15 modes
solve
fini
/post1
set,list ! Frequencies obtained from Modal Solve
rsys,solu
/show,jpeg
/eshape,0
/efacet,2
/view,1,1,1,1
/graphics,power
eplot
/replot
set,1,1
plnsol,u,sum
set,1,2
plnsol,u,sum
set,1,3
plnsol,u,sum
set,1,4
plnsol,u,sum
set,1,5
plnsol,u,sum
set,1,6
plnsol,u,sum
set,1,7
plnsol,u,sum
set,1,8
plnsol,u,sum
set,1,9
plnsol,u,sum
set,1,10
plnsol,u,sum
set,1,11
plnsol,u,sum
set,1,12
plnsol,u,sum
set,1,13
plnsol,u,sum
set,1,14
plnsol,u,sum
set,1,15
plnsol,u,sum
finish
/com,----------------------------------------------------------------------------
/com,
/com,================
/com, Spectrum Solve
/com,================
/com,
/solution
antype,spectr ! Perform Spectrum Analysis
spopt,sprs ! Single Point Excitation Response Spectrum
dmprat,0.02 ! Constant Damping Ratio
grp,0.0 ! Group Modes based on significance level
svtyp,2 ! Seismic Acceleration Response Loading
sed,1 ! Excitation in X direction
freq
freq,3.1,4,5,5.81,7.1,8.77,10.99,14.08,17.24
freq,25,28.5,30,34.97,55,80,140,162,588.93
sv,0.02,400,871,871,700,1188,1188,440,775,775
sv,0.02,533.2,467.2,443.6,380,289,239.4,192.6,184.1,145
solve
sed,,1 ! Excitation in Y direction
freq
freq,3.1,4,5,5.81,7.1,8.77,10.99,14.08,17.24
freq,25,28.5,30,34.97,55,80,140,162,588.93
sv,0.02,266.7,580.7,580.7,466.7,792,792,293.3,516.7,516.7
sv,0.02,355.5,311.5,295.7,253.3,192.7,159.6,128.4,122.7,96.7
solve
sed,,,1 ! Excitation in Z direction
freq
freq,3.1,4,5,5.81,7.1,8.77,10.99,14.08,17.24
freq,25,28.5,30,34.97,55,80,140,162,588.93
sv,0.02,400,871,871,700,1188,1188,440,775,775
sv,0.02,533.2,467.2,443.6,380,289,239.4,192.6,184.1,145
solve
fini
/com,----------------------------------------------------------------------------
/post1
/input,281,mcom
/com, *========================================================
/com, * Displacement Sum and Equivalent Stress
/com, *==========================================================
plnsol,u,sum
*get,umax,plnsol,0,max
*get,umin,plnsol,0,min
plnsol,s,eqv
*get,smax,plnsol,0,max
*get,smin,plnsol,0,min
/show,close
*stat,umax
*stat,umin
*stat,smax
*stat,smin
finish
/exit,nosave
