A method and system for calculating the oil film force of an extrusion oil film damper

Abstract

This application discloses a method and system for calculating the oil film force of an extrusion oil film damper. The method includes: dimensionlessly transforming the oil film pressure control equation to obtain a dimensionless oil film pressure control equation; meshing the oil film region of the extrusion oil film damper to obtain multiple oil film mesh elements; performing a weak-form transformation of the dimensionless oil film pressure control equation using the finite element method based on the oil film region to obtain a weak-form oil film pressure control equation; determining the element stiffness matrix and element load vector corresponding to each oil film mesh element; determining the global stiffness matrix and global load vector by constructing an element-global mapping matrix; and determining the oil film force of the extrusion oil film damper based on the global stiffness matrix and global load vector, thereby improving the calculation efficiency of the oil film force.

Description

Technical Field

[0001] This application relates to the technical field of oil film force calculation for extrusion oil film dampers, and in particular to a method and system for calculating the oil film force of extrusion oil film dampers. Background Technology

[0002] Rotor-support systems are widely used in modern industrial production. As rotors develop towards higher speeds, greater compactness (smaller structures), and higher parameters (high speed, high pressure, large flow), the integrity and reliability of their structures become increasingly important. Vibration and stability of rotor-support systems are often the root causes of these problems.

[0003] Squeeze film dampers (SFDs) suppress vibration by adding an external damping term to the shaft support system. For elastically supported rotor systems, a well-designed SFD can effectively reduce the external excitation of the rotor, significantly lower the vibration peak value when the rotor is overcritical, and thus reduce the impact on the support structure, making it an effective and economical technical means to solve the vibration and noise problems of high-speed rotating machinery. Therefore, the calculation of the SFD film force is particularly important for the design of the SFD-rotor system and the analysis of the system's dynamic response.

[0004] Currently, the SFD oil film force is mainly calculated using the over-relaxation iterative method. However, the over-relaxation iterative method can only handle working conditions with regular boundary conditions, and it usually requires hundreds of iterations to meet the convergence tolerance requirements, resulting in low computational efficiency. Summary of the Invention

[0005] This application aims to at least address the technical problems existing in the prior art. To this end, this application proposes a method and system for calculating the oil film force of an extrusion oil film damper, which is unaffected by boundary conditions and improves the calculation efficiency of the oil film force.

[0006] A first aspect of this application provides a method for calculating the oil film force of a squeeze oil film damper, comprising the following steps: In constructing the oil film pressure control equation, the oil film pressure control equation is dimensionless to obtain the dimensionless oil film pressure control equation. Having obtained the oil film region of the extrusion oil film damper, the oil film region is divided into grids to obtain multiple oil film grid units; Based on the oil film region, the dimensionless oil film pressure control equation is transformed into a weak form using the finite element method to obtain the weak form of the oil film pressure control equation. Based on the weak form of the oil film pressure control equation, the element stiffness matrix and element load vector corresponding to each oil film mesh element are determined. In constructing the element-global mapping matrix, the global stiffness matrix is ​​determined based on the element-global mapping matrix and all the element stiffness matrices; the global load vector is determined based on the element-global mapping matrix and all the element load vectors. Under the condition of applying a preset boundary condition, the oil film force of the extrusion oil film damper is determined based on the global stiffness matrix and the global load vector.

[0007] The oil film force calculation method for the extrusion oil film damper according to the embodiments of this application has at least the following beneficial effects: This application achieves dimensionless oil film pressure control equations by constructing them, obtaining dimensionless oil film pressure control equations; it then meshes the oil film region of the extrusion oil film damper to obtain multiple oil film mesh elements; based on the oil film region, it performs a weak-form transformation of the dimensionless oil film pressure control equations using the finite element method to obtain weak-form oil film pressure control equations; based on the weak-form oil film pressure control equations, it determines the element stiffness matrix and element load vector corresponding to each oil film mesh element; it determines the global stiffness matrix based on the element-global mapping matrix and all element stiffness matrices, and the global load vector based on the element-global mapping matrix and all element load vectors; and it determines the oil film force of the extrusion oil film damper based on the global stiffness matrix and global load vector, under the condition of applying preset boundary conditions. This achieves independence from the influence of boundary conditions and improves the calculation efficiency of oil film force.

[0008] According to some embodiments of this application, the oil film pressure control equation is dimensionless using the following formula to obtain the dimensionless oil film pressure control equation:

[0009] in, The oil film thickness at any dimensionless location. The pressure distribution value is the dimensionless value. This represents the horizontal displacement of the journal. This represents the vertical displacement of the journal. The oil film thickness at any location. Let be the radius of the outer wall surface of the journal. Let be the radius of the inner wall of the bearing housing. This represents the initial radial clearance of the oil film. This represents the pressure distribution value. For dimensionless parameters, For intermediate variable values, The angle between the journal and the horizontal axis in the counterclockwise direction. The x-axis velocity component is the dimensionless value. The velocity component in the y-direction after dimensionless scaling is denoted as . Here, θ is the rheological index, θ is the angle between the counterclockwise direction and the reference line, Z is the dimensionless axial coordinate, and L is the oil film length. This represents the dimensionless displacement of the journal in the horizontal direction. This represents the dimensionless displacement of the journal in the vertical direction. is the consistency coefficient of the fluid.

[0010] According to some embodiments of this application, based on the oil film region, the dimensionless oil film pressure control equation is transformed into a weak form using the finite element method through the following formula to obtain the weak form oil film pressure control equation:

[0011] in, This is the value of the first intermediate parameter. This is the value of the second intermediate parameter. This is the value of the third intermediate parameter. This is the value of the fourth intermediate parameter. For the oil film area, For the trial function, Let be the set of trial functions.

[0012] According to some embodiments of this application, determining the element stiffness matrix and element load vector corresponding to each oil film mesh element based on the weak form of the oil film pressure control equation includes: When the trial function is a 4-node bilinear function, through 2 2. The Gaussian integration method constructs the partial derivative at the Gaussian point corresponding to each oil film grid cell using the following formula:

[0013] in, Let be the shape function of the first node among the four nodes corresponding to the oil film mesh element. Let be the shape function of the second node out of the four nodes corresponding to the oil film mesh element. Let be the shape function of the third node out of the four nodes corresponding to the oil film mesh element. Let be the shape function of the fourth node out of the four nodes corresponding to the oil film mesh element. The x-axis of the reference unit is... The ordinate of the reference unit; Based on the partial derivatives at the Gaussian point corresponding to each oil film mesh element and the weak form of the oil film pressure control equation, the element stiffness matrix corresponding to each oil film mesh element is determined by the Picard iteration method and the following formula:

[0014] in, The two-dimensional integral is And the grid discretization function is The true value of stress, The two-dimensional integral is And the grid discretization function is The trial function within the oil film mesh element, The fourth node among the four nodes corresponding to the oil film mesh element The stress value of the node, The fourth node among the four nodes corresponding to the oil film mesh element The node's trial function, The two-dimensional integral is The The shape function of the node. The two-dimensional integral is The The value of the fifth intermediate parameter in the next iteration. For the first Pressure distribution value in the next iteration This is the set of four nodes corresponding to the oil film mesh element. For Jacobian matrices, For the first Oil film mesh unit Dimensional value in direction, For the first Oil film mesh unit Dimensional value in direction, For the first The first oil film mesh unit The element stiffness matrix of the next iteration The two-dimensional integral is The The value of the fifth intermediate parameter in the next iteration. For the oil film area corresponding to the oil film grid cell, For transpose; The element load vector corresponding to each oil film mesh element is determined by the following formula:

[0015] in, The fourth node among the four nodes corresponding to the oil film mesh element The weight value of the node. The fourth node among the four nodes corresponding to the oil film mesh element The weight value of the node. For the first The first oil film mesh unit The element load vector of the next iteration. The two-dimensional integral is The The fourth intermediate parameter of the node.

[0016] According to some embodiments of this application, the global stiffness matrix is ​​determined by the following formula:

[0017] in, For the first The global stiffness matrix of the next iteration. For the first The cell-global mapping matrix of each oil film mesh cell; The global load vector is determined using the following formula:

[0018] in, For the first The global load vector for the next iteration.

[0019] According to some embodiments of this application, the oil film force of the extrusion oil film damper includes the oil film reaction force in the horizontal direction under the Cartesian coordinate system and the oil film reaction force in the vertical direction under the Cartesian coordinate system. Determining the oil film force of the extrusion oil film damper based on the global stiffness matrix and the global load vector includes: Based on the global stiffness matrix and the global load vector, the first... Pressure distribution value in the next iteration:

[0020] in, For the first Pressure distribution values ​​in the next iteration; In the first If the pressure distribution value of the next iteration meets the preset convergence condition, the oil film force of the extrusion oil film damper is calculated using the following formula:

[0021] in, This represents the true pressure value. For the oil film region of the e-th oil film unit, The two-dimensional integral is The true value of stress, The two-dimensional integral is And the grid discretization function is The dimensionless pressure within the oil film mesh cell The finite element approximation, For the e-th oil film unit The x-coordinate of the reference element of the node. For the e-th oil film unit The reference element ordinate of the node, For the e-th oil film unit The weight value of the node. Let e ​​be the radial oil film reaction force of the e-th oil film unit. Let e ​​be the oil film reaction force of the e-th oil film unit in the tangential direction. The two-dimensional integral is And the grid discretization function is The trial function within the oil film mesh element, for Total number of oil film mesh cells in the direction, for Total number of oil film mesh cells in the direction, The radial reaction force across the entire oil film region. The tangential reaction force across the entire oil film region. The oil film reaction force in the horizontal direction in the Cartesian coordinate system. Let be the oil film reaction force in the vertical direction in the Cartesian coordinate system.

[0022] According to some embodiments of this application, the preset convergence condition is set by the following formula:

[0023] in, To preset the nonlinear convergence tolerance value, This is the preset linear convergence tolerance value.

[0024] A second aspect of this application provides a system for calculating the oil film force of an extrusion oil film damper, the system comprising: The dimensionless module is used to make the oil film pressure control equation dimensionless when constructing the oil film pressure control equation, so as to obtain the dimensionless oil film pressure control equation. The mesh generation module is used to divide the oil film region of the extrusion oil film damper into multiple oil film mesh units after obtaining the oil film region. The weak form conversion module is used to perform a weak form conversion on the dimensionless oil film pressure control equation based on the oil film region using the finite element method, so as to obtain a weak form oil film pressure control equation. The element stiffness matrix and element load vector determination module is used to determine the element stiffness matrix and element load vector corresponding to each oil film mesh element based on the weak form of the oil film pressure control equation. A global stiffness matrix and global load vector determination module is used to determine the global stiffness matrix based on the element-global mapping matrix and all the element stiffness matrices, and to determine the global load vector based on the element-global mapping matrix and all the element load vectors, given the construction of the element-global mapping matrix. The oil film force determination module is used to determine the oil film force of the extrusion oil film damper based on the global stiffness matrix and the global load vector under the condition of applying preset boundary conditions.

[0025] This system achieves dimensionless oil film pressure control equations by constructing them first, then subdividing them into multiple oil film mesh elements. Based on these oil film regions, the system performs a weak-form transformation of the dimensionless oil film pressure control equations using the finite element method, resulting in a weak-form oil film pressure control equation. Based on this weak-form equation, the system determines the element stiffness matrix and element load vector for each oil film mesh element. With the element-global mapping matrix constructed, the system determines the global stiffness matrix based on the element-global mapping matrix and all element stiffness matrices. Similarly, the system determines the global load vector based on the element-global mapping matrix and all element load vectors. Finally, with preset boundary conditions applied, the system determines the oil film force of the extrusion oil film damper based on the global stiffness matrix and global load vector. This process is independent of boundary conditions and improves the computational efficiency of the oil film force.

[0026] A third aspect of this application provides an electronic device for calculating the oil film force of an extrusion oil film damper, including at least one control processor and a memory for communicatively connecting to the at least one control processor; the memory stores instructions that can be executed by the at least one control processor, and the instructions are executed by the at least one control processor to enable the at least one control processor to perform the above-described method for calculating the oil film force of the extrusion oil film damper.

[0027] In a fourth aspect, this application provides a computer-readable storage medium storing computer-executable instructions for causing a computer to execute the above-described method for calculating the oil film force of a squeeze oil film damper.

[0028] It should be noted that the beneficial effects of the second to fourth aspects of this application compared with the prior art are the same as the beneficial effects of the oil film force calculation system for the extrusion oil film damper compared with the prior art, and will not be described in detail here.

[0029] Additional aspects and advantages of this application will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of this application. Attached Figure Description

