/batch,list
/verify,demonstration_problem3_289_290
JPGPRF,500,100,1
/title,Piping model meshed with PIPE289 and ELBOW290 elements
/filnam,289-290
/PREP7
YoungModulus1 = .258e+8 ! Young's Modulus
Nu = 0.3 ! Poissons ratio
ShearModulus1 = YoungModulus1/(2*(1+Nu)) ! Shear Modulus
WMass = 1.042868e-03 ! Density
WTick=0.216 ! Wall Thickness
OD=3.5 ! Outer Diameter
RADCUR=48.003 ! Radius curvature
temp=60 ! Temperature
maxm=15 ! Number of modes to extract
et,1,pipe289,,, ! Straight pipe elements
et,2,elbow290,,6 ! Curved pipe elements
et,3,combin14 ! Spring-damper elements
keyopt,3,2,1 ! UX Degree Of Freedom
et,4,combin14, ! Spring-damper elements
keyopt,4,2,2 ! UY Degree Of Freedom
et,5,combin14 ! Spring-damper elements
keyopt,5,2,3 ! UZ Degree Of Freedom
et,6,combin14 ! Spring-damper elements
keyopt,6,2,1 ! UX Degree Of Freedom
et,7,combin14 ! Spring-damper elements
keyopt,7,2,2 ! UY Degree Of Freedom
/com,
/com, Real Constants
/com,****************
sectype,1,PIPE,ctube ! Pipe Section Definition
secdata,OD,WTick,24
r,3,0.2e+8 ! Stiffness
r,4,0.2e+8 ! Stiffness
r,5,0.2e+8 ! Stiffness
r,6,0.2e+5 ! Stiffness
r,7,0.2e+5 ! Stiffness
/com,----------------------------------------------------------------------------
/com,
/com, Material Properties
/com,*********************
mp,ex,1,YoungModulus1
mp,nuxy,1,Nu
mp,gxy,1,ShearModulus1
mp,dens,1,WMass
mp,ex,2,YoungModulus1
mp,nuxy,2,Nu
mp,gxy,2,ShearModulus1
mp,dens,2,WMass
/com,----------------------------------------------------------------------------
/com, Keypoints
/com,***********
k,1,0,0,0
k,2,0,12,0
k,3,35.687,60,32.110
k,4,55,60,49.5
k,5,74.329,60,66.882
k,6,110,12,99
k,7,110,0,99
k,8,110,-24,99,
k,9,110,-48,99,
k,10,110,-72,99
k,11,110,-96,99
k,12,110,-120,99
k,13,110,-144,99
k,14,110,-168,99
k,15,110,-198,99
k,16,110,-228,99
k,17,110,-252,99
k,18,110,-276,99
k,19,110,-300,99
k,20,110,-324,99
k,21,99.6,-349.4,99
k,22,89.2,-374.8,99
k,23,78.8,-400,99
k,24,68.4,-425.6,99
k,25,58,-451,99
k,26,58,-475,99
k,27,58,-487,99
k,28,103.537,-535,114.179
k,29,124.269,-535,121.1
k,30,145,-535,128
k,31,184.975,-535,123.615
k,32,214.8,-536,102.8
k,33,254.585,-535,81.849
k,34,279.312,-535,75
k,35,331,-535,75
k,36,383,-535,75
/com,
/com, Elastic support Keypoints
/com,***************************
k,37,10,0,0
k,38,0,10,0
k,39,0,0,10
k,40,55,70,49.5
k,41,110,0,109
k,42,120,0,99
k,43,110,-168,109
k,44,120,-168,109
k,45,110,-324,109
k,46,120,-324,99
k,47,58,-475,109
k,48,68,-475,99
k,49,103.537,-545,114.179
k,50,103.537,-535,104.179
k,51,393,-535,75
k,52,383,-545,75
k,53,383,-535,85
/com,----------------------------------------------------------------------------
/com,
/com, Modeling of Straight Pipe (Tangent)
/com,************************************
l, 1, 2
l, 3, 4
l, 4, 5
l, 6, 7
l, 7, 8
l, 8, 9
l, 9,10
l,10,11
l,11,12
l,12,13
l,13,14
l,14,15
l,15,16
l,16,17
l,17,18
l,18,19
l,19,20
l,20,21
l,21,22
l,22,23
l,23,24
l,24,25
l,25,26
l,26,27
l,28,29
l,29,30
l,31,32
l,32,33
l,34,35
l,35,36 ! line number 30
/com,
/com, Modeling of Pipe Bend
/com,***********************
larch,2,3,4,RADCUR
larch,5,6,4,RADCUR
larch,27,28,26,RADCUR
larch,30,31,29,RADCUR
larch,33,34,32,RADCUR ! line number 35
/com, Elastic supports and anchors
/com,******************************
l,1,37
l,36,51
l,4,40
l,7,41
l,26,47
l,28,49
l,1,38
l,36,52
l,7,42
l,26,48
l,28,50
l,1,39
l,36,53
l,14,43
l,20,45
l,14,44
l,20,46
/com, **********************************
/com, Meshing for Straight pipe
/com, **********************************
type,1
secnum,1
mat,1
lsel,s,,,2,29
allsel,below,line
lesize,all,,,4
lmesh,all
allsel,all,all
/com, **********************************
/com, Meshing for bend pipe
/com, **********************************
type,2
secnum,1
mat,2
/com, *****************************************
/com, Note: The end elements are
/com, modeled using elbow elements
/com, in order to compare with the shell281
/com, results.
/com, With elbow290 in addition to constraining
/com, the nodal DOF you can also constrain the
/com, section constraints which is not possible
/com, with pipe289 elements.
lsel,s,line,,1
lsel,a,line,,30
lesize,all,,,6
lmesh,all
lsel,all
lsel,s,,,31,35
allsel,below,line
lesize,all,,,4
lmesh,all
allsel,all,all
/com, ************************************************************
/com, Using ELBOW, to convert some PIPE289 into ELBOW290 near
/com, the pipe ends
/com, ************************************************************
elbow,on,,,sect
esel,s,ename,,290
nsle,s
esln,s
nsle,s
esln,s
nsle,s
esln,s
esel,u,ename,,290
emodif,all,type,2
allsel,all
/com, ****************************
/com, Spring - damper elements
/com, ******************************
type,3
real,3
lsel,s,,,36,41
allsel,below,line
lesize,all,,,1
lmesh,all
allsel,all,all
type,4
real,4
lsel,s,,,42,46
allsel,below,line
lesize,all,,,1
lmesh,all
allsel,all,all
type,5
real,5
lsel,s,,,47,48
allsel,below,line
lesize,all,,,1
lmesh,all
allsel,all,all
type,6
real,6
lsel,s,,,49,50
allsel,below,line
lesize,all,,,1
lmesh,all
allsel,all,all
type,7
real,7
lsel,s,,,51,52
allsel,below,line
lesize,all,,,1
lmesh,all
allsel,all,all
n1 = node(55,60,49.500)
n2 = node(55,70,49.500)
n3 = node(35.687,60,32.110)
n4 = node(110,0,99.00)
n5 = node(110,0,109.00)
n6 = node(120,0,99.00)
n7 = node(110,-168.0,99.00)
n8 = node(110,-168.0,109.00)
n9 = node(120,-168.0,109.00)
n10 = node(110,-324.0,99.00)
n11 = node(110,-324,109.00)
n12 = node(120,-324,99.00)
n13 = node(58,-475,99.00)
n14 = node(58,-475,109.00)
n15 = node(68,-475,99)
n16 = node(103.54,-535,114.18)
n17 = node(103.54,-545,114.18)
n18 = node(103.54,-535,104.18)
n19 = node(10,0,0)
n20 = node(393,-535,75)
n21 = node(55,70,49.500)
n22 = node(110.0,0,109.0)
n23 = node(58.00,-475,109.00)
n24 = node(103.54,-545,114.18)
n25 = node(0,10,0)
n26 = node(0,0,0)
n27 = node(383,-535,75)
n28 = node(383,-545,75)
n29 = node(120,0,99)
n30 = node(68,-475,99)
n31 = node(103.54,-535,104.18)
n32 = node(0,0,10)
n33 = node(383.00,-535,85)
n34 = node(110,-168,109)
n35 = node(110,-324,109)
n36 = node(120,-168,109)
n37 = node(120,-324,99)
allsel,all
/com, rotate nodes with less than 3 supports
/com,
/com, rotate nodes with less than 3 supports
/com,
n1 = 2
n2 = 292
n3 = 1
ics = 11
wplane,,nx(n1),ny(n1),nz(n1),nx(n2),ny(n2),nz(n2),nx(n3),ny(n3),nz(n3)
cswplane,ics,0
nrotat,n1
nrotat,n2
csys,0
n1 = 19
n2 = 293
n3 = 298
ics = ics + 1
wplane,,nx(n1),ny(n1),nz(n1),nx(n2),ny(n2),nz(n2),nx(n3),ny(n3),nz(n3)
cswplane,ics,0
nrotat,n1
nrotat,n2
nrotat,n3
csys,0
n1 = 75
n2 = 303
n3 = 305
ics = ics + 1
wplane,,nx(n1),ny(n1),nz(n1),nx(n2),ny(n2),nz(n2),nx(n3),ny(n3),nz(n3)
cswplane,ics,0
nrotat,n1
nrotat,n2
nrotat,n3
csys,0
n1 = 123
n2 = 304
n3 = 306
ics = ics + 1
wplane,,nx(n1),ny(n1),nz(n1),nx(n2),ny(n2),nz(n2),nx(n3),ny(n3),nz(n3)
cswplane,ics,0
nrotat,n1
nrotat,n2
nrotat,n3
csys,0
n1 = 171
n2 = 294
n3 = 299
ics = ics + 1
wplane,,nx(n1),ny(n1),nz(n1),nx(n2),ny(n2),nz(n2),nx(n3),ny(n3),nz(n3)
cswplane,ics,0
nrotat,n1
nrotat,n2
nrotat,n3
csys,0
n1 = 187
n2 = 295
n3 = 300
ics = ics + 1
wplane,,nx(n1),ny(n1),nz(n1),nx(n2),ny(n2),nz(n2),nx(n3),ny(n3),nz(n3)
cswplane,ics,0
nrotat,n1
nrotat,n2
nrotat,n3
csys,0
allsel,all,all
/com,----------------------------------------------------------------------------
/com,
/com, Constraints
/com,*************
nsel,,node,,290,306
d,all,all
allsel
d,230,rotx,,,,,roty,rotz
d,243,rotx,,,,,roty,rotz
d,230,sect,0
d,243,sect,0
allsel,all,all
finish
/com,----------------------------------------------------------------------------
/com
/com,==============
/com, Modal solve
/com,==============
/com,
/solution
antype,modal ! Perform Modal Analysis
modopt,lanb,maxm
lumpm,on ! Use Lumped Mass Approximation
mxpand,maxm,,,yes ! Expand solution with Element Calculations ON
solve
finish
/post1
set,list ! Frequencies obtained from Modal solve
/show,jpeg
/eshape,5
/efacet,2
/view,1,1,2,3
/ang,1
/auto,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,spectrum ! Perform Spectrum Analysis
spopt,sprs,15 ! Single Point Excitation Response Spectrum
srss,0.0 ! SRSS mode combination
gval = 386.4
svtyp, 2, gval ! Seismic Acceleration Response Loading
freq, 2.5 , 5.0 , 8.0, 12.35, 13.51, 16.95, 20.83, 23.25, 29.41
sv,, 2.275, 2.275, 1.0, 0.8, 0.925, 0.925, 0.8 , 1.0 , 1.0
freq, 34.48
sv,, 0.875
sed,1,0,0 ! Excitation in X direction
SOLVE
svtyp, 2, gval ! Seismic Acceleration Response Loading
freq
freq, 2.5 , 5.0 , 8.0, 12.35, 13.51, 16.95, 20.83, 23.25, 29.41
sv,, 1.517, 1.517, 0.667, 0.534, 0.617, 0.617, 0.534, 0.667, 0.667
freq, 34.48
sv,, 0.584
sed,0,1,0 ! Excitation in Y direction
SOLVE
fini
/com,----------------------------------------------------------------------------
/post1
/input,289-290,mcom
/com, *========================================================
/com, * Displacement Sum and Equivalent Stress
/com, *==========================================================
/auto,1
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
