CN116484755A - Inner curve hydraulic motor roller optimal design method based on elastic flow lubrication - Google Patents

Inner curve hydraulic motor roller optimal design method based on elastic flow lubrication Download PDF

Info

Publication number
CN116484755A
CN116484755A CN202310057286.6A CN202310057286A CN116484755A CN 116484755 A CN116484755 A CN 116484755A CN 202310057286 A CN202310057286 A CN 202310057286A CN 116484755 A CN116484755 A CN 116484755A
Authority
CN
China
Prior art keywords
roller
film thickness
oil film
logarithmic
contact pressure
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
Application number
CN202310057286.6A
Other languages
Chinese (zh)
Other versions
CN116484755B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202310057286.6A priority Critical patent/CN116484755B/en
Publication of CN116484755A publication Critical patent/CN116484755A/en
Application granted granted Critical
Publication of CN116484755B publication Critical patent/CN116484755B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Geometry (AREA)
  • Fluid Mechanics (AREA)
  • Computing Systems (AREA)
  • Evolutionary Computation (AREA)
  • Operations Research (AREA)
  • Computer Hardware Design (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Support Of The Bearing (AREA)
  • Rolling Contact Bearings (AREA)

Abstract

The invention discloses an inner curve hydraulic motor roller optimal design method based on elastic flow lubrication. An improved roller busbar logarithmic modification equation is introduced first, the singularity of the traditional equation at the roller end is eliminated, and then two variables are introduced into the modification equation as optimal design parameters. And calculating the contact pressure distribution between the roller and the guide rail according to the elastic mechanics theory, then dispersing the roller into a plurality of slices along the axial direction, and calculating the thickness of an oil film between each slice and the guide rail based on the elastic flow lubrication theory. And finally, optimizing the logarithmic contour of the roller with the aim of reducing the contact pressure, increasing the oil film thickness, improving the contact pressure and the oil film thickness distribution uniformity and improving the bearing capacity, and finally obtaining the optimal logarithmic repair type roller design parameters.

Description

Inner curve hydraulic motor roller optimal design method based on elastic flow lubrication
Technical Field
The invention belongs to the technical fields of tribology and fluid transmission, and particularly relates to an inner curve hydraulic motor roller optimization design method based on elastohydrodynamic lubrication.
Background
The roller in the inner curve hydraulic motor is one of the most important transmission parts of the motor, is in a line contact state with the inner curve guide rail, and converts hydraulic pressure into mechanical torque for pushing the motor to rotate through rolling. However, the straight busbar roller is easy to generate stress concentration at two ends of the roller under heavy load, and the fatigue life of the roller is affected. Lundberg proposed a roller profile logarithmic modification method, which has been considered as one of the most optimal roller modification methods in the dry friction state. However, in practical application, the roller and the guide rail thereof are often in a spring-flow lubrication state, and the edge effect at two ends of the roller cannot be reduced to the greatest extent by the conventional modification mode. In addition, the edge stress cannot be completely eliminated due to the fact that the roller is too small in modification amount, and the contact stress at the center of a roller bus is increased due to the fact that the roller is too large in modification amount. Therefore, there is a problem of optimizing design in the modification amount design.
The calculation of the film thickness of the elastic lubricating oil of the roller and the guide rail thereof usually adopts an analytical empirical formula or a two-dimensional Reynolds equation numerical method. However, the analysis of the empirical formula treats the roller as a wireless length, the axial oil film thickness is constant, and the axial oil film distribution under the influence of stress concentration cannot be described. The numerical solution of the two-dimensional Reynolds equation is too complex, a great deal of time is consumed in the optimization process, and the optimization calculation efficiency is low.
In summary, a simple and effective method is provided to calculate the thickness of the lubricant film between the roller and the guide rail, and to optimally design the bus profile.
Disclosure of Invention
The invention aims to provide a simple and effective method for calculating the thickness of the lubricating oil film between a roller and a guide rail, optimally designing the contour of a bus, improving the bearing capacity and the fatigue life of the roller, and reducing friction and abrasion. Has the advantages of high efficiency, low cost, strong practicability and the like. Can be combined with a test method to guide the test, and save the test cost.
The invention is realized by the following technical scheme: the method is used for optimizing the design of the contour of the bus of the limited long-line contact friction pair of the inner curve hydraulic motor, and specifically comprises the following steps:
(1) According to the logarithmic modification equation, two variables are introduced into the coefficient and logarithmic term of the logarithmic modification equation, so that the singularity of the roller end is eliminated;
(2) According to the elastic mechanics theory, a conjugate gradient method is adopted to calculate the contact pressure distribution between the roller and the guide rail thereof, the bus profile after the modification of the modified roller pair modification equation in the step (1) is considered during calculation, and the working condition of the roller is considered;
(3) Dispersing the roller into a plurality of slices along the axial direction, and calculating the contact pressure between each slice and the guide rail through the step (2); according to the elastic fluid dynamic pressure lubrication theory, the load parameter in the oil film thickness calculation formula is calculated through the contact pressure of each slice, and then the oil film thickness of each slice in contact with the guide rail is calculated;
(4) Taking two variables introduced in the improved logarithmic modification equation in the step (1) as optimization design parameters of the inner curve hydraulic motor roller, so as to reduce contact pressure, increase oil film thickness, improve uniformity of contact pressure and oil film thickness distribution and improve bearing capacity as targets, optimizing the logarithmic profile of the roller, wherein the contact pressure and the oil film thickness distribution are respectively obtained through calculation in the step (2) and the step (3); and finally obtaining the optimal logarithmic repair type roller design parameters.
Further, according to the logarithmic modification proposed by Lundberg, the modified logarithmic modification equation is:
wherein C is k (y) is the amount of modification of the roller in the axial direction, y is the roller axisThe axial coordinates are directed, and the origin of the coordinates is the axial center of the roller; p is p H The Hertz contact pressure between the roller and the guide rail under the current working condition is represented by E, the equivalent elastic modulus is represented by Σρ, the sum of contact curvatures of the roller and the guide rail is represented by l, the length of the roller is represented by k 1 And k 2 The initial values of the optimal design variables of the roller shaping profile are all set to be 1.
Further, the minimum oil film thickness calculation formula is:
wherein H is m Is a dimensionless film thickness parameter, h m For actual oil film thickness, R x Is equivalent radius of curvature; g * 、U * The non-dimensional material parameters and the non-dimensional speed parameters are respectively.
Further, in the step (4), the logarithmic contour optimization objective function of the roller is:
wherein p is max And p c Maximum contact pressure and center pressure, h, respectively, of the axial contact length of the roller min And h c The minimum oil film thickness and the center oil film thickness of the axial contact length of the roller are respectively.
The invention has the beneficial effects that:
according to the non-Hertz contact stress distribution between the roller and the guide rail, an oil film thickness calculation formula between the roller and the guide rail is deduced, so that a complex numerical iteration process of an oil film thickness numerical calculation method is avoided, and the defect that the axial oil film distribution cannot be described by a traditional analysis method is overcome. The optimal logarithmic repair type roller design parameters are obtained by optimizing the logarithmic contour of the roller, so that the effects of reducing contact pressure, increasing oil film thickness, improving the uniformity of contact pressure and oil film thickness distribution and improving bearing capacity are achieved.
Drawings
FIG. 1 is a flow chart of an inner curve motor roller optimization design.
Fig. 2 is a partial schematic view of a roller-rail in an inner curve hydraulic motor.
FIG. 3 is a graph comparing contact pressure between the optimized front and rear roller tracks.
FIG. 4 is a graph of oil film thickness comparisons between the optimized front and rear roller tracks.
Detailed Description
The invention will be described in detail with reference to the accompanying drawings and specific examples, as shown in fig. 1 and 2, comprising the following steps:
(1) The roller end singularity is eliminated by improving the logarithmic modification equation according to Lundberg. Two variables were introduced into the modified logarithmic modification equation, with the initial value set to 1.
(2) According to the elastic mechanics theory, a conjugate gradient method is adopted to calculate the contact pressure distribution between the roller and the guide rail thereof, and the roller modification equation established in the step (1) is considered during calculation, and meanwhile, the working condition of the roller is considered.
(3) The roller is axially discretized into a plurality of slices, and the contact pressure between each slice and the guide rail is calculated in the step (2). According to the elastic fluid dynamic pressure lubrication theory, the load parameter in the oil film thickness calculation formula is calculated through the contact stress of each slice, and then the oil film thickness of each slice in contact with the guide rail is calculated.
(4) And (3) optimizing the logarithmic contour of the roller with the two variables determined in the step (1) as optimization design parameters and aiming at reducing the contact pressure, increasing the oil film thickness, improving the uniformity of the contact pressure and the oil film thickness distribution and improving the bearing capacity, wherein the contact pressure and the oil film thickness distribution are calculated by the step (2) and the step (3), respectively. And finally obtaining the optimal logarithmic repair type roller design parameters.
To eliminate the singularities at both ends of the logarithmic repair roller, the improved logarithmic repair equation is:
wherein C is k And (y) is the axial modification quantity of the roller, y is the axial coordinate of the roller, and the origin of the coordinate is the axial center of the roller. P is p H The Hertz contact pressure between the roller and the guide rail under the current working condition is represented by E, the equivalent elastic modulus is represented by Σρ, the sum of contact curvatures of the roller and the guide rail is represented by l, the length of the roller is represented by k 1 And k 2 The design variable is optimized for the roller shaping profile.
To calculate the contact pressure between the roller and the rail, the Boussineq equation describing the relationship between contact stress and contact deformation in infinite half space needs to be solved, namely:
where w (x, y) is the contact deformation at the coordinates (x, y), p (ζ, η) is the contact pressure at the coordinates (ζ, η), v is the poisson's ratio of the material, and E is the elastic modulus. Furthermore, the gap h (x, y) between the roller and the rail can be described as:
h(x,y)=B(x,y)+w(x,y)-δ
where δ is the relative displacement of the roller and rail under external load and B (x, y) is the initial gap between the roller and rail:
wherein R is 1 And R is 2 The radii of the rollers and the guide rail are respectively, and beta is the relative deflection angle between the rollers and the guide rail. In order to balance with the external load W, the contact pressure between the roller and the rail needs to satisfy the balance equation:
wherein A is c Is the actual contact area of the roller and the rail. There are two additional conditions:
when h (x, y) is greater than 0, the rollers are not in contact with the guide rail, and the contact pressure is 0; when h (x, y) is less than 0, the rollers and the guide rail are in contact, and the contact pressure is greater than 0. According to the equation, the contact pressure is calculated iteratively by adopting a conjugate gradient method, and the contact deformation calculation in iteration adopts a discrete convolution-fast Fourier transform (DC-FFT) method.
The minimum oil film thickness calculation formula is:
H m =h m /R x =2.65G *0.54 U *0.7 W *(-0.13)
wherein H is m Is a dimensionless film thickness parameter, h m For actual oil film thickness, R x Is the equivalent radius of curvature. G * 、U * 、W * The non-dimensional material parameter, the non-dimensional speed parameter and the non-dimensional load parameter are respectively. Non-dimensional oil film thickness parameter W * The method comprises the following steps: w (W) * =W/ER x l. According to the Hertz contact theory, the relationship between the external load and the maximum contact stress of the roller is as follows:
substituting to calculate dimensionless oil film thickness parameter W * . Substituting the contact pressure of each slice of the roller axial direction and the guide rail into p H The minimum oil film thickness of the slice and the guide rail can be calculated.
The roller modification parameter optimization aims at reducing contact pressure, increasing oil film thickness, improving contact pressure and oil film thickness distribution uniformity and improving bearing capacity, and the objective function is as follows:
wherein p is max And p c Maximum contact pressure and center pressure, h, respectively, of the axial contact length of the roller min And h c The minimum oil film thickness and the center oil film thickness of the axial contact length of the roller are respectively.
Finally with f (k) 1 ,k 2 ) Minimum value for objective function, k 1 ,k 2 And for designing variables, adopting a genetic algorithm to solve the optimal logarithmic modification parameters, and improving the bearing capacity and the lubrication state of the roller.
Further described in one embodiment, the roller radius R 1 7mm, a roller length l of 26mm, a rail diameter R 2 45mm, the external load of the roller is 10kN, the relative deflection angle between the roller and the guide rail is 0.1 DEG, a logarithmic modification mode is adopted, and the initial modification parameter k is the initial modification parameter 1 ,k 2 All set to 0. Optimized modification parameter k 1 ,k 2 1.38 and 0.84, respectively. The initial parameters gave an objective function result of 0.0628 and the optimized parameters gave an objective function result of 0.0254. The contact pressure and the oil film thickness of the two groups of parameters are respectively shown in fig. 3 and 4, and compared with an unmodified roller, the initial modification parameters can effectively reduce the contact pressure concentration caused by the edge effect and improve the oil film thickness of the edge of the roller, but the contact pressure and the oil film thickness are uneven. The optimized modification parameters not only reduce the edge contact pressure, but also make the pressure distribution and oil film thickness distribution of the whole contact length more uniform, and improve the bearing capacity of the roller.
The above-described embodiments are intended to illustrate the present invention, not to limit it, and any modifications and variations made thereto are within the spirit of the invention and the scope of the appended claims.

