Method for evaluating axial load distribution and wear life of connecting part

Abstract

The invention relates to the technical field of connecting part load distribution and life analysis, in particular to a method for evaluating axial load distribution and wear life of a connecting part, which comprises the following steps of: generating a three-dimensional geometric model by inputting key structure parameters of a spline coupling so as to identify a tooth surface contact area, a transition area and a non-key area; the method comprises the following steps of: constructing a combination assembly model through rigid body transformation, identifying a physical contact area to generate a node pair set, and identifying a potential stress concentration and contact instability area after applying a load working condition in the assembly model. The method comprises the following steps of: selecting a local mesh as a candidate region of local mesh refinement, integrating to generate a complete finite element input file, calculating the wear depth of a contact node based on a wear model, and determining whether to terminate simulation, so that the whole process automation of the connecting part from modeling, assembling to wear prediction is realized, and the analysis efficiency and the result accuracy are remarkably improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of connection component load distribution and life analysis, in particular to a method for evaluating axial load distribution and wear life of a connection component. BACKGROUND

[0002] With the continuous progress of high-end equipment manufacturing technology, typical connection components such as splines and gears have been widely used in various mechanical transmission systems. The axial load distribution, assembly registration accuracy, tooth surface contact state and wear evolution process of the connection components have important influence on the stability and life of the whole machine;

[0003] However, the existing finite element analysis method has many limitations in actual engineering. It usually requires engineering personnel to manually complete steps such as three-dimensional modeling, contact definition and loading setting. The operation process is tedious and time-consuming, and it is difficult to achieve high-precision modeling and truly reflect the mechanical behavior when dealing with complex structural features such as assembly eccentricity, tooth profile gradual change and transition chamfer. This analysis method relying on manual setting is inefficient when facing variable working conditions and batch evaluation requirements, and has been difficult to meet the refinement and automation requirements of connection component performance analysis of complex equipment.

[0004] In order to solve the above defects, a technical scheme is provided. SUMMARY

[0005] In view of the above-mentioned shortcomings of the prior art, the present application provides a method for evaluating the axial load distribution and wear life of a connection component, which can effectively solve the problem that the complex geometric features such as assembly eccentricity, tooth profile gradual change and transition chamfer are difficult to model accurately, and the analysis process relies on manual setting, is low in efficiency and poor in result consistency.

[0006] To achieve the above purpose, the present application can be realized by the following technical scheme:

[0007] The present application provides a method for evaluating the axial load distribution and wear life of a connection component, comprising the following steps:

[0008] Through the preset Python script interface, the user inputs the key structural parameters of the spline coupling;

[0009] After receiving the key structural parameters input by the user, the three-dimensional modeling process is automatically executed to construct a three-dimensional geometric model;

[0010] Based on the three-dimensional geometric model, the tooth surface contact area, transition area and non-key area in the model are automatically identified, and a zoning and grading mesh division strategy is adopted according to the characteristics of each area;

[0011] Based on the set assembly parameters, perform geometric registration between assemblies, apply rigid body transformation to the inner spline body using the coordinate transformation matrix, complete the registration, and generate the assembly model;

[0012] In the assembly model, identify the area where physical contact occurs, and construct a node pair set;

[0013] Apply typical working condition load to the assembly model, identify potential stress concentration areas and contact instability areas as candidate areas for local mesh refinement;

[0014] Integrate three-dimensional geometric model, structured grid, assembly model, node pair set and candidate area into complete finite element input file;

[0015] Based on the Archard wear model, the wear depth of the tooth surface node is obtained, and if the wear depth of any point on the tooth surface is greater than the preset wear depth threshold, the simulation is terminated, and the life prediction result is output.

[0016] Further, construct a three-dimensional geometric model, the specific process is as follows:

[0017] According to the standard involute tooth profile formula:

[0018]

[0019] Where m represents the modulus, a represents the pressure angle, and theta represents the polar angle parameter;

