Analytical study of the spherical hydrostatic bearing dynamics through a unique technique

Because of their self-alignment property and design simplicity, the spherical hydrostatic bearings have the advantage over the other bearing configurations. Their static and dynamic performances have been intensively studied. Focusing on the bearing dynamic performance, it could be realized that the researchers used to mechanically excite the bearing in the experimental studies and perturb the rotor (spatial finite displacement) in the theoretical studies, observing its behavior and expressing it by dynamic stiffness and damping coefficient. Owing to a lack of information on bearing oscillation, this study adopts a new method to analyze this bearing behavior theoretically and grasp its nature. Unusually, the bearing vibration is studied hydro-dynamically rather than mechanically, showing the effect of eccentricity, inertia, restrictor type, seat configuration, and the supply pressure on the performance. New and unique formulas have been derived to predict the frequency, stiffness, and dampness, explaining how these parameters create the bearing self-alignment property. Some thrust bearings, designed with different techniques, have been examined where the results show vibration even in the bearing’s stationary state. Also, the newly derived formulas could precisely predict the rotation and rotational speed of the Kugel ball (recently proved and designed as a spherical hydrostatic bearing).


Analytical study of the spherical hydrostatic bearing dynamics through a unique technique Ahmad Waguih Elescandarany
Because of their self-alignment property and design simplicity, the spherical hydrostatic bearings have the advantage over the other bearing configurations.Their static and dynamic performances have been intensively studied.Focusing on the bearing dynamic performance, it could be realized that the researchers used to mechanically excite the bearing in the experimental studies and perturb the rotor (spatial finite displacement) in the theoretical studies, observing its behavior and expressing it by dynamic stiffness and damping coefficient.Owing to a lack of information on bearing oscillation, this study adopts a new method to analyze this bearing behavior theoretically and grasp its nature.Unusually, the bearing vibration is studied hydro-dynamically rather than mechanically, showing the effect of eccentricity, inertia, restrictor type, seat configuration, and the supply pressure on the performance.New and unique formulas have been derived to predict the frequency, stiffness, and dampness, explaining how these parameters create the bearing self-alignment property.Some thrust bearings, designed with different techniques, have been examined where the results show vibration even in the bearing's stationary state.Also, the newly derived formulas could precisely predict the rotation and rotational speed of the Kugel ball (recently proved and designed as a spherical hydrostatic bearing).

