# Path Transmissibility Analysis Considering Two Types of Correlations in Hydropower Stations

Authors: Baoping Zhi; Zhenyue Ma
Journal: Journal of Applied Mathematics
 (2013)

---

## Abstract

A new vibration model is built by introducing the head-cover vibration transfer path based on a previous analysis of the vertical vibration model for hydropower station units and powerhouses. This research focuses on disturbance- and parameter-related transfer paths in a practical situation. In a complex situation, the application of the stochastic perturbation method is expanded using an algebra synthesis method  the Hadamard product, and theoretical analyses, and numerical simulations of transfer paths in the new vibration model are carried out through the expanded perturbation method. The path transfer force, the path transmissibility, and the path disturbance ranges in the frequency domain are provided. The results indicate that the methods proposed in this study can efficiently reduce the disturbance range and can accurately analyze the transfer paths of hydraulic-source vertical vibration in hydropower stations.


---

## Body

## 1. Introduction

Hydraulic vibration is the main vibration source in hydropower station units and powerhouses. The vertical vibration in units is usually caused by hydraulic pressure fluctuations or other loads on the water turbine’s flow passage components. Field and model tests have shown that there are three main vibration transfer paths running from the water turbine to the powerhouse [1]: (1) runner-shaft-bearing-fixed components (machine frame, head-cover-powerhouse;) (2) flow pressure-spiral case-powerhouse; and (3) runner-runner negative pressure region-head cover-powerhouse. Previous studies concerning the vertical vibration produced by a hydraulic source mainly focus on path (1), while the effects of paths (2) and (3) are usually ignored [2]. However, as the scale and capacity of hydropower stations increase, the flow passage area of the head-cover system continually increases, with corresponding increases in the vibration of the head-cover system. Therefore, the influence of the head-cover system becomes more important in the hydraulic vibration transfer path, and ignoring path (3) will produce a larger error. Therefore, it is imperative that the contribution of path (3) should be analyzed; more specifically, the contribution of the vibration of the head-cover system of the hydropower station vertical vibration transfer should be analyzed. However, due to the presence of randomness, it is difficult to clearly and accurately describe the contribution of this transfer path to the structural vibration.Theoretical analyses of the vibration transfer path can be classified as a stochastic structural system problem. At present, the Monte-Carlo numerical simulation method (MCSM) [3, 4] and the perturbation method [5–7] are the most popular analysis methods. MCSM is used less frequently because of the large amount of computation required when dealing with a large-scale structure. Conversely, the perturbation method is applied by many researchers in various fields. Collins and Thompson [8] initially employed the perturbation method to analyze stochastic dynamical systematic characteristics in 1969; the perturbation method was later employed in a static analysis by Hisada and Nakagiri [9] and in a dynamic analysis by Liu et al. [10]. Kronecker algebra was introduced to the expansion of the perturbation method by Vetter [11]. After several decades of development, relevant studies on perturbation theory were quite abundant. The primary methods included the L-P method [12], the multiple scale method [13, 14], the average method [15], the KBM (Krylov-Bogoliubov-Mitropolsky) method [15, 16] and the singular perturbation method [17, 18]. Recent developments include the homotopy perturbation method [19]. In the analysis of transfer paths combined with the perturbation method, Zhang et al. built a theoretical model for vibration transfer path analysis by implementing the perturbation method with Kronecker algebra in the static analysis [20], the dynamic analysis [21] and the reliability analysis [22, 23]. This model has been used for the analysis and design of mechanical components. Computations of the transfer path considering multi-vibration sources have been carried out by Zhao and Zhang [24] based on Zhang’s research. These results have been used for isolation vibration analysis. The above research studies focus on additive disturbance analysis. Using multiplicative disturbance analysis, Gao et al. separately analyzed the dynamical characteristics of a truss structure by using the interval factor method [25], the random factor method [26], and nonstationary random excitation [27]. Ma et al. conducted a dynamical characteristic analysis of a linear [28] and nonlinear [29] truss structure using fuzzy variables. In 2010, a two-factor method was proposed by Ma et al. [30], which considers two multiplicative disturbances. However, because of the limitations of these methods, they have not been applied to transfer path analysis with respect to a multiplicative disturbance. Furthermore, previous studies have not jointly analyzed these two types of disturbances. In practice, the two types of disturbances always exist together, and there are some correlations between these two disturbances. There are some correlations among the parameters of the structure as well. Thus, the previous studies were not based on an accurate analysis because they only considered a single disturbance, neglected the correlations between disturbances, or neglected the correlations among parameters. Therefore, it is necessary to jointly analyze the problem of the vibration transfer path of time-invariable parameters for the two types of disturbances and their related parameters. Only a few studies have reported on the two types of correlations in perturbation theory. Previous studies on correlations ignored the two types of correlation or converted relevant variables into irrelevant variables. These analyses did not consider the effects of correlations and thus were not rigorous. At present, only a few researchers have attempted to explore this aspect. Some examples include studies on molecular chemistry [31], nuclear physics [32], and vocal vibration [33, 34]. Of these studies, only the study on vocal vibration considered the correlation parameters, but it did not involve the perturbation method. At present, there are few available studies on structural analysis. Pirrotta [35] conducted a perturbation analysis of delta-correlated processes, but he did not study the correlation between parameters; Ambrogio [36] described the correlations by considering the correlation coefficient between additive and multiplicative disturbances for one parameter, but he also did not refer to the correlation between parameters. Husain et al. [37] solved the problem of parameter correlation, but they only considered computational results and did not extend their results to improve the perturbation method. Khodaparast et al. [38] corrected the structure variation model, but they ignored the correlation between modified parameters and measurement values. This study focuses on the application of perturbation theory with two types of correlation from the perspective of structure analysis.Hydropower station units and powerhouses are large-scale structures. Most of the parameters’ disturbances should be less than 10% of their mean values in such a large structure. Therefore, problems regarding hydropower stations can be solved by the perturbation method. First, the head-cover system is introduced, which has the same vibration source as a path(1); the elastic foundation constraint is selected; and a new vibration model is built on the basis of the previous vertical vibration model. Furthermore, a method for solving the path transfer force is proposed using the general method of dynamic analysis. Second, for the test signal, correlations between the two types of disturbances and correlations between the parameters are considered based on the single disturbance vibration path analysis. By using the coefficient of variation algebra synthesis method [39], Kronecker algebra [11], and the Hadamard product [40], gradient-sorting estimations of the transfer paths in the frequency domain are carried out, and methods for determining the transfer force and transmissibility and their disturbances are proposed. Finally, the method described in this paper is verified using the model of a large hydropower station.


## 2. Analysis Model with the Introduction of a Head-Cover System

For a hydroelectric generating unit, regardless of whether it has a suspension or umbrella structure, the weight of its rotating parts is successively transferred to the reinforced concrete machine foundation through the thrust bearing, frame (suspension units containing the stator frame), and sole screw. The head-cover system is always fixed on the base ring strengthening plate, and the head-cover system and the strengthening plate are considered to be one part. Taking a vertical vibration characteristic analysis of the umbrella unit as an example, the model contains a shaft system, thrust bearing, and lower bracket (see Figure1). The heave shaft can be simplified as a massless elastic continuous beam, and then its mass can be regarded as a mass attached to three nodes, m1, m2, and m3.  m1 can be defined as the mass of the excitation rotor and the shaft, the half shafting mass, which is measured from the heavy shaft top to the rotor frame, and another mass added on top of the heavy shaft;  m2 can be defined as the mass of the central body of the rotor frame, the half mass of the whole gate arm, and the half mass of the whole shaft; and m3 can be defined as the mass of the water, turbine runner, the additional mass of water and the half shafting mass, which is measured from the rotor frame to the hydraulic turbine. The rotor gate arm can be simplified as a massless elastic continuous rod, and then its mass can be assigned to the runner margin and the central body of the rotor frame. k4 can be defined as the sum of the vertical stiffness of the whole gate arm, and m4 can be defined as the lumped mass of the runner margin. The outer end of the lower bracket is fixed to the concrete foundation. The lower bracket gate arm can then be simplified as a gravity-free beam if the coupling effect of the foundation is ignored. k52 can be defined as the vertical stiffness of the lower bracket gate arm. m5 can be defined as the lumped mass of one end of the lower bracket, which is close to the heavy shaft, and the half mass of the gate arm; m5 can be connected to m2 by the thrust bearing, which is simplified by an equivalent stiffness k51.Figure 1
Simplified model of the coupling system between the umbrella unit and the powerhouse.With the introduction of the head-cover system, the vibration transfer path is as follows: first, the hydraulic vertical vibration is transferred to the head cover by the heavy shaft seal and guide bearing in the water turbine runner chamber; next, the vibration is transferred from the head cover to the spiral case base ring strengthening plate, which is connected to the outer end of the head cover; and finally, the vibration is transferred to the machine foundation by the wrapped concrete outside the spiral case. Ignoring the coupling effect, the control parts and other additional parts on the head cover can be regarded as the attached mass of the head-cover system. As the lumped mass is close to the heavy shaft,m6 can be defined as the mass of the central body and the half mass of the whole head-cover system. The head-cover system can be simplified as a gravity-free beam, and k62 can be defined as the vertical stiffness. m6 can be connected to m3 by the sealing spring, which lies between the head-cover structure and the water turbine runner. The vertical stiffness of the connection can be simplified by an equivalent stiffness k61.k
51 is the series stiffness of the elastic oil tank stiffness (thrust bearing support system) and the oil film stiffness. k61 is the series stiffness of the sealing structure stiffness and the clearance water stiffness. These two parameters exhibit a linear relationship because the unit’s axial water thrust varies with the unit’s conditions. However, k51 and k61 are simplified as a single stochastic variable in this study.In the process of examining the unit’s vertical vibration, the machine foundation pier can be regarded as an elastic foundation. The structure can be treated as a single node,m7, which ranges from the machine foundation pier to the concrete floor in the turbine layer. The structure can be treated as a single node, m8, which ranges from the turbine layer to the foundation. The stiffness of this structure is replaced by the equivalent stiffness, which is based on the strengthening plate, spiral structure, and concrete structure. Thus, the hydropower station base is treated by dividing it into two nodes at the connection of the head-cover system and the hydropower station base, where k62 is only connected to m8.For the entire model, the form and meaning of the dampness matrix are similar to those of the stiffness matrix. Assuming that the system is linear, the differential equation of the vibration is found from the Lagrange equation:(1)Mu¨+Cu˙+Ku=F(t).By merging these dynamical balance equations for the shaft, the rotor, the lower bracket, the head-cover system, and the machine foundation pier, an equation with 8 degrees of freedom and 26 parameters can be obtained, and the total stiffness matrix can be written as(2)K=(k1-k1000000-k1k1+k3+k4+k51-k3-k4-k510000-k3k3+k6100-k61000-k40k400000-k5100k51+k520-k52000-k6100k61+k620-k620000-k520k52+k7-k700000-k62-k7k62+k7+k8).The total mass matrix is obtained from the lumped mass:(3)M=diag{m1,m2,m3,m4,m5,m6,m7,m8}.Generally, the characteristics of the vertical vibration source in the hydropower station are unique position, simple contact surface with structure, and few ingredients. Conversely, the characteristics of the lateral (radial) vibration are multipositions, strong nonlinear contact surface, and complex excitation. For simplifying the model and the calculation, the total structure of the unit’s powerhouse is considered as suffering vertical harmonic excitation. The vertical harmonic excitation is located on the water turbine runner, and the response of the total structure is in a steady state. In addition, the initial phase remains constant throughout the whole process. Setting the steady-state response asui(t)=U-iei(ωt+φi) and Ui=U-ieiφi yields(4a)ui(t)=Uieiωt,u˙i(t)=iωUieiωt,u¨i(t)=-ω2Uieiωt,(4b)U={u1,u2,u3,u4,u5,u6,u7,u8}T,(4c)F(t)={0,0,F0eiωt,0,0,0,0,0}T.Substituting (4a)–(4c) into (1) gives the dynamical balance equation of the hydropower station’s vertical vibration with the head-cover system in the frequency domain:
(5)(-ω2M+iωC+K)U=F(t).The response vectorU at each node of the structure was obtained by numerical calculations. Fzhou and Fding denote the forces that are transferred to the machine foundation pier from the water turbine by the shaft system and the head-cover system, respectively, and can be written as(6a)Fzhou=kzhou(u3-u7)+czhou(u˙3-u˙7)=(kzhou+iωczhou)(U3-U7)eiωt,(6b)Fding=kding(u3-u7)+cding(u˙3-u˙7)=(kding+iωcding)(U3-U7)eiωt.In (6a)–(6b), kzhou and czhou denote the path stiffness and the path dampness of the vibration path through the shaft system, respectively. Similarly, kding and cding denote the path stiffness and the path dampness through the head-cover system, respectively. When calculating the path stiffness and dampness of the shaft system, m1 and m4 are treated as dynamic vibration absorbers with dampness.


## 3. Vibration Path Sorting Considering Two Types of Correlations

In the analysis of a hydropower station, parameters such as the stiffness, mass, and dampness are attributed to the multiplicative disturbances arising from the material properties, manufacturing technology, and other factors. These parameters, obtained by measurement, are attributed to the additive disturbances due to the testing noise and environmental noise. Furthermore, because the testing method and environment are the same for every parameter, the two disturbances of each parameter are interrelated. These parameters are also interrelated by means of continuous structure discretization.Given the factors mentioned above, the mass, stiffness, and dampness are described as a random vectora with n random variables, where every random variable involves two types of disturbances. When the variable disturbance is lower than 15% of the mean value, the stochastic variable can be expressed as
(7)ai=ai1aid+ai2.In (7), ai denotes the ith element in the random vector a. ai1 and ai2 denote the multiplicative and additive disturbance of the random variable ai, respectively. Assuming that their mean values are 1 and 0, respectively, aid denotes the deterministic component of ai, which represents the mean value after multiple samplings.According to the coefficient of variation algebra synthesis method [40], calculating the mathematical expectation and variance of (7) yields
(8)Eai=E[ai1aid+ai2]=E(ai1aid)+E(ai2)=aid,σai2=Var(ai)=E[(ai-Eai)2]=E[((ai1-ai)aid+ai2)2]=(aid)2σai12+σai22+2aidCov(ai1,ai2),
where
(9)Cov(ai1,ai2)=ρai1,ai2σai1σai2,
where ρa1,a2 is the correlation coefficient of the multiplicative and additive disturbances contained in the parameter a. If these parameters are interrelated, then
(10)Cov(ai,aj)=ρai,ajσaiσaj,
where ρai,aj is the correlation coefficient of the random variables ai and aj. Generally, if the random variables follow a normal distribution, their linear transformations and multiplication will also follow a normal distribution. For example, if the vector a follows a normal distribution, then the function Fi(a) will follow a normal distribution. Using the Taylor expansion, expanding the transfer force Fi at the mean value Fid yields
(11)Fi=Fid+∂Fid∂aT(a-ad)+O(ap)=Fid+∂Fid∂aTap+O(ap),
where
(12)ap=a-ad=(a1-1)ad+a2.Neglecting components of the second order and above, the disturbance of the transfer force is(13)Fip=∑k=1m∂Fid∂akakp.In (13), ∂Fid/∂ak is the partial derivative of Fi with respect to the random variable ak, which is also the first-order sensitivity of ak. The covariance of Fip and Fjp is
(14)Cov(Fi,Fj)=E[FipFjp]=E[(∑k=1m∂Fid∂akakp)(∑l=1m∂Fjd∂alalp)].That is,(15)Cov(Fi,Fj)=∑k=1m∑l=1m∂Fid∂ak∂Fjd∂alE(akpalp)=∑k=1m∑l=1m∂Fid∂ak∂Fjd∂alCov(ak,al).Equation (16) shows that Cov(Fi,Fj) of the transfer force Fi and Fj can be expressed by Cov(ak,al) of the random structure parameter:
(16)Cov(Fi,Fj)=∑i=1m∑j=1m∂Fid∂ak∂Fjd∂alρak,alσakσal.Ifi=j,
(17)σFi2=Cov(Fi,Fi)=∑k=1m∑l=1m∂Fid∂ak∂Fid∂alρak,alσakσal.According to Kronecker algebra [11], the corresponding stochastic analysis theory and Hadamard product [41] yield
(18)σFi2=|∂Fid∂aT|[2](ρ⃑∘(σa)[2]),
where the subscript [2] denotes the Kronecker power; that is, a[2]=a⊗a; if a is of order n×1, then a[2] is of order n2×1. ρ⃑ is the stacking vector of the matrix of correlation coefficients, which is also of order n2×1. The symbol ∘ denotes the Hadamard product. Therefore, σFi2 is of order 1×1.The sensitivity matrix|∂Fid(a)/∂aT| of the transfer force of each parameter is
(19)|∂Fid∂aT|=|∂Fid∂a1⋯∂Fid∂an|.Substituting (19) into (18), the variance of the transfer force Fi for each path can be calculated. Equation (11) only involves the first-order Taylor expansion. A higher-order Taylor expansion will improve the accuracy, but it involves complicated mathematical calculations. Equation (18) shows that the transfer force variance can be directly obtained from the random variables’ numerical characteristics. The calculations can thus be simplified because the sample is not included in the mathematical operations. The correlations between parameters are merely modified in (18) and do not increase the number of calculations required. The transmissibility is defined as the ratio of the amplitudes between the transfer force and the vibration source excitation force:
(20)βi=|FiF0|.Neglecting the disturbance of the excitation forceF0, based on the random variable algebra synthesis method, the expectation, variance, and transfer coefficient of the transmissibility can be written as(21a)Eβ=E[β(a)]=E(Ft)E(F0),(21b)σβ2=E[Fi]2E[F0]2|Var[Fi][E(Fi)]2|,(21c)θ=Eβσβ.The result of a random variable that follows a normal distribution divided by a constant also follows a normal distribution.E[F0] and E[Fi] are known in the deterministic system, so the variance of the transmissibility will follow a normal distribution. The transfer coefficient is the transfer efficiency of the excitation force. The gradient sorting of the vibration path transmissibility in the frequency domain can be obtained by changing the frequency of the excitation force.A solution of the vibration transfer path for the two types of disturbances and correlations is proposed based on the above methods. The path transmissibility and its probabilistic characteristic are provided. In this study, the solution only involves the first-order sensitivity of the random parameter and the probabilistic characteristic of the random variable. Additionally, the introduction of the two types of correlations does not involve excessive calculations. Therefore, the calculation accuracy is improved for practical problems.


