1. Introduction
Hydropower, as a green, renewable, and clean energy, is of great significance in promoting economic development, reducing pollution, and reducing greenhouse gas emissions [
1,
2,
3], and has been widely used worldwide [
4]. However, hydropower development in high-altitude areas is faced with severe conditions such as high velocity and large head variation, and a large number of high-hardness sand particles contained in the water flow can easily lead to the wear failure of turbine runner blades operating in this area [
5]. When the runner blades are worn, the fluid exciting force on the turbine runner system will change, which easily induces abnormal vibration of the main shaft system of the hydropower units, leading to security risks and even economic losses [
6]. In order to avoid the failure of hydropower units caused by severe blade wear, it is necessary to judge the wear degree of turbine runner blades in time. Therefore, the vibration characteristics of the shaft system of mixed-flow hydropower units are analyzed in this paper, considering the different wear degrees of the runner blades; the wear degree of underwater blades can be judged from the changes of shafting vibration characteristics that can be monitored.
The wear of turbine runner blades is a complicated problem involving hydraulics, mechanical engineering, and material science. Kang Min-Woo et al. [
7] found that the wear degree near the runner outlet was serious, and the wear rate was almost linearly positively correlated with the amount of silt. Chen et al. [
8] found that when cavitation corrosion acts on the surface of ductile metal materials, it will lead to plastic deformation and fatigue damage, which will further aggravate the wear of the materials. Mack R et al. [
9] found that the runner wear degree is highly correlated with the solid particle size, and a change in the shape parameter of the movable guide vanes will affect the fluid pressure pulsation. According to the studies of various scholars, wear is usually affected by a variety of factors, and the wear degree is most serious at the water edge of the blade [
10,
11]. Therefore, based on the above research results and practical engineering problems, it is necessary to take the wear degree of the runner blade into consideration when analyzing the fluid exciting force of the turbine runner during its long life cycle, so as to obtain a more accurate change trend of the fluid exciting force.
The dynamic behavior of the main shaft system is greatly affected by the fluid exciting force of the turbine runner system, and the changing trend of the force is the premise of the analysis of the unit dynamic characteristics. At present, the research on the internal flow field of hydraulic turbines mainly focus on the pressure pulsation of flow fields and the fluid exciting force. Chiappa et al. [
12] found that the vibration characteristics of the Karman vortex band are the key factors leading to blade erosion and fatigue failure. Grein and his team [
13] carried out detailed numerical simulation and experimental research on the blade vortex, analyzed the vibration causes of the turbine, and proposed the vibration reduction methods of the unit. Vu et al. [
14,
15] found that there is a large deviation between the CFD calculation results and the test results when deviating from the design working conditions. Therefore, proper analysis and restriction of the working conditions should be carried out in the numerical simulation. Rodriguez et al. [
16] found that the frequency components of static and static interference are mainly blade passing frequency and its harmonics. Zhai et al. [
17] found that in the low-load working zone, the pressure distribution of the gap between the upper crown and the lower ring was obviously uneven, resulting in a large gap imbalance force. Saeed et al. [
18] found that the pressure difference between the inlet and outlet of the runner and the pressure difference between the pressure side and the suction side of the blade was large, which led to the unbalanced force and torque on the turbine main shaft. Although the pressure pulsation and exciting force of the internal flow field of the hydraulic turbine have been studied deeply, the influence of the blade wear degree on the flow field is rarely taken into account. Therefore, this paper will analyze the flow field of runner blades under different wear degrees, explore the relationship between the fluid exciting force and runner blade wear degrees, and provide a basis for the subsequent analysis of the shafting vibration characteristics of hydropower units in this paper.
The shafting vibration of hydropower units belongs to the category of rotor dynamics and is a research field involving nonlinear vibration caused by complex exciting forces in the operation of the main shaft system. Bai et al. [
19] found that the random fluid exciting force has a great impact on the water guide bearing and runner system, but it is mainly manifested in vibration form and stability, rather than vibration amplitude. Zhang et al. [
20] analyzed the bend-torsional coupled vibration response of the shafting of hydropower units, but only considered the linear steady fluid exciting force on the runner, without taking into account the transient flow characteristics of the turbine. Based on the Lagrange equation, Mokhtar et al. [
21] studied the influence of parameters such as clearance and the friction coefficient at different rotational speeds and found that some parameter values would produce special friction-related features. Based on the Lagrange equation, Li et al. [
22] established a mathematical model of coupling the adjustment system and the spindle system and analyzed the influence of the gyroscopic effect and axial deviation on the stability of the system. Zhang et al. [
23] found that when the nonlinear sealing force model is adopted and there is no other excitation influence, the axis locus of the runner does not exceed the seal clearance even if the runner is unstable, and the friction probability between the runner system and the seal system is small. An et al. [
24] derived the energy equation of the runner and blade of a hydraulic turbine and found that with the continuous increase in the mass of the runner blade, the bifurcation phenomenon would occur in the runner system. Zhuang et al. [
25] considered the combined effects of hydraulic, mechanical, and electromagnetic forces and found that the hydraulic instability caused by changes in the blade outlet flow angle, outlet diameter, and guide blade opening determines the overall trend of the dynamic characteristics of shafting. However, the influence of the fluid exciting force change caused by the blade wear degree on spindle vibration characteristics was not considered in the present study.
Therefore, the relationship between the turbine blade wear degree and the vibration response of the shaft system was studied in this paper. A turbine model with blades exhibiting various degrees of wear was developed to analyze the effects of blade wear on the fluid exciting force. The exciting force was subsequently integrated into a dynamic model of the hydropower units shaft system to analyze the vibration response under different wear conditions. A mathematical relationship between the blade wear and the vibration response of the units was established. This relationship was intended to utilize measurable vibration characteristics to assess the wear level of turbine runner blades, providing a foundation for future fault prediction in hydropower units.
3. Dynamics Modeling of Shafting of Hydropower Units
The main shaft-system structure of the hydropower-generating units is shown in
Figure 8a.
O1 and
O2, respectively, represent the geometric centers of the generator rotor and the turbine runner, and
B1,
B2, and
B3, respectively, represent the geometric centers of the upper guide bearing, the lower guide bearing, and the water guide bearing. The centroid of the rotating part is shown in
Figure 8b.
is used to represent the centroid coordinates of the generator rotor and the turbine rotor, respectively, and
represents the centroid coordinates.
In order to facilitate the analysis, this paper adopts appropriate simplified assumptions: (1) In the modeling stage, the torsional motion and axial vibration of shafting are ignored, and only the radial displacement and inclination of the spindle system are analyzed. (2) Assume that all components of the spindle system behave isotropic during operation. In the rotation process of the disk, the angle of its rotation is defined as, where ω is the angular speed of the disk, and t is the time variable. Based on the coordinate function of the center of mass, the velocity coordinate function can be obtained by obtaining the first derivative of time, and the general expression of the translational kinetic energy of the center of mass of a disk can be derived:
When considering the tilt of the shafting of hydropower units, the simplified expression of the rotational kinetic energy can be obtained according to the rotor dynamics and the principle of conservation of the momentum moment [
22]:
By combining Equations (7) and (8), the total kinetic energy equation of the shafting of hydropower units can be obtained, as follows:
where
m1 and
m2 are the masses of the generator rotor and the turbine runner, respectively.
Jp is the polar moment of inertia.
Jd is the diameter moment of inertia.
φ is the rotation angle of the rotor around the axis of rotation.
The potential energy equation of hydropower units’ shafting is [
27]:
In Equation (10), there is:
where
A1 = (
a +
b)(
b +
c +
d),
A2 =
ab,
B =
b(
b +
c +
d),
k1,
k2,
k3 are the stiffness coefficients of the upper guide bearing, the lower guide bearing, and the water guide bearing, respectively.
By combining Equations (9) and (10) and following the principle of Lagrange mechanics, the Lagrange function of the main shaft system in hydropower generating units can be constructed, and its expression is:
Through the previous analysis, the lateral displacement and deflection angle of the main shaft system of the hydropower generating units are selected as generalized displacement coordinates, namely
qi = (
x1,
y1,
x2,
y2,
θx1,
θy1,
θx2,
θy2). The Lagrange equation can be expressed as:
where
L is the Lagrange function,
qi is the generalized coordinate, and
Qi is the generalized force.
The external forces acting on the shafting of hydropower units can be described as:
where
c1,
c2, and
c3 are the rotor damping coefficient, turbine runner damping coefficient, and structure bending damping coefficient, respectively.
k4 and
k5 stiffness coefficients of the large axis inclination caused by force and torque, respectively, represent the ratio of torque change caused by the change in unit inclination when the main axis is tilted.
Fx and
Fy are the fluid excitation forces acting on the x and y directions of the turbine runner.
Fx-ump and
Fy-ump are the unbalanced magnetic forces acting on the x and y directions of the generator rotor.
In Equation (14), the nonlinear unbalanced magnetic tension [
28] can be expressed as:
where
R is the radius of the generator rotor.
L is the length of the generator rotor.
kj is the fundamental wave coefficient of the magnetomotive force.
Ij is the generator excitation current.
μ0 is the permeability of air.
γ is the rotation angle, and
.
e is the radial displacement of the generator rotor,
.
ε is relative rotor eccentricity.
δ0 is the average air gap length of the rotor.
is the air gap permeability, which can be expressed as:
By substituting Equations (12) and (14) into the Lagrange equation, the differential equations of the dynamic motion of the spindle system can be obtained:
where
mi is the rotor mass,
xi and
yi are the radial displacement coordinates,
θi is the inclination angle,
J is the moment of inertia,
K11 and
K22 can be understood as the generalized stiffness, and the expressions are as follows: