A heavy load sliding bearing self-adaptive wear calculation method and system
By using an adaptive wear calculation method, the problem of inaccurate wear prediction for sliding bearings under heavy load conditions has been solved, enabling more accurate wear process simulation and performance evaluation, and improving the equipment's maintenance and life prediction capabilities.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- XI AN JIAOTONG UNIV
- Filing Date
- 2026-02-10
- Publication Date
- 2026-06-02
Smart Images
Figure CN122133275A_ABST
Abstract
Description
Technical Field
[0001] This invention belongs to the field of sliding bearings, and specifically relates to an adaptive wear calculation method and system for heavy-duty sliding bearings. Background Technology
[0002] Sliding bearings are widely used in high-load and harsh-condition fields such as shipbuilding, aviation, aerospace, and nuclear power. These fields have extremely high requirements for the performance and reliability of mechanical equipment, and the performance of sliding bearings directly affects the stability, safety, and service life of the entire system. Currently, the operating conditions faced by sliding bearings in major equipment are becoming increasingly extreme. Under conditions such as start-up, shutdown, and heavy loads, bearings inevitably suffer wear, which may lead to a continuous decline in bearing lubrication performance and ultimately failure. Therefore, wear prediction of sliding bearings under heavy load conditions has become particularly important. Accurate wear prediction allows for timely understanding of the bearing's wear condition. On the one hand, it helps to identify potential problems early, avoiding unexpected downtime and increased maintenance costs; on the other hand, it provides a scientific basis for equipment maintenance and replacement.
[0003] Under heavy loads or start-stop conditions, wear easily occurs in the contact area between the journal and bearing shell of a sliding bearing. Currently, wear is mainly calculated using the Archard model. However, the Archard model is only applicable to steady-state wear conditions and is not suitable for wear conditions such as heavy loads, running-in under extreme environments, adhesive wear, and three-body wear. The accuracy of its results heavily depends on the wear coefficient input by the user. The accuracy of wear predictions is compromised; furthermore, the Archard model assumes that the material roughness is constant, but the roughness of actual materials may vary with time, wear distance, and load, making the wear prediction results unreliable. Summary of the Invention
[0004] The purpose of this invention is to provide an adaptive wear calculation method and system for heavy-duty sliding bearings. This invention iteratively calculates the bearing's wear process, gradually updating information such as wear depth, roughness, film thickness, pressure, and temperature until a preset stopping condition is reached. The final simulation results are output to evaluate the bearing's wear behavior and performance changes, as well as its wear life.
[0005] To achieve the above objectives, the present invention adopts the following technical solution: A method for calculating adaptive wear of heavy-duty sliding bearings includes the following steps: 1) Read the input parameters, including bearing geometry parameters, lubricating oil parameters, operating condition parameters, hardness coefficient, and current wear parameters; 2) Based on the input parameters, initialize the relevant variables, including the roughness of the journal and bearing, the wear depth array, and the time step array; 3) Calculate the combined roughness based on the initial roughness of the journal and bearing bush. ; 4) Based on oil film thickness h and combined roughness Determine whether to ignore subsequent wear and tear. h / If the value is ≤4, proceed to the next step; if not, jump to step 13). 5) Based on the initialized relevant variables, and following the idea that the oil film and roughness peaks jointly bear the load, the contact pressure is calculated using the Jane roughness contact model. ; 6) Based on contact pressure Adaptive calculation of maximum sliding distance L ; 7) Based on maximum sliding distance L Adaptive calculation of maximum time step ; 8) Based on the initialized relevant variables and contact pressure and maximum time step Update real-time wear coefficient K ; 9) Based on real-time wear coefficient K and maximum sliding distance L Calculate the current wear depth based on the Archard model. ; 10) Determine the shaft hardness coefficient H j With tile hardness coefficient H b If the relationship, H j ≈ H b If the shaft and bearing wear simultaneously, the wear depth array for both will be adjusted synchronously; if H j < H b Then adjust the journal wear depth array; if H j > H b Then adjust the bearing wear depth array; 11) Update the wear depth array and calculate the roughness change; 12) Determine if the bearing wear depth or maximum time step is greater than the planned depth or time step. If it is greater, proceed to the next step; if it is less, update the time. t , update film thickness h Then return to step 2); 13) Determine the final wear depth, film thickness, and contact pressure, and output the results.
[0006] A further improvement of this invention is that, in step 3), the formula for calculating the combined roughness is:
[0007] In the formula, and This refers to the surface roughness of the mating surfaces of the journal and the bearing bush.
[0008] A further improvement of this invention is that, in step 5), the contact pressure is calculated using Jane's rough contact model. ,include: According to Jane's rough contact model, the relationship between the dimensionless average clearance and the roughness parameters, material parameters, and dimensionless average pressure is determined by the following formula:
[0009]
[0010]
[0011]
[0012] In the formula, The gap is dimensionless average. L The dimensionless average contact pressure is defined as follows: ; The root mean square of the surface roughness; The rough peak autocorrelation coefficient is defined as follows: ,express , The ratio of directional feature lengths, Indicates vertical stripes, and vice versa; the default value is 1. The rough peak shape factor is defined as follows: This represents the ratio of the rough peak feature length to the feature height, with a default value of 40. For rough peak skewness, This indicates that the peak is lower than the valley, and vice versa; the default value is 0; K represents the rough peak kurtosis. This indicates a high peak density, and vice versa; the default value is 3. The characteristic length of the rough peak in the x-direction; The characteristic length of the rough peak in the y-direction is M; M is the relative Young's modulus of the material, defined as... ; E The Young's modulus of the material. The equivalent Young's modulus of the material; Y The yield strength of the material; Based on the dimensionless average gap and dimensionless average contact pressure L The relationship between them depends on the film thickness. h Interpolation to contact pressure .
[0013] A further improvement of the present invention is that, in step 6), based on contact pressure... Adaptive calculation of maximum sliding distance L ,include: Maximum sliding distance L：
[0014] In the formula, The oil film thickness over a specific time period. For maximum contact pressure, It is a stable term in the denominator.
[0015] A further improvement of the present invention is that, in step 7), based on the maximum sliding distance... L Adaptive calculation of maximum time step ,include: Time step :
[0016] In the formula, L The sliding distance, For the bearing angular velocity, Where is the bearing radius.
[0017] A further improvement of the present invention is that, in step 8), based on the initialized relevant variables and contact pressure... and maximum time step Update real-time wear coefficient K ,include: Calculate the real-time wear coefficient K : The initial calculation uses the recommended values obtained from the following experimental measurements, and subsequent iterations use the values calculated using the following formula; K Determined through material wear testing; Subsequent iteration calculation formula:
[0018]
[0019]
[0020]
[0021] In the formula, For wear depth,H It refers to the hardness of the softer side of the journal and bearing. To contact pressure, v For sliding speed, a and b These represent the steepness of the support ratio curve based on the roughness meter and the height of the midpoint of the curve, respectively. and These are the material ratios corresponding to the highest peak height and the lowest valley depth of the roughness profile, respectively. The total roughness height, For core roughness height, k It is a dimensionless wear coefficient.
[0022] A further improvement of the present invention is that, in step 9), based on the real-time wear coefficient... K and maximum sliding distance L Calculate the current wear depth based on the Archard model. ,include: Wear depth Calculation formula:
[0023] In the formula, For wear depth, For surface roughness, L The sliding distance, K This represents the real-time wear coefficient.
[0024] An adaptive wear calculation system for heavy-duty sliding bearings includes: The reading module reads input parameters, including bearing geometric parameters, lubricating oil parameters, operating condition parameters, hardness coefficient, and current wear parameters; The initialization module initializes relevant variables based on the input parameters, including the roughness of the journal and bearing, the wear depth array, and the time step array. The first calculation module calculates the combined roughness based on the initial roughness of the journal and bearing. ; The first judgment module determines the oil film thickness. h and combined roughness Determine whether to ignore subsequent wear and tear. h / If the value is ≤4, proceed to the next step; if not, jump to the output module. The second calculation module, based on the initialized relevant variables and following the idea that the oil film and roughness peaks jointly bear the load, calculates the contact pressure using the Jane rough contact model. ; The third calculation module is based on contact pressure. Adaptive calculation of maximum sliding distanceL ; The fourth calculation module is based on the maximum sliding distance. L Adaptive calculation of maximum time step ; The first update module is based on the initialized relevant variables and contact pressure. and maximum time step Update real-time wear coefficient K ; The fifth calculation module is based on the real-time wear coefficient. K and maximum sliding distance L Calculate the current wear depth based on the Archard model. ; The second judgment module determines the shaft hardness coefficient. H j With tile hardness coefficient H b If the relationship, H j ≈ H b If the shaft and bearing wear simultaneously, the wear depth array for both will be adjusted synchronously; if H j < H b Then adjust the journal wear depth array; if H j > H b Then adjust the bearing wear depth array; The second update module updates the wear depth array and calculates the roughness change. The third judgment module determines whether the bearing wear depth or the maximum time step is greater than the planned depth or time step. If it is greater, proceed to the next step; if it is less, update the time. t , update film thickness h Then return to the initialization module; The output module determines the final wear depth, film thickness, and contact pressure, and outputs the results.
[0025] A further improvement of this invention is that, in the first calculation module, the formula for calculating the combined roughness is:
[0026] In the formula, and This refers to the surface roughness of the mating surfaces of the journal and the bearing bush.
[0027] A computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of the adaptive wear calculation method for heavy-duty sliding bearings.
[0028] Compared with the prior art, the present invention has at least the following beneficial technical effects: 1) This invention provides an adaptive wear calculation method for heavy-duty sliding bearings, which can realize the evolution simulation of sliding bearing wear with information such as wear depth, roughness, film thickness, pressure and temperature, and calculate bearing performance and wear life.
[0029] 2) This invention derives the contact pressure of the sliding bearing using the Jane rough contact model, taking into account the elastoplastic deformation of the rough peaks, the mutual influence between different rough peaks, and the influence of statistical parameters of the rough peaks, such as root mean height, skewness, and kurtosis, on the results.
[0030] 3) This invention significantly improves the accuracy of wear process prediction by continuously updating the wear coefficient in real time. This makes the prediction results closer to actual application conditions, thus providing more reliable data support for equipment maintenance and life prediction. It can effectively prevent unexpected failures, optimize maintenance plans, extend equipment life, and ensure operational continuity and efficiency.
[0031] 4) The design method of the present invention has the advantages of simplicity and efficiency. It only requires continuous analysis and iteration of the nominal film thickness to determine the maximum time step in the wear process until the preset stopping condition is reached, which saves tedious calculation steps and is more efficient. Attached Figure Description
[0032] Figure 1 This is a flowchart of an adaptive wear calculation method for heavy-duty sliding bearings according to the present invention.
[0033] Figure 2 This is a schematic diagram of the roughness of the peak skewness.
[0034] Figure 3 This is a schematic diagram of the rough peak warp.
[0035] Figure 4 The surface profile and load-bearing ratio curves are shown.
[0036] Figure 5 The figure shows an example of the final implementation effect of the adaptive wear calculation method for heavy-duty sliding bearings of the present invention.
[0037] Figure 6 This is a structural block diagram of an adaptive wear calculation system for heavy-duty sliding bearings according to the present invention. Detailed Implementation
[0038] Exemplary embodiments of the present disclosure will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided to enable a more thorough understanding of the present disclosure and to fully convey the scope of the disclosure to those skilled in the art. It should be noted that, unless otherwise specified, the embodiments and features described herein can be combined with each other. The present invention will now be described in detail with reference to the accompanying drawings and embodiments.
[0039] Example 1 This invention provides a method for calculating adaptive wear of heavy-duty sliding bearings, comprising the following steps: 1) Read the input parameters, including bearing geometry parameters, lubricating oil parameters, operating condition parameters, hardness coefficient, and current wear parameters; 2) Based on the input parameters, initialize the relevant variables, including the roughness of the journal and bearing, the wear depth array, and the time step array; 3) Calculate the combined roughness based on the initial roughness of the journal and bearing bush. ; 4) Based on oil film thickness h and combined roughness Determine whether to ignore subsequent wear and tear. h / If the value is ≤4, proceed to the next step; if not, jump to step 13). 5) Based on the initialized relevant variables, and following the idea that the oil film and roughness peaks jointly bear the load, the contact pressure is calculated using the Jane roughness contact model. ; 6) Based on contact pressure Adaptive calculation of maximum sliding distance L ; 7) Based on maximum sliding distance L Adaptive calculation of maximum time step ; 8) Based on the initialized relevant variables and contact pressure and maximum time step Update real-time wear coefficient K ; 9) Based on real-time wear coefficient K and maximum sliding distance L Calculate the current wear depth based on the Archard model. ; 10) Determine the shaft hardness coefficient H j With tile hardness coefficient H b If the relationship,H j ≈ H b If the shaft and bearing wear simultaneously, the wear depth array for both will be adjusted synchronously; if H j < H b Then adjust the journal wear depth array; if H j > H b Then adjust the bearing wear depth array; 11) Update the wear depth array and calculate the roughness change; 12) Determine if the bearing wear depth or maximum time step is greater than the planned depth or time step. If it is greater, proceed to the next step; if it is less, update the time. t , update film thickness h Then return to step 2); 13) Determine the final wear depth, film thickness, and contact pressure, and output the results.
[0040] Example 2 A flowchart of an adaptive wear calculation method for heavy-duty sliding bearings according to an embodiment of the present invention includes: reading input parameters: bearing properties, initial state, and current wear parameters; initializing relevant variables: roughness, wear depth array, and time step array; calculating combined roughness; calculating contact pressure; adaptively calculating the maximum sliding distance; adaptively calculating the maximum time step; calculating the real-time wear coefficient; calculating the current wear depth according to the Archard model; determining whether the maximum wear depth or time step of the bearing is greater than the planned depth or time step; determining the final wear depth, film thickness, contact pressure, temperature, and roughness, and outputting the results. The design method of the present invention has the advantages of simplicity and efficiency, eliminating tedious calculation steps and making it more efficient.
[0041] like Figure 1 The diagram shows a flowchart of an adaptive wear calculation method for heavy-duty sliding bearings, including the following steps: The technical solution of the present invention is as follows: A method for calculating adaptive wear of heavy-duty sliding bearings, comprising the following steps: Step 1: Read the input parameters This includes bearing geometry parameters, lubricant parameters, operating condition parameters, hardness coefficient, and current wear parameters; Step 2: Initialize relevant variables Initialize relevant variables based on the input parameters, including roughness, wear depth array, time step array, etc. Step 3: Calculate the combined roughness Two rough surfaces The mean square value is the combined roughness:
[0042] In the formula, and This refers to the surface roughness of the mating surfaces of the journal and the bearing bush.
[0043] Step 4: Based on the oil film thickness h and combined roughness The ratio determines whether subsequent wear should be ignored. h / If the value is ≤4, then subsequent wear cannot be ignored, proceed to the next step; otherwise, skip to step thirteen.
[0044] Step 5: Based on the initialized variables, and following the idea that the oil film and roughness peaks jointly bear the load, calculate the contact pressure using the Jane roughness contact model. ; According to Jane's contact model, the relationship between the dimensionless average clearance and the roughness parameters, material parameters, and dimensionless average pressure can be determined by the following formula:
[0045]
[0046]
[0047]
[0048] In the formula, The gap is dimensionless average. L The dimensionless average contact pressure is defined as follows: ; The root mean square of the surface roughness; The rough peak autocorrelation coefficient is defined as follows: ,express , The ratio of directional feature lengths, Indicates vertical stripes, and vice versa; the default value is 1. The rough peak shape factor is defined as follows: This represents the ratio of the rough peak feature length to the feature height, with a default value of 40. Roughness peak skewness (see schematic diagram) Figure 1 (as shown) This indicates that the peak is lower than the trough, and vice versa; the default value is 0. K Rough peak kurtosis (see diagram) Figure 2 (as shown) This indicates a high peak density, and vice versa; the default value is 3. The characteristic length (autocorrelation length) of the rough peak in the x-direction; is the characteristic length (autocorrelation length) of the rough peak in the y-direction; M is the relative Young's modulus of the material, defined as... ; E The Young's modulus of the material. The equivalent Young's modulus of the material; Y The yield strength of the material.
[0049] Based on the dimensionless average gap and dimensionless average contact pressure L The relationship between them can be determined based on the film thickness. h Interpolation to contact pressure .
[0050] T For temperature, k The thermal conductivity coefficient, for y Thermal conductivity in the direction, This is the equivalent viscosity.
[0051] Step 6: Adaptively calculate the maximum sliding distance L。
[0052] Maximum sliding distance L：
[0053] In the formula, The oil film thickness over a specific time period. For maximum contact pressure, It is a stable term in the denominator.
[0054] Step 7: Adaptively calculate the maximum time step .
[0055] Time step :
[0056] In the formula, L The sliding distance, For the bearing angular velocity, Where is the bearing radius.
[0057] Step 8: Calculate the real-time wear coefficient K .
[0058] The initial calculations used the recommended values obtained from the following experimental measurements, while subsequent iterations used the values calculated using the following formulas. K The wear coefficient needs to be determined through material wear tests, as different materials have significantly different wear coefficients, ranging from 5 to 50,000. It is also related to operating conditions such as load and linear velocity.
[0059] Tin bronze can be taken from 5~10
[0060] Babbitt alloy can be used in the range of 500~2000.
[0061] Modified PEEK can be 200~400.
[0062] Subsequent iteration calculation formula:
[0063]
[0064]
[0065]
[0066] In the formula, For wear depth, H It refers to the hardness of the softer side of the journal and bearing. To contact pressure, v For sliding speed, a and b These represent the steepness of the support ratio curve based on the roughness meter and the height of the midpoint of the curve, respectively. and These are the material ratios corresponding to the highest peak height and the lowest valley depth of the roughness profile, respectively. The total roughness height, For the core roughness height, such as Figure 3 As shown, k It is a dimensionless wear coefficient.
[0067] Step 9: Calculate the current wear depth based on the Archard model. .
[0068]
[0069] In the formula, For wear depth, For surface roughness, L The sliding distance, K This represents the real-time wear coefficient.
[0070] Step 10: Determine the shaft stiffness coefficient H j With tile hardness coefficient H b If the relationship, H j ≈H b If the shaft and bearing wear simultaneously, the wear depth array for both will be adjusted synchronously; if H j < H b Then adjust the journal wear depth array; if H j > H b Then adjust the bearing wear depth array.
[0071] Step 11: Update the wear depth array and calculate the roughness change.
[0072] Step 12: Determine if the maximum wear depth or time step of the bearing is greater than the planned wear depth or time step. If it is greater, proceed to the next step; otherwise, update the time. t , update film thickness h Then return to step two.
[0073] Step 13: Determine the final wear depth, film thickness, contact pressure, temperature, and roughness, and output the results.
[0074] Figure 2 This is a schematic diagram of the roughness of the peak skewness. Figure 3 The diagram above illustrates the roughness peak curvature. The Jane contact model calculates the contact performance of rough surfaces with a large number of different parameter combinations. Through multi-parameter regression analysis, it provides a series of fitting formulas for parameters such as average gap, contact area, and plastic deformation rate as a function of contact pressure. The two diagrams above show the parameters required for calculation as follows: roughness peak skewness, roughness peak kurtosis, roughness peak autocorrelation coefficient, and roughness peak shape coefficient.
[0075] like Figure 5 The image shows an example of the final implementation effect of an adaptive wear calculation method for heavy-duty sliding bearings. This method can accurately calculate the film thickness, wear depth, and contact pressure of heavy-duty sliding bearings at different time points during actual operation. It can also clearly determine the maximum wear depth and highest contact pressure under the operating conditions, and predict the wear trend of sliding bearings. This provides strong data support for bearing design optimization, maintenance strategies, and life management.
[0076] This invention provides an adaptive wear calculation method for heavy-duty sliding bearings. It iteratively calculates the bearing's wear process, gradually updating information such as wear depth, roughness, film thickness, pressure, and temperature until a preset stopping condition is reached. The final simulation results are output to evaluate the bearing's wear behavior and performance changes.
[0077] Example 3 like Figure 6As shown, the present invention provides an adaptive wear calculation system for heavy-duty sliding bearings, comprising: The reading module reads input parameters, including bearing geometric parameters, lubricating oil parameters, operating condition parameters, hardness coefficient, and current wear parameters; The initialization module initializes relevant variables based on the input parameters, including the roughness of the journal and bearing, the wear depth array, and the time step array. The first calculation module calculates the combined roughness based on the initial roughness of the journal and bearing. ; The first judgment module determines the oil film thickness. h and combined roughness Determine whether to ignore subsequent wear and tear. h / If the value is ≤4, proceed to the next step; if not, jump to the output module. The second calculation module, based on the initialized relevant variables and following the idea that the oil film and roughness peaks jointly bear the load, calculates the contact pressure using the Jane rough contact model. ; The third calculation module is based on contact pressure. Adaptive calculation of maximum sliding distance L ; The fourth calculation module is based on the maximum sliding distance. L Adaptive calculation of maximum time step ; The first update module is based on the initialized relevant variables and contact pressure. and maximum time step Update real-time wear coefficient K ; The fifth calculation module is based on the real-time wear coefficient. K and maximum sliding distance L Calculate the current wear depth based on the Archard model. ; The second judgment module determines the shaft hardness coefficient. H j With tile hardness coefficient H b If the relationship, H j ≈ H b If the shaft and bearing wear simultaneously, the wear depth array for both will be adjusted synchronously; if H j < H b Then adjust the journal wear depth array; if H j > Hb Then adjust the bearing wear depth array; The second update module updates the wear depth array and calculates the roughness change. The third judgment module determines whether the bearing wear depth or the maximum time step is greater than the planned depth or time step. If it is greater, proceed to the next step; if it is less, update the time. t , update film thickness h Then return to the initialization module; The output module determines the final wear depth, film thickness, and contact pressure, and outputs the results.
[0078] Example 4 The present invention provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of the adaptive wear calculation method for heavy-duty sliding bearings.
[0079] Those skilled in the art will understand that embodiments of this application can be provided as methods, systems, or computer program products. Therefore, this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this application can take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.
[0080] This application is described with reference to flowchart illustrations and / or block diagrams of methods, systems, and computer program products according to embodiments of this application. It will be understood that each block of 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, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart... Figure 1 One or more processes and / or boxes Figure 1 A system that specifies functions in one or more boxes.
[0081] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.
[0082] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.
[0083] The foregoing has shown and described the basic principles, main features, and advantages of the present invention. It will be apparent to those skilled in the art that the invention is not limited to the details of the exemplary embodiments described above, and that the invention can be implemented in other specific forms without departing from its spirit or essential characteristics. Therefore, the embodiments should be considered illustrative and non-limiting in all respects, and the scope of the invention is defined by the appended claims rather than the foregoing description. Thus, all variations falling within the meaning and scope of equivalents of the claims are intended to be included within the scope of the invention. No reference numerals in the claims should be construed as limiting the scope of the claims.
[0084] Furthermore, it should be understood that although this specification describes embodiments, not every embodiment contains only one independent technical solution. This narrative style is merely for clarity. Those skilled in the art should consider the specification as a whole, and the technical solutions in each embodiment can be appropriately combined to form other embodiments that can be understood by those skilled in the art. The above content is only for illustrating the technical concept of the present invention and should not be construed as limiting the scope of protection of the present invention. Any modifications made based on the technical concept proposed in this invention shall fall within the scope of protection of the claims of this invention.
Claims
1. A method for calculating adaptive wear of heavy-duty sliding bearings, characterized in that, Includes the following steps: 1) Read the input parameters, including bearing geometry parameters, lubricating oil parameters, operating condition parameters, hardness coefficient, and current wear parameters; 2) Based on the input parameters, initialize the relevant variables, including the roughness of the journal and bearing, the wear depth array, and the time step array; 3) Calculate the combined roughness based on the initial roughness of the journal and bearing bush. ; 4) Based on oil film thickness h and combined roughness Determine whether to ignore subsequent wear and tear. h / If the value is ≤4, proceed to the next step; if not, jump to step 13). 5) Based on the initialized relevant variables, and following the idea that the oil film and roughness peaks jointly bear the load, the contact pressure is calculated using the Jane roughness contact model. ; 6) Based on contact pressure Adaptive calculation of maximum sliding distance L ; 7) Based on maximum sliding distance L Adaptive calculation of maximum time step ; 8) Based on the initialized relevant variables and contact pressure and maximum time step Update real-time wear coefficient K ; 9) Based on real-time wear coefficient K and maximum sliding distance L Calculate the current wear depth based on the Archard model. ; 10) Determine the shaft hardness coefficient H j With tile hardness coefficient H b If the relationship, H j ≈ H b If the shaft and bearing wear simultaneously, the wear depth array for both will be adjusted synchronously; if H j < H b Then adjust the journal wear depth array; if H j > H b Then adjust the bearing wear depth array; 11) Update the wear depth array and calculate the roughness change; 12) Determine if the bearing wear depth or maximum time step is greater than the planned depth or time step. If it is greater, proceed to the next step; if it is less, update the time. t , update film thickness h Then return to step 2); 13) Determine the final wear depth, film thickness, and contact pressure, and output the results.
2. The adaptive wear calculation method for heavy-duty sliding bearings according to claim 1, characterized in that, In step 3), the formula for calculating the combined roughness is: In the formula, and This refers to the surface roughness of the mating surfaces of the journal and the bearing bush.
3. The adaptive wear calculation method for heavy-duty sliding bearings according to claim 1, characterized in that, In step 5), the contact pressure is calculated using Jane's rough contact model. ,include: According to Jane's rough contact model, the relationship between the dimensionless average clearance and the roughness parameters, material parameters, and dimensionless average pressure is determined by the following formula: In the formula, The gap is dimensionless average. L The dimensionless average contact pressure is defined as follows: ; The root mean square of the surface roughness; The rough peak autocorrelation coefficient is defined as follows: ,express , The ratio of directional feature lengths, Indicates vertical stripes, and vice versa; the default value is 1. The rough peak shape factor is defined as follows: This represents the ratio of the rough peak feature length to the feature height, with a default value of 40. For rough peak skewness, This indicates that the peak is lower than the valley, and vice versa; the default value is 0; K represents the rough peak kurtosis. This indicates a high peak density, and vice versa; the default value is 3. The characteristic length of the rough peak in the x-direction; The characteristic length of the rough peak in the y-direction is M; M is the relative Young's modulus of the material, defined as... ; E The Young's modulus of the material. The equivalent Young's modulus of the material; Y The yield strength of the material; Based on the dimensionless average gap and dimensionless average contact pressure L The relationship between them depends on the film thickness. h Interpolation to contact pressure .
4. The adaptive wear calculation method for heavy-duty sliding bearings according to claim 1, characterized in that, In step 6), based on contact pressure Adaptive calculation of maximum sliding distance L ,include: Maximum sliding distance L： In the formula, The oil film thickness over a specific time period. For maximum contact pressure, It is a stable term in the denominator.
5. The adaptive wear calculation method for heavy-duty sliding bearings according to claim 1, characterized in that, In step 7), based on the maximum sliding distance L Adaptive calculation of maximum time step ,include: Time step : In the formula, L The sliding distance, For the bearing angular velocity, Where is the bearing radius.
6. The adaptive wear calculation method for heavy-duty sliding bearings according to claim 1, characterized in that, In step 8), based on the initialized relevant variables and contact pressure... and maximum time step Update real-time wear coefficient K ,include: Calculate the real-time wear coefficient K : The initial calculation uses the recommended values obtained from the following experimental measurements, and subsequent iterations use the values calculated using the following formula; K Determined through material wear testing; Subsequent iteration calculation formula: In the formula, For wear depth, H It refers to the hardness of the softer side of the journal and bearing. To contact pressure, v For sliding speed, a and b These represent the steepness of the support ratio curve based on the roughness meter and the height of the midpoint of the curve, respectively. and These are the material ratios corresponding to the highest peak height and the lowest valley depth of the roughness profile, respectively. The total roughness height, For core roughness height, k It is a dimensionless wear coefficient.
7. The adaptive wear calculation method for heavy-duty sliding bearings according to claim 1, characterized in that, In step 9), based on the real-time wear coefficient K and maximum sliding distance L Calculate the current wear depth based on the Archard model. ,include: Wear depth Calculation formula: In the formula, For wear depth, For surface roughness, L The sliding distance, K This represents the real-time wear coefficient.
8. A heavy-duty sliding bearing adaptive wear calculation system, characterized in that, include: The reading module reads input parameters, including bearing geometric parameters, lubricating oil parameters, operating condition parameters, hardness coefficient, and current wear parameters; The initialization module initializes relevant variables based on the input parameters, including the roughness of the journal and bearing, the wear depth array, and the time step array. The first calculation module calculates the combined roughness based on the initial roughness of the journal and bearing. ; The first judgment module, based on the oil film thickness... h and combined roughness Determine whether to ignore subsequent wear and tear. h / If the value is ≤4, proceed to the next step; if not, jump to the output module. The second calculation module, based on the initialized relevant variables and following the idea that the oil film and roughness peaks jointly bear the load, calculates the contact pressure using the Jane rough contact model. ; The third calculation module is based on contact pressure. Adaptive calculation of maximum sliding distance L ; The fourth calculation module is based on the maximum sliding distance. L Adaptive calculation of maximum time step ; The first update module is based on the initialized relevant variables and contact pressure. and maximum time step Update real-time wear coefficient K ; The fifth calculation module is based on the real-time wear coefficient. K and maximum sliding distance L Calculate the current wear depth based on the Archard model. ; The second judgment module determines the shaft hardness coefficient. H j With tile hardness coefficient H b If the relationship, H j ≈ H b If the shaft and bearing wear simultaneously, the wear depth array for both will be adjusted synchronously; if H j < H b Then adjust the journal wear depth array; if H j > H b Then adjust the bearing wear depth array; The second update module updates the wear depth array and calculates the roughness change. The third judgment module determines whether the bearing wear depth or the maximum time step is greater than the planned depth or time step. If it is greater, proceed to the next step; if it is less, update the time. t , update film thickness h Then return to the initialization module; The output module determines the final wear depth, film thickness, and contact pressure, and outputs the results.
9. The adaptive wear calculation system for heavy-duty sliding bearings according to claim 8, characterized in that, In the first calculation module, the formula for calculating combined roughness is: In the formula, and This refers to the surface roughness of the mating surfaces of the journal and the bearing bush.
10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements the steps of the adaptive wear calculation method for heavy-duty sliding bearings according to any one of claims 1-7.