vm-nr1677-02-3c Input Listing


/COM,ANSYS MEDIA REL. 2024R2 (05/10/2024) REF. VERIF. MANUAL: REL. 2024R2
/verify,vm-nr1677-02-3c
/title,vm-nr1677-02-3c,NRC piping benchmarks problems,Volume II,Problem 3c

/com, *******************************************************************************
/com, Reference: Piping benchmark problems,Dynamic analysis independant support
/com,            motion response spectrum method, P. Bezler, M. Subudhi and
/com,            M.Hartzman, NUREG/CR--1677-Vol.2, August 1985.
/com, 
/com, Elements used: PIPE289, ELBOW290, COMBIN14
/com, 
/com, Results:
/com, The following results are outputted
/com, 1. Frequencies obtained from modal solution.
/com, 2. Maximum nodal displacements and rotations obtained from spectrum solution.
/com, 3. Element forces/moments obtained from spectrum solution.
/com, 4. Reaction forces obtained from spectrum solution.
/com,
/com, *******************************************************************************

/out,scratch

/prep7

YoungModulus1 = 0.277e+08			! Young's Modulus
Nu = 0.3						! Minor Poisson's Ratio
ShearModulus1 = YoungModulus1/(2*(1+Nu))	! Shear Modulus
WMass = 1.546e-03					! Density
WTick=0.3750					! Wall Thickness
OD=12.750						! Outer Diameter

Temperature = 400
Pressue = 615
maxm=15						! No. of Modes to Extract

radcur1 = -60              ! radius of curvature for bend pipes
radcur2 = -18              ! radius of curvature for bend pipes             


/com,------------------------------------------------------------------------------------

et, 1,pipe289,,,,2				! Element 1 - PIPE289
et, 2,elbow290,,6				! Element 2 - ELBOW290
et, 3,elbow290,,6				! Element 3 - ELBOW290

et, 4,combin14					! Element 4 - COMBIN14
keyopt,4,2,2					! Y Degree of Freedom
et, 5,combin14					! Element 5 - COMBIN14
keyopt,5,2,1					! X Degree of Freedom
et, 6,combin14					! Element 6 - COMBIN14
keyopt,6,2,2					! Y Degree of Freedom
et, 7,combin14					! Element 7 - COMBIN14
keyopt,7,2,3					! Z Degree of Freedom
et, 8,combin14					! Element 8 - COMBIN14
keyopt,8,2,4					! ROT-X Degree of Freedom
et, 9,combin14					! Element 9 - COMBIN14
keyopt,9,2,5					! ROT-Y Degree of Freedom
et,10,combin14					! Element 10 - COMBIN14
keyopt,10,2,6					! ROT-Z Degree of Freedom

/com,------------------------------------------------------------------------------------

/com, Real Constants
/com,****************

sectype,1,pipe
secdata,OD,WTick,24

r, 4, 0.1e+2
r, 5, 0.1e+13
r, 6, 0.1e+13
r, 7, 0.1e+13
r, 8, 0.1e+13
r, 9, 0.1e+13
r,10, 0.1e+13

/com, -------------------------------------------------------------------------------------

/com, Material Properties
/com,*********************

mp,ex,  1, YoungModulus1			
mp,nuxy,1, Nu					
mp,gxy ,1, ShearModulus1			
mp,dens,1, WMass					

		

/com,------------------------------------------------------------------------------------

/com, key points
/com,*******

k,1,0,1226.875,0 
k,2,30.021,1226.875,30.550
k,3,60.042,1226.875,61.100 
k,4,90.064,1226.875,91.651   
k,5,105.154,1226.875,102.817 
k,6,139.302,1226.875,120.593 
k,7,181.878,1226.875,142.757 
k,8,224.435,1226.875,164.922 
k,9,243.383,1226.875,174.774 
k,10,262.311,1226.875,184.628
k,11,292.798,1226.875,191.342

k,12,334.171,1226.875,189.421 
!k,120,334.171,1226.875,189.421

k,13,375.543,1226.875,187.500

k,14,405.511,1226.875,186.110  
k,140,431.483,1226.875,184.904
 
k,15,501.172,1226.875,181.669 
k,16,570.860,1226.875,178.433  
k,17,579.777,1226.875,178.683 
k,18,615.118,1226.875,182.316 
k,20,633.028,1226.875,184.156

k,21,678.227,1226.875,188.802 
k,22,723.426,1226.875,193.448 
k,23,768.625,1226.875,198.095 
k,24,809.602,1226.875,187.256 
k,25,814.057,1226.875,184.079 
k,26,852.626,1226.875,156.568 
k,27,891.195,1226.875,129.058 
k,28,929.764,1226.875,101.547 

k,29,968.332,1226.875,74.036  
k,290,978.101,1226.875,67.067  

k,31,1012.600,1226.875,42.430
!k,310,1012.600,1226.875,42.430

k,32,1047.098,1226.875,17.793 
k,34,1061.752,1244.875,7.340   
k,35,1061.752,1272.375,7.340   
k,36,1072.214,1290.375,-7.307 
k,37,1081.623,1290.375,-20.48 
k,38,1108.85,1290.375,-58.399 
k,39,1136.077,1290.375,-96.317
k,40,1163.304,1290.375,-134.236 
k,41,1190.531,1290.375,-172.154

/com,
/com, Elastic Support 
/com,***********************

k,410,1190.531,1290.375,-172.154
   
k,43,1197.006,1290.375,-182.019   
k,44,1207.729,1290.375,-209.536   
k,45,1211.63,1290.375,-241.111
k,46,1215.531,1290.375,-272.687   
k,47,1219.432,1290.375,-304.262   
k,48,1223.333,1290.375,-335.873   
k,49,1227.234,1290.375,-367.413   