[0030] The above and / or additional aspects and advantages of this application will become apparent and readily understood from the description of the embodiments taken in conjunction with the following drawings, in which: Figure 1 This is a schematic flowchart of an embodiment of the oil film force calculation method for the extrusion oil film damper provided in this application; Figure 2 This is a schematic diagram of the oil film coordinate system and symbols for the oil film force calculation method of the extrusion oil film damper provided in this application; Figure 3 This is a schematic diagram of journal whirl in the oil film force calculation method of the extrusion oil film damper provided in this application; Figure 4 This is a schematic diagram of the discretized oil film region using bilinear quadrilateral units for the oil film force calculation method of the extrusion oil film damper provided in this application; Figure 5 This is a schematic diagram of the bilinear isoparametric mapping of the oil film force calculation method for the extrusion oil film damper provided in this application; Figure 6 This is a comparative schematic diagram of the SBA (short bearing approximate analytical solution) and the FEM (finite element method) of the oil film force calculation method of the extrusion oil film damper provided in this application; Figure 7 This is a schematic diagram of the oil film force distribution along the circumferential direction using the LBA (long bearing approximate analytical solution) of the oil film force calculation method for the extrusion oil film damper provided in this application and the FEM (finite element method) of this application. Figure 8 This is a comparative schematic diagram of the SOR (super-relaxation iterative method) and the FEM (finite element method) of the oil film force calculation method of the extrusion oil film damper provided in this application; Figure 9 This is a comparative schematic diagram of the FDM (Finite Difference Method) and the FEM (Finite Element Method) of the present application for calculating the oil film force of the extrusion oil film damper. Figure 10 This is a schematic diagram of an embodiment of the oil film force calculation system for the extrusion oil film damper provided in this application; Figure 11 This is a schematic diagram of the structure of an embodiment of the electronic device provided in this application. Detailed Implementation

[0031] The embodiments of this application are described in detail below. Examples of the embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the accompanying drawings are exemplary and are only used to explain this application, and should not be construed as limiting this application.

[0032] In the description of this application, the use of terms such as "first," "second," etc., is for the purpose of distinguishing technical features only and should not be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated or the order of the technical features indicated.

[0033] In the description of this application, it should be understood that the orientation descriptions, such as up, down, etc., are based on the orientation or positional relationship shown in the accompanying drawings, and are only for the convenience of describing this application and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation, and therefore should not be construed as a limitation of this application.

[0034] In the description of this application, it should be noted that, unless otherwise explicitly defined, terms such as "setup," "installation," and "connection" should be interpreted broadly, and those skilled in the art can reasonably determine the specific meaning of the above terms in this application in conjunction with the specific content of the technical solution.

[0035] Rotor-support systems are widely used in modern industrial production. As rotors develop towards higher speeds, greater compactness (smaller structures), and higher parameters (high speed, high pressure, large flow), the integrity and reliability of their structures become increasingly important. Vibration and stability of rotor-support systems are often the root causes of these problems.

[0036] Squeeze film dampers (SFDs) suppress vibration by adding an external damping term to the shaft support system. For elastically supported rotor systems, a well-designed SFD can effectively reduce the external excitation of the rotor, significantly lower the vibration peak value when the rotor is overcritical, and thus reduce the impact on the support structure, making it an effective and economical technical means to solve the vibration and noise problems of high-speed rotating machinery. Therefore, the calculation of the SFD film force is particularly important for the design of the SFD-rotor system and the analysis of the system's dynamic response.

[0037] Currently, the SFD oil film force is mainly calculated using the over-relaxation iterative method. However, the over-relaxation iterative method can only handle working conditions with regular boundary conditions, and it usually requires hundreds of iterations to meet the convergence tolerance requirements, resulting in low computational efficiency.

[0038] To address the aforementioned technical deficiencies, this application provides a method and system for calculating the oil film force of an extrusion oil film damper.

[0039] Please see Figure 1 This is a flowchart illustrating a method for calculating the oil film force of a squeeze oil film damper, provided in an embodiment of this application. This method is applied to electronic devices, such as servers. Figure 1 As shown, the calculation method for the oil film force of this extrusion oil film damper includes: Step S101: Under the condition of constructing the oil film pressure control equation, the oil film pressure control equation is dimensionless to obtain the dimensionless oil film pressure control equation. Step S102: After obtaining the oil film region of the extrusion oil film damper, the oil film region is divided into meshes to obtain multiple oil film mesh units. Step S103: Based on the oil film region, the dimensionless oil film pressure control equation is transformed into a weak form using the finite element method to obtain the weak form oil film pressure control equation. Step S104: Based on the weak form of the oil film pressure control equation, determine the element stiffness matrix and element load vector corresponding to each oil film mesh element; Step S105: With the element-global mapping matrix constructed, determine the global stiffness matrix based on the element-global mapping matrix and all element stiffness matrices; determine the global load vector based on the element-global mapping matrix and all element load vectors. Step S106: Under the condition of applying preset boundary conditions, determine the oil film force of the extrusion oil film damper based on the global stiffness matrix and the global load vector.

[0040] The oil film force of the aforementioned extrusion oil film damper includes the oil film reaction force in the horizontal direction under the Cartesian coordinate system and the oil film reaction force in the vertical direction under the Cartesian coordinate system.

[0041] This application achieves dimensionless oil film pressure control equations by constructing them, obtaining dimensionless oil film pressure control equations; it then meshes the oil film region of the extrusion oil film damper to obtain multiple oil film mesh elements; based on the oil film region, it performs a weak-form transformation of the dimensionless oil film pressure control equations using the finite element method to obtain weak-form oil film pressure control equations; based on the weak-form oil film pressure control equations, it determines the element stiffness matrix and element load vector corresponding to each oil film mesh element; it determines the global stiffness matrix based on the element-global mapping matrix and all element stiffness matrices, and the global load vector based on the element-global mapping matrix and all element load vectors; and it determines the oil film force of the extrusion oil film damper based on the global stiffness matrix and global load vector, under the condition of applying preset boundary conditions. This achieves independence from the influence of boundary conditions and improves the calculation efficiency of oil film force.

[0042] In some embodiments, the oil film pressure control equation is dimensionless using the following formula to obtain the dimensionless oil film pressure control equation:

[0043] in, The oil film thickness at any dimensionless location. The pressure distribution value is the dimensionless value. This represents the horizontal displacement of the journal. This represents the vertical displacement of the journal. The oil film thickness at any location. Let be the radius of the outer wall surface of the journal. Let be the radius of the inner wall of the bearing housing. This represents the initial radial clearance of the oil film. This represents the pressure distribution value. For dimensionless parameters, For intermediate variable values, The angle between the journal and the horizontal axis in the counterclockwise direction. The x-axis velocity component is the dimensionless value. The velocity component in the y-direction after dimensionless scaling is denoted as . Here, θ is the rheological index, θ is the angle between the counterclockwise direction and the reference line, Z is the dimensionless axial coordinate, and L is the oil film length. This represents the dimensionless displacement of the journal in the horizontal direction. This represents the dimensionless displacement of the journal in the vertical direction. is the consistency coefficient of the fluid.

[0044] Specifically, refer to Figure 2 , Figure 2 In the diagram, "Upper wall" represents the upper surface, "Lower wall" represents the lower surface, and "Oil-film area" represents the oil film region. This oil film region is formed between two relatively moving rigid surfaces, filled with a non-Newtonian fluid. A Cartesian coordinate system is established with the X-axis representing the horizontal direction of the oil film, the Y-axis representing the oil film thickness direction, and the Z-axis representing the main flow direction (axial direction) of the oil film. Indicates oil film thickness, outer wall surface of journal Inner wall of bearing housing , Points on the outer wall surface The velocity component along the horizontal axial direction, Points on the outer wall surface The velocity component along the vertical direction, Points on the outer wall surface The velocity component along the axial direction, Represents a point on the inner wall surface The velocity component along the horizontal direction, Represents a point on the inner wall surface The velocity component along the vertical direction, Represents a point on the inner wall surface The velocity component along the axial direction.

[0045] The journal can only translate and not rotate in the XY plane. The bearing housing is fixed and exists. .

[0046] Before deriving the oil film pressure control equation, the following assumptions need to be made: (1) the oil film thickness is much smaller than its length and width; (2) the fluid is incompressible and its density is low. (3) Neglecting volume forces; (4) No slippage on the wall, and the fluid velocity at the boundary is the same as that on the journal and bearing housing surface; (5) Neglecting inertial forces; (6) Constant pressure along the film thickness direction. (7) Laminar flow; (8) The velocity gradient in the y-direction is much greater than the velocity gradients in the x and z directions; The constitutive equation of the power-law model can be expressed as:

[0047] in, For shear stress, The consistency coefficient of the fluid. Shear strain rate It is a dimensionless rheological index.

[0048] Based on the above assumptions (3), (5) and (6), the simplified momentum equation is obtained:

[0049] Where P represents the pressure distribution. Let be the shear stress component along the x-direction in the Y-plane. Let be the shear stress component along the z-direction in the Y-plane.

[0050] make The volumetric flow rate per unit width in the z-direction is defined as: ,in, This represents the volumetric flow rate per unit width in the z-direction.

[0051] Define the volumetric flow rate per unit width in the x-direction as: ,in, This represents the volumetric flow rate per unit width in the x-direction.

[0052] The continuity equation for incompressible fluids can be written as: Along the oil film thickness direction from 0 to Integrating and simplifying, we can obtain After substituting into the flow equation, the two-dimensional Reynolds equation using the finite difference method for modeling non-Newtonian fluids with a power-law model can be simplified to:

[0053] Reference Figure 3 , Center of journal The center of the bearing housing, and The line connecting them serves as the reference baseline. This refers to the eccentricity of the journal during the vortex process. Indicates the journal eccentricity and the oil film thickness at any position. .

[0054] The velocity components of the journal surface in the horizontal and vertical directions can be expressed as follows:

[0055] according to The Reynolds equation in cylindrical coordinates can be obtained as follows:

[0056] Set dimensionless parameters , , , , , The dimensionless oil film pressure control equation is: .

[0057] This application reduces calculation errors and improves calculation efficiency by making the oil film pressure control equation dimensionless.

[0058] In some embodiments, when obtaining the oil film region of the extrusion oil film damper, the oil film region is meshed to obtain multiple oil film mesh elements, which can be: Reference Figure 4 The oil film area is expanded into a rectangular area, in direction and Draw the directions separately and Equally spaced lines are used to construct a discrete space using Q4 bilinear quadrilateral elements. Each element has 4 nodes, and the total number of grid cells is [number missing]. The size of a single grid is The total number of nodes is .

[0059] In some embodiments, the dimensionless oil film pressure control equation is transformed into a weak form using the finite element method based on the oil film region, according to the following formula:

[0060] in, This is the value of the first intermediate parameter. This is the value of the second intermediate parameter. This is the value of the third intermediate parameter. This is the value of the fourth intermediate parameter. For the oil film area, For the trial function, Let be the set of trial functions.

[0061] Specifically, in the expanded rectangular area The pressure satisfies the following formula: ; Multiply both sides of the formula that the pressure satisfies by any... And in Integrating further, we transform the equation into a weak form, yielding the weak form of the oil film pressure control equation: .

[0062] The finite element method is used to perform a weak form transformation on the dimensionless oil film pressure control equation, which improves the accuracy of pressure calculation in the boundary region.

[0063] In some embodiments, based on the weak form of the oil film pressure control equation, the element stiffness matrix and element load vector corresponding to each oil film mesh element are determined, including: When the trial function is a 4-node bilinear function, through 2 2. The Gaussian integration method constructs the partial derivative at the Gaussian point corresponding to each oil film grid cell using the following formula:

[0064] in, Let be the shape function of the first node among the four nodes corresponding to the oil film mesh element. Let be the shape function of the second node out of the four nodes corresponding to the oil film mesh element. Let be the shape function of the third node out of the four nodes corresponding to the oil film mesh element. Let be the shape function of the fourth node out of the four nodes corresponding to the oil film mesh element. The x-axis of the reference unit is... The ordinate of the reference unit; Based on the partial derivatives at the Gaussian point corresponding to each oil film mesh element and the weak form of the oil film pressure control equation, the element stiffness matrix corresponding to each oil film mesh element is determined by the Picard iteration method and the following formula:

[0065] in, The two-dimensional integral is And the grid discretization function is The true value of stress, The two-dimensional integral is And the grid discretization function is The trial function within the oil film mesh element, The fourth node among the four nodes corresponding to the oil film mesh element The stress value of the node, The fourth node among the four nodes corresponding to the oil film mesh element The node's trial function, The two-dimensional integral is The The shape function of the node. The two-dimensional integral is The The value of the fifth intermediate parameter in the next iteration. For the first Pressure distribution value in the next iteration This is the set of four nodes corresponding to the oil film mesh element. For Jacobian matrices, For the first Oil film mesh unit Dimensional value in direction, For the first Oil film mesh unit Dimensional value in direction, For the first The first oil film mesh unit The element stiffness matrix of the next iteration The two-dimensional integral is The The value of the fifth intermediate parameter in the next iteration. For the oil film area corresponding to the oil film grid cell, For transpose; The element load vector corresponding to each oil film mesh element is determined by the following formula:

[0066] in, The fourth node among the four nodes corresponding to the oil film mesh element The weight value of the node. The fourth node among the four nodes corresponding to the oil film mesh element The weight value of the node. For the first The first oil film mesh unit The element load vector of the next iteration. The two-dimensional integral is The The fourth intermediate parameter of the node.

