1. Introduction
A great amount of energy is consumed in the course of iron and steel manufacture. To recover the waste heat and gas in the production process and to convert them into reusable energy, the TRT is used in the blast furnace for energy recovery [
1,
2]. The TRT uses the blast furnace gas generated in the operation process of the blast furnace to push the turbine, and the turbine drives the generator to generate electricity. The aerodynamic force is converted into mechanical energy by using the aforesaid method, and the generated electric energy is supplied to in-plant equipment to recover energy.
In an axial turbine, the risk of blade damage is considerable. In addition to reducing the efficiency of the turbine, the failed blade may hit other blades and cause greater losses. Since the TRT system is an axial turbine, it is important to avoid the failure of the blades during the operation of the system. During the operation of the system, in addition to pure metal fatigue failure, the damage may also be caused by metal corrosion defects, and failure due to the expansion of defects due to resonance [
3]. In practice, the locations where axial turbine blades are prone to failure are the bladed disk blades and the upper edges of the dovetail grooves of the blades. Wei et al. [
4] found that a large number of cracks occurred in the root of the blade in the steam turbine, and even the blades that failed due to this. After destructive and non-destructive inspections, it was found that the initial defects of the blades were caused by corrosion, which eventually caused cracks to grow under the influence of various excitations during operation, and finally caused the blades to fail. Poursaeidi et al. [
5] found a failed blade in a gas turbine, and the damage location was at the center of the blade. After the modal analysis by the FEM, it is found that the blade has a chance to excite the third mode of the blade under working conditions, which will cause the fatigue failure of the blade. Shevkhlari et al. [
6] tested the equivalent working conditions of the blades of the steam turbine. After an equivalent test of 500 working hours, the blades failed at the dovetail slot. The reason is that when the blade is manufactured, the geometry of the dovetail groove of the first section of the blade is different from the design, and the stress concentration occurs there, which eventually leads to the occurrence of fatigue failure. Katinic et al. [
7] found blade failure resulted from defects at the dovetail groove of a steam turbine blade. The defects are caused by the corrosion effect of the dovetail groove of the blade. Metal fatigue occurred due to structural vibration during operation. Zhao et al. [
8] found signs of fretting wear at the dovetail groove of the failed blade, and, at the operating speed, the second mode of the blade aggravated the crack growth at the wear part, and finally caused structural failure. In a TRT system, blast furnace gas must be received from the blast furnace, and its working conditions are more complicated. Liu et al. [
9] discussed the causes of TRT blade failure in their research and found that fretting and corrosion occurred at the dovetail groove of the blade at the same time. Finally, under structural vibration, cracks grew until the blade failed.
The centrifugal force generated during the operation of the rotor is also one of the keys to the damage of the blade structure. Cano et al. [
10] used the FEM to analyze the cenrtifugal force of blade members, and the assessment of corresponding blades service life can then be carried out. The results of the study show that the centrifugal force during the operation of the structure contributes to the generation of cracks on the blades, and more closely leads to the growth of cracks. The fatigue failure of turbine blades is gradual. The initial blade defects cause cracks to grow gradually, causing blade failure during the process. Rodriguee et al. [
11] used a vibration exciter to analyze the natural frequency of a cracked blade and found that when the blade has a crack, the natural frequency of the structure will change, and the position of the crack will change, and the blade vibration amplitude will also change. The longer the crack length, the greater the drop in the natural frequency of the structure. In addition to using vibration exciters for operating modal analysis of the blade structure, the EMA is also a common method to explore the dynamic characteristics of the structure. Zhang et al. [
12] used experimental modal analysis to explore the dynamic characteristics of Kaplan Turbine blades, and compared them with the frequency results of the Theoretical Modal Analysis (TMA) based on FEM to improve the reliability of the finite element model. Wu [
13] used experimental modal analysis methods to explore the dynamic characteristics of large turbines. After knowing the dynamic behavior of the system, using the disturbance behavior of the jet airflow to the blades, the frequency that may be excited by the external force is established, and is helpful to the discussion of blade failure analysis.
From previous studies, it can be known that the failure behavior of blades is mostly caused by initial defects. Under the external excitation of the structure, cracks are generated at the defects, and the cracks grow gradually. Therefore, understanding the structural dynamic characteristics of the blade is an important point that cannot be ignored in the rotor structure. The Campbell diagram is a method used to evaluate the dynamic characteristics of the axial flow turbine blade disc structure. The Campbell diagram can effectively determine whether the structure will resonate at the working speed. To explore the structural mode of the system resonance, Singh et al. [
14] proposed a Singh’s Advanced Frequency Evaluation (SAFE) diagram, which integrates the structural modes and resonance frequency information to determine the structural modes that may be excited. Bertini et al. [
15] used Campbell diagrams and SAFE diagrams to discuss the dynamic characteristics of the impeller. The results indicated that for the conditions under which the structural mode is excited, except for the same natural frequencies of the structure, the excitation phase of the external force needs to be consistent with the structural modes; otherwise, the structure will not resonate with the excitation. To monitor the moving blades of the turbine effectively, some rotor systems will be equipped with a Blade Tip Timing (BTT) system. The BTT technology can be used to measure the displacement of the tip of the rotor blade, and then analyze the dynamic characteristics of the rotor blade during operation. Madhavan et al. [
16] used the BTT system to monitor the blade dynamics and found abnormal blade vibration during the experiment. After shutting down the inspection, they found that the blade heel was cracked, which indicated that the BTT system was effective in monitoring the dynamic characteristics of the blade.
In the past studies, more attention was paid to the dynamic characteristics of the blades, and the simplified equivalent blisk/blade model was used for discussion. In this paper, a more complete model verified by EMA is used to implement strength analysis of the bladed disks of the turbine in rated working conditions. We firstly performed EMA on the practical TRT rotor, and then used modal verification technology to verify and update the finite element model. The centrifugal force analysis is employed to explore the centrifugal force that the TRT rotor structure bears during operation, evaluate the strength of the structure, and verify the stress generated by the centrifugal force at the blade and the analytical solution. In addition to the dynamic characteristics of the blade itself, we also pay attention to the contact behavior at the dovetail groove of the blade. After the verification is completed, the aerodynamic force of the TRT rotor structure is taken into consideration to evaluate the stress of the structure under actual working conditions. In the end, we used Campbell and SAFE diagrams to explore the dynamic characteristics of the TRT bladed disks to estimate the mode shapes of the TRT rotor that may be excited under working conditions.
2. Materials and Methods
The analysis of this study is the steelworks of a two-stage dry-type TRT comprising the stationary blades of two stages and moving blades of two stages. When the blast furnace top pressure gas flows into the TRT system, it passes by the 1st stage stationary blades with variable guide vanes at first. The purpose is to control the flow velocity of gas by adjusting the guide vanes to adjust the generating power of the TRT system. After passing by the 1st stage stationary blades, the air flows by the 1st stage moving blade and generates a lift to push the TRT rotor. When the air passes by the 1st stage blades, the flow velocity and pressure decrease to make the 2nd stage moving blade generate the torque equal to that of the 1st stage blade. The air flows by the 2nd stage stationary blades at first to increase the flow velocity. As the gas pressure decreases, to achieve the torque identical to that of the 1st stage, the 2nd stage rotor needs a larger flow area to reduce energy loss. Therefore, the 2nd stage moving blade is longer than the 1st stage moving blade.
This study aims to analyze the strength of the TRT moving blade in working conditions. To improve the reliability of the finite element model of moving blades, we perform modal verification to check the consistency of the actual structure and the corresponding finite element model. We use the MAC to evaluate the correlation between the mode shapes obtained from Theoretical and Experimental modal analysis to check the similarity of the mode shapes. It is difficult to perform in-situ measurements for TRT under working conditions; therefore, we use Computer-Aided-Engineering (CAE) software for analyzing centrifugal force and an aerodynamic force to discuss the stress distribution on the structure when the TRT rotor is in operation condition, and to evaluate the probable initial region if the structure fails.
2.1. Modal Verification of TRT Blades
To investigate the dynamic characteristics of TRT rotor blades, we used commercial software ANSYS to perform modal analysis, centrifugal force analysis, and aerodynamic analysis of TRT rotors through finite element analysis. To improve the reliability of the finite element models, the mode shapes and natural frequencies in the modal analysis of the TRT rotor will be captured and compared with those of the practical structures. To verify the agreement between the mode shapes of the finite element analysis and those of the actual structure, the MAC is introduced here to confirm the reliability of the finite element model.
2.1.1. Theoretical Modal Analysis of TRT Rotors Based on FEM
To implement modal verification, we used commercial software ANSYS to create a finite element model of the TRT rotor, which is employed in conjunction with appropriate boundary conditions and load conditions for computer-aided engineering. For the theoretical modal analysis of the TRT rotor, the free-free boundary condition is set, and an appropriate impulse force is applied at the chosen excitation point. Its purpose is to simulate the response of each measurement point of the TRT rotor system subjected to a broadband excitation, and the natural frequencies and mode shapes of the structure can then be obtained from the modal analysis based on finite element method.
2.1.2. Experimental Modal Analysis of TRT Rotors
To obtain the modal parameters of the TRT bladed disk structures, an impact hammer is employed to give a pulse excitation to the TRT rotor, the purpose of which is to excite the dynamic characteristics of the TRT blades. We set up an accelerometer at the measuring point, and use the accelerometer to collect the response data generated by the excitation of the TRT rotor blade. After obtaining the excitation and response signals, the excitation and response data are incorporated into the Brüel & Kjær RT Pro Photon 7.0 data acquisition system, and the Fast Fourier Transform (FFT) of correlation function operations are performed, and then the Frequency Response Function (FRF) of the system is obtained. The experimental procedure is as shown in
Figure 1.
The response of the system contains the information of excitation and the dynamic characteristics of the system. When the excitation frequency is close to the natural frequency of the system, the system has the largest magnification, which is regarded as resonance. However, the inner product of the mode shape and the excitation vectors is equally important to the contribution of resonance, and its physical meaning represents the similarity between the excitation phase (distribution in space) and the phase of the mode shape. When the two phases are similar (or the same), the inner product is larger up to unity; when it is completely orthogonal, the inner product is zero [
12]. Therefore, to excite the structural modes with mutually orthogonal direction components effectively, we apply excitations in the axial and radial directions, respectively, to excite the structural modes with axial or radial direction components effectively. In the experiment, since the signal measurement is often susceptible to environmental vibration or contaminated with other noise, it is impossible to guarantee whether the measured response is all induced from the excitation, so the coherence function
is employed as the importance index for a combination of the FFT of correlation functions between the response and excitation data, which is defined as
where
,
and
are, respectively, the auto power spectral density function of the system responses, the auto power spectral density function of the input excitation, and the cross-power spectral density function between the response and the excitation.
is close to unity, which means that the input and output data have a fairly good correlation. On the contrary, the closer the value is to zero, it means that the correlation between input and output data is not good. After evaluating
to estimate the correlation between the response and the excitation data, we evaluate the frequency response function (FRF)
of the system by using the combination of the power spectrum density function
and
defined as follows:
The FRF combined with the auto- and cross-power spectrum used in this study is one of the frequency-domain modal identifications of the structure under ambient vibration. This method is based on the assumption that the input signal is stationary white noise under the condition of unknown excitation, and the characteristics of the correlation function are employed to estimate the modal parameters of the system by evaluating the auto-power spectrum of the structural response data and the magnitude and phase of the cross-power spectrum with the corresponding response data of a reference channel.
2.1.3. Modal Assurance Criterion (MAC)
To ensure the consistency of the mode shapes obtained from the theoretical modal analysis and the experimental modal analysis, we use the modal assurance criteria (MAC) proposed by Allemang [
17] and defined as follows:
where
and
are the mode-shape vectors obtained from theoretical modal analysis and experimental modal analysis, respectively.
and
are, respectively, transpose and conjugate complex operators of the vector or matrix. Using the
indicator can verify the agreement between the mode-shape vectors obtained from finite element analysis and experimental modal analysis. If the
value approaches zero, it means that the correlation between the two modes is low; on the contrary, if the
is close to unity, it means that the correlation between the two modes is high.
2.2. Analysis of Centrifugal Force of TRT Rotor
When the TRT rotor is in operation, it must bear the centrifugal force resulting from rotation. The centrifugal force will occur in the radial direction of the rotating shaft, so that the blade bears the tensile stress resulted from centrifugal force during rotation, as shown in
Figure 2.
To discuss the influence of centrifugal force on the TRT rotor blade, two techniques are used for analysis and validation. First, the analytic solution is used for the determination of tensile stress of blade, and the tensile stress caused from centrifugal force of blade is analyzed by using the FEM. The analysis results are compared with the analytic solution, and then verified the reliability of centrifugal force analysis.
2.2.1. Analytical Solution
The TRT rotor working speed is 3600 rpm (60 rps). To obtain the tensile stress distribution on the blade due to centrifugal force, the centrifugal force (
) shall be calculated which is expressed as follows:
where
is the object mass;
is the angular velocity of an object rotating along the shaft;
is the shortest distance between the object centroid and rotation axis. According to the above equation, the magnitude of centrifugal force is correlated with the rotational speed of the object and related to the mass and centroid position of the object. As the length direction of the blade is normal to the direction of centrifugal force, there are different centrifugal forces in different length positions of the blade, expressed as follows:
where
is the distance from the center of circle
O to the blade tip;
is the distance from the center of circle
O to the blade root;
is the distance from the center of circle
O to the blade section of the desired centrifugal force;
is the blade material density;
is the angular velocity of rotation of the rotor;
is the function of the locational relation between the cross-section area and blade length;
is the function of the locational relation between the centroid position and blade length. The aforesaid symbols and blade diagram are shown in
Figure 3. The centrifugal force on the section of arbitrary blade length positions can be obtained through the aforesaid procedure.
When the centrifugal force in each length position of the blade is obtained, the tensile stress (
) in various positions of the blade can be calculated which is expressed as follows:
where
is the tensile stress resulting from centrifugal force;
is the centrifugal force; and
is the cross-sectional area. The tensile stress resulting from the centrifugal force in any radius position of the blade can be obtained by Equations (4)–(6).
2.2.2. Centrifugal Force Analysis of TRT Rotor Based on FEM
To discuss the stress intensity resulting from the centrifugal force in various positions of the TRT rotor blade, the CAE software is used for analyzing the centrifugal force of the TRT rotor. In the centrifugal force analysis, the load and boundary conditions are set up according to actual working conditions. The shaft speed is set as 3600 rpm, and the axial and radial constraints are given where the shaft bearing is fixed. After analysis, the stress of each section is extracted and compared with the analytic solution; the reliability of analysis is enhanced; and the influence of maximum stress on the structure is evaluated.
2.3. Analysis of Aerodynamic Force of TRT Rotor
The TRT system uses blast furnace top pressure to push the rotor blade, and the shaft drives the generator to generate electricity. When the rotor is in operation, the blast furnace gas flows across the flow channel between stationary blades, and the blast furnace gas is transferred to the moving blade. In this working condition, the moving blade moves among different flow channels due to the rotation of the shaft, so the external force on the moving blade is a periodic function. In actuation conditions, the aerodynamic force on the moving blade is a periodic function, and the aerodynamic load change can be expressed as follows.
where
is the external force on the moving blade;
~
are the magnitude of aerodynamic force flowing through the flow channel between different stationary blades when the rotor is in operation;
is the frequency of the moving blade passing through the flow channel when the rotor is in operation; and
~
are the phase difference of aerodynamic force between flow channels. According to Equation (7), the aerodynamic external force of various moving blades in different time conditions can be calculated as long as the magnitude of the aerodynamic force, the frequency of the rotor blade passing through the flow channel, and the phase difference of aerodynamic force between flow channels can be obtained.
The phase difference of aerodynamic force between flow channels has resulted from unequal numbers of moving blades and flow channels of the TRT rotor. Therefore, the phase difference is related to the geometry of the moving blade and flow channel. The law of phase difference repetition can be obtained by calculating the greatest common factor between the number of moving blades and the number of stationary blades, and the phase difference is calculated in the minimum unit. The minimum phase difference
is calculated as the following equation:
where
is the number of moving blades;
is the number of stationary blades; gcf(
,
) is the greatest common factor of moving blades and the number of stationary blades. After the angle of minimum phase difference and the law of phase difference repetition are obtained, the phase difference of each blade can be calculated by the following equation.
where
is the aerodynamic phase difference of the m-th moving blade;
is the minimum phase difference;
is the repetition number of the phase difference. To obtain the ordinal number between the blade and flow channel, the following equation is used for computation.
Equation (11) uses the geometrical relationship between the moving blade and flow channel to calculate the ordinal number between the moving blade and flow channel. is the parameter of discriminant; m is the desired blade number; gcf(,) is the greatest common factor of the moving blade and the number of stationary blades; is the number of moving blades; mod is the remainder operator; the calculated constituent ordinal number is substituted in Equation (9) to obtain the aerodynamic phase difference of moving blade m.
In the aerodynamic load change (Equation (7), besides the phase difference of aerodynamic force between flow channels (
~
), the frequency
of the moving blade passing through the flow channel when the rotor is in operation must be known. The frequency of moving blades passing through the flow channel can be obtained through the following procedure.
where
is the frequency of the blade passing by the nozzle when the rotor is in operation;
is the number of flow channels of stationary blades on the circumference; and
is the working speed. The calculation results of Equations (9) and (12) are substituted in Equation (4) to obtain the aerodynamic force parameter on the moving blade.
The finite element aerodynamic analysis is used for discussing the stress behavior of the shaft and blade of the TRT rotor under the effect of the aerodynamic force in working conditions. In the analysis, the loading conditions are set up according to actual working conditions; a simple harmonic external force is applied to the face of each blade as per the result of Equation (4). Besides loading conditions, the boundary condition setting for the TRT rotor in practice is of great importance, given the boundaries’ aim to simulate the constraints of the TRT rotor in the actual structure. The main referential boundary is the interface of the connection between the bearing on the rotor and the generator.
2.4. Campbell and SAFE Diagrams of TRT Bladed Disk
When the TRT rotor is under operation, the frequency in rotation obtained from order analysis will be generated due to the movement of the shaft. The frequency in rotation obtained from order analysis will vary with the number of blades. When the order frequency is close to the natural frequency of the structure, the mode shape of the structure system will be excited, causing the structure to vibrate. To understand the possibility that the rotor is excited by the frequency in rotation obtained from order analysis when the rotor is running, the theoretical modal analysis is employed to explore the dynamic behavior of the rotor at different speeds, and we can then obtain the Campbell diagram containing the relationship between the modal frequency and speed of the system at different speeds. The purpose of the Campbell diagram is to evaluate whether the shaft and blades will be excited by the frequency in rotation obtained from order analysis under the working speed of the system, which will cause the vibration of the structural system.
The Campbell diagram can be used to explore the relationship between the natural frequencies of the TRT rotor and the frequency in rotation obtained from order analysis, but it is unavailable to obtain the mode shapes of the excited structural modes. To understand the information of mode shapes at the intersection of the order line and the modal frequency of structural systems further, the SAFE diagram is applied to explore the types of nodal diameters to corresponding structural modes by drawing the structural modes with different nodal diameters of all moving blades at different speeds. The information contained in the SAFE diagram includes the system frequency where resonance occurs and the number of nodal diameters of its mode shapes; SAFE-diagram-related information can be used to estimate the durability of the blades further.
4. Conclusions
This study discussed the aerodynamic force, centrifugal force, and maximum stress on the structure and the positions of occurrence using finite element analysis (FEA) when the first-stage and second-stage moving blades of the TRT rotor system are in rated working conditions. The dynamic characteristics and vibration behaviors of rotor blades are also investigated through Campbell and SAFE diagrams. To validate the effectiveness of the finite element models, the mode shapes and natural frequencies in the FEA-based modal analysis of the TRT rotor are captured and compared with those of the practical structures through experimented modal analysis (EMA). To verify the agreement between the mode shapes of the finite element analysis and those of the actual structure, the modal assurance criterion (MAC) is introduced here to confirm the reliability of the finite element model. We use the modal verification method to explore the reliability of the TRT rotor finite element model. Since the TRT rotor is a large structure and the blade surface attached to the accelerometer is curved, some modal recognition results are poor, but the key modal guarantee indexes of the modal shapes are all higher than 0.6, so the finite element model is regarded as reliable. The CAE software is used for centrifugal force analysis and aerodynamic analysis to discuss the stress behavior and distribution when the structure is in rated conditions. In centrifugal force analysis, the analytic solution of the average stress of the blade section was compared with simulation results. Their errors are less than 5%, and thus confirm the reliability of simulation from CAE. It is found in the simulation that the first section dovetail slot of the blade bears maximum stress when the rotor blade is in a rated working condition of 3600 rpm. The rotor system may have metal fatigue failure in the first section dovetail slot of the blade. To evaluate the possibility of structural vibration of the rotor blades under the working conditions of the TRT rotor structure, we gathered the modal parameters at each speed to obtain the Campbell and SAFE diagrams, and determined the modal frequencies and the corresponding modal shapes that may be excited. From the dynamic analysis of the TRT rotor structure, it is pointed out that the key structural resonance frequencies of the primary bladed disk are 1630.50 Hz, 1504.00 Hz, 1575.10 Hz, and 1543.50 Hz, all of which are the torsional modal behavior of the blade; the key structural resonance frequency of the secondary blade is 438.86 Hz and 944.43 Hz, which are the bending and torsion modes of the blade, respectively.