Dynamics Modeling and Numerical Analysis of Rotor with Elastic Support/Dry Friction Dampers

，，，

School of Power and Energy，Northwestern Polytechnical University，Xi′an 710072，P.R.China

0 Introduction

In some cases dry friction is disadvantageous in mechanical systems，such as in rubbing faults of rotor systems［1－4］.However，dry friction has been widely used to increase the stability of mechanical systems，such as the dry friction damping blade of aero－engines and turbines［5－8］.

In order to restrain the vibration of rotor system，series of measures have been taken by researchers，such as magnetorheological fluid damper and dry friction damper［9－10］.It has been proved theoretically and experimentally that the elastic support/dry friction damper can significantly attenuate the vibration amplitude of a rotor system in the critical speed region［11－14］.Fan proposed that friction damping could be applied in rotor system［12］.He designed an elastic support/dry friction damper using a spring to provide the pressing force and experimentally studied the damping performance of this damper.The experimental results showed a great damping performance.Based on his work，Wang Siji designed an active elastic support/dry friction damper using an electro－magnetic actuator，which can conveniently be actively controlled by adjusting the control voltage of the electromagnet［15－17］.For different conditions，this active damper can be operated under different control strategies.Wang Siji also studied its performance of the damper for some rotor faults such as a crack in the elastic support，sudden blade－out event and unstable vibration of the rotor and found that this damper has a good damping effect［18－21］.

For theoretical study of the elastic support/dry friction damper，only a simple mechanism model was established in which the rotor was simplified to a single－degree－of－freedom system or a two－degree－of－freedom system，and the classic one－dimensional hysteretic dry friction model was used［12］.This type of analysis model was sufficient for illustrating the friction damping mechanism of the dry friction damper，but it is limited and insufficient for the next study on the damping performance with variable design parameters.The limitation is mainly reflected in the following two points：（1）The axial location of the damper influences the damping performance of the elastic support/dry friction damper，which cannot be studied in the current analysis model.（2）The relative motion between the friction pair is a twodimensional motion，which is considered one－dimensionally in the current analysis model.

In this paper，a two－dimensional friction model－ball/plate model based on the classic onedimensional hysteretic dry friction model is proposed.By using the ball/plate model，a rotordynamic model for rotordynamic analysis of a rotor with elastic support/dry friction dampers was established and experimentally verified.Finally，the damping performance with respect to some variable parameters has been studied numerically，laying the foundation for the dynamic design of the elastic support/dry friction damper.

1 Basic Operation Principles

Fig.1shows the basic operation principles of a rotor with elastic support/dry friction dampers.The rotor is supported by two elastic supports，and at the end cross section of each one，the dry friction damper is affixed.Each damper consists of three key components：The elastic support，the friction pairs（stationary disk and moving disk）and the actuator.The elastic supports redistribute the strain energy of the whole rotor－support system and concentrate the vibration，which is then dissipated by the dry friction between the friction pairs.The moving disk，which is fixed to the end cross section of the elastic support，is connected to the bearing outer ring and vibrates with the elastic support but does not rotate with the rotor.The stationary disk is fixed to the casing and can be moved in the axial direction by the actuator.The frictional force between the two disks can be changed by adjusting the pressing force from the actuator.So the actuator is the core component of the damper.

If the moving disk fixed to the end cross section of the elastic support vibrates with the rotor vibration，a relative motion between the moving disk and the stationary disk will take place.When the actuator provides an appropriate pressing force，this will lead to dry friction damping dissipating the vibration energy.

Fig.1 Basic operation principles of the elastic support/dry friction damper

2 Rotordynamic Modeling

2.1 Friction model

Fan and Wang［12，16］studied the elastic support/dry friction damper on a rotor by using the hysteretic dry friction model shown in Fig.2，where k is the stiffness of the contact interface.

The whole rotor system with dampers was simplified as a mass－spring system shown in Fig.3，where M is the equivalent mass，k the equivalent stiffness，c the equivalent damp，f（τ）the load on the system，andxandzrepresent the displacement.The basic characteristics of the damper can be explained by this simplified model.However，for further study，this simplified model is inadequate.

The inadequacy of the model in Fig.2is that in a rotor system，a two－dimensional translation may occur between the moving disk and the stationary disk，and the friction force may be a circumferential force，but in the model in Fig.2，only one－dimensional motion was considered.The rotor model in Fig.3was oversimplified for further study of the dynamics.

Fig.2 Hysteretic dry friction model

Fig.3 Simplified model of a rotor system with an elastic support/dry friction damper［12］

Based on the above considerations，a two－dimensional friction model－ball/plate model（shown in Fig.4）is proposed.This ball/plate model was developed from the hysteretic dry friction model，and it is clear and simple to describe the rotor system with elastic support/dry friction dampers.

As shown in Fig.4，in the ball/plate dry friction model，the stationary disk is represented by a rectangle.The moving disk consists of a ball and a plate.The plate （without considering its mass）represents the contact interface between the moving disk and the stationary disk，and the ball represents the moving disk.The ball and the plate are connected with ideal springs and linear damping in two directions.The ideal springs represent the tangential contact stiffness of the contact interface.The displacement between the ball and the plate represents the microscopic sliding in the state of static friction.So if the applied force on the ball is greater than maximum static friction，the ball will touch the edge of the plate，and the plate will begin to move.

Fig.4 Ball/plate dry friction model

Both the stationary disk and the moving disk have two degrees of freedom in the ball/plate dry friction model.When the stationary disk is fixed，its displacement remains zero in each degree of freedom.But when there is a stiffness between the stationary disk and ground，the displacement will depend on the friction force of the contact interface.

The ball/plate dry friction model is actually a type of two－dimensional hysteretic dry friction model.The frictional force between the stationary disk and the plate is Coulomb′s friction force.The frictional force at any time depends on the motion state of the plate and the applied force with which the ball acts on the plate.

Without considering the mass of the plate，the resultant force acting on the plate will be 0at any time.When the plate is stationary relative to the stationary disk，namely the friction pair is in the state of static friction，the frictional force acting on the plate is

In Eqs.（1，2）：r1＝x1＋iy1is the displacement of the ball；r2＝x2＋iy2is the displacement of the plate；rj＝xj＋iyjis the displacement of the stationary disk，its value is 0when the stationary disk remains at rest；k is the stiffness coefficient between the ball and the plate；d is the damping coefficient between the ball and the plate；μis the friction coefficient；and N is the pressing force.

By mechanical analysis of this two－dimensional friction model，the motion equations of the friction pair can be established conveniently.

2.2 Rotor system

A single－disk flexible rotor with elastic support/dry friction dampers is shown in Fig.5.The system consists of two parts：A rotor and two elastic support/dry friction dampers（in the dashed boxes）.The rotor is a single offset disk with a flexible shaft that is supported by two elastic supports at both ends.The two elastic support/dry friction dampers can be set up at each elastic support，shown in the dashed boxes of Fig.5.In Fig.5，L is the length between the two supports；a is the length between the left support and the rotor disk；b is the length between the right support and the rotor disk；r is the diameter of the rotor disk；εis the eccentricity of the rotor disk；t is the diameter of the shaft.

Fig.5 A single－disk flexible rotor with elastic support/dry friction dampers

For the rotor disk，four coordinates are required to describe its motion.Two of them （x，y）describe its two－dimensional translation；the other two（φx，φy）describe its two－dimensional swing.For the left friction pair，based on Section 2.1，six coordinates are required to describe the two－dimensional motion of the moving disk （the ball and the plate）and the stationary disk：xb1and yb1describe the motion of the ball；xd1and yd1describe the motion of the plate；and xj1and yj1describe the motion of the stationary disk.In the same way，another six coordinates describe the motion of the right friction pair.They are xb2，yb2，xd2，yd2，xj2，yj2.So in total，the whole rotor system is a system with 16degrees of freedom（16－DOF system）.

