# Dynamic Model for Axial Motion of Horizontal Pelton Turbine and Validation in Actual Failure Case

Authors: C. G. Rodriguez; D. Zambrano; S. Reyes; J. Tapia; M. Egusquiza; E. Egusquiza
Journal: Shock and Vibration
 (2020)

---

## Abstract

Pelton turbines are important machines for power generation from a renewable energy source such as water. For power rates below 20 MW, the rotor of Pelton turbines is usually in horizontal position. Considering ideal mounting and operating conditions, there are no axial forces acting on the rotor. In practice, there is an hydraulic force due to the difference between nozzle centerline and bucket centerline, and there is a magnetic force due to the difference between axial position of stator and rotor magnetic field centers. These forces are supported by bearings. In this article, a nonlinear dynamic model considering these axial forces and bearings behavior is presented and solved for two different actual Pelton turbines. The nonlinear dynamic model allows determining and evaluating the source of axial motion and therefore provides valuable information in order to reduce it when the axial displacement is high enough to produce damage.


---

## Body

## 1. Introduction

Pelton turbines are widely used all over the world for electric power generation. This type of turbines with horizontal rotors is commonly used for power rates below 20 MW for small hydro power plant which allows environmentally friendly, local, and stable power supply.There are several studies about the dynamic behavior of Pelton turbines and its failures. These studies are mainly focused on the performance calculation [1, 2]; the bucket erosion [3–6]; and the bucket fracture [7–9]. Computer fluid dynamics (CFD) simulations and most experimental studies are carried out considering that jet centerline coincides with bucket centerline. Taking into account this consideration, there are no axial forces and therefore there is no need to study the axial motion of horizontal Pelton turbines. However, in practice, it is found that, under some conditions, excessive axial motion can occur, producing, in some extreme cases, catastrophic failure of the bearing. The axial motion of an horizontal Pelton turbine is originated due to axial forces on the rotor, which can be originated due to two main causes (Figure 1):(i)
the axial force due to water jet on the bucket,fh(ii)
the axial force due to the magnetic field on generator rotor,fm.Figure 1
Horizontal Pelton turbine.Axial force on the bucket due to water jet is zero if the jet centerline coincides with bucket centerline. If both centerlines do not coincide, then there will be an axial forcefh due to uneven separation of water jet as schematically shown in Figure 2. Another axial force is due to the magnetic force on the rotor generator. This force should be zero if the axial position of the stator magnetic field center coincides with the axial position of rotor magnetic field center. If both axial positions do not coincide, then there will be an axial force fm on the rotor as schematically shown in Figure 3. These undesirable resultant axial forces on the rotor are supported by bearing reaction which can be degraded over the time. In the following study, a model for axial motion considering these forces will be presented and then it will be used in two actual cases of horizontal Pelton turbines: one in normal operating condition and another presenting excessive axial motion with bearing failure as a consequence.Figure 2
Force on bucket due to the difference between water jet centerline and bucket centerline.Figure 3
Force on rotor due to the difference between axial position of stator and rotor magnetic field centers.


## 2. The Model

### 2.1. Axial Force on the Bucket due to Water Jet

To determine the axial force generated by water jet on the bucket, let us consider that the centerlines of the jet and bucket are parallel but displaced from each other by a distancez−zn0, where z is the moving axial position of bucket plane and zn0 is the fixed axial position of nozzle centerline. Considering the diameter of water jet as d and that the jet is separated by bucket cut-out at z, then the amount of flow rate in each side of bucket exit will be QL at the left and QR at the right, as shown in Figure 4.Figure 4
Uneven separation of water jet in Pelton turbine bucket.Figure5 shows the jet separation where θ is the angle of intersection between water jet diameter and bucket cut-out that is given by θ=2cos−12z−zn0/d. Considering that v is the mean inlet relative velocity of water into bucket, then QL and QR are determined by the following equation:(1)QL=vπd24−θd28+z−zn0d2sinθ2,QR=vθd28−z−zn0d2sinθ2.Figure 5
Geometry of water jet separation in bucket.
(a)(b)Then, the axial force on bucket due to change in momentum can be determined as follows:(2)fh=14ρv2d2sinβπ−2cos−12z−zn0d+4z−zn0dsincos−12z−zn0d.In equation (2), ρ is the water density and β is the angle of water exit.The axial hydraulic forcefh in equation (2) is shown in Figure 6. It can be observed that for a positive displacement in z−zn0, there will be a force on the bucket in the same direction. This behavior is equivalent to a negative stiffness giving instability to axial rotor position.Figure 6
Axial force on the bucket due to difference between nozzle and bucket centerlines.The bucket’s relative position with respect to jet impact position changes as rotor rotates (Figure7); thus, there is a variation of force magnitude every time a bucket passes, which leads into a magnitude variation at the bucket passing frequency ZBΩ, with ZB being the number of buckets and Ω the rotating speed. This amplitude modulation is incorporated to the model using the modulation function Ah in equation (3). The shape of modulation function is based on CFD results of force and torque variation due to rotation as [10] (solid line in Figure 8 for single bucket as rotor rotates). The experimental data of [11] shows the same behavior but with a short duration peak when the bucket is first reached by jet (dashed line in Figure 8). The existence of this peak is not explained with CFD simulations and it is supposed to be originated in the way pressure is measured [12–14]. This peak in the modulation function would produce several harmonics of blade passing frequency. During prototype measurements, it is found that there are no high harmonics of bucket passing frequency >5ZBΩ in axial motion; therefore, peak consideration does not provide significant difference in axial calculated motion, and for that reason, the modulation function considers the first four harmonics of blade passing frequency ZBΩ,2ZBΩ,3ZB,Ωand4ZBΩ.Figure 7
Bucket relative position to water jet.
(a)(b)(c)(d)(3)Ah=1+∑n=14ahncosnZBΩt+∑n=14bhnsinnZBΩt.Figure 8
Torque as a function of angular position due to one bucket from CFD calculation and experiments (schematically from [11, 14]).In equation (3), ahn and bhn are the Fourier coefficients that are calculated to describe the hydraulic force modulation function (obtained from Figure 8 repeated for each bucket). After these considerations, the axial force on the buckets is obtained as follows:(4)fh=14ρv2d2sinβπ−2cos−12z−zn0d+4z−zn0dsincos−12z−zn0d∗1+∑n=14ahncosnZBΩt+∑n=14bhnsinnZBΩt.