k,51,1232.114,1290.375,-407.115   
k,52,1233.704,1295.647,-419.787   
k,53,1234.945,1305.772,-429.836
   
k,55,1254.329,1318.500,-439.952 
k,56,1279.579,1318.500,-436.387 
k,57,1304.829,1318.500,-432.823 
k,58,1330.078,1318.500,-429.258 
k,59,1355.328,1318.500,-425.693 

!k,61,431.943,1226.875,194.899 
!k,62,616.14,1226.875,172.368  
!k,63,974.139,1226.875,82.176  
k,65,1227.234,1300.375,-367.413   
!k,66,1255.726,1318.500,-449.852

!k,67,1182.401,1290.375,-177.966   
k,68,105.154,1236.875,102.817 
k,69,224.435,1236.875,164.922 
k,70,405.511,1236.875,186.110  
k,71,633.028,1236.875,184.156 
k,72,814.057,1236.875,184.079 
k,73,978.101,1236.875,67.067  
k,74,1081.623,1300.375,-20.48 
k,75,1190.531,1300.375,-172.154   

!k,101,10,1226.875,0
!k,102,0,1236.875,0 
!k,103,0,1226.875,10

!k,591,1345.328,1318.5,-425.693 
!k,592,1355.328,1328.5,-425.693
!k,593,1355.328,1318.5,-415.693 

k,601,93.495,1226.9,94.879

k,701,90.618,1226.9,92.207
k,702,91.180,1226.9,92.757
k,703,91.748,1226.9,93.298
k,704,92.324,1226.9,93.833
k,705,92.906,1226.9,94.360

k,711,94.091,1226.9,95.390
k,712,94.694,1226.9,95.894
k,713,95.304,1226.9,96.390
k,714,95.919,1226.9,96.877
k,715,96.542,1226.9,97.357

k,602,97.170,1226.9,97.828

k,721,97.804,1226.9,98.291
k,722,98.443,1226.9,98.746
k,723,99.089,1226.9,99.192
k,724,99.740,1226.9,99.630
k,725,100.40,1226.9,100.06

k,731,101.73,1226.9,100.89
k,732,102.40,1226.9,101.29
k,733,103.08,1226.9,101.69
k,734,103.77,1226.9,102.07
k,735,104.46,1226.9,102.45

k,603,101.06,1226.9,100.48

k,741,267.08,1226.9,186.85
k,742,272.03,1226.9,188.65
k,743,277.11,1226.9,190.00
k,744,282.29,1226.9,190.91
k,745,287.54,1226.9,191.36

k,751,572.35,1226.9,178.38
k,752,573.84,1226.9,178.37
k,753,575.32,1226.9,178.39
k,754,576.81,1226.9,178.45
k,755,578.29,1226.9,178.55

k,761,775.83,1226.9,198.40
k,762,783.03,1226.9,197.84
k,763,790.10,1226.9,196.41
k,764,796.95,1226.9,194.15
k,765,803.49,1226.9,191.09

k,771,1199.7,1290.4,-186.18
k,772,1202.1,1290.4,-190.55
k,773,1204.1,1290.4,-195.11
k,774,1205.7,1290.4,-199.81
k,775,1206.9,1290.4,-204.63

k,781,1050.9,1227.5,15.088
k,782,1054.4,1229.3,12.567
k,783,1057.5,1232.1,10.402
k,784,1059.8,1235.9,8.7405
k,785,1061.3,1240.2,7.6962

k,791,1062.1,1277.0,6.8409
k,792,1063.2,1281.4,5.3777
k,793,1064.8,1285.1,3.0499
k,794,1067.0,1288.0,0.16395E-01
k,795,1069.5,1289.8,-3.5163

k,801,1232.4,1290.5,-409.45
k,802,1232.7,1291.0,-411.75
k,803,1233.0,1291.7,-413.97
k,804,1233.2,1292.8,-416.08
k,805,1233.5,1294.1,-418.02

k,811,1236.0,1309.0,-432.99
k,812,1238.2,1312.1,-435.76
k,813,1241.3,1314.7,-437.97
k,814,1245.3,1316.7,-439.46
k,815,1249.7,1318.0,-440.13

/com,------------------------------------------------------------------------------------

/com,
/com, Straight Pipe (Tangent) Elements
/com,**********************************

mat,1						! Material ID 1
type,1						! Element Type 1
secnum,1					! Section 1

l,1,2                       ! line number 1 
l,2,3   
l,3,4   

l,5,6
l,6,7   
l,7,8  
 
l,8,9   
l,9,10
  
l,11,12 
l,12,13 
l,13,14 
l,14,140
l,140,15 
l,15,16 

l,17,18
 
l,18,20  
l,20,21 
l,21,22 
l,22,23

l,24,25 
l,25,26 
l,26,27 
l,27,28 
l,28,29 
l,29,290
l,290,31 
l,31,32
 
l,34,35
 
l,36,37 
l,37,38 
l,38,39 
l,39,40 
l,40,41 
l,41,43
 
l,44,45 
l,45,46 
l,46,47 
l,47,48
l,48,49 
l,49,51

l,52,53

l,55,56
l,56,57
l,57,58
l,58,59             ! line number 45

lesize,all,,,1
lmesh,all
allsel,all

/com,
/com, Pipe Bend Elements
/com,********************

type,2
secnum,1
mat,1

