Research Article  Open Access
Fractal Theory and Contact Dynamics Modeling Vibration Characteristics of Damping Blade
Abstract
The contact surface structure of dry friction damper is complicate, irregular, and selfsimilar. In this paper, contact surface structure is described with the fractal theory and damping blade is simplified as 2DOF cantilever beam model with lumped masses. By changing the position of the damper, lacing and shroud structure are separately simulated to study vibration absorption effect of damping blade. The results show that both shroud structure and lacing could not only dissipate energy but also change stiffness of blade. Under the same condition of normal pressure and contact surface, the damping effect of lacing is stronger than that of shroud structure. Meanwhile, the effect on changing blade stiffness of shroud structure is stronger than that of lacing. This paper proposed that there is at least one position of the blade, at which the damper dissipates the most vibration energy during a vibration cycle.
1. Introduction
Lacing and shroud structure are widely adopted in turbine blades to effectively restrain blade resonance and have features such as simple structure and easy installation. Lacing and shroud structure are based on the mechanism of dry friction between adjacent blades to achieve lower blade vibration. The blades are pretwistedly installed so that the lacing and shroud structure among adjacent blades come into contact and generate some positive pressure. In operation, blades twisted under the centrifugal force and the positive pressure caused by torsion provide adequate positive pressure for adjacent blades.
While forced by the exciting force, blades vibrate and the relative displacements among adjacent blades are formed. The dry friction structures of lacing and shroud, on the one hand, dissipate vibration energy. On the other hand, the structure not only changes the mass of the blade but also changes the stiffness of the blade so that the inherent frequency and the exciting force frequency of the blade could avoid each other well. Thus the purpose of reducing the vibration level of the blade is achieved.
The contact surface structure of dry friction damper can be described with the fractal theory accurately. Since the fractal theory was proposed by Benoit B. Mandelbrot in 1975, it has been widely applied. Based on fractal dimension, Ritchie and Olff [1] proposed biodiversity calculation model by researching the relationship between space scaling relation and biodiversity. Taylor et al. [2] studied Pollock’s watercolors by using the fractal theory.
Using the dry friction model to inhibit vibration is a very effective method. The famous scholars Den Hartog [3] proposed the Coulomb model. After fully developing the Coulomb model and intensively studying dry friction, Meng et al. did further research on twodimensional sliding movement on the basis of the microscopic sliding study in 1991 [4, 5]. Muszynska et al. [6–8] simplified the blade and damper structure using the centralized quality model of multidegrees of freedom. Csaba [9] then combined the onebar microslip model with the centralized quality model of two degrees of freedom to study the blade with a damper. Lui et al. [10] and others proposed a new model to research the dry friction under different roughness and positive pressure condition.
In order to accurately describe contact interface and reveal the vibration reduction mechanisms of lacing and shroud structure, a centralized quality cantilever model of two degrees of freedom is built using the fractal theory in this paper. By placing the damper on different masses, the vibration reduction effect of lacing and shroud structure is simulated. By using the onebar microslip model, the constitutive equation of the force among friction surface is established. The timefrequency domain interactive method is adopted in the solving process of lacing and shroud structure. Therefore, the corresponding results at the steady state are obtained.
2. The OneBar Microslip Model
For damping blade, the damper structure and boundary conditions are very complex. The contact interface, especially, consists of a series of concave and convex structure in different sizes. It has high selfsimilarity and it is difficult to describe the structure using the same size. In order to develop a model as simple as possible but also complete enough to show the most important properties of friction interface, the contact interface is described using fractal dimension and fractal length in this paper. Fractal dimension can be calculated by a variety of methods. And similarity dimension method is adopted here.
Assuming objects or geometric shapes can be divided into parts, each part is similar to the whole object with a likelihood ratio . The fractal dimension can be expressed as
The fractal dimension does not have to be an integer.
The fractal curve length estimation model given by Mandelbrot is expressed as where is the fractal curve’s Euclidean length; is the initial operation length of fractal curve; is the scale of fractal curve; is the fractal dimension.
The damper is modelled by a rectangular bar pressed against a rigid surface with a normal load and subjected to a force , shown in Figure 1. The normal load on the bar is assumed to be constant over the width of the bar and defined by a quadratic normal load function in the lengthwise direction:
The force can be expressed as The bar has a modulus of elasticity and a cross section area . The length of the bar is , and the coefficient of friction is assumed to be constant across the contact zone and independent whether the bar is sliding or not.
The bar starts to slip, as the force is applied. Then the bar may be divided into two zones, one that is slipping and one that is stuck. The length of the slip zone is defined as the slip length and the slip length corresponding to the maximum force is , shown in Figure 2(a).
(a) Stuck and slip zones when initially loading the bar
(b) Illustrative plot of force and displacement when is decreasing
(c) Illustrative plot of force and displacement when is increasing
The zone that is slipping is stretched after the initial loading with . As the force decreases from to , the bar may be divided into three zones in this analysis, which is shown in Figure 2(b). Zone is stuck and has zero strain. Zone is stuck and stretched. Zone is slipping and compressed.
The length of the compression zone is denoted as the slip length . Zone will increase and zone will decrease as decreases. This will continue until . Zone is then eliminated and equals .
Here we have the opposite situation as in the previous part. As the force increases from to , the bar may be divided into three zones in this analysis, which is shown in Figure 2(c). Zone is stuck and has zero strain. Zone is stuck and compressed. Zone is slipping and stretched.
The length of stretched zone is denoted as the slip length . Zone will increase and zone will decrease as increases. This will continue until . Zone is then eliminated and equals .
The length of the slip zone is defined as the slip length . The slip length corresponding to the maximum force is . The displacement of the bar is defined as . According to the mechanics analysis, it can be seen that force and displacement are the function of slip length .
As the force is applied, the force at the right bar end and displacement can be expressed as As the force decreases from to , the force at the right bar end and displacement can be expressed as As the force decreases from to , the force at the right bar end and displacement can be expressed as Tow linearization techniques are discussed to transform the nonlinear properties of friction into equivalent damping and stiffness. Two criteria are established in Lazan’s linearization method, for equivalence of the two loops.(a)The same value of loop area is the same as damping energy per cycle .(b)The amplitude of force and displacement should be the same.
The first criterion gives the equivalent viscous damping , while the second criterion together with gives the equivalent stiffness :
3. The Forced Vibration Response Analysis of Damped Blade
A twodegree freedom system simplifying the blade is used for forced vibration response analysis in this paper. The 2DOF system is built up of a massless beam, two concentrated masses, and the damper. Compared with the previous study, this paper simulates lacing and shroud structure by locating damper on and , respectively, and the structure is shown in Figures 3 and 4. The damper is described by a viscous damper and spring, and and and , respectively. Displacements where the damper is attached and the force is applied are and , respectively.
The following parameters of the whole system are as follows, assuming that : beam system N/m, N/m, kg, and kg; damper system N, m, and .
The equations of motion for the system are
Assuming that harmonic motion yields we define the complex stiffness as
Solving equations for and yields
The equations of motion for the system are
Assuming that harmonic motion yields we define the complex stiffness as
Solving equations for and yields
4. Results and Discussion
Through blades vibration analysis with lacing in Figure 3, amplitudefrequency response curves of and in the range of 0~3000 Hz are achieved. The results are shown in Figures 5 and 6. Different positive pressure and damper parameters under different resonance states are shown in Table 1.