[0020] Generate polar coordinate point set of tooth profile by traversing angle interval, and convert it to two-dimensional point array in Cartesian coordinate system;

[0021] Subsequently, combined with the geometric features of involute, circular arc transition section and chamfer area, through multi-segment curve splicing and high-density sampling fitting, a complete tooth profile is constructed, which meets the accuracy requirements of finite element analysis;

[0022] Then use the modeling API to automatically draw the two-dimensional tooth profile section according to the single tooth profile, and according to the set number of teeth Z, copy the equiangular array to form a complete internal / external spline profile structure;

[0023] Based on the above two-dimensional profile, further perform stretching operation to construct three-dimensional geometric model, and according to the user set axial length, step structure, keyway parameter to supplement the shaft class detail structure, and according to the tolerance requirement to adjust the geometric size.

[0024] Further, a zoning and grading meshing strategy is adopted according to the characteristics of each region, the specific process is as follows:

[0025] For tooth surface contact area, use high-precision eight-node or twenty-node three-dimensional solid element;

[0026] For the transition region, medium-density ten-node tetrahedral elements are used;

[0027] For non-critical areas, linear elements with reduced degrees of freedom are used;

[0028] Automatically calculate the theoretical contact radius r of the tooth surface based on the typical Hertz contact theory, which serves as a reference scale for mesh refinement;

[0029] The contact radius calculation formula is: Among them, T in It is represented by the currently applied input torque, R is represented by the equivalent curvature radius of the tooth surface (obtained by synthesizing the curvature of the internal and external spline tooth profiles), and E * Expressed as the equivalent elastic modulus of the material;

[0030] The grid density is arranged in levels according to the contact radius r, and the specific arrangement is:

[0031] Divide the radius into a dense grid with a size less than r / 3;

[0032] Set a gradient transition grid in the range of r to 3r, gradually transitioning to the background coarse grid;

[0033] Set a sparse grid larger than 3r for the remaining areas;

[0034] The entire meshing process also includes a quality control mechanism to ensure that the unit distortion is controlled within 0.8 and the maximum unit size in the contact area does not exceed the minimum contact width. If the local stress gradient exceeds a preset threshold, the mesh refinement or reconstruction process in that area will be automatically triggered.

[0035] Furthermore, the assembly model of the assembly is generated. The specific process is as follows:

[0036] The transformation matrix is ​​expressed as: T = R(θ0)·T(δe);

[0037] in,

[0038] The offset δe refers to the offset of the internal spline relative to the external spline, and the initial rotation angle θ0 refers to the initial rotation angle of the internal spline around its own axis. e ) refers to the translation matrix, and R(θ0) refers to the rotation matrix around the z-axis;

[0039] Then, the bounding box detection algorithm is used to automatically detect the interference after registration. If there are penetration or overlapping areas, the system will prompt you to adjust the parameters. If the assembly is successful, the assembly model of the combination will be generated.

[0040] Furthermore, a node pair set is constructed. The specific process is as follows:

[0041] Let Γ in be the inner tooth surface and Γ out be the outer tooth surface. in For any point X ∈ Γ out , define the minimum distance from X to the nearest point Y ∈ Γ

[0042] Set the contact tolerance ε, define the points that satisfy d min (X) ≤ ε as the potential contact area, and construct the contact area set Ω c = {X ∈ Γ in | d min (X) ≤ ε}.

[0043] For each X ∈ Ω c , extract its normal vector n in (X) and the normal vector n out (Y) of the corresponding nearest point Y, and construct the contact direction relationship.

[0044] In the identified contact area Ω c , extract the most adjacent node pairs on the inner and outer tooth surfaces, and establish a one-to-one mapping relationship, which is as follows:

[0045] If X i ∈ Γ in , Y j ∈ Γ out , and satisfy , then (X i , Y j ) is taken as a pair of contact node pairs, and each node pair is assigned a unique identifier, where the indices i, j, and k represent the inner tooth number, the outer tooth number, and the node sequence number on the tooth surface, respectively.

