Theoretical Study of the Static and Dynamic Characteristics of a Slotted Adaptive Hydrostatic Thrust Bearing with a Regulator of the Lubricant Output Flow

: This manuscript considers the design of a slotted adaptive hydrostatic thrust bearing with a regulator of the lubricant output ﬂow. A theoretical study of its static and dynamic characteristics was carried out. The aim of the study was to test the reliability of the hypothesis concerning the possibility of obtaining a stable-to-oscillation adaptive thrust bearing of negative compliance, avoiding the need for a complex system of external combined throttling by replacing it with a simple slotted throttle. Mathematical modeling of the thrust bearing movement was carried out. The possibility of reducing the compliance to negative values is shown, providing the bearing with an adaptive function, consisting of using the structure as a bearing and as an active deformation compensator of the elastic system of a metal-cutting machine in order to improve the quality of the metalwork. Analysis of load static characteristics showed that negative compliance is provided over a wide range of loads, which can be up to 75% or more of the range of permissible bearing loads. Based on the study of dynamic characteristics, it was concluded that with a targeted selection of parameters that have a major effect on the dynamics of the structure, the considered adaptive hydrostatic thrust bearing can attain a very high quality of dynamics. It is shown that the viscous damping of the lubricating ﬁlm enclosed in the gaps and the damping of the material of the elastic ring of the regulator are important resources for ensuring the optimal dynamics of the thrust bearing.


Introduction
Sliding bearings with hydrostatic lubrication (journal bearings, thrust bearings, and guides) used in modern metal-cutting machines usually work as passive supporting structures [1,2]. At the same time, these devices can act as adaptive units, which, in addition to their load-carrying function, are capable of performing the role of automatic deformation compensators of the technological system of the machine, thereby helping to reduce errors and, consequently, reduce the time of metalworking [3,4]. The adaptive function of such bearings arises only in those designs that are able to compensate for the positive compliance of the technological system of the machine, operating in the modes of negative compliance and thereby providing the total zero compliance of the machine [5,6]. In the latter cases, the direction of the reaction vector of the lubricating film under pressure is opposite to that of the cutting force.
In hydrostatic and gas-static bearings with flow regulators, which are us stalled at the inlet of the lubricating flow, compliance reduction is obtained wit tomatic supply of an additional volume of lubricant, which contributes to an a increase in pressure in the bearing gap and, consequently, the desired complian being achieved [13,14]. Studies have shown that the inflow of additional energy utes to a significant deterioration in the dynamic quality of the structure up to th stability [15,16]. It is possible to improve the dynamics of the bearings due to th cation of the design through the use of an external combined throttling system whereby an additional hydraulic resistance is installed between the bearing gap regulator, which makes it possible to put in place an inter-throttle chamber, the of lubrication of which acts as a viscous damper for energy dissipation [17,18].
Bearings with displacement compensators are divided into two types: those compensation and those with partial compensation [19,20]. In contrast to beari active compensation of the flow rate, passive throttles are installed at the inlet of t which provides better dynamics for structures [11]. However, even such beari negative compliance have low dynamics; therefore, the ECTS is also used in the first, the displacement compensator is installed between the input passive throttl bearing gap, which are separated by an intermediate throttle of low resistance ECTS chamber [21]. A decrease in compliance in such structures is achieved b the positive compliance of the bearing gap and the negative compliance of the co tor. Partial displacement compensation is the result of a combination of a disp compensator and a lubricant output flow regulator (LOFR) [10,22]. The improv the dynamics of such structures is provided by the use of an ECTS.
In work [4], it is shown that the stability of a step hydrostatic thrust bearing LOFR can be provided without the use of an ECTS, which simplifies the design. T ematical model of the bearing takes into account the compliance of the elastic ment compensator, but what is not taken into account is that, at the same time, as a damper, which improves the dynamics of the bearing.
In this paper, we consider the design of a slotted hydrostatic thrust bearing LOFR and study its static and dynamic characteristics, taking into account the co and damping of the regulator. Figure 1 shows a diagram of a slotted adaptive hy thrust bearing with LOFR. The structure contains a shaft, 1, and a base, 2, to which a rigid ring, 4, of radius r1 and an outer radius r0 is hermetically connected by means of an elastic r an inner radius r2. The thrust bearing is powered by a lubricant injection source to the slotted throttle, 5, under pressure ps. In working conditions, a lubricatin thickness h is created between the surfaces of ring 4 and shaft 1. Between base 2 4 there is a blind annular gap of thickness ht, at the entrance to which the lubrican The structure contains a shaft, 1, and a base, 2, to which a rigid ring, 4, of an inner radius r 1 and an outer radius r 0 is hermetically connected by means of an elastic ring, 3, of an inner radius r 2 . The thrust bearing is powered by a lubricant injection source supplied to the slotted throttle, 5, under pressure p s . In working conditions, a lubricating gap of thickness h is created between the surfaces of ring 4 and shaft 1. Between base 2 and ring 4 there is a blind annular gap of thickness h t , at the entrance to which the lubricant is under pressure p t < p s , formed as a result of overcoming the hydraulic resistance to the lubricant Mathematics 2022, 10, 355 3 of 18 flow in slotted throttle 5. In the operating state, with increasing load f on the thrust bearing, the gap h decreases and its hydraulic resistance increases, as a result of which the pressure p t increases. In this case, the force reactions w t and w e of the lubricant film (which is under pressure in the gaps h t and h) both increase. At w t > w e , under the action of a force equal to the difference between these reactions, the axial deformation e of elastic ring 3 increases, limiting the gap h, again increasing its hydraulic resistance and pressure p t , which brings about a decrease in the compliance of the thrust bearing where is the thickness of the lubricating gap enclosed between the working surfaces of base 2 and shaft 1, e is the deformation of ring 3. The force w t − w e is balanced by the force of resistance to deformation of ring 3. By varying the size of the elastic ring and its axial compliance k e , it is possible to reduce the compliance of the thrust bearing k to negative values. The paper considers a mathematical model of the unsteady state of the thrust bearing; static and dynamic characteristics of the bearing are calculated and studied.