2.3 Motion equations

By means of Euler′s laws of motion，the motion equations of the rotor disk，moving disk and stationary disk can be obtained as follows.

The rotor disk

The moving disk （the ball in the ball/plate model）

In Eqs.（3—10），r ＝x＋iyis the displacement of the rotor disk；φ ＝φx＋iφyis the swing angle of the rotor disk；rb1＝ xb1＋iygbg1is the displacement of the left moving disk；rb2＝xb2＋iyb2is the displacement of the right moving disk；rd1＝xd1＋iyd1is the displacement of the left plate of the moving disk；rd2＝ xd2＋iyd2is the displacement of the right plate of the moving disk；rj1＝ xj1＋iyj1is the displacement of the left stationary disk；rj2＝ xj2＋iyj2is the displacement of the right stationary disk；m，mb1，mb2，mj1，and mj2represent the masses of the rotor disk，the moving disks and the stationary disks，respectively；s11，s12，s21and s22represent the stiffness coefficients of the shaft at the rotor disk，where s11is the displacement stiffness，s22the angle stiffness，s12and s21are cross stiffness；sb1，sb2，sj1，and sj2represent the stiffness coefficients of the elastic support and the stationary disk respectively；db1，db2，dj1，and dj2represent the damping coefficients of the elastic support and the stationary disk，respectively；Jpand Jdrespectively represent the polar moment of inertia and the moment of inertia about a diameter of the rotor disk；N1and N2represent the pressing force between the moving disk and stationary disk of the two dampers；k and d are the tangential contact stiffness coefficient and damping coefficient of the contact interface introduced by the ball/plate model；andΩis the rotational speed of the rotor.

To solve the equations by numerical methods，the following variables are introduced

Eq.（13）can be solved conveniently by numerical method for ordinary differential equations.

3 Model Verification

In the following numerical analysis，the model is verified using some test data，which are obtained from experiments with two types of elastic support/dry friction dampers：One is uncontrollable and the other one is controllable.

3.1 Time domain response of a rotor with elastic support/dry friction dampers under on－off control

The time domain responses were tested by a rotor with two active elastic support/dry friction dampers under on－off control.And the parameters of the rotor test rig had been given in Ref.［16］.This type of damper can be controlled by electricity.The two dampers were fixed to each elastic support.The test results are shown in Figs.6，7.

Figs.6，7are the time domain responses of the rotor under on－off control.The control voltagesUare 8Vand 10V.WhenU＝8V，the pressing forces between the moving disk and the stationary disk of the two dampers on both ends are 58Nand 70N，and whenU＝10V，the pressing forces are 87Nand 104N.

Numerical simulations were carried out by using the parameters of the rotor test rig.In the simulations，the eccentricityεis 13×10－5m，the rotation speedΩis 1 410r/min，the friction coefficientμis 0.25，the tangential contact stiffness of the contact interface k is 3×105N/m，and the pressing forces of two dampers on both ends are 58，70Nand 87，104N.The results are shown in Figs.8，9.

From Fig.8（b）and Fig.9（b），the peak amplitude shows a fluctuation when the damper is switched off.This phenomenon is mainly caused by transient excitation.The transient excitation contains the natural frequency of the rotor system，which leads to a free vibration.Because the natural frequency is close to the rotation speed，it causes a beat vibration.Due to the system damping，the free vibration decays with time and the peak amplitude of rotor is gradually stabilized.What′s more，the larger the damping is，the faster the free vibration decays.So when the damper is switched on，the beat vibration phenomenon is slight in Fig.8（a）and almost nonexistent in Fig.9（a）.

Fig.6 Test results of time domain responses of the rotor（Ω ＝1 407r/min，U＝8V）

Fig.7 Test results of time domain responses of the rotor（Ω ＝1 407r/min，U＝10V）