### 2.2. Axial Force on the Rotor Generator due to Magnetic Forces

Magnetic fields in stator and rotor generator can be considered as symmetrical with respect to an axial position in both stator and rotor. The axial position of stator magnetic field center is fixed, while the axial position of rotor magnetic field center moves with the rotor. If both axial positions do not coincide, then there will be a magnetic force or magnetic pull that tries to maintain both magnetic fields centers at the same axial position.Figure9 shows the axial rotor magnetic field center at z+zrs that does not coincide with axial position of stator magnetic field center which is at zs0, with zrs being the axial distance between bucket centerline and rotor magnetic field center. The axial distance between both magnetic centers is z+zrs−zs0. The diameter of stator is ds, and the diameter of rotor is dr.Figure 9
Axial position of stator and rotor magnetic field centers.The magnetic force on rotor due to interaction of theZp poles as a consequence of different axial position is given by the following equation [15]:(5)fm=−12i2N2Zp∂Pz∂z.In equation (5), i is the current, N is the number of turns, and P is the permeance. Permeance is given by the following equation:(6)P=P1+P2.In equation (6), P1 and P2 are the permeances of the two parts in the generator air gap as shown in Figure 10, which are given by the following equations, where μ0 is the magnetic permeability of air.(7)P1=πμ0ds+drL−z+zrs−zs0ds−dr,(8)P2=μ0ds+drln1+πz+zrs−zs0ds−dr.Figure 10
Axial magnetic force calculation.Using equations (6) to (8) into (5), the axial magnetic force in equation (9) is obtained. Figure 11 shows magnetic force magnitude considering typical values of ds, dr, and L for a 1000 rev/min and 10 MW generator. References [16–18] obtain a similar behavior for axial force.(9)fm=−12μ0π2i2ZpN2ds+drds−drz+zrs−zs0ds−dr+πz+zrs−zs0.Figure 11
Axial force magnitude on rotor generator.As in the case of axial force due to jet, there will be a change in the magnitude of the magnetic force as rotor rotates due to relative position of poles. This change in magnetic force will be considered as the modulation functionAm in equation (10) at twice pole passing frequency (2ZPΩ, ZP being the number of poles). The shape of modulation function shown in Figure 12 is taken from electrical simulations of magnetic pull during rotation as [19]. The modulation function considers the first four harmonics of twice pole passing frequency 2ZPΩ,4ZPΩ,6ZPΩ,and8ZPΩ.Figure 12
Shape of magnetic force (schematically from Kang et al., 2017).(10)Am=1+∑n=14amncos2nZPΩt+∑n=14bmnsin2nZPΩt.In equation (10), amn and bmn are the Fourier coefficients that are calculated to describe time waveform of magnetic force modulation function in Figure 12.After these considerations, the axial magnetic force on rotor is given by(11)fm=−12μ0π2i2ZpN2ds+drds−drz+zrs−zs0ds−dr+πz+zrs−zs0∗1+∑n=14amncos2nZPΩt+∑n=14bmnsin2nZPΩt.


### 2.3. Axial Reaction Force on Bearing

The bearings are usually two and are mounted with the rotor centered by the nozzle centerline and bucket centerline at the same position. One of the bearings will be mounted in order to support axial motion in one direction and the other bearing will be mounted in order to support axial motion in the other axial direction. In each one of the bearings, there will be a working side and a nonworking side. The working side in each bearing is obtained giving an adequate axial clearance to form lubricant film between rotor and bearing surfaces as shown in Figure13. In some turbine designs, a single bearing has both working sides.Figure 13
Axial working and nonworking side of two bearings.If the clearance in the working side is higher than the one necessary to obtain an adequate lubricant behavior, then there will be an interval of axial motionΔzb0 where the bearing is not able to provide a reaction force. To incorporate this condition to model, a zero stiffness interval is incorporated to bearing behavior. If the clearance in the working side is necessary to obtain an adequate lubricant behavior, then the bearing behaves as a constant stiffness spring. To incorporate this condition to model, a stiffness of k1 is incorporated to bearing behavior. The value of k1 is obtained making use of bearing geometry and the calculations based on [20]. If the clearance in the working side is lower than the one necessary to obtain an adequate lubricant behavior, then there is a high axial force (in comparison to the one considered by design) and bearing could be damaged. To incorporate this condition to model, a high stiffness k2 is incorporated to bearing behavior. This behavior is equivalent to a single bearing with maximum displacement of shaft in axial direction lb=l2−l1 centered at an axial distance zb0 which is schematically shown in Figures 14 and 15.Figure 14
Equivalent bearing.Figure 15
Bearing axial stiffness.Taking the previous analysis into account, the bearing axial stiffness is considered as follows:(12)kb=0ifz+zrb−zb0<Δzb02,k1,ifΔzb≤z+zrb−zb0<lb2,k2=∞,ifz+zrb−zb0≥lb2.In equation (12), zrb is the distance between bucket plane and axial position of bearing in rotor. To take into account the bearing stiffness, it is not necessary to add a modulation function because bearings are flat, and no significant variation of stiffness is found in a single rotation of rotor.


### 2.4. Equation of Motion

Due to high torque and low rotating speed in hydraulic turbines, the rotor has a higher diameter in relation to its length when compared to other types of turbines as gas or steam turbines. For that reason, the rotor first natural frequency in flexion is usually much higher than its rotating speed. Moreover, the first natural frequency in axial compression and traction is much higher than the first natural frequency in flexion. For these reasons, the rotor is considered as a rigid body in axial direction. Simulations of horizontal Pelton turbines [21] show that axial natural frequencies are much higher than those of excitation forces. By international standards as API 670, the rotor should be considered as a rigid body if its main excitation frequency is below 0.7 times rotor natural frequency. In the case of the proposed model, the excitation in axial direction with higher frequency is the bucket passing frequency ZBΩ which must be lower than 0.7 times its first axial natural frequency fn,a to be considered as rigid. For the simulated cases, this condition is verified before computing the model that is represented by(13)∑f=mz¨,fh+fm-kbz+zrb-zb0-cbz˙=mz¨.In equation (13), m is rotor mass, and cb is an introduced bearing damping, considered in order to eliminate initial transient response. Considering the previous axial forces and bearing axial stiffness, the model equation of motion is given by equation (14) and schematically shown in Figure 16.(14)14ρv2d2sinβπ−2cos−12z−zn0d+4z−zn0dsincos−12z−zn0d∗1+∑n=14ahncosnZBΩt+∑n=14bhnsinnZBΩt+−12μ0π2i2ZpN2ds+drds−drz+zrs−zs0ds−dr+πz+zrs−zs0∗1+∑n=14amncosnZPΩt+∑n=14bmnsinnZPΩt−kbz+zrb−zb0−cbz˙=mz¨.Figure 16
Model of Pelton turbine axial motion.In Figure16, the conditions that are important for adequate axial behavior are as follows:(i)
The positions for zero axial force are equal forfh and fm, which means that zs0−zn0=zrs(ii)
The position for zero axial forces coincides with the axial position of zero reaction at bearing, which means thatzb0−zn0=zrb(iii)
Axial clearance of bearing is adequateThe nonlinear differential equation of motion in equation (14) was solved making use of MATLAB, using ode113 with a variable time step.


### 2.5. Model Variables to be Determined

Rotor dimensions, rotor material, flow rate, power rate, and generator current and voltage are known. Then, the model can be solved for different values ofzn0,zb0,zs0,zrb,zrs,Δzb0,andlb. These values are adjusted in order to obtain an axial behavior similar to the actual axial behavior.The valuezn0=0 is considered in simulations as the reference for axial position. For the first simulation for adjustment, it was considered that (i) the distance zrb is equal to the distance between bucket centerline and the axial position of bearing center Zrb determined by turbine specifications; (ii) the distance zrs is equal to the distance between bucket centerline and the axial position of generator rotor geometric center Zrs determined by turbine specifications; and (iii) Δzb0 is considered as zero which means that the axial clearance is supposed to be equal or lower than the one necessary to obtain an adequate lubricant film behavior. Taking into account these considerations for the first simulation, then the values of zs0,zb0,andlb are determined. If these values cannot provide a good agreement of model to experimental data, then the values of zrb,zrs, and Δzb0 are modified and adjusted. After this, the modulation functions for hydraulic and magnetic forces are modified until the best agreement with measurements is obtained. This procedure is schematically shown in Figure 17.Figure 17
Adjustment of model parameters.The adjustment was carried out making use of different measurements of turbine. It is useful to count on the axial motion measured rotating at constant speed, zero current, and zero voltage in generator: full speed no load (FSNL). This condition allows focusing onzb0 and lb because there is not magnetic force on rotor at this operating condition. Then, it is necessary to count on the axial motion measured at different power in order to focus on zs0.