From Figures 5 and 6, it is obvious that, when the damper is placed on the mass to simulate lacing structure, the resonance amplitude of decreases with the increase of the positive pressure. On the contrary, the resonance frequency of increases with the increase of positive pressure. It could also be seen that, with the increase of positive pressure, the equivalent stiffness increases. From the vibration mechanics, the resonant frequency is positively correlated with the system stiffness as the resonance amplitudes negatively correlated with the stiffness. As a result, the resonance amplitude of decreases and the resonant frequency increases.
The resonance amplitude of decreases first and then increases along with the increase of positive pressure. Moreover, the changing trend of the resonance amplitude in Figure 6 is consistent with the conclusion as above. While the positive pressure equals 320, the resonance amplitude reaches the minimum. This indicates that there exists optimal positive pressure which makes blade resonance response reach minimum for a blade with lacing structure. As the same as mass , the resonance frequency of increases with the increase of positive pressure.
Comparing the amplitudes of the mass and , the amplitude of is found greater than that of in calculation of all frequencies. In this case, the vibration inhibition effect towards that damping structure brought is much greater than the inhibition effect towards .
The damper is placed on the mass to simulate shroud structure. Figures 7 and 8 present the amplitudefrequency response curves in the frequency range of 0~3200 Hz of and with shroud structure. Table 2 shows the resonance parameters of the damper under various positive pressures.

As shown in Figures 7 and 8, when the damper is placed on the mass , the resonance amplitude of decreases first and then increases along with the increase of positive pressure. While the positive pressure equals 160, the resonance amplitude reaches the minimum. And the resonance frequency of increases gradually with the increase of positive pressure. On the other hand, the resonance amplitude of decreases with the increase of positive pressure. The resonance frequency of increases along with the increase of positive pressure.
Comparing the amplitude of and , it is found that, when the positive pressure is smaller than or equal to 80, the amplitude of is greater than the amplitude of and the difference between the two amplitudes is smaller than the difference obtained when the damper is placed on . When the positive pressure is greater than 80, the amplitude of is greater than that of while the frequency is near the resonance frequency. However, the amplitude of remains smaller than that of when the frequency is beyond this range. In this case, the vibration inhibition effect towards that damping structure brought is much greater than the inhibition effect towards . When the positive pressure is relatively small, the inhibition effect could be small accordingly. But if the positive pressure is relatively large, particularly under the resonance state, the equivalent stiffness of the damper is very large and the inhibition effect towards would be great as well. This results in the fact that the amplitude of stays smaller than that of within certain range around the resonance frequency.
Comparing the results obtained when the damper is placed separately on and , it is found that the resonance amplitude changing trend of in the first condition is the same as the changing trend of in the second condition. And the resonance amplitude changing trend of is the same as the changing trend of in the second condition. But, under the same positive pressure, the resonance frequency of with the damper is higher than that of with the damper. By the contrast of Tables 1 and 2, it is obvious that the equivalent stiffness of with the damper is much greater than that of with the damper. So, in the same situation, the resonance frequency of increases drastically.
Figures 9 and 10 show hysteresis loops in an excitation cycle of the damper at resonance state. Table 3 gives energy dissipation in a vibration cycle of the damper at resonance state.

