Rub-Impact Detection in Rotor Systems with Pedestal Looseness Using a Nonlinearity Evaluation

In this paper, a nonlinearity evaluation is proposed in order to identify the rub-impact in rotor systems with pedestal looseness. Nonlinear mathematical models have been established for bearing-rotor systems with single pedestal looseness and pedestal looseness coupled with rub-impact. Piecewise linear stiffness and damping are considered regarding the position of pedestal looseness, while radial impact forces are defined using the Colulomb type of frictional relationship during rub-impact. +e nonlinearity evaluation is employed to quantify the nonlinearity of the dynamics of bearing-rotor systems, which are calculated at different looseness clearances.+e experiments for rotor systems with pure pedestal looseness and pedestal looseness coupled with rub-impact are conducted respectively to collect the vibration signals on different looseness clearances. Two different curves are obtained using the nonlinear fitting method for the values of nonlinearity evaluation. +e rub-impact within rotor systems with pedestal looseness can then be identified by comparing the curves that denote the trend of nonlinearity evaluation for the measured vibration responses.


Introduction
Pedestal looseness is a common issue in rotating machinery.
e excessive vibration caused by pedestal looseness often results in a secondary phenomenon, rub-impact [1].Rotorto-stator rub remains one of the primary causes of serious malfunctions that occur in rotating machinery.In recent years, there has been a focus on determining whether rubbing is occurring in a rotor system [2][3][4][5][6][7].However, the partial rub and pedestal looseness have the same contact and separation process, and the dynamic characteristics and signatures of the measured vibration signals of a rotor system will appear similar.us, distinguishing between partial rub and pedestal looseness by measured signals is challenging [8].Using signal-based methods alone, it is difficult to identify rub-impact from rotor systems with pedestal looseness in rotating machinery.
ere has been extensive research in the past decades regarding the detection of rubbing faults in a rotor system.e methods of detection can be divided into signals-based identification methods [9][10][11][12] and model-based identification methods [13][14][15][16].e former uses the measured vibrations to describe the dynamics of the rotor system to detect the possible underlying reason for rubbing.e various signal processing techniques such as FFT, STFT, Hilbert-Huang transform (HHT) [8], ensemble local means decomposition [17], and variational mode decomposition (VMD) [18] are widely used and have been proven to be powerful.To overcome the difficulty in detecting the rubbing location, the model-based identification method is employed to provide quantitative information regarding the faults, including locations and severity.is method combines the dynamic characteristics and the dynamic model of the rotor systems with the signal processing techniques [19][20][21][22].Despite a lot of research focusing on the analysis of the rotor systems with pedestal looseness [23][24][25][26][27][28], these methods have not been applied to identify the rub-impact in rotor systems with pedestal looseness.
Rotor systems with pedestal looseness can be described by a nonlinear mathematical model containing stiffness and damping forces with trilinear forms [28].Due to the existence of nonlinear bearing forces within the system, a bearing-rotor system with pedestal looseness demonstrates nonlinear dynamic behaviors, with complex and nonstationary vibration characteristics.When rub-impact occurs in rotor systems with pedestal looseness, the nonlinear rubbing forces with a piecewise linear form of stiffness are also considered in the system [29,30].Rotor systems with pedestal looseness coupled with rub-impact will become much more complicated, and its nonlinearity effects on the dynamics of the system will become increasingly more serious.Numerically, the severity of the nonlinearity can be quantified by a nonlinearity evaluation (also called a nonlinearity measure or a nonlinearity test), which represents an approach to systematically quantify the degree of nonlinearity of a dynamical system [31,32].
erefore, a nonlinearity evaluation has the capacity to identify the variation of nonlinearity severity on the dynamics of the rotor systems.
In this paper, a nonlinearity evaluation is proposed to identify the rub-impact in rotor systems with pedestal looseness.Nonlinear mathematical models are established for rotor systems with single pedestal looseness as well as pedestal looseness coupled with rub-impact.e elastic force and the rubbing force are considered in regards to the position of pedestal looseness and rubbing.
en, the severity of the nonlinearity of the rotor systems is quantified using a nonlinearity evaluation.e definition of nonlinearity measure is introduced and a special case is used to quantify the nonlinearity of dynamics.e dynamics of rotor systems with pedestal looseness coupled with rub-impact is analyzed by integration using the adaptive fourth-order Runge-Kutta method for the mathematical models.Different looseness clearances are used in order to calculate the degree of nonlinearity of the dynamics in rotor systems.e experiments are conducted to verify the effectiveness of the proposed method.Two kinds of experiments are used to collect the vibration responses on the different setting looseness clearances.Two curves are obtained by nonlinear fitting method for the values of nonlinearity measures on different looseness clearances.en, rub-impact can be identified from the rotor systems with pedestal looseness by comparing the curves that denote the trend of the nonlinearity evaluation.