Mathematical Modeling
When studying the characteristics of the thrust bearing, we take a plane-parallel movement model of the working surfaces of the moving elements in which the thickness of the lubricant layer in the gaps does not depend on the radius. In addition, we assume that the loaded thrust bearing has reached a steady-state temperature regime at which the dynamic viscosity of the lubricant is constant.
It is convenient to study the characteristics of the thrust bearing in a dimensionless form. This allows us to significantly reduce the number of variable parameters that affect the characteristics of the thrust bearing. Following [4], we take the scales of values: p s for pressures, r 0 for radii, t 0 for the current time, h 0 for gaps and deformation e, πh 3 0 p s 6µ for volumetric lubricant flow rate, 2πr 2 0 p s for axial forces, where h 0 corresponds to the gap h in the bearing that carries the design load f 0 , and µ is the lubricant viscosity. Dimensionless quantities are signified by Greek and capital Latin letters.
The reactions of pressure forces in mentioned areas are determined by the general formula [4] where P(R,τ) is a function of the pressure in the clearances of the thrust bearing. The formula for the lubricant flow rate at the inlet and outlet of the corresponding gaps has the form [4] The pressure function satisfies the unsteady Reynolds equation for an incompressible fluid [22] ∂ ∂R where H i (τ) is the function of the gap, R is the current radius, σ = 12µr 2 0 h 2 0 p s t 0 is the so-called "compression number" of the film of the liquid in the gap, and τ is the dimensionless time.
On the ring R ∈ [R 1 , 1] in the gap H region, the boundary value problem for Equation (5) satisfies the boundary conditions For the blind gap region, the boundary conditions for the boundary value problem with Equation (5) are The solution of the boundary value problem (5), (6) is the function On the ring in the region of the gap H, for the boundary value problem (5), (7) has the solution On the ring R ∈ [R 1 , R 2 ], in the region of the blind gap H t , the solution to the boundary value problem (5), (8) is the function Substituting (8), (9), and (10) into (3), after integration, we find expressions for the force reactions where The bearing capacity of the thrust bearing is determined by the expression The formula for the force equal to the difference between the reactions causing the deformation of ring 3 has the form The flow through the slotted throttle is determined by the formula where A s is the criterion for the similarity of the throttle [12].
The flow rate at the outlet from the gap H s in accordance with formula (4) has the form where B 3 = σR 2 1 2 . At the entrance to the gap H, the flow rate is where At the entrance to the gap H t , the flow rate is determined by the formula where B 5 = σ 2 R 2 2 − R 2 1 . Formulae for hydrostatic forces and flow rates, derived from the analytical solutions to boundary value problems for the Reynolds equation, are verified using an alternative numerical finite difference method [23]. In a dynamic system, the force of inertia of the mass of shaft 1 acts as where M = mh 0 2πr 2 0 p s t 2 0 , m is the mass of the shaft.
The bearing capacity of the thrust bearing is Hydrostatic force acts on the working surfaces of ring 3 The mathematical model of the thrust bearing dynamics includes two equations for the force balance of shaft 1 and ring 3, one equation for the balance of lubricant flow rate, and formulae for determining the gap H t and the total gap H s where ε is the dimensionless axial deformation of ring 3, and H t0 is the thickness of the blind gap in the absence of load on ring 3. Formula (23) is obtained on the basis of Hooke's law [24] for small axial deformations of elastic ring 3, where K e is the coefficient of its axial compliance (elasticity), and D e is the damping coefficient of the elastic ring material. Substituting (21)- (23) and (24) into (25), (26), and (28), and taking into account (27) and the fact that we represent the mathematical model of the thrust bearing dynamics in the form of the following combined system of nonlinear differential and algebraic equations The closed system of three nonlinear differential Equations (29)-(31) and one algebraic Equation (32) can be considered as a continuous automatic control system [25] containing an input-disturbing function of the external load F and four output-controlled functions H, H s , ε, and P t .
First, we conduct a study of the static characteristics of the thrust bearing in order to detect modes of negative compliance that provide the thrust bearing with an adaptive function. Then, we conduct a study of the stability criteria and the stability margin of the dynamic system, determine the optimal dynamic modes under which the structure has the maximum speed and stability margin in the negative compliance modes, and establish the particular parameters that have a decisive effect on the quality of the thrust bearing dynamics.