Claims (4)

1. The optimization design method for the inner curve hydraulic motor roller based on the elastic flow lubrication is characterized by being used for optimizing the contour of the limited long line contact friction pair bus of the inner curve hydraulic motor, and specifically comprising the following steps:
(1) According to the logarithmic modification equation, two variables are introduced into the coefficient and logarithmic term of the logarithmic modification equation, so that the singularity of the roller end is eliminated;
(2) According to the elastic mechanics theory, a conjugate gradient method is adopted to calculate the contact pressure distribution between the roller and the guide rail thereof, the bus profile after the modification of the modified roller pair modification equation in the step (1) is considered during calculation, and the working condition of the roller is considered;
(3) Dispersing the roller into a plurality of slices along the axial direction, and calculating the contact pressure between each slice and the guide rail through the step (2); according to the elastic fluid dynamic pressure lubrication theory, the load parameter in the oil film thickness calculation formula is calculated through the contact pressure of each slice, and then the oil film thickness of each slice in contact with the guide rail is calculated;
(4) Taking two variables introduced in the improved logarithmic modification equation in the step (1) as optimization design parameters of the inner curve hydraulic motor roller, so as to reduce contact pressure, increase oil film thickness, improve uniformity of contact pressure and oil film thickness distribution and improve bearing capacity as targets, optimizing the logarithmic profile of the roller, wherein the contact pressure and the oil film thickness distribution are respectively obtained through calculation in the step (2) and the step (3); and finally obtaining the optimal logarithmic repair type roller design parameters.
2. The optimization design method for the inner curve hydraulic motor roller based on the elastic flow lubrication according to claim 1 is characterized in that according to the logarithmic modification proposed by Lundberg, the improved logarithmic modification equation is as follows:
wherein C is k (y) is the amount of modification of the roller in the axial direction, y is the roller axisThe axial coordinates are directed, and the origin of the coordinates is the axial center of the roller; p is p H The Hertz contact pressure between the roller and the guide rail under the current working condition is represented by E, the equivalent elastic modulus is represented by Σρ, the sum of contact curvatures of the roller and the guide rail is represented by l, the length of the roller is represented by k 1 And k 2 The initial values of the optimal design variables of the roller shaping profile are all set to be 1.
3. The optimal design method for the inner curve hydraulic motor roller based on the elastic flow lubrication according to claim 2 is characterized in that a minimum oil film thickness calculation formula is as follows:
wherein H is m Is a dimensionless film thickness parameter, h m For actual oil film thickness, R x Is equivalent radius of curvature; g * 、U * The non-dimensional material parameters and the non-dimensional speed parameters are respectively.
4. The optimization design method for the inner curve hydraulic motor roller based on the elastic flow lubrication according to claim 2, wherein in the step (4), the logarithmic contour optimization objective function of the roller is as follows:
wherein p is max And p c Maximum contact pressure and center pressure, h, respectively, of the axial contact length of the roller min And h c The minimum oil film thickness and the center oil film thickness of the axial contact length of the roller are respectively.
CN202310057286.6A 2023-01-18 2023-01-18 Inner curve hydraulic motor roller optimal design method based on elastic flow lubrication Active CN116484755B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310057286.6A CN116484755B (en) 2023-01-18 2023-01-18 Inner curve hydraulic motor roller optimal design method based on elastic flow lubrication

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310057286.6A CN116484755B (en) 2023-01-18 2023-01-18 Inner curve hydraulic motor roller optimal design method based on elastic flow lubrication

Publications (2)

Publication Number Publication Date
CN116484755A true CN116484755A (en) 2023-07-25
CN116484755B CN116484755B (en) 2023-10-24

Family

ID=87225639

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310057286.6A Active CN116484755B (en) 2023-01-18 2023-01-18 Inner curve hydraulic motor roller optimal design method based on elastic flow lubrication

Country Status (1)

Country Link
CN (1) CN116484755B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102644663A (en) * 2012-04-26 2012-08-22 北京石油化工学院 Cylindrical roller multi-circular-arc variable curvature profile engineering simulation method
CN110765617A (en) * 2019-10-25 2020-02-07 常州市乾憬轴承科技有限公司 Lubrication theory-based cylindrical roller bearing roller logarithmic modification design method
CN111919040A (en) * 2018-03-28 2020-11-10 Ntn株式会社 Tapered roller bearing

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102644663A (en) * 2012-04-26 2012-08-22 北京石油化工学院 Cylindrical roller multi-circular-arc variable curvature profile engineering simulation method
CN111919040A (en) * 2018-03-28 2020-11-10 Ntn株式会社 Tapered roller bearing
CN110765617A (en) * 2019-10-25 2020-02-07 常州市乾憬轴承科技有限公司 Lubrication theory-based cylindrical roller bearing roller logarithmic modification design method