The area surrounded by a hysteresis loop represents the energy dissipation amount in a vibration cycle of the damper. From Figures 9 and 10 and Table 3, it could be seen that, whether the damper is placed on or , the energy dissipation amount of the damper decreases overall along with the increase of positive pressure. Based on formula (8), the energy dissipation amount in a vibration cycle of the damper is relevant to the positive pressure and the sliding part length . The energy dissipation amount is positively correlated with the positive pressure and negatively correlated with the sliding part length.
Under the same condition, the consumed energy of vibration when the damper is placed on lacing is greater than that when the damper is placed on shroud structure. Comparing Tables 1 and 2 and Figures 9 and 10, it is found that, when the damper is placed on lacing structure, the maximum excitation force on the damper presents linear variation with the increase of pressure. And the excitation force on the damper increases as the positive pressure increases. When damper is placed on shroud structure, the maximum excitation force on the damper presents linear variation with the increase of pressure. The variation amplitude is pretty small and decreases as the positive pressure increases.
Assuming the placing of the damper on the blade root, as the vibration displacement amplitude at the root is close to zero, the consumed energy in a vibration cycle is less than the consumed energy at the lacing structure. Meanwhile, from Figures 9 and 10 and Table 3, it could be seen that, under the same positive pressure, the consumed energy when the damper is placed on lacing wire is greater than that while the damper is placed on shroud structure. Combining mathematical knowledge, it could be obtained that there existed one position so that, when the damper is placed on this position, the consumed energy in a vibration cycle reaches maximum.
5. Conclusions
(1)Because of the complexity and selfsimilarity of contact surface geometry, contact surface can be described accurately using fractal dimension and fractal length.(2)Whether shroud or lacing structure, it not only has the effect of dissipating energy but also changes blade stiffness. Meanwhile, in the same condition of pressure and contact area, the damping effect of lacing structure is stronger than that of shroud structure, but the stiffness effect is opposite.(3)When damper is located on mass , as pressure increases, the resonance frequencies of masses and increase. The resonance amplitude of mass decreases and the resonance amplitude of mass decreases first and then increases. It is obvious that there is an optimal pressure, making resonance amplitude minimum. Moreover, the resonance amplitude of is less than that of .(4)When damper is located on mass , as pressure increases, the resonance frequencies of masses and increase. The resonance amplitude of mass decreases first and then increases while the resonance amplitude of mass decreases. There is an optimal pressure that makes resonance amplitude of mass minimum. When pressure is relatively small, the resonance amplitude of is smaller than that of . When pressure is relatively large, near the resonance frequency, there is a frequency range that makes resonance amplitude of bigger than that of .(5)Under the same circumstance, resonance frequency is lower when damper is located on than when it is located on , but the vibration amplitude is opposite.(6)The damper consumes more energy when it is located on mass than the situation when it is located on . No matter which mass damper is located, the energy consumed in a vibration cycle declined as the pressure increases.(7)Under the same circumstance, there is at least one position on the blade, making damper energy consumption maximum in one vibration cycle.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 M. E. Ritchie and H. Olff, “Spatial scaling laws yield a synthetic theory of biodiversity,” Nature, vol. 400, no. 6744, pp. 557–560, 1999. View at: Publisher Site  Google Scholar
 R. P. Taylor, A. P. Micolich, and D. Jonas, “Fractal analysis of Pollock's drip paintings,” Nature, vol. 399, no. 6735, p. 422, 1999. View at: Google Scholar
 J. P. Den Hartog, “Forced vibrations with combined Coulomb and viscous friction,” Transactions of the American Society of Mechanical Engineers, vol. 53, no. 9, pp. 107–115, 1931. View at: Google Scholar
 C.H. Meng, J. Bielak, and J. H. Griffin, “The influence of microslip on vibratory response, part I: a new microslip model,” Journal of Sound and Vibration, vol. 107, no. 2, pp. 279–293, 1986. View at: Google Scholar
 C.H. Meng and P. Chidamparam, “Friction damping of twodimensional motion and its application in vibration control,” Journal of Sound and Vibration, vol. 144, no. 3, pp. 427–447, 1991. View at: Google Scholar
 A. Muszyńska and D. I. G. Jones, “On tuned bladed disk dynamics: some aspects of friction related mistuning,” Journal of Sound and Vibration, vol. 86, no. 1, pp. 107–128, 1983. View at: Google Scholar  Zentralblatt MATH
 A. Muszynska, D. I. G. Jones, T. Lagnese, and L. Whitford, “On nonlinear response of multiple blade systems,” Shock and Vibration Bulletin, vol. 51, no. 3, pp. 89–110, 1981. View at: Google Scholar
 D. I. G. Jones and A. Muszynska, “Design of turbine blades for effective slip damping at high rotational speeds,” Shock and Vibration Bulletin, vol. 49, pp. 87–96, 1979. View at: Google Scholar
 G. Csaba, Modelling Microslip Friction Damping and Its Influence on Turbine Blade Vibrations, Division of Machine Design Department of Mechanical Engineering, Linköping University, 1998.
 Y. Liu, B. Shangguan, and Z. Xu, “A friction contact stiffness model of fractal geometry in forced response analysis of a shrouded blade,” Nonlinear Dynamics, vol. 70, no. 3, pp. 2247–2257, 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Ruishan Yuan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.