The natural frequencies of the 2D axisymmetric model without rotation are evaluated and compared with the results of 3D solid model in the following table.
Mode # | 2D Axisymmetric Model (Hz) | 3D Solid Model (Hz) | Error (%) |
1 | 189.72 | 189.00 | 0.38 |
2 | 208.96 | 208.04 | 0.44 |
3 | 639.91 | 639.59 | 0.05 |
4 | 654.54 | 653.25 | 0.20 |
5 | 733.26 | 732.77 | 0.07 |
6 | 807 | 805.64 | 0.17 |
7 | 990.36 | 991.73 | 0.14 |
8 | 1780.5 | 1782.50 | 0.11 |
9 | 1781.1 | 1785.20 | 0.23 |
10 | 2016.9 | 2009.60 | 0.36 |
11 | 2092.6 | 2086.90 | 0.27 |
12 | 3291.7 | 3287.80 | 0.12 |
The natural frequencies of the 2D axisymmetric model in rotation (50,000 rpm) also show good agreement with the 3D solid model results, as shown in the following table.
Mode # | 2D Axisymmetric Model (Hz) | 3D Solid Model (Hz) | Error (%) |
1 | 169.07 | 168.25 | 0.49 |
2 | 232.54 | 231.75 | 0.34 |
3 | 627.97 | 627.19 | 0.12 |
4 | 652.32 | 651.39 | 0.14 |
5 | 752.14 | 751.59 | 0.07 |
6 | 808.51 | 807.25 | 0.16 |
7 | 990.36 | 991.73 | 0.14 |
8 | 1763.00 | 1766.00 | 0.17 |
9 | 1798.90 | 1802.00 | 0.17 |
10 | 1931.30 | 1923.50 | 0.41 |
11 | 2192.50 | 2187.60 | 0.22 |
12 | 3291.70 | 3287.80 | 0.12 |
The /POST1 results of the Campbell diagram analysis are shown in the figure that follows.
With the help of the Campbell diagram analysis, one can identify the forward (FW) and backward (BW) whirls, as well as possible unstable frequencies (though none are present in this example). In the table below, the whirls and natural frequencies of the 2D axisymmetric model at maximum rotational speed (100,000rpm) are compared with the 3D solid model results.
Frequency (Hz) | ||||
---|---|---|---|---|
Whirl | 100000 RPM | |||
Mode # | 2D Axisymmetric Model | 3D Solid Model | 2D Axisymmetric Model | 3D Solid Model |
1 | BW | BW | 144.98 | 144.17 |
2 | FW | FW | 264.59 | 263.92 |
3 | BW | BW | 602.73 | 601.48 |
4 | BW | BW | 649.45 | 648.72 |
5 | FW | FW | 792.94 | 791.73 |
6 | FW | FW | 814.48 | 813.99 |
7 | FW | BW | 990.36 | 991.73 |
8 | BW | BW | 1745.08 | 1748.14 |
9 | FW | FW | 1817.30 | 1820.05 |
The Campbell diagram analysis helps to determine the critical speeds of the rotating structure (PRCAMP ). Critical speeds are compared in the table below. For a synchronous excitation, the critical speeds correspond to the intersection points between the frequency curves and the 1.0 slope line. The critical speeds of the 2D axisymmetric and 3D solid models show strong agreement.
Critical Speeds (RPM) | |||
Mode # | 2D Axisymmetric Model | 3D Solid Model | Error (%) |
1 | 11107.97 | 11064.65 | 0.39 |
2 | 12902.64 | 12847.70 | 0.43 |
3 | 37852.19 | 37812.80 | 0.10 |
4 | 39167.83 | 39107.83 | 0.15 |
5 | 45015.50 | 44982.13 | 0.07 |
6 | 48507.73 | 48431.91 | 0.16 |
7 | 59421.77 | 59503.56 | 0.14 |
8 | none | none | - |
9 | none | none | - |
The results of the unbalance response analysis post processed in /POST26 are shown in the following figure. The logarithmic plots show the variation of the displacement amplitudes of two selected nodes with respect to the frequency of excitation. The first node is located near the rigid disk, and it corresponds to the light blue curve. The second node is located near the bearing location, and it corresponds to the purple curve.
The critical frequencies appear where the amplitudes are largest, and correspond to the critical speeds.
The orbits after a full harmonic analysis can be plotted in POST1 (PLORB) as shown in the figure below. For the solid element and the axisymmetric element model, it is necessary to add massless line elements using BEAM188 on the rotational velocity axis to plot these orbits. The orbits of the 2D axisymmetric model at a frequency of 1666.67 Hz are shown in the figure below. The rotor line is in dark blue, while the orbits are in light blue.
The following input fragment shows the steps to produce an orbits plot at a given frequency:
/POST1 esel,r,ename,, 188 ! Select BEAM188 elements to produce orbits set,1, 200 ! Visualize orbits at frequency 1666.67 Hz /view,,1,1,1 plorb ! Displays the orbital motion of a rotating structure
The memory and CPU usage of the 2D model is shown in the following table.
2D Axisymmetric Model | |||||
---|---|---|---|---|---|
Elements # | Nodes # | No. of Equations | Memory required for in-core (MB) | CPU Time (Sec) | |
Campbell Diagram Analysis | 2208 | 6751 | 20225 | 53.283 | 17.30 |
Unbalance Response Analysis | 118.325 | 1347.53 |
The memory and CPU usage of the 3D model is shown in the following table.
3D Solid Model | |||||
---|---|---|---|---|---|
Elements # | Nodes # | No. of Equations | Memory required for in-core (MB) | CPU Time (Sec) | |
Campbell Diagram Analysis | 9239 | 15123 | 45341 | 186.141 | 45.56 |
Unbalance Response Analysis | 605.464 | 4645.95 |
The CPU times for the unbalance response analysis are represented in the following bar graph.