[0046] Sort the node pairs in the geometric order of the tooth surface contact path, and construct the node pair set.

[0047] Further, identify the potential stress concentration area and the contact instability area, and the specific process is as follows:

[0048] Let the contact pressure corresponding to the node X i be p i , and its rate of change over time be

[0049] The stress gradient of the node position is

[0050] The pre-set stress gradient judgment threshold is , and the contact pressure change rate threshold is ε p .

[0051] Then define the candidate area A adaptA node set satisfying any of the following conditions: Wherein, Indicates the gradient size of the current node stress in space, which is used to capture the area of stress change, Indicates the trend of pressure change per unit time during contact, which is used to identify the dynamic contact instability area.

[0052] Further, the candidate area is subjected to local mesh refinement, and the specific process is as follows:

[0053] For the area determined as A adapt , automatic mesh local refinement is performed, and the operation includes: local mesh refinement reconstruction, contact patch mesh alignment and error-driven mesh updating mechanism.

[0054] Further, the life prediction result is output as: the axial load distribution of the tooth surface, the target wear value, the predicted life, and the fitting parameter, which can be saved as a txt report.

[0055] The technical scheme provided by the application has the following beneficial effects compared with the known prior art:

[0056] 1、The application realizes full-process automation from structure parameter input to three-dimensional modeling, mesh division, assembly registration and finite element solution file generation through the preset Python script interface, without manual intervention throughout, greatly improving modeling efficiency and consistency;

[0057] 2、The application automatically identifies and processes complex structural features such as assembly eccentricity, tooth profile gradual change and transition chamfer, adapts to registration and load application under various actual working conditions, and significantly improves the adaptability to complex engineering scenes;

[0058] 3、The application proposes a multi-scale partitioning and grading mesh division strategy based on contact mechanics theory, combines Hertz contact radius and stress gradient distribution, realizes high-density division in key areas, and balances overall solving efficiency while ensuring calculation accuracy;

[0059] 4、The application automatically constructs the node mapping relationship between the inner and outer tooth surfaces through the nearest neighbor search and contact tolerance control mechanism, accurately restores the physical contact boundary, and lays a stable foundation for subsequent wear prediction and contact analysis;

[0060] 5、The application realizes intelligent identification and local mesh adaptive reconstruction of the contact instability area through the dual threshold mechanism of stress gradient and contact pressure change rate, significantly enhances the accuracy and stability of the simulation result. BRIEF DESCRIPTION OF DRAWINGS

[0061] In order to make the technical solutions in the embodiments of the present application or the prior art clearer, the accompanying drawings needed in the embodiments or prior art description will be briefly introduced. Obviously, the accompanying drawings in the following description only aim to some of the embodiments of the present application, and for those skilled in the art, other drawings can be obtained without creative effort on the basis of these drawings.

[0062] Figure 1 The overall flowchart of the present application. DETAILED DESCRIPTION

[0063] In order to make the technical solutions in the embodiments of the present application or the prior art clearer, the accompanying drawings needed in the embodiments or prior art description will be briefly introduced. Obviously, the accompanying drawings in the following description only aim to some of the embodiments of the present application, and for those skilled in the art, other drawings can be obtained without creative effort on the basis of these drawings.

[0064] As shown in Figure 1 A method for evaluating the axial load distribution and wear life of a connecting component, comprising:

[0065] Step one, parameter input: through a preset Python script interface, a user inputs key structural parameters of a spline coupling, the input key structural parameters include but are not limited to tooth number, module, pressure angle, chamfer radius, diameter of each circle, tolerance, assembly eccentricity / rotation amount, material parameter, torque, friction coefficient, wear parameter, the key structural parameters constitute a complete set of modeling and simulation input, which serves as a basis for subsequent geometric modeling, assembly definition and finite element analysis.

[0066] Step two, three-dimensional geometric modeling: after receiving the key structural parameters input by the user, a three-dimensional modeling process is automatically executed to complete the whole process from tooth profile curve generation to assembly geometry construction, the specific implementation steps are as follows:

[0067] According to the standard involute tooth profile formula:

[0068]

[0069] Wherein, m represents the module, a represents the pressure angle, and θ represents the polar angle parameter.

[0070] The polar coordinate point set of the tooth profile is generated by traversing the angle interval, and is converted into a two-dimensional point array in the Cartesian coordinate system.

[0071] Subsequently, combined with the geometric characteristics of involute, circular arc transition section and chamfer area, the complete tooth profile contour is constructed by multi-segment curve splicing and high-density sampling fitting, which meets the accuracy requirements of finite element analysis.

[0072] Then, the two-dimensional tooth profile section is automatically drawn according to the single tooth profile by using modeling API (such as Abaqus / Python or SolidWorks API), and the complete internal / external spline contour structure is formed by equiangular array replication according to the set number of teeth Z.

[0073] Based on the above two-dimensional contour, further stretching operation is performed to construct a three-dimensional geometric model, and the shaft details such as axial length, step structure and keyway parameters are supplemented according to user settings, and the geometric size is fine-tuned according to the tolerance requirements to ensure the assembly accuracy.

[0074] After completing the construction of the three-dimensional geometric model, it is exported to a standard neutral geometric format (such as.step or.sat), and automatically completes the geometric topology check (such as face closure, body consistency) to ensure that the model can be directly imported into the finite element software (such as Abaqus) and used for meshing.

[0075] Step three, automatic structured meshing: based on the three-dimensional geometric model, the tooth surface contact area, transition area and non-critical area in the model are automatically identified, and a zoning and grading meshing strategy is adopted according to the characteristics of each area.

[0076] Among them, the tooth surface contact area mainly refers to the contact tooth profile of the internal and external spline meshing surface;

[0077] The transition area mainly refers to the parts with relatively large stress gradient such as hub connection and shaft sleeve root;

[0078] The non-critical area mainly refers to the geometric area far away from the contact area of the rim;

[0079] According to different areas, appropriate mesh element types are automatically matched, which are as follows:

[0080] For the tooth surface contact area, high-precision eight-node or twenty-node three-dimensional solid elements are used to enhance the contact accuracy;

[0081] For the transition area, medium-density ten-node tetrahedral elements are used to balance accuracy and efficiency;

[0082] For the non-critical area, linear elements with reduced degrees of freedom are used to compress the calculation resources;

[0083] In order to realize the automatic optimization of contact accuracy and local stress capture effect, the theoretical contact radius r of the tooth surface is automatically calculated based on the typical Hertz contact theory, which is used as the reference scale for mesh refinement.

[0084] The contact radius calculation formula is: wherein T in represents the current applied input torque, R represents the equivalent curvature radius of the tooth surface (obtained by synthesizing the curvatures of the inner and outer spline tooth profiles), E * represents the equivalent elastic modulus of the material, and is calculated as follows: wherein E1, E2 and v1, v2 are the elastic modulus and Poisson's ratio of the materials on both sides of the contact;

[0085] The grid density is arranged according to the contact radius r, and the specific arrangement is as follows:

[0086] Within the radius r, dense grids with a size less than r / 3 are divided;

[0087] In the range of r to 3r, a gradually transitioned grid is set, which gradually transitions to a background coarse grid;

[0088] For the remaining area, sparse grids with a size greater than 3r are set;

[0089] The entire grid division process also includes a quality control mechanism to ensure that the unit twist is controlled within 0.8, the maximum unit size in the contact area does not exceed the minimum contact width, and if the local stress gradient exceeds a preset threshold (such as 5% / mm), the grid refinement or reconstruction process in that area will be automatically triggered;

[0090] In addition, for high curvature or geometric mutation sites, high-order units with surface fitting are preferentially used for approximation modeling to ensure modeling accuracy and convergence stability.