## 2.1. Axial Force on the Bucket due to Water Jet

To determine the axial force generated by water jet on the bucket, let us consider that the centerlines of the jet and bucket are parallel but displaced from each other by a distancez−zn0, where z is the moving axial position of bucket plane and zn0 is the fixed axial position of nozzle centerline. Considering the diameter of water jet as d and that the jet is separated by bucket cut-out at z, then the amount of flow rate in each side of bucket exit will be QL at the left and QR at the right, as shown in Figure 4.Figure 4
Uneven separation of water jet in Pelton turbine bucket.Figure5 shows the jet separation where θ is the angle of intersection between water jet diameter and bucket cut-out that is given by θ=2cos−12z−zn0/d. Considering that v is the mean inlet relative velocity of water into bucket, then QL and QR are determined by the following equation:(1)QL=vπd24−θd28+z−zn0d2sinθ2,QR=vθd28−z−zn0d2sinθ2.Figure 5
Geometry of water jet separation in bucket.
(a)(b)Then, the axial force on bucket due to change in momentum can be determined as follows:(2)fh=14ρv2d2sinβπ−2cos−12z−zn0d+4z−zn0dsincos−12z−zn0d.In equation (2), ρ is the water density and β is the angle of water exit.The axial hydraulic forcefh in equation (2) is shown in Figure 6. It can be observed that for a positive displacement in z−zn0, there will be a force on the bucket in the same direction. This behavior is equivalent to a negative stiffness giving instability to axial rotor position.Figure 6
Axial force on the bucket due to difference between nozzle and bucket centerlines.The bucket’s relative position with respect to jet impact position changes as rotor rotates (Figure7); thus, there is a variation of force magnitude every time a bucket passes, which leads into a magnitude variation at the bucket passing frequency ZBΩ, with ZB being the number of buckets and Ω the rotating speed. This amplitude modulation is incorporated to the model using the modulation function Ah in equation (3). The shape of modulation function is based on CFD results of force and torque variation due to rotation as [10] (solid line in Figure 8 for single bucket as rotor rotates). The experimental data of [11] shows the same behavior but with a short duration peak when the bucket is first reached by jet (dashed line in Figure 8). The existence of this peak is not explained with CFD simulations and it is supposed to be originated in the way pressure is measured [12–14]. This peak in the modulation function would produce several harmonics of blade passing frequency. During prototype measurements, it is found that there are no high harmonics of bucket passing frequency >5ZBΩ in axial motion; therefore, peak consideration does not provide significant difference in axial calculated motion, and for that reason, the modulation function considers the first four harmonics of blade passing frequency ZBΩ,2ZBΩ,3ZB,Ωand4ZBΩ.Figure 7
Bucket relative position to water jet.
(a)(b)(c)(d)(3)Ah=1+∑n=14ahncosnZBΩt+∑n=14bhnsinnZBΩt.Figure 8
Torque as a function of angular position due to one bucket from CFD calculation and experiments (schematically from [11, 14]).In equation (3), ahn and bhn are the Fourier coefficients that are calculated to describe the hydraulic force modulation function (obtained from Figure 8 repeated for each bucket). After these considerations, the axial force on the buckets is obtained as follows:(4)fh=14ρv2d2sinβπ−2cos−12z−zn0d+4z−zn0dsincos−12z−zn0d∗1+∑n=14ahncosnZBΩt+∑n=14bhnsinnZBΩt.


## 2.2. Axial Force on the Rotor Generator due to Magnetic Forces

