4.1. Initial Finite Element Modeling and Structural Analysis
A certain dual-rotor turbofan engine with a high bypass ratio was used to establish the engine dynamics model. The engine weighed 2.80 T, had a thrown ratio of 3.80, and maximum thrust of 86.7 KN. The low-pressure rotor system consisted of 1 stage fan, 3 stage compressors, and 4 stage low-pressure turbines. The high-pressure rotor system consisted of 9 stage high-pressure compressors and 1 stage high-pressure turbine. The stator system is divided into the outer bypass casing and the inner bypass casing. The geometric model of the whole machine is established, and the semi-sectional view of the model is shown in
Figure 4a, while the FE model of the whole engine is shown in
Figure 4b.
The double-rotor turbofan engine used in this study had five supports and two mountings. Its overall supports and mounting plan are shown in
Figure 5. The main mounting section is located in the intermediate casing, and the auxiliary mounting section is located at the upper end of the turbine rear casing. The support stiffness value is taken as listed in
Table 1, and the values for the main and auxiliary mounting are set as 3 × 10
9 N/m.
When the twin-rotor aero engine is operating, there is a certain speed relationship between the high-pressure rotor and the low-pressure rotor; usually, the method of controlling the high-pressure rotor speed and the low-pressure rotor speed adaptation is adopted to achieve speed coordination. In the whole speed range, the speed ratio of high- and low-pressure rotors is not constant, and the speed of high- and low-pressure rotors conforms to the law of common working line of high- and low-pressure rotors. The common working line of high- and low-pressure rotors of a certain type of high-bypass ratio twin-rotor aero-engine established in this paper is shown in
Figure 6. In the entire speed range, there are three operating speeds, namely ground slow speed (low pressure rotor speed 1151 rpm, high-pressure rotor speed 9879 rpm), cruising speed (low-pressure rotor speed 4358 rpm, high-pressure rotor speed 14,022 rpm) and takeoff speed, that is, the maximum speed (low-pressure rotor speed 4779 rpm, high-pressure rotor speed 14,511 rpm).
Referring to the common working line of the high-pressure rotor (HPR) and low-pressure rotor (LPR) of a certain type of aero-engine, seven groups of high- and low-pressure speeds were selected on the common working line (
Table 2), and the critical speed of the system was analyzed by the Campbell diagram method.
The Campbell plot is plotted with the LPR speed as the abscissa and the precession frequency as the ordinate, as shown in
Figure 7. It can be seen from the figure that the critical speed of low-pressure rotor excitation in the operating speed range have two orders, 2755 rpm and 2881 rpm, respectively. Obviously, when the low-pressure rotor speed is 0 rpm, the forward and reverse precession speeds of each stage are not the same, and as the order increases, the difference between the forward and reverse precession speeds increases. This phenomenon is mainly caused by the asymmetry of rotor support stiffness in the horizontal and vertical directions.
The Campbell diagram is drawn with the HPR speed as the abscissa and the precession frequency as the ordinate, as shown in
Figure 8; the intersection point of the iso-speed line and the high-pressure rotor forward precession line is the critical speed of the high-pressure rotor excitation. It can be seen from the figure that the critical speed of the high-pressure rotor excitation within the operating speed range is only 1 order, which is 12,669 rpm. Due to the dense mode and clear performance, only the forward precession lines of each order are drawn for the high-pressure rotor.
The first 3 critical speed risk coefficients, rotor strain energy risk coefficients and cross-section rotor–stator rubbing risk coefficients of the engine are calculated, respectively. The structural vibration of the engine is evaluated and analyzed using three defined coefficients.
According to the calculated critical speed of the engine, the critical speed margin (that is,
), and then the critical speed risk coefficient of each order is obtained. Among them, the subscript L means low-pressure rotor and H means high-pressure rotor. For example,
CiL is the low-pressure rotor critical speed risk coefficient of the
i-th order, and
CiH is the high-pressure rotor critical speed risk coefficient of the
i-th order.
W1L represents the first operating speed of the low-pressure rotor, which is 1151 rpm,
W2L represents the second operating speed of the low-pressure rotor, which is 4358 rpm,
W3L represents the third operating speed of the low-pressure rotor, which is 4794 rpm, and
W1H,
W2H and
W3H are the corresponding the first, the second and the third high pressure rotor operating speeds, which are 9879 rpm, 14,022 rpm and 14,511 rpm respectively. The calculation process and results of the low-pressure critical speed risk coefficient and the high pressure critical speed risk coefficient of each order are shown in
Table 4 and
Table 5.
It can be shown in
Table 4 and
Table 5, the critical speed risk coefficients of the first 3 order low-pressure rotor are 0.632, 0.661 and 0.696, respectively, and the critical rotor risk coefficients of the first 3 order high-pressure rotor are 0.883, 0.892 and 0.904, respectively. Among them, the critical speed of the 3rd order high-pressure rotor is the most dangerous, and the critical speed risk coefficient reaches 0.9. Its safety should be further verified.
The strain energy analysis is carried out under the mode corresponding to each critical speed. The percentage of strain energy of each support, high- and low-pressure rotor and mountings corresponding to the three order critical speed are shown in
Table 6, where
BiL represents the percentage of low pressure rotor strain energy of the modal shape at the critical speed of the
i-th order, which is the rotor strain energy risk coefficient of the low pressure rotor of
i-th order. Similarly,
BiH represents the high-pressure rotor strain energy risk coefficient of modal shape at the
i-th order critical speed.
- 1)
The modal shape at the 1st critical speed (N
1 = 2755 rpm, N
2 = 12,384 rpm) of low-pressure excitation is characterized by the 1st order bend of the low-pressure rotor shaft, low-pressure turbine pitch, and high-pressure rotor pitch, as shown in
Figure 9. In this modal shape, 73.40% of the strain energy is concentrated on the low-pressure rotor, and this modal shape will be very sensitive to the imbalance of the low-pressure rotor, and its vibration amplitude needs to be strictly controlled. Secondly, in this modal shape, the strain energy ratio of the support 5 reaches 13.99%, which is mainly caused by the low-pressure turbine pitch of the order modal shape, and the fatigue life of the bearing should also be noted.
- 2)
The modal shape at the 2nd critical speed (N
1 = 2881 rpm, N
2 = 12,512 rpm) excitation by low-pressure rotor is mainly manifested as fan pitch, as shown in
Figure 10. 51.44% of the strain energy is concentrated in the low-pressure rotor, mainly distributed on the low-pressure fan rotor, in this order modal shape, the fan will be very sensitive to imbalance, and its dynamic balance and vibration amplitude should be strictly controlled. Support 1 withstands 24.08%, and support 2 withstands 12.13%. In this modal shape, the structural strength of the connection between the fan rotor and the low-pressure turbine rotor should attract sufficient attention. Notably, the stator withstands 10.13% of the strain energy; mainly due to support 1 and support 2, located on the stator structure, some energy transfer occurs, which alleviates the energy concentration of the low-pressure fan rotor and support 1 and support 2 to a certain extent, which is conducive to the stability of the structural vibration of the whole engine.
- 3)
The modal shape at the 3rd critical speed of high-pressure excitation (N
1 = 3034 rpm, N
2 = 12,669 rpm) is manifested as low-pressure rotor second-order bending, low-pressure turbine pitching, and high-pressure rotor pitching, as shown in
Figure 11. In total, 24.07% of the strain energy is concentrated in the low-pressure rotor and 21.89% in the high-pressure rotor. In addition, 12.27% and 20.06% are distributed in support 3 and support 4, respectively. In this order modal shape, attention needs to be paid to the loads of support 3 and support 4. In this modal shape, the stator structure bears 14.69% of the strain energy, which is distributed in the inner and outer casing, which alleviates the load of the rotor and support to a certain extent.
According to the first three order modal shape, the five sections which have the greatest rubbing risk are selected as the key sections. Cross sections 1–5 are the fan section, the 1st stage low-pressure compressor section, the 1st stage high-pressure compressor section, high-pressure turbine section and the 4th stage low-pressure turbine section; the location of each section is shown in
Figure 5. The modal displacement of the rotor tip of each key section and the corresponding stator casing position are extracted, and the modal displacement of the top node of the low-pressure rotor rectifier cone is set as 1, and each modal shape is normalized to calculate the minimum relative displacement of the rotor blade and stator casing at each section. According to Formula (3), the cross-section rotor–stator rubbing risk coefficients are calculated, as shown in
Table 7.
- 1)
Under the 1st modal shape, the 1st order bend of the low-pressure rotor shaft, low-pressure turbine pitch, and high-pressure rotor pitch. The cross-section rotor–stator rubbing risk coefficient of
Section 4 and
Section 5 is large, indicating that under the modal shape, the high-pressure turbine section and the 4th stage low-pressure turbine section have a higher degree of rubbing, which is in line with the characteristics of the mode shape.
- 2)
Under the fan pitch modal corresponding to the 2nd order critical speed, the cross-section rotor–stator rubbing risk coefficients of
Section 1 and
Section 2 are the highest. It shows that under this modal shape, the fan section and the 1st stage low-pressure compressor section have a high degree of rubbing risk, which is in line with the characteristics of this modal shape.
- 3)
Under the 3rd modal shape, low-pressure rotor second-order bend, low-pressure turbine pitches, and high-pressure rotor pitches,
Section 3 and
Section 4 have the highest cross-section rotor–stator rubbing risk coefficient. It shows that under the modal shape, the 1st stage high-pressure compressor section and the high-pressure turbine section have a high degree of rubbing, which is in line with the characteristics of this modal shape.
4.2. Determining the Optimization Model
A multi-objective optimization method is used to solve the aeroengine support stiffness optimization design problem. The first three orders’ critical speed risk coefficient, rotor strain energy risk coefficient, and rotor-stator rubbing risk coefficient for five cross sections (cross section locations shown in
Figure 5) are determined as the optimization objectives. Since the engine has two rotors and the three order critical speeds in the working speed range have to be analyzed, there are 27 index parameters; thus, the optimization problem becomes complicated. To solve the aeroengine support stiffness optimization problem efficiently and accurately, the calculated indicators are selected and processed. As shown in
Figure 12, the most dangerous index in each group of stiffness is selected as the optimization target. The selected optimization targets are the most dangerous critical speed risk factor
f1, the highest rotor strain energy risk factor
f2 and the most dangerous rotor–stator rubbing risk coefficient
f3.
Aeroengine bearing 4 is an intermediate support. It has high stiffness and is subjected to complex working conditions including various factors, and the scope of adjustment is limited in the late design stage. Hence, it was not considered for optimization. The influence law of the stiffness of other supports on each optimization objective was analyzed. The influence laws of each support stiffness on
f1,
f2 and
f3 are shown in
Figure 13,
Figure 14 and
Figure 15, respectively.
As shown in
Figure 7, with the increase in the stiffness of support 1,
f1 first remains unchanged and then increases. The influence of support 2 on
f1 is small. Meanwhile,
f1 increases with an increase in the stiffness of support 3. With the increase in the stiffness of support 5,
f1 increases linearly.
As shown in
Figure 13, with the increase in the stiffness of support 1 and support 2,
f2 decreases. Meanwhile,
f1 increases with an increase in the stiffness of support 3 and support 5.
As shown in
Figure 14, with the increase in the stiffness of support 1 and support 2,
f3 decreases. The influence of support 3 on
f3 is small. With the increase in the stiffness of support 5,
f3 first decreases and then remains almost unchanged.
Based on the analysis of the influence of the stiffness of each support on each optimization target, the stiffness values of supports 1, 2, 3, and 5 were determined to be the parameters to be optimized. The stiffness values of supports 1, 2, and 5 vary from 2 × 10
7 N/m to 1 × 10
8 N/m, and the stiffness values of support 3 range from 6 × 10
7 N/m to 2 × 10
8 N/m. The optimization model is finally determined as follows:
4.5. Computing Agent Model Obtaining Based on SVM
Support vector regression analysis is performed on the obtained sample data, and the mapping relationships between
f1,
f2,
f3 and
k1,
k2,
k3 and
k5 are obtained. The mapping relationship reflects the functional relationship between the input (
k1,
k2,
k3,
k5) and the output (
f1,
f2,
f3) of the FE model, which can be shown as follows:
The NAGA-II parameters are set to cross probability 0.9, probability of variation 0.1, population size 1000, maximum evolutionary algebra 1000, tournament size 2, cross distribution index 20, and variation distribution index 20.
To verify the accuracy of the training model, using the cross-validation method, 350 sets of samples are selected as training samples and 50 groups of samples are selected as test samples, and three models are verified respectively; the results are shown in
Figure 18. It can be used to predict the objective function values of each stiffness combination of a certain type of engine in place of the FE calculation.
As shown in the figure, the predicted objective function values are in good agreement with the actual objective function values, and the model is considered to have good generalization capability.
It can be used to predict the objective function values of each stiffness combination of a certain type of engine in place of FE calculation.