Fig.8 Calculation results of time domain responses of the rotor（Ω＝1 410r/min，N1＝58N，N2＝70N）

Fig.9 Calculation results of time domain responses of the rotor（Ω＝1 410r/min，N1＝87N，N2＝104N）

Comparing each plot in Fig.8with Fig.6and that in Fig.9with Fig.7，the calculation results show good agreement with the test results.And there are some differences between the experimental damping effect and the calculation results.The main reason is that the pressing force can be affected by the assembly parameters of the damper，which may cause differences in calculation results.

3.2 Amplitude－frequency characteristics of a rotor with elastic support/dry friction dampers under variable pressing force

With the test rig given in Ref.［12］，the amplitude－frequency characteristics of a rotor with elastic support/dry friction dampers under variable pressing force were tested，and two friction pairs were used in the experiments.One was brass/steel and the other was steel/steel.The test results are shown in Fig.10（a）and Fig.11（a）.Numerical simulations were carried out using the parameters of the rotor test rig.In the simulations，the eccentricities ε are 5.3 × 10－5m and 2×10－5m，the friction coefficientsμare 0.19（brass/steel）and 0.1（steel/steel），the tangential contact stiffness k of the contact interfaces are 1×106N/m （brass/steel）and 1×108N/m （steel/steel），and the rotational speeds are 500—3 000 r/min and 500—4 000r/min.The results are shown in Fig.10（b）and Fig.11（b）.

Fig.10 Amplitude－frequency characteristics of the rotor（brass/steel）

Fig.10shows the amplitude－frequency characteristics when the friction pair is brass/steel，and Fig.11shows the amplitude－frequency characteristics when the friction pair is steel/steel.The specific values of the peak amplitude in Figs.10，11are shown in Table 1.The difference between the calculation results and the test results is less than 10%.So comparing test results with calculation results in both Fig.10and Fig.11，it is found that the calculation results show good agreement with the test results.

Fig.11 Amplitude－frequency characteristics of the rotor（steel/steel）

Table 1 Peak amplitude of the rotor system in test and calculation results

By the comparisons above，the dynamic model of the rotor with elastic support/dry friction dampers has been verified.The comparisons also show that the variable parameters are very important for the damping performance of the elastic support/dry friction damper.Numerical analysis of some variable parameters will be carried out in the next section.

4 Numerical Analysis

4.1 Basic parameters of the rotor

The geometric dimensions of the rotor for numerical analysis are shown in Fig.5.The values are as follows：L＝700mm，a＝250mm，b＝450mm，t＝28mm，r＝120mm，h＝40mm.

It is assumed that the material of the rotor is steel.Its densityρis 7.8×103kg/m3，and its elastic modulus E is 2.1×1011N/m2.

Based on the geometric dimensions and material parameters above，the mass，stiffness，moment of inertia，etc.，of the rotor can be obtained.In addition，in order to investigate the effect of mode shape on damping effect，mb2is valued much larger than mb1to make the vibration amplitude at the right elastic support higher than left on the second mode.

All of the basic parameters are listed in Table 2.

4.2 Mode and unbalance response of the rotor

The mode and unbalance response of the rotor can be obtained by solving the homogeneous equations corresponding to Eqs.（3—6）analytically or numerically.The results are shown in Fig.12，in which Fig.12（a）is the mode shape of the rotor and Fig.12（b）is the unbalance response of the rotor.

As shown in Fig.12，the vibration amplitudes of the left and right supports are approximately equal when the rotor is at the first mode，while at the second mode，the vibration of the rotor is concentrated in the right support.

Table 2 Basic parameters of the rotor for numerical analysis

Fig.12 Mode shape and unbalance response of the rotor

4.3 Damping performance

In this section，the damping performance of the elastic support/dry friction damper with respect to variable parameters such as the location of the damper，stiffness of the elastic support，pressing force between the moving and stationary disks，etc.，will be discussed with numerical analysis.