Magnetic fields in stator and rotor generator can be considered as symmetrical with respect to an axial position in both stator and rotor. The axial position of stator magnetic field center is fixed, while the axial position of rotor magnetic field center moves with the rotor. If both axial positions do not coincide, then there will be a magnetic force or magnetic pull that tries to maintain both magnetic fields centers at the same axial position.Figure9 shows the axial rotor magnetic field center at z+zrs that does not coincide with axial position of stator magnetic field center which is at zs0, with zrs being the axial distance between bucket centerline and rotor magnetic field center. The axial distance between both magnetic centers is z+zrs−zs0. The diameter of stator is ds, and the diameter of rotor is dr.Figure 9
Axial position of stator and rotor magnetic field centers.The magnetic force on rotor due to interaction of theZp poles as a consequence of different axial position is given by the following equation [15]:(5)fm=−12i2N2Zp∂Pz∂z.In equation (5), i is the current, N is the number of turns, and P is the permeance. Permeance is given by the following equation:(6)P=P1+P2.In equation (6), P1 and P2 are the permeances of the two parts in the generator air gap as shown in Figure 10, which are given by the following equations, where μ0 is the magnetic permeability of air.(7)P1=πμ0ds+drL−z+zrs−zs0ds−dr,(8)P2=μ0ds+drln1+πz+zrs−zs0ds−dr.Figure 10
Axial magnetic force calculation.Using equations (6) to (8) into (5), the axial magnetic force in equation (9) is obtained. Figure 11 shows magnetic force magnitude considering typical values of ds, dr, and L for a 1000 rev/min and 10 MW generator. References [16–18] obtain a similar behavior for axial force.(9)fm=−12μ0π2i2ZpN2ds+drds−drz+zrs−zs0ds−dr+πz+zrs−zs0.Figure 11
Axial force magnitude on rotor generator.As in the case of axial force due to jet, there will be a change in the magnitude of the magnetic force as rotor rotates due to relative position of poles. This change in magnetic force will be considered as the modulation functionAm in equation (10) at twice pole passing frequency (2ZPΩ, ZP being the number of poles). The shape of modulation function shown in Figure 12 is taken from electrical simulations of magnetic pull during rotation as [19]. The modulation function considers the first four harmonics of twice pole passing frequency 2ZPΩ,4ZPΩ,6ZPΩ,and8ZPΩ.Figure 12
Shape of magnetic force (schematically from Kang et al., 2017).(10)Am=1+∑n=14amncos2nZPΩt+∑n=14bmnsin2nZPΩt.In equation (10), amn and bmn are the Fourier coefficients that are calculated to describe time waveform of magnetic force modulation function in Figure 12.After these considerations, the axial magnetic force on rotor is given by(11)fm=−12μ0π2i2ZpN2ds+drds−drz+zrs−zs0ds−dr+πz+zrs−zs0∗1+∑n=14amncos2nZPΩt+∑n=14bmnsin2nZPΩt.


## 2.3. Axial Reaction Force on Bearing

The bearings are usually two and are mounted with the rotor centered by the nozzle centerline and bucket centerline at the same position. One of the bearings will be mounted in order to support axial motion in one direction and the other bearing will be mounted in order to support axial motion in the other axial direction. In each one of the bearings, there will be a working side and a nonworking side. The working side in each bearing is obtained giving an adequate axial clearance to form lubricant film between rotor and bearing surfaces as shown in Figure13. In some turbine designs, a single bearing has both working sides.Figure 13
Axial working and nonworking side of two bearings.If the clearance in the working side is higher than the one necessary to obtain an adequate lubricant behavior, then there will be an interval of axial motionΔzb0 where the bearing is not able to provide a reaction force. To incorporate this condition to model, a zero stiffness interval is incorporated to bearing behavior. If the clearance in the working side is necessary to obtain an adequate lubricant behavior, then the bearing behaves as a constant stiffness spring. To incorporate this condition to model, a stiffness of k1 is incorporated to bearing behavior. The value of k1 is obtained making use of bearing geometry and the calculations based on [20]. If the clearance in the working side is lower than the one necessary to obtain an adequate lubricant behavior, then there is a high axial force (in comparison to the one considered by design) and bearing could be damaged. To incorporate this condition to model, a high stiffness k2 is incorporated to bearing behavior. This behavior is equivalent to a single bearing with maximum displacement of shaft in axial direction lb=l2−l1 centered at an axial distance zb0 which is schematically shown in Figures 14 and 15.Figure 14
Equivalent bearing.Figure 15
Bearing axial stiffness.Taking the previous analysis into account, the bearing axial stiffness is considered as follows:(12)kb=0ifz+zrb−zb0<Δzb02,k1,ifΔzb≤z+zrb−zb0<lb2,k2=∞,ifz+zrb−zb0≥lb2.In equation (12), zrb is the distance between bucket plane and axial position of bearing in rotor. To take into account the bearing stiffness, it is not necessary to add a modulation function because bearings are flat, and no significant variation of stiffness is found in a single rotation of rotor.


## 2.4. Equation of Motion

