/batch,list
/verify,demonstration_problem3_281
JPGPRF,500,100,1
/title,Piping model meshed with SHELL281 elements
/filname,281
/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
out_rad = OD/2 ! Outer Radius
in_rad = OD/2 - WTick ! Inner Radius
hf_thick = WTick/2 ! Half Thicknes
midd = in_rad + hf_thick
et,1,shell281 ! SHELL281 elements
sectype,1,shell ! Shell Section Definition
secdata,WTick,1,0,3
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
et,8,conta175 ! Contact175 element
keyopt,8,2,2 ! Multipoint Constraint
keyopt,8,4,2 ! Contact normal direction
keyopt,8,12,6 ! Bonded
et,9,targe170 ! Target170 element
keyopt,9,2,1 ! Boundary conditions for target nodes
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, Material Properties
/com,*********************
mp,ex,1,YoungModulus1
mp,nuxy,1,Nu
mp,gxy,1,ShearModulus1
mp,dens,1,WMass
/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, Forming the area
/com, *************************
k,101,midd,0,0
k,102,0,0,midd
k,103,-midd,0,0
k,104,0,0,-midd
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,2,2
adrag,1,2,3,4,,,5
larc,2,3,4,RADCUR
adrag,6,13,11,9,,,14
l,3,4
adrag,18,15,22,20,,,23
l,4,5
adrag,24,31,29,27,,,32
larc,5,6,4,RADCUR
adrag,33,40,38,36,,,41
l,6,7
adrag,42,45,47,49,,,50
l,7,8
adrag,51,54,56,58,,,59
l,8,9
adrag,60,63,65,67,,,68
l,9,10
adrag,69,72,74,76,,,77
l,10,11
adrag,78,81,83,85,,,86
l,11,12
adrag,87,90,92,94,,,95
l,12,13
adrag,96,99,101,103,,,104
l,13,14
adrag,105,108,110,112,,,113
l,14,15
adrag,114,117,119,121,,,122
l,15,16
adrag,123,126,128,130,,,131
/com, ***************************************
l,16,17
adrag,132,135,137,139,,,140
l,17,18
adrag,141,144,146,148,,,149
l,18,19
adrag,150,153,155,157,,,158
l,19,20
adrag,159,162,164,166,,,167
l,20,21
adrag,168,171,173,175,,,176
l,21,22
adrag,177,180,182,184,,,185
l,22,23
adrag,186,189,191,193,,,194
l,23,24
adrag,195,198,200,202,,,203
l,24,25
adrag,204,207,209,211,,,212
l,25,26
adrag,213,216,218,220,,,221
l,26,27
adrag,222,225,227,229,,,230
larch,27,28,26,RADCUR
adrag,231,234,236,238,,,239
l,28,29
adrag,240,243,245,247,,,248
l,29,30
adrag,249,252,254,256,,,257
larch,30,31,29,RADCUR
adrag,258,261,263,265,,,266
l,31,32
adrag,267,270,272,274,,,275
l,32,33
adrag,276,279,281,283,,,284
larch,33,34,32,RADCUR ! line number 35
adrag,285,288,290,292,,,293
l,34,35
adrag,294,297,299,301,,,302
l,35,36
adrag,303,306,308,310,,,311
allsel,all
lesize,all,,,10
/com,
/com, Meshing the area
/com, ***************************
type,1
mat,1
secnum,1
shpp,off
amap,1,54,55,101,102
amap,2,55,56,102,103
amap,3,56,57,103,104
amap,4,54,57,101,104
amap,5,54,55,58,59
amap,6,55,56,59,60
amap,7,56,57,60,61
amap,8,54,57,58,61
amap,9,59,60,62,63
amap,10,60,61,63,64
amap,11,58,61,64,65
amap,12,58,59,62,65
amap,13,62,63,66,67
amap,14,63,64,67,68
amap,15,64,65,68,69
amap,16,62,65,66,69
amap,17,66,67,70,71
amap,18,67,68,71,72
amap,19,68,69,72,73
amap,20,66,69,70,73
amap,21,70,71,74,75
amap,22,71,72,75,76
amap,23,72,73,76,77
amap,24,70,73,74,77
amap,25,74,75,78,79
amap,26,75,76,79,80
amap,27,76,77,80,81
amap,28,74,77,78,81
amap,29,78,79,82,83
amap,30,79,80,83,84
amap,31,80,81,84,85
amap,32,78,81,82,85
amap,33,82,83,86,87
amap,34,83,84,87,88
amap,35,84,85,88,89
amap,36,82,85,86,89
amap,37,86,87,90,91
amap,38,87,88,91,92
amap,39,88,89,92,93
amap,40,86,89,90,93
amap,41,90,91,94,95
amap,42,91,92,95,96
amap,43,92,93,96,97
amap,44,90,93,94,97
amap,45,94,95,98,99
amap,46,95,96,99,100
amap,47,96,97,100,105
amap,48,94,97,98,105
amap,49,98,99,106,107
amap,50,99,100,107,108
amap,51,100,105,108,109
amap,52,98,105,106,109
amap,53,106,107,110,111
amap,54,107,108,111,112
amap,55,108,109,112,113
amap,56,106,109,110,113
amap,57,110,111,114,115
amap,58,111,112,115,116
amap,59,112,113,116,117
amap,60,110,113,114,117
amap,61,114,115,118,119
amap,62,115,116,119,120
amap,63,116,117,120,121
amap,64,114,117,118,121
amap,65,118,119,122,123
amap,66,119,120,123,124
amap,67,120,121,124,125
amap,68,118,121,122,125
amap,69,122,123,126,127
amap,70,123,124,127,128
amap,71,124,125,128,129
amap,72,122,125,126,129
amap,73,126,127,130,131
amap,74,127,128,131,132
amap,75,128,129,132,133
amap,76,126,129,130,133
amap,77,130,131,134,135
amap,78,131,132,135,136
amap,79,132,133,136,137
amap,80,130,133,134,137
amap,81,134,135,138,139
amap,82,135,136,139,140
amap,83,136,137,140,141
amap,84,134,137,138,141
amap,85,138,139,142,143
amap,86,139,140,143,144
amap,87,140,141,144,145
amap,88,138,141,142,145
amap,89,142,143,146,147
amap,90,143,144,147,148
amap,91,144,145,148,149
amap,92,142,145,146,149
amap,93,146,147,150,151
amap,94,147,148,151,152
amap,95,148,149,152,153
amap,96,146,149,150,153
amap,97,150,151,154,155
amap,98,151,152,155,156
amap,99,152,153,156,157
amap,100,150,153,154,157
amap,101,154,155,158,159
amap,102,155,156,159,160
amap,103,156,157,160,161
amap,104,154,157,158,161
amap,105,158,159,162,163
amap,106,159,160,163,164
amap,107,160,161,164,165
amap,108,158,161,162,165
amap,109,162,163,166,167
amap,110,163,164,167,168
amap,111,164,165,168,169
amap,112,162,165,166,169
amap,113,166,167,170,171
amap,114,167,168,171,172
amap,115,168,169,172,173
amap,116,166,169,170,173
amap,117,170,171,174,175
amap,118,171,172,175,176
amap,119,172,173,176,177
amap,120,170,173,174,177
amap,121,174,175,178,179
amap,122,175,176,179,180
amap,123,176,177,180,181
amap,124,174,177,178,181
amap,125,178,179,182,183
amap,126,179,180,183,184
amap,127,180,181,184,185
amap,128,178,181,182,185
amap,129,182,183,186,187
amap,130,183,184,187,188
amap,131,184,185,188,189
amap,132,182,185,186,189
amap,133,186,187,190,191
amap,134,187,188,191,192
amap,135,188,189,192,193
amap,136,186,189,190,193
amap,137,190,191,194,195
amap,138,191,192,195,196
amap,139,192,193,196,197
amap,140,190,193,194,197
shpp,on
/com,
/com, Defining MPC contacts
/com, **********************
/com, support at node 1
tshap,pilot
n,1000000,0,0,0
type,9
real,9
e,1000000
type,8
real,9
lsel,s,line,,1,4,1
nsll,s,1
esln,s
esurf
allsel,all
/com, support at end node
tshap,pilot
n,1000001,383,-535,75
type,9
real,10
e,1000001
type,8
real,10
lsel,s,line,,312
lsel,a,line,,315
lsel,a,line,,317
lsel,a,line,,319
nsll,s,1
esln,s
esurf
allsel,all
/com, support at key point 7
tshap,pilot
n,1000002,110,0,99
type,9
real,11
e,1000002
type,8
real,11
lsel,s,line,,51
lsel,a,line,,54
lsel,a,line,,56
lsel,a,line,,58
nsll,s,1
esln,s
esurf
allsel,all
/com, support at key point 14
tshap,pilot
n,1000003,110,-168,99
type,9
real,12
e,1000003
type,8
real,12
lsel,s,line,,114
lsel,a,line,,117
lsel,a,line,,119
lsel,a,line,,121
nsll,s,1
esln,s
esurf
allsel,all
/com, support at key point 20
tshap,pilot
n,1000004,110,-324,99
type,9
real,13
e,1000004
type,8
real,13
lsel,s,line,,168
lsel,a,line,,171
lsel,a,line,,173
lsel,a,line,,175
nsll,s,1
esln,s
esurf
allsel,all
/com, support at key point 26
tshap,pilot
n,1000005,58,-475,99
type,9
real,14
e,1000005
type,8
real,14
lsel,s,line,,222
lsel,a,line,,225
lsel,a,line,,227
lsel,a,line,,229
nsll,s,1
esln,s
esurf
allsel,all
/com, support at key point 28
tshap,pilot
n,1000006,103.537,-535,114.179
type,9
real,15
e,1000006
type,8
real,15
lsel,s,line,,240
lsel,a,line,,243
lsel,a,line,,245
lsel,a,line,,247
nsll,s,1
esln,s
esurf
allsel,all
/com, support at key point 4
tshap,pilot
n,1000007,55,60,49.5
type,9
real,16
e,1000007
type,8
real,16
lsel,s,line,,24
lsel,a,line,,31
lsel,a,line,,29
lsel,a,line,,27
nsll,s,1
esln,s
esurf
allsel,all
/com, Nodes for elastic support
/com, ******************************
n,2000000,2,0,0
n,2000001,0,2,0
n,2000002,0,0,2
n,2000003,55,62,49.5
n,2000004,110,0,101
n,2000005,112,0,99
n,2000006,110,-168,101
n,2000007,112,-168,101
n,2000008,110,-324,101
n,2000009,112,-324,99
n,2000010,58,-475,101
n,2000011,60,-475,99
n,2000012,103.537,-537,114.179
n,2000013,103.537,-535,112.179
n,2000014,385,-535,75
n,2000015,383,-537,75
n,2000016,383,-535,77
/com, ****************************
/com, Spring - damper elements
/com, ******************************
type,3
real,3
e,1000000,2000000
e,1000001,2000014
e,1000007,2000003
e,1000002,2000004
e,1000005,2000010
e,1000006,2000012
allsel,all,all
type,4
real,4
e,1000000,2000001
e,1000001,2000015
e,1000002,2000005
e,1000005,2000011
e,1000006,2000013
allsel,all,all
type,5
real,5
e,1000000,2000002
e,1000001,2000016
allsel,all,all
type,6
real,6
e,1000003,2000006
e,1000004,2000008
allsel,all,all
type,7
real,7
e,1000003,2000007
e,1000004,2000009
allsel,all,all
/com, rotate nodes with less than 3 supports
/com,
n,3000000,35.687,60,32.110
n1 = 1000007
n2 = 2000003
n3 = 3000000
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 = 1000002
n2 = 2000004
n3 = 2000005
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 = 1000003
n2 = 2000006
n3 = 2000007
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 = 1000004
n2 = 2000008
n3 = 2000009
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 = 1000005
n2 = 2000010
n3 = 2000011
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 = 1000006
n2 = 2000012
n3 = 2000013
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, Constraints
/com,*************
nsel,,node,,2000000,2000016
d,all,all
allsel
nsel,s,node,,1000000
d,all,rotx,0
d,all,roty,0
d,all,rotz,0
allsel,all
nsel,s,node,,1000001
d,all,rotx,0
d,all,roty,0
d,all,rotz,0
allsel,all
finish
/com,----------------------------------------------------------------------------
/com
/com,==============
/com, Modal solve
/com,==============
/com,
/solution
antype,modal ! Perform ModalAnalysis
modopt,lanb,maxm ! Use LANB eigensolver
lumpm,on ! Use Lumped Mass Approximation
mxpand,maxm,,,yes ! Expand all modes
solve
finish
/post1
set,list ! List Frequencies obtained from Modal Solve
/show,jpeg
/eshape,0
/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 Response Spectrum Analysis
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
finish
/com,----------------------------------------------------------------------------
/post1
/input,281,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
