Research Article  Open Access
Zhijun Li, Wenjun Yang, Huiqun Yuan, "Vibration Analysis of Aeroengine Blisk Structure Based on a Prestressed CMS SuperElement Method", Shock and Vibration, vol. 2016, Article ID 1021402, 10 pages, 2016. https://doi.org/10.1155/2016/1021402
Vibration Analysis of Aeroengine Blisk Structure Based on a Prestressed CMS SuperElement Method
Abstract
For vibration analysis of aeroengine blisk structure, a prestressed component modal synthesis (CMS) superelement method is put forward with the fixed interface prestressing and free interface superelement approach. Based on this method, natural vibration characteristics of blisk structure are calculated at different modal truncation numbers. Comparing with the accurate result of global method, the selection principle of modal truncation number is obtained which affects the accuracy of prestressed CMS superelement method. Vibration response of twostage blisk structure is calculated by this method, and the effects of different blade aspect ratios have been discussed on vibration characteristics. The results show that prestressed CMS superelement method is in the high accuracy and efficiency on blisk vibration analysis. Resonant frequencies in vibration response are nearly the same between the firststage blisk and the secondstage blisk, and they are both approximately located in the range 588 Hz–599 Hz. The maximum displacement and dynamic stress are at blade tip and root of the firststage blisk, respectively. Blade aspect ratio is a key factor of blisk vibration; the effects of blade aspect ratio on natural frequencies are different in the conditions of fixed width and fixed length. This research provides the theoretical basis for dynamic design of aeroengine compressor rotor system.
1. Introduction
With the increasing development of aeronautical manufacturing technology, a whole system of bladed disk (blisk) has been widely used in newly developed aeroengines. Blade and disk are integrated to the whole system of bladed disk by the advanced technology. Traditional joint of tenon and mortise is removed out between blade and disk; it can greatly simplify the structure to achieve the lightweighting design. Multiple crack failure is avoided in mortise bottom, and the reliability of rotor system is improved significantly.
Many scholars have carried out extensive research on the blisk structure. Ferria et al. [1] developed a blisk numerical model on the flutter stability of subsonic turbine and pointed out the conditions which affected the stability of turbine system. Di Maio and Ewins [2] implemented a practical method to measure vibration response of simplified blisk structure by Scanning Laser Doppler Vibrometer systems. Ji et al. [3] put forward a multidisciplinary optimization method for the design of compressor structure, and this method was verified that it could maintain the harmony and consistency of blisk structure well. Bhaumik et al. [4] conducted a theoretical research on the failure mechanism of turbine blisk; the criteria of avoiding failure were given according to the performance of turbine material. S. Lu and F.J. Lu [5] presented a weightlighting optimized design of blisk structure with a guarantee of structural safety. In the field of crack growth, Xu et al. [6] developed some measures to improve the reliability of axial compressor blisk system.
Due to the limitation of experimental conditions and lack of appropriate experimental methods, major approach is the finite element analysis (FEA) in the analysis of structural dynamics. However, for the large complex structures such as aeroengine blade disk system, the number of degrees of freedom (DOFs) can be up to one million in discrete model. Corresponding kinematic equations cannot be solved efficiently, and it is difficult to develop dynamic analysis with the finite element method. Although it can be analyzed with FEA, a lot of time would be eventually consumed. Thus, calculation efficiency cannot be guaranteed. From the review of related studies, it can be found that computational method of blisk dynamics remains for further research. So the method of dynamic substructure becomes an appropriate solution method on the basis of modal reduction technology.
Theoretical method and practical application of modal synthesis technology have been researched by related scholars [7–10]. Hurty firstly established the concepts of modal coordinates and modal synthesis method [11], which laid the foundation of fixed interface modal synthesis method. Then, Bampton and Craig Jr. [12] proposed an improved method to make the fixed interface method more simple and practical. This method eliminated the boundary of rigid mode and constraint mode and no longer distinguished between both of them. Hou [13] and Goldman [14] explored the modal synthesis method with free interface. While the effect of higherorder substructure modes is ignored, the accuracy of this method is challenged. MacNeal [15] and Rubin [16] introduced the residual stiffness to consider the effect of higherorder truncated modes; global precision was improved. For further improving Rubin’s method, Wang and Du [17] put forward the dynamic substructure technique with double coordination. In this method, residual stiffness was regarded as the Ritz base; modified free interface method was put into the orbit of Ritz analysis. Comprehensive precision and efficiency were improved greatly. With the reduction technology, Leung [18] concentrated internal coordinates of each substructure to the coordinates of each substructure interface. Motion equation of the system was established according to the conditions of displacement coordination and force equilibrium. Based on Leung’s method, Yun et al. [19] developed a superelement modal synthesis method by the frequency conversion of dynamic modes. Although there is a certain progress in the technique of dynamic substructure, its application still needs to be further researched on dynamic analysis.
For the analysis of dynamic characteristics in an aeroengine blisk structure, in our research a prestressed component modal synthesis (CMS) method is proposed with the fixed interface prestressing and free interface superelement approach. Based on this method, natural vibration characteristics of blisk structure are calculated at different modal truncation numbers. Comparing with the accurate result of global method, the selection principle of modal truncation number is obtained which affects the accuracy of prestressed CMS superelement method. Vibration response of twostage blisk structure is calculated by this method, and the effects of different blade aspect ratios have been discussed on vibration characteristics.
2. Dynamics Modeling of Aeroengine Blisk Structure
2.1. The Method of Prestressed CMS SuperElement Method
Blisk structure of aeroengine system is very complex. It is a challenge for numerical simulation, as the number of finite elements is much enormous. In this research, the prestressed CMS superelement method is proposed to analyze dynamic characteristics of blisk structure. By this method, it can solve the structure with large numbers of elements. Moreover, the prestressed CMS superelement method has higher precision and efficiency compared with traditional method.
The prestressed CMS superelement method is based on modal synthesis technology, and it is a method which utilizes matrix reduction technology to reduce model order. For the th substructure finite element model, the general undamped free vibration equation can be expressed aswhere is the mass matrix, is the stiffness matrix, is the displacement vector, and is the number of substructures.
For components such as blades, they are very thin on the direction of one or two degrees of freedom. At the action of centrifugal force, stress state may affect structural natural frequency and dynamic response. So in the analysis of rotor dynamic characteristics, the effect of centrifugal rigidification should be considered. In our research, the linear stress analysis is developed under the static state; centrifugal load is transformed into structure prestressed matrix . Then, dynamic equation related to centrifugal rigidification can be expressed aswhere is the stiffness matrix which considers the matrix of prestressed effect. Displacement vector and coefficient matrix in (2) are divided into the master DOF on the boundary and the slaver DOF beyond the boundary. Here subscripting m denotes the master DOF and subscripting s denotes the slaver DOF. The transformed form is
In the equation, is the interfacial force. Coordinates of substructure are transformed by the following formula:
Here is the coordinate of substructure mode, is the transformed coordinate matrix in the fixed interface [12], is substructure eigenvector in the condition of fixed boundary nodes, and is the unit matrix. According to formula (4), dynamic equation in formula (2) can be expressed with modal coordinates:
Here . Based on the DOF reduction method [11], the following expression can be obtained from formula (4):
In formula (6), is the reduced modal coordinate which is transformed from DOF generalized coordinate . is the new transformed coordinate matrix: . Here is the loworder modal set with higherorder truncation. Thus, modal coordinate vector in (5) is transformed as .
To consider the rigid connection between substructures, the following coordinate transformation [12] is utilized to transform nonindependent coordinate into the generalized coordinate of global structure.
Thus, the free vibration equation is established as follows:
By the reduction of DOF, and are much less than mass matrix and stiffness matrix in the generalized coordinates of original system. Thus, natural frequency and modal shape can be obtained from (8).
2.2. Physical Model and Analysis Process
Blisk structure is a new manufacturing technology of aeroengine. The joint between tenon and mortise is removed out, so the weight of tenonmortise connection and supporting structure is reduced greatly. In addition, bolts, nuts, locking plate, and other connectors are no longer needed. Compared with traditional blade disk, the hub becomes thinner in blisk structure, and bore diameter becomes larger. Figure 1 shows the structure of first twostage blisk in aeroengine compressor. Number of blades in the first stage is 38, and number of blades in the second stage is 53.
(a)
(b)
In the finite element model of blisk structure, parts of blade and disk are established with the element Solid 185. As the blades on disk are leaned, the joint zone of blade and disk is meshed with the element Solid 187. Here total number of elements is 400124 and total number of nodes is 561878. Blade material is titanium alloy TA11, the density is 4400 kg·m^{−3}, elasticity modulus is 114 GPa, and Poisson’s ratio is 0.3. Disk material is titanium alloy TC17, the density is 8200 kg·m^{−3}, elasticity modulus is 166 GPa, and Poisson’s ratio is 0.3.
For applying the prestressed CMS superelement method, substructure models of twostage blisk are established in this research. The firststage blisk and the secondstage blisk are packed as a substructure, respectively, as shown in Figure 2. In the firststage substructure, it contains 194825 elements and 275406 nodes. In the secondstage substructure, it contains 205299 elements and 287432 nodes.
(a)
(b)
After establishing FEA model of substructure, the DOFs in hub tube are constrained for further analysis. Considering the effect of centrifugal force, operating rotational speed is applied. Prestressed analysis with fixed interface is performed in each substructure. Then, prestressed option is set as open; constrained master DOFs are released. Generation process of substructure mode synthesis is performed in the free interface, and superelements are created. After these, application process is performed. The superelements of substructures are connected to take the vibration analysis of whole model. Finally, dynamic response of superelement master DOFs is extended into each inner DOF of superelement orderly. Complete solution of dynamic response is obtained, and expansion process of the prestressed CMS method is accomplished.
Analysis process of prestressed CMS superelement method is shown in Figure 3.
3. Vibration Characteristics of Blisk Structure
3.1. Accuracy Verification of the Prestressed CMS Method
The modal truncation number is a key factor of calculation accuracy in prestressed CMS method. For obtaining more suitable modal truncation number in this research, natural vibration characteristics of blisk structure are calculated at different modal truncation numbers. Substructure models of twostage blisk have been established, and operating speed of blisk structure is 11383 rpm. The DOFs of nodes in hub section are constrained as the boundary conditions. By the prestressed CMS method, natural frequencies of twostage blisk are solved, respectively, at different modal truncation numbers. Here the result of global method is assumed as the accurate value. and are, respectively, the modal truncation number in the firststage blisk and the secondstage blisk.
In Figure 4, it is shown that calculation result of the prestressed CMS method at and is far away from the frequency of global method. With the increase of and , calculation results of the prestressed CMS method are much closer to the accurate solution. When and , the result of prestressed CMS method is basically consistent with the result of global method. Continuing to increase the modal truncation number such as and , calculation accuracy of prestressed CMS method is much higher. Therefore, it illustrates that the requirement of calculation precision can be satisfied.
For further examining the calculation efficiency of prestressed CMS method, computing time of prestressed CMS method is compared with that of global method at different modal truncation numbers, as is shown in Table 1.