[0091] Step four, automatic assembly and coordinate transformation: based on the set assembly parameters, such as offset δ e and initial rotation angle θ0, the geometric registration between the assembled bodies is performed, and a rigid body transformation is applied to the inner spline body using the coordinate transformation matrix to complete the registration, realize the precise assembly alignment between the inner and outer splines, and ensure that the model meets the definition requirements of the mechanical boundary conditions in the initial state;

[0092] The transformation matrix is expressed as: T = R(θ0) · T(δe);

[0093] wherein,

[0094] The offset δe refers to the offset of the inner spline relative to the outer spline, the initial rotation angle θ0 refers to the initial rotation angle of the inner spline around its own axis, and T(δ e ) refers to the translation matrix, and R(θ0) refers to the rotation matrix around the z-axis;

[0095] The bounding box detection algorithm is used to automatically detect the interference after registration. If there is penetration or overlapping area, the parameters will be adjusted. If the assembly is successful, the assembly model of the combination is generated.

[0096] Step five, automatic identification of contact surface and tooth surface nodes: In the assembly model of the combination, the area where physical contact occurs is identified, and a set of node pairs is constructed. The specific steps are as follows:

[0097] Let the inner tooth surface be Γ in and the outer tooth surface be Γ out . For any point X∈Γ in , define the minimum distance from it to the nearest point Y∈Γ out on the outer tooth surface as:

[0098] Set the contact tolerance ε, and define the points that satisfy d min (X)≤ε as the potential contact area, and construct the contact area set: Ω c ={X∈Γ in |d min (X)≤ε};

[0099] Where, for each X∈Ω c , extract its normal vector n in (X) and the normal vector n out (Y) of the corresponding nearest point Y, and construct the contact direction relationship.

[0100] In the identified contact area Ω c , the nearest node pairs on the inner and outer tooth surfaces are extracted, and a one-to-one mapping relationship is established, which is as follows:

[0101] If X i ∈Γ in , Y j ∈Γ out , and satisfy , then (X i , Y j ) is taken as a pair of contact nodes, and each node pair is assigned a unique identifier, where the indices i, j, k represent the inner tooth number, the outer tooth number, and the node sequence number on the tooth surface, respectively.

[0102] The node pairs are sorted in the geometric order of the tooth surface contact path, and a set of node pairs is constructed, which serves as the input basis for subsequent finite element contact analysis and wear iteration.

[0103] Step six, adaptive mesh region identification: Apply typical working condition loads (such as axial force, torque, etc.) to the assembly model of the combination, identify potential stress concentration areas and contact instability areas, and use them as candidate areas for local mesh refinement. The specific identification process is as follows:

[0104] Set node X i The corresponding contact pressure is p i The rate of change of which over time is

[0105] The stress gradient of the node position is

[0106] The preset stress gradient determination threshold is The contact pressure change rate threshold is ε p ;

[0107] Then define the candidate area A adapt is a node set that meets any of the following conditions: Wherein, represents the gradient size of the current node stress in space, which is used to capture the area of rapid stress change, characterizes the rapid change trend of pressure per unit time during contact, which is used to identify the dynamic contact unstable area;

[0108] For the area determined as A adapt , automatically perform local mesh refinement, including:

[0109] Local mesh refinement reconstruction: introduce a size control function h(X), adaptively adjust the unit size according to the spatial position X of the target area, generate a denser grid structure, and improve the simulation accuracy of the area;

[0110] Contact patch element mesh alignment: for the master and slave contact surface areas, adopt the normal fitting strategy (i.e. rearrange the nodes along the contact normal), realize the coordinated alignment of the grid nodes of the two contact surfaces, and guarantee the stability and accuracy of the contact pressure and wear mechanics calculation;

[0111] Error-driven mesh updating mechanism: introduce an error estimation function e(X), and update the real-time error of A adapt according to the local error estimation control the next iteration step and the area, so as to dynamically respond to the evolution of physical field variables.

[0112] Step seven, automatic generation of finite element input file: integrate three-dimensional geometric model, structured grid, assembly model, node pair set and candidate area, etc. Into a complete finite element input file (such as Abaqus inp file), so as to support automatic iteration operation and wear update, without manual operation in CAE software;

[0113] Step eight, wear life analysis and termination criterion control: based on the Archard wear model, the wear depth of the tooth surface node is obtained, and if the wear depth of any point on the tooth surface is greater than the preset wear depth threshold, the simulation is terminated, and the life prediction result is output;

[0114] The output life prediction result is a tooth surface axial load distribution, a target wear value, a predicted life, a fitting parameter, etc., and can be saved as a txt report.

[0115] A computer device comprises a memory and a processor, the memory stores a computer program, and the processor implements the method of any one of the above when executing the computer program.

[0116] A computer readable storage medium stores a computer program, and the computer program is executed by a processor to implement the method of any one of the above.

[0117] The preferred embodiments of the application disclosed above are only used to help explain the application. The preferred embodiments do not describe all the details and limit the application to the specific embodiments. Obviously, many modifications and variations can be made according to the content of the present application. The present application selects and specifically describes these embodiments in order to better explain the principles and practical applications of the application, so that those skilled in the art can well understand and use the application. The application is limited by the claims and their full scope and equivalents.

Claims

1. A method for evaluating the axial load distribution and wear life of a connection component, characterized in that: The following steps are involved: Through the preset Python script interface, users input the key structural parameters of the spline coupling; After receiving the key structural parameters input by the user, the 3D modeling process is automatically executed to construct a 3D geometric model; Based on the 3D geometric model, the tooth contact area, transition area and non-critical area in the model are automatically identified, and a partitioned and graded meshing strategy is adopted according to the characteristics of each area; Based on the set assembly parameters, perform geometric registration between assemblies, apply rigid body transformation to the internal spline body using the coordinate transformation matrix, complete the registration, and generate the assembly model of the combination; In the assembly model of the composite body, the areas where physical contact occurs are identified and a set of node pairs is constructed; Apply typical working condition loads to the assembly model and identify potential stress concentration areas and contact instability areas as candidate areas for local mesh refinement. Integrate 3D geometric models, structured meshes, assembly models, node pair sets, and candidate regions into a complete finite element input file; The wear depth of the tooth surface nodes is obtained based on the Archard wear model. If the wear depth of any point on the tooth surface is greater than the preset wear depth threshold, the simulation is terminated and the life prediction result is output.

2. A method for evaluating the axial load distribution and wear life of a connection component according to claim 1, characterized in that: Construct a 3D geometric model. The specific process is as follows: According to the standard involute tooth profile formula: Among them, m represents the modulus, α represents the pressure angle, and θ represents the polar angle parameter; Generate the polar coordinate point set of the tooth profile by traversing the angle interval and convert it into a two-dimensional point matrix in the Cartesian coordinate system; Then, by combining geometric features such as involutes, arc transition sections, and chamfered areas, a complete tooth profile with the accuracy required for finite element analysis was constructed through multi-segment curve splicing and high-density sampling and fitting. Then, using the modeling API, a two-dimensional tooth profile is automatically drawn based on the single tooth profile, and an equiangular array is replicated according to the set number of teeth Z to form a complete internal / external spline profile structure; Based on the above two-dimensional contour, a stretching operation is further performed to construct a three-dimensional geometric model. The shaft details are supplemented according to the axial length, step structure, and keyway parameters set by the user, and the geometric dimensions are fine-tuned according to the tolerance requirements.