Static Characteristics of the Thrust Bearing
The mathematical model of the stationary state of the thrust bearing is a special case of the model of dynamics in the absence of oscillations, that is, with zero derivatives of the output functions H, H s , ε, and P t . The model can be represented as the following nonlinear system of four algebraic equations The study of the static characteristics is preceded by the adjustment of the thrust bearing parameters to the mode of the so-called "design point". This mode corresponds to a point on the load curve that in turn corresponds to a unit dimensionless gap H = 1, at which the dimensional gap is equal to the design gap h = h 0 .

Calculation of the Parameters of the "Design Point" Mode
It is convenient to adjust the thrust bearing to this mode by setting the pressure at the outlet of slot 5 using the normalized pressure setting factor χ ∈ [0, 1]. At the "design point", the values of pressure P t and this coefficient coincide so that P t = χ. We also set the radii R 1 and R 2 and the coefficient of elasticity K e . Using Formulae (31)-(34), we calculate the deformation ε, the criterion for the similarity of the gap A d , the gap H s , the bearing capacity W, and the external force F It was earlier mentioned that to ensure the adaptive function, it is necessary that the condition W e > 0 be satisfied. It follows from (21) that this takes place at A 5 > 0. This condition allows us to find the smallest radius of elastic ring 3, at which can be ensured a decrease in the compliance of the thrust bearing Let us find the value of the elasticity coefficient K e = K e0 , at which the thrust bearing at the "design point" will have zero compliance K = 0.
The compliance of the structure is determined by the formula Substituting (32) and (34) into (35), and performing the differentiation of this expression, taking into account (33) we obtain the formula for the compliance of the thrust bearing at the "design point" It follows that the value of the elasticity coefficient K e0 = K e , at which zero compliance of the thrust bearing K 0 = 0 is achieved at the "design point", can be calculated using the formula Formula (37) is convenient to use for calculating the static characteristics of bearing capacity, compliance, and lubricant flow rate, setting the elasticity coefficient K e in proportion to the value of K e0 . With K e = K e0 , the thrust bearing at the "design point" will have zero compliance K = 0, with 0 ≤ K e < K e0 -positive and K e > K e0 -negative.