Due to high torque and low rotating speed in hydraulic turbines, the rotor has a higher diameter in relation to its length when compared to other types of turbines as gas or steam turbines. For that reason, the rotor first natural frequency in flexion is usually much higher than its rotating speed. Moreover, the first natural frequency in axial compression and traction is much higher than the first natural frequency in flexion. For these reasons, the rotor is considered as a rigid body in axial direction. Simulations of horizontal Pelton turbines [21] show that axial natural frequencies are much higher than those of excitation forces. By international standards as API 670, the rotor should be considered as a rigid body if its main excitation frequency is below 0.7 times rotor natural frequency. In the case of the proposed model, the excitation in axial direction with higher frequency is the bucket passing frequency ZBΩ which must be lower than 0.7 times its first axial natural frequency fn,a to be considered as rigid. For the simulated cases, this condition is verified before computing the model that is represented by(13)∑f=mz¨,fh+fm-kbz+zrb-zb0-cbz˙=mz¨.In equation (13), m is rotor mass, and cb is an introduced bearing damping, considered in order to eliminate initial transient response. Considering the previous axial forces and bearing axial stiffness, the model equation of motion is given by equation (14) and schematically shown in Figure 16.(14)14ρv2d2sinβπ−2cos−12z−zn0d+4z−zn0dsincos−12z−zn0d∗1+∑n=14ahncosnZBΩt+∑n=14bhnsinnZBΩt+−12μ0π2i2ZpN2ds+drds−drz+zrs−zs0ds−dr+πz+zrs−zs0∗1+∑n=14amncosnZPΩt+∑n=14bmnsinnZPΩt−kbz+zrb−zb0−cbz˙=mz¨.Figure 16
Model of Pelton turbine axial motion.In Figure16, the conditions that are important for adequate axial behavior are as follows:(i)
The positions for zero axial force are equal forfh and fm, which means that zs0−zn0=zrs(ii)
The position for zero axial forces coincides with the axial position of zero reaction at bearing, which means thatzb0−zn0=zrb(iii)
Axial clearance of bearing is adequateThe nonlinear differential equation of motion in equation (14) was solved making use of MATLAB, using ode113 with a variable time step.


## 2.5. Model Variables to be Determined

Rotor dimensions, rotor material, flow rate, power rate, and generator current and voltage are known. Then, the model can be solved for different values ofzn0,zb0,zs0,zrb,zrs,Δzb0,andlb. These values are adjusted in order to obtain an axial behavior similar to the actual axial behavior.The valuezn0=0 is considered in simulations as the reference for axial position. For the first simulation for adjustment, it was considered that (i) the distance zrb is equal to the distance between bucket centerline and the axial position of bearing center Zrb determined by turbine specifications; (ii) the distance zrs is equal to the distance between bucket centerline and the axial position of generator rotor geometric center Zrs determined by turbine specifications; and (iii) Δzb0 is considered as zero which means that the axial clearance is supposed to be equal or lower than the one necessary to obtain an adequate lubricant film behavior. Taking into account these considerations for the first simulation, then the values of zs0,zb0,andlb are determined. If these values cannot provide a good agreement of model to experimental data, then the values of zrb,zrs, and Δzb0 are modified and adjusted. After this, the modulation functions for hydraulic and magnetic forces are modified until the best agreement with measurements is obtained. This procedure is schematically shown in Figure 17.Figure 17
Adjustment of model parameters.The adjustment was carried out making use of different measurements of turbine. It is useful to count on the axial motion measured rotating at constant speed, zero current, and zero voltage in generator: full speed no load (FSNL). This condition allows focusing onzb0 and lb because there is not magnetic force on rotor at this operating condition. Then, it is necessary to count on the axial motion measured at different power in order to focus on zs0.


## 3. Actual Cases

### 3.1. Measurements

Measurements were carried out in two different power stations. In both stations, power rate, stator and rotor currents and voltages, and nozzle spear position were obtained from power station control system. Axial motion was measured making use of displacement sensor Baumer IWRM 18U9501 with a measuring range from 2 to 5 mm mounted in a support with a magnetic ferrite base (Figure18). The displacement sensor was powered by a 24 V power supply. The rotating speed was determined using the voltage output of a Monarch PLT200 optical tachometer which detected a reflective tape in rotor using a laser light source. The voltage signal from displacement sensor and optical tachometer was recorded using a NI 9229 module installed in a CompactDAQ 9174 USB chassis both from National Instruments. The analog to digital converter of NI 9229 allowed measuring a voltage range of ±60 V at 24 bits. The sampling frequency used in measurements was 25600 Hz.Figure 18
Axial movement measurement.


### 3.2. Case I: Actual Case with Normal Axial Motion