4.3.1 Location of the damper

The damping of the elastic support/dry friction damper is caused by the friction between the moving and stationary disks.The moving disk is fixed at the end cross section of the elastic support，so the damping performance of the damper is closely related to the vibration of the elastic support.

Fig.12（b）shows that the first and second critical speeds of the rotor are 1 748and 2 655r/min，so the numerical analysis of the damper′s performance with respect to the location of the damper was carried out at 1 748and 2 655r/min.

Fig.13shows the vibration displacement of the rotor disk at 1 748r/min（Fig.13（a））and 2 655r/min（Fig.13（b））.In Fig.13，all three lines represent the vibration of the disk in the rotor system.What makes these lines different is that the dotted line represents the rotor without dry friction damping；the solid line represents the rotor with dry friction damping on the left elastic support（N1＝100N，N2＝0N）；and the dashed line represents the rotor with dry friction damping on the right elastic support（N1＝0N，N2＝100N）.As shown in Fig.13（a）and Fig.13（b），for the first critical speed，the damping performances of the dampers on both elastic supports are almost equal.However，for the second critical speed，the damping performance of the damper on the right elastic support is better than that on the left elastic support.

Fig.13 Damping performance with respect to the location of the damper

The vibration energy of the rotor system is dissipated by the friction between the friction pairs.At 1 748r/min，the first mode shape of the rotor in Fig.12（a）shows that the vibrations of the two supports are almost equal，so in each individual rotation cycle，the vibration energy dissipated by the dampers on the left and right supports is almost the same，which means that the dampers on the two elastic supports show the same damping performance.At 2 655r/min，the second mode shape of the rotor in Fig.12（a）shows that the vibration amplitude of the right support is larger than that of the left support，so the vibration energy dissipated by the damper on the right elastic support is larger than that dissipated by the left elastic support，which means that the damper on the right elastic support shows a better performance than the damper on the left elastic support.

4.3.2 Stiffness of the elastic support

The critical speed and mode shape can be affected by the elastic support.To study the performance of the elastic support/dry friction damper with respect to the stiffness of the elastic support，let the stiffness of the elastic support be the variable parameters（sb1＝sb2＝3×103，7×103，and 12×105N/m）.The friction force is applied simultaneously on the left and right elastic supports（N1＝ N2＝20N）.

Fig.14is the mode shape of the rotor with the parameters（sb1＝sb2＝3×105，7×105，and 12×105N/m）.From Fig.14，we can see that as the stiffness coefficient of the elastic supports increases，the mode shapes of both the first and second modes become more and more bent，which means that the smaller the stiffness coefficient of the elastic supports，the more vibration energy is concentrated in the elastic support.

Fig.14 Mode shape of the rotor

Fig.15is the unbalance response of the rotor with the parameters（sb1＝sb2＝3×105，7×105，and 12×105N/m）.In this figure，the dotted line represents the response without friction（N1＝ N2＝0N），and the solid line represents the response with friction（N1＝ N2＝20N）.From the figure，we can see that the damper performance is best when the stiffness coefficient of the support is the smallest，3×105N/m；the damper performance is the worst when the stiffness coefficient of the support is the largest，12×105N/m.

Fig.15 Unbalance response of the rotor（sb1＝ sb2＝3×105，7×105，and 12×105 N/m）

As mentioned above，the smaller the stiffness coefficient of the elastic support，the more vibration energy is concentrated in the elastic support.Additionally，the more vibration energy is concentrated in the elastic support，the more conducive is for the damping performance of the elastic support/dry friction damper.So the damper performs better and better as the stiffness coefficient of the elastic support decreases.For the stiffness coefficient of the elastic support，the op－erating speed，static deformation，strength and reliability should be taken into consideration.

The above analysis indicates that the damping performance of the elastic support/dry friction damper is closely related to the operating speed of the rotor and the characteristics of the rotor′s mode.Even the same dampers fixed on different rotors or different support locations will perform differently.