Modelling for Rotor Systems with Pedestal
Looseness Coupled with Run-Impact 2.1.Model for Rotor Systems with Pedestal Looseness.A simple rotor system with two identical oil film bearings at both sides is shown in Figure 1. e shaft sections in the system between the disk and the bearings are considered to be elastic and massless.e equivalent lumped mass in the position of the disk is m 2 , which is m 1 in the position of the bearings.It is assumed that the left support has the single pedestal looseness, the maximum static gap of the looseness is δ, and the lumped mass of the pedestal involving the looseness is m 3 .c 1 , c 2 are the equivalent damping coefficients in the positions of bearing and disk, respectively, while k is the stiffness coefficient of the shaft.e foundation, or the joint to the pedestal, is equivalent to a spring-damping system with the damping coefficients k b and c b .
It is assumed that the horizontal and vertical displacements in the right-bearing position are x 1 , y 1 , in the disk position x 2 , y 2 , in the left-bearing position x 3 , y 3 .e small movement of the left pedestal in the horizontal direction is considered to be negligible, and its displacement in the vertical direction is y 4 .e differential equations for the bearing-rotor system with single pedestal looseness can then be written as follows: where e b is the unbalance, ω is the rotating speed.
When the pedestal looseness occurs, the equivalent stiffness and the damping coefficients of the foundation, or the joint to the pedestal, can be expressed using the following piecewise linear structure [28]: where stiffness and damping actions are considered in three parts.
If the bearings are oil film bearings, 1) denote the force components of the oil film in the right bearing in the horizontal and vertical directions, while ) are that of the left bearing.e nonlinear forces from the journal bearing are obtained under the short bearing theory [33,34].For the sake of convenience, F x , F y are used to denote the components of the oil film force for two bearings in the horizontal and vertical directions, which is calculated using  Shock and Vibration where the adjustment factor s � μωRL[R/c] 2 [L/(2R)] 2 , ω is the rotating speed, R is the radius of bearing, L is the length of bearing, c is the bearing clearance, and μ is the oil viscosity.e force components f x , f y in the x and y directions can be obtained using the following equations: where V, E, S, β can be represented as If there is no pedestal looseness existed in the rotor systems, the dynamics of rotor systems can be described by the following differential equations: where denote the force components of the oil film in the right bearing in the horizontal and vertical directions, while F x (x 3 , y 3 , _ x 3 , _ y 3 ), F y (x 3 , y 3 , _ x 3 , _ y 3 ) are that of the left bearing.

Model for Rotor Systems with Pedestal Looseness Coupled
with Rub-Impact.If rub-impact occurs in the rotor system (1) with single pedestal looseness, the previously described model with identical oil film bearings with rotor-to-stator rub is shown in Figure 2. Partial rub between the rotor and the stator is taken into consideration in the rotor-bearing system with pedestal looseness in this paper.An initial clearance of δ 1 between the rotor and the stator is assumed for rub-impact.When comparing with one complete rotation period, the duration time of the rub-impact is very short.An elastic impact model is then used, neglecting the thermal effect of friction.Additionally, the Colulomb type of frictional relationship is assumed in the analysis of the dynamics.
When rub-impact occurs, as shown in Figure 3, (x, y) denotes the coordinates of the rotor center, φ is the angular displacement, and ω is the rotating speed.e radial impact force G n and the tangential rub force G r can thus be expressed as [6] where f is the friction coefficient between the rotor and the stator, k c is the radial stiffness of the stator, and is the radial displacement of the rotor.ese two forces can be written in the x − y coordinate as follows: which can be rewritten as Because sin φ � y/R, cos φ � x/R, the forces G x , G y can be expressed as follows: When comparing to the differential equation ( 1) for the rotor system with single pedestal looseness, rub forces are taken into consideration for this situation.en, the differential equation for rotor systems with pedestal looseness coupled with rub-impact can be obtained as follows: Figure 2: Bearing-rotor system with pedestal looseness coupled with rub-impact.