[0067] Specifically, refer to Figure 5 Considering the computational domain The circumferential-axial rectangular feature, using a reference square To the physical quadrilateral Bilinear isoparametric mapping is used to handle geometry and gradients, and 2 is used within the oil film mesh element. 2. By calculating each term using Gaussian integrals, a linear algebraic system that can be directly solved is assembled. Figure 5 In Let G1 be the normal vector, G2 be the Gaussian integration points, and N1 be the shape functions corresponding to the four nodes.

[0068] Reference unit Mapping to The first coordinate alignment Oil film mesh unit At this time, the mapping can be written as a one-dimensional linear relationship that is independent of each other in two directions:

[0069] in, Let be the geometric center coordinates of the oil film mesh element.

[0070] Calculate the corresponding Jacobian matrix.

[0071] When the trial function is a 4-node bilinear function, through 2 2. Gaussian integration method to construct the partial derivative at the Gaussian point corresponding to each oil film grid cell.

[0072] Based on the partial derivatives at the Gaussian point corresponding to each oil film mesh element and the weak form of the oil film pressure control equation, the element stiffness matrix corresponding to each oil film mesh element is determined by the Picard iteration method.

[0073] remember , According to the chain rule, it exists , , , The product in the weak form can be rewritten as a vector quadratic form: for Xiang You ,right Xiang You ,in, Define parameter one for the intermediate stage. Define parameter two for the intermediate stage. For the first Pressure value of each oil film mesh cell For the first Trial function for each oil film mesh element.

[0074] Given the first The result of the second iteration In the The next iteration Evaluate This transforms the problem into a forced linear problem and updates it. Continue until the residuals satisfy the convergence criterion.

[0075] No. The expression for the weak form of the oil film pressure control equation in the next iteration is as follows:

[0076] For the first iteration, if the initial value If we take 0, then This can lead to ill-conditioned linear systems. To ensure the stability of the numerical solution, a micro-regularization strategy is employed during the computation process, namely... ,in, , The value can be .

[0077] According to the The left-hand side of the weak form of the oil film pressure control equation in the next iteration is defined as follows: The element stiffness matrix corresponding to each oil film mesh element in each iteration.

[0078] According to the The right-hand side of the weak form of the oil film pressure control equation in the next iteration defines the... The element load vector corresponding to each oil film mesh element in each iteration.

[0079] After obtaining the element stiffness matrix and element load vector Then, using the unit-global mapping matrix Will and Embed the global stiffness matrix respectively and global load vector In the middle, where DOF equals When numbering global nodes, a row-major principle is adopted, that is, numbering from left to right and from bottom to top. Here, DOF represents the global degrees of freedom. It represents the set of real numbers.

[0080]

[0081] Where (s, r) represent the abscissa and ordinate values ​​of the straight lines along the θ and Z directions when meshing the oil film region, respectively. The e-th element is numbered sequentially using Q4 local nodes starting from the lower left corner; the element numbering is as follows:

[0082] Mapping matrix A sparse matrix containing only elements 0 and 1 is specifically defined as:

[0083] in, During the indexing process The OK, During the indexing process The List.

[0084] Based on the weak form of the oil film pressure control equation, this application determines the element stiffness matrix and element load vector corresponding to each oil film grid element, providing more accurate data for subsequent calculation of the oil film force of the extrusion oil film damper, thereby improving the calculation efficiency and accuracy of the oil film force.

[0085] In some embodiments, the global stiffness matrix is ​​determined using the following formula:

[0086] in, For the first The global stiffness matrix of the next iteration. For the first The cell-global mapping matrix of each oil film mesh cell; The global load vector is determined using the following formula:

[0087] in, For the first The global load vector for the next iteration.

[0088] In some embodiments, the oil film force of the extrusion oil film damper is determined based on the global stiffness matrix and the global load vector, including: Based on the global stiffness matrix and global load vector, the first... Pressure distribution value in the next iteration:

[0089] in, For the first Pressure distribution values ​​in the next iteration; In the If the pressure distribution value of the next iteration meets the preset convergence condition, the oil film force of the extrusion oil film damper is calculated using the following formula:

[0090] in, This represents the true pressure value. For the oil film region of the e-th oil film unit, The two-dimensional integral is The true value of stress, The two-dimensional integral is And the grid discretization function is The dimensionless pressure within the oil film mesh cell The finite element approximation, For the e-th oil film unit The x-coordinate of the reference element of the node. For the e-th oil film unit The reference element ordinate of the node, For the e-th oil film unit The weight value of the node. Let e ​​be the radial oil film reaction force of the e-th oil film unit. Let e ​​be the oil film reaction force of the e-th oil film unit in the tangential direction. The two-dimensional integral is And the grid discretization function is The trial function within the oil film mesh element, for Total number of oil film mesh cells in the direction, for Total number of oil film mesh cells in the direction, The radial reaction force across the entire oil film region. The tangential reaction force across the entire oil film region. The oil film reaction force in the horizontal direction in the Cartesian coordinate system. Let be the oil film reaction force in the vertical direction in the Cartesian coordinate system.

[0091] This application determines the oil film force of the extrusion oil film damper based on the global stiffness matrix and global load vector, thereby improving the calculation accuracy of the oil film force.

[0092] In some embodiments, a preset convergence condition is set using the following formula:

[0093] in, To preset the nonlinear convergence tolerance value, This is the preset linear convergence tolerance value.

[0094] Specifically, Denote the Euclidean norm, take , .

[0095] This application improves the convergence efficiency by setting preset convergence conditions, thereby improving the calculation efficiency of oil film force.

[0096] Specifically, refer to Figure 6 ,when , hour, For an oil film length L = 10 mm, journal radius R = 50 mm, oil film radial clearance c = 0.2 mm, horizontal displacement x = 0.1 mm, vertical displacement y = 0 mm, and horizontal velocity... m / s, vertical velocity For operating conditions with a flow rate of m / s and a dynamic viscosity of µ = 0.02 Pa·s, the length-to-diameter ratio L / R = 0.2 < 1 / 4 satisfies the approximate solution for short bearings. The oil film pressure of the SFD is calculated using both the approximate analytical solution for short bearings and the proposed finite element method. Figure 6 As shown, the oil film pressure calculated by the proposed finite element method is basically consistent with the oil film pressure calculated by the approximate solution of the short bearing, whether in the axial or circumferential direction.