3. The method for evaluating the axial load distribution and wear life of a connection component according to claim 1, characterized in that: A zoning and hierarchical grid division strategy is adopted based on the characteristics of each region. The specific process is as follows: For the tooth surface contact area, high-precision eight-node or twenty-node three-dimensional solid elements are used; For the transition region, medium-density ten-node tetrahedral elements are used; For non-critical areas, linear elements with reduced degrees of freedom are used; Automatically calculate the theoretical contact radius r of the tooth surface based on the typical Hertz contact theory, which serves as a reference scale for mesh refinement; The contact radius calculation formula is: Among them, T in It is represented by the currently applied input torque, R is represented by the equivalent curvature radius of the tooth surface (obtained by synthesizing the curvature of the internal and external spline tooth profiles), and E * Expressed as the equivalent elastic modulus of the material; The grid density is arranged in levels according to the contact radius r, and the specific arrangement is: Divide the radius into a dense grid with a size less than r / 3; Set a gradient transition grid in the range of r to 3r, gradually transitioning to the background coarse grid; Set a sparse grid larger than 3r for the remaining areas; The entire meshing process also includes a quality control mechanism to ensure that the unit distortion is controlled within 0.8 and the maximum unit size in the contact area does not exceed the minimum contact width. If the local stress gradient exceeds a preset threshold, the mesh refinement or reconstruction process in that area will be automatically triggered.

4. The method for evaluating the axial load distribution and wear life of a connection component according to claim 1, characterized in that: Generate an assembly model of the assembly. The specific process is as follows: The transformation matrix is ​​expressed as: T = R(θ0)·T(δe); in, The offset δe refers to the offset of the internal spline relative to the external spline, and the initial rotation angle θ0 refers to the initial rotation angle of the internal spline around its own axis. e ) refers to the translation matrix, and R(θ0) refers to the rotation matrix around the z-axis; Then, the bounding box detection algorithm is used to automatically detect the interference after registration. If there are penetration or overlapping areas, the system will prompt you to adjust the parameters. If the assembly is successful, the assembly model of the combination will be generated.

5. The method for evaluating the axial load distribution and wear life of a connection component according to claim 1, characterized in that: Construct a node pair set. The specific process is as follows: Assume that the inner tooth surface is Γ in , the outer tooth surface is Γ out , for any point X∈Γ in , define its closest point to the outer tooth surface Y∈Γ out The minimum distance is: Set the contact tolerance ε and define the contact tolerance d min The points where (X)≤ε are potential contact areas, forming the contact area set: Ω c ={X∈Γ in |d min (X)≤ε}; Among them, for each X∈Ω c , extract its normal vector n in (X) and the normal vector n corresponding to the nearest point Y out (Y), construct contact direction relationship; In the identified contact area Ω c In the , the nearest node pairs on the inner and outer tooth surfaces are extracted and a one-to-one mapping relationship is established, specifically: If X i ∈Γ in ,Y j ∈Γ out , and satisfies Then (X i , Y j ) as a pair of contact node pairs, and assign a unique identifier to each node pair, where the indexes i, j, and k represent the internal tooth number, the external tooth number, and the node number on the tooth surface, respectively; The node pairs are sorted according to the geometric order of the tooth surface contact path and a node pair set is constructed.

6. The method for evaluating the axial load distribution and wear life of a connection component according to claim 1, characterized in that: Identify potential stress concentration areas and contact instability areas. The specific process is as follows: Set node X i The corresponding contact pressure is p i , whose rate of change over time is The stress gradient at the node is The preset stress gradient judgment threshold is The contact pressure change rate threshold is ε p ; Then define the candidate area A adapt A set of nodes that meet any of the following conditions: in, Indicates the gradient of the current node stress in space, which is used to capture the area where stress changes dramatically. Characterizes the dramatic change trend of pressure per unit time during the contact process and is used to identify the unstable area of ​​dynamic contact.

7. The method for evaluating the axial load distribution and wear life of a connection component according to claim 6, characterized in that: Perform local mesh refinement on the candidate area. The specific process is as follows: The pair is judged as A adapt In the area of ​​​​the mesh, local mesh refinement is automatically performed, and the operations include: local mesh refinement and reconstruction, contact element mesh alignment and error-driven mesh update mechanism.

8. The method for evaluating the axial load distribution and wear life of a connection component according to claim 1, characterized in that: The output life prediction results are: tooth surface axial load distribution, target wear value, predicted life, fitting parameters, which can be saved as a txt report.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method according to any one of claims 1 to 8 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 8 are implemented.

Images