The first actual case, using the model, considers the Pelton turbine described in Table1. This unit was operating for twenty-two months presenting no bearing failure and normal working temperatures. After parameter adjustment, the values of axial position of stator and rotor magnetic field centers, bearing center, bearing axial clearance, and bearing nonworking interval are shown in Table 2. Figure 19 shows the waveform of model prediction and measurements at different operating conditions. Table 3 shows the maximum, minimum, and averaged axial position obtained with model predictions and measurements. A good agreement is observed between simulated and measured axial motion.Table 1
Horizontal Pelton turbine, Case I.
Power rate10MWNet head600mFlow1.89m3/sNumber of buckets19Rotating speed1000rpmRotor length1.240mRotor diameter1.480mStator diameter1.500mTotal rotor weight11800kgNumber of poles6Table 2
Calculated parameters, Case I.
zrs−zs0−zn01.70mmzrb−zb0−zn0−0.30mmΔzb00.00mmlb0.90mmMáxahn,bhn0.00Máxamn,bmn0.00Figure 19
Simulated and measured waveform, Case I. (a) Simulation. (b) Experimental.
(a)(b)Table 3
Model prediction and experimental measurements, Case I.
Power rate in MWAxial position of rotor in mmModel predictionExperimental measurementDifference0Maximum0.450.450.00Minimum−0.45−0.450.002Maximum−0.27−0.280.01Minimum−0.29−0.280.014Maximum−0.25−0.260.01Minimum−0.26−0.260.006Maximum−0.22−0.230.01Minimum−0.23−0.230.008Maximum−0.19−0.200.01Minimum−0.20−0.200.0010Maximum−0.17−0.170.00Minimum−0.18−0.170.01After this analysis, and taking into account the values in Table2, it was concluded that the axial motion was normal and adequate for continuous operation with no restrictions.


### 3.3. Case II: Actual Case with Excessive Axial Motion and Bearing Failure

The second actual case, using the model, considers the Pelton turbine described in Table4. The stator poles of this Pelton turbine were replaced due to the end of its useful life. After startup, there was a bearing failure, showing Babbitt erosion and contact in axial direction. After this failure, bearing was repaired and mounted. The measurements shown in this section were taken after reparation. After parameter adjustment, the values of axial position of stator and rotor magnetic field centers, bearing center, bearing axial clearance, and bearing nonworking interval are shown in Table 5. Figure 20 shows the waveform of model prediction and measurements at different operating conditions. Table 6 shows the maximum, minimum, and averaged axial position obtained with model predictions and measurements. A good agreement is observed between simulated and measured axial motion.Table 4
Horizontal Pelton turbine, Case II.
Power rate8MWNet head475mFlow1.91m3/sNumber of buckets19Rotating speed300rpmRotor length0.500mRotor diameter3.330mStator diameter3.460mTotal rotor weight20800kgNumber of poles20Table 5
Calculated parameters, Case II.
zrs−zs0−zn09.65mmzrb−zb0−zn010.28mmΔzb00.00mmlb0.80mmMáxahn,bhn0.00Máxamn,bmn0.00Figure 20
Simulated and measured waveform, Case II. (a) Simulation. (b) Experimental.
(a)(b)Table 6
Model prediction and experimental measurements, Case II.
Power rate in MWAxial position of rotor in mmModel predictionExperimental measurementDifference2Maximum0.040.030.01Minimum0.030.030.004Maximum0.180.180.00Minimum0.170.180.016Maximum0.300.300.00Minimum−0.50−0.500.008Maximum0.300.300.00Minimum−0.50−0.500.00In this case, there was a nonstationary behavior for power rates over 6 MW. Over 6 MW, there is an axial bounce of rotor against bearing. After this analysis, and taking into account the values in Table5, it was concluded that the axial motion was mainly due to the difference between axial position of stator and rotor magnetic field centers and due to the difference between axial positions of bearing. Until corrective action, it was suggested not to operate over 4 MW. A corrective action was performed moving bearing in order to correct axial position. It was not considered to move stator frame or to move nozzle because it is more expensive in time and costs when comparing to bearing adjustment. After this correction, no bounces in axial motion were detected at 6 and 8 MW, and the unit was able to operate at all its power ranges.


## 3.1. Measurements

Measurements were carried out in two different power stations. In both stations, power rate, stator and rotor currents and voltages, and nozzle spear position were obtained from power station control system. Axial motion was measured making use of displacement sensor Baumer IWRM 18U9501 with a measuring range from 2 to 5 mm mounted in a support with a magnetic ferrite base (Figure18). The displacement sensor was powered by a 24 V power supply. The rotating speed was determined using the voltage output of a Monarch PLT200 optical tachometer which detected a reflective tape in rotor using a laser light source. The voltage signal from displacement sensor and optical tachometer was recorded using a NI 9229 module installed in a CompactDAQ 9174 USB chassis both from National Instruments. The analog to digital converter of NI 9229 allowed measuring a voltage range of ±60 V at 24 bits. The sampling frequency used in measurements was 25600 Hz.Figure 18
Axial movement measurement.


## 3.2. Case I: Actual Case with Normal Axial Motion