Shock and Vibration 3
N r (ω) : Because of the partial rub between rotor and stator, the radial impact force and the tangential rub force are considered in the equation (11).Comparing with the dynamical model ( 1), this will give more complicate nonlinear dynamics than the rotor system with pure pedestal looseness.It can be concluded that the nonlinearity effects on the dynamics of the rotor system will become more serious, which is verified by the finite element analysis [30].If the partial rub occurs in rotor systems with pure pedestal looseness, the severity of nonlinearity in rotor systems will increase significantly.Mathematically, this variation can be identified by using a nonlinearity measure to quantify the nonlinearity degree of the dynamics.

Definition of Nonlinearity Evaluation.
e evaluation of the nonlinearity represents an approach to systematically quantify the degree of nonlinearity for dynamical systems [31,32], which is also referred to as the nonlinearity measure or nonlinearity test.e fundamental idea underlying the evaluation of nonlinearity is to compare the dynamic behaviors of the nonlinear systems with a certain linear system in an appropriate setup.
e most common setup for nonlinearity measures is depicted in Figure 4 [32].
Assuming that a nonlinear stable dynamical system N has the transfer operator N : u↦y � N(u), this is compared to a linear model that has the linear transfer operator G : u↦y � G(u), which is used to approximate the behaviors of the system N. u, y, y represent the input and output trajectories of the systems N and G, respectively.Without a loss of generality, it is assumed that N[0] � 0. e error is defined as the difference between the output y of the nonlinear system N and the output y of the linear system G.
is error measurement provides information regarding how well the dynamical behaviors of a nonlinear system resemble that of a certain linear model.To evaluate the nonlinearity of the dynamics, the following definition is employed to compute the nonlinearity degree, which has been proposed by T. Schweickhardt and F. Allgower [31]: where N(u) is the output of the nonlinear system N, G(u) is the output of the linear model G, and ϕ u N denotes the value of the nonlinearity quantification using the input u.It is evident that the nonlinearity evaluation depends on the linear system G and the set of considered inputs U.
e set U typically describes the region of operation in which the nonlinearity of the system N is to be assessed.In this paper, the constant rotational speed ω is selected as the input u for systems N and G. en, the following definition is used to quantify the nonlinearity of rotor-bearing systems: e most important step in nonlinearity evaluation is to identify a proper linear model to compare its dynamic behaviors with that of the original nonlinear model.Computation of ϕ ω N in equation ( 13) can be extremely complicated, because the minimization must be performed over all linear operators G in the set G. Due to the computation requirement for engineering applications, it is appropriate to use a simple definition of nonlinearity quantification.In some situations [35,36], Taylor expansions are used to obtain the linear approximations on the vicinity of the equilibrium position.Shock and Vibration Linearizing the nonlinear terms on the static equilibrium position allows for a linear system L(ω) to be derived to calculate the dynamical differences between the outputs.e details of the procedure for Taylor expansions of nonlinear terms are shown in Appendix A. As a special case of linear approximated system, then e equation ( 14) is a special case of the definition of nonlinearity measure, which is used to assess the pedestal looseness of rotor systems in References [35,36].However, different linear model will influence on the nonlinearity quantification on a real rotor system by using the directly measured data.To improve the applicability of the method, the model N 0 (ω) in equation ( 6) for rotor system with no pedestal looseness is introduced in equation ( 14) to obtain an upper bound to measure the severity of nonlinearity for dynamics: where d � ‖N 0 (ω) − L(ω)‖/‖N(ω)‖, and it can be computed according to the model ( 6) and its linear approximation by Taylor expansion.en, the following definition of nonlinearity measure is used to quantify the nonlinearity of rotor systems (11): Similarly, the following definition ( 17) is used to quantify the nonlinearity of rotor systems (1).
In the above equations, ‖(•)‖ is the norm, with the following definition: where x(t) denotes the nonlinear vibration signals that describe the dynamical outputs of the dynamic systems.e energy contained in the signal x(t) over the infinite time interval t ≥ 0 is described by equation (18).However, in principle, any norm can be used for the following considerations.e evaluation of norm in equation (18) requires the following numerical approximation: where t i , i � 1, 2, . . ., t nim   are equally spaced mesh points in the interval [0, T max ], and Δt is the average sampling time.e maximum integration time T max can be considered as a characteristic time of the system, based on a certain degree of physical knowledge of rotor systems.
e result of the nonlinearity evaluation ( 16) and ( 17) will satisfy ϕ p , ϕ r ≥ d.Discussing the nonlinear equations of motions in an analytical way is quite challenging; therefore, numerical methods have to be resorted for the calculations of the nonlinearity evaluation.e parameter values used for the nonlinearity measure and simulations are given in Table 1 of Appendix B, while the values of the parameters for oil film bearings and rub-impact used in the numerical examples are given in Tables 2 and 3, respectively.Some parameters used in the simulations are chosen as the same as Tables 1 and 2 in References [35,36].
e adaptive fourth-order Runge-Kutta method is used to integrate the dynamic systems, while a smaller marching step is chosen to ensure a stable solution and to avoid the numerical divergence at the point where the damping and stiffness parameters are discontinuous.To maintain the accuracy of integrations, nondimensionalization is used to transform the systems into new ones, which can be implemented by using where x 2 , y 2 , x 3 , y 3 ], c is the average thickness of oil film, and ω is the rotating speed of the bearing-rotor systems.

Nonlinearity Evaluations for the Dynamics of Rotor
Systems.A nonlinear vibrating system with seven degrees of freedom and piecewise linear stiffness and damping is represented by equation (11).In comparison to equation ( 1), the difference is that there are two rub forces in the position of the disk.e force components of the oil film and the rubimpact forces have nonlinear characteristics.To calculate the value of d in equation (11), Taylor expansion is used to Shock and Vibration 5 linearize these nonlinear terms of equation ( 6) on the static equilibrium position, and a linear model ( 21) can be developed as follows: L(ω) : , and  F y (x 3 , y 3 , _ x 3 , _ y 3 ) denote the linear approximations of the oil firm forces of bearings.
e calculations for the procedure are given in (A.1) and (A.2) of Appendix A. e values for the parameters used in the nonlinearity measure and the simulations for equations ( 6), (11), and ( 21) are given in Appendix B. In order to maintain the accuracy of integration by fourth-order Runge-Kutta method, nondimensionalization is used to transform these systems to new systems using equation (20).Looseness clearances of the pedestal is used as a control parameter in order to perform a detailed investigation on the nonlinear dynamics of the bearing-rotor system, which uses the values 0 to 0.0035 m for simulations.
e vibration amplitudes at the positions of pedestal looseness and disk are collected for investigations; the amplitude spectrums in the frequency domain are also given for analysis purposes.
In the process of nonlinearity evaluation for pedestal looseness coupled with rub-impact, the practical process of generation for rub-impact in the rotor system is simulated.A critical point of looseness clearance is chosen for the generation of rub-impact, and this process is divided into two stages.When the looseness clearance in the pedestal is smaller than the critical point, only single pedestal looseness without rub-impact occurs in the rotor system.Otherwise, the rotor system will have pedestal looseness coupled with rub-impact.Assuming that the critical point of looseness clearance is set to be 0.0002 m, rubbing between the rotor and the stator starts to occur in the rotor system.When the looseness clearance is smaller than 0.0002 m, the equations for rotor systems with single pedestal looseness are used to investigate the dynamic behaviors and for the nonlinearity evaluation.e vibration amplitudes and its spectrums for a bearing-rotor system at the positions of pedestal and disk with a looseness clearance of 0.00002 m are shown in Figures 5 and 6, while the orbit at the position of pedestal looseness is shown in Figure 7.
When the looseness clearance is greater than 0.0002 m, the rotor system will have a pedestal looseness coupled with rub-impact.e vibration amplitudes and its spectrums for a bearing-rotor system at the positions of pedestal and disk with a looseness clearance of 0.0002 m are shown in Figures 8 and 9, respectively, while the orbit at the positions of pedestal looseness is shown in Figure 10.
e vibration amplitudes and its spectrums for a bearing-rotor system at the positions of pedestal and disk with a looseness clearance of 0.002 m are shown in Figures 11 and 12, respectively, while the orbit at the positions of pedestal looseness is shown in Figure 13.
It can be seen in Figures 8 and 11 that the vibration amplitudes will decrease with a reduction in looseness clearances; the responses exhibit more harmonic components with a decreasing in looseness clearances.
e vibration amplitudes and the orbits of a rotor system exhibit increasingly nonlinear dynamics with an increase in looseness clearances.e match of looseness clearances and the nonlinearity evaluation is shown in Figure 14, where dot ( * ) denotes the results of nonlinearity quantification at a given looseness clearance.A nonlinear fitting method with exponential functions, using the same time interval, is used to obtain the trend for the nonlinearity evaluation with the increasing looseness clearances.It can also be found that the values of nonlinearity measure in Figure 14 are larger than d, which is calculated according the model of the rotor system with no pedestal looseness.