[0097] Specifically, refer to Figure 7For an oil film length L = 400 mm, journal radius R = 50 mm, oil film radial clearance c = 0.2 mm, horizontal displacement x = 0.1 mm, vertical displacement y = 0 mm, and horizontal velocity... m / s, vertical velocity For operating conditions with m / s and dynamic viscosity µ = 0.02 pa·s, the length-to-diameter ratio L / R = 4 satisfies the applicable range of the approximate solution for long bearings. For example... Figure 7 As shown, the calculation results of the two methods agree well in the circumferential direction. Under the condition of small length-to-diameter ratio, the SFD oil film distribution pressure calculated by the proposed finite element method is consistent with the result calculated by the analytical solution of the short bearing, and under the condition of large length-to-diameter ratio, it is consistent with the result calculated by the approximate analytical solution of the long bearing. This verifies the effectiveness of the proposed finite element method.

[0098] Specifically, refer to Figure 8 For an oil film length L = 20 mm, journal radius R = 50 mm, oil film radial clearance c = 0.2 mm, horizontal displacement x = 0.1 mm, vertical displacement y = 0 mm, and horizontal velocity... m / s, vertical velocity Under operating conditions of m / s and dynamic viscosity µ = 0.02 Pa·s, such as Figure 8 As shown, the oil film distribution pressure calculated by the over-relaxation iterative method and the finite element method remains consistent in both the axial and circumferential directions, but the solution efficiency differs significantly. This computer has an i5-8300H CPU and 8 GB of RAM, and is programmed using Julia with a mesh density set to (40...). 20) For SOR, with a relaxation factor of 1, the equation converged after 70 iterations in 0.323 seconds; for the proposed FEM, convergence was achieved in just 2 iterations in only 0.048 seconds. The solution efficiency of the finite element method is approximately 6.7 times that of the over-relaxation iterative method. Furthermore, the over-relaxation iterative method, based on the finite difference method, can only handle cases with regular boundary conditions, while the finite element method can adapt to arbitrary boundary conditions, thus having a wider range of applications.

[0099] Specifically, refer to Figure 9 For an oil film length L = 20 mm, journal radius R = 50 mm, oil film radial clearance c = 0.2 mm, horizontal displacement x = 0.1 mm, vertical displacement y = 0 mm, and horizontal velocity... m / s, vertical velocity Under operating conditions of m / s and dynamic viscosity µ = 0.02 Pa·s, such as Figure 9As shown, the oil film distribution pressure calculated using both the finite difference method and the finite element method remains consistent in both the axial and circumferential directions. However, the finite difference method can only handle conditions with regular boundary conditions, while the finite element method can adapt to arbitrary boundary conditions and has a wider range of applications. The finite difference method takes 0.207 seconds, which is 4.31 times that of this invention.

[0100] Additionally, refer to Figure 10 One embodiment of this application provides a system for calculating the oil film force of a squeeze oil film damper, including a data dimensionless conversion module 1100, a mesh generation module 1200, a weak form transformation module 1300, an element stiffness matrix and element load vector determination module 1400, a global stiffness matrix and global load vector determination module 1500, and an oil film force determination module 1600, wherein: The dimensionless transformation module 1100 is used to transform the oil film pressure control equation into a dimensionless form when constructing the oil film pressure control equation, so as to obtain the dimensionless oil film pressure control equation. The mesh generation module 1200 is used to divide the oil film region into multiple oil film mesh units after obtaining the oil film region of the extrusion oil film damper. The weak form conversion module 1300 is used to perform a weak form conversion on the dimensionless oil film pressure control equation based on the oil film region using the finite element method, so as to obtain the weak form oil film pressure control equation. The element stiffness matrix and element load vector determination module 1400 is used to determine the element stiffness matrix and element load vector corresponding to each oil film mesh element based on the weak form of the oil film pressure control equation. The global stiffness matrix and global load vector determination module 1500 is used to determine the global stiffness matrix based on the element-global mapping matrix and all element stiffness matrices, and to determine the global load vector based on the element-global mapping matrix and all element load vectors, given the element-global mapping matrix has been constructed. The oil film force determination module 1600 is used to determine the oil film force of the extrusion oil film damper based on the global stiffness matrix and the global load vector under the condition of applying preset boundary conditions.

[0101] This system achieves dimensionless oil film pressure control equations by constructing them first, then subdividing them into multiple oil film mesh elements. Based on these oil film regions, the system performs a weak-form transformation of the dimensionless oil film pressure control equations using the finite element method, resulting in a weak-form oil film pressure control equation. Based on this weak-form equation, the system determines the element stiffness matrix and element load vector for each oil film mesh element. With the element-global mapping matrix constructed, the system determines the global stiffness matrix based on the element-global mapping matrix and all element stiffness matrices. Similarly, the system determines the global load vector based on the element-global mapping matrix and all element load vectors. Finally, with preset boundary conditions applied, the system determines the oil film force of the extrusion oil film damper based on the global stiffness matrix and global load vector. This process is independent of boundary conditions and improves the computational efficiency of the oil film force.

[0102] It should be noted that the system embodiments described above are based on the same inventive concept as the method embodiments described above. Therefore, the relevant content of the method embodiments described above is also applicable to the system embodiments described above, and will not be repeated here.

[0103] Figure 11 A schematic diagram of the hardware structure for calculating the oil film force of the extrusion oil film damper provided in an embodiment of this application is shown.

[0104] The oil film force calculation device in the extrusion oil film damper may include a processor 301 and a memory 302 storing computer program instructions.

[0105] Specifically, the processor 301 may include a central processing unit (CPU), an application-specific integrated circuit (ASIC), or one or more integrated circuits that can be configured to implement the embodiments of this application.

[0106] Memory 302 may include mass storage for data or instructions. For example, and not limitingly, memory 302 may include a hard disk drive (HDD), floppy disk drive, flash memory, optical disk, magneto-optical disk, magnetic tape, or Universal Serial Bus (USB) drive, or a combination of two or more of these. Where appropriate, memory 302 may include removable or non-removable (or fixed) media. Where appropriate, memory 302 may be internal or external to the integrated gateway disaster recovery device. In a particular embodiment, memory 302 is non-volatile solid-state memory.

[0107] In some embodiments, memory 302 may include read-only memory (ROM), random access memory (RAM), disk storage media device, optical storage media device, flash memory device, electrical, optical, or other physical / tangible memory storage device. Thus, generally, memory includes one or more tangible (non-transitory) computer-readable storage media (e.g., memory devices) encoded with software including computer-executable instructions, and when the software is executed (e.g., by one or more processors), it is operable to perform the operations described with reference to the method according to one aspect of this disclosure.