4.3.3 Pressing force and friction coefficient

The damping provided by an elastic support/dry friction damper to a rotor is derived from the sliding friction force between the stationary and moving disks.So，as the product of the pressing forces N1and N2and the friction coefficientμ，the sliding friction forcesμN(yùn)1andμN(yùn)2directly determine the damping performance of the damper.

Let the sliding friction forcesμN(yùn)1andμN(yùn)2be the variable parameters，withμN(yùn)1＝ μN(yùn)2＝ 0，4.5，10.5，13.5，37.5，75，120，and 195N，applied simultaneously.Fig.17is the unbalance response of the rotor for the parameterμN(yùn)1＝μN(yùn)2＝0，4.5，10.5，13.5，37.5，75，120，and 195N.

Fig.16illustrates the sliding friction forces withμN(yùn)1andμN(yùn)2increasing.It is found that the critical speed region of the rotor system moves upward，and the peak amplitude of the rotor at the critical speed first decreases and then increases.When the sliding friction forces are large enough，the peak amplitude of the rotor even exceeds the peak amplitude without friction.There must be an optimal sliding friction force under which the unbalance response of the rotor system will be smallest for all rotational speeds，and the rotor can pass through the critical speed smoothly.For this rotor，the optimal sliding friction force is between 24Nand 37.5N.

Fig.16 Unbalance response of the rotor（μN(yùn)1＝ μN(yùn)2＝0，4.5，10.5，13.5，37.5，75，120，and 195N）

This phenomenon can be explained by analyzing the friction hysteresis loops of ball/plate dry friction model in Fig.17.It shows the relationship between the friction and displacement in a period，wherexandyrepresent the two－dimensional displacement，F(xiàn)frepresents the absolute value of friction whose direction is tangential to the circle.There are two moving trajectories in each period，bigger one is for the ball and smaller one for the plate.Both of them are moving at the rotation speed，and the ball′s phase is ahead of the plate.Their distance eyrepresents the microscopic sliding.Then，the flank area of the cylinder formed by the plate circle andx－yplane represents the friction work in a period.When the pressure gets smaller，the trajectory of friction goes to theμN(yùn)′－A′circle，on the contrary，the trajectory of friction goes to theμN(yùn)″－A″circle.Therefore，there is a maximum value of the area when the pressing force changes，also it′s the same for damping effect of damper.Furthermore，the elastic support/dry friction damper not only provides external damping to the rotor system，but also introduces extra stiffness into the rotor system.By calculating the response of the rotor system like Fig.16，the optimum value can be found.

Fig.17 Hysteresis loops of friction with displacement

4.3.4 Stiffness of stationary disk and tangential contact stiffness of contact interface

Between the moving disk and the mounting base of the stationary disk，there is a combined stiffness that consists of the stiffness of the stationary disk and the tangential contact stiffness of the contact interface.

（1）The stiffness of the stationary disk

Under the traction of the moving disk，the stationary disk will move slightly.Let the stiffness of the stationary disks sj1and sj2be the variable parameters，sj1＝ sj2＝ 0.05×107，0.2×107，1×107，and 15×107N/m，and the pressing forces of the two dampers，applied simultaneously，be N1＝150Nand N2＝150N.

Fig.18is the unbalance response of the rotor with the parameters sj1＝sj2＝0.05×107，0.2×107，1×107，and 15×107N/m，and N1＝150N，N2＝150N.The figure shows that as the stiffness coefficient of the stationary disk increases，the peak amplitude of the rotor decreases，while the damping performance of the elastic support/dry friction damper improves.When the stiffness coefficients sj1and sj2increase to some extent，the unbalance response curve is nearly constant and the damping performance of the damper no longer changes.

Fig.18 Unbalance response of the rotor（sj1＝sj2＝5×105，2×106，1×107，and 1.5×108 N/m；N1＝150N，N2＝150N）

Fig.19is the time domain waveform of the left support at 1 800r/min，with the stiffness coefficients sj1＝sj2＝0.05×107N/m and sj1＝sj2＝1×107N/m.As shown in Fig.19（a），when the stiffness coefficient is small，sj1＝sj2＝0.05×107N/m，the motion of the stationary disk under the traction of the moving disk is obvious，which makes the relative motion between the moving and stationary disks smaller，so it is unfavourable for the damping performance of the elastic support/dry friction damper.In Fig.19（b），when the stiffness coefficient is large，sj1＝sj2＝1×107N/m，the stationary disk barely moves，which is very favourable for the damper.

Fig.19 Time domain waveform of the left support（Ω＝1 800r/min，N1 ＝ 150N， N2 ＝150N）

（2）The tangential contact stiffness of the contact interface

In the ball/plate model（Fig.4），the tangential contact stiffness of the contact interface is represented by an ideal spring with stiffness coefficient k.Let this stiffness coefficient k be the variable parameter，k＝0.8×105，3×105，and 10×105N/m，and the pressing forces of the two dampers，applied simultaneously，N1＝ 150N and N2＝150N.

Fig.20is the unbalance response of the rotor under the parameter k＝0.8×105，3×105，and 10×105N/m，and N1＝150N，N2＝ 150N.From Fig.20，we can see that as the stiffness coefficient k increases，the peak amplitude of the rotor decreases and the damping performance of the elastic support/dry friction damper improves.

In the ball/plate model（Fig.4），the plate moves under the traction of the ball.The damping comes from the relative motion between the plate and the stationary disk.When the tangential contact stiffness k is small，the relative motion between the plate and the stationary disk is small，so the damping will be small.

Fig.20 Unbalance response of the rotor（k＝8×104，3×105，and 1.0×106 N/m；N1＝150N，N2＝150N

Based on the analysis in this section，it′s clear that the stiffness of the stationary disk sjand the tangential contact stiffness of the contact interface k are connected in series between the moving disk and the mounting base of the stationary disk.The larger this combined stiffness is，the better the damper′s damping performance is.

In addition，the analysis in this paper is mainly based on the circular whirl assumption.In this situation，the trajectory of the moving disk is a circle.So the value of friction does not change with the response of the rotor system in each period，but its direction changes with time.If the circular whirl assumption is invalid，there will be transformation between kinetic and static friction，which may leads to nonlinear phenomenon.Further research is needed in this situation.

5 Conclusions

In this paper，a two－dimensional friction model－ball/plate model was proposed by which the rotordynamic model of a rotor with elastic support/dry friction dampers was established and experimental verified.Lastly，the damping performance with some variable parameters was studied numerically.The main conclusions are summarized as follows：

（1）A two－dimensional friction model－ball/plate model was proposed.This model is clear and simple，and it can be used to establish the rotordynamic model of a rotor with elastic support/dry friction dampers.

（2）The damping performance of the elastic support/dry friction damper is closely related to the characteristics of the rotor′s mode.To obtain the damper′s best performance，the damper should be located at the elastic support in which the vibration energy is concentrated.

（3）For vibration control，the friction forces should not be maximized，as there is an optimal value at which the elastic support/dry friction damper performs best.The damper not only provides external damping to the rotor system，but also introduces extra stiffness into the rotor system.

（4）The stiffness of the stationary disk and the tangential contact stiffness of the contact interface are connected in series between the moving disk and the mounting base of the stationary disk.The larger this combined stiffness is，the better the damper′s damping performance is.

Acknowledgement

This work was supported by the National Natural Science Foundation of China（No.51405393）.

［1］ ABU－MAHHOUZ I，BANERJEE A.On the investigation of nonlinear dynamics of a rotor with rub－impact using numerical analysis and evolutionary algorithms［J］.Procedia Computer Science，2013，20（1）：140－147.