As shown in Table 1, computing time of prestressed CMS method increases with the growth of modal truncation number. When the modal truncation number rises to a certain value, computing time of prestressed CMS method is nearly equal to computing time of global method. It has been known that calculation accuracy of prestressed CMS method is improved significantly with the increasing of modal truncation number. While computing time is much longer, calculation efficiency of prestressed CMS method cannot be guaranteed. So the calculation accuracy and efficiency should be both considered to determine the modal truncation number.
From the above analysis, it can be known that prestressed CMS method has high enough precision at and . Compared with the global method, calculation efficiency is improved by 36%. For dynamic analysis of twostage blisk structure, hence and are selected as the modal truncation number. The results are compared with global method, as can be seen from Table 2.

Table 2 exhibits that 38order vibration frequencies of the firststage blisk are in the range of 434 Hz–443 Hz and 53order vibration frequencies of the secondstage blisk are in the range of 546 Hz–600 Hz. It can be found that natural frequencies perform the bending shape of blade in loworder modes. The number of vibration modes is equal to the number of blades. The results of global method and prestressed CMS method have very small difference; it illustrates that prestressed CMS method is credible. In order to ensure the accuracy of prestressed CMS method, truncated frequency of the substructure is required to be greater than the corresponding frequency of solved system. For dynamic analysis of blisk structure, the modal truncation number is required to be greater than the number of blades. This selection principle of modal truncation number has been verified according to Figure 4 and Table 2.
3.2. Vibration Response of TwoStage Blisk Structure
In the operation of blisk structure, main source of vibration is caused by the uneven air flow on blade pressure and suction surfaces. Airflow exciting force can be estimated on the basis of aerodynamic calculation, and the velocity and pressure on rotor blades can be determined by the experiment. Airflow exciting force at average radius of air passage can be expressed as the form of Fourier series; it is a series of harmonic superposition. For the system of blisk structure, practical exciting form is very complex; it has a great relationship with the working conditions. Assume that airflow exciting force on blade surface is simplified as the single point form. Aerodynamic force is applied at each blade tip of leading edge; the load in space is required to meet the form of traveling wave.
For the system of blisk structure, motion equation of forced vibration is as follows:
Here is the vector of displacement and is the vector of exciting force. The parameters , , and are mass matrix, viscous damping matrix, and stiffness matrix, respectively. In the forced response of blisk structure, exciting force is usually expressed as
In formula (10), is the amplitude of exciting force in th blade, is the frequency of exciting force, and is the number of blades. is the phase of exciting force in th blade, and it is defined as
Here is the order of exciting force.
3.2.1. Response Analysis of Vibration Displacement
Figure 5 shows the relationship between maximum displacement and exciting frequency, and the responses of vibration displacement in twostage blisk are compared.
(a)
(b)
In Figure 5, it is exhibited that resonant frequency range of twostage blisk is at 588 Hz–599 Hz. The frequency is around 594 Hz at the maximum response. It is clear to note that the maximum displacement of the firststage blisk is much higher than that of the second stage.
Figure 6 shows the displacement contour maps of twostage blisk at resonant frequency of 594 Hz.
(a)
(b)
In Figure 6, it can be seen that when the frequency of exciting force is at 594 Hz, the vibration of twostage blisk is mainly focused on the firststage blisk. Additionally, the maximum displacement appears at the position of blade tip, and the displacement vibration shows a nodal diameter mode. The amplitude of displacement vibration is relatively small in each region of the secondstage blisk.
3.2.2. Response Analysis of Dynamic Stress
The load acting on the element has a remarkable change with time, or each node on the component has significant acceleration under the load. The stress is generated by dynamic load in the component, and it is called dynamic stress. Dynamic stress analysis is the basis of solving the problem of component dynamic failure. According to the analysis of harmonic response, dynamic stress of twostage blisk is obtained under aerodynamic excitation force.
Figure 7 shows the relationship between maximum dynamic stress and exciting frequency; the responses of vibration stress in twostage blisk are compared.
(a)
(b)
From Figure 7, it is found that the response of dynamic stress is basically consistent with the response of vibration displacement. The range of resonant frequency is at 588 Hz–599 Hz, and the response peak occurs at the frequency of 594 Hz. Moreover, the maximum dynamic stress of the firststage blisk is much higher than that of the secondstage blisk.
Figure 8 shows the stress contour maps of twostage blisk at resonant frequency of 594 Hz.
(a)
(b)
In Figure 8, it can be shown that dynamic stress occurs mainly at blade root position of the firststage blisk, and the stress vibration shows a nodal diameter mode. Besides, the stress amplitude of the secondstage blisk is relatively small in each region.
In order to avoid resonant response of blisk structure, the frequency of external load should be far away from the resonant frequency of 594 Hz. In addition, it can be found that the maximum vibration displacement is 6.21 mm and the maximum dynamic stress is 749 MPa in the limit working condition of 594 Hz, while vibration amplitude and material strength of blisk structure are both in the safe working condition.
4. Discussion of Blisk Vibration at Different Aspect Ratios
Blade aspect ratio is the ratio of blade length and blade width; it stands for blade relative length or blade relative width. Blade aspect ratio is one of key factors on blisk vibration characteristics. In order to explore the effects of blade aspect ratios, natural frequencies of blisk structure are discussed at different aspect ratios.
4.1. Design and Modeling of Blisk Structure
For discussing the effects of aspect ratio on blisk vibration, blisk models are established at different aspect ratios. Relevant parameters including blade inclination and wheel size are kept constant in the process of modeling; blade aspect ratio is only adjusted accordingly. Thus, analysis results can be deemed to be credible.
According to practical experiences, blade aspect ratio is set as 1.50, 1.75, 2.00, 2.25, and 2.50, respectively, as shown in Tables 3 and 4. At the conditions of fixed width and fixed length, vibration characteristics of blisk structure are discussed.


