Study on Impeller Fracture Model Based on Vibration Characteristics and Fractal Theory
During the operation of centrifugal compressor, failure easily occurs in the presence of complicated external forces. The failure process characterizes with strong nonlinearity, and hence it is difficult to be described by conventional methods. In this paper, firstly, the cracks in different positions are described using crack fractal theory. The basic failure modes of the impeller are summarized. Secondly, a three-dimensional finite element model of the impeller is constructed. Then the von Mises stress under the centrifugal force is calculated, and the corresponding impeller failure process is simulated by “element life and death technology” in ANSYS. Finally, the impeller failure mechanism is analyzed. It can be found that the static stress is not the main cause of the impeller failure, and the dynamic characteristics of the impeller are not perfect because of the pitch vibration modes which appeared in the investigated frequency range. Meanwhile, the natural frequency of the impeller also cannot avoid the frequency of the excitation force.
Centrifugal compressors are widely used in various fields such as aviation, metallurgy, petrifaction, environmental protection, and pharmacy industry. Centrifugal impeller is a core component of centrifugal compressors, which is composed of hub, blades, and shrouds . With the development of processing technology, welding impeller has been extensively used. Particularly, more and more high-speed centrifugal impellers in recent years are welded. The failure of centrifugal impeller with high-speed rotation is very serious. It not only affects the normal function of machinery, but also threats the safety of operator.
Scholars at home and aboard have carried out significative research on the impeller failure modes, summing up the following two main reasons for the cause of the impeller failure fracture: one is its own defects or material processing defects and the other is the alternating stress fatigue. The material and manufacturing defects are well understood [2–8]. It is shown that alternating stress in internal compressor is mainly from unsteady flow and aerodynamic loads [5, 9–12]. Combining the preceding research, the fatigue life design theory is proposed [3, 4, 13, 14].
Fracture mechanics is developed in the recent decades. Based on the macroscopic continuum mechanics, fracture mechanics can be applied to describe propagation and unstability of crack.
Fractal theory can accurately describe the crack process. Since 1975, Benoit B. Mandelbrot put forward fractal theory; the fractal theory has been widely applied. Ritchie and Olff  studied the relationship between the space scaling and biodiversity, providing a model to calculate the species diversity by using habitat fractal dimension. Lizzi et al.  studied the emergency management on multiband fractal antenna wireless communication systems. Taylor et al.  researched Pollock watercolor by using fractal theory. Mirny  used the fractal theory to describe chromosome. The fractal structures are also widely observed in dynamical systems [19, 20] and can induce long irregular transients .
In order to accurately reveal the failure mechanisms of centrifugal compressor impeller, we use two basic fundamental concepts of fractal theory, that is, fractal dimension and fractal length, and to describe the impeller fracture cracks, based on the sufficient analysis of centrifugal compressor impeller failure modes, the centrifugal compressor impeller common failure modes are summarized. Combining with the actual situation, the three-dimensional finite element model of a centrifugal compressor impeller is constructed, and the failure mechanisms are simulated.
2. Impeller Failure Mode Analysis
2.1. Description on Impeller Cracks
Impeller suffered from a variety of force including centrifugal force, air force, and other exciting force which lead to failure. The process has strong nonlinearity, which is mainly caused by two reasons. First, the external load of impeller is nonlinear; second, impeller failure process is nonlinear, which is difficult to be described by the conventional methods. Since in the failure process, a variety of different scales of fracture, which is distributed around the failure surface, displays a strong similarity. Therefore, to use fractal theory to describe the failure process is an effective method.
The specific process can be described by two basic concepts, fractal dimension and fractal length. Assuming an impeller cracks into pars and one part is similar to another one with similarity ratio β, then the similar dimension can be expressed asand similar dimension does not need to be integers.
Mandelbrot gives the estimation model to calculate length of the fractal curvewhere is fractal curves’ Euclidean length, is the initial operation length of the fractal curve, is fractal curve’ scaling, and is the fractal dimension.
Based on the above fractal theory, with specific forms of centrifugal compressor, it is easy to present the impeller failure modes.
2.2. Centrifugal Compressor Impeller Failures Modes
Figure 1 is the sectional view for a centrifugal compressor. The rotor is constituted by shaft and three centrifugal impellers (shown in Figure 2). Centrifugal impeller is welded together by three bodies, that is, hub, blades, and shroud, which are made of carbon quenched and tempered steel. By a period of operation time, this centrifugal impeller failure occurs frequently (as shown in Figures 3(a), 3(b), and 3(c)). In some serious spots, the impeller could be broken into pieces (shown in Figure 4). According to the collected data from the scene, its life ranges from a few days to more than a year. Figuring out the mechanism and finally solving this problem are extremely important.
(a) The eyes of the shroud weld
(b) Hub weld
(c) Weld between the hub and the blade
There are three main failure modes based on the failure characteristics of the scene macrocracks.
Failure Mode 1. Macrocracks are distributed in the vicinity of blades and shroud welding arc zone (hereinafter referred to as point A), and the shroud radially cracks in the blade inlet side (shown in Figure 5(a)).
(a) Failure mode 1
(b) Failure mode 2
(c) Failure mode 3
Failure Mode 2. Macrocracks are distributed in the vicinity of blades and hub welding arc zone (hereinafter referred to as point B) and break circumferentially (shown in Figure 5(b)).
Failure Mode 3. Macrocracks are distributed in the vicinity of blades and hub welding arc zone (hereinafter referred to as point C) and crack outwardly from the blade along the direction of the weld (shown in Figure 5(c)).
By comparing the three failure modes mentioned above, the cracks show the feature of self-similarity and can be described by fractal theory, although the cracks locations are different. Fractal dimension and fractal length can be calculated using formulas (1) and (2).
3. Numerical Simulations of Impeller Failure Modes
3.1. Elastic Strain Energy Density Unit
According to [22–27], the research focuses on the deformation of a unit volume element under external force. Assuming that there is no heat exchange with the outside, that is, a closed system, the total energy generated from the external power is . According to the first law of thermodynamics, the total strain energy generated in the unit is where is strain energy dissipation and is releasable elastic strain energy of a unit.
Figure 6 shows the value of the element energy dissipation and releasable strain energy per unit volume. The area of indicates the energy consumed when the unit damage and the plastic deformation occur. The shaded area represents the releasable strain energy stored in elements, and is unloaded elastic modulus. From the thermodynamic view, the energy dissipation is unidirectional and irreversible, but the energy release is bidirectional, as long as certain conditions are satisfied.
We then define the amount of energy damage of per unit volume:where is material critical strain dissipation energy and it is a constant which is irrelevant to the state of stress and is determined by the material test.
Under any stress state, , when , the material loses strength; that is,
Strain energy density is defined as , and because the above theoretical formulas are derived based on unit volume, (3), (4), and (5) evolve as follows:where is the strain energy dissipation density of a unit, is releasable elastic strain energy density of a unit, and is the critical strain energy density dissipated when a unit loses strength.
According to (6) we can obtain
Alloy material failure usually goes through several stages as shown in Figure 7: elastic stage, yield stage, strength hardening stage, and local plastic deformation stage. In other words, before the material enters the yield stage, reaches the maximum value . can be obtained based on material test and the material loses strength when ; so formula (7) becomes
Equation (8) is the criteria of unit failure (strength loss). That is, the unit failure occurs when the total unit strain energy density reaches the maximum value for the sum of material critical strain energy dissipation density and releasable elastic strain energy density of a unit.
3.2. Finite Element Model of the Impeller
The three-dimensional finite element model of the impeller is built according to the design drawings, with a total of 19 blades, uniformly distributed in the circumferential direction, and can be described as cycle symmetry structure. In order to reduce the computation cost, symmetric cycle model is used. We can focus on 1/19 of the original model. Figure 8(a) is the 1/19 impeller segment finite element model that contains only one blade. Figure 8(b) is the whole finite element model. Assuming that there is no heat exchange with the outside of the impeller during operation, only the effect of the centrifugal at working speed of 5130 r/min and impact of airflow on import and export are considered.
(a) The 1/19 model
(b) The whole model
Impeller material parameters are shown in Table 1.
3.3. Finite Element Numerical Simulation
ANSYS provides the function of calculating element volume and strain energy , where
The average strain energy density of the element is .
Based on the failure criteria of element elastic strain energy density theory in Section 3.1, an extension of the impeller crack is simulated in this paper, and the elements located in the cracks remain in failure state.
“Element life and death technology” in ANSYS is not really adding or deleting unit again, and the stiffness matrix is multiplied by a small factor. Dead element’s load will be zero; it is ineffective against the load vector.
We only consider the steady centrifugal stress corresponding to the operating speed. Figure 9 shows the segments von Mises stress distribution during the operating speed of 5130 r/min. The maximum von Mises stress appears in a section near the mouth ring which is point A. Using “element life and death technology” in ANSYS, the largest strain energy density unit is gradually killed, and the crack extension path of the part section is simulated.
Firstly, the finite element model shown in Figure 6 is established based on the solid model and the finite element model has elements. After loading and calculating, the strain energy and volume per unit are obtained. Then average strain energy density of each element is calculated and arranged finally. The maximum is identified, and the number of the element is recorded and killed. Repeating the previous calculating and reloading process, a sequence of unit numbers is obtained. The units killed in the finite element model are reproduced and remarked with the red color. In this way, the crack initiation and extension are simulated. Numerical simulation process of the whole crack extension is shown in Figure 10, and, in order to simulate crack extension perfectly, the finite element model should be finely meshed.
Figure 11 to Figure 15 show the impeller cracking process using the finite element method. The red elements represent the killed elements, and in calculation represents the number of units cracked. From the simulation results, Crack initiation point is located on inlet side of shroud and blades weld (point A) (as shown in Figure 11); crack gradually expands along surface to the mouth turn and depth directions until the ring is disconnected completely (as shown in Figure 12). When the cracking unit number reaches 2000, the shroud is totally disconnected (as shown in Figure 13). It is the failure mode 1.
When the cracks first appear (as shown in Figure 14) at the hub inlet and blade welding joint (point B), cracks extend to both sides along the weld zone, and finally cracks in the circumferential direction adjacent blades are penetrating along cracks extension (as shown in Figure 15). It is the failure mode 2.
(a) The shroud
(b) The hub
Cracks extension on hub may also have another failure mode. When the weld quality of blades and hub is poor (lack of penetration, such as shown in Figure 16) or there are obvious welding defects near the weld joint of hub inlet and blades, the blade would be gradually torn off inside-out from the disc. It is the failure mode 3.
At the same time, from numerical simulation results, it can be found that the cracks have strong self-similarity, which is consistent with the reality. Therefore it is feasible to adopt the fractal theory to describe crack.
3.4. Impeller Failure Mechanism Analysis
3.4.1. Criterion for the Avoiding Rate
The impellers suffer centrifugal force due to high-speed rotation and the airflow of the impact load. Airflow pressure is constantly changing when compressor is operating which is regarded as the dynamic load for the impeller. A large number of studies show that there are two types of working stress that affect fatigue life of the impeller and the blade: steady and alternating stress. Any excessively high stress will result in reduced fatigue life. The steady stress can be reduced through structural optimization design, and the alternating stress is mainly reduced by modulating frequency.
For the excitation frequency that needs to meet certain avoiding rate as , the adjacent frequency differences and need to meet the following relationship:
Equation (10) is the criterion for the avoiding rate between frequency of exciting force and natural frequency. At the same time, it can be obtained that the maximum investigated frequency is not less than .
3.4.2. Case Study of the Impeller
During the impeller operation, there is always unsteady alternating stress adding to the steady work stress. The steady work stress is mainly due to the centrifugal force of the impeller and has been fully considered during the design. The unsteady alternating stress is mainly caused by random load in actual operation. It can be seen in Figure 9 that the impeller maximum von Mises stress is 673 MPa, and yield strength of the impeller material is 850 MPa in Table 1. The impeller von Mises stress at steady-state operation does not reach yield strength of material, and it can be seen that the static strength design is reasonable.
As can be seen from the sectional view of the compressor, there is a backflow structure in front of the third-stage impeller inlet and has 14 blades. The whole impeller has 19 blades. The natural frequency needs to meet at least 8% avoiding rate with the exciting force frequency. There are two types of airflow exciting force of the impeller as follows.
The exciting force caused by the backflow structure blades:
The exciting force caused by the impeller blades:
From Section 3.4.1, we at least need to consider avoiding rate when the natural frequency of impeller and exciting frequency range 0–1755 Hz.
The natural frequency value can be obtained by modal analysis. From Figure 17, it can be seen that, for exciting force 1, Hz and only Hz meet the requirements, but Hz < = 1394.3 Hz, which does not meet the avoiding guidelines formula (10); for exciting force 2, Hz; does not meet the avoiding guidelines formula (10). Original design natural frequency did not avoid the two types of impeller flow excitation force; it is the main reason of impeller failure.
By investigating the impeller vibration modes, there is a pitch vibration (shown in Figure 18) in the 17th mode with frequency of 1599.8 Hz. Reference  shows that umbrella vibration and pitch vibration will occur when impeller rigidity is insufficient only. So, it can be seen that the rigidity of impeller is insufficient. There is a design flaw of impeller.
(a) Shroud side
(b) Hub side
The failure process of a third-stage impeller is investigated in this paper. The cracks on different locations were described by the fractal theory. The crack extension is simulated using “element life and death technology,” together with the element strain energy density theory. The principle for the avoiding rate between frequency of exciting force and natural frequency is established, and the failure mechanism of the impeller is studied. The main conclusions are given as follows.
(1) The impeller cracks have different sizes but show a very strong self-similarity. The fractal theory can be used to describe the impeller crack, which lays the foundation for the analysis.
(2) By analyzing the strength of the impeller, the maximum von Mises stress is only 673 MPa as shown in Figure 9. It is far away from the yield strength of the impeller material. The design of impeller strength can meet the related requirements. The static stress is not the main cause of impeller failure.
(3) By analyzing vibration characteristics as shown in Figure 17, it can be found that the dynamic characteristics of the impeller are not perfect. The fact that the natural frequency of the impeller does not avoid the frequency of the exciting force effectively would cause vibration risk.
(4) By analyzing the impeller model as shown in Figure 18, the rigidity of impeller structure is insufficient, as the pitch vibration mode appears in the investigated frequency range.
(5) Comparing the impeller macroscopic fracture in Figure 16 with failure mode 3, it is found that the impeller has welding defects, such as imperfect welding, not penetration.
(6) It can be found that, by using strain energy density, “element life and death technology” together, the crack extension of the impeller can be well established.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
H. Y. Wu, Structural Parts of Turbine and Strength Calculation, Mechanical Industry Press, Beijing, China, 1982.
C. W. Wu, “Reliability analysis of large centrifugal compressor impeller blade fatigue,” Equipment and Manufacturing Technology, pp. 1–3, 2008.View at: Google Scholar
X.-Y. Yang, Z.-X. Geng, X.-H. Cai, and H.-Q. Shi, “Fracture analysis of the 2nd stage compressor blade in an aeroengine,” Journal of Aerospace Power, vol. 16, no. 4, pp. 327–330, 2001.View at: Google Scholar
L. Wei and Q. Yuan, “Centrifugal compressor surge and control,” Fan Technology, vol. 1, pp. 43–47, 2004.View at: Google Scholar
Y. J. Mao, D. T. Qi, and Q. Xu, “Numerical study on high cycle fatigue failure of a centrifugal compressor impeller blades,” Journal of Xi'an Jiaotong University, vol. 42, no. 11, pp. 1336–1339, 2008.View at: Google Scholar
C. W. Wu, C. K. Guan, and X. L. Guo, “Large centrifugal compressor impeller blade fatigue reliability analysis,” Equipment Manufacturing Technology, vol. 8, no. 1, article 3, 2008.View at: Google Scholar
X. Y. Yang and Z. X. Geng, “Study of a two engine compressor blade failure analysis following the crack,” Aerospace, vol. 16, no. 4, pp. 327–330, 2001.View at: Google Scholar
R. P. Taylor, A. P. Micolich, and D. Jonas, “Fractal analysis of Pollock's drip paintings,” Nature, vol. 399, no. 6735, p. 422, 1999.View at: Google Scholar
F. Erdogan and G. C. Sih, “On crack extension in plane loading and transverse shear,” Journal of Basic Engineering, vol. 85, no. 2, pp. 519–527, 1963.View at: Google Scholar
G. C. Sih, “Energy-density concept in fracture mechanics,” Engineering Fracture Mechanics, vol. 5, pp. 1037–1040, 1973.View at: Google Scholar
L. P. Gan, “Based on the total strain energy density roulette low cycle fatigue life and reliability of the model study,” Chinese Scientific Papers, vol. 7, no. 8, pp. 616–621, 2012.View at: Google Scholar
Q. Sun, S.-C. Li, X.-D. Feng, W.-T. Li, and C. Yuan, “Study of numerical simulation method of rock fracture based on strain energy density theory,” Rock and Soil Mechanics, vol. 32, no. 5, pp. 1575–1582, 2011.View at: Google Scholar
H. P. Xie, “Based on the energy dissipation and the rock strength and overall failure criterion release principle,” Rock Mechanics and Engineering, vol. 24, no. 17, pp. 3003–3010, 2005.View at: Google Scholar