The first actual case, using the model, considers the Pelton turbine described in Table1. This unit was operating for twenty-two months presenting no bearing failure and normal working temperatures. After parameter adjustment, the values of axial position of stator and rotor magnetic field centers, bearing center, bearing axial clearance, and bearing nonworking interval are shown in Table 2. Figure 19 shows the waveform of model prediction and measurements at different operating conditions. Table 3 shows the maximum, minimum, and averaged axial position obtained with model predictions and measurements. A good agreement is observed between simulated and measured axial motion.Table 1
Horizontal Pelton turbine, Case I.
Power rate10MWNet head600mFlow1.89m3/sNumber of buckets19Rotating speed1000rpmRotor length1.240mRotor diameter1.480mStator diameter1.500mTotal rotor weight11800kgNumber of poles6Table 2
Calculated parameters, Case I.
zrs−zs0−zn01.70mmzrb−zb0−zn0−0.30mmΔzb00.00mmlb0.90mmMáxahn,bhn0.00Máxamn,bmn0.00Figure 19
Simulated and measured waveform, Case I. (a) Simulation. (b) Experimental.
(a)(b)Table 3
Model prediction and experimental measurements, Case I.
Power rate in MWAxial position of rotor in mmModel predictionExperimental measurementDifference0Maximum0.450.450.00Minimum−0.45−0.450.002Maximum−0.27−0.280.01Minimum−0.29−0.280.014Maximum−0.25−0.260.01Minimum−0.26−0.260.006Maximum−0.22−0.230.01Minimum−0.23−0.230.008Maximum−0.19−0.200.01Minimum−0.20−0.200.0010Maximum−0.17−0.170.00Minimum−0.18−0.170.01After this analysis, and taking into account the values in Table2, it was concluded that the axial motion was normal and adequate for continuous operation with no restrictions.


## 3.3. Case II: Actual Case with Excessive Axial Motion and Bearing Failure

The second actual case, using the model, considers the Pelton turbine described in Table4. The stator poles of this Pelton turbine were replaced due to the end of its useful life. After startup, there was a bearing failure, showing Babbitt erosion and contact in axial direction. After this failure, bearing was repaired and mounted. The measurements shown in this section were taken after reparation. After parameter adjustment, the values of axial position of stator and rotor magnetic field centers, bearing center, bearing axial clearance, and bearing nonworking interval are shown in Table 5. Figure 20 shows the waveform of model prediction and measurements at different operating conditions. Table 6 shows the maximum, minimum, and averaged axial position obtained with model predictions and measurements. A good agreement is observed between simulated and measured axial motion.Table 4
Horizontal Pelton turbine, Case II.
Power rate8MWNet head475mFlow1.91m3/sNumber of buckets19Rotating speed300rpmRotor length0.500mRotor diameter3.330mStator diameter3.460mTotal rotor weight20800kgNumber of poles20Table 5
Calculated parameters, Case II.
zrs−zs0−zn09.65mmzrb−zb0−zn010.28mmΔzb00.00mmlb0.80mmMáxahn,bhn0.00Máxamn,bmn0.00Figure 20
Simulated and measured waveform, Case II. (a) Simulation. (b) Experimental.
(a)(b)Table 6
Model prediction and experimental measurements, Case II.
Power rate in MWAxial position of rotor in mmModel predictionExperimental measurementDifference2Maximum0.040.030.01Minimum0.030.030.004Maximum0.180.180.00Minimum0.170.180.016Maximum0.300.300.00Minimum−0.50−0.500.008Maximum0.300.300.00Minimum−0.50−0.500.00In this case, there was a nonstationary behavior for power rates over 6 MW. Over 6 MW, there is an axial bounce of rotor against bearing. After this analysis, and taking into account the values in Table5, it was concluded that the axial motion was mainly due to the difference between axial position of stator and rotor magnetic field centers and due to the difference between axial positions of bearing. Until corrective action, it was suggested not to operate over 4 MW. A corrective action was performed moving bearing in order to correct axial position. It was not considered to move stator frame or to move nozzle because it is more expensive in time and costs when comparing to bearing adjustment. After this correction, no bounces in axial motion were detected at 6 and 8 MW, and the unit was able to operate at all its power ranges.


## 4. Conclusions

A model of axial motion of Pelton horizontal turbine is proposed. The model needs a set of experimental data under different operating conditions in order to determine the unknown values of axial position of stator and rotor magnetic field centers, bearing position, and bearing clearance. After the adjustment of these variables, the model accurately predicts axial motion of Pelton horizontal turbine. Knowing these variables allows evaluating the influence of different sources of axial motion: due to the difference in centerlines of nozzle and bucket and due to the difference in axial position of stator and rotor magnetic field centers; due to bearing positioning and due to bearing axial clearance.The model is tested in two cases: with a turbine under normal axial motion and with a turbine presenting excessive axial motion and damage in one of the bearings. In the case of excessive axial motion, it was possible to determine the main causes of axial motion. Taking this into account, a corrective action was suggested. After corrective action, the excessive axial motion was reduced and no longer axial failure in bearing occurred.



---

*Source: 8811961-2020-11-25.xml*