After establishing blisk models at different aspect ratios, substructure models need to be divided. Substructure can be natural component of global structure, and it also can be a certain part of manual separation. As shown in Figure 9, blisk structure is divided into sectors. Each sector of blisk is regarded as a substructure.
(a)
(b)
4.2. Blisk Vibration Characteristics at Different Aspect Ratios
Considering the effect of centrifugal force, blisk models at different aspect ratios are analyzed with the prestressed CMS superelement method. At the conditions of fixed width and fixed length, natural vibration frequencies of blisk structure are obtained. As blisk system has many blades, frequencies of the same vibration shapes are similar. Therefore, the first 120order natural frequencies of blisk structure are solved. With the increase of blade aspect ratio, variation trends of natural frequencies and description of vibration shapes are illustrated as shown in Table 5.
 
“—” stands for no obvious change, “↗” stands for rise, and “↘” stands for decline. 
For observing the relationship between natural frequencies and aspect ratios, one typical order of the similar frequencies is selected as the representative. The frequencies of 1st order, 39th order, 77th order, 84th order, and 120th order have been extracted as shown in Tables 6 and 7.


According to the data in Tables 6 and 7, effect curves of aspect ratios on natural frequencies are drawn at the conditions of fixed width and fixed length, as shown in Figure 10.
(a)
(b)
From the analysis of effect curves in Figure 10, it can be found that each order frequency of blisk structure declines with the increasing of blade aspect ratio in the condition of fixed width. And the curve of high order frequency is much steep; it illustrates that the effect of blade aspect ratio on high order frequency is more obvious than the effect on loworder frequency. In the condition of fixed length, blade distorted frequencies between 40th order and 77th order show a certain rise with the increasing of blade aspect ratio. However, other order frequencies have no clear change.
5. Conclusions
In this research a prestressed CMS superelement method is put forward for the vibration analysis of aeroengine blisk structure. Based on this method, dynamic characteristics of blisk structure are calculated at different modal truncation numbers. And the effects of different blade aspect ratios have been discussed on blisk vibration characteristics. Through the above analysis, we can draw a conclusion.(1)Compared with the result of global method, the accuracy of prestressed CMS method can meet the requirement of blisk dynamic analysis. For the selection principle of modal truncation number, natural frequency of the substructure is required to be greater than the corresponding frequency of solved system.(2)Resonant frequencies of the firststage blisk and the secondstage blisk are basically consistent; they are mainly at 588 Hz–599 Hz. The maximum displacement and maximum dynamic stress appear at blade tip and blade root of the firststage blisk, respectively, and show the vibration mode of nodal diameter.(3)Effects of aspect ratio on blisk vibration are different at the conditions of fixed width and fixed length. Natural frequencies of blisk structure decline with the increasing of blade aspect ratio in the condition of fixed width, and the effect of blade aspect ratio is more obvious on high order frequency, while blade distorted frequencies show a certain rise in the condition of fixed length.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
The work is supported by the National Science Foundation of China (Grant no. 51275081) and the National Science Foundation of China (Grant no. 51335003).
References
 H. Ferria, P. Ferrand, F. Pacull, and S. Aubert, “Numerical investigation of flutter stability in subsonic space turbine blisk with emphasis on cuton/cutoff modes and interblade phase angles,” Journal of Thermal Science, vol. 21, no. 6, pp. 492–499, 2012. View at: Publisher Site  Google Scholar
 D. Di Maio and D. J. Ewins, “Experimental measurements of outofplane vibrations of a simple blisk design using blade tip timing and scanning LDV measurement methods,” Mechanical Systems and Signal Processing, vol. 28, pp. 517–527, 2012. View at: Publisher Site  Google Scholar
 J. Ji, D. Zhang, S. Li, and B. Chen, “A master model approach to multidisciplinary design and simulation for compressor blisk,” Advanced Materials Research, vol. 308–310, pp. 209–212, 2011. View at: Publisher Site  Google Scholar
 S. K. Bhaumik, T. A. Bhaskaran, R. Rangaraju, M. A. Venkataswamy, M. A. Parameswara, and R. V. Krishnan, “Failure of turbine rotor blisk of an aircraft engine,” Engineering Failure Analysis, vol. 9, no. 3, pp. 287–301, 2002. View at: Publisher Site  Google Scholar
 S. Lu and F.J. Lu, “Structure optimization design for blisk based on ANSYS,” Journal of Aerospace Power, vol. 27, no. 6, pp. 1218–1224, 2012. View at: Google Scholar
 F. Xu, C.S. Li, Q.S. Luo, B.H. Zhou, and C.B. Wu, “Study on the blade crack of an axial flow compressor blisk,” Gas Turbine Experiment and Research, vol. 25, no. 4, pp. 25–29, 2012. View at: Google Scholar
 Y. G. Tsuei and E. K. L. Yee, “Direct component modal synthesis technique for system dynamic analysis,” Journal of Biological Chemistry, vol. 171, no. 2, pp. 639–640, 1947. View at: Google Scholar
 T. Mochizuki and I. Hagiwara, “A comparison between modal differential substructure method and conventional component modal synthesis methods,” Journal of System Design & Dynamics, vol. 5, no. 2, pp. 320–331, 2011. View at: Publisher Site  Google Scholar
 J.B. Qiu, Z.G. Ying, and L. H. Yam, “New modal synthesis technique using mixed modes,” AIAA Journal, vol. 35, no. 12, pp. 1869–1875, 1997. View at: Publisher Site  Google Scholar
 P. Ren, L. Wei, and G. Gong, “Analysis for modal characteristics of crane boom based on modal synthesis method,” Hoisting & Conveying Machinery, vol. 7, pp. 59–64, 2014. View at: Google Scholar
 W. C. Hurty, “Vibrations of structural systems by component mode synthesis,” Journal of the Engineering Mechanics, vol. 86, pp. 51–70, 2015. View at: Google Scholar
 M. C. C. Bampton and R. R. Craig Jr., “Coupling of substructures for dynamic analyses,” AIAA Journal, vol. 6, no. 7, pp. 1313–1319, 1968. View at: Publisher Site  Google Scholar
 S. N. Hou, “Review of modal synthesis techniques and a new approach,” Shock and Vibration Bulletin, vol. 40, no. 4, pp. 25–39, 1969. View at: Google Scholar
 R. L. Goldman, “Vibration analysis by dynamic partitioning,” AIAA Journal, vol. 7, no. 6, pp. 1152–1154, 1969. View at: Publisher Site  Google Scholar
 R. H. MacNeal, “A hybrid method of component mode synthesis,” Computers & Structures, vol. 1, no. 4, pp. 581–601, 1971. View at: Publisher Site  Google Scholar
 S. Rubin, “Improved componentmode representation for structural dynamic analysis,” AIAA Journal, vol. 13, no. 8, pp. 995–1006, 1975. View at: Publisher Site  Google Scholar
 W. L. Wang and Z. R. Du, Structural Vibration and Substructuring Method, Fudan University Press, Shanghai, China, 1985.
 Y. T. Leung, “An accurate method of dynamic substructuring with simplified computation,” International Journal for Numerical Methods in Engineering, vol. 14, no. 8, pp. 1241–1256, 1979. View at: Publisher Site  Google Scholar
 W. Yun, G. Duan, and Z. Hu, “The superelement method of modal synthesis and its application to dynamlc calculation of ship structure,” Journal of Shanghai Scientific Research Institute of Shipping, no. 1, pp. 1–18, 1981. View at: Google Scholar
Copyright
Copyright © 2016 Zhijun Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.