larc,4,702,701, radcur1         ! line number 46
larc,702,704,703, radcur1
larc,704,601,705, radcur1
larc,601,712,711, radcur1
larc,712,714,713, radcur1
larc,714,602,715, radcur1
larc,602,722,721, radcur1
larc,722,724,723, radcur1
larc,724,603,725, radcur1
larc,603,732,731, radcur1
larc,732,734,733, radcur1
larc,734,5,735, radcur1
larc,10,742,741, radcur1
larc,742,744,743, radcur1
larc,744,11,745, radcur1
larc,16,752,751, radcur1
larc,752,754,753, radcur1
larc,754,17,755, radcur1
larc,23,762,761, radcur1
larc,762,764,763, radcur1
larc,764,24,765, radcur1
larc,43,772,771, radcur1
larc,772,774,773, radcur1
larc,774,44,775, radcur1            ! line number 69

lsel,s,line,,46,69,1
lesize,all,,,1
lmesh,all
allsel,all


type,3
secnum,1
mat,1

larc,32,782,781, radcur2            ! line number 70
larc,782,784,783, radcur2
larc,784,34,785, radcur2
larc,35,792,791, radcur2
larc,792,794,793, radcur2
larc,794,36,795, radcur2
larc,51,802,801, radcur2
larc,802,804,803, radcur2
larc,804,52,805, radcur2
larc,53,812,811, radcur2
larc,812,814,813, radcur2
larc,814,55,815, radcur2            ! line number 81


lsel,s,line,,70,81,1
lesize,all,,,1
lmesh,all
allsel,all
/com,

/com, Spring Elements
/com,*****************

type,4
real,4
l,49,65                         ! line number 82

lesize,82,,,1
lmesh,all
allsel,all

type,6
real,6
l,5,68                          ! line number 83
l,8,69
l,14,70
l,20,71
l,25,72
l,290,73
l,37,74
l,410,75                        ! line number 90

lsel,s,line,,83,90,1
lesize,all,,,1
lmesh,all
allsel,all


n1 = node(1254.329,1318.500,-439.952) 

n,66000,1255.726,1318.500,-449.852

n3 = node(1279.579,1318.500,-436.387) 
n4 = node(1190.531,1290.375,-172.154)

n, 67000,1182.401,1290.375,-177.966

n6 = node(1163.304,1290.375,-134.236)
n7 = node(375.543,1226.875,187.500)

n, 120000,334.171,1226.875,189.421

n9 = node(1355.328,1318.500,-425.693)
n10 = node(431.483,1226.875,184.904)

n,61000,431.943,1226.875,194.899

n12 = node(501.172,1226.875,181.669)
n13 = node(615.118,1226.875,182.316)

n,62000,616.14,1226.875,172.368

n15 = node(579.777,1226.875,178.683)
n16 = node(968.332,1226.875,74.036)

n,63000,974.139,1226.875,82.176

n18 = node(929.764,1226.875,101.547)
n19 = node(1047.098,1226.875,17.793)

n,310000,1012.600,1226.875,42.430

n21 = node(1061.752,1244.875,7.340)
n22 = node(1190.531,1290.375,-172.154)

!n,410000,1190.531,1290.375,-172.154

n,61000,431.943,1226.875,194.899
 
n,62000,616.14,1226.875,172.368

n,63000,974.139,1226.875,82.176

n27 = node(1227.234,1300.375,-367.413)   

n,66000,1255.726,1318.500,-449.852

n,67000,1182.401,1290.375,-177.966
     

n30 = node(105.154,1236.875,102.817) 
n31 = node(224.435,1236.875,164.922) 
n32 = node(405.511,1236.875,186.110)  
n33 = node(633.028,1236.875,184.156) 
n34 = node(814.057,1236.875,184.079) 
n35 = node(978.101,1236.875,67.067)  
n36 = node(1081.623,1300.375,-20.48) 
n37 = node(1190.531,1300.375,-172.154)
   
n,101000,10,1226.875,0

n,102000,0,1236.875,0 

n,103000,0,1226.875,10   

n,120000,334.171,1226.875,189.421

n,310000,1012.600,1226.875,42.430

n,591000,1345.328,1318.5,-425.693
 
n,592000,1355.328,1328.5,-425.693

n,593000,1355.328,1318.5,-415.693

n46 = node(0,1226.875,0) 
 

type,5
real,5
n1 = 92
n2 = 66000
n3 = 93

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
e,n1,n2

n1 = 72
n2 = 67000
n3 = 70

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 
csys,0
e,n1,n2

/com,
/com, Snubber Elements
/com,******************

n1 = 22
n2 = 120000
n3 = 99

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 
csys,0
e,n1,n2

n1 = 26
n2 = 61000
n3 = 28

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 
csys,0
e,n1,n2

n1 = 33
n2 = 62000
n3 = 32

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 
csys,0
e,n1,n2

n1 = 52
n2 = 63000
n3 = 50

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 
csys,0
e,n1,n2

n1 = 58
n2 = 310000
n3 = 60
!n3 = 59

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 
csys,0
e,n1,n2

/com,
/com, 3D Support at both ends
/com,*************************

type,5
real,5
e,1,101000
e,99,591000

type,6
real,6
e,1,102000
e,99,592000

type,7
real,7
e,1,103000
e,99,593000

type,8
real,8
e,1,101000
e,99,591000

type,9
real,9
e,1,102000
e,99,592000

type,10
real,10
e,1,103000
e,99,593000

/com,------------------------------------------------------------------------------------

/com,
/com, Model Rigid Region
/com,********************

cerig,72,172,uy

/com,------------------------------------------------------------------------------------

/com,
/com, Convert some PIPE289 into ELBOW290 using ELBOW command
/com, 

elbow,on,,,sect
allsel,all


/com,
/com, Constraints
/com,*************

nsel,s,node,,61000,63000
nsel,a,node,,164
nsel,a,node,,66000,67000,1
nsel,a,node,,165,173,1
nsel,a,node,,101000,103000
nsel,a,node,,120000
nsel,a,node,,310000
nsel,a,node,,591000,593000
d,all,all,0