## 4. Example Analysis

The main structure of the umbrella unit in a large hydropower station is shown in Figure2, and the simplified model of transfer path is shown in Figure 1. In this example, the effect of the spiral case and substructure is ignored [41, 42] because their effect is far less than that of the upper structure. The excitation is assumed to be a simple harmonic excitation. The mean values of the random parameters can be obtained from the hydropower station design diagrams. These mean values are m1=8.28×104, m2=1.042×106, m3=3.29×105, m4=9×105, m5=1.2×105, m6=1.15×105, m7=1.39×105 and m8=8.92×105, in units of kg; k1=7.26×1010, k3=5.72×1010, k4=2.32×1010, k51=2.20×1012, k52=9.41×109, k61=1.73×108, k62=1.73×1010, k7=7.70×109, and k8=4.26×108, where the unit of the stiffness k is N/m; and c1=5.48×106, c3=4.11×106, c4=1.02×107, c51=2.57×107, c52=7.51×105, c61=2.23×105, c62=9.99×104, c7=1.64×106, and c8=9.74×105, where the unit of the dampness c is N·s/m. Each parameter comprises two disturbances that follow a normal distribution. The multiplicative disturbance is related to the difficulty of obtaining the parameters. The multiplicative variance coefficients of k51 (including the vertical stiffness of the thrust bearing) and k61 (the sealing equivalent vertical stiffness between the head cover and runner) are set to 0.10. The equivalent bending rigidity of k62s multiplicative variance coefficient is set to 0.075 because there are many uncertain factors in the control components and other attached components on the head cover. The multiplicative random variance coefficients of the other parameters are set at 0.05. The additive disturbance is related to the measurement range. The standard deviations of the additive random variable with respect to mass and stiffness are set at 104 and 108, respectively, according to each parameter’s mean value. The variance coefficient of damping is similar to that of the mass and stiffness. In this example, the parameters’ disturbances are generated by a function in the MATLAB software, and the sample size is 10000 for each disturbance. All errors are determined to be smaller than 0.1% by comparing the variances and mean values of the generated samples with the corresponding set values. Therefore, the generated sample variances and mean values are used in the study. The correlation function is determined by the generated function. The partial derivative is calculated by the software Mathematica. The coefficient correlation of the generated sample is 0.0329.Figure 2
Cutaway view of the main structure of hydropower.Figure3 displays the transmissibility variances of the head-cover transfer path obtained using different methods. Method 1 does not involve the two types of correlations. The Monte-Carlo simulation is not used in method 1 (due to the inclusion of 26 random variables and because the calculation time is too long, this method only focuses on how each group of variables changes with time, that is, the number of calculations is equal to the sample size). Method 2 (following the method developed by W. J. Vetter in 1973) only involves the correlation with respect to one parameter, which involves the correlations between the disturbances of the parameters. Method 3 (similar to method 2) involves the correlations between parameters instead of the correlations between the disturbances of the parameters. Method 4 involves both types of correlations. The results for method 4 are obtained from (18). Figure 3 displays the characteristic curve of the head-cover system transmissibility variances in the frequency band of 0–1000 rad/s. Figure 4 presents the optimal curve of the two types of correlations for the transmissibility variance range as the correlation coefficient changes. This figure contains two curves: the disturbance curve, which describes the mean values of the results of method 2 subtracted from method 1 and the results of method 4 subtracted from method 3, and the parameter curve, which is similar to the disturbance curve.Figure 3
Head-cover transmissibility variances in the frequency domain as predicted by the four methods.Figure 4
Optimization range of the correlation coefficient characteristic curve.Figures3 and 4 show the following results: (1) compared to the other methods, the variance range for the method that contains two types of correlations is the smallest. This finding indicates that this method significantly decreases the disturbance range under theoretical calculations. The path transmissibility variance decreases by 72.51%. (2) In optimizing the method considering correlations between the disturbances and variances, the range merely decreases by 0.644%. This small decrease may have arisen because the difference between the additive and multiplicative disturbances is large, and the optimizing function is weak. (3) The optimization effect is significant for the method considering the correlations between parameters. The variance range decreases by 75.71% when the correlation coefficient is 0.0085. Thus, neglecting this type of correlation leads to an inaccurate analysis. (4) As the correlation coefficient increases, the optimization of the transmissibility variance is more efficient, and the effect on the disturbance range is more significant. This result demonstrates that the relationship between the optimization range and the correlation coefficients for the parameters is linear and that the relationship between the optimization range and the correlation coefficients for the disturbances is proportional. However, the latter relationship does exhibit some fluctuations.Figure5 compares the characteristic curve for each path transmissibility, βi, to the excitation frequency. This curve shows the ratio between the forces that are transmitted to the machine pier by the corresponding paths and the excitation force of the water turbine hydraulic source. This ratio is equivalent to the amplification coefficient. The curve indicates two results: (1) resonance occurs when the frequency of the excitation force is the same as the natural frequency. The transmissibility reaches its maximum value for each path at this frequency. This result indicates that the path has a large effect on the structural vibration at this frequency. (2) Considering the whole process of vertical vibration transfer, the transmissibility of the head cover is significantly smaller than that of the shafting system. Throughout the entire frequency range, the mean value of the path transmissibility ratio is 32.37. This value indicates that the vibration effect of the head-cover system can be neglected in accurate calculations.Figure 5
Path transmissibility of the excitation frequency curve.Figure6 shows that the characteristic curve of the transmissibility variance varies with the excitation frequency for the two types of correlations. This curve reflects the diversion of the transmissibility and indicates the following findings. (1) The path transmissibility variance of the shaft system is greater than that of the head-cover system due to the large number of components and the randomness of these components’ parameters. Thus, the variance of a structure with a high level of parameter randomness is large in the transfer path. (2) The variance maximum occurs at the natural frequency. This frequency is determined by the structural parameter, which accounts for the parameter contribution of every transfer path, enabling the transfer force to be obtained from the expectation and variance of the transmissibility.Figure 6
Variance of the path transmissibility of the excitation frequency characteristic curve.Figure7 shows that the characteristic curve of the transfer coefficients varies with the excitation frequency for the two types of correlations. This curve reflects the transfer efficiency of the excitation force. The figure shows that the transfer efficiency of the head-cover system is more efficient than that of the shafting system. Therefore, more attention should be given to the assigned vibration proportion of the shaft and head-cover system in a design.Figure 7
Transfer coefficient of the paths of the excitation frequency characteristic curve.