[0108] The processor 301 reads and executes the computer program instructions stored in the memory 302 to implement any of the oil film force calculation methods for the extrusion oil film damper in the above embodiments.

[0109] In one example, the oil film force calculation device for the squeeze oil film damper may further include a communication interface 303 and a bus 310. Wherein, as Figure 11 As shown, the processor 301, memory 302, and communication interface 303 are connected through bus 310 and complete communication with each other.

[0110] The communication interface 303 is mainly used to realize communication between various modules, devices, units and / or equipment in the embodiments of this application.

[0111] Bus 310 includes hardware, software, or both, that couples components of the oil film force calculation device for the squeeze film damper together. For example, and not limitingly, the bus may include an Accelerated Graphics Port (AGP) or other graphics bus, an Enhanced Industry Standard Architecture (EISA) bus, a Front Side Bus (FSB), HyperTransport (HT) interconnect, an Industry Standard Architecture (ISA) bus, an Unlimited Available Network Bandwidth Interconnect, a Low Pin Count (LPC) bus, a memory bus, a Microchannel Architecture (MCA) bus, a Peripheral Component Interconnect (PCI) bus, a PCI-Express (PCI-X) bus, a Serial Advanced Technology Attachment (SATA) bus, a Video Electronics Standards Association Local (VLB) bus, or other suitable buses, or combinations of two or more of these. Where appropriate, bus 310 may include one or more buses. Although specific buses are described and illustrated in embodiments of this application, any suitable bus or interconnect is contemplated herein.

[0112] The oil film force calculation device for the extrusion oil film damper can execute the oil film force calculation method of the extrusion oil film damper in this application embodiment based on a three-dimensional design model, thereby achieving a combination of Figure 1 and Figure 10 The method and system for calculating the oil film force of the extrusion oil film damper are described.

[0113] Furthermore, in conjunction with the oil film force calculation method for the extrusion oil film damper in the above embodiments, this application embodiment can provide a computer storage medium for implementation. This computer storage medium stores computer program instructions; when these computer program instructions are executed by a processor, they implement any of the oil film force calculation methods for the extrusion oil film damper in the above embodiments.

[0114] It should be clarified that this application is not limited to the specific configurations and processes described above and shown in the figures. For the sake of brevity, detailed descriptions of known methods are omitted here. In the above embodiments, several specific steps are described and shown as examples. However, the method process of this application is not limited to the specific steps described and shown. Those skilled in the art can make various changes, modifications, and additions, or change the order of steps, after understanding the spirit of this application.

[0115] The functional blocks shown in the above block diagram can be implemented as hardware, software, firmware, or a combination thereof. When implemented in hardware, they can be, for example, electronic circuits, application-specific integrated circuits (ASICs), appropriate firmware, plug-ins, function cards, etc. When implemented in software, the elements of this application are programs or code segments used to perform the required tasks. Programs or code segments can be stored on a machine-readable medium or transmitted over a transmission medium or communication link via data signals carried on a carrier wave. "Machine-readable medium" can include any medium capable of storing or transmitting information. Examples of machine-readable media include electronic circuits, semiconductor memory devices, ROM, flash memory, erasable ROM (EROM), floppy disks, CD-ROMs, optical disks, hard disks, fiber optic media, radio frequency (RF) links, etc. Code segments can be downloaded via computer networks such as the Internet, intranets, etc.

[0116] It should also be noted that the exemplary embodiments mentioned in this application describe methods or systems based on a series of steps or apparatus. However, this application is not limited to the order of the above steps; that is, the steps can be performed in the order mentioned in the embodiments, or in a different order, or several steps can be performed simultaneously.

[0117] The aspects of this disclosure have been described above with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this disclosure. It should be understood that each block in the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing apparatus to produce a machine such that these instructions, executable via the processor of the computer or other programmable data processing apparatus, enable the implementation of the functions / actions specified in one or more blocks of the flowchart illustrations and / or block diagrams. Such a processor can be, but is not limited to, a general-purpose processor, a special-purpose processor, a special application processor, or a field-programmable logic circuit. It is also understood that each block in the block diagrams and / or flowcharts, and combinations of blocks in the block diagrams and / or flowcharts, can also be implemented by special-purpose hardware performing the specified functions or actions, or can be implemented by a combination of special-purpose hardware and computer instructions.

[0118] The above description is merely a specific implementation of this application. Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working processes of the systems, modules, and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here. It should be understood that the protection scope of this application is not limited thereto. Any person skilled in the art can easily conceive of various equivalent modifications or substitutions within the technical scope disclosed in this application, and these modifications or substitutions should all be covered within the protection scope of this application.

Claims

1. A method for calculating the oil film force of a squeeze oil film damper, characterized in that, The method for calculating the oil film force of the extrusion oil film damper includes: In constructing the oil film pressure control equation, the oil film pressure control equation is dimensionless to obtain the dimensionless oil film pressure control equation. Having obtained the oil film region of the extrusion oil film damper, the oil film region is divided into grids to obtain multiple oil film grid units; Based on the oil film region, the dimensionless oil film pressure control equation is transformed into a weak form using the finite element method to obtain the weak form of the oil film pressure control equation. Based on the weak form of the oil film pressure control equation, the element stiffness matrix and element load vector corresponding to each oil film mesh element are determined. In constructing the element-global mapping matrix, the global stiffness matrix is ​​determined based on the element-global mapping matrix and all the element stiffness matrices; the global load vector is determined based on the element-global mapping matrix and all the element load vectors. Under the condition of applying a preset boundary condition, the oil film force of the extrusion oil film damper is determined based on the global stiffness matrix and the global load vector.

2. The method for calculating the oil film force of a squeeze oil film damper according to claim 1, characterized in that, The oil film pressure control equation is dimensionlessly transformed using the following formula to obtain the dimensionless oil film pressure control equation: in, The oil film thickness at any dimensionless location. The pressure distribution value is the dimensionless value. This represents the horizontal displacement of the journal. This represents the vertical displacement of the journal. The oil film thickness at any location. Let be the radius of the outer wall surface of the journal. Let be the radius of the inner wall of the bearing housing. This represents the initial radial clearance of the oil film. This represents the pressure distribution value. For dimensionless parameters, For intermediate variable values, The angle between the journal and the horizontal axis in the counterclockwise direction. The x-axis velocity component is the dimensionless value. The velocity component in the y-direction after dimensionless scaling is denoted as . Here, θ is the rheological index, θ is the angle between the counterclockwise direction and the reference line, Z is the dimensionless axial coordinate, and L is the oil film length. This represents the dimensionless displacement of the journal in the horizontal direction. This represents the dimensionless displacement of the journal in the vertical direction. is the consistency coefficient of the fluid.