nsel,all

/com,------------------------------------------------------------------------------------

/com,
/com, Loads
/com,*******

/com, Temperature Input
/com,*******************

bf,all,temp,Temperature

esel,r,type,,1
esel,a,type,,2
esel,a,type,,3

/com, Pressure Input
/com,****************

sfe,all,1,pres,,Pressue,,,

allsel,all,all
finish

/com,------------------------------------------------------------------------------------

/com,
/com,=============
/com,	Modal Solve
/com,=============
/com,

/solution
antype,modal							! Perform modal solve
modopt,lanb,maxm
lump,on								! Lumped mass formulation
mxpand,maxm,,,yes							! Expand the modes with stress calculation
solve
save,
finish


/post1
/out,
/com, ********************************************
/com,  Frequencies obtained from modal solution
/com, *********************************************
set,list
finish

/com,
/com,================
/com,	Spectrum Solve
/com,================
/com,
/out,scratch

/solution
sfedele,all,1,pres,,,
bfdele,all,temp,,,

antype,spectrum					
spopt,mprs,maxm					! Multi-point response spectrum analysis

srss,,,						
gval = 386.4

/com,------------------------------------------------------------------------------------
/com,

/com,
/com, Spectrum 1 (Group 1 - X - Direction Excitation)
/com,*************************************************

id = 1
spunit, id,accg, gval

spfrq,id, 1.0,1.0428,1.1025,1.1905,1.2270,1.2739,1.2937
spval,id,,  1.162,1.282,1.399,1.549,1.606,1.676,1.704

spfrq,id, 1.3423,1.3889,1.4104,1.4347,1.5552,1.6949,1.7825
spval,id,,  1.774,1.839,1.869,1.904,2.084,2.246,2.304

spfrq,id, 1.9305,2.0747,2.2779,2.4752,2.6042,2.6596,2.9499
spval,id,,  2.382,2.479,2.592,2.644,2.64,2.639,2.782

spfrq,id, 3.2362,3.3898,3.4965,3.5714,3.6101,3.663,3.7313
spval,id,,  2.951,3.066,3.215,3.384,3.532,3.766,4.189

spfrq,id, 3.8168,3.8911,3.9216,4.218,4.6948,4.7847,5.0505
spval,id,, 4.793,5.189,5.224,5.232,5.227,5.152,3.002

spfrq,id, 5.0761,5.3476,5.7471,5.9524,5.9880,6.6225,7.4627
spval,id,,  2.923,2.90,2.873,2.849,2.844,2.761,2.667

spfrq,id, 7.8125,7.874,7.9365,8.3333,8.9286,9.5238,9.6154
spval,id,,  2.635,2.685,2.755,2.087,2.797,2.744,2.674

spfrq,id, 9.7087,10.4167,10.8696,11.6279,11.7647,12.1951,12.50
spval,id,,  2.627,2.781,2.931,3.077,3.112,3.134,3.134

spfrq,id, 12.8205,13.1579,13.3333,13.4953,13.5135,13.8889,14.2857
spval,id,, 3.116,2.975,2.687,2.56,2.399,2.064,1.855

spfrq,id, 15.3846,15.625,17.8571,18.8679,22.7273,23.8095,24.3902
spval,id,,  1.524,1.512,1.472,1.335,1.09,1.073,1.07

spfrq,id, 25.641,26.3158,27.027,27.7778,28.5714,40.0,76.9231
spval,id,,  1.049,1.004,0.9823,0.9669,0.956,0.893,0.83

spfrq,id, 1000.0
spval,id,, 0.771

/com,
/com, Spectrum 2 (Group 1 - Y - Direction Excitation)
/com,*************************************************

id = 2
spunit, id,accg, gval

spfrq,id, 0.5,1.4993,1.6207,1.9011,2.0704,3.8023,4.2553
spval,id,,  0.35,1.45,1.8,1.8,2.61,2.61,2.78

spfrq,id, 5.1813,5.4054,7.8125,8.1301,9.901,11.5207,14.1044
spval,id,,  2.78,2.58,2.58,3.25,3.25,3.62,3.62

spfrq,id, 14.4928,17.6991,23.9981,59.988
spval,id,,  3.05,3.05,1.20,0.75

/com,
/com, Spectrum 3 (Group 1 - Z - Direction Excitation)
/com,*************************************************

id = 3
spunit, id,accg, gval

spfrq,id, 1.0,1.0428,1.1025,1.1905,1.2270,1.2739,1.2937
spval,id,,  1.162,1.282,1.399,1.549,1.606,1.676,1.704

spfrq,id, 1.3423,1.3889,1.4104,1.4347,1.5552,1.6949,1.7825
spval,id,,  1.774,1.839,1.869,1.904,2.084,2.246,2.304

spfrq,id, 1.9305,2.0747,2.2779,2.4752,2.6042,2.6596,2.9499
spval,id,,  2.382,2.479,2.592,2.644,2.64,2.639,2.782

spfrq,id, 3.2362,3.3898,3.4965,3.5714,3.6101,3.663,3.7313
spval,id,,  2.951,3.066,3.215,3.384,3.532,3.766,4.189

spfrq,id, 3.8168,3.8911,3.9216,4.218,4.6948,4.7847,5.0505
spval,id,, 4.793,5.189,5.224,5.232,5.227,5.152,3.002

spfrq,id, 5.0761,5.3476,5.7471,5.9524,5.9880,6.6225,7.4627
spval,id,,  2.923,2.90,2.873,2.849,2.844,2.761,2.667

spfrq,id, 7.8125,7.874,7.9365,8.3333,8.9286,9.5238,9.6154
spval,id,,  2.635,2.685,2.755,2.087,2.797,2.744,2.674

spfrq,id, 9.7087,10.4167,10.8696,11.6279,11.7647,12.1951,12.50
spval,id,,  2.627,2.781,2.931,3.077,3.112,3.134,3.134

spfrq,id, 12.8205,13.1579,13.3333,13.4953,13.5135,13.8889,14.2857
spval,id,, 3.116,2.975,2.687,2.56,2.399,2.064,1.855

spfrq,id, 15.3846,15.625,17.8571,18.8679,22.7273,23.8095,24.3902
spval,id,,  1.524,1.512,1.472,1.335,1.09,1.073,1.07

spfrq,id, 25.641,26.3158,27.027,27.7778,28.5714,40.0,76.9231
spval,id,,  1.049,1.004,0.9823,0.9669,0.956,0.893,0.83

spfrq,id, 1000.0
spval,id,, 0.771

/com,------------------------------------------------------------------------------------
/com,

/com,
/com, Spectrum 4 (Group 2 - X - Direction Excitation)
/com,*************************************************

id = 4
spunit, id,accg, gval

spfrq,id, 1.0,1.0428,1.1025,1.1905,1.2270,1.2739,1.2937
spval,id,,  1.162,1.282,1.399,1.549,1.606,1.676,1.704

spfrq,id, 1.3423,1.3889,1.4104,1.4347,1.5552,1.6949,1.7825
spval,id,,  1.774,1.839,1.869,1.904,2.084,2.246,2.304

spfrq,id, 1.9305,2.0747,2.2779,2.4752,2.6042,2.6596,2.9499
spval,id,,  2.382,2.479,2.592,2.644,2.64,2.639,2.782

spfrq,id, 3.2362,3.3898,3.4965,3.5714,3.6101,3.663,3.7313
spval,id,,  2.951,3.066,3.215,3.384,3.532,3.766,4.189

spfrq,id, 3.8168,3.8911,3.9216,4.218,4.6948,4.7847,5.0505
spval,id,, 4.793,5.189,5.224,5.232,5.227,5.152,3.002

spfrq,id, 5.0761,5.3476,5.7471,5.9524,5.9880,6.6225,7.4627
spval,id,,  2.923,2.90,2.873,2.849,2.844,2.761,2.667

spfrq,id, 7.8125,7.874,7.9365,8.3333,8.9286,9.5238,9.6154
spval,id,,  2.635,2.685,2.755,2.087,2.797,2.744,2.674

spfrq,id, 9.7087,10.4167,10.8696,11.6279,11.7647,12.1951,12.50
spval,id,,  2.627,2.781,2.931,3.077,3.112,3.134,3.134

spfrq,id, 12.8205,13.1579,13.3333,13.4953,13.5135,13.8889,14.2857
spval,id,, 3.116,2.975,2.687,2.56,2.399,2.064,1.855

spfrq,id, 15.3846,15.625,17.8571,18.8679,22.7273,23.8095,24.3902
spval,id,,  1.524,1.512,1.472,1.335,1.09,1.073,1.07

spfrq,id, 25.641,26.3158,27.027,27.7778,28.5714,40.0,76.9231
spval,id,,  1.049,1.004,0.9823,0.9669,0.956,0.893,0.83

spfrq,id, 1000.0
spval,id,, 0.771

/com,
/com, Spectrum 5 (Group 2 - Y - Direction Excitation)
/com,*************************************************

id = 5
spunit, id,accg, gval

spfrq,id, 0.50,1.4556,1.9011,2.0704,2.907,4.065,4.9603
spval,id,,  1.45,1.8,1.8,2.68,2.68,3.17,3.17

spfrq,id, 5.0,7.5988,8.5034,10.9051,11.5207,14.0845,16.0
spval,id,,  3.03,3.03,4.82,4.82,5.95,5.95,4.49

spfrq,id, 19.1205,21.0084,50.0
spval,id,,  4.49,1.85,1.05

/com,
/com, Spectrum 6 (Group 2 - Z - Direction Excitation)
/com,*************************************************

id = 6
spunit, id,accg, gval

spfrq,id, 1.0,1.0428,1.1025,1.1905,1.2270,1.2739,1.2937
spval,id,,  1.162,1.282,1.399,1.549,1.606,1.676,1.704

spfrq,id, 1.3423,1.3889,1.4104,1.4347,1.5552,1.6949,1.7825
spval,id,,  1.774,1.839,1.869,1.904,2.084,2.246,2.304

spfrq,id, 1.9305,2.0747,2.2779,2.4752,2.6042,2.6596,2.9499
spval,id,,  2.382,2.479,2.592,2.644,2.64,2.639,2.782

spfrq,id, 3.2362,3.3898,3.4965,3.5714,3.6101,3.663,3.7313
spval,id,,  2.951,3.066,3.215,3.384,3.532,3.766,4.189

spfrq,id, 3.8168,3.8911,3.9216,4.218,4.6948,4.7847,5.0505
spval,id,, 4.793,5.189,5.224,5.232,5.227,5.152,3.002

spfrq,id, 5.0761,5.3476,5.7471,5.9524,5.9880,6.6225,7.4627
spval,id,,  2.923,2.90,2.873,2.849,2.844,2.761,2.667

spfrq,id, 7.8125,7.874,7.9365,8.3333,8.9286,9.5238,9.6154
spval,id,,  2.635,2.685,2.755,2.087,2.797,2.744,2.674

spfrq,id, 9.7087,10.4167,10.8696,11.6279,11.7647,12.1951,12.50
spval,id,,  2.627,2.781,2.931,3.077,3.112,3.134,3.134

spfrq,id, 12.8205,13.1579,13.3333,13.4953,13.5135,13.8889,14.2857
spval,id,, 3.116,2.975,2.687,2.56,2.399,2.064,1.855

spfrq,id, 15.3846,15.625,17.8571,18.8679,22.7273,23.8095,24.3902
spval,id,,  1.524,1.512,1.472,1.335,1.09,1.073,1.07

spfrq,id, 25.641,26.3158,27.027,27.7778,28.5714,40.0,76.9231
spval,id,,  1.049,1.004,0.9823,0.9669,0.956,0.893,0.83

spfrq,id, 1000.0
spval,id,, 0.771

/com,------------------------------------------------------------------------------------
/com,

/com,
/com, Spectrum 7 (Group 3 - X - Direction Excitation)
/com,*************************************************

id = 7
spunit, id,accg, gval

spfrq,id, 1.0,1.0428,1.1025,1.1905,1.2270,1.2739,1.2937
spval,id,,  1.162,1.282,1.399,1.549,1.606,1.676,1.704

spfrq,id, 1.3423,1.3889,1.4104,1.4347,1.5552,1.6949,1.7825
spval,id,,  1.774,1.839,1.869,1.904,2.084,2.246,2.304

spfrq,id, 1.9305,2.0747,2.2779,2.4752,2.6042,2.6596,2.9499
spval,id,,  2.382,2.479,2.592,2.644,2.64,2.639,2.782

spfrq,id, 3.2362,3.3898,3.4965,3.5714,3.6101,3.663,3.7313
spval,id,,  2.951,3.066,3.215,3.384,3.532,3.766,4.189

spfrq,id, 3.8168,3.8911,3.9216,4.218,4.6948,4.7847,5.0505
spval,id,, 4.793,5.189,5.224,5.232,5.227,5.152,3.002

spfrq,id, 5.0761,5.3476,5.7471,5.9524,5.9880,6.6225,7.4627
spval,id,,  2.923,2.90,2.873,2.849,2.844,2.761,2.667

spfrq,id, 7.8125,7.874,7.9365,8.3333,8.9286,9.5238,9.6154
spval,id,,  2.635,2.685,2.755,2.087,2.797,2.744,2.674

spfrq,id, 9.7087,10.4167,10.8696,11.6279,11.7647,12.1951,12.50
spval,id,,  2.627,2.781,2.931,3.077,3.112,3.134,3.134

spfrq,id, 12.8205,13.1579,13.3333,13.4953,13.5135,13.8889,14.2857
spval,id,, 3.116,2.975,2.687,2.56,2.399,2.064,1.855

spfrq,id, 15.3846,15.625,17.8571,18.8679,22.7273,23.8095,24.3902
spval,id,,  1.524,1.512,1.472,1.335,1.09,1.073,1.07

spfrq,id, 25.641,26.3158,27.027,27.7778,28.5714,40.0,76.9231
spval,id,,  1.049,1.004,0.9823,0.9669,0.956,0.893,0.83

spfrq,id, 1000.0
spval,id,, 0.771

/com,
/com, Spectrum 8 (Group 3 - Y - Direction Excitation)
/com,*************************************************

id = 8
spunit, id,accg, gval

spfrq,id, 0.5,1.9011,2.0704,3.003,4.0486,4.9505,7.1942
spval,id,,  0.4,1.88,2.72,2.72,3.42,3.42,3.82

spfrq,id, 8.1301,9.5238,10.352,12.6422,13.5135,15.4083,15.7978
spval,id,, 7.21,7.21,8.18,8.18,6.39,6.39,5.92    

spfrq,id, 17.6991,21.0084,50.0
spval,id,,  5.92,2.25,1.55

/com,
/com, Spectrum 9 (Group 3 - Z - Direction Excitation)
/com,*************************************************

id = 9
spunit, id,accg, gval

spfrq,id, 1.0,1.0428,1.1025,1.1905,1.2270,1.2739,1.2937
spval,id,,  1.162,1.282,1.399,1.549,1.606,1.676,1.704

spfrq,id, 1.3423,1.3889,1.4104,1.4347,1.5552,1.6949,1.7825
spval,id,,  1.774,1.839,1.869,1.904,2.084,2.246,2.304

spfrq,id, 1.9305,2.0747,2.2779,2.4752,2.6042,2.6596,2.9499
spval,id,,  2.382,2.479,2.592,2.644,2.64,2.639,2.782

spfrq,id, 3.2362,3.3898,3.4965,3.5714,3.6101,3.663,3.7313
spval,id,,  2.951,3.066,3.215,3.384,3.532,3.766,4.189

spfrq,id, 3.8168,3.8911,3.9216,4.218,4.6948,4.7847,5.0505
spval,id,, 4.793,5.189,5.224,5.232,5.227,5.152,3.002

spfrq,id, 5.0761,5.3476,5.7471,5.9524,5.9880,6.6225,7.4627
spval,id,,  2.923,2.90,2.873,2.849,2.844,2.761,2.667

spfrq,id, 7.8125,7.874,7.9365,8.3333,8.9286,9.5238,9.6154
spval,id,,  2.635,2.685,2.755,2.087,2.797,2.744,2.674

spfrq,id, 9.7087,10.4167,10.8696,11.6279,11.7647,12.1951,12.50
spval,id,,  2.627,2.781,2.931,3.077,3.112,3.134,3.134

spfrq,id, 12.8205,13.1579,13.3333,13.4953,13.5135,13.8889,14.2857
spval,id,, 3.116,2.975,2.687,2.56,2.399,2.064,1.855

