14.7. Results and Discussion

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.

Figure 14.15: Campbell Diagram

2D Axisymmetric Model
Campbell Diagram
3D Solid Model
Campbell Diagram

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)
Whirl100000 RPM
Mode #2D Axisymmetric Model3D Solid Model2D Axisymmetric Model3D Solid Model
1BWBW144.98144.17
2FWFW264.59263.92
3BWBW602.73601.48
4BWBW649.45648.72
5FWFW792.94791.73
6FWFW814.48813.99
7FWBW990.36991.73
8BWBW1745.081748.14
9FWFW1817.301820.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.

Figure 14.16: Unbalance Response

2D Axisymmetric Model
Unbalance Response
3D Solid Model
Unbalance Response

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

Figure 14.17: Orbits Plot of 2D Axisymmetric Model

Orbits Plot of 2D Axisymmetric Model

14.7.1. Performance Benefits of the 2D Axisymmetric Model

The memory and CPU usage of the 2D model is shown in the following table.

2D Axisymmetric Model
Elements #Nodes #No. of EquationsMemory required for in-core (MB)CPU Time (Sec)
Campbell Diagram Analysis220867512022553.28317.30
Unbalance Response Analysis118.3251347.53

The memory and CPU usage of the 3D model is shown in the following table.

3D Solid Model
Elements #Nodes #No. of EquationsMemory required for in-core (MB)CPU Time (Sec)
Campbell Diagram Analysis92391512345341186.14145.56
Unbalance Response Analysis605.4644645.95

The CPU times for the unbalance response analysis are represented in the following bar graph.

Figure 14.18: CPU Time of 2D Axisymmetric and 3D Solid Models (Unbalance Response Analysis)

CPU Time of 2D Axisymmetric and 3D Solid Models (Unbalance Response Analysis)