3. The method for calculating the oil film force of a squeeze oil film damper according to claim 2, characterized in that, Based on the oil film region, the dimensionless oil film pressure control equation is transformed into a weak form using the finite element method, yielding the following formula: in, This is the value of the first intermediate parameter. This is the value of the second intermediate parameter. This is the value of the third intermediate parameter. This is the value of the fourth intermediate parameter. For the oil film area, For the trial function, Let be the set of trial functions.

4. The method for calculating the oil film force of a squeeze oil film damper according to claim 3, characterized in that, The oil film pressure control equation based on the weak form determines the element stiffness matrix and element load vector corresponding to each oil film mesh element, including: When the trial function is a 4-node bilinear function, through 2 2. The Gaussian integration method constructs the partial derivative at the Gaussian point corresponding to each oil film grid cell using the following formula: in, Let be the shape function of the first node among the four nodes corresponding to the oil film mesh element. Let be the shape function of the second node out of the four nodes corresponding to the oil film mesh element. Let be the shape function of the third node out of the four nodes corresponding to the oil film mesh element. Let be the shape function of the fourth node out of the four nodes corresponding to the oil film mesh element. The x-axis of the reference unit is... The ordinate of the reference unit; Based on the partial derivatives at the Gaussian point corresponding to each oil film mesh element and the weak form of the oil film pressure control equation, the element stiffness matrix corresponding to each oil film mesh element is determined by the Picard iteration method and the following formula: in, The two-dimensional integral is And the grid discretization function is The dimensionless pressure within the oil film mesh cell The finite element approximation, The two-dimensional integral is And the grid discretization function is The trial function within the oil film mesh element, The fourth node among the four nodes corresponding to the oil film mesh element The stress value of the node, The fourth node among the four nodes corresponding to the oil film mesh element The node's trial function, The two-dimensional integral is The The shape function of the node. The two-dimensional integral is The The value of the fifth intermediate parameter in the next iteration. For the first Pressure distribution value in the next iteration This is the set of four nodes corresponding to the oil film mesh element. For Jacobian matrices, For the first Oil film mesh unit Dimensional value in direction, For the first Oil film mesh unit Dimensional value in direction, For the first The first oil film mesh unit The element stiffness matrix of the next iteration The two-dimensional integral is The The value of the fifth intermediate parameter in the next iteration. For the oil film area corresponding to the oil film grid cell, For transpose; The element load vector corresponding to each oil film mesh element is determined by the following formula: in, The fourth node among the four nodes corresponding to the oil film mesh element The weight value of the node. The fourth node among the four nodes corresponding to the oil film mesh element The weight value of the node. For the first The first oil film mesh unit The element load vector of the next iteration. The two-dimensional integral is The The fourth intermediate parameter of the node.

5. The method for calculating the oil film force of a squeeze oil film damper according to claim 4, characterized in that, The global stiffness matrix is ​​determined using the following formula: in, For the first The global stiffness matrix of the next iteration. For the first The cell-global mapping matrix of each oil film mesh cell; The global load vector is determined using the following formula: in, For the first The global load vector for the next iteration.

6. The method for calculating the oil film force of a squeeze oil film damper according to claim 4, characterized in that, The oil film force of the extrusion oil film damper includes the oil film reaction force in the horizontal direction under Cartesian coordinates and the oil film reaction force in the vertical direction under Cartesian coordinates. Determining the oil film force of the extrusion oil film damper based on the global stiffness matrix and the global load vector includes: Based on the global stiffness matrix and the global load vector, the first... Pressure distribution value in the next iteration: in, For the first Pressure distribution values ​​in the next iteration; In the first If the pressure distribution value of the next iteration meets the preset convergence condition, the oil film force of the extrusion oil film damper is calculated using the following formula: in, This represents the true pressure value. For the oil film region of the e-th oil film unit, The two-dimensional integral is The true value of stress, The two-dimensional integral is And the grid discretization function is The dimensionless pressure within the oil film mesh cell The finite element approximation, For the e-th oil film unit The x-coordinate of the reference element of the node. For the e-th oil film unit The reference element ordinate of the node, For the e-th oil film unit The weight value of the node. Let e ​​be the radial oil film reaction force of the e-th oil film unit. Let e ​​be the oil film reaction force of the e-th oil film unit in the tangential direction. The two-dimensional integral is And the grid discretization function is The trial function within the oil film mesh element, for Total number of oil film mesh cells in the direction, for Total number of oil film mesh cells in the direction, The radial reaction force across the entire oil film region. The tangential reaction force across the entire oil film region. The oil film reaction force in the horizontal direction in the Cartesian coordinate system. Let be the oil film reaction force in the vertical direction in the Cartesian coordinate system.

7. The method for calculating the oil film force of a squeeze oil film damper according to claim 6, characterized in that, The preset convergence condition is set using the following formula: in, To preset the nonlinear convergence tolerance value, This is the preset linear convergence tolerance value.

8. A system for calculating the oil film force of a squeeze oil film damper, characterized in that, The oil film force calculation system of the extrusion oil film damper includes: The dimensionless module is used to make the oil film pressure control equation dimensionless when constructing the oil film pressure control equation, so as to obtain the dimensionless oil film pressure control equation. The mesh generation module is used to divide the oil film region of the extrusion oil film damper into multiple oil film mesh units after obtaining the oil film region. The weak form conversion module is used to perform a weak form conversion on the dimensionless oil film pressure control equation based on the oil film region using the finite element method, so as to obtain a weak form oil film pressure control equation. The element stiffness matrix and element load vector determination module is used to determine the element stiffness matrix and element load vector corresponding to each oil film mesh element based on the weak form of the oil film pressure control equation. A global stiffness matrix and global load vector determination module is used to determine the global stiffness matrix based on the element-global mapping matrix and all the element stiffness matrices, and to determine the global load vector based on the element-global mapping matrix and all the element load vectors, given the construction of the element-global mapping matrix. The oil film force determination module is used to determine the oil film force of the extrusion oil film damper based on the global stiffness matrix and the global load vector under the condition of applying preset boundary conditions.

9. A device for calculating the oil film force of a squeeze oil film damper, characterized in that, It includes at least one control processor and a memory for communicatively connecting to the at least one control processor; the memory stores instructions executable by the at least one control processor, which, when executed by the at least one control processor, enable the at least one control processor to perform a method for calculating the oil film force of a squeeze oil film damper as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer-executable instructions for causing a computer to perform a method for calculating the oil film force of a squeeze oil film damper as described in any one of claims 1 to 7.

Images