CN110083924B - Oil film lubrication performance simulation method under static pressure thrust bearing unbalance loading working condition - Google Patents

Oil film lubrication performance simulation method under static pressure thrust bearing unbalance loading working condition Download PDF

Info

Publication number
CN110083924B
CN110083924B CN201910332243.8A CN201910332243A CN110083924B CN 110083924 B CN110083924 B CN 110083924B CN 201910332243 A CN201910332243 A CN 201910332243A CN 110083924 B CN110083924 B CN 110083924B
Authority
CN
China
Prior art keywords
oil film
oil
block
setting
thrust bearing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910332243.8A
Other languages
Chinese (zh)
Other versions
CN110083924A (en
Inventor
于晓东
王宇
伍广鹏
周文凯
毕宏伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin University of Science and Technology
Original Assignee
Harbin University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Priority to CN201910332243.8A priority Critical patent/CN110083924B/en
Publication of CN110083924A publication Critical patent/CN110083924A/en
Application granted granted Critical
Publication of CN110083924B publication Critical patent/CN110083924B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Sliding-Contact Bearings (AREA)

Abstract

The invention discloses a method for simulating oil film performance under a static pressure thrust bearing unbalance loading working condition, and mainly relates to a method for simulating oil film lubricating performance under a double rectangular cavity static pressure thrust bearing unbalance loading working condition. The invention aims to solve the problem that the static pressure thrust bearing does not have an actual effective method for simulating the lubricating property of an oil film under the condition of unbalance loading. The method comprises the following steps: 1. establishing 6 symmetrical double-rectangular oil pad models through UG; 2. leading the model into ANSYS ICEM CFD software, creating an oil film model part, and creating a topological structure division structure grid; 3. checking whether the points and lines of the model and the points and lines of the block are all associated; 4. dividing the number of grid nodes, and encrypting the number of oil inlet nodes; 5. setting boundary oil film conditions in ANSYS CFX software; 6. and adopting the default standard of ANSYS CFX software for convergence, and judging whether the fluctuation is stable. The method is applied to the field of lubrication performance simulation under the condition of unbalance loading of the hydrostatic thrust bearing.

Description

Oil film lubrication performance simulation method under static pressure thrust bearing unbalance loading working condition
Technical Field
The invention relates to a method for simulating lubricating performance of an oil film of a hydrostatic thrust bearing, in particular to a method for simulating lubricating performance of an oil film when a double rectangular cavity is in unbalanced load.
Background
The hydrostatic thrust bearing is one of the core parts of heavy numerical control machining equipment, and the better the performance of the hydrostatic thrust bearing is, the higher the machining quality and the operating efficiency of the equipment are. Under the ideal state, the center of mass of the workpiece is positioned at the center of the rotary worktable to bear the central load, and the working state of the hydrostatic thrust bearing is in the optimal performance state at the moment. However, in the actual machining process, the mass center of the workpiece may not be at the center of the rotary table, so that the eccentric load condition occurs, the oil film of the hydrostatic thrust bearing is locally thinned, and even dry friction occurs, so that the hydrostatic thrust bearing fails, and the problem restricts the development of heavy numerical control machining equipment in China to the direction of high precision and high efficiency. And the simulation method for the oil film lubrication performance of the hydrostatic thrust bearing under the condition of unbalance loading is less researched at home and abroad, so that the research on the new simulation method for the lubrication performance of the hydrostatic thrust bearing under the condition of unbalance loading is necessary, and a basis is provided for the further optimization design of the hydrostatic thrust bearing.
Disclosure of Invention
The invention aims to solve the problems that in the actual working process of a hydrostatic thrust bearing, the center of mass of a workpiece is not coincident with the rotation center of a rotary table, unbalance loading occurs, the thickness of an oil film is uneven, the lubricating performance of the hydrostatic thrust bearing is greatly reduced, even dry friction occurs, and an effective method is lacked for simulating the phenomenon. In order to research the lubricating performance of an oil film of a hydrostatic thrust bearing during unbalance loading, a simulation method of the lubricating performance of the oil film of the hydrostatic thrust bearing during unbalance loading is provided.
The simulation method of the lubricating performance of the hydrostatic thrust bearing during unbalance loading is realized by the following technical scheme.
(1) Establishing 6 symmetrical double-rectangular oil pad three-dimensional models through UG three-dimensional modeling software;
(2) the model is led into ANSYS ICEM CFD software, the required surface is split, an oil film model part is created, and a topological structure division structure grid is created;
(3) checking whether the points and lines of the model are all associated with the points and lines of the block, if all the lines of the block are green, all the lines are associated, otherwise, all the lines are not associated, and therefore, the association is carried out again according to the steps;
(4) dividing the number of grid nodes; the number of nodes of the oil inlet is encrypted, and the number of adjacent nodes is not more than two, so that the grid quality is ensured;
(5) setting boundary oil film conditions in ANSYS CFX software, wherein the boundary oil film conditions comprise setting of viscosity-temperature relationship and setting of domains;
(6) the convergence adopts the default standard of ANSYS CFX software, namely the residual value is lower than 10-4And judging whether the fluctuation is stable or not and whether the fluctuation has an obvious descending trend or not.
The invention can simulate the actual working condition of the hydrostatic thrust bearing when the workpiece is subjected to unbalance loading through a computer, and obtain oil film lubrication parameters of the hydrostatic thrust bearing under different unbalance loading degrees. The method can save a large amount of time and obtain oil film lubrication parameters which accord with experimental results. By researching lubricating oil film parameters, the influence rule of different unbalance loading degrees on the lubricating property of the oil film is disclosed. The method provides an effective reference basis for further optimization research of the hydrostatic thrust bearing, and lays a solid foundation for realizing high efficiency and high precision of heavy equipment using the hydrostatic thrust bearing.
Drawings
FIG. 1 is a flow chart of a method for simulating oil film lubrication performance under a static pressure thrust bearing unbalance loading condition.
FIG. 2 is an oil film model under the unbalanced load condition mentioned in step A.
Fig. 3 is a post-Part definition list referred to in step B1.
Fig. 4 is a topological model diagram of the oil film mesh proposed in step B2.
Fig. 5 is a diagram of oil film grid division referred to in step C.
Fig. 6 is a diagram illustrating the oil film mesh quality in step D.
Fig. 7 shows the oil film position settings in step E.
Fig. 8 is a graph of the iterative residuals referred to in step F.
Detailed Description
The method can be realized by the following technical scheme,
step A: and (4) modeling the oil pad under the unbalanced loading condition by utilizing UG software, and creating models with different shapes and sizes according to different conditions.
Step B1: the established oil film models under different working conditions are exported into a. x.t file by utilizing UG software, and the exported. x.t file is imported into ANSYS ICEM CFD software, wherein the unit is set to be mm; because the CFD software cannot automatically divide the surface established in UG into two surfaces, the first step after model import is to split the surface in UG into required surfaces; secondly, part of an oil film model is created, which is mainly used for leading IN the setting of boundary conditions after CFX, and the part is set with the definition surfaces of IN1, IN2, OUT1, OUT2, OUT3, OUT4, INTERFACE1, INTERFACE2, ROTATE and WALL.
Step B2: and a top-down partitioning method is adopted when the topological structure partitioning structure grid is created. Firstly clicking an option Blocking, selecting a Check Block option, selecting a sub-topology structure in split Block, firstly dividing the X-direction Block, then dividing the Y, Z-direction Block, and deleting unnecessary topology structures. When the point mapping is created, all vertexes of the three-dimensional model are sequentially mapped to the vertexes of the block, and the Auto association automatic association option in the block can be directly selected for line mapping, so that the mapping of oil film and oil cavity points and lines is completed.
Step B3: the oil inlet is cylindrical, an O-shaped grid is used for the grid, the Block of the O-shaped grid is divided independently, two points are added to the circularly symmetrical part, the Block in the X, Y direction is divided, then a Check Block option is selected to stretch out the lower rectangular Block, the rectangular Block is further obtained, and during mapping, a quadrangle is mapped to a circle. And finally clicking the Select face option in the Split Block, selecting the upper surface and the lower surface of the cylinder, and clicking to determine, thus finishing all Block division and point-line mapping.
Step C: and checking whether the points and lines of the model are all related to the points and lines of the block, hiding both points and edges under Geometry, and checking whether the lines of the block are all green. If yes, all the relations are carried out; if the place is not green, it indicates that there is no correlation, so the correlation is performed again according to the above steps.
Step D: and clicking a Mesh Quality module to check the oil film grid and the angle grid Quality. And when the grids are Output, an Output/Select Solver is adopted to Select an ANSYS CFX Solver, and the default Output file format is CFX 5.
Step E1: ANSYS CFX software is selected to carry out simulation analysis on the oil film. Setting the FLUID domain as FLUID, selecting 32# lubricating oil as a material, selecting a FLUID calculation model as total energy, setting the temperature to be 293K, and keeping other settings as default settings.
And E2, setting an oil Inlet Inlet. The inlet IN1 and IN2 were set at a flow rate of 0.035, giving a room temperature of 293K, i.e. 20 ℃.
Step E3: setting of an oil Outlet, four oil outlets are provided, namely OUT1, OUT2, OUT3 and OUT4, and the relative atmospheric pressure is 0 Pa.
Step E4: clicking Interface, selecting Fluid by Interface Type, selecting INTERFACE1 by Interface side1 option, selecting INTERFACE2 by Interface side2 option, selecting the Rtta Axis rotation Axis as Z Axis, selecting GGI grid splicing mode from Mesh connecion option, and setting the periodic axial and radial tolerance as 0.0001.
Step E5: and setting a Wall surface Wall, wherein the upper surface of the oil film is a rotary working surface and is set to be roll Wall, and the surface of the lower surface of the oil film, which is in contact with the oil pad, and the inner Wall surface of the oil cavity are both set to be Wall.
Step E6: setting of ROTATE: and selecting rolling Wall by the Option, Rotating to rotate axis to select Z axis global Z, setting the Rotating speed of the rotary worktable around the Z axis to be 100r/min, and finishing the setting of the oil film boundary condition.
Step F: the convergence criterion used is the default criterion of ANSYS CFX software, i.e. the residual value is already lower than 10-4And the fluctuation condition is stable without obvious downward trend.

Claims (1)

1. The simulation method for the lubricating performance of the oil film under the condition of the unbalance loading of the hydrostatic thrust bearing is characterized in that the simulation method under the condition of the unbalance loading of the hydrostatic thrust bearing is realized according to the following steps:
step A: modeling the oil pad under the unbalanced loading condition by utilizing UG software, and creating models with different shapes and sizes according to different conditions;
step B1: the established oil film models under different working conditions are exported into a. x.t file by utilizing UG software, and the exported. x.t file is imported into ANSYS ICEM CFD software, wherein the unit is set to be mm; because the CFD software cannot automatically divide the surface established in UG into two surfaces, the first step after model import is to split the surface in UG into required surfaces; secondly, creating part of an oil film model, wherein the part is set to introduce the setting of boundary conditions after CFX, and the part is set with definition surfaces of IN1, IN2, OUT1, OUT2, OUT3, OUT4, INTERFACE1, INTERFACE2, ROTATE and WALL;
step B2: the method comprises the steps that a top-down partitioning method is adopted when a topological structure partitioning structure grid is created, firstly, an option Blocking is clicked, a Check Block option is selected, a sub-topological structure in a split Block is selected, an X-direction Block is partitioned firstly, then an Y, Z-direction Block is partitioned, unnecessary topological structures are deleted, when point mapping is created, all vertexes of a three-dimensional model are sequentially mapped to the vertexes of the Block, line mapping can directly select an Auto association option in the Blocking, and mapping of oil films, oil cavity points and lines is completed;
step B3: the oil inlet is cylindrical, an O-shaped grid is used for the grid, the Block of the grid is divided independently, two points are added to a circularly symmetric part, the Block in the X, Y direction is divided, then a Check Block option is selected to stretch out the rectangular Block below the Block, the rectangular Block is further obtained, a quadrangle is mapped to a circle during mapping, finally a Select face option in a Split Block is clicked, the upper surface and the lower surface of the cylinder are selected, and clicking determination is carried out, so that all Block division and point line mapping are completed;
and C: checking whether the points and lines of the model are all associated with the points and lines of the block, hiding the points and edges under Geometry, checking whether the lines of the block are all green, and if so, all associating; if a certain place is not green, indicating that no association exists, and performing association again according to the steps;
step D: clicking a Mesh Quality module, checking the Quality of an oil film grid and an angle grid, selecting an ANSYS CFX Solver by adopting an Output/Select Solver when the grid is Output, and defaulting the Output file to be CFX 5;
step E1: selecting ANSYS CFX software to perform analog simulation analysis on the oil film, setting a FLUID domain as FLUID, selecting a material as 32# lubricating oil, selecting a FLUID calculation model as total energy, setting the temperature to be 293K, and keeping other settings to be default settings;
step E2, setting the flow rate of inlets IN1 and IN2 to be 0.035 and setting the temperature to be room temperature 293K, namely 20 ℃ for setting an oil Inlet Inlet;
step E3: setting an oil Outlet Outlet, namely setting four oil outlets of OUT1, OUT2, OUT3 and OUT4 respectively, wherein the relative atmospheric pressure is 0 Pa;
step E4: clicking Interface, selecting Fluid by Interface Type, selecting INTERFACE1 by Interface side1 option, selecting INTERFACE2 by Interface side2 option, selecting GGI grid splicing mode by Rtta Axis rotation Axis as Z Axis, and selecting periodic axial and radial tolerance as 0.0001;
step E5: setting a Wall surface Wall, wherein the upper surface of an oil film is a rotary working surface and is set as a Rotate Wall, and the surface of the lower surface of the oil film, which is in contact with an oil pad, and the inner Wall surface of an oil cavity are both set as walls;
step E6: setting of ROTATE: selecting rolling Wall by an Option, Rotating to rotate axis to select Z axis global Z, setting the Rotating speed of the rotary worktable around the Z axis to be 100r/min, and finishing the setting of oil film boundary conditions;
step F: the convergence criterion used is the default criterion of the ANSYS CFX software, i.e. the residual value is already below 10-4And the fluctuation condition is stable without obvious downward trend.
CN201910332243.8A 2019-04-23 2019-04-23 Oil film lubrication performance simulation method under static pressure thrust bearing unbalance loading working condition Active CN110083924B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910332243.8A CN110083924B (en) 2019-04-23 2019-04-23 Oil film lubrication performance simulation method under static pressure thrust bearing unbalance loading working condition

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910332243.8A CN110083924B (en) 2019-04-23 2019-04-23 Oil film lubrication performance simulation method under static pressure thrust bearing unbalance loading working condition

Publications (2)

Publication Number Publication Date
CN110083924A CN110083924A (en) 2019-08-02
CN110083924B true CN110083924B (en) 2022-05-27

Family

ID=67416395

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910332243.8A Active CN110083924B (en) 2019-04-23 2019-04-23 Oil film lubrication performance simulation method under static pressure thrust bearing unbalance loading working condition

Country Status (1)

Country Link
CN (1) CN110083924B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111666644A (en) * 2020-06-18 2020-09-15 哈尔滨工业大学 Performance analysis method for parameterized CFD aerostatic bearing
CN113158369B (en) * 2021-04-19 2023-11-28 哈尔滨理工大学 Oil film flow simulation monitoring method for oil sealing edge of oil pad of hydrostatic thrust bearing
CN113742978A (en) * 2021-11-08 2021-12-03 哈尔滨理工大学 Friction failure prediction method for oil pad inclinable hydrostatic thrust bearing under extreme working condition

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102141084A (en) * 2011-03-30 2011-08-03 哈尔滨理工大学 Numerical simulation method for temperature and thickness relation of static thrust bearing gap oil film based on film thickness variation
US9169728B2 (en) * 2011-12-08 2015-10-27 General Electric Company Dynamic load reduction system
CN106594076A (en) * 2016-12-07 2017-04-26 北京工业大学 Dynamic and static pressure mixed fan-shaped oil pad
CN105135197B (en) * 2015-09-02 2017-09-05 哈尔滨理工大学 Double square chamber hydrostatic thrust bearing greasy property forecasting procedure
CN109117522A (en) * 2018-07-24 2019-01-01 北京工业大学 A kind of calculation method of the static-pressure slide carriage fluid structurecoupling working performance based on MATLAB-ANSYS software

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9371835B2 (en) * 2013-07-19 2016-06-21 Praxair Technology, Inc. Coupling for directly driven compressor

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102141084A (en) * 2011-03-30 2011-08-03 哈尔滨理工大学 Numerical simulation method for temperature and thickness relation of static thrust bearing gap oil film based on film thickness variation
US9169728B2 (en) * 2011-12-08 2015-10-27 General Electric Company Dynamic load reduction system
CN105135197B (en) * 2015-09-02 2017-09-05 哈尔滨理工大学 Double square chamber hydrostatic thrust bearing greasy property forecasting procedure
CN106594076A (en) * 2016-12-07 2017-04-26 北京工业大学 Dynamic and static pressure mixed fan-shaped oil pad
CN109117522A (en) * 2018-07-24 2019-01-01 北京工业大学 A kind of calculation method of the static-pressure slide carriage fluid structurecoupling working performance based on MATLAB-ANSYS software

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Yu Xiaodong 等."Research on the pressure field of hydrostatic center rest".《Proceedings of 2011 International Conference on Electronic & Mechanical Engineering and Information Technology》.2011, *
于晓东 等."极端工况双矩形腔静压推力轴承动态特性".《力学学报》.2018,第50卷(第4期), *
于晓东 等."极端工况静压推力轴承承载性能动压补偿".《推进技术》.2018,第39卷(第5期), *
于晓东等.环形腔多油垫静压推力轴承热变形.《吉林大学学报(工学版)》.2015,第45卷(第2期), *

Also Published As

Publication number Publication date
CN110083924A (en) 2019-08-02

Similar Documents

Publication Publication Date Title
CN110083924B (en) Oil film lubrication performance simulation method under static pressure thrust bearing unbalance loading working condition
CN102141084B (en) Numerical simulation method for temperature and thickness relation of static thrust bearing gap oil film based on film thickness variation
CN107066750B (en) Method for solving electromechanical comprehensive pipeline arrangement
US11200355B2 (en) 3D geometry generation for computer aided design considering subtractive manufacturing forces
CN107341288B (en) Optimization method for controlling vibration of combined cycle unit by adjusting elevation of bearing
CN102063549B (en) Method for stimulating fluid-solid-heat coupling value of machine tool spindle
CN101515168B (en) Method for generating ring-cutting roughing NC tool path of product triangular mesh model
CN102194018A (en) System and method for predictive analysis of energy consumption of building based on cloud computing
CN115937467B (en) Method and system for dividing random three-dimensional fracture in upscale model grid
CN103869755A (en) Method for guaranteeing smoothness of machine tool output power by adjusting machining code feeding speed
CN114082988A (en) Method for repairing aero-engine cold and hot end blade
CN112949112B (en) Rotor-sliding bearing system lubrication basin dynamic grid parallel computing method
CN111209639A (en) Efficient quantitative modeling method for impeller-bearing-rotor system
CN109558614B (en) Simulation method and system for gas flow in shale gas reservoir multi-scale fracture
Kalbfleisch et al. Computational valve plate design in axial piston pumps/motors
CN112711889B (en) Dynamic heat-fluid-solid simulation method based on LS-DYNA software
CN111428359B (en) Anisotropic rock modeling method considering grain occlusion
CN115130352B (en) Automatic modeling and analyzing method for complex pipeline system of rocket engine
CN113254525B (en) Collaborative manufacturing engine system constructed based on industrial PAAS and implementation method
CN114996917A (en) Single-disc crack rotor dynamic characteristic evaluation method based on similarity theory
CN111985095A (en) Assembly process parameter optimization method
Zloto et al. Modeling the load of the kinematic pair piston-cylinder in an axial piston pump by means of FEA
CN116522717B (en) Variable thickness thin layer unit modeling method for bolt connection structure of mounting edge of casing
CN112100755B (en) Numerical simulation method for asymmetric supporting structure containing rolling bearing gap
CN108628250A (en) A kind of optimization method of machine tool chief axis and handle of a knife fit structure size

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant