CN110083924A - Film lubrication performance simulation method under hydrostatic thrust bearing unbalance loading operating condition - Google Patents
Film lubrication performance simulation method under hydrostatic thrust bearing unbalance loading operating condition Download PDFInfo
- Publication number
- CN110083924A CN110083924A CN201910332243.8A CN201910332243A CN110083924A CN 110083924 A CN110083924 A CN 110083924A CN 201910332243 A CN201910332243 A CN 201910332243A CN 110083924 A CN110083924 A CN 110083924A
- Authority
- CN
- China
- Prior art keywords
- block
- oil
- option
- setting
- film
- 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.)
- Granted
Links
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
Oil film performance simulation method under hydrostatic thrust bearing unbalance loading operating condition, the invention mainly relates to the analogy methods of film lubrication performance under double square chamber hydrostatic thrust bearing unbalance loading operating condition.The invention solves hydrostatic thrust bearings under unbalance loading operating condition without actually active method come the problem of simulating film lubrication performance.This method is to pass through: 1, establishing 6 symmetrical double square lubricating pad models by UG；2, model is imported in ANSYS ICEM CFD software, creates Fuel film model part, create topological structure partition structure grid；3, check whether the Points And lines of model are all associated with the point of block with line；4, grid division number of nodes, the encryption of oil inlet number of nodes；5, boundary film condition is set in ANSYS CFX software；6, convergence uses the standard of ANSYS CFX software default, judges whether fluctuation is stable.The present invention is applied to greasy property simulation field under hydrostatic thrust bearing unbalance loading operating condition.
Description
Technical field
It is inclined in particular for double square chamber the present invention relates to a kind of analogy method of hydrostatic thrust bearing film lubrication performance
The analogy method of film lubrication performance when load.
Background technique
Hydrostatic thrust bearing is one of the kernel component of heavy type numerical control process equipment, its performance is better, and equipment adds
Working medium amount and operational efficiency are higher.In the ideal situation, the mass center of workpiece should carry center at the center of rotary table
Load, hydrostatic thrust bearing working condition is in optimum performance state at this time.But in the actual processing process, the mass center of workpiece can
The case where unbalance loading can not occur at the center of rotary table, keep hydrostatic thrust bearing oil film locally thinning, or even occur dry
Friction, so as to cause hydrostatic thrust bearing failure, this problem restricts China's heavy type numerical control process equipment in high precision, efficiently
The direction of rate is developed.And both at home and abroad be directed to unbalance loading when hydrostatic thrust bearing film lubrication behavior pattern analogy method research compared with
It is few, therefore it is very necessary for studying new hydrostatic thrust bearing analogy method of greasy property in unbalance loading operating condition, is also quiet
The design that advanced optimizes of pressure thrust bearing provides foundation.
Summary of the invention
The purpose of the present invention is to solve hydrostatic thrust bearing is in practical work process, workpiece mass center and revolution
Work table rotation center, which is not overlapped, there is unbalance loading, and oil film thickness becomes uneven distribution, so as to cause hydrostatic thrust bearing lubricity
It can be greatly reduced or even be occurred dry friction, lack a kind of effective method come the problem of simulation to this phenomenon.In order to
The greasy property for probing into hydrostatic thrust bearing oil film in unbalance loading proposes a kind of hydrostatic bearing film lubrication performance in unbalance loading
Analogy method.
The analogy method of greasy property when hydrostatic thrust bearing unbalance loading is achieved through the following technical solutions.
(1), 6 symmetrical double square lubricating pad threedimensional models are established by UG 3 d modeling software；
(2), model is imported in ANSYS ICEM CFD software, face required for splitting creates Fuel film model part, and creation is opened up
Flutter structure partition structure grid；
(3), check whether the Points And lines of model are all associated with the point of block with line, if block all lines are all green
Color then all in association, otherwise without all associations, therefore need to be associated by above-mentioned steps again；
(4), grid division number of nodes；The encryption of oil inlet number of nodes, adjacent segments points are no more than twice, guarantee grid matter with this
Amount；
(5), boundary film condition, the setting of setting and domain including viscous kelvin relation are set in ANSYS CFX software；
(6), convergence uses the standard of ANSYS CFX software default, i.e. residual values are lower than 10-4, judge whether fluctuation is stable, has
Without apparent downward trend.
The present invention can simulate actual condition of the hydrostatic thrust bearing when unbalance loading occurs for workpiece by computer, obtain
The film lubrication parameter of hydrostatic thrust bearing under different unbalance loading degree.This method can save the plenty of time, obtain and tie with experiment
The film lubrication parameter that fruit is consistent.By studying lubricating oil film parameter, different unbalance loading degree are disclosed to film lubrication characteristic
Affecting laws.This provides effective reference frame for the further optimizing research of hydrostatic thrust bearing, also to use static pressure
The heavy equipment of thrust bearing realizes high efficiency, has established solid foundation in high precision.
Detailed description of the invention
Fig. 1 is film lubrication performance simulation method flow diagram under hydrostatic thrust bearing unbalance loading operating condition.
Fig. 2 is the Fuel film model under the unbalance loading operating condition that step A is mentioned.
Fig. 3 is list after the definition Part that step B1 is mentioned.
Fig. 4 is the oil film network topology illustraton of model that step B2 is mentioned.
Fig. 5 is the oil film grid dividing figure that step C is mentioned.
Fig. 6 is the oil film mesh quality schematic diagram that step D is mentioned.
Fig. 7 is the oil film each position setting that step E is mentioned.
Fig. 8 is the iteration residual plot that step F is mentioned.
Specific embodiment
It can be achieved through the following technical solutions,
Step A: the lubricating pad under unbalance loading operating condition is modeled using UG software, is created that different shape size according to different operating conditions
Model.
Step B1: the Fuel film model under established different operating conditions is exported as into .x.t file using UG software, will be exported
.x.t file imported into ANSYS ICEM CFD software, setting unit mm；Since CFD software can not will be established in UG
Face is divided into two faces automatically, so what is done first after model importing is that the face in UG is split as required face；Secondly
Be create Fuel film model part, the step for mainly for import CFX back boundary condition setting, part be set with IN1,
IN2, OUT1, OUT2, OUT3, OUT4, INTERFACE1, INTERFACE 2, these definition faces ROTATE, WALL.
Step B2: top-down division methods are used when creation topological structure partition structure grid.First click on option
Blocking selects Check Block option, selects the subtopology structure in split block, first draws to X-direction block
Point, then Y, Z-direction block are divided, delete unnecessary topological structure.When creation point mapping, by each top of threedimensional model
Point is successively mapped in the vertex of block, and the Auto association that line mapping can directly select in Blocking is closed automatically
Join option, completes the mapping of oil film and oil pocket Points And lines.
Step B3: the shape of oil inlet is cylinder, then grid should use O-shaped grid, individually be divided to its block,
Two points are dosed at the symmetrical position of circle first, the block of X, Y-direction is divided, then select Check Block option
By following rectangular block it is stretched come, and then obtain the block of rectangular shape, when mapping, quadrangle be mapped to
Circle.The Select face option in Split Block is finally clicked, the top and bottom of cylinder are selected, clicks and determines, so far complete
All block are divided and dotted line mapping.
Step C: checking whether the Points And lines of model are all associated with the Points And lines of block, will be under Geometry
Point and edge hides, and checks whether the line of block is all green.If it is, all in association；If somewhere is not
Be it is green, show herein without that therefore need to be associated again by above-mentioned steps in association.
Step D: Mesh Quality module is clicked, oil film grid and angle mesh quality are checked.Export grid
Shi Caiyong Output/Select Solver selects ANSYS CFX solver, and default output file format is .CFX5.
Step E1: selection ANSYS CFX software carries out simulation analysis to oil film.Fluid domain is set as FLUID, is selected
Material is 32# lubricating oil, and fluid calculation model is selected as total energy, temperature setting 293K, other settings are kept
Default setting.
Step E2: the setting of oil inlet Inlet.The flow of inlet porting IN1 and IN2 are 0.035, and temperature is given as room temperature
293K i.e. 20 DEG C.
Step E3: four oil outlets, respectively OUT1, OUT2, OUT3, OUT4, phase is arranged in the setting of oil outlet Outlet
It is 0Pa to atmospheric pressure.
Step E4: clicking Interface, and Interface Type selects Fluid, interface Interface side1 choosing
Item selection INTERFACE1, interface Iterfance side2 select option INTERFACE2, and Rtation Axis rotary shaft is
Z axis, selects GGI grid connecting method in Mesh connecion option, and periodical axial-radial tolerance is 0.0001.
Step E5: the setting of wall surface Wall, oil film upper surface are rotating face, are set as Rotate Wall, under oil film
The face and oil pocket inner wall that surface is in contact with lubricating pad are both configured to Wall.
The setting of step E6:ROTATE: Option option selects Rotating Wall, rotates to rotate axis selection
Z axis global Z, revolving speed is 100r/min to setting rotary table about the z axis, completes the setting of Oil Boundary condition.
Step F: the convergence of use is the standard of ANSYS CFX software default, i.e., residual values are already below 10-4, and
And fluctuation situation is stablized, without apparent downward trend.
Claims (1)
1. the analogy method of film lubrication performance under hydrostatic thrust bearing unbalance loading operating condition, it is characterised in that hydrostatic thrust bearing unbalance loading
Operating condition Imitating method is realized according to the following steps:
Step A: the lubricating pad under unbalance loading operating condition is modeled using UG software, is created that different shape size according to different operating conditions
Model；
Step B1: exporting as .x.t file for the Fuel film model under established different operating conditions using UG software, will be derived
.x.t file is imported into ANSYS ICEM CFD software, setting unit mm；The face that can not will be established in UG due to CFD software
Automatically two faces are divided into, so what is done first after model importing is that the face in UG is split as required face；Followed by
Create the part of Fuel film model, the step for be setting in order to import CFX back boundary condition, part be set with IN1, IN2,
OUT1, OUT2, OUT3, OUT4, INTERFACE1, INTERFACE 2, these definition faces ROTATE, WALL；
Step B2: top-down division methods are used when creation topological structure partition structure grid, first click on option
Blocking selects Check Block option, selects the subtopology structure in split block, first draws to X-direction block
Point, then Y, Z-direction block are divided, unnecessary topological structure is deleted, when creation point maps, by each top of threedimensional model
Point is successively mapped in the vertex of block, and the Auto association that line mapping can directly select in Blocking is closed automatically
Join option, completes the mapping of oil film and oil pocket Points And lines；
Step B3: the shape of oil inlet is cylinder, then grid should use O-shaped grid, individually be divided to its block, first
Two points are dosed at the symmetrical position of circle, the block of X, Y-direction is divided, then select Check Block option will under
The rectangular block in face it is stretched come, and then obtain the block of rectangular shape, when mapping, quadrangle be mapped to circle, most
The Select face option in Split Block is clicked afterwards, selects the top and bottom of cylinder, is clicked and is determined, so far completes all
Block is divided and dotted line mapping；
Step C: checking whether the Points And lines of model are all associated with the Points And lines of block, by under Geometry point and
Edge hides, and checks whether the line of block is all green, if it is, all in association；If somewhere is not green
, show in no association herein, therefore need to be associated again by above-mentioned steps；
Step D: clicking Mesh Quality module, check oil film grid and angle mesh quality, adopts when exporting grid
ANSYS CFX solver is selected with Output/Select Solver, default output file format is .CFX5；
Step E1: selection ANSYS CFX software carries out simulation analysis to oil film, sets fluid domain as FLUID, selects material
For 32# lubricating oil, fluid calculation model is selected as total energy, temperature setting 293K, other settings keep default
Setting；
Step E2: the setting of oil inlet Inlet, the flow of inlet porting IN1 and IN2 are 0.035, and temperature is given as room temperature 293K
I.e. 20 DEG C；
Step E3: four oil outlets, respectively OUT1, OUT2, OUT3, OUT4 are arranged in the setting of oil outlet Outlet, relatively large
Air pressure is 0Pa；
Step E4: clicking Interface, and Interface Type selects Fluid, the choosing of interface Interface side1 option
INTERFACE1 is selected, interface Iterfance side2 selects option INTERFACE2, and Rtation Axis rotary shaft is Z axis,
GGI grid connecting method is selected in Mesh connecion option, periodical axial-radial tolerance is 0.0001；
Step E5: the setting of wall surface Wall, oil film upper surface are rotating face, are set as Rotate Wall, oil film lower surface
The face and oil pocket inner wall being in contact with lubricating pad are both configured to Wall；
The setting of step E6:ROTATE: Option option selects Rotating Wall, rotates to rotate axis selection Z axis
Global Z, revolving speed is 100r/min to setting rotary table about the z axis, completes the setting of Oil Boundary condition；
Step F: the convergence of use is the standard of ANSYS CFX software default, i.e., residual values are already below 10-4, and wave
Emotionally condition is stablized, without apparent downward trend.
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 true CN110083924A (en) | 2019-08-02 |
| CN110083924B 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) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111666644A (en) * | 2020-06-18 | 2020-09-15 | 哈尔滨工业大学 | Performance analysis method for parameterized CFD aerostatic bearing |
| CN113158369A (en) * | 2021-04-19 | 2021-07-23 | 哈尔滨理工大学 | Oil film flow simulation monitoring method for oil seal edge of hydrostatic thrust bearing oil pad |
| CN113742978A (en) * | 2021-11-08 | 2021-12-03 | 哈尔滨理工大学 | Friction failure prediction method for oil pad inclinable hydrostatic thrust bearing under extreme working condition |
Citations (6)
| 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 |
| US20150023785A1 (en) * | 2013-07-19 | 2015-01-22 | Michael J. Stanko | Coupling for directly driven compressor |
| US9169728B2 (en) * | 2011-12-08 | 2015-10-27 | General Electric Company | Dynamic load reduction system |
| CN105135197A (en) * | 2015-09-02 | 2015-12-09 | 哈尔滨理工大学 | Lubricating property forecasting method of double-rectangular-cavity hydrostatic thrust bearing |
| 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 |
-
2019
- 2019-04-23 CN CN201910332243.8A patent/CN110083924B/en active Active
Patent Citations (6)
| 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 |
| US20150023785A1 (en) * | 2013-07-19 | 2015-01-22 | Michael J. Stanko | Coupling for directly driven compressor |
| CN105135197A (en) * | 2015-09-02 | 2015-12-09 | 哈尔滨理工大学 | Lubricating property forecasting method of double-rectangular-cavity hydrostatic thrust bearing |
| 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)
| Title |
|---|
| YU XIAODONG 等: ""Research on the pressure field of hydrostatic center rest"", 《PROCEEDINGS OF 2011 INTERNATIONAL CONFERENCE ON ELECTRONIC & MECHANICAL ENGINEERING AND INFORMATION TECHNOLOGY》, 19 September 2011 (2011-09-19), pages 458 - 461 * |
| 于晓东 等: ""极端工况双矩形腔静压推力轴承动态特性"", 《力学学报》, vol. 50, no. 4, 13 May 2018 (2018-05-13), pages 1085 - 1091 * |
| 于晓东 等: ""极端工况静压推力轴承承载性能动压补偿"", 《推进技术》, vol. 39, no. 5, 2 April 2018 (2018-04-02), pages 899 - 907 * |
| 于晓东等: "环形腔多油垫静压推力轴承热变形", 《吉林大学学报(工学版)》, vol. 45, no. 2, 15 March 2015 (2015-03-15), pages 125 - 130 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111666644A (en) * | 2020-06-18 | 2020-09-15 | 哈尔滨工业大学 | Performance analysis method for parameterized CFD aerostatic bearing |
| CN113158369A (en) * | 2021-04-19 | 2021-07-23 | 哈尔滨理工大学 | Oil film flow simulation monitoring method for oil seal edge of hydrostatic thrust bearing oil pad |
| 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 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110083924B (en) | 2022-05-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110083924A (en) | Film lubrication performance simulation method under hydrostatic thrust bearing unbalance loading operating condition | |
| CN102141084B (en) | Numerical simulation method for temperature and thickness relation of static thrust bearing gap oil film based on film thickness variation | |
| CN103020393B (en) | Aircraft assembly data visualization method based on Web | |
| CN109855238A (en) | A kind of modeling of central air-conditioning and efficiency optimization method and device | |
| CN107341288B (en) | Optimization method for controlling vibration of combined cycle unit by adjusting elevation of bearing | |
| CN112487734B (en) | CFD numerical simulation method for air system structure | |
| CN109471377A (en) | Based on two twin axis X-Y leather and fur products Cutting machine control methods of number | |
| CN105046432A (en) | Building dimensional guide construction method | |
| WO2023138013A1 (en) | Method and system for computing power gridding in multi-cloud scene | |
| CN105893695A (en) | Through-flow clearance error compensation method for small-size turbine in cylinder combination state | |
| CN104615801B (en) | The heavy critical lubricating status tachometer value of hydrostatic bearing determines method | |
| US6526550B1 (en) | Analyzing characteristics of geometries | |
| CN107526914B (en) | Variable-watershed flow field calculation method of tilting-pad sliding bearing based on structured dynamic grid | |
| CN113111453B (en) | Numerical simulation simplifying method for drag reduction performance of micro-texture of rotating blade | |
| CN112949112B (en) | Rotor-sliding bearing system lubrication basin dynamic grid parallel computing method | |
| CN109063223B (en) | BIM model light weight method and device and BIM model processing method and system | |
| CN114155341A (en) | Cross-platform consistency method for three-dimensional model of large-scale digital workshop based on Unity platform | |
| CN116702374A (en) | Turbine pump sealing structure design method, device and equipment | |
| CN102902769A (en) | Database benchmark test system of cloud computing platform and method thereof | |
| CN104598666A (en) | Large-scale hydrostatic support critical load parameter acquiring method | |
| CN112711889B (en) | Dynamic heat-fluid-solid simulation method based on LS-DYNA software | |
| CN114840946B (en) | Dynamic pressure type mechanical seal reliability calculation method based on liquid film vaporization phase transition | |
| CN115601548A (en) | Cutting method for CFD analysis rotating domain of rotating machine | |
| CN114528657A (en) | Energy loss dimensionality reduction coupling method and system, computer equipment and readable storage medium | |
| CN104731051A (en) | Method for achieving sensor point position objectification |
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 |