VM-NR1677-01-1-a
VM-NR1677-01-1-a
NUREG/CR-1677: Volume 1, Benchmark Problem No. 1
Overview
Test Case
This benchmark problem contains three straight sections, two bends, and two fixed anchors (refer to Figure 619: FE Model of Benchmark Problem). The total mass of the system is represented by structural mass element (MASS21) specified at individual nodes. Modal and response spectrum analysis is performed on the piping model. Frequencies obtained from modal solve and the nodal/element solution obtained from spectrum solve are compared against reference results.
| Material Properties | Geometric Properties | Loading | |||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Pipe Elements:
3D mass (lb-sec2/in):
| Straight Pipe:
Bend Pipe:
| Acceleration response spectrum curve defined by SV and FREQ commands. |
Results Comparison
Table 34: Frequencies Obtained from Modal Solution
| Mode | Target | Mechanical APDL | Ratio |
|---|---|---|---|
| 1 | 28.530 | 28.534 | 1.00 |
| 2 | 55.770 | 55.771 | 1.00 |
| 3 | 81.500 | 81.499 | 1.00 |
| 4 | 141.700 | 141.739 | 1.00 |
| 5 | 162.800 | 162.816 | 1.00 |
Table 35: Maximum Displacements and Rotations Obtained from Spectrum Solve
| Result Node | Target | Mechanical APDL | Ratio |
|---|---|---|---|
| UX at node5 | 0.0078 | 0.0078 | 0.999 |
| UY at node7 | 0.0025 | 0.0025 | 1.00 |
| UZ at node4 | 0.0174 | 0.0174 | 1.00 |
| ROTX at node3 | 0.0002 | 0.0002 | 1.00 |
| ROTY at node7 | 0.0002 | 0.0002 | 1.00 |
| ROTZ at node3 | 0.0001 | 0.0001 | 1.00 |
Table 36: Element Forces and Moments Obtained from Spectrum Solve
| Result | Target | Mechanical APDL | Ratio |
|---|---|---|---|
| Element 1 | |||
| PX(I) | 4.964 | 4.958 | 1.001 |
| VY(I) | 17.879 | 17.880 | 1.000 |
| VZ(I) | 36.436 | 36.430 | 1.000 |
| TX(I) | 629.715 | 629.600 | 1.000 |
| MY(I) | 3226.871 | 3227.000 | 1.000 |
| MZ(I) | 1393.841 | 1394.000 | 1.000 |
| PX(J) | 4.964 | 4.958 | 1.001 |
| VY(J) | 17.879 | 17.880 | 1.000 |
| VZ(J) | 36.436 | 36.430 | 1.000 |
| TX(J) | 629.715 | 629.600 | 1.000 |
| MY(J) | 1259.879 | 1260.000 | 1.000 |
| MZ(J) | 473.244 | 474.200 | 0.998 |
| Element 10 | |||
| PX(I) | 24.022 | 24.020 | 1.000 |
| VY(I) | 7.478 | 7.472 | 1.001 |
| VZ(I) | 34.800 | 34.780 | 1.001 |
| TX(I) | 112.957 | 113.000 | 1.000 |
| MY(I) | 1871.654 | 1871.000 | 1.000 |
| MZ(I) | 650.273 | 650.100 | 1.000 |
| PX(J) | 24.022 | 24.020 | 1.000 |
| VY(J) | 7.478 | 7.472 | 1.001 |
| VZ(J) | 34.800 | 34.780 | 1.001 |
| TX(J) | 112.957 | 113.000 | 1.000 |
| MY(J) | 2477.415 | 2477.000 | 1.000 |
| MZ(J) | 774.763 | 774.500 | 1.000 |
| Element 4 | |||
| PX(I) | 9.304 | 9.300 | 1.000 |
| VY(I) | 10.635 | 10.630 | 1.000 |
| VZ(I) | 9.245 | 9.239 | 1.001 |
| TX(I) | 142.101 | 142.100 | 1.000 |
| MY(I) | 289.991 | 289.900 | 1.000 |
| MZ(I) | 828.470 | 828.400 | 1.000 |
| PX(J) | 12.388 | 12.38 | 1.001 |
| VY(J) | 10.635 | 10.63 | 1.000 |
| VZ(J) | 4.310 | 4.305 | 1.001 |
| TX(J) | 423.720 | 423.7 | 1.000 |
| MY(J) | 261.350 | 261.3 | 1.000 |
| MZ(J) | 541.933 | 541.9 | 1.000 |
Note: PX (I) and PX (J) = Section axial force at node I and J.
VY (I) and VY (J) = Section shear forces along Y direction at node I and J.
VZ (I) and VZ (J) = Section shear forces along Z direction at node I and J.
TX (I) and TX (J) = Section torsional moment at node I and J.
MY (I) and MY (J) = Section bending moments along Y direction at node I and J.
MZ (I) and MZ (J) = Section bending moments along Z direction at node I and J.
The element forces and moments along Y and Z directions are flipped between Mechanical APDL and NRC results.