Experiments and Verifications
In this section, experiments are conducted to collect the vibration signals of a rotor system for nonlinearity evaluation and demonstrate the effect of the proposed method.
e motions of the rotor system are described by the analytical model (1), which the values of the parameters for simulations are given in Tables 1-3 of Appendix B. Some parameters used in simulations are kept the same within experiments, while the others are obtained from the references (using the similar experiment setup to collect the vibration signals).
e analytical models show the increasingly nonlinear dynamics of the rotor system with an increase in looseness clearances, which can guide the implementation of experiments and data analysis.Due to the restrictions of experimental conditions, less number of looseness clearances can be used to conduct the process of signals collection, dynamics analysis, and nonlinearity evaluation.It is shown that the trend of nonlinearity evaluation obtained from experiments remains to be largely consistent with that from simulations.Two kinds of experiments are conducted for the rotor system with pure pedestal looseness and with pedestal looseness coupled with rub-impact.e trend curves with the increasing looseness clearances show that the rubimpact will have important in uences on the nonlinearity degrees of the dynamics, which can be used to identify the existence of rub-impact in rotor systems with pure pedestal looseness.

Shock and Vibration
As shown in Figure 15, the experimental setup that produced by SpectraQuest company consists of a rigid cylindrical shaft supported by two cylindrical journal bearings (an oil lm journal bearing at the left end near the driving motor, and an oil lm journal bearing in the position of pedestal looseness at the right end).e signal collector is B&K Pulse.ere are ve accelerometers in the rig, which positions are supporters in the horizontal and vertical directions on the position of pedestal looseness, supporter in the horizontal and vertical directions at the left end, and on the position of fastening bolt.
e experimental setup used in this paper is the same as in Reference [35].e only di erence in this paper is that no rubber gasket is used for the fastening bolt in the position of pedestal looseness, while a rubber gasket is used in Reference [35].8

Shock and Vibration
In the experiment mentioned above, the rotor speed is set to be 2100 rpm and the sample frequency is 3.2 kHz.e pedestal looseness is simulated by adjusting the fastening bolt ⑧ according to the dial ⑨.
e distance between threads is measured in advance, and the bolt can be adjusted according to the scale on the dial.e looseness clearances will change 1 × 10 −4 accordingly, when the M10 bolt is adjusted 36 °.With the same conditions, the experiments for vibration signals acquisition are conducted into three stages.First, the vibration signals for the rotor system without pedestal looseness are collected as the basis of nonlinearity measure, and d in equations ( 16) and ( 17) is computed numerically.Second, the experiments for the rotor system with single pedestal looseness are conducted to collect the vibration signals with di erent looseness clearances.Lastly, the experiments for rotor systems with pedestal looseness  Because the results of nonlinearity evaluation should be obtained at every setting looseness clearance for two kinds of experiments, it needs a huge amount of computations.For simpli cation, only vibration signals and their amplitude spectrums for experiments of rotor systems without pedestal looseness and pedestal looseness coupled with rubbing can be exhibited here.In the rst stage of the experiment, vibration signals and their amplitude spectrums of the supporter in the vertical direction for rotor system without pedestal looseness are given in Figure 17.
In the experiments for the rotor system with pedestal looseness coupled with rub-impact, Figures 18 and 19 are vibration signals and their amplitude spectrums of the supporter in the vertical direction on the position of pedestal looseness, where the looseness clearances are set to be 0.0002 m and 0.002 m.
In order to compare this scenario to the scenario of single pedestal looseness, the same number of looseness clearances with equal intervals is used in order to calculate the dynamical behaviors of this rotor system, and the nonlinearity evaluation    16) and ( 17), and the inner product is calculated using the equation (19).Two curves are separately obtained using the nonlinear tting method with exponential functions according to the values of the nonlinearity measure, which are shown in Figure 20.It can be found that the results of the nonlinearity evaluation are all equal or larger than d.
Each curve displays the trend for the results of the nonlinearity evaluation on di erent looseness clearances.When the looseness clearance is set to values greater than 0.0002 m, the rub-impact is taken into consideration in the simulations, and the trend and the amplitude of the curve begins to change.Additionally, it is obvious that when small values are chosen for the looseness clearance, the results of the nonlinearity measure have a slow increasing trend.If   values larger than the threshold are chosen for the looseness clearance, the nonlinearity evaluation of the rotor system with pedestal looseness coupled with rub-impact increases more dramatically than that of the rotor system with single pedestal looseness.It comes to the same conclusions with the numerical example.
On the other hand, the curve that denotes the trend of the nonlinearity evaluation for the rotor system with single pedestal looseness is close to a straight line, although it is also given by a tting method with exponential functions.However, the curve that denotes the trend of the nonlinearity evaluation for the rotor system with pedestal looseness coupled with rub-impact displays more characteristics other than a straight line.As shown in Figure 20, the values of nonlinearity measure for the rotor system with pedestal looseness coupled with rub-impact is obviously lager than that with pure pedestal looseness on the same looseness clearances.