Also Published As

Publication number Publication date
CN116484755B (en) 2023-10-24

Similar Documents

Publication Publication Date Title
CN108984933B (en) Boundary element method for calculating load and pressure of rolling bearing under elastohydrodynamic lubrication condition
CN110765617B (en) Design method for cylindrical roller bearing roller logarithmic modification based on lubrication theory
CN106649982B (en) Method for calculating friction torque of large-cone-angle conical main bearing of wind driven generator
CN103962390B (en) Tension force and draught pressure synthetic setting method in VC roller planisher wet jetting piles process
CN109753723B (en) Radial rolling bearing fatigue life calculation method
CN103810354A (en) Optimal design method for logarithm shaping curve of cylindrical roller bearing
CN116484755B (en) Inner curve hydraulic motor roller optimal design method based on elastic flow lubrication
CN107563082B (en) Bearing parameter optimization method based on cylindrical roller bearing contact deformation and load distribution
CN116000106A (en) Rolling force setting method in cold continuous rolling speed increasing and decreasing stage
CN110688715A (en) Variable-crown roller type self-adaptive design method based on genetic algorithm
CN111140598A (en) Shape-modifying roller of self-aligning bearing for wind power equipment
CN111291455A (en) Shape modification design method of self-aligning bearing roller for wind power equipment
CN113094629A (en) Rolling force setting method for asynchronous rolling of steel strip
Jiqiang et al. Geometric modification on mixed lubrication performance of low-speed cylindrical roller bearing
Hrček et al. Increase of durability in the area of line contact with curved body geometry
CN112836445B (en) Method for solving dynamic contact mixed lubrication parameters of rolling bearing
CN115510735A (en) Rail grinding target profile optimization design method based on wheel-rail contact parameters
CN111709098B (en) Method for reducing friction coefficient of bearing by changing lubricating oil under variable load
Jiang et al. Asymmetric loading multi-roller planetary traction drive: Modeling and performance analysis
CN207111258U (en) The attachment structure of piston rod
Ishmuratov et al. Research wear resistance teeth of gears at rolling
CN112069706A (en) Tribology and dynamics coupling analysis method for tank turret seat ring
CN109241695B (en) Rolling bearing reliability multi-objective optimization design method based on novel PSO-SADE
CN112347581B (en) Matching method and matching system for RV reducer ceramic ball main bearing and storage medium
CN115630554B (en) Wheel-rail two-dimensional rolling contact force calculation method considering material nonuniformity

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