Static Characteristics and Their Discussion
It is convenient to calculate the static characteristics of the thrust bearing in a parametric form, taking pressure as the parameter P t ∈ [0, 1]. Using Formula (31), one can determine the external load F and the bearing capacity of the thrust bearing W, which are equal in the static mode. Equation (32) allows us to find the deformation of elastic ring 3 with the formula ε = K e A 5 P t .
Transforming (33), we obtain the formula for calculating the gap Finally, using (34), we obtain a formula for calculating the total clearance Taking into account (31), (32), and (39)-(41), we differentiate (30) and (33) with respect to the pressure P t and substitute its results in (36). In this case, we obtain a formula for calculating the static compliance at any permissible load F on the thrust bearing Analysis of Formula (42) shows that, at extreme values of the load, Therefore, there are certain load values at which the compliance reaches a minimum. To find these values, let us determine the derivative of the compliance with respect to the load, which can be obtained by differentiating (42) with respect to the pressure P t It can be easily observed that this function has a single zero P t = 2/3, which is always constant and does not depend on the parameters of the thrust bearing, including the coefficient of elasticity K e . Thus, the optimal mode of least compliance is achieved at χ = 2/3, at which the minima of the dependences K(F) and K(χ) coincide. Figure 2 shows the graph of the load characteristics for various values of the coefficient of elasticity K e . The adjustment factor is χ = 2/3, radii R 1 = 0.5, R 2 = 0.9 > R 2,min = 0.736, at which, according to (35), the adaptive function of the thrust bearing can be provided. The coefficient of elasticity of ring 3, according to (38), is K e0 = 11.2. With this value, the thrust bearing at the "design point" has zero compliance. For each curve in Figure 2, the coefficient of elasticity is set by the formula K e = i K e0 . It can be easily observed that this function has a single zero Pt = 2/3, which is always constant and does not depend on the parameters of the thrust bearing, including the coefficient of elasticity Ke. Thus, the optimal mode of least compliance is achieved at χ = 2/3, at which the minima of the dependences K(F) and K(χ) coincide. Figure 2 shows the graph of the load characteristics for various values of the coefficient of elasticity Ke. The adjustment factor is χ = 2/3, radii R1 = 0.5, R2 = 0.9 > R2,min = 0.736, at which, according to (35), the adaptive function of the thrust bearing can be provided. The coefficient of elasticity of ring 3, according to (38), is Ke0 = 11.2. With this value, the thrust bearing at the "design point" has zero compliance. For each curve in Figure 2, the coefficient of elasticity is set by the formula Ke = i Ke0. The thrust bearing has a maximum permissible load Fmax = 0.271. The "design point" corresponds to the load F0 = χFmax = 0.181. The curve for i = 0 corresponds to the load characteristic of a conventional thrust bearing, for which the coefficient of elasticity is Ke = 0. The graph shows that with an increase in i and, consequently, an increase in the coefficient Ke, the nature of the curves changes. In this case, in the region of small and moderate values of the load F, the compliance K of the thrust bearing decreases, which can be determined by the change in the slope of the curves.
At Ke = Ke0 (i = 1), at the "design point", the thrust bearing acquires zero compliance K = 0. A further increase in the coefficient for Ke > Ke0 (i > 1) contributes to the establishment of a negative compliance mode. This mode is characterized by the presence of areas of simultaneous increases in load and clearance. It is in these areas that the thrust bearing acquires an adaptive function as it is able to perform not only the load bearing but also adaptive functions in order to actively compensate for the passive deformation of the elastic system of the machine due to the opposite direction of the force deformation of the lubricating gap of the thrust bearing.
The thrust bearing has a maximum permissible load F max = 0.271. The "design point" corresponds to the load F 0 = χF max = 0.181. The curve for i = 0 corresponds to the load characteristic of a conventional thrust bearing, for which the coefficient of elasticity is K e = 0. The graph shows that with an increase in i and, consequently, an increase in the coefficient K e , the nature of the curves changes. In this case, in the region of small and moderate values of the load F, the compliance K of the thrust bearing decreases, which can be determined by the change in the slope of the curves.
At K e = K e0 (i = 1), at the "design point", the thrust bearing acquires zero compliance K = 0. A further increase in the coefficient for K e > K e0 (i > 1) contributes to the establishment of a negative compliance mode. This mode is characterized by the presence of areas of simultaneous increases in load and clearance. It is in these areas that the thrust bearing acquires an adaptive function as it is able to perform not only the load bearing but also adaptive functions in order to actively compensate for the passive deformation of the elastic system of the machine due to the opposite direction of the force deformation of the lubricating gap of the thrust bearing.
The effect of negative compliance is illustrated by the curves of dependence of the compliance K on the external force F, which are shown in Figure 3. With an increase in the coefficient of elasticity K e , the compliance K decreases in the range of working loads. At K e = K e0 , the thrust bearing reaches zero compliance for χ = 2/3 only at the "design point". At K e > K e0 , the load range-where the thrust bearing has zero and negative compliance-becomes wider. In the range of loads F corresponding to negative compliance K, as mentioned above the thrust bearing is able to carry the load and compensate for the elastic deformations o the machine, that is, it is able to perform an adaptive function. At Ke = 2Ke0, when at the "design point" the negative compliance in absolute value coincides with the compliance of a conventional thrust bearing (Ke = 0), the width of the adaptive load range in negative compliance modes is approximately 75% of the allowable load range width. With Ke > 2Ke0 the adaptive load range will be wider still.
The graph in Figure 4 shows the dependences ε(F) of the deformation of elastic ring 4 on the external load on the bearing. These dependences are linear. To achieve negative compliance at the "design point", the deformation of elastic ring 3 must meet the condi tion ε = Ke/Ke0 > 1. In the range of loads F corresponding to negative compliance K, as mentioned above, the thrust bearing is able to carry the load and compensate for the elastic deformations of the machine, that is, it is able to perform an adaptive function. At K e = 2K e0 , when at the "design point" the negative compliance in absolute value coincides with the compliance of a conventional thrust bearing (K e = 0), the width of the adaptive load range in negative compliance modes is approximately 75% of the allowable load range width. With K e > 2K e0 , the adaptive load range will be wider still.
The graph in Figure 4 shows the dependences ε(F) of the deformation of elastic ring 4 on the external load on the bearing. These dependences are linear. To achieve negative compliance at the "design point", the deformation of elastic ring 3 must meet the condition ε = K e /K e0 > 1.
It was mentioned above that as compliance decreases, the dynamic quality of the thrust bearing usually deteriorates, often to the point of loss of stability. Therefore, an important condition for the operability of an adaptive thrust bearing is its dynamic performance, which must meet certain criteria for response speed and stability margin. The speed of response is important so that, with a minimum duration of the transient response, the thrust bearing reaches the calculated negative compliance mode, and the stability margin is important from the point of view of guaranteeing the accident-free operation of the structure.
Thus, to ensure the operability of the adaptive thrust bearing, it is necessary to optimize its parameters based on the study of its dynamic characteristics in order to deliver the structure with sufficient dynamic quality. compliance modes is approximately 75% of the allowable load range width. With Ke > 2Ke0 the adaptive load range will be wider still.
The graph in Figure 4 shows the dependences ε(F) of the deformation of elastic ring 4 on the external load on the bearing. These dependences are linear. To achieve negative compliance at the "design point", the deformation of elastic ring 3 must meet the condition ε = Ke/Ke0 > 1. It was mentioned above that as compliance decreases, the dynamic quality of the thrust bearing usually deteriorates, often to the point of loss of stability. Therefore, an important condition for the operability of an adaptive thrust bearing is its dynamic performance, which must meet certain criteria for response speed and stability margin. The