［2］ DENG X，LIAO M F，LIEBICH R，et al.Experimental research of bending and torsional vibrations of a double disc rotor due to rotor－to－stator contacts［J］.Journal of Aerospace Power，2002，17（2）：205－211.

［3］ OBERST S，LAI J C S，MARBURG S.Guidelines for numerical vibration and acoustic analysis of disc brake squeal using simple models of brake systems［J］.Journal of Sound ＆ Vibration，2013，332（9）：2284－2299.

［5］ GRIFFIN J H，MENQ C H.Friction damping of circular motion and its implications to vibration control［J］.Journal of Vibration and Acoustics，1991，113（2）：225－229.

［6］ QI W K，ZHANG Y J.Reduced vibration characteristics of turbine blade with platform damper［J］.Joernal of Nanjing University of Aeronautics ＆ Astronautics.2014，46（2）：280－284.（in Chinese）

［7］ HE B B，OUYANG H J，REN X M，et al.Dynamic response of a simplified turbine blade model with under－platform dry friction dampers considering normal load variation ［J］.Applied Sciences，2017，7（3）：228.

［8］ DING Q，CHEN Y.Analyzing resonant response of a system with dry friction damper using an analytical method［J］.Journal of Vibration ＆ Control，2008，14（8）：1111－1123.

［9］ WANG L Q，LI W Z，GU L，et al.Reducing vibration with friction－damping in high－speed rotor system［J］.Transactions of Nanjing University of Aeronautics ＆ Astronautics，2007，24（1）：48－53.

［10］XING J，HE L D，WANG K.Optimizing control for rotor vibration with magnetorheological fluid damper［J］.Transactions of Nanjing University of Aeronautics ＆ Astronautics，2014，31（5）：538－545.

［11］FAN T Y，LIAO M F.Dynamic behavior of a rotor with dry friction dampers［J］.Mechanical Science and Technology，2003，22（5）：743－745.（in Chinese）

［12］FAN T Y.Vibration reduction by elastic support dry friction damper［D］.Xi′an：Northwestern Polytechnical University，2006.（in Chinese）

［13］WANG S J，LIAO M F，YANG S J.Experimental investigation on vibration control by elastic support/dry friction damper［J］.Journal of Aerospace Power，2007，22（11）：1893－1897.（in Chinese）

［14］WANG S J，LIAO M F.Experimental investigation of an active elastic support/dry friction damper on vibration control of rotor systems［J］.International Journal of Turbo ＆Jet Engines，2014，31（1）：13－17.

［15］WANG S J，LIAO M F.Control strategy and methods of rotor systems by an elastic support/dry friction damper［J］.Journal of Aerospace Power，2011，26（10）：2014－2219.（in Chinese）

［16］ WANG S J.Vibration control techniques for rotor systems by an active elastic support/dry friction damper［D］.Xi′an：Northwestern Polytechnical University，2008.（in Chinese）

［17］WANG S J，LIAO M F.Online Control of rotor system instability by elastic support/dry friction damper［J］.Journal of Vibration Measurement ＆ Diagnosis，2012，32（2）：323－327.（in Chinese）

［18］WANG S J，LIAO M F.Application of an active elastic support/dry friction damper to control sudden unbalance response of rotor systems［J］.Mechanical Science and Technology，2008，27（5）：667－672.（in Chinese）

［19］LIAO M F，SONG M B，WANG S J.Active elastic support/dry friction damper with piezoelectric ceramic actuator［J］.Shock and Vibration，2014，2014（2）：1－10.

［20］SONG M B.Dynamic design of elastic support/dry friction damper matching rotor［D］.Xi′an：Northwestern Polytechnical University，2016.（in Chinese）

［21］WANG S J，LIAO M F.Protection of a rotor against violent vibrations caused by a crack in elastic supports by using active elastic support/dry friction dampers［J］.Journal of Aerospace Power，2008，23（11）：2026－2030.（in Chinese）