Bearing dynamic parameter determination method and device, storage medium and terminal

Abstract

The invention discloses a bearing power parameter determination method and device, a storage medium and a terminal, relates to the technical field of compressors, and mainly aims to solve the problem that the bearing operation condition analysis effectiveness is affected due to poor bearing power parameter identification accuracy of an existing gear assembly type compressor. Comprising the following steps: determining a static characteristic initial value based on a bearing load, a bearing geometric parameter, a lubrication parameter and a rotating speed, and constructing a bearing three-dimensional model based on the static characteristic initial value and the bearing geometric parameter; based on the bearing three-dimensional model, determining the bearing tile torque and journal resultant force deviation of the bearing; and when the deviations of the bearing pad torque and the shaft neck resultant force are both smaller than a first preset deviation threshold value and the static characteristic parameter deviation is smaller than a second preset deviation threshold value, bearing dynamic characteristic analysis processing is carried out based on a small displacement vortex motion analysis method, and bearing dynamic parameters are obtained.

Description

Technical Field

[0001] This application relates to the field of compressor technology, and in particular to a method and apparatus for determining bearing dynamic parameters, a storage medium, and a terminal. Background Technology

[0002] Multi-shaft assembled compressors have wide applications in fluid machinery. Compared to single-shaft compressors, multi-shaft assembled compressors face more complex load conditions. The loads faced by the bearings include not only the mass of the shaft itself, but also the cantilever mass and constantly changing meshing forces. With the continuous increase in overall power, the bearing operating load range becomes wider, and the load magnitude can vary by more than ten times under different operating conditions. Therefore, the analysis and prediction of bearing dynamic characteristics have a profound impact on the stable operation of the rotor of the assembled compressor.

[0003] Currently, existing bearing dynamic parameters are usually predicted using engineering analysis software through numerical calculation methods. However, since the bearing dynamic parameters for different load sizes and directions can produce large errors through numerical calculation alone, it affects the judgment of bearing operation and results in a large difference between the actual performance of the unit in practice. Therefore, there is an urgent need for a method to determine the bearing dynamic parameters to solve the above problems. Summary of the Invention

[0004] In view of this, this application provides a method, apparatus, storage medium, and terminal for determining bearing dynamic parameters, with the main purpose of solving the problem that the accuracy of existing bearing dynamic parameter identification is poor, which in turn affects the effectiveness of bearing operation analysis.

[0005] According to one aspect of this application, a method for determining bearing dynamic parameters is provided, comprising: The initial values ​​of static characteristics are determined based on the bearing load, bearing geometric parameters, lubrication parameters, and rotational speed, and a three-dimensional model of the bearing is constructed based on the initial values ​​of static characteristics and the bearing geometric parameters. The bearing pad torque and journal resultant force deviation of the bearing are determined based on the three-dimensional model of the bearing. When the deviations between the bearing pad torque and the journal resultant force are both less than the first preset deviation threshold, and the deviation of the static characteristic parameters is less than the second preset deviation threshold, the bearing dynamic characteristic analysis is performed based on the small displacement eddy analysis method to obtain the bearing dynamic parameters.

[0006] Furthermore, the method also includes: When the deviation between the bearing pad torque and the journal resultant force is greater than or equal to a first preset deviation threshold, the bearing pad position, journal position, and static characteristic parameter deviation are iteratively determined based on the adjusted shaft center position and the bearing pad swing angle. When the deviations of the bearing pad torque and the journal resultant force after iteration are both less than the first preset deviation threshold, and the deviations of the static characteristic parameters after iteration are greater than or equal to the second preset deviation threshold, the initial static characteristic value is re-determined, so as to re-execute the step of constructing a three-dimensional bearing model based on the initial static characteristic value and the bearing geometric parameters.

[0007] Furthermore, when the deviations between the bearing pad torque and the journal resultant force are both less than a first preset deviation threshold, and the deviation of the static characteristic parameters is less than a second preset deviation threshold, the bearing dynamic characteristic analysis is performed based on the small displacement eddy analysis method to obtain the bearing dynamic parameters, including: When the deviation between the bearing pad torque and the journal resultant force is less than a first preset deviation threshold, static characteristic parameters are obtained, and the difference between the static characteristic parameters and the initial static characteristic value is calculated to obtain the static characteristic parameter deviation. When the deviation of the static characteristic parameter is less than the second preset deviation threshold, the dynamic characteristic calculation of the bearing is carried out. The small displacement eddy analysis method is adopted, the eddy displacement signal is applied, and the oil film excitation force is collected in real time. Based on data processing of the oil film excitation force and eddy displacement, the bearing dynamic parameters are obtained.

[0008] Furthermore, determining the bearing pad torque and journal resultant force deviation based on the bearing's three-dimensional model includes: The target bearing pad position and the target journal position are determined based on the three-dimensional model of the bearing. Based on the fluid dynamics components, the bearing pressure field characteristics of the target bearing pad position and the target journal position are calculated to obtain the bearing pad torque and the oil film resultant force. The journal resultant force deviation is obtained by calculating based on the oil film resultant force and bearing load.

[0009] Furthermore, the construction of the bearing three-dimensional model based on the initial static characteristic value and the bearing geometric parameters includes: The 3D modeling component is activated, and based on the initial static characteristic values ​​and the bearing geometric parameters, the 3D modeling component is driven to generate a 3D model of the bearing; The bearing three-dimensional model is a three-dimensional oil film flow field model corresponding to the connected domain composed of the oil inlet, oil bladder area and bearing oil film area. The oil bladder area is formed between each bearing pad, and the bearing oil film area is formed between the bearing pad and the journal.

[0010] Furthermore, the determination of the initial static characteristic values ​​based on bearing load, bearing geometric parameters, lubrication parameters, and rotational speed includes: Construct the Reynolds equation and balance equation for the bearing load, bearing geometric parameters, lubrication parameters, and rotational speed. The bearing load includes the bearing load magnitude and bearing load direction. The bearing geometric parameters include the bearing diameter, pad width, number of pads, pad arc surface geometric parameters, pad swing angle, shaft center position and arrangement, clearance, and preload. The lubrication parameters include the lubricating oil type and inlet temperature. The initial values ​​of the static characteristics are obtained by jointly solving the Reynolds equation and the equilibrium equation.

[0011] Furthermore, before determining the initial values ​​of static characteristics based on bearing load, bearing geometric parameters, and lubrication parameters, the method further includes: Obtain compressor operating condition data, and determine the bearing load magnitude and bearing load direction based on the operating condition data.

[0012] According to another aspect of this application, a device for determining bearing dynamic parameters is provided, comprising: The module is used to determine the initial values ​​of static characteristics based on bearing load, bearing geometric parameters, lubrication parameters and rotational speed, and to construct a three-dimensional model of the bearing based on the initial values ​​of static characteristics and the bearing geometric parameters. The determination module is used to determine the bearing pad torque and journal resultant force deviation of the bearing based on the three-dimensional model of the bearing. The processing module is used to perform bearing dynamic characteristic analysis processing based on the small displacement eddy analysis method to obtain bearing dynamic parameters when the deviation between the bearing pad torque and the journal resultant force is less than a first preset deviation threshold and the deviation of static characteristic parameters is less than a second preset deviation threshold.

[0013] Furthermore, the device also includes: an adjustment module, The adjustment module is used to iteratively determine the bearing pad position, journal position, and static characteristic parameter deviation based on the adjusted shaft center position and the bearing pad swing angle when the deviation between the bearing pad torque and the journal resultant force is greater than or equal to a first preset deviation threshold. The determining module is further configured to redetermine the initial static characteristic value when the deviation between the bearing pad torque and the journal resultant force after iteration is less than the first preset deviation threshold, and the deviation of the static characteristic parameter after iteration is greater than or equal to the second preset deviation threshold, so as to re-execute the step of constructing a three-dimensional bearing model based on the initial static characteristic value and the bearing geometric parameters.

[0014] Furthermore, The processing module is specifically used to: acquire static characteristic parameters when the deviations of the bearing pad torque and the journal resultant force are both less than a first preset deviation threshold; and calculate the difference between the static characteristic parameters and the initial static characteristic values ​​to obtain the static characteristic parameter deviations; when the static characteristic parameter deviations are less than a second preset deviation threshold, perform bearing dynamic characteristic calculations, apply eddy displacement signals using a small displacement eddy analysis method, and collect oil film excitation forces in real time; and obtain bearing dynamic parameters based on data processing of the oil film excitation forces and eddy displacements.

[0015] Furthermore, The determining module is specifically used to determine the target bearing pad position and the target journal position based on the bearing three-dimensional model; to calculate the bearing pressure field characteristics of the target bearing pad position and the target journal position based on the fluid dynamics component, and to obtain the bearing pad torque and oil film resultant force; and to calculate the journal resultant force deviation based on the oil film resultant force and the bearing load.

[0016] Furthermore, The construction module is specifically used to start the three-dimensional modeling component and drive the three-dimensional modeling component to generate a three-dimensional model of the bearing based on the initial static characteristic value and the bearing geometric parameters; wherein, the three-dimensional model of the bearing is a three-dimensional oil film flow field model corresponding to the connected domain composed of the oil inlet, oil bladder area and bearing oil film area, the oil bladder area is formed between each bearing pad, and the bearing oil film area is formed between the bearing pad and the journal.

[0017] Furthermore, the construction module is specifically used to construct the Reynolds equation and the balance equation for the bearing load, the bearing geometric parameters, the lubrication parameters, and the rotational speed. The bearing load includes the bearing load magnitude and bearing load direction. The bearing geometric parameters include the bearing diameter, bearing pad width, number of bearing pads, bearing pad arc surface geometric parameters, bearing pad swing angle, shaft center position and arrangement, clearance, and preload. The lubrication parameters include the lubricating oil type and inlet temperature. The Reynolds equation and the balance equation are solved jointly to obtain the initial values ​​of the static characteristics.

[0018] Furthermore, the device also includes: The acquisition module is used to acquire the operating condition data of the compressor and determine the bearing load size and bearing load direction based on the operating condition data.

[0019] According to another aspect of this application, a storage medium is provided, wherein at least one executable instruction is stored therein, the executable instruction causing a processor to perform an operation corresponding to the above-described method for determining bearing dynamic parameters.

[0020] According to another aspect of this application, a terminal is provided, comprising: a processor, a memory, a communication interface, and a communication bus, wherein the processor, the memory, and the communication interface communicate with each other through the communication bus; The memory is used to store at least one executable instruction, which causes the processor to perform the operation corresponding to the method for determining the bearing dynamic parameters described above.

[0021] By employing the above technical solutions, the technical solutions provided in the embodiments of this application have at least the following advantages: This application provides a method, apparatus, storage medium, and terminal for determining bearing dynamic parameters. Compared with the prior art, the embodiments of this application determine the initial value of static characteristics based on bearing load, bearing geometric parameters, lubrication parameters, and rotational speed, and construct a three-dimensional bearing model based on the initial value of static characteristics and the bearing geometric parameters; determine the bearing pad torque and journal resultant force deviation of the bearing based on the three-dimensional bearing model; when the bearing pad torque and the journal resultant force deviation are both less than a first preset deviation threshold, and the static characteristic parameter deviation is less than a second preset deviation threshold, perform bearing dynamic characteristic analysis processing based on the small displacement eddy analysis method to obtain the bearing dynamic parameters. This realizes the calculation of bearing dynamic parameters by iteratively correcting the static characteristic parameters in combination with the three-dimensional model, which meets the calculation accuracy of bearing dynamic parameters of different load sizes and directions, thereby improving the effectiveness of bearing dynamic parameter calculation.

[0022] The above description is only an overview of the technical solution of this application. In order to better understand the technical means of this application and to implement it in accordance with the contents of the specification, and to make the above and other objects, features and advantages of this application more obvious and understandable, specific embodiments of this application are given below. Attached Figure Description

[0023] Various other advantages and benefits will become apparent to those skilled in the art upon reading the following detailed description of preferred embodiments. The accompanying drawings are for illustrative purposes only and are not intended to limit the scope of this application. Furthermore, the same reference numerals denote the same parts throughout the drawings. In the drawings: Figure 1 A flowchart illustrating a method for determining bearing dynamic parameters according to an embodiment of this application is shown; Figure 2 A flowchart illustrating another method for determining bearing dynamic parameters provided in an embodiment of this application is shown; Figure 3 This paper shows a schematic diagram of the distribution of different regions of an oil film flow field according to an embodiment of this application; Figure 4 This illustration shows a schematic diagram of an oil film flow field structure provided in an embodiment of this application; Figure 5 A schematic diagram of an ideal bearing and shaft diameter structure provided in an embodiment of this application is shown; Figure 6 This illustration shows a schematic diagram of a bearing and shaft diameter structure under load according to an embodiment of this application; Figure 7 This illustration shows a schematic diagram of a three-dimensional bearing model provided in an embodiment of this application; Figure 8 A block diagram of a bearing dynamic parameter determination device provided in an embodiment of this application is shown; Figure 9 A schematic diagram of the structure of a terminal provided in an embodiment of this application is shown. Detailed Implementation

[0024] Exemplary embodiments of the present disclosure will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.

[0025] This application provides a method for determining bearing dynamic parameters, such as... Figure 1 As shown, the method includes: 101. Determine the initial values ​​of static characteristics based on bearing load, bearing geometric parameters, lubrication parameters, and rotational speed, and construct a three-dimensional model of the bearing based on the initial values ​​of static characteristics and the bearing geometric parameters.

[0026] In this embodiment, the bearing can be a full-flow-field tilting pad bearing in a compressor. The current execution end, as the main body for analyzing the dynamic characteristics of the bearing, can be a terminal processor or a remote server; this embodiment does not impose any specific limitations. The bearing load includes the bearing load magnitude and bearing load direction. The bearing geometric parameters include the bearing diameter, pad width, number of pads, pad arc surface geometric parameters, pad tilt angle, shaft center position and arrangement, clearance, and preload. The lubrication parameters include the lubricating oil type and inlet temperature. The bearing load, bearing geometric parameters, lubrication parameters, and rotational speed (rotor speed) can be acquired and calculated using compressor engineering software. Furthermore, the bearing's static characteristics are used to characterize features under steady-state conditions, including eccentricity, offset angle, fulcrum oil film thickness, and pad tilt angle relative to the fulcrum. The initial static characteristic values ​​are calculated based on the bearing load, bearing geometric parameters, lubrication parameters, and rotational speed. The calculation method can be to construct mathematical equations through mechanical analysis of the bearing pads to obtain the initial static characteristic parameter values; this embodiment does not impose any other limitations.

[0027] It should be noted that, in order to perform characteristic analysis using model calibration, the current execution end constructs a three-dimensional model of the bearing based on the initial static characteristic values ​​and bearing geometric parameters. At this point, the three-dimensional bearing model is a three-dimensional model describing the spatial position of components such as the bearing bearing pads. That is, due to the high-speed rotation of the rotor and the complex structure of the tilting pad bearing, the flow of fluid inside the bearing is a very complex three-dimensional flow. Furthermore, in the Computational Fluid Dynamics (CFD) software used to analyze the dynamic characteristics of the bearing, in order to obtain the oil film excitation force on the shaft diameter surface, a full three-dimensional model of the tilting pad bearing with a 360° fluid domain can be established. This application embodiment does not impose specific limitations on this. Additionally, the compressor in this application embodiment can be a compressor in a single-shaft centrifugal compressor unit, and the bearing can include, but is not limited to, different types of support bearings, such as round bearings, elliptical bearings, and multi-blade bearings. This application embodiment does not impose specific limitations on this.

[0028] 102. Based on the three-dimensional model of the bearing, determine the bearing pad torque and journal resultant force deviation of the bearing.

[0029] In this embodiment, the bearing pad torque is the torque of the bearing pad, and the journal resultant force deviation is the difference between the resultant force of the journal and the load in the bearing, so as to characterize the force situation in the oil film flow field in the bearing.

[0030] It should be noted that the bearing 3D model is constructed using 3D modeling components, which can perform force analysis and spatial distance analysis for different spatial locations. Therefore, relevant force data can be obtained from the bearing 3D model to determine the bearing pad torque and journal resultant force deviation. This application embodiment does not make specific limitations.

[0031] 103. When the deviations of the bearing pad torque and the journal resultant force are both less than the first preset deviation threshold, and the deviation of the static characteristic parameters is less than the second preset deviation threshold, the bearing dynamic characteristic analysis is performed based on the small displacement eddy analysis method to obtain the bearing dynamic parameters.

[0032] In this embodiment, after obtaining the bearing pad torque and bearing resultant force deviation, the current execution end compares them with a first preset deviation threshold to determine whether the bearing is in a balanced position. The first preset deviation threshold is any small integer, preferably 0.01, but this embodiment does not impose a specific limitation. Simultaneously, when both the bearing pad torque and journal resultant force deviation are less than the first preset deviation threshold, it is determined whether the static characteristic parameter deviation is less than a second preset deviation threshold. In this case, the static characteristic parameter deviation characterizes the difference between the initial static characteristic parameters and the static characteristic parameters obtained through iterative calculation. The static characteristic parameters may include optimal parameters selected for the execution engineering project, such as heart rate, offset angle, fulcrum oil film thickness, and pad tilt angle relative to the fulcrum; this embodiment does not impose a specific limitation. The second preset deviation threshold can be any small integer, preferably 0.01. When the deviations of the bearing pad torque and the journal resultant force are both less than the first preset deviation threshold, and the deviation of the static characteristic parameters is less than the second preset deviation threshold, the current execution end sets the eddy displacement signal based on the small displacement eddy analysis method to perform bearing dynamic characteristic analysis and obtain the bearing dynamic parameters. That is, the bearing dynamic parameters can be obtained by analysis using computational fluid dynamics software CFD. At this time, the bearing dynamic parameters are used to characterize the relevant characteristic parameters of bearing dynamics. The small displacement eddy theory can be used, and the transient calculation method of CFD software can be used to calculate the bearing dynamic characteristic parameters. This application embodiment does not limit this.

[0033] In another embodiment of this application, for further definition and explanation, the steps also include: When the deviation between the bearing pad torque and the journal resultant force is greater than or equal to a first preset deviation threshold, the bearing pad position, journal position, and static characteristic parameter deviation are iteratively determined based on the adjusted shaft center position and the bearing pad swing angle. When the deviations of the bearing pad torque and the journal resultant force after iteration are both less than the first preset deviation threshold, and the deviations of the static characteristic parameters after iteration are greater than or equal to the second preset deviation threshold, the initial static characteristic value is re-determined, so as to re-execute the step of constructing a three-dimensional bearing model based on the initial static characteristic value and the bearing geometric parameters.

[0034] To improve the accuracy of dynamic characteristic parameter calculations and ensure the effectiveness of data modeled at the static equilibrium position for precise calculation of dynamic characteristic parameters, the current execution end first compares the bearing pad torque and journal resultant force deviation with a first preset deviation threshold. After both the bearing pad torque and the journal resultant force deviation are less than the first preset deviation threshold, static characteristic parameters are extracted and subtracted from their initial values ​​to obtain the static characteristic parameter deviation, which is then compared with a second preset deviation threshold. If the bearing pad torque and journal resultant force deviation are greater than or equal to the first preset deviation threshold, the pad position, journal position, and static characteristic parameter deviation are iteratively determined based on the adjusted shaft center position and the pad swing angle. Figure 2 As shown, after readjusting the shaft center position and bearing pad swing angle, the static characteristic parameters are recalculated using the bearing 3D model. This leads to the next iteration's step comparing the bearing pad torque and journal resultant force deviation with a first preset deviation threshold. If the deviations of the bearing pad torque and journal resultant force after iteration are both less than the first preset deviation threshold, and the deviations of the static characteristic parameters after iteration are greater than or equal to a second preset deviation threshold, the initial static characteristic values ​​are redefined. This allows for the re-execution of the step of constructing the bearing 3D model based on the initial static characteristic values ​​and bearing geometric parameters, thus entering the next iteration. If the deviations of the bearing pad torque and journal resultant force after iteration are both less than the first preset deviation threshold, and the deviations of the static characteristic parameters after iteration are less than the second preset deviation threshold, then the step of performing bearing dynamic characteristic analysis by setting a whirl displacement signal to obtain the bearing dynamic parameters is executed.

[0035] In some embodiments, the iterative process following the adjustment of the shaft center position and tile swing angle can begin by iterating over the shaft center position, adjusting the eccentricity and offset angle. This causes a deviation in the torque acting on the tile. Therefore, the tile swing angle is iterated over. When all tile torques meet the aforementioned conditions, the resultant force on the shaft diameter is compared. If the aforementioned conditions are not met, the shaft center position can be iterated over again until the static characteristic parameter deviation is satisfied. The iterative formula is expressed as: ; ; in, The axis position, and The iteration step size, The number of iterations. For the corners of the tiles, For the oil film's combined force, This is to calculate the bearing torque so that the corresponding oil film resultant force and bearing torque can be calculated after each adjustment of the bearing swing angle and shaft center position. Additionally, in this embodiment, the journal resultant force deviation is expressed as ΔF = |F p -W|, where W represents the bearing load.

[0036] It should be noted that the adjustment of the shaft position and the bearing swing angle can be made manually through the bearing 3D model, or the adjustment coefficient can be preset and adjusted in the bearing 3D model. This application embodiment does not make specific limitations.

[0037] In another embodiment of this application, for further definition and explanation, when the deviations of the bearing pad torque and the journal resultant force are both less than a first preset deviation threshold, and the deviation of the static characteristic parameters is less than a second preset deviation threshold, the bearing dynamic characteristic analysis is performed based on the small displacement eddy analysis method to obtain the bearing dynamic parameters including: When the deviation between the bearing pad torque and the journal resultant force is less than a first preset deviation threshold, static characteristic parameters are obtained, and the difference between the static characteristic parameters and the initial static characteristic value is calculated to obtain the static characteristic parameter deviation. When the deviation of the static characteristic parameter is less than the second preset deviation threshold, the dynamic characteristic calculation of the bearing is carried out. The small displacement eddy analysis method is adopted, the eddy displacement signal is applied, and the oil film excitation force is collected in real time. Based on data processing of the oil film excitation force and eddy displacement, the bearing dynamic parameters are obtained.

[0038] To achieve the goal of calculating bearing dynamic parameters based on iteratively corrected loads, thereby improving the effectiveness and accuracy of determining bearing dynamic characteristic parameters, the current execution end compares the relationship between the bearing pad torque, journal resultant force deviation, and the first preset deviation threshold, as follows: Figure 2 As shown, when the deviations of the bearing pad torque and the journal resultant force are both less than the first preset deviation threshold, static characteristic parameters are obtained, and the difference between the static characteristic parameters and the initial static characteristic values ​​is calculated to obtain the static characteristic parameter deviation. At this time, the static characteristic parameters can be pre-configured, and this application embodiment does not impose specific limitations. Furthermore, when the static characteristic parameter deviation is less than the second preset deviation threshold, it indicates that the bearing's stable operating state meets engineering requirements. Therefore, based on the small displacement eddy analysis method, an eddy displacement signal is set to perform bearing dynamic characteristic analysis processing to obtain the bearing dynamic parameters. That is, the bearing dynamic parameters can be obtained by analysis using computational fluid dynamics (CFD) software, and this application embodiment does not impose specific limitations. In addition, the bearing dynamic parameters in this application embodiment represent relevant parameters generated when the bearing oil film rotates during compressor operation, including but not limited to dynamic bearing stiffness and damping coefficients, and this application embodiment does not impose specific limitations.

[0039] In some embodiments, the resultant force Fp of the oil film on the bearing pad surface and the torque M on the pad, when the torque M on the pad is less than a first preset deviation threshold δ1, the resultant force deviation ΔF on the journal is |F|. p -W|<δ1 indicates that the equilibrium position has been reached; otherwise, iterative calculations are performed. When the static characteristic deviation is less than the second preset deviation threshold δ2, an eddy displacement signal is applied, and bearing CFD analysis is carried out to extract the oil film excitation force and analyze the bearing dynamic characteristic coefficients.

[0040] In another embodiment of this application, for further definition and explanation, the step of determining the bearing pad torque and journal resultant force deviation based on the bearing three-dimensional model includes: The target bearing pad position and the target journal position are determined based on the three-dimensional model of the bearing. Based on the fluid dynamics components, the bearing pressure field characteristics of the target bearing pad position and the target journal position are calculated to obtain the bearing pad torque and the oil film resultant force. The journal resultant force deviation is obtained by calculating based on the oil film resultant force and bearing load.

[0041] To calculate the bearing pad torque and journal resultant force deviation using a 3D bearing model, thereby characterizing the bearing's steady-state stress and improving the accuracy of using bearing pad torque and journal resultant force deviation as iterative judgment criteria, the 3D bearing model is constructed based on initial static characteristics and bearing geometric parameters using a 3D modeling component. The target pad and journal positions are determined using the 3D bearing model. Then, the bearing pressure field characteristics at the target pad and journal positions are calculated using CFD analysis with a fluid dynamics component, yielding the bearing pad torque and oil film resultant force. This involves pressure integration calculations on the inner surfaces of each pad to extract the oil film resultant force F on the bearing pad surface. p And the bearing moment M. Simultaneously, based on the resultant force of the oil film and the bearing load M, the journal resultant force deviation is calculated, expressed as ΔF = |F p -W|.

[0042] In another embodiment of this application, for further definition and explanation, the step of constructing a three-dimensional bearing model based on the initial static characteristic value and the bearing geometric parameters includes: The 3D modeling component is activated, and based on the initial static characteristics and the bearing geometric parameters, the 3D modeling component is driven to generate a 3D model of the bearing.

[0043] To improve the effectiveness of bearing dynamic characteristic parameter calculation through iterative correction, the current execution end constructs a 3D bearing model based on initial static characteristic values ​​(initial set bearing pad positions and shaft diameter positions) and bearing geometric parameters (bearing diameter, bearing pad width, bearing pad wrap angle, etc.). Specifically, the 3D modeling component is first activated; however, this embodiment does not specify a particular 3D modeling component. The generated bearing 3D model is a 3D oil film flow field model corresponding to the connected domain composed of the oil inlet, oil bladder area, and bearing oil film area. The oil bladder area is formed between the bearing pads, and the bearing oil film area is formed between the bearing pad and the journal. Figure 3 As shown, the oil inlet 1, the oil sac area 2 between the bearings, and the bearing oil film area 3 between the bearing and the journal are interconnected, forming the entire oil film flow field. Figure 4 As shown, the solid structure of tile 5 does not need to be modeled. In this case, the shaded area represents the flow field obtained after modeling.

[0044] In some embodiments, such as Figure 5The ideal bearing and shaft diameter shown are illustrated. Under load, the shaft diameter exhibits eccentricity and misalignment angle, the bearing pads have a swing angle relative to the fulcrum, and a lubricating oil film exists between the pads and the shaft diameter. Due to differences in bearing structure, the oil inlet area varies. This embodiment uses an oil inlet designed between two adjacent pads as an example for modeling, constructing a three-dimensional bearing model of the oil film domain, i.e., the lubricating oil flow field region. The model does not include the solid bearing region. Figure 6 As shown, the static equilibrium position circle of the shaft diameter is established based on the eccentricity and offset angle calculated from the aforementioned initial static characteristics, i.e., as shown in the figure. Figure 6 The toroidal surface centered at 01 shown can be used to create the inner surface of the bearing pads based on the pad angle, pad size, and arrangement. Simultaneously, the bearing oil inlet structure is simplified to create oil inlet channels and inlets between the pads, resulting in the final 3D bearing model. Figure 7 As shown.

[0045] In another embodiment of this application, for further definition and explanation, the step of determining the initial value of the static characteristics based on the bearing load, bearing geometric parameters, lubrication parameters, and rotational speed includes: Construct the Reynolds equations and balance equations for the bearing load, bearing geometry, lubrication parameters, and rotational speed; The initial values ​​of the static characteristics are obtained by jointly solving the Reynolds equation and the equilibrium equation.

[0046] To analyze the dynamic characteristics of assembled compressor bearings under different loads by combining engineering algorithms and CFD technology, and to achieve the goal of calculating bearing dynamic characteristic parameters using iterative correction, the current execution end calculates the initial values ​​of static characteristics. Specifically, the operator can use operating conditions such as gravity, gear meshing force, and bearing load as input parameters for the current execution end. Then, using the input load, bearing geometric parameters, lubrication parameters, and rotational speed, Reynolds equations and equilibrium equations are constructed. By simultaneously solving the Reynolds equations and equilibrium equations (forces and torques), the initial values ​​of the bearing's static characteristics are obtained, including but not limited to initial eccentricity, offset angle, bearing oil film thickness, and bearing pad tilt angle relative to the bearing. This application does not specifically limit these values ​​in its embodiments.

[0047] In some embodiments, the constructed Reynolds equation is expressed as: ; in, The circumferential position of the axis, The axial position of the axis. For oil film pressure, For oil film thickness, For linear velocity, The turbulence coefficient is... , For preload, For -eccentricity, It is an off-center angle. Let the arc length of the tile be _____. The time-shifting angle of the tile. The angle represents the position of the tile's fulcrum.

[0048] In some embodiments, the constructed force balance equations are expressed as: ,in, For the oil film force of the i-th bearing, For external load.

[0049] The torque balance equation is expressed as: ;in, The distance from the pressure point to the fulcrum. Let be the integration domain, and be the surface of the i-th tile.

[0050] It should be noted that, in the embodiments of this application, the bearing load includes the bearing load size and bearing load direction, the bearing geometric parameters include the bearing diameter, pad width, number of pads, pad arc surface geometric parameters, pad swing angle, shaft center position and arrangement, as well as clearance and preload, and the lubrication parameters include the lubricating oil type and inlet temperature. The above parameters can be entered by the operator as input parameters for solving equations simultaneously. The embodiments of this application do not impose specific limitations.

[0051] In another embodiment of this application, for further definition and explanation, before determining the initial values ​​of the static characteristics based on the bearing load, bearing geometric parameters, and lubrication parameters, the method further includes: Obtain compressor operating condition data, and determine the bearing load magnitude and bearing load direction based on the operating condition data.

[0052] In order to meet the power characteristic analysis of different loads, the embodiments of this application also need to obtain the operating condition data of the compressor so as to determine the load and bearing direction of the bearing based on the engineering condition data. At this time, the embodiments of this application can perform calculations based on existing calculation components such as gear meshing force. The embodiments of this application do not make specific limitations.

[0053] This application provides a method for determining bearing dynamic parameters. The method determines initial static characteristic values ​​based on bearing load, bearing geometric parameters, lubrication parameters, and rotational speed. A three-dimensional bearing model is then constructed based on these initial static characteristic values ​​and the bearing geometric parameters. The bearing pad torque and journal resultant force deviation are determined based on the three-dimensional bearing model. When both the bearing pad torque and the journal resultant force deviation are less than a first preset deviation threshold, and the static characteristic parameter deviation is less than a second preset deviation threshold, the bearing dynamic characteristics are analyzed using a small displacement eddy analysis method to obtain the bearing dynamic parameters. This method calculates bearing dynamic parameters by iteratively correcting the static characteristic parameters using a three-dimensional model, ensuring the accuracy of bearing dynamic parameter calculations for different load sizes and directions, thereby improving the effectiveness of bearing dynamic parameter calculation.

[0054] Furthermore, as a response to the above Figure 1 The implementation of the method shown in this application provides a device for determining bearing dynamic parameters, such as... Figure 8 As shown, the device includes: The construction module 21 is used to determine the initial value of static characteristics based on the bearing load, bearing geometric parameters, lubrication parameters and rotational speed, and to construct a three-dimensional model of the bearing based on the initial value of static characteristics and the bearing geometric parameters. The determination module 22 is used to determine the bearing pad torque and journal resultant force deviation of the bearing based on the three-dimensional model of the bearing. The processing module 23 is used to perform bearing dynamic characteristic analysis processing based on the small displacement eddy analysis method to obtain bearing dynamic parameters when the deviation between the bearing pad torque and the journal resultant force is less than a first preset deviation threshold and the deviation of static characteristic parameters is less than a second preset deviation threshold.

[0055] Furthermore, the device also includes: an adjustment module, The adjustment module is used to iteratively determine the bearing pad position, journal position, and static characteristic parameter deviation based on the adjusted shaft center position and the bearing pad swing angle when the deviation between the bearing pad torque and the journal resultant force is greater than or equal to a first preset deviation threshold. The determining module is further configured to redetermine the initial static characteristic value when the deviation between the bearing pad torque and the journal resultant force after iteration is less than the first preset deviation threshold, and the deviation of the static characteristic parameter after iteration is greater than or equal to the second preset deviation threshold, so as to re-execute the step of constructing a three-dimensional bearing model based on the initial static characteristic value and the bearing geometric parameters.