Dynamic Characteristics of the Thrust Bearing
Let us consider the dynamics of the thrust bearing, which, under the influence of small disturbances of the external load, bring about small fluctuations in the output values relative to their values in the steady (stationary) state. Let us take a unit dimensionless mass M = 1 of shaft 1. This allows us to determine the scale of the current time.

Linear Model of Thrust Bearing Dynamics
We represent dynamic quantities in the form where the first term in the expressions represents the static constituent and the second the small dynamic constituent. After substituting (44) into the nonlinear system of Equations (27)- (30), with its linearization and application to the resulting linear system of equations of the integral Laplace transform [25], we obtain a linear model of the dynamics of the studied thrust bearing in Laplace transforms. Eliminating ∆H, we obtain a system of three equations. The matrix form of this system is where ∆F(s) is the Laplace transform of a small deviation of the disturbing external force, ∆H s (s), ∆ε(s), ∆P t (s) are the transforms of small deviations of the corresponding output quantities, and s is the Laplace transform variable, Applying Cramer's rule [26] and the algorithm of symbolic expansion of the determinants of polynomial matrices to the solution of the system of linear Equation (45), we determined the transfer function (TF) of the dynamic compliance of the thrust bearing, which is a rational function of the variable s where The denominator of the TF (46) is the characteristic polynomial of the system, by which it is possible to judge the stability and quality of transients in a linear dynamic system that describes the dynamics of an adaptive hydrostatic thrust bearing for small deviations of non-stationary quantities from their stationary state.

Some Considerations Regarding the Methodology for Applying the Linear Dynamics Model to the Study of Sliding Bearings with Liquid and Gas Lubrication
When studying the quality of the dynamics of linear oscillatory systems, the model of the simplest damped harmonic oscillator (DHO) with a second-order characteristic equation is often used j + cs + ms 2 = 0, where j, c, and m refer to the stiffness, damping, and mass of the oscillator [27], respectively. According to the Hurwitz conditions, such a system is stable if all the coefficients of Equation (47) are nonzero and have the same sign [28]. As the mass and damping of the generator are always positive, the stiffness should also be positive in stable systems. However, in slide bearings with regulators, which include the studied design, the stiffness, which is the opposite of the compliance, can be infinite or negative. On this basis, it could be assumed that the adaptive bearings of negative compliance are always unstable. However, this assumption is wrong.
In accordance with the basic provisions of the theory of automatic control, the characteristic polynomial of the control object is the denominator of the TF, which for DHO has the form [25] Φ(s) = b 0 a 0 + a 1 s + a 2 s 2 .
The characteristic polynomial of the TF (48) can be reduced to the form of the polynomial (47) by multiplying the first by the constant b 0 . Studies of adaptive slip bearings show that the sign of static compliance is determined by the coefficient b 0 [18,19]. In this case, the coefficients of the denominator polynomial for stable systems always remain positive, which indicates the stability of such dynamical systems. This also applies to the bearing under study. If we turn to TF (46), then only the coefficient b 0 depends on the compliance K e of elastic ring 3, and the coefficient a 0 does not. Thus, an oscillator, which is an adaptive hydrostatic thrust bearing, can, like similar structures, have negative static compliance and still be stable.
Note that PFs of the form (48) and (46), which are the Laplace image of dynamic compliance, are rational functions; therefore, from the point of view of the theory of automatic control, the preferred term is "compliance" rather than "stiffness", the latter term often being used for the analysis of the characteristics of sliding bearings, including adaptive bearings. A particular case of a TF at s = 0 is the static compliance K(0) = b 0 /a 0 . Therefore, in this work and other works of the authors on the study of bearings with lubricant flow controllers and displacement compensators that provide the bearings with an adaptive function, the term "compliance" is used [4,5,10,11].

Dynamic Characteristics of the Thrust Bearing and Their Discussion
To assess the quality of the thrust bearing dynamics, we used the root criteria-the degree of stability and damping of oscillations over the period [25]. The degree of stability η characterizes the speed of the system, that is, the speed of damping of its free oscillations. By this criterion, one can determine the degree of stability of the system. For η < 0 the system is unstable, for η > 0 it is stable, and for η = 0 the system is on the boundary between stable and unstable. The criterion ξ of the damping of oscillations over a period can be applied to the assessment of the stability margin of the system. The smaller the value of ξ, the greater the oscillation of the transient response and the lower the stability margin of the dynamic system. It is considered that the system is well damped if ξ ≥ 90% [25].
The calculation of the dynamic criteria was carried out for the same parameter values that were used to calculate the static characteristics of the thrust bearing. Three new parameters were added to the dynamics model: - Unit mass of the shaft M = 1, taking into account the influence of the force of its inertia; -"Compression number" σ, which is a measure of the viscous damping of the lubricant film under pressure; -Damping coefficient D e of the elastic ring material.
As the first parameter is constant, the other two parameters, as well as the parameter H t0 , are varied.
Calculations have shown that an increase in the parameter H t0 negatively affects the dynamics of the thrust bearing. The best value is H t0 = 0.5.
Let us consider the influence of the parameter σ on the dynamic characteristics of the structure. Figure 5 shows the curves of dependence of the degree of stability η on the parameter σ for different values of the elasticity coefficient K e without taking into account the damping coefficient of elastic ring 3 (D e = 0). The graph is constructed for the values of the elasticity coefficient K e , which provide the thrust bearing negative static compliance in the "design point" mode.
value of ξ, the greater the oscillation of the transient response and the lower the stability margin of the dynamic system. It is considered that the system is well damped if ξ ≥ 90% [25].
The calculation of the dynamic criteria was carried out for the same parameter values that were used to calculate the static characteristics of the thrust bearing. Three new pa rameters were added to the dynamics model:

−
Unit mass of the shaft M = 1, taking into account the influence of the force of its iner tia; − "Compression number" σ, which is a measure of the viscous damping of the lubri cant film under pressure; − Damping coefficient De of the elastic ring material.
As the first parameter is constant, the other two parameters, as well as the parameter Ht0, are varied.
Calculations have shown that an increase in the parameter Ht0 negatively affects the dynamics of the thrust bearing. The best value is Ht0 = 0.5.
Let us consider the influence of the parameter σ on the dynamic characteristics of the structure. Figure 5 shows the curves of dependence of the degree of stability η on the parameter σ for different values of the elasticity coefficient Ke without taking into accoun the damping coefficient of elastic ring 3 (De = 0). The graph is constructed for the values of the elasticity coefficient Ke, which provide the thrust bearing negative static compliance in the "design point" mode.
Comparing the graphs of Figures 2 and 3 with Figure 5, we find that the stability (η > 0) of the thrust bearing of negative compliance is provided only for σ > σmin, where σmi is the value of the parameter σ at which the bearing reaches the stability boundary (η = 0) Therefore, for a bearing of negative compliance Ke = 2Ke0, the stability limit corresponds to σmin = 42. At lower values of σ, the thrust bearing is unstable (η < 0). It can be seen that the larger the value of Ke, that is, the lower the compliance K, the larger the σmin corresponding to the stability boundary. This result can be interpreted in such a way that the lower the compliance of the thrust bearing, the smaller the calculated dimensional clearance of the bearing lubricant gap must be to ensure the stability of the thrust bearing at the "design point".  Comparing the graphs of Figures 2 and 3 with Figure 5, we find that the stability (η > 0) of the thrust bearing of negative compliance is provided only for σ > σ min , where σ min is the value of the parameter σ at which the bearing reaches the stability boundary (η = 0). Therefore, for a bearing of negative compliance K e = 2K e0 , the stability limit corresponds to σ min = 42. At lower values of σ, the thrust bearing is unstable (η < 0). It can be seen that the larger the value of K e , that is, the lower the compliance K, the larger the σ min corresponding to the stability boundary. This result can be interpreted in such a way that the lower the compliance of the thrust bearing, the smaller the calculated dimensional clearance of the bearing lubricant gap must be to ensure the stability of the thrust bearing at the "design point".
It is also noteworthy that the dependences η(σ) have a pronounced extreme character, which indicates the existence of an optimal value of the parameter σ = σ opt . At this value, the dependence η(σ) reaches a maximum, providing the thrust bearing with the maximum response rate, at which the fastest decay of transient processes is observed. For K e ≥ K e0 , the optimal performance takes place at approximately σ opt = 60.
A comparison (presented in Figure 6) of the graph of the dependence of the damping criterion ξ of oscillations for the period on the parameter σ with the graph in Figure 5 shows that, in the stability region at σ < σ opt , the transient response has a pronounced oscillatory character (ξ < 80%). When approaching the point σ = σ opt , the criterion ξ grows rapidly, which indicates a significant weakening of the oscillation of the transient characteristics and an increase in the stability margin of the thrust bearing. At the same time, for thrust bearings of low compliance (K e > 2K e0 ), the transient characteristics even at σ > σ opt remain somewhat oscillatory (ξ < 90%).
Mathematics 2022, 9, x FOR PEER REVIEW 13 of 18 It is also noteworthy that the dependences η(σ) have a pronounced extreme character, which indicates the existence of an optimal value of the parameter σ = σopt. At this value, the dependence η(σ) reaches a maximum, providing the thrust bearing with the maximum response rate, at which the fastest decay of transient processes is observed. For Ke ≥ Ke0, the optimal performance takes place at approximately σopt = 60.
A comparison (presented in Figure 6) of the graph of the dependence of the damping criterion ξ of oscillations for the period on the parameter σ with the graph in Figure 5 shows that, in the stability region at σ < σopt, the transient response has a pronounced oscillatory character (ξ < 80%). When approaching the point σ = σopt, the criterion ξ grows rapidly, which indicates a significant weakening of the oscillation of the transient characteristics and an increase in the stability margin of the thrust bearing. At the same time, for thrust bearings of low compliance (Ke > 2Ke0), the transient characteristics even at σ > σopt remain somewhat oscillatory (ξ < 90%). The graph in Figure 7 shows the dependences of the degree of stability η on the viscous damping coefficient De for different values of the elasticity coefficient Ke for a fixed value of σ = 60. The graph in Figure 7 shows the dependences of the degree of stability η on the viscous damping coefficient D e for different values of the elasticity coefficient K e for a fixed value of σ = 60.
The graph shows that the damping of the elastic ring material 3 has a significant impact on the dynamics of the thrust bearing. This influence is especially noticeable with the thrust bearing operating in the modes of negative compliance when the performance at optimal D e values can be increased several times. Therefore, for K e = 2K e0 without taking into account the coefficient D e , the optimal criterion is η = 0.03. As can be seen from Figure 7, with the optimal D e = 11, the criterion η = 0.15. That is, the choice of the elastic compensator material, which is optimal from the point of view of damping, enables the speed of the thrust bearing to be increased by a factor of 5.
As seen from Figure 8, the parameter D e has a particularly significant effect on the criterion ξ of damping of oscillations over a period. For D e > 3, the criterion is virtually equal to ξ = 100%. This demonstrates that the corresponding transient processes are aperiodic in nature, indicating a minimum oscillation or even the complete absence of one. This result shows that the negative compliance thrust bearing, designed to take into account the optimal recommendations, has an absolute margin of stability to fluctuations. The graph in Figure 7 shows the dependences of the degree of stability η on the vis cous damping coefficient De for different values of the elasticity coefficient Ke for a fixed value of σ = 60.  The graph shows that the damping of the elastic ring material 3 has a significant impact on the dynamics of the thrust bearing. This influence is especially noticeable with the thrust bearing operating in the modes of negative compliance when the performance at optimal De values can be increased several times. Therefore, for Ke =2Ke0 without taking into account the coefficient De, the optimal criterion is η = 0.03. As can be seen from Figure  7, with the optimal De = 11, the criterion η = 0.15. That is, the choice of the elastic compensator material, which is optimal from the point of view of damping, enables the speed of the thrust bearing to be increased by a factor of 5.
As seen from Figure 8, the parameter De has a particularly significant effect on the criterion ξ of damping of oscillations over a period. For De > 3, the criterion is virtually equal to ξ = 100%. This demonstrates that the corresponding transient processes are aperiodic in nature, indicating a minimum oscillation or even the complete absence of one This result shows that the negative compliance thrust bearing, designed to take into account the optimal recommendations, has an absolute margin of stability to fluctuations. It is interesting to consider the quality of the thrust bearing dynamics at various values of the load across the entire range of its variation. The graphs of Figures 9 and 10 show the curves of the degree of stability η and damping ξ of oscillations for the period on the load F for various values of the coefficient Ke at σ = 60 and De = 11.
The graph in Figure 9 shows that the thrust bearing remains stable for all values of the Ke coefficient over the entire range of load variation. In this case, the maximum speed at which the function η(F) reaches its maximum value falls on loads close to the calculated value (F = 0.181). On this basis, a significant conclusion can be reached-that the adaptive bearing can be operable not only in the "design point" mode but also, importantly, across the entire range of permissible loads.
The stability margin of the thrust bearing across the range of permissible loads is represented by the graph in Figure 10. As can be seen from the graph, the thrust bearing has an ideal stability margin over almost the entire range of permissible loads (ξ = 100%) It is interesting to consider the quality of the thrust bearing dynamics at various values of the load across the entire range of its variation. The graphs of Figures 9 and 10 show the curves of the degree of stability η and damping ξ of oscillations for the period on the load F for various values of the coefficient K e at σ = 60 and D e = 11.  The study of the static and dynamic characteristics of the slotted adaptive hydrostatic thrust bearing shows that an almost ideal quality of its dynamics can be achieved withou the use of an external combined throttling system that would over-complicate the design Viscous damping of the pressurized lubricating film in the gaps of the structure and damping of the material of the elastic ring of the regulator are important resources fo ensuring optimal dynamics of the thrust bearing.

Checking the Correctness of the Obtained Data
To confirm the correctness of the obtained data, they were compared with the theo retical and experimental data of other authors. The comparison was carried out for hydro static thrust bearings, in which ring 4 was motionless, which corresponds to the coefficien of elasticity of the compensator Ke = 0 (absolutely rigid ring 3). The data were compared with the theoretical data from [29] and the experimental results from [30]. Table 1 shows comparative theoretical data on the static load capacity of the thrus bearing for radii r1 = 30 mm, r0 = 40 mm for various values of external pressure ps at pt =  The study of the static and dynamic characteristics of the slotted adaptive hydrostatic thrust bearing shows that an almost ideal quality of its dynamics can be achieved without the use of an external combined throttling system that would over-complicate the design Viscous damping of the pressurized lubricating film in the gaps of the structure and damping of the material of the elastic ring of the regulator are important resources for ensuring optimal dynamics of the thrust bearing.

Checking the Correctness of the Obtained Data
To confirm the correctness of the obtained data, they were compared with the theoretical and experimental data of other authors. The comparison was carried out for hydrostatic thrust bearings, in which ring 4 was motionless, which corresponds to the coefficient of elasticity of the compensator Ke = 0 (absolutely rigid ring 3). The data were compared with the theoretical data from [29] and the experimental results from [30]. Table 1 shows comparative theoretical data on the static load capacity of the thrust bearing for radii r1 = 30 mm, r0 = 40 mm for various values of external pressure ps at pt = 0.5ps. The graph in Figure 9 shows that the thrust bearing remains stable for all values of the K e coefficient over the entire range of load variation. In this case, the maximum speed at which the function η(F) reaches its maximum value falls on loads close to the calculated value (F = 0.181). On this basis, a significant conclusion can be reached-that the adaptive bearing can be operable not only in the "design point" mode but also, importantly, across the entire range of permissible loads.
The stability margin of the thrust bearing across the range of permissible loads is represented by the graph in Figure 10. As can be seen from the graph, the thrust bearing has an ideal stability margin over almost the entire range of permissible loads (ξ = 100%). A decrease in the values of the criterion ξ is observed only in the vicinity of the extreme points at which the thrust bearing loses its load-carrying function.
The study of the static and dynamic characteristics of the slotted adaptive hydrostatic thrust bearing shows that an almost ideal quality of its dynamics can be achieved without the use of an external combined throttling system that would over-complicate the design. Viscous damping of the pressurized lubricating film in the gaps of the structure and damping of the material of the elastic ring of the regulator are important resources for ensuring optimal dynamics of the thrust bearing.

Checking the Correctness of the Obtained Data
To confirm the correctness of the obtained data, they were compared with the theoretical and experimental data of other authors. The comparison was carried out for hydrostatic thrust bearings, in which ring 4 was motionless, which corresponds to the coefficient of elasticity of the compensator K e = 0 (absolutely rigid ring 3). The data were compared with the theoretical data from [29] and the experimental results from [30]. Table 1 shows comparative theoretical data on the static load capacity of the thrust bearing for radii r 1 = 30 mm, r 0 = 40 mm for various values of external pressure p s at p t = 0.5p s . The presented theoretical results are in good agreement. A slight error is explained by the fact that in [29] a simplified method for calculating the bearing capacity was used, assuming that the pressure of the lubricating film on the surface of the outer ring is equal to half the pressure at the outlet of their supply gap. Table 2 shows comparative theoretical and experimental data on the static bearing capacity of thrust bearings for radii r 1 = 24 mm and r 0 = 32 mm for different values of the gap h at an external pressure p s = 0.8 MPa. A comparison of theoretical and experimental data provides a more significant difference. However, such a result can be considered quite satisfactory as the indicated data on the difference between theory and experiment can be explained by many reasons known to researchers, therefore the experiment as a whole confirms theoretical calculations. This allows us to conclude that for the considered particular case of the mathematical model of the thrust bearing, confirmation of its adequacy has been obtained.
As for the model of the movement of ring 4 on elastic suspension 3, for such a bearing the conduction of an experimental study is being planned, in the process of preparation and implementation of which several important practical problems will have to be solved, which are complex and beyond the scope of this paper. In this regard, the present theoretical study is fundamental to a certain extent, and the data obtained from the verification of a particular case provide reason to trust the calculated data and general conclusions about the performance of the considered design in the context of assessing its quantitative performance characteristics.

Conclusions
This paper considered the design of a slotted adaptive hydrostatic thrust bearing with a regulator of the lubricant output flow and conducted a theoretical study of its static and dynamic characteristics. The aim of the study was to test the reliability of the hypothesis that it is possible to obtain a vibration-resistant adaptive thrust bearing of negative compliance without having to use a complex system of external combined throttling by replacing it with a simple slotted throttle. Mathematical modeling of the movement of the thrust bearing and its structural elements was carried out. The paper presented the possibility of a significant decrease in the compliance of the structure, including negative compliance in order to provide an adaptive function, which consists of using the bearing not only as a support but also as an active compensator for the deformation of the elastic system of a metal-cutting machine in order to improve the quality of the metalwork. It was found that negative compliance is provided across a wide range of loads, which can be up to 75% or more of the range of permissible bearing loads. The study showed that with a purposeful selection of parameters that have the main influence on the dynamics of the structure, the considered adaptive hydrostatic thrust bearing can attain an almost ideal quality of its dynamics. Viscous damping of the lubricating film enclosed in the gaps and damping of the elastic ring material of the regulator are important resources for ensuring optimal thrust bearing dynamics.