spfrq,id, 15.3846,15.625,17.8571,18.8679,22.7273,23.8095,24.3902
spval,id,,  1.524,1.512,1.472,1.335,1.09,1.073,1.07

spfrq,id, 25.641,26.3158,27.027,27.7778,28.5714,40.0,76.9231
spval,id,,  1.049,1.004,0.9823,0.9669,0.956,0.893,0.83

spfrq,id, 1000.0
spval,id,, 0.771

/com,------------------------------------------------------------------------------------
/com,

/com,
/com, Spectrum 10 (Group 4 - X - Direction Excitation)
/com,**************************************************

id = 10
spunit, id,accg, gval

spfrq,id, 1.0,1.0428,1.1025,1.1905,1.2270,1.2739,1.2937
spval,id,,  1.162,1.282,1.399,1.549,1.606,1.676,1.704

spfrq,id, 1.3423,1.3889,1.4104,1.4347,1.5552,1.6949,1.7825
spval,id,,  1.774,1.839,1.869,1.904,2.084,2.246,2.304

spfrq,id, 1.9305,2.0747,2.2779,2.4752,2.6042,2.6596,2.9499
spval,id,,  2.382,2.479,2.592,2.644,2.64,2.639,2.782

spfrq,id, 3.2362,3.3898,3.4965,3.5714,3.6101,3.663,3.7313
spval,id,,  2.951,3.066,3.215,3.384,3.532,3.766,4.189

spfrq,id, 3.8168,3.8911,3.9216,4.218,4.6948,4.7847,5.0505
spval,id,, 4.793,5.189,5.224,5.232,5.227,5.152,3.002

spfrq,id, 5.0761,5.3476,5.7471,5.9524,5.9880,6.6225,7.4627
spval,id,,  2.923,2.90,2.873,2.849,2.844,2.761,2.667

spfrq,id, 7.8125,7.874,7.9365,8.3333,8.9286,9.5238,9.6154
spval,id,,  2.635,2.685,2.755,2.087,2.797,2.744,2.674

spfrq,id, 9.7087,10.4167,10.8696,11.6279,11.7647,12.1951,12.50
spval,id,,  2.627,2.781,2.931,3.077,3.112,3.134,3.134

spfrq,id, 12.8205,13.1579,13.3333,13.4953,13.5135,13.8889,14.2857
spval,id,, 3.116,2.975,2.687,2.56,2.399,2.064,1.855

spfrq,id, 15.3846,15.625,17.8571,18.8679,22.7273,23.8095,24.3902
spval,id,,  1.524,1.512,1.472,1.335,1.09,1.073,1.07

spfrq,id, 25.641,26.3158,27.027,27.7778,28.5714,40.0,76.9231
spval,id,,  1.049,1.004,0.9823,0.9669,0.956,0.893,0.83

spfrq,id, 1000.0
spval,id,, 0.771

/com,
/com, Spectrum 11 (Group 4 - Y - Direction Excitation)
/com,**************************************************

id = 11
spunit, id,accg, gval

spfrq,id, 0.5,1.3004,2.0,3.003,4.0486,4.9505,6.993
spval,id,,  0.5,1.4,2.75,2.75,3.5,3.5,4.5 

spfrq,id, 8.1301,10.4167,11.4943,14.845,17.5439,22.2222,50.0
spval,id,,  12.1,12.1,10.7,10.7,2.7,1.8,1.5 

/com,
/com, Spectrum 12 (Group 4 - Z - Direction Excitation)
/com,**************************************************

id = 12
spunit, id,accg, gval

spfrq,id, 1.0,1.0428,1.1025,1.1905,1.2270,1.2739,1.2937
spval,id,,  1.162,1.282,1.399,1.549,1.606,1.676,1.704

spfrq,id, 1.3423,1.3889,1.4104,1.4347,1.5552,1.6949,1.7825
spval,id,,  1.774,1.839,1.869,1.904,2.084,2.246,2.304

spfrq,id, 1.9305,2.0747,2.2779,2.4752,2.6042,2.6596,2.9499
spval,id,,  2.382,2.479,2.592,2.644,2.64,2.639,2.782

spfrq,id, 3.2362,3.3898,3.4965,3.5714,3.6101,3.663,3.7313
spval,id,,  2.951,3.066,3.215,3.384,3.532,3.766,4.189

spfrq,id, 3.8168,3.8911,3.9216,4.218,4.6948,4.7847,5.0505
spval,id,, 4.793,5.189,5.224,5.232,5.227,5.152,3.002

spfrq,id, 5.0761,5.3476,5.7471,5.9524,5.9880,6.6225,7.4627
spval,id,,  2.923,2.90,2.873,2.849,2.844,2.761,2.667

spfrq,id, 7.8125,7.874,7.9365,8.3333,8.9286,9.5238,9.6154
spval,id,,  2.635,2.685,2.755,2.087,2.797,2.744,2.674

spfrq,id, 9.7087,10.4167,10.8696,11.6279,11.7647,12.1951,12.50
spval,id,,  2.627,2.781,2.931,3.077,3.112,3.134,3.134

spfrq,id, 12.8205,13.1579,13.3333,13.4953,13.5135,13.8889,14.2857
spval,id,, 3.116,2.975,2.687,2.56,2.399,2.064,1.855

spfrq,id, 15.3846,15.625,17.8571,18.8679,22.7273,23.8095,24.3902
spval,id,,  1.524,1.512,1.472,1.335,1.09,1.073,1.07

spfrq,id, 25.641,26.3158,27.027,27.7778,28.5714,40.0,76.9231
spval,id,,  1.049,1.004,0.9823,0.9669,0.956,0.893,0.83

spfrq,id, 1000.0
spval,id,, 0.771