List of symbols A
6µ q e 3 π p i = Q A s Bearing projected area ( A w sin 2 θ i ) A w Wet projected area ( πR Bearings constitute a fundamental component of any rotating machine, where a good understanding of their behavior allows the prediction of the system function.High speeds and high power with reliable operation are necessary.Accurate prediction and control of bearing behavior is required.It is commonly known that more than 40% of machine failures are the result of bearing problems.The most efficient approach to solving the problem is through analysis conducted in numerous studies, which leads to better performance.The potential characteristics of low running friction, small viscous dissipation, high load carrying capacity, and high stiffness have made hydrostatic bearings indispensable for industrial applications in present-day heavy-duty and highspeed devices such as machine tools, precision measuring instruments, hydraulic piston pumps, motors, telescopes, gyroscope gimbals, dynamometers, radar tracking units, and craft engines.Spherical bearings have advantages over other configurations due to their self-alignment property and their capacity to accept radial and thrust loads.Therefore, the operation of the spherical bearings is not affected by angular misalignments.Mayer 1 and Elescandarany 2 derived the Reynolds governing equation for spherical bearing applications.Kazama 3 presented an optimal design for the types of fitted and clearance thrust spherical bearings with capillary and orifice restrictors.This design was based on minimizing power losses and maximizing stiffness.The study concluded that the central pressure ratio is two-thirds for the fitted type and 0.69 for the clearance type.Yacout [4][5][6][7][8] analytically studied both types after developing the Reynolds equation to apply it to the thrust spherical bearing, considering the effects of surface roughness, centripetal inertia, viscosity, and seat configuration on bearing performance.The central pressure ratio was calculated to be two-thirds for both types, leading to consistent bearing performance.Examples of optimal designs based on minimum friction and flow rate were provided.Elescandarany [9][10][11] developed two new techniques to design a fitted type, with and without a restrictor.Based on the steady and unsteady states, Gosh et al. [12][13][14] and San Ander 15 numerically studied the dynamic behavior of the journal bearing.This behavior is represented through stiffness and damping coefficients by exciting the bearing journal while also considering the fluid compressibility in the bearing recess.Applying the first-order perturbation method, Yacout 16,17 studied the dynamic behavior (dynamic stiffness and damping coefficient) of the thrust spherical bearing, considering the compressibility of the fluid in the bearing recess.Sharma et al. 18,19 numerically studied thrust and journal spherical bearings with Newtonian and non-Newtonian fluids, considering the effect of surface roughness on bearing performance.Gupta and Kumar 20,21 examined the behavior of a hydrodynamic squeeze film between a stationary spherical surface and a hemispherical bearing.The research focused on understanding the impact of surface roughness on bearing performance when subjected to a constant load.Based on self-excited vibration, the studies generally are concerned with techniques to control and suppress the vibration.
Cameron 22 stated that one of the most troublesome features of high-speed bearings could be oil whirl.It is a vibration that occurs slightly below half the speed of the shaft and is usually found during acceptance trials, which is the most awkward time for serious faults to appear.The main features of the problem were analyzed and understood a few years ago.There are now a few guidelines for the solution to the problem.A bearing that completely suppresses oil whirl has not yet been produced, though new anti-whirl bearing configurations are patented yearly.Rowe 23 addressed the whirl phenomenon, stating that the puzzling question is the following: How is it possible for a heavily damped system to act as though it has zero dampness, and replying that the explanation is not altogether simple and lies like the hydrodynamic stiffness causing a positive reaction force in the direction of the velocity, in turn, this force acts in opposition to the squeeze force and cancels out the dampness.Jinhao et al. 24 presented an experimental study on the active control of self-excited vibrations in a rotor bearing system supported on a pair of externally pressurized thrust bearings in the axial direction and on actively controlled journal gas bearings in the radial direction.Feedback control systems were constructed with gap sensors to measure the vibration of the piezoelectric actuators embedded in the rotor and thee PID (proportional-integral-derivative) controllers.The experimental results show that self-excited vibration could be effectively suppressed with the designed feedback control system if the gains of the PID controllers are appropriately tuned.An experimental study offered by Talukder and Stowell 25 was made of the occurrence of a pneumatic hammer in an externally pressurized air journal with admission into double plane.Various parameters that affect the onset of the pneumatic hammer were investigated, such as the depth of the recess, the diameter of the orifice and the bearing mass, stating that the pneumatic hammer is relatively easily avoided in the journal bearings.Taiwari et al. 26 presented research reviewing the bearings' dynamic parameters in rotating machines.Attention is given to vibration-based methods where the review covers descriptions of experimental measurement techniques, mathematical modeling, parameter extraction algorithms, and uncertainty in the estimates applied to a variety of bearings and suggests some articles as the need for more experimental work in the field of rotor dynamics to study the bearing and support influence upon the rotor response, especially for full-scale rotor systems.Ales Tondl 27 presented a survey of publications dealing with vibration suppression methods in different types of self-excitation.The attention is given to passive and active means using parametric excitation, focusing on the latter, and the necessary steps for answering the main questions are formulated to initiate further investigation, especially application to actual systems.Feng and Hahn 28 presented a study on circular and elliptical bearings to examine and compare their impact on the rotor support of industrial machines.The findings showed significant variations in natural frequencies and critical speeds depending on the type of bearing used.Charki et al. 29 provided a numerical simulation and an experimental study to assess the stiffness and dampness characteristics of thrust air bearings with multiple orifices using finite element modeling to solve the Reynolds equation.
The results obtained numerically showed that performance is related to bearing design, and the experimental investigation allowed analyzing bearing behavior.Peiji et al. 30 studied the impacts of cylindrical bearings and four-lobe bearings on the vibration of the rotor system, combining theoretical analysis and experimental results, concluding that different types of bearing have effects on the critical characteristics of the rotor system and the vibration of the rotor supported by the cylindrical bearing is more stable than the four-lobe bearing, where the results showed good agreement between theoretical and experimental analysis.Zhifeng et al. 31 presented a review of research articles related to the introduction of developments in hydrostatic bearings (the basic theory contains equations and analysis methods including analytic, numerical, and experimental methods) where typical applications were based on rectangular oil pads, circular oil pads, and journal bearings.Polach et al. 32 presented a simulation for a system consisting of a rotor supported on two journal bearings and a controller, which controls the speed of the rotor, to study self-excited and flow-induced vibrations of the rotor, reporting that subsynchronous flow-induced vibrations can occur only if the rotor is well balanced and is radially supported on poorly lubricated journal bearings, whose supply bores are located in the lower bearings half.Davorka and Marko 33 presented a modern method for testing the condition of bearings in operation using vibration analysis and demonstrating the importance of this method by changing the parameters related to two different bearing structures grouped in the so-called Sommerfeld number.Dmytro et al. 34 presented a study to improve machining accuracy in precision cutting and grinding by designing a developed spindle with water-lubricated hydrostatic bearings.Dynamic characteristics such as stiffness, dampness, and frequency of the machine spindle were investigated by simulation and experiments, where it is reported that the simulation results were in good agreement with the actual spindle dynamics obtained experimentally.Yingjie Wang et al. 35 studied the characteristics of hydrostatic bearings by analyzing load capacity, stiffness, and dampness using CFD and 3D meshing, hence Compared to traditional empirical formulae, 3D meshing and CFD analysis were more accurate and could be adapted to the analysis of various shapes of hydrostatic bearings.Runlin Chen et al. 36 studied a motorized hydrostatic spindle system of a grinding machine using a two-degree-of-freedom stiffness model, considering the rotation of the rotor, where stiffness was analyzed under different excitation frequencies, and the vibration test in the spindle system was implemented through the hammering method.The stiffness of the spindle system was found to agree well with the theoretical calculation, with an average error of about 14%.Xin Qin et al. 37 presented a novel hydrostatic squeeze film metal mesh journal bearing, which uses both a hydrostatic squeeze film damper and a metal mesh damper to suppress vibration of rotor bearing systems where lubrication equations were introduced to calculate the dynamic characteristics.Experiments were conducted to study the reduction in vibration of the bearing rotor system with the bearing; the theoretical and experimental results showed that the bearings exhibit excellent dampness and vibration attenuation characteristics.Despite the attention given to the journal bearing through numerous studies handling its static and dynamic behavior, in addition to the recent studies on its self-excitation vibration and endless research on eliminating, suppressing, and controlling this undesired disturbance by active or passive means, nothing has been offered about the vibratory behavior or/and the self-alignment property of the spherical hydrostatic type of bearings.Therefore, as the next step to designing an ideal model of such a bearing, this paper presents a new study investigating the natural properties of this bearing by applying a new method depending on the Reynolds equation and the hydrodynamics of the bearing.

Derivation of mathematical expressions
Arranging and integrating Eq. (2) give: (1) Relation between ( − θ ) and the frequency (f) Equation (1b) could be put in the form: Hence: Then:

The characteristic equation solution
Hint: It would be useful to replace (e) with (E) in the equation to prevent the confusion between (e) the eccentricity and (e) the exponential number and after the solution it would be returned to its original form.
Hence it would be: Then, the solution is: Using the boundary conditions: Returning the equation to its original form will be:

Bearing pressure
Arranging Eq. (3) gives: The integration could be found in the appendix as: (21) www.nature.com/scientificreports/

Bearing load capacity
From 4 : The integration of Eq. ( 25) could be found in the appendix: The bearing dynamic stiffness From Eqs. (23, 26), the dynamic stiffness could be calculated as (Fig. 2):

The bearing dampness
The dampness could be calculated as: (24)

Bearing vibration inspection
Rotating machines must be subjected to vibration inspection to evaluate condition and efficiency where defects could be revealed.

The vibration measurement
Vibration could be measured using one of the following measurement techniques: displacement, velocity or acceleration.

The measurement position
The probe position on the bearing depends on the measuring technique.

The vibration remedy
The vibration could be suppressed via passive or active means.

Results
The goal of this paper, as stated previously, is to investigate the vibratory behavior of the fitted thrust spherical hydrostatic bearing.Applying a new method, new and unique formulas have been derived to predict the vibratory bearing characteristics such as the vibration frequency, the stiffness, the dampness etc.The effect of eccentricity, supply pressure, restrictor type and its dimensions, inertia, and seat configuration on the bearing characteristics is studied.Seven cases have been checked to assess the new method's validity.Figures

Discussion
The conception of this research arises from the author's references 6-8, 16, 17 which point to the vibratory phenomenon of this type of bearings promising to be studied in detail in future separate research.However, a new method is used to study the bearing vibratory behavior avoiding the assumptions of infinitesimal oscillation or disturbance or/and applying the first-order perturbation method as done before.The new method adopted is mainly dependent on bearing hydrodynamics rather than reliance on mechanical analysis as in 2,11 .

The new derived formulas
Pressure and load capacity New pressure and load capacity Eqs.(24, 26) have been developed [Appendix A1, A2]. Figure 3a depicts the pressure and load distributions that coincide with that obtained from 4 and the pressure has the same feature as in 19 despite the use of non-Newtonian lubricant.The wave parameters Figure 3b, c, d show the parameters distribution and their mean values.Figure 3b, c express the 1st and 2ndtime derivatives of the variable parameter (ꝋ).The 1st derivative is the angular wave velocity, while the 2nd is the angular wave acceleration.Figure 3d is the superimposing of the angular speed on acceleration, which gives an idea about their relationship, which may not be easily realized from the mathematical equations.The optical examination of Fig. 3d provides a direct impression of the mathematical nature of both curves.It could be mathematically realized that the angular acceleration is a form of the sine function, so the angular velocity is a cosine function (as its time integration) lagged by ninety degrees.Hence, it should be predicted that the displacement (e) in the spatial domain must be in a sinusoidal form (as an integration of the velocity).Their maximum, mean, and minimum values are used in the mathematical calculations, while their positions on the seat could be helpful for industrial inspection.

The axial wave velocity
The axial wave velocity is represented by Eq. ( 17), and its distribution is shown in Fig. 4a.The velocity nonlinearity is revealed by its distribution.The nonlinearity of the velocity produces two equal averages at two different positions on the bearing seat, which could be interpreted mathematically as inconsistency, which exists from the energy difference (Kinetic energy) between these two positions.A mechanical moment is produced by the difference between these two kinetic energy values.This moment causes the bearing rotor to rotate if it has zero degrees  of freedom, as in the case of the Kugel ball 2,11 , and to oscillate if it has one degree of freedom, as in the case of the current bearing. 2,11 alculated this moment mechanically by force analysis and determined the frequency.The optical examination of the curve reveals that it is of a parabolic form (two equal roots at two different positions), which mathematically means that this axial velocity is a 2nd-degree function.Hence, it should be expected (mathematically) that the displacement (e) in the time domain must be a 3rd-degree function.
The wavelength Equation (18) connects Eqs.(12, 17), yielding the wavelength.Figure 4b depicts its distribution, as well as its mean value and position.

The wave frequency
The unique frequency Eq. ( 12) calculates and predicts the wave frequency (vibration), while Fig. 4c shows the frequency distribution and the mean value, and its position.

The wave power
Equation (19) calculates the wave power per unit wavelength to be (0.2517 N/s).The wave power is determined to provide only an idea about the energy which could be lost.The displacement characteristic equation The unique characteristic Eq. ( 22) of the rotor's displacement represents mathematically a partial differential equation of an oscillating damped motion.This demonstrates that the bearing eccentricity oscillates, resulting in rotor vibration in the case of a stationary bearing seat (without a damper).

The (PDE) solution
Equation ( 23) is the mathematical solution to the partial differential Eq. ( 22).The solution demonstrates that displacement is both a time and a spatial function.Figure 5a depicts the synchronous function distributions in both domains.Figure 5b shows the time domain function distributions (at various spatial values) and the distribution of their mean values.Figure 5c shows the spatial domain function distributions (at different time values) and the distribution of their mean values.Figure 5d depicts the superimposing of the mean value distributions in the two domains and the average of each distribution representing the function mean values in the two domains.From the optical examination, it could be realized that Fig. 5b is a 3rd-degree function as predicted before.The function's continuous increase in the time domain leads to the unbalance.Figure 5c looks sinusoidal with severe dampness.This suppresses the function inconsistency in the spatial domain.Hence the function behavior removes the mystery of the resistance to angular displacement property of this type of bearing.The superimposing of the function components in the two domains, Fig. 5d, shows two equal averages at two different positions.This leads to the bearing inconsistency, as was explained previously, by creating a moment  that results in rotating or oscillating the bearing rotor depending on its degree of freedom.Hence it could be confirmed that this generated moment is a result exerted due to the potential energy difference between the two averages of the displacement (e) or/and the kinetic energy difference between the two averages of the axial wave velocity (e o ) as declared before.
The bearing stiffness Equation (27) relates between Eqs. (23, 26) existing the bearing stiffness.The bearing is represented by a loaded vibrating spring in Fig. 2. The stiffness is also a time and spatial function, and Fig. 6a-d represent its distribution as done with displacement (e). Figure 6a shows the synchronous stiffness function distributions in both domains.Figure 6b, c shows that the stiffness in the time domain decreases while increasing in the spatial domain.This makes the bearing more stiffened at the lubricant exit than at the inlet.Figure 6d shows two equal averages of the stiffness function at different times.The difference between the two stiffness averages in the time domain creates internal dampness (stiffness × time = dampness).This internal dampness reveals the mystery of the self-alignment property characterizing this type of bearing.

The bearing dampness
Equation (28) relates Eqs.(17, 26) existing the bearing dampness.Figure 7 shows its distribution where it reaches its maximum at the lubricant exit, which enhances the bearing self-alignment property.The effects of these factors have been applied to the main bearing parameters, the displacement, the frequency, the stiffness, and the dampness.Equations (6, 12, 23, 27, 28) are the main equations that govern this behavior.To ease the grasp of this complicated behavior, attention should be paid to the simple and basic formula used in deriving Eqs.(6).

The eccentricity effect on the behavior
As the eccentricity increases ( θ e ) decreases, in turn, the time (t) must also reduce to keep the (r.h.s) of the formula constant.Hence, the periodic time (T) automatically decreases due to Eq. ( 6).In turn, the frequency increases, Fig. 8b.Looking at Eq. ( 23), it could be realized, despite its complexity that its exponents (K1, K2) increase due to the frequency increase leading to increasing the displacement despite the time reduction Fig. 8a.Also, the eccentricity increase leads to a decrease in the load.In turn, load and displacement increase leads to an increase in the stiffness due to Eq. ( 27) Fig. 8c.
The reduction in time causes increasing the axial velocity (e o ).Hence, the decreased load and the increased axial velocity cause decreasing dampness due to Eq. ( 28), Fig. 8c.
The eccentricity increase severely harms the bearing where the vibration increases and the stiffness and dampness decrease.The configuration effect on the behavior Figure 9 a-d show the bearing behavior under different configurations, where the seat arc length is expressed by ( η ).Referring to the basic formula and the bearing configuration in Fig. 1, it could be seen that the increase in ( η ) increases ( θ e ), in turn, the time (t) to balance the (r.h.s.) of the formula, subsequently the periodic time (T) increases leading to decrease the frequency as could be shown in Fig. 9b.With the increase in (η), the load increases, and with the time increase, the axial velocity (e o ) decreases; hence from Eq. ( 28), the dampness increases as in Fig. 9d.With the frequency and time increase, the exponents of Eq. ( 23) increase leading to the displacement increase, Fig. 9a.With a quick look, the stiffness Fig. 9c shows a contradiction between Eq. ( 23) and   the results.This untrue contradiction arises from the relation between the load and the displacement, where both load and displacement are increasing.However, this means that the stiffness parameter of each case increases as the seat arc length increase, and it is a fact without any doubt about it.The seat arc length increase is crucial for improving bearing behavior.This results in decreased vibration and increased dampness and stiffness.
The rotor speed effect on the behavior Figure 10 a-d show the bearing behavior under different rotational speeds.The rotational speed (N) results in existing inertia, which increases the load 4 .This case could be treated as the previous item.Like in the case of ( η ), the rotational speed (N) increases the load, which means that the (r.h.s) of the basic formula increases.To balance this increase, the (l.h.s) should also increase; hence, because the ( θ e ) is constant, the time (t) should be decreased; in turn, the periodic time (T) decreases, resulting in the frequency increase, Fig. 10b.With the decrease in the time (t), the axial velocity (e o ) increases; hence, the dampness decreases as predicted by Eq. ( 28), where both of load and the axial velocity increase, Fig. 10d.
Figure 10a shows the decrease in the displacement due to the time decrease despite the frequency increase, where the complexity of Eq. ( 23) doesn't allow the behavior prediction without calculation.However, in Fig. 10c, the increase in both load and displacement causes the stiffness increase, as explained before.Two important finds were observed; the 1st is that the rotational inertia could severely harm the bearing where the increase in (N) leads to an increase in the frequency and decrease in the dampness, despite the stiffness increase, while the second observation is that the bearing generates vibration in the stationary state (N = 0), this vibration is recognized from the figures and the mathematical calculations confirm it.What is worthy is the ability to internally control and fix the inertia effect by designing the bearing with certain conditions 4 .
The supply pressure effect on the behavior Figure 11a-d show the bearing behavior under different supply pressures.
Figure 11b shows the unpleasant effect of the supply pressure on the bearing frequency.It could be observed from Eqs. (13, 14).
Figure 11c, d show the pleasant effect of the supply pressure on the bearing stiffness and dampness, respectively.It could be realized from Eqs. (27, 28), where the increased pressure increases the load, and at fixed eccentricity (e) and axial velocity (e o ), the stiffness and dampness increase.
Figure 11a shows the undesired effect of the supply pressure on the bearing displacement.The complicated relation of displacement with pressure could not be easily predicted without calculation, as seen from Eq. ( 23).
Hence, despite improving the stiffness and dampness, the supply pressure increases the rotor vibration and displacement.Therefore the designer could be guided to compromise between the supply pressure gain and drawback.Figure 12a, b show the frequency distribution using the capillary and the orifice at different diameter ratios, the Increasing diameters of the restrictors negatively affect the bearing vibration, and Eqs.(13, 14) predict the frequency increase at increasing the restrictor diameter.

The restrictor effect on the behavior
Figure 12c compares the mean frequencies of the capillary and the orifice when varying their diameters; it is clear that at a specific dimension, it could be possible to get the equivalent restrictor which replaces each other.
Figure 12d shows the effect of the capillary length on the frequency, where the increased length highly decreases the frequency, which could be predicted by Eq. (13).Another observation that could be noted from Fig. 12a-c is that the self-restriction bearing 10 , where (m o = 1), generates relatively high vibration.

Tested examples
Seven cases have been tested.The first two cases were selected from 4,5 , respectively, and the rest were chosen from 2 .

The first two examples
Figures 13 and 14 show that the two bearings have good performances where the frequency and the stiffness are good while the dampness is excellent.

The rest of the examples
The five cases are significant because their vibration data could be directly compared with the present ones predicted by the new technique.Figures 15b and 19b and Table 1 below show the coincidence between the results despite their completely different scientific basics.This confirms the validity and reliability of the previous and present techniques.
Despite their different dimensions, the rest of the Figs.15, 16, 17, 18 and 19 describing the hydrosphere (Kugel ball) show excellent characteristics of high stiffness, high dampness, and low displacement.It may be thought that these characteristics are not significant to be presented here, and the author finds it necessary to get more information about this Kugel ball as a type of hydrostatic thrust spherical bearing.

Figure 5 .
Figure 5. Bearing vibratory factors distribution.(a) The normalized eccentricity as a function of (ꝋ) and (t) domains, (b) The eccentricity in the time domain at different (ꝋ), (c) The eccentricity in the (ꝋ) domain at different time (t), (d) Mean eccentricity in (ꝋ) & (t) domains.

Figure 6 .
Figure 6.Bearing vibratory factors distribution.(a) The stiffness as a function of (ꝋ) and (t) domains, (b) The stiffness in the time domain at different (ꝋ), (c) The stiffness in the (ꝋ) domain at different time (t), (d) The mean stiffness in (ꝋ) & (t) domains.

Figure 8 .
Figure 8. Eccentricity effect on parameters.(a) The mean displacement, (b) The mean frequency, (c) The mean stiffness, (d) The dampness.

Figure 9 .
Figure 9.The configuration effect on parameter.(a) The mean displacement, (b) The mean frequency, (c) The mean stiffness, (d) The dampness.

Figure 10 .
Figure 10.Rotor speed effect on parameters.(a) The mean displacement, (b) The mean frequency, (c) The mean stiffness, (d) The dampness.

Figure 11 .
Figure 11.Supply pressure effect on parameter.(a) The mean displacement, (b) The mean frequency, (c) The mean stiffness, (d) The dampness.

Figure 12 .
Figure 12.The restrictors effect on the frequency.(a) The capillary diameter effect, (b) The capillary length effect on the mean frequency, (c) The orifice effect, (d) The restrictors type effect on the mean frequency.
Figure 12a-d show the restrictor-type effect on the bearing vibration.

Table 1 .
Results of the frequency calculations.Frequency comparison of the examined cases.