Conclusions
In this paper, a nonlinearity evaluation was proposed to identify the rub-impact in rotor systems with pedestal looseness.First, nonlinear mathematical models were established for bearing-rotor systems with single pedestal looseness and pedestal looseness coupled with rubimpact.e elastic force and rubbing force on the positions of pedestal looseness and rubbing points were taken into consideration.Secondly, the de nition of nonlinearity evaluation is given, and a special case is applied to quantify the nonlinearity of dynamics for rotor systems.
is evaluation is computed numerically based on the dynamical models on di erent looseness clearances.Finally, experiments are conducted to verify the e ectiveness of the proposed method.
e experiments for the rotor systems with pure pedestal looseness and pedestal looseness coupled with rub-impact are used to collect the vibration signals on di erent looseness clearances.Two di erent curves were obtained separately using a nonlinear tting method with exponential functions, according to the values of the nonlinearity measure.It can be observed that if the rotor systems with pedestal looseness have a rub-impact between the stator and the rotor, the nonlinearity degree will have a larger value than rotor systems with pure pedestal looseness.
is will provide a new idea to identify the rub-impact in rotor systems with pedestal looseness.

Figure 4 :
Figure 4: Setup for the comparison of a nonlinear system with a linear system.

Figure 5 :
Figure 5: Vertical vibration amplitude signals and its frequency diagram on left bearing.

Figure 6 :
Figure 6: Vertical vibration amplitude signals and its frequency diagram on the disk.

Figure 9 :Figure 10 : 4 Figure 8 :
Figure 9: Vertical vibration amplitude signals and its frequency diagram on the disk.

Figure 12 :Figure 11 :
Figure 12: Vertical vibration amplitude signals and its frequency diagram on the disk.

Figure 14 :
Figure 14: e relationship of results of nonlinearity measure and looseness clearances.

Figure 17 :
Figure 17: Vertical vibration amplitude signals and its frequency diagram on bearing without pedestal looseness.(a) Vibration amplitude signals and its frequency diagram.(b) Denosing vibration amplitude signals and its frequency diagram.

Figure 16 :
Figure 16: e schematic diagram of rubbing institutions.

Figure 18 :
Figure 18: Vertical vibration amplitude signals and their frequency diagram on bearing with pedestal looseness coupled with rub-impact (looseness clearance 0.0002 m).(a) Vibration amplitude signals and its frequency diagram.(b) Denosing vibration amplitude signals and its frequency diagram.

Figure 19 :
Figure 19: Vertical vibration amplitude signals and their frequency diagram on bearing with pedestal looseness coupled with rub-impact (looseness clearance 0.002 m).(a) Vibration amplitude signals and its frequency diagram.(b) Denosing vibration amplitude signals and its frequency diagram.

Table 1 :
Values of the simulation parameters.

Table 2 :
Values of the parameters for oil film bearing.

Table 3 :
Values of the parameters for rub-impact.