/com,------------------------------------------------------------------------------------
/com,
/com,
/com, Nodal Components for Excitation Points (non rotated)
/com,******************************************************

allsel,all,all

nsel,s,node,,165
nsel,a,node,,61000
nsel,a,node,,168
nsel,a,node,,63000
nsel,a,node,,170
nsel,a,node,,310000
nsel,a,node,,591000,593000
cm,group_1,node

allsel,all,all

nsel,s,node,,167
nsel,a,node,,120000
nsel,a,node,,62000
nsel,a,node,,67000
nsel,a,node,,173	
nsel,a,node,,164   
cm,group_2,node

allsel,all,all
nsel,s,node,,166
nsel,a,node,,169  
nsel,a,node,,171
cm,group_3,node

allsel,all,all
nsel,s,node,,101000,103000
nsel,a,node,,66000
cm,group_4,node

allsel,all,all
eplo
 
/com,
/com, *****************************************************
! -- Support Group 1 - spectrum 1 (Along X - Direction)

sed,1,,,group_1
pfact,1
sed,0,,,group_1

! -- Support Group 1 - spectrum 2 (Along Y - Direction)

sed,,1,,group_1
pfact,2
sed,,0,,group_1

! -- Support Group 1 - spectrum 3 (Along Z - Direction)

sed,,,1,group_1
pfact,3
sed,,,0,group_1

! -- Support Group 2 - spectrum 4 (Along X - Direction)

sed,1,,,group_2
pfact,4
sed,0,,,group_2

! -- Support Group 2 - spectrum 5 (Along Y - Direction)

sed,,1,,group_2
pfact,5
sed,,0,,group_2

! -- Support Group 2 - spectrum 6 (Along Z - Direction)

sed,,,1,group_2
pfact,6
sed,,,0,group_2


! -- Support Group 3 - spectrum 8 (Along Y - Direction)

sed,,1,,group_3
pfact,8
sed,,0,,group_3


! -- Support Group 4 - spectrum 10 (Along X - Direction)

sed,1,,,group_4
pfact,10
sed,0,,,group_4

! -- Support Group 4 - spectrum 11 (Along Y - Direction)

sed,,1,,group_4
pfact,11
sed,,0,,group_4

! -- Support Group 4 - spectrum 12 (Along Z - Direction)

sed,,,1,group_4
pfact,12
sed,,,0,group_4


srss,0,,YES			! activate Absolute Sum for MPRS

solve

fini
/com, *****************************************************
/post1
/input,,mcom

/out,

/com,
/com,===============================================================
/com, Maximum nodal displacements and rotations from spectrum solution
/com,===============================================================
/com,

/out,scratch

*GET,AdisX,NODE,63,U,X
*GET,AdisY,NODE,87,U,Y
*GET,AdisZ,NODE,63,U,Z
*GET,ArotX,NODE,58,ROT,X
*GET,ArotY,NODE,52,ROT,Y
*GET,ArotZ,NODE,61,ROT,Z

/out,
*stat,AdisX
*stat,AdisY
*stat,AdisZ
*stat,ArotX
*stat,ArotY
*stat,ArotZ

/com,==============================================================
/com, Element Forces and Moments obtained from spectrum solution
/com,==============================================================

/out,scratch

/com,==========
/com,	 Node I
/com,==========

/com, Element #1 (Pipe289 element)
/com,************************


esel,s,elem,,1
etable,pxi_1,smisc,1
etable,vyi_1,smisc,6
etable,vzi_1,smisc,5
etable,txi_1,smisc,4
etable,myi_1,smisc,2
etable,mzi_1,smisc,3
esel,all

/out,
/com, ************************************************
/com,  Element forces and moments at element1, node i
/com, ************************************************

pretab,pxi_1,vyi_1,vzi_1,txi_1,myi_1,mzi_1


/out,scratch
/com,==========
/com,  Node J
/com,==========

/com, Element #1 (Pipe289 element)  
/com,************************

esel,s,elem,,1
etable,pxj_1,smisc,14
etable,vyj_1,smisc,19
etable,vzj_1,smisc,18
etable,txj_1,smisc,17
etable,myj_1,smisc,15
etable,mzj_1,smisc,16
esel,all

/out,
/com, **********************************************
/com,  Element forces and moments at element1, node j
/com, ************************************************

pretab,pxj_1,vyj_1,vzj_1,txj_1,myj_1,mzj_1



/out,scratch


/com,==========
/com,  Node I
/com,==========

/com, Element #50 (Elbow 290 element)
/com,*************************

esel,s,elem,,50
etable,pxi_50,smisc,1
etable,vyi_50,smisc,6
etable,vzi_50,smisc,5
etable,txi_50,smisc,4
etable,myi_50,smisc,2
etable,mzi_50,smisc,3
esel,all

/out,
/com, *************************************************
/com,  Element forces and moments at element 50, node i
/com, ************************************************

pretab,pxi_50,vyi_50,vzi_50,txi_50,myi_50,mzi_50


/out,scratch
/com,==========
/com,  Node J
/com,==========



/com, Element #50 (Elbow290 element)
/com,********************************
esel,s,elem,,50

etable,pxj_50,smisc,36
etable,vyj_50,smisc,41
etable,vzj_50,smisc,40
etable,txj_50,smisc,39
etable,myj_50,smisc,37
etable,mzj_50,smisc,38
esel,all

allsel,all

/out,
/com, **********************************************
/com,  Element forces and moments at element 50, node j
/com, ************************************************

pretab,pxj_50,vyj_50,vzj_50,txj_50,myj_50,mzj_50



/com,----------------------------------------------------------------------------

/com, ***************************************
/com,  Reaction forces from spectrum solution
/com, ****************************************

prrsol

finish