/COM,ANSYS MEDIA REL. 2024R2 (05/10/2024) REF. VERIF. MANUAL: REL. 2024R2
/verify,vm-nr1677-2-4a-a
/title,vm-nr1677-2-4a-a,NRC piping benchmarks problems,Volume II,Problem 4a
/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,
/com, Element used: Pipe16, Pipe18, Combin14 and Mass21
/com,
/com, Results comparsion:
/com, The following results are compared against NRC piping benchmark values
/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
YoungModulus = 0.283e+8 ! Young's Modulus
Nu = 0.3 ! Minor Poisson's Ratio
ShearModulus = YoungModulus/(2*(1+Nu)) ! Shear Modulus
K = 0.911e-5 ! Thermal Expansion
maxm = 50 ! No. of Modes Extracted
/com,
/com, Section Property
/com,******************
/com,
/com, Wall Thickness
/com,----------------
WTick_1 = 2.25
WTick_2 = 3.44
WTick_3 = 1.0
WTick_4 = 2.64
WTick_5 = 74.775
WTick_6 = 1.0
WTick_7 = 2.64
WTick_8 = 0.7180
WTick_9 = 1.62
WTick_10 = 46.035
WTick_11 = 0.7180
WTick_12 = 0.906
WTick_13 = 0.365
/com,
/com, Outer Diameter
/com,----------------
OD_1 = 32.25
OD_2 = 15.625
OD_3 = 10.75
OD_4 = 16.03
OD_5 = 160.3
OD_6 = 10.75
OD_7 = 16.03
OD_8 = 6.625
OD_9 = 9.87
OD_10 = 98.7
OD_11 = 6.625
OD_12 = 8.625
OD_13 = 10.75
/com,
/com, Bend Radius
/com,--------------
RADCUR_1 = 15.0
RADCUR_2 = 14.9
RADCUR_3 = 9.0
RADCUR_4 = 12.0
RADCUR_5 = 40.0
RADCUR_6 = 8.0
/com,-----------------------------------------------------------------------------------
/com,
/com, Element Types
/com,***************
et,1,pipe16
et,2,pipe16
et,3,pipe16
et,4,pipe16
et,5,pipe16
et,6,pipe16
et,7,pipe16
et,8,pipe16
et,9,pipe16
et,10,pipe16
et,11,pipe16
et,12,pipe16
et,13,pipe16
et,14,pipe18
et,15,pipe18
et,16,pipe18
et,17,pipe18
et,18,pipe18
et,19,pipe18
et,20,pipe18
et,21,pipe18
et,22,pipe18
et,23,mass21 ! 3D Mass without rotatary inertia
keyopt,23,3,2
et,24,mass21
keyopt,24,3,2
et,25,mass21
keyopt,25,3,2
et,26,mass21
keyopt,26,3,2
et,27,mass21
keyopt,27,3,2
et,28,mass21
keyopt,28,3,2
et,29,mass21
keyopt,29,3,2
et,30,mass21
keyopt,30,3,2
et,31,mass21
keyopt,31,3,2
et,32,mass21
keyopt,32,3,2
et,33,mass21
keyopt,33,3,2
et,34,mass21
keyopt,34,3,2
et,35,mass21
keyopt,35,3,2
et,36,mass21
keyopt,36,3,2
et,37,mass21
keyopt,37,3,2
et,38,mass21
keyopt,38,3,2
et,39,mass21
keyopt,39,3,2
et,40,mass21
keyopt,40,3,2
et,41,mass21
keyopt,41,3,2
et,42,mass21
keyopt,42,3,2
et,43,mass21
keyopt,43,3,2
et,44,mass21
keyopt,44,3,2
et,45,mass21
keyopt,45,3,2
et,46,mass21
keyopt,46,3,2
et,47,mass21
keyopt,47,3,2
et,48,combin14,,1 ! Spring Elements Types and Real constants
keyopt,48,2,1 ! X Degree of Freedom
et,49,combin14,,2
keyopt,49,2,2 ! Y Degree of Freedom
et,50,combin14,,3
keyopt,50,2,3 ! Z Degree of Freedom
et,51,combin14,,2
keyopt,51,2,2 ! Y Degree of Freedom
et,52,combin14,,2
keyopt,52,2,2 ! Y Degree of Freedom
et,53,combin14,,3
keyopt,53,2,3 ! Z Degree of Freedom
et,54,combin14,,1
keyopt,54,2,1 ! X Degree of Freedom
et,55,combin14,,2
keyopt,55,2,2 ! Y Degree of Freedom
et,56,combin14,,1
keyopt,56,2,1 ! X Degree of Freedom
et,57,combin14,,1
keyopt,57,2,1 ! X Degree of Freedom
et,58,combin14,,2
keyopt,58,2,2 ! Y Degree of Freedom
et,59,combin14,,1
keyopt,59,2,1 ! X Degree of Freedom
et,60,combin14,,2
keyopt,60,2,2 ! Y Degree of Freedom
et,61,combin14,,1
keyopt,61,2,1 ! X Degree of Freedom
et,62,combin14,,2
keyopt,62,2,2 ! Y Degree of Freedom
et,63,combin14,,2
keyopt,63,2,2 ! Y Degree of Freedom
et,64,combin14,,1
keyopt,64,2,1 ! X Degree of Freedom
et,65,combin14,,2
keyopt,65,2,2 ! Y Degree of Freedom
et,66,combin14,,3
keyopt,66,2,3 ! Z Degree of Freedom
et,67,combin14,,2
keyopt,67,2,2 ! Y Degree of Freedom
et,68,combin14,,1
keyopt,68,2,1 ! X Degree of Freedom
et,69,combin14,,1
keyopt,69,2,1 ! X Degree of Freedom
et,70,combin14,,2
keyopt,70,2,2 ! Y Degree of Freedom
et,71,combin14,,1
keyopt,71,2,1 ! X Degree of Freedom
et,72,combin14,,2
keyopt,72,2,2 ! Y Degree of Freedom
et,73,combin14,,1
keyopt,73,2,1 ! X Degree of Freedom
et,74,combin14,,1
keyopt,74,2,1 ! X Degree of Freedom
et,75,combin14,,3
keyopt,75,2,3 ! Z Degree of Freedom
et,76,combin14,,2
keyopt,76,2,2 ! Y Degree of Freedom
et,77,combin14,,2
keyopt,77,2,2 ! Y Degree of Freedom
et,78,combin14,,1
keyopt,78,2,1 ! X Degree of Freedom
et,79,combin14,,2
keyopt,79,2,2 ! Y Degree of Freedom
et,80,combin14,,2
keyopt,80,2,2 ! Y Degree of Freedom
et,81,combin14,,1
keyopt,81,2,1 ! X Degree of Freedom
et,82,combin14,,2
keyopt,82,2,2 ! Y Degree of Freedom
et,83,combin14,,3
keyopt,83,2,3 ! Z Degree of Freedom
/com,----------------------------------------------------------------------------------
/com
/com, Real Constants
/com,****************
r, 1,OD_1,WTick_1
r, 2,OD_2,WTick_2
r, 3,OD_3,WTick_3
r, 4,OD_4,WTick_4
r, 5, OD_5,WTick_5
r, 6, OD_6,WTick_6
r, 7, OD_7,WTick_7
r, 8, OD_8,WTick_8
r, 9, OD_9,WTick_9
r, 10, OD_10,WTick_10
r, 11, OD_11,WTick_11
r, 12, OD_12,WTick_12
r, 13, OD_13,WTick_13
r, 14, OD_3, WTick_3,RADCUR_1
r, 15, OD_6, WTick_6,RADCUR_1
r, 16, OD_6, WTick_6,RADCUR_2
r, 17, OD_13, WTick_13,RADCUR_1
r, 18, OD_13, WTick_13,RADCUR_2
r, 19, OD_8, WTick_8,RADCUR_3
r, 20, OD_12, WTick_12,RADCUR_4
r, 21, OD_12, WTick_12,RADCUR_5
r, 22, OD_12, WTick_12,RADCUR_6
r,23,1.69306
r,24,5.07505
r,25,4.96894
r,26,1.20212
r,27,1.42495
r,28,1.88768
r,29,2.18323
r,30,2.4397
r,31,2.98188
r,32,1.41874
r,33,1.04943e+1
r,34,9.30124e-1
r,35,1.6118
r,36,6.74431e-1
r,37,6.43116e-1
r,38,1.06962
r,39,1.20549
r,40,1.05642
r,41,1.25388
r,42,1.3543
r,43,6.66149e-1
r,44,2.27769
r,45,1.15217
r,46,1.23214
r,47,1.52976
r,48, 0.1e+9
r,49, 0.1e+9
r,50, 0.1e+9
r,51, 0.1080e+4
r,52, 0.6001e+5
r,53, 0.6001e+5
r,54, 0.6001e+5
r,55, 0.7541e+6
r,56, 0.7541e+6
r,57, 0.6001e+5
r,58, 0.6000e+3
r,59, 0.6001e+5
r,60, 0.7601e+5
r,61, 0.6001e+5
r,62, 0.8000e+3
r,63, 0.6001e+5
r,64, 0.1000e+9
r,65, 0.1000e+9
r,66, 0.1000e+9
r,67, 0.2600e+3
r,68, 0.5901e+5
r,69, 0.2400e+5
r,70, 0.7601e+5
r,71, 0.2801e+5
r,72, 0.2460e+6
r,73, 0.6001e+5
r,74, 0.6001e+5
r,75, 0.7501e+5
r,76, 0.4660e+6
r,77, 0.3400e+3
r,78, 0.6001e+5
r,79, 0.5000e+6
r,80, 0.5000e+6
r,81, 0.1000e+9
r,82, 0.1000e+9
r,83, 0.1000e+9
/com,------------------------------------------------------------------------------------
/com, Material Properties
/com,*********************
MP,EX, 1, YoungModulus
MP,NUXY,1, Nu
MP,GXY ,1, ShearModulus
MP,ALPX,1, K
MP,EX, 2, YoungModulus
MP,NUXY,2, Nu
MP,GXY ,2, ShearModulus
MP,ALPX,2, K
MP,EX, 3, YoungModulus
MP,NUXY,3, Nu
MP,GXY ,3, ShearModulus
MP,ALPX,3, K
MP,EX, 4, YoungModulus
MP,NUXY,4, Nu
MP,GXY ,4, ShearModulus
MP,ALPX,4, K
MP,EX, 5, YoungModulus
MP,NUXY,5, Nu
MP,GXY ,5, ShearModulus
MP,ALPX,5, K
MP,EX, 6, YoungModulus
MP,NUXY,6, Nu
MP,GXY ,6, ShearModulus
MP,ALPX,6, K
MP,EX, 7, YoungModulus
MP,NUXY,7, Nu
MP,GXY ,7, ShearModulus
MP,ALPX,7, K
MP,EX, 8, YoungModulus
MP,NUXY,8, Nu
MP,GXY ,8, ShearModulus
MP,ALPX,8, K
MP,EX, 9, YoungModulus
MP,NUXY,9, Nu
MP,GXY ,9, ShearModulus
MP,ALPX,9, K
MP,EX, 10, YoungModulus
MP,NUXY,10, Nu
MP,GXY ,10, ShearModulus
MP,ALPX,10, K
MP,EX, 11, YoungModulus
MP,NUXY,11, Nu
MP,GXY ,11, ShearModulus
MP,ALPX,11, K
MP,EX, 12, YoungModulus
MP,NUXY,12, Nu
MP,GXY ,12, ShearModulus
MP,ALPX,12, K
MP,EX, 13, YoungModulus
MP,NUXY,13, Nu
MP,GXY ,13, ShearModulus
MP,ALPX,13, K
MP,EX, 14, YoungModulus
MP,NUXY,14, Nu
MP,GXY ,14, ShearModulus
MP,ALPX,14, K
MP,EX, 15, YoungModulus
MP,NUXY,15, Nu
MP,GXY ,15, ShearModulus
MP,ALPX,15, K
MP,EX, 16, YoungModulus
MP,NUXY,16, Nu
MP,GXY ,16, ShearModulus
MP,ALPX,16, K
MP,EX, 17, YoungModulus
MP,NUXY,17, Nu
MP,GXY ,17, ShearModulus
MP,ALPX,17, K
MP,EX, 18, YoungModulus
MP,NUXY,18, Nu
MP,GXY ,18, ShearModulus
MP,ALPX,18, K
MP,EX, 19, YoungModulus
MP,NUXY,19, Nu
MP,GXY ,19, ShearModulus
MP,ALPX,19, K
MP,EX, 20, YoungModulus
MP,NUXY,20, Nu
MP,GXY ,20, ShearModulus
MP,ALPX,20, K
MP,EX, 21, YoungModulus
MP,NUXY,21, Nu
MP,GXY ,21, ShearModulus
MP,ALPX,21, K
MP,EX, 22, YoungModulus
MP,NUXY,22, Nu
MP,GXY ,22, ShearModulus
MP,ALPX,22, K
/com,------------------------------------------------------------------------------------
/com, Nodes
/com,*******
n,1,-203.808,4715.952,139.224
n,2,-215.808,4715.952,139.224
n,3,-203.808,4703.952,139.224
n,4,-203.808,4715.952,127.224
n,5,-203.268,4719.852,154.860
n,6,-202.308,4726.800,182.748
n,7,-202.008,4728.948,191.340
n,8,-201.876,4729.392,194.916
n,9,-201.888,4729.392,194.988
n,10,-205.908,4729.392,205.740
n,11,-210.036,4729.392,210.168
n,12,-212.160,4729.392,212.436
n,13,-213.324,4729.392,213.684
n,14,-220.968,4729.392,221.868
n,15,-220.968,4734.492,221.868
n,16,-231.888,4729.392,233.568
n,17,-233.952,4729.392,235.776
n,18,-233.952,4741.392,235.776
n,19,-235.980,4729.392,237.936
n,20,-246.204,4714.392,248.904
n,21,-246.204,4705.860,248.904
n,22,-246.204,4686.792,248.904
n,23,-246.204,4676.952,248.904
n,24,-246.204,4665.900,248.904
n,25,-246.204,4638.876,248.904
n,26,-246.204,4631.628,248.904
n,27,-246.204,4628.952,248.904
n,28,-256.812,4613.952,238.296
n,29,-263.004,4613.952,232.104
n,30,-263.004,4625.952,232.104
n,31,-264.048,4613.952,231.060
n,32,-275.352,4613.952,219.744
n,33,-275.352,4619.052,219.744
n,34,-286.668,4613.952,208.440
n,35,-291.996,4613.952,203.112
n,36,-302.604,4613.952,198.720
n,37,-314.988,4613.952,198.720
n,38,-314.988,4613.952,186.720
n,39,-340.380,4613.952,198.720
n,40,-366.000,4613.952,198.720
n,41,-381.000,4628.952,198.720
n,42,-381.000,4630.560,198.720
n,43,-381.000,4655.952,198.720
n,44,-369.035,4655.952,197.800
n,45,-381.000,4689.588,198.720
n,46,-381.000,4756.872,198.720
n,47,-381.000,4760.952,198.720
n,48,-396.000,4775.952,198.720
n,49,-399.000,4775.952,198.720
n,50,-399.000,4787.952,198.720
n,51,-402.600,4775.952,198.720
n,52,-402.600,4775.872,210.685
n,53,-440.904,4775.952,198.720
n,54,-450.612,4775.952,198.720
n,55,-465.612,4790.952,198.720
n,56,-465.612,4796.592,198.720
n,57,-453.643,4796.592,197.860
n,58,-465.612,4840.068,198.720
n,59,-465.612,4927.032,198.720
n,60,-465.612,4943.952,198.720
n,61,-459.864,4958.952,212.568
n,62,-458.712,4958.952,215.340
n,63,-458.712,4970.952,215.340
n,64,-455.748,4958.952,222.480
n,65,-463.860,4958.952,242.076
n,66,-484.044,4958.952,250.452
n,67,-533.124,4958.952,270.828
n,68,-537.724,4958.952,259.745
n,69,-552.516,4958.952,278.880
n,70,-552.600,4958.952,278.916
n,71,-560.652,4958.952,298.392
n,72,-559.464,4958.952,301.248
n,73,-559.464,4970.952,301.248
n,74,-558.312,4958.952,304.020
n,75,-553.866,4961.957,293.286
n,76,-554.868,4958.952,312.324
n,77,-548.736,4958.952,327.108
n,78,-548.736,4981.752,327.108
n,79,-546.708,4958.952,331.992
n,80,-542.604,4958.952,341.880
n,81,-538.152,4958.952,352.620
n,82,-532.404,4973.952,366.468
n,83,-532.404,4988.628,366.468
n,84,-532.404,5039.952,366.468
n,85,-532.404,5040.048,366.468
n,86,-546.168,5054.952,372.180
n,87,-551.904,5054.952,374.568
n,88,-551.904,5066.952,374.568
n,89,-560.988,5054.952,378.336
n,90,-560.988,5066.952,378.336
n,91,-565.128,5054.952,380.052
n,92,-578.988,5069.952,385.812
n,93,-578.988,5091.924,385.812
n,94,-578.988,5141.940,385.812
n,95,-590.988,5141.940,385.812
n,96,-578.988,5129.940,385.812
n,97,-578.988,5141.940,373.812
n,98,-242.400,4676.952,252.708
n,99,-242.382,4676.952,252.720
n,100,-238.434,4676.952,256.680
n,101,-232.056,4667.952,263.052
n,102,-232.056,4665.936,263.052
n,103,-232.056,4622.952,263.052
n,104,-225.696,4613.952,256.692
n,105,-218.688,4613.952,249.684
n,106,-216.048,4613.952,243.312
n,107,-216.048,4613.952,208.800
n,108,-216.048,4613.952,181.644
n,109,-216.048,4625.952,181.644
n,110,-216.048,4613.952,155.748
n,111,-216.048,4613.952,103.944
n,112,-216.048,4613.952,78.048
n,113,-204.048,4613.952,78.048
n,114,-216.048,4613.952,34.980
n,115,-216.048,4613.952,27.048
n,116,-225.048,4613.952,18.048
n,117,-231.048,4613.952,18.048
n,118,-231.048,4616.447,29.785
n,119,-246.048,4613.952,18.048
n,120,-246.048,4625.952,18.048
n,121,-267.048,4613.952,18.048
n,122,-276.048,4613.952,27.048
n,123,-276.048,4613.952,40.452
n,124,-276.048,4613.952,137.544
n,125,-276.048,4613.952,168.948
n,126,-276.048,4622.952,177.948
n,127,-276.048,4625.952,177.948
n,128,-274.814,4625.952,189.884
n,129,-276.048,4668.492,177.948
n,130,-276.048,4753.560,177.948
n,131,-276.048,4766.952,177.948
n,132,-285.048,4775.952,177.948
n,133,-288.048,4775.952,177.948
n,134,-288.048,4787.952,177.948
n,135,-300.048,4775.952,177.948
n,136,-300.048,4781.441,188.620
n,137,-350.532,4775.952,177.948
n,138,-451.512,4775.952,177.948
n,139,-472.020,4775.952,177.948
n,140,-481.020,4784.952,177.948
n,141,-481.020,4800.792,177.948
n,142,-469.580,4800.792,181.571
n,143,-481.020,4855.320,177.948
n,144,-481.020,4964.388,177.948
n,145,-481.020,4979.592,177.948
n,146,-481.020,4979.592,189.948
n,147,-481.020,5000.592,177.948
n,148,-489.336,5009.952,181.392
n,149,-492.876,5009.952,182.868
n,150,-492.876,5021.952,182.868
n,151,-520.044,5009.952,194.148
n,152,-558.600,5009.952,210.156
n,153,-566.916,5018.952,213.600
n,154,-566.916,5021.892,213.600
n,155,-566.916,5048.952,213.600
n,156,-563.472,5057.952,221.912
n,157,-562.320,5057.952,224.688
n,158,-562.320,5069.952,224.688
n,159,-558.864,5057.952,233.004
n,160,-556.716,5057.952,238.188
n,161,-553.884,5057.952,245.004
n,162,-553.884,5062.044,245.004
n,163,-548.904,5057.952,257.016
n,164,-546.600,5057.952,262.560
n,165,-545.436,5057.952,265.368
n,166,-556.148,5060.878,260.821
n,167,-534.588,5057.952,291.528
n,168,-512.880,5057.952,343.824
n,169,-502.020,5057.952,369.984
n,170,-502.020,5069.952,369.984
n,171,-500.988,5057.952,372.468
n,172,-496.392,5045.952,383.544
n,173,-496.392,5037.192,383.544
n,174,-496.392,5033.952,383.544
n,175,-506.700,5021.952,389.676
n,176,-539.784,5021.952,409.358
n,177,-596.976,5059.452,495.504
n,178,-588.540,5021.952,438.372
n,179,-588.624,5021.952,438.420
n,180,-595.416,5029.848,442.464
n,181,-595.416,5035.452,442.464
n,182,-605.748,5047.452,448.596
n,183,-610.896,5047.452,451.668
n,184,-615.072,5047.452,468.120
n,185,-613.716,5047.452,470.388
n,186,-603.960,5047.452,486.792
n,187,-597.504,5047.452,494.988
n,188,-596.976,5047.452,495.504
n,189,-570.048,5047.452,521.736
n,190,-516.180,5047.452,574.200
n,191,-462.312,5047.452,626.676
n,192,-435.384,5047.452,652.920
n,193,-447.384,5047.452,652.920
n,194,-435.384,5035.452,652.920
n,195,-435.384,5047.452,640.920
n,271,-249.31,4618.4,245.79
n,272,-560.66,4959,286.99
n,273,-574.94,5059.4,384.11
n,274,-536.41,5050.6,368.12
n,275,-233.93,4674.4,261.18
n,276,-483.33,5007.2,178.89
n,277,-564.48,5012.6,212.59
n,278,-598.45,5044.0,444.25
n,279,-616.39,5047.5,459.03
/com,------------------------------------------------------------------------------------
/com,
/com, Straight Pipe (Tangent) Elements
/com,**********************************
mat,1
type,1
real,1
e,1,5
type,2
real,2
e,5,6
type,3
real,3
e,6,7
e,8,9
e,10,11
type,4
real,4
e,11,12
e,12,13
e,13,14
e,14,16
e,31,32
e,32,34
type,5
real,5
e,14,15
e,32,33
e,77,78
type,6
real,6
e,16,17
e,17,19
e,20,21
e,21,22
e,22,23
e,23,24
e,24,25
e,25,26
e,26,27
e,28,29
e,29,31
e,34,35
e,36,37
e,37,39
e,39,40
e,41,42
e,42,43
e,43,45
e,45,46
e,46,47
e,48,49
e,49,51
e,51,53
e,53,54
e,55,56
e,56,58
e,58,59
e,59,60
e,61,62
e,62,64
e,65,66
e,66,67
e,67,69
e,69,70
e,71,72
e,72,74
e,74,76
type,7
real,7
e,76,77
e,77,79
e,79,80
type,13
real,13
e,80,81
e,82,83
e,83,84
e,84,85
e,86,87
e,87,89
e,89,91
e,92,93
e,93,94
id = 8
type,8
real,8
e,23,98
e,98,99
e,99,100
e,101,102
e,102,103
e,104,105
e,106,107
e,107,108
e,108,110
e,110,111
e,111,112
e,112,114
e,114,115
e,116,117
e,117,119
e,119,121
e,122,123
e,123,124
e,124,125
e,126,127
e,127,129
e,129,130
e,130,131
e,132,133
e,133,135
e,135,137
e,137,138
e,138,139
e,140,141
e,141,143
e,143,144
e,144,145
e,145,147
e,148,149
e,149,151
e,151,152
e,153,154
e,154,155
e,156,157
e,157,159
e,163,164
type,9
real,9
e,159,160
e,160,161
e,161,163
type,10
real,10
e,161,162
type,12
real,12
e,164,165
e,165,167
e,167,168
e,168,169
e,169,171
e,172,173
e,173,174
e,175,176
e,176,178
e,178,179
e,180,181
e,182,183
e,184,185
e,185,186
e,187,188
e,188,189
e,189,190
e,190,191
e,191,192
/com,
/com, Pipe Bend Elements
/com,********************
mat,1
type,14
real,14
e,7,8,6
e,9,10,11
type,15
real,15
e,19,20,21
e,27,271,26
e,271,28,29
e,35,36,37
e,40,41,42
e,47,48,49
e,54,55,53
e,60,61,62
e,64,65,66
type,16
real,16
e,70,272,69
e,272,71,72
type,17
real,17
e,81,82,83
e,91,273,89
e,273,92,93
type,18
real,18
e,85,274,83
e,274,86,87
type,19
real,19
e,100,275,99
e,275,101,102
e,103,104,105
e,105,106,107
e,115,116,114
e,121,122,123
e,125,126,124
e,131,132,133
e,139,140,138
e,147,276,145
e,276,148,149
e,152,277,151
e,277,153,154
e,155,156,157
type,20
real,20
e,171,172,173
e,174,175,176
e,181,278,180
e,278,182,183
e,183,279,182
e,279,184,185
type,21
real,21
e,186,187,188
type,22
real,22
e,179,180,181
/com, Point Mass without considering rotatary inertia
/com,*************************************************
type,23
real,23
e,12
e,21
e,25
type,24
real,24
e,15
type,25
real,25
e,33
type,26
real,26
e,34
type,27
real,27
e,39
e,42
type,28
real,28
e,45
e,46
type,29
real,29
e,53
type,30
real,30
e,58
e,59
type,31
real,31
e,66
type,32
real,32
e,69
type,33
real,33
e,78
type,34
real,34
e,79
e,83
e,84
type,35
real,35
e,93
type,36
real,36
e,102
e,104
e,107
type,37
real,37
e,110
e,111
type,38
real,38
e,114
type,39
real,39
e,123
e,124
type,40
real,40
e,129
e,130
type,41
real,41
e,137
e,138
type,42
real,42
e,143
e,144
type,43
real,43
e,151
e,154
e,160
type,44
real,44
e,162
type,45
real,45
e,167
e,168
type,46
real,46
e,173
e,176
e,178
e,185
type,47
real,47
e,189
e,190
e,191
/com,
/com, Elastic supports and anchors
/com,******************************
type,48
real,48
e,1,2
type,49
real,49
e,1,3
type,50
real,50
e,1,4
id = 51
type,51
real,51
e,17,18
type,52
real,52
e,29,30
type,53
real,53
e,37,38
type,55
real,55
e,49,50
type,58
real,58
e,62,63
type,60
real,60
e,72,73
type,62
real,62
e,87,88
type,63
real,63
e,89,90
type,64
real,64
e,94,95
type,65
real,65
e,94,96
type,66
real,66
e,94,97
type,67
real,67
e,108,109
type,68
real,68
e,112,113
type,70
real,70
e,119,120
type,72
real,72
e,133,134
type,75
real,75
e,145,146
type,76
real,76
e,149,150
type,77
real,77
e,157,158
id = 79
type,79
real,79
e,169,170
type,80
real,80
e,188,177
type,81
real,81
e,192,193
type,82
real,82
e,192,194
type,83
real,83
e,192,195
/com,
/com, rotate nodes with less than 3 supports
/com,****************************************
wplane,,nx(43),ny(43),nz(43),nx(44),ny(44),nz(44),nx(45),ny(45),nz(45)
cswplane,11,0
nrotat,43
nrotat,44
csys,0
real,54
type,54
e,43,44
wplane,,nx(51),ny(51),nz(51),nx(52),ny(52),nz(52),nx(49),ny(49),nz(49)
cswplane,12,0
nrotat,51
nrotat,52
csys,0
real,56
type,56
e,51,52
wplane,,nx(56),ny(56),nz(56),nx(57),ny(57),nz(57),nx(55),ny(55),nz(55)
cswplane,13,0
nrotat,56
nrotat,57
csys,0
real,57
type,57
e,56,57
wplane,,nx(67),ny(67),nz(67),nx(68),ny(68),nz(68),nx(66),ny(66),nz(66)
cswplane,14,0
nrotat,67
nrotat,68
csys,0
real,59
type,59
e,67,68
wplane,,nx(74),ny(74),nz(74),nx(75),ny(75),nz(75),nx(72),ny(72),nz(72)
cswplane,15,0
nrotat,74
nrotat,75
csys,0
real,61
type,61
e,74,75
wplane,,nx(117),ny(117),nz(117),nx(118),ny(118),nz(118),nx(116),ny(116),nz(116)
cswplane,16,0
nrotat,117
nrotat,118
csys,0
real,69
type,69
e,117,118
wplane,,nx(127),ny(127),nz(127),nx(128),ny(128),nz(128),nx(126),ny(126),nz(126)
cswplane,17,0
nrotat,127
nrotat,128
csys,0
real,71
type,71
e,127,128
wplane,,nx(135),ny(135),nz(135),nx(136),ny(136),nz(136),nx(133),ny(133),nz(133)
cswplane,18,0
nrotat,135
nrotat,136
csys,0
real,73
type,73
e,135,136
wplane,,nx(141),ny(141),nz(141),nx(142),ny(142),nz(142),nx(140),ny(140),nz(140)
cswplane,19,0
nrotat,141
nrotat,142
csys,0
real,74
type,74
e,141,142
wplane,,nx(165),ny(165),nz(165),nx(166),ny(166),nz(166),nx(164),ny(164),nz(164)
cswplane,20,0
nrotat,165
nrotat,166
csys,0
real,78
type,78
e,165,166
/com,------------------------------------------------------------------------------------
/com,
/com, Constraints
/com,*************
nsel,s,node,,2,4
nsel,a,node,,18
nsel,a,node,,30
nsel,a,node,,38
nsel,a,node,,44
nsel,a,node,,50
nsel,a,node,,52
nsel,a,node,,57
nsel,a,node,,63
nsel,a,node,,68
nsel,a,node,,73
nsel,a,node,,75
nsel,a,node,,88
nsel,a,node,,90
nsel,a,node,,95,97
nsel,a,node,,109
nsel,a,node,,113
nsel,a,node,,118
nsel,a,node,,120
nsel,a,node,,128
nsel,a,node,,134
nsel,a,node,,136
nsel,a,node,,142
nsel,a,node,,146
nsel,a,node,,150
nsel,a,node,,158
nsel,a,node,,166
nsel,a,node,,170
nsel,a,node,,177
nsel,a,node,,193,195
d,all,all
allsel,all,all
d,1,rotx,,,,,roty,rotz
d,94,rotx,,,,,roty,rotz
d,192,rotx,,,,,roty,rotz
finish
/com,------------------------------------------------------------------------------------
/com,
/com,=============
/com, Modal Solve
/com,=============
/com,
/solution
antype,modal ! Perform Modal Analysis
modopt,lanb,maxm
lumpm,on ! Use Lumped Mass Matrix Approximation
mxpand,maxm,,,yes ! Expand solutions with Element Calculations turned ON
solve
save
/com,
/com,===========================
/com, Compare Modal Frequencies
/com,===========================
/com,
*dim,Amode,ARRAY,maxm
*dim,Emode,ARRAY,maxm
*dim,ERmode,ARRAY,maxm
*dim,moden,ARRAY,maxm
*do,i,1,maxm
*GET, Amode(i), MODE, i, FREQ
*enddo
*VFILL,Emode,DATA,2.612,2.914,4.337,4.66,5.734,5.833,7.359,7.769,9.952,10.329
*VFILL,Emode(11),DATA,10.679,10.943,12.03,12.286,13.251,13.407,14.429,14.72,15.253,15.553
*VFILL,Emode(21),DATA,16.172,16.797,17.23,17.275,17.453,18.71,18.898,19.993,21.46,21.523
*VFILL,Emode(31),DATA,22.736,23.281,24.067,24.593,25.117,26.516,26.935,27.509,28.662,29.542
*VFILL,Emode(41),DATA,30.596,31.274,32.283,35.484,36.022,36.394,36.769,38.0,38.42,40.185
*do,i,1,maxm
ERmode(i) = ABS(Amode(i)/Emode(i))
moden(i) = i
*enddo
save,table_1
finish
/com,------------------------------------------------------------------------------------
/com,
/com,================
/com, Spectrum Solve
/com,================
/com,
/solution
antype,spectr ! Perform Spectrum Analysis
spopt,sprs,maxm ! Single Point Response Spectrum
SRSS,0.0 ! Combine modes using SRSS mode combination
/com,------------------------------------------------------------------------------------
/com,
/com,
/com, spectrum 1 (X - Direction)
/com,****************************
svtyp, 2, 386.20
freq,0.2,0.5,0.5263,0.578,0.7752,1.7241,2.0,2.1978,2.2472,
SV,,0.16,0.16,0.16,0.23,0.339,0.55,0.555,0.555,0.65,
freq,2.9851,3.125,4.1667,5.0,5.7143,7.6923,9.0909,9.5238,10.4167
sv,,0.65,1.5,1.5,1.36,0.45,0.935,0.935,0.995,1.31,
freq,14.0845,17.8571,19.6078,27.027,33.3333,50.0,100.0
sv,,1.31,1.0,1.222,1.222,0.5,0.252,0.252,
sed,1,0,0
solve
/com,
/com, spectrum 2 (Y - Direction)
/com,****************************
svtyp, 2, 386.20
freq,,
FREQ,0.2,0.6061,0.667,1.9231,2.439,3.7037,4.0,5.3476,5.5866,
SV,,0.085,0.115,0.125,0.36,0.40,0.525,0.54,0.54,0.506,
FREQ,8.333,10.989,12.5,14.7059,20.0,21.2766,22.2222,27.7778,28.5714,
SV,,2.6,2.6,0.9,0.875,0.875,0.76,0.85,0.85,0.522,
FREQ,30.303,33.3333,50.0,100.0
SV,,0.255,0.295,0.194,0.194
sed,0,1,0
solve
/com,
/com, spectrum 3 (Z - Direction)
/com,***************************
svtyp, 2, 386.20
freq,,,
FREQ,0.2,0.5,0.5814,1.7544,2.1739,2.8571,3.0303,4.1667,5.0,
SV,,0.16,0.16,0.23,0.56,0.56,0.80,1.5,1.5,1.15,
FREQ,5.2632,5.7143,5.8824,6.2893,6.5789,7.4074,9.6154,10.4167,14.0845,
SV,,0.64,0.64,0.7,0.7,0.875,1.05,1.05,1.31,1.31,
FREQ,15.873,17.2414,25.0,27.027,32.2581,50.0,100.0
SV,,0.85,1.15,1.15,1.1,0.5,0.225,0.225,
sed,0,0,1
solve
finish
/com,------------------------------------------------------------------------------------
/com,
/post1
/input,,mcom
/com,-----------------------------------
/com, *Labels*
*dim,label2,char,1,6
*dim,label3,char,6,1
*dim,label4,char,6,1
*dim,label5,char,22,1
/com,-------------------------
label2(1,1) = 'ux_81'
label2(1,2) = 'uy_155'
label2(1,3) = 'uz_61'
label2(1,4) ='rotx_143'
label2(1,5) ='roty_149'
label2(1,6) ='rotz_84'
/com,-----------------------
label3(1,1)='PX(I)'
label3(2,1)='VY(I)'
label3(3,1)='VZ(I)'
label3(4,1)='TX(I)'
label3(5,1)='MY(I)'
label3(6,1)='MZ(I)'
/com,-----------------------
label4(1,1)='PX(J)'
label4(2,1)='VY(J)'
label4(3,1)='VZ(J)'
label4(4,1)='TX(J)'
label4(5,1)='MY(J)'
label4(6,1)='MZ(J)'
/com,-----------------------------------
/com,------------------------------------------------------------------------------------
/com,
/com,========================================================
/com, Maximum nodal displacements and rotations comparsion
/com,========================================================
/com,
/com, Solution obtained from Mechanical APDL
/com, ****************************
*GET,AdisX,NODE,81,U,X
*GET,AdisY,NODE,155,U,Y
*GET,AdisZ,NODE,61,U,Z
*GET,ArotX,NODE,143,ROT,X
*GET,ArotY,NODE,149,ROT,Y
*GET,ArotZ,NODE,84,ROT,Z
/com,
/com, Expected results from NRC manual
/com, *********************************
*SET,EdisX,9.29050e-01
*SET,EdisY,3.19026e-01
*SET,EdisZ,6.18342e-01
*SET,ErotX,5.98220e-03
*SET,ErotY,1.00078e-02
*SET,ErotZ,9.17224e-03
/com,
/com, Error computation
/com,********************
ERdisX=ABS(AdisX/EdisX)
ERdisY=ABS(AdisY/EdisY)
ERdisZ=ABS(AdisZ/EdisZ)
ERrotX=ABS((ArotX)/(ErotX))
ERrotY=ABS((ArotY)/(ErotY))
ERrotZ=ABS((ArotZ)/(ErotZ))
*dim,value,,6,3
*vfill,value(1,1),data,EdisX
*vfill,value(1,2),data,AdisX
*vfill,value(1,3),data,ERdisX
*vfill,value(2,1),data,EdisY
*vfill,value(2,2),data,AdisY
*vfill,value(2,3),data,ERdisY
*vfill,value(3,1),data,EdisZ
*vfill,value(3,2),data,AdisZ
*vfill,value(3,3),data,ERdisZ
*vfill,value(4,1),data,ErotX
*vfill,value(4,2),data,ArotX
*vfill,value(4,3),data,ERrotX
*vfill,value(5,1),data,ErotY
*vfill,value(5,2),data,ArotY
*vfill,value(5,3),data,ERrotY
*vfill,value(6,1),data,ErotZ
*vfill,value(6,2),data,ArotZ
*vfill,value(6,3),data,ERrotZ
save,table_2
/com,------------------------------------------------------------------------------------
/com,========================================================
/com, Element Forces and Moments Comparison
/com,========================================================
/com, Solution obtained from Mechanical APDL
/com,******************************
*dim,elem_res_I,,1,6
*dim,elem_res_J,,1,6
*dim,pxi,,1
*dim,vyi,,1
*dim,vzi,,1
*dim,txi,,1
*dim,myi,,1
*dim,mzi,,1
*dim,pxj,,1
*dim,vyj,,1
*dim,vzj,,1
*dim,txj,,1
*dim,myj,,1
*dim,mzj,,1
esel,s,ename,,16
esel,a,ename,,18
/com,==========
/com, Node I
/com,==========
/com, Element #1
/com,***********
*get,pxi(1,1),elem,1,smisc,1
*get,vyi(1,1),elem,1,smisc,2
*get,vzi(1,1),elem,1,smisc,3
*get,txi(1,1),elem,1,smisc,4
*get,myi(1,1),elem,1,smisc,5
*get,mzi(1,1),elem,1,smisc,6
*vfill,elem_res_I(1,1),data,pxi(1,1)
*vfill,elem_res_I(1,2),data,vyi(1,1)
*vfill,elem_res_I(1,3),data,vzi(1,1)
*vfill,elem_res_I(1,4),data,txi(1,1)
*vfill,elem_res_I(1,5),data,myi(1,1)
*vfill,elem_res_I(1,6),data,mzi(1,1)
/com,==========
/com, Node J
/com,==========
/com, Element #1
/com,************
*get,pxj(1,1),elem,1,smisc,7
*get,vyj(1,1),elem,1,smisc,8
*get,vzj(1,1),elem,1,smisc,9
*get,txj(1,1),elem,1,smisc,10
*get,myj(1,1),elem,1,smisc,11
*get,mzj(1,1),elem,1,smisc,12
*vfill,elem_res_J(1,1),data,pxj(1,1)
*vfill,elem_res_J(1,2),data,vyj(1,1)
*vfill,elem_res_J(1,3),data,vzj(1,1)
*vfill,elem_res_J(1,4),data,txj(1,1)
*vfill,elem_res_J(1,5),data,myj(1,1)
*vfill,elem_res_J(1,6),data,mzj(1,1)
/com,----------------------------------------------------------------------------
/com, Results from NRC benchmarks
/com, ***************************
*dim,exp_I,,1,6
*dim,exp_J,,1,6
/com, Element #1
/com,************
*vfill,exp_I(1,1),data,2.273e+03
*vfill,exp_I(1,2),data,3.719e+03
*vfill,exp_I(1,3),data,2.287e+03
*vfill,exp_I(1,4),data,1.050e+05
*vfill,exp_I(1,5),data,2.299e+05
*vfill,exp_I(1,6),data,3.516e+05
*vfill,exp_J(1,1),data,2.273e+03
*vfill,exp_J(1,2),data,3.719e+03
*vfill,exp_J(1,3),data,2.287e+03
*vfill,exp_J(1,4),data,1.050e+05
*vfill,exp_J(1,5),data,1.987e+05
*vfill,exp_J(1,6),data,2.920e+05
/com,---------------------------------------------------------------------------
/com, Error computation
/com, *********************
*dim,elem_error_I,,1,6
*dim,elem_error_J,,1,6
*dim,elem_tab,,12,3
/com,============
/com, Node I
/com,============
i = 1
*do,j,1,6
*vfill,elem_error_I(i,j),data,abs(elem_res_I(i,j)/exp_I(i,j))
*enddo
/com,============
/com, Node J
/com,============
i = 1
*do,j,1,6
*vfill,elem_error_J(i,j),data,abs(elem_res_J(i,j)/exp_J(i,j))
*enddo
/com,--------------------------------------------------------------------------
i = 1
cs=(i-1)*6
*do,j,1,6
n=cs+j
*vfill,elem_tab(n,1),data,exp_I(i,j)
*vfill,elem_tab(n,2),data,elem_res_I(i,j)
*vfill,elem_tab(n,3),data,elem_error_I(i,j)
*enddo
*do,j,1,6
m=cs+j+6
*vfill,elem_tab(m,1),data,exp_J(i,j)
*vfill,elem_tab(m,2),data,elem_res_J(i,j)
*vfill,elem_tab(m,3),data,elem_error_J(i,j)
*enddo
save,table_3
/com,------------------------------------------------------------------------------------
/com,*****************************
/com, Reaction forces comparision
/com,******************************
*dim,Areac,,36
*dim,Ereac,,36
*dim,ERreac,,36
*dim,Nreac,STRING,8,36
/com, Labels
/com,********
Nreac(1,1) = 'FX1'
Nreac(1,2) = 'FY1'
Nreac(1,3) = 'FZ1'
Nreac(1,4) = 'FY17'
Nreac(1,5) = 'FY29'
Nreac(1,6) = 'FZ37'
Nreac(1,7) = 'FX43'
Nreac(1,8) = 'FY49'
Nreac(1,9) = 'FZ51'
Nreac(1,10) = 'FX56'
Nreac(1,11) = 'FY62'
Nreac(1,12) = 'FZ67'
Nreac(1,13) = 'FY72'
Nreac(1,14) = 'FZ74'
Nreac(1,15) = 'FY87'
Nreac(1,16) = 'FY89'
Nreac(1,17) = 'FX94'
Nreac(1,18) = 'FY94'
Nreac(1,19) = 'FZ94'
Nreac(1,20) = 'FY108'
Nreac(1,21) = 'FX112'
Nreac(1,22) = 'FZ117'
Nreac(1,23) = 'FY119'
Nreac(1,24) = 'FZ127'
Nreac(1,25) = 'FY133'
Nreac(1,26) = 'FZ135'
Nreac(1,27) = 'FX141'
Nreac(1,28) = 'FZ145'
Nreac(1,29) = 'FY149'
Nreac(1,30) = 'FY157'
Nreac(1,31) = 'FX165'
Nreac(1,32) = 'FY169'
Nreac(1,33) = 'FY188'
Nreac(1,34) = 'FX192'
Nreac(1,35) = 'FY192'
Nreac(1,36) = 'FZ192'
/com, Solution obtained from Mechanical APDL
/com,******************************
*GET,Areac(1),NODE,2,RF,FX
*GET,Areac(2),NODE,3,RF,FY
*GET,Areac(3),NODE,4,RF,FZ
*GET,Areac(4),NODE,18,RF,FY
*GET,Areac(5),NODE,30,RF,FY
*GET,Areac(6),NODE,38,RF,FZ
*GET,Areac(7),NODE,44,RF,FX
*GET,Areac(8),NODE,50,RF,FY
*GET,Areac(9),NODE,52,RF,FX
*GET,Areac(10),NODE,57,RF,FX
*GET,Areac(11),NODE,63,RF,FY
*GET,Areac(12),NODE,68,RF,FX
*GET,Areac(13),NODE,73,RF,FY
*GET,Areac(14),NODE,75,RF,FX
*GET,Areac(15),NODE,88,RF,FY
*GET,Areac(16),NODE,90,RF,FY
*GET,Areac(17),NODE,95,RF,FX
*GET,Areac(18),NODE,96,RF,FY
*GET,Areac(19),NODE,97,RF,FZ
*GET,Areac(20),NODE,109,RF,FY
*GET,Areac(21),NODE,113,RF,FX
*GET,Areac(22),NODE,118,RF,FX
*GET,Areac(23),NODE,120,RF,FY
*GET,Areac(24),NODE,128,RF,FX
*GET,Areac(25),NODE,134,RF,FY
*GET,Areac(26),NODE,136,RF,FX
*GET,Areac(27),NODE,142,RF,FX
*GET,Areac(28),NODE,146,RF,FZ
*GET,Areac(29),NODE,150,RF,FY
*GET,Areac(30),NODE,158,RF,FY
*GET,Areac(31),NODE,166,RF,FX
*GET,Areac(32),NODE,170,RF,FY
*GET,Areac(33),NODE,177,RF,FY
*GET,Areac(34),NODE,193,RF,FX
*GET,Areac(35),NODE,194,RF,FY
*GET,Areac(36),NODE,195,RF,FZ
/com, Expected results from NRC manual
/com,**********************************
*VFILL,Ereac,DATA,3724,2390,2156,42,2466,4850,4765,3835,3482,2101
*VFILL,Ereac(11),DATA,61,6860,2669,6554,109,5015,3334,4739,861,64
*VFILL,Ereac(21),DATA,2312,2079,1153,1829,886,889,1858,2571,1349,106
*VFILL,Ereac(31),DATA,4370,1340,1170,970,749,2952
/com, Error computation
/com,*******************
*do,i,1,36
ERreac(i) = abs(Areac(i)/Ereac(i))
*enddo
save,table_4
/com,------------------------------------------------------------------------------------
/com,
/out,
/com,
/com, ------------------------vm-nr1677-2-4a-a Results Verification-------------------------
/com,
/nopr
resume,table_1
/gopr
/out,vm-nr1677-2-4a-a,vrt
/com,
/com, ===========================================
/com, COMPARISON OF MODAL FREQUENCY
/com, WITH EXPECTED RESULTS
/com, ===========================================
/com,
/com, Mode | Expected | Mechanical APDL | Ratio
/com,
*VWRITE,moden(1),Emode(1),Amode(1),ERmode(1)
(1X,F3.0,2X,F8.4,3X,F8.4,3X,F4.2,' ')
/com,
/com, -------------------------------------------------------------------------
/nopr
resume,table_2
/gopr
/com,
/com,====================================================
/com, COMPARISON OF NODAL DISPLACEMENTS AND ROTATIONS
/com, WITH EXPECTED RESULTS
/com,====================================================
/com,
/com, Result_Node | Expected | Mechanical APDL | Ratio
/com,
*vwrite,label2(1,1),value(1,1),value(1,2),value(1,3)
(1x,a8,' ',f10.4,' ',f10.4,' ',f5.3)
*vwrite,label2(1,2),value(2,1),value(2,2),value(2,3)
(1x,a8,' ',f10.4,' ',f10.4,' ',f5.3)
*vwrite,label2(1,3),value(3,1),value(3,2),value(3,3)
(1x,a8,' ',f10.4,' ',f10.4,' ',f5.3)
*vwrite,label2(1,4),value(4,1),value(4,2),value(4,3)
(1x,a8,' ',f10.4,' ',f10.4,' ',f5.3)
*vwrite,label2(1,5),value(5,1),value(5,2),value(5,3)
(1x,a8,' ',f10.4,' ',f10.4,' ',f5.3)
*vwrite,label2(1,6),value(6,1),value(6,2),value(6,3)
(1x,a8,' ',f10.4,' ',f10.4,' ',f5.3)
/com,
/com,-------------------------------------------------------------------------
/com,
/nopr
resume,table_4
/gopr
/com,
/com, ===========================================
/com, COMPARISON OF REACTION FORCES
/com, WITH EXPECTED RESULTS
/com, ===========================================
/com,
/com, Node | Expected | Mechanical APDL | Ratio
/com,
*VWRITE,Nreac(1),Ereac(1),Areac(1),ERreac(1)
(5X,a,2X,F12.4,3X,F12.4,3X,F8.4,' ')
/com,
/com,--------------------------------------------------------------------------
/com,
/nopr
resume,table_3
/gopr
/com,
/com,===============================================
/com, COMPARISON OF ELEMENT FORCES AND MOMENTS
/com, WITH EXPECTED RESULTS
/com,===============================================
/com,
/com,--------------------------------------------------------
/com, Note: Element Forces and Moments for some elements
/com, along Y & Z directions are flipped between Mechanical APDL
/com, and NRC results
/com,--------------------------------------------------------
/com, Result | Expected | Mechanical APDL | Ratio
/com,
/com,===============
/com, Element 1
/com,===============
/com,
*vwrite,label3(1,1),elem_tab(1,1),elem_tab(1,2),elem_tab(1,3)
(1x,a8,' ',f12.4,' ',f12.4,' ',f5.3)
/com,
*vwrite,label4(1,1),elem_tab(7,1),elem_tab(7,2),elem_tab(7,3)
(1x,a8,' ',f12.4,' ',f12.4,' ',f5.3)
/com,
/com,
/com,*******************************************************************
/com,*******************************************************************
/com,
/com,
/out,
*list,vm-nr1677-2-4a-a,vrt
finish