## 5. Conclusions

(1)
The multivibration source, multipath vibration model, is improved with respect to the coupling effect between hydropower station units and powerhouses by the introduction of a head-cover system. The scalar expression of each vibration transfer path is provided in this study.
(2)
A perturbation analysis of the transfer path is carried out by using a test signal. A complex situation with two types of correlations is considered. Kronecker algebra, the Hadamard product, and probabilistic statistics are employed to develop an analysis method for the transfer path. In the solution process, excessive additional computations are not required because only the numerical characteristics of the random variables are used. The application of the stochastic perturbation method is expanded, and a method for analyzing the vibration transfer path is developed. The analysis method efficiently decreases the disturbance range of the path’s contribution. The disturbance ranges of the path transmissibility and the contribution rate are efficiently reduced by considering the two types of correlations. This reduction is significant for the optimization function with respect to the results.
(3)
The simulation results indicate that when the effect of the spiral case and substructure is ignored, the influence on the disturbance range continually increases as the correlation coefficient increases in the hydropower station model. The optimization function of the correlations between parameters is too important to be neglected when calculating the disturbance range.
(4)
In conclusion, the analysis of the proposed model shows that the effect of the water turbine head cover is not evident in the vertical vibration transfer, but its transfer efficiency is significant.The analysis of vibration transfer paths of hydropower station units and powerhouses is complicated. Based on the complex disturbances and parameters, the contribution of each transfer path can be calculated by analyzing the sensitivity of the transfer force and the transmissibility of the vibration model in the frequency domain. This study provides a reference for future comprehensive research on the transfer paths of hydropower station units and powerhouses.



---

*Source: 802546-2013-09-08.xml*