[0056] Furthermore, The processing module is specifically used to: acquire static characteristic parameters when the deviations of the bearing pad torque and the journal resultant force are both less than a first preset deviation threshold; and calculate the difference between the static characteristic parameters and the initial static characteristic values ​​to obtain the static characteristic parameter deviations; when the static characteristic parameter deviations are less than a second preset deviation threshold, perform bearing dynamic characteristic calculations, apply eddy displacement signals using a small displacement eddy analysis method, and collect oil film excitation forces in real time; and obtain bearing dynamic parameters based on data processing of the oil film excitation forces and eddy displacements.

[0057] Furthermore, The determining module is specifically used to determine the target bearing pad position and the target journal position based on the bearing three-dimensional model; to calculate the bearing pressure field characteristics of the target bearing pad position and the target journal position based on the fluid dynamics component, and to obtain the bearing pad torque and oil film resultant force; and to calculate the journal resultant force deviation based on the oil film resultant force and the bearing load.

[0058] Furthermore, The construction module is specifically used to start the three-dimensional modeling component and drive the three-dimensional modeling component to generate a three-dimensional model of the bearing based on the initial static characteristic value and the bearing geometric parameters; wherein, the three-dimensional model of the bearing is a three-dimensional oil film flow field model corresponding to the connected domain composed of the oil inlet, oil bladder area and bearing oil film area, the oil bladder area is formed between each bearing pad, and the bearing oil film area is formed between the bearing pad and the journal.

[0059] Furthermore, the construction module is specifically used to construct the Reynolds equation and the balance equation for the bearing load, the bearing geometric parameters, the lubrication parameters, and the rotational speed. The bearing load includes the bearing load magnitude and bearing load direction. The bearing geometric parameters include the bearing diameter, bearing pad width, number of bearing pads, bearing pad arc surface geometric parameters, bearing pad swing angle, shaft center position and arrangement, clearance, and preload. The lubrication parameters include the lubricating oil type and inlet temperature. The Reynolds equation and the balance equation are solved jointly to obtain the initial values ​​of the static characteristics.

[0060] Furthermore, the device also includes: The acquisition module is used to acquire the operating condition data of the compressor and determine the bearing load size and bearing load direction based on the operating condition data.

[0061] This application provides a device for determining bearing dynamic parameters. This device determines initial static characteristic values ​​based on bearing load, bearing geometric parameters, lubrication parameters, and rotational speed, and constructs a three-dimensional bearing model based on these initial static characteristic values ​​and the bearing geometric parameters. Based on the three-dimensional bearing model, it determines the bearing pad torque and journal resultant force deviation. When both the bearing pad torque and the journal resultant force deviation are less than a first preset deviation threshold, and the static characteristic parameter deviation is less than a second preset deviation threshold, it performs bearing dynamic characteristic analysis using a small displacement eddy analysis method to obtain the bearing dynamic parameters. This achieves the calculation of bearing dynamic parameters by iteratively correcting the static characteristic parameters using a three-dimensional model, ensuring the accuracy of bearing dynamic parameter calculations for different load sizes and directions, thereby improving the effectiveness of bearing dynamic parameter calculation.

[0062] According to one embodiment of this application, a storage medium is provided, the storage medium storing at least one executable instruction, the computer-executable instruction being able to execute the method for determining bearing dynamic parameters in any of the above method embodiments.

[0063] Figure 9 The diagram shows a structural schematic of a terminal according to one embodiment of the present application. The specific embodiments of the present application do not limit the specific implementation of the terminal.

[0064] like Figure 9 As shown, the terminal may include: a processor 302, a communications interface 304, a memory 306, and a communications bus 308.

[0065] The processor 302, communication interface 304, and memory 306 communicate with each other via communication bus 308.

[0066] Communication interface 304 is used to communicate with other network elements such as clients or other servers.

[0067] The processor 302 is used to execute program 310, specifically the relevant steps in the above-described method embodiment for determining bearing dynamic parameters.

[0068] Specifically, program 310 may include program code that includes computer operation instructions.

[0069] Processor 302 may be a central processing unit (CPU), an application-specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of this application. The terminal includes one or more processors, which may be processors of the same type, such as one or more CPUs; or they may be processors of different types, such as one or more CPUs and one or more ASICs.

[0070] Memory 306 is used to store program 310. Memory 306 may include high-speed RAM memory, and may also include non-volatile memory, such as at least one disk storage device.

[0071] Specifically, program 310 can be used to cause processor 302 to perform the following operations: The initial values ​​of static characteristics are determined based on the bearing load, bearing geometric parameters, lubrication parameters, and rotational speed, and a three-dimensional model of the bearing is constructed based on the initial values ​​of static characteristics and the bearing geometric parameters. The bearing pad torque and journal resultant force deviation of the bearing are determined based on the three-dimensional model of the bearing. When the deviations between the bearing pad torque and the journal resultant force are both less than the first preset deviation threshold, and the deviation of the static characteristic parameters is less than the second preset deviation threshold, the bearing dynamic characteristic analysis is performed based on the small displacement eddy analysis method to obtain the bearing dynamic parameters.

[0072] Obviously, those skilled in the art should understand that the modules or steps of this application described above can be implemented using general-purpose computing devices. They can be centralized on a single computing device or distributed across a network of multiple computing devices. Optionally, they can be implemented using computer-executable program code, thereby storing them in a storage device for execution by a computing device. In some cases, the steps shown or described can be performed in a different order than those presented here, or they can be fabricated as separate integrated circuit modules, or multiple modules or steps can be fabricated as a single integrated circuit module. Thus, this application is not limited to any particular combination of hardware and software.

[0073] The above description is merely a preferred embodiment of this application and is not intended to limit this application. Various modifications and variations can be made to this application by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the scope of protection of this application.

Claims

1. A method for determining bearing dynamic parameters, characterized in that, include: The initial values ​​of static characteristics are determined based on the bearing load, bearing geometric parameters, lubrication parameters, and rotational speed, and a three-dimensional model of the bearing is constructed based on the initial values ​​of static characteristics and the bearing geometric parameters. The bearing pad torque and journal resultant force deviation of the bearing are determined based on the three-dimensional model of the bearing. When the deviations between the bearing pad torque and the journal resultant force are both less than the first preset deviation threshold, and the deviation of the static characteristic parameters is less than the second preset deviation threshold, the bearing dynamic characteristic analysis is performed based on the small displacement eddy analysis method to obtain the bearing dynamic parameters.

2. The method according to claim 1, characterized in that, The method further includes: When the deviation between the bearing pad torque and the journal resultant force is greater than or equal to a first preset deviation threshold, the bearing pad position, journal position, and static characteristic parameter deviation are iteratively determined based on the adjusted shaft center position and the bearing pad swing angle. When the deviations of the bearing pad torque and the journal resultant force after iteration are both less than the first preset deviation threshold, and the deviations of the static characteristic parameters after iteration are greater than or equal to the second preset deviation threshold, the initial static characteristic value is re-determined, so as to re-execute the step of constructing a three-dimensional bearing model based on the initial static characteristic value and the bearing geometric parameters.

3. The method according to claim 1, characterized in that, When the deviations between the bearing pad torque and the journal resultant force are both less than a first preset deviation threshold, and the deviations in static characteristic parameters are less than a second preset deviation threshold, the bearing dynamic characteristic analysis is performed based on the small displacement eddy analysis method to obtain the bearing dynamic parameters, including: When the deviation between the bearing pad torque and the journal resultant force is less than a first preset deviation threshold, static characteristic parameters are obtained, and the difference between the static characteristic parameters and the initial static characteristic value is calculated to obtain the static characteristic parameter deviation. When the deviation of the static characteristic parameter is less than the second preset deviation threshold, the dynamic characteristic calculation of the bearing is carried out. The small displacement eddy analysis method is adopted, the eddy displacement signal is applied, and the oil film excitation force is collected in real time. Based on data processing of the oil film excitation force and eddy displacement, the bearing dynamic parameters are obtained.

4. The method according to claim 1, characterized in that, The determination of the bearing pad torque and journal resultant force deviation based on the bearing three-dimensional model includes: The target bearing pad position and the target journal position are determined based on the three-dimensional model of the bearing. Based on the fluid dynamics components, the bearing pressure field characteristics of the target bearing pad position and the target journal position are calculated to obtain the bearing pad torque and the oil film resultant force. The journal resultant force deviation is obtained by calculating based on the oil film resultant force and bearing load.

5. The method according to claim 1, characterized in that, The construction of the bearing 3D model based on the initial static characteristic value and the bearing geometric parameters includes: The 3D modeling component is activated, and based on the initial static characteristic values ​​and the bearing geometric parameters, the 3D modeling component is driven to generate a 3D model of the bearing; The bearing three-dimensional model is a three-dimensional oil film flow field model corresponding to the connected domain composed of the oil inlet, oil bladder area and bearing oil film area. The oil bladder area is formed between each bearing pad, and the bearing oil film area is formed between the bearing pad and the journal.

6. The method according to claim 1, characterized in that, The determination of initial static characteristic values ​​based on bearing load, bearing geometric parameters, lubrication parameters, and rotational speed includes: Construct the Reynolds equation and balance equation for the bearing load, bearing geometric parameters, lubrication parameters, and rotational speed. The bearing load includes the bearing load magnitude and bearing load direction. The bearing geometric parameters include the bearing diameter, pad width, number of pads, pad arc surface geometric parameters, pad swing angle, shaft center position and arrangement, clearance, and preload. The lubrication parameters include the lubricating oil type and inlet temperature. The initial values ​​of the static characteristics are obtained by jointly solving the Reynolds equation and the equilibrium equation.

7. The method according to claim 6, characterized in that, Before determining the initial values ​​of static characteristics based on bearing load, bearing geometric parameters, and lubrication parameters, the method further includes: Obtain compressor operating condition data, and determine the bearing load magnitude and bearing load direction based on the operating condition data.

8. A device for determining bearing dynamic parameters, characterized in that, include: The module is used to determine the initial values ​​of static characteristics based on bearing load, bearing geometric parameters, lubrication parameters and rotational speed, and to construct a three-dimensional model of the bearing based on the initial values ​​of static characteristics and the bearing geometric parameters. The determination module is used to determine the bearing pad torque and journal resultant force deviation of the bearing based on the three-dimensional model of the bearing. The processing module is used to perform bearing dynamic characteristic analysis processing based on the small displacement eddy analysis method to obtain bearing dynamic parameters when the deviation between the bearing pad torque and the journal resultant force is less than a first preset deviation threshold and the deviation of static characteristic parameters is less than a second preset deviation threshold.

9. A computer-readable storage medium having a computer program / instructions stored thereon, characterized in that, When the computer program / instructions are executed by the processor, they implement the steps of the method of claim 1.

10. A terminal, comprising a memory, a processor, and a computer program stored in the memory, characterized in that, The processor executes the computer program to implement the steps of the method of claim 1.

Images