Method for calculating time-varying meshing stiffness of internal meshing spur cylindrical gear pair considering bolt loosening failure influence

Abstract

This invention discloses a method for calculating the time-varying meshing stiffness of an internal meshing spur gear pair considering the influence of bolt loosening faults. The method is characterized by obtaining the tooth body compliance, base compliance, and tooth body compliance of the planetary gears and the gear ring using analytical methods; establishing a finite element static analysis model of the gear ring with bolt loosening faults; training an MLP neural network model with the extracted data to predict the base compliance of the gear ring under bolt loosening faults; and then establishing a gear bearing contact analysis model to solve for the time-varying meshing stiffness value. The advantage is that it achieves a deep integration of computational efficiency and prediction accuracy for time-varying meshing stiffness, significantly improving the efficiency and intelligence of dynamic simulation of internal meshing gear transmission systems while ensuring the accuracy of time-varying meshing stiffness calculation.

Description

Technical Field

[0001] This invention relates to the field of fault diagnosis technology in mechanical transmission dynamics, and in particular to a method for calculating the time-varying meshing stiffness of an internal meshing spur gear pair that takes into account the effects of bolt loosening faults. Background Technology

[0002] As a core power component, the efficiency of gear transmission systems directly affects the overall mechanical operation efficiency. When faced with adverse factors such as overload, lubrication loss, and loose bolts, gear transmission systems may experience faults such as cracks, pitting, galling, and even tooth breakage, severely hindering their smooth operation. Given the widespread use of gears in most transmission mechanisms, fault diagnosis of gears is particularly important. Currently, planetary gearbox components widely use bolted connections, and the tightness of these bolts has a significant impact on the time-varying meshing stiffness of the overall gear transmission system. Therefore, in experimental or engineering applications, real-time monitoring of the bolt connection status of planetary gearbox components is a crucial measure. This effectively optimizes the operating conditions of the gear transmission system, promptly identifies and addresses potential problems, and prevents situations that may cause permanent damage to the system. However, the core prerequisite for achieving accurate monitoring and early warning of bolt loosening is establishing a system dynamics model that can accurately characterize fault features, and then deeply analyzing the system response evolution under fault conditions through this model.

[0003] Accurately obtaining the time-varying meshing stiffness curve of a planetary gear transmission system with loose bolts is crucial when constructing a dynamic model of the system, as it determines the accuracy and reliability of the simulation. However, current methods for obtaining this stiffness curve have significant limitations: Firstly, while high-precision finite element simulation models can simulate complex boundaries, they are extremely computationally expensive (requiring significant computation time and hardware resources), inefficient, complex in modeling, have high professional barriers, and are prone to convergence problems in nonlinear solutions, making them unsuitable for the needs of rapid analysis, parametric research, or large-scale dynamic simulation in practical engineering applications. Secondly, while traditional analytical calculation methods are relatively efficient, they are mostly based on the assumption that the gear ring of the planetary gear set is externally rigidly fixed. However, under the condition of loose bolts, the gear ring is fixedly connected to the gearbox housing by bolts, and bolt loosening directly changes the boundary flexibility of the gear ring matrix, which is fundamentally different from the rigid boundary of traditional analytical models. Therefore, existing analytical models cannot accurately characterize the impact of changes in the gear ring matrix flexibility caused by bolt loosening on the time-varying meshing stiffness. Summary of the Invention

[0004] The technical problem to be solved by the present invention is to provide a method for calculating the time-varying meshing stiffness of an internal meshing spur gear pair that takes into account the effect of bolt loosening failure, which can realize the rapid and accurate acquisition of the time-varying meshing stiffness under bolt loosening failure conditions.

[0005] The technical solution adopted by this invention to solve the above-mentioned technical problems is: a time-varying meshing stiffness calculation method for internal meshing spur gear pairs considering the influence of bolt loosening faults, including the following specific steps:

[0006] (1) The tooth body compliance, base compliance and tooth body compliance of the planetary gear are obtained by analytical calculation;

[0007] (2) Establish a finite element static analysis model of the gear ring with bolt loosening fault;

[0008] (3) Apply a unit load to the finite element static analysis model of the gear ring and extract the displacement data of the loaded tooth and the adjacent tooth;

[0009] (4) Calculate the matrix compliance of the gear ring under bolt loosening fault using displacement data;

[0010] (5) Establish an MLP neural network model, take the measured radius of the load-bearing tooth and the force-bearing point on the adjacent tooth of the gear ring as the input feature, and take the calculated matrix compliance of the gear ring under bolt loosening fault as the output feature. Train the MLP neural network. After training, predict the matrix compliance of the gear ring when the measured radius of the load-bearing point of the load-bearing tooth under bolt loosening fault is taken as the input feature through the MLP neural network.

[0011] (6) Establish a gear bearing contact analysis model, and input the tooth body flexibility, base flexibility and tooth body flexibility of the planetary gear obtained by analytical method, and the base flexibility of the gear ring under bolt loosening fault into the gear bearing contact analysis model, solve for the time-varying meshing stiffness value of the planetary gear and the gear ring at different meshing positions of the loaded teeth in one meshing cycle, and thus obtain the time-varying meshing stiffness curve.

[0012] Furthermore, in step (1), the analytical method employs the potential energy method.

[0013] Furthermore, in step (1), the matrix compliance of the planetary gear includes matrix-induced tooth body local compliance and matrix-induced tooth body coupling compliance.

[0014] Furthermore, the specific establishment process of step (2) is as follows:

[0015] (2.1) Establish the finite element static analysis model of the gear ring;

[0016] (2.2) Set the elastic modulus of the gear body material of the gear ring to 10% of the elastic modulus of the base material of the gear ring. 5 times;

[0017] (2.3) Apply full-degree-of-freedom constraints to all bolt holes on the base of the gear ring to simulate the assembly boundary conditions under gear ring fastening;

[0018] (2.4) Release the full degree of freedom of at least one bolt hole to obtain the finite element static analysis model of the gear ring with bolt loosening fault.

[0019] Furthermore, in step (3), the specific extraction process of the displacement data of the loaded tooth and adjacent teeth of the gear ring is as follows:

[0020] (3.1) Apply along one of the teeth of the finite element static analysis model of the gear ring with bolt loosening fault established in step (2). x direction and y Unit load in the direction;

[0021] (3.2) Extracting data using a Python script x direction and y Displacement data of the loaded tooth and adjacent teeth under a unit load applied in the direction are generated into a text file.

[0022] Furthermore, in step (4), the matrix compliance of the gear ring includes matrix-induced local compliance of the tooth body and matrix-induced coupled compliance of the tooth body. The relationship between the matrix-induced local compliance of the gear ring under bolt loosening fault is as follows:

[0023] ,

[0024] in: Indicates unit composite load f The matrix-induced local compliance of the tooth body when acting on the loaded tooth of the gear ring. This represents the unit composite load acting on the loaded tooth. f The angle between its direction and the direction of its tangent on its circumference. and Correspondingly, this indicates the application of force to the loaded teeth of the gear ring. x Under a unit load in the direction, the point of force application on the loaded tooth is at x direction and y Displacement component generated by direction; and Correspondingly, this indicates that an application is made at the same position on the loaded tooth of the gear ring. y Under a unit load in the direction, the point of force application on the loaded tooth is at x direction and y Displacement component generated by direction;

[0025] The relationship between the matrix-induced tooth coupling compliance of the gear ring under bolt loosening fault is as follows:

[0026] ,

[0027] in: Indicates unit composite load f When acting on the loaded tooth of the gear ring, the adjacent tooth induces tooth body coupling compliance with the matrix of the loaded tooth; This indicates the angle between the direction of the load on an adjacent tooth and the direction of its circumferential tangent; z Indicates the number of teeth on the gear ring. and Correspondingly, this indicates the application of force to the loaded teeth of the gear ring. x Under a unit load in the direction, the force-bearing point of adjacent teeth is... x direction and y Displacement component generated by direction; and Correspondingly, this indicates that an application is made at the same position on the loaded tooth of the gear ring. y Under a unit load in the direction, the force-bearing point of adjacent teeth is... x direction and y The displacement component generated by the direction.

[0028] Furthermore, in step (5), when training the MLP neural network, a single-input single-output network based on the measurement point radius of the loaded tooth is established for the matrix-induced local flexibility of the tooth body of the gear ring, and a dual-input single-output network based on the dual measurement point radius of the loaded tooth and the adjacent tooth is established for the matrix-induced coupling flexibility of the tooth body of the gear ring.

[0029] Compared with the prior art, the advantages of the present invention are:

[0030] (1) Achieving a deep integration of computational efficiency and prediction accuracy for time-varying meshing stiffness; This method breaks through the limitations of traditional analytical methods, establishes a matrix compliance mapping framework for gear rings with bolt holes through finite element simulation, and introduces an MLP neural network to intelligently predict the matrix compliance of the gear ring under bolt loosening fault, effectively capturing the complex boundary deformation characteristics caused by bolt loosening. This fusion scheme of "finite element simulation + analysis" not only significantly improves the efficiency and intelligence level of dynamic simulation of internal meshing gear transmission system, but also ensures the computational accuracy of time-varying meshing stiffness;

[0031] (2) In the finite element static analysis model of the gear ring, this invention innovatively sets the elastic modulus of the gear body material to 10% of the elastic modulus of the matrix material. 5This mechanical processing method physically equates the tooth body of the gear ring to an "absolute rigid body," thereby ensuring that after applying a unit load, all displacements at the stress point are purely caused by the elastic deformation of the base body, achieving a complete decoupling of tooth deformation and base body deformation in terms of mechanical mechanism; moreover, it completely locks the flexibility of the tooth body itself, ensuring that the extracted network prediction target (i.e., the base flexibility of the gear ring) is absolute. This allows the prediction results of the neural network to achieve perfect closure of the mechanical boundary when combined with the calculation results of the analytical method, fundamentally guaranteeing the theoretical rigor and numerical reliability of the calculation of the time-varying meshing stiffness curve;

[0032] (3) The present invention applies a unit load to the rigid tooth body in the orthogonal directions of x and y respectively, and extracts the displacement data of the loaded tooth and the adjacent tooth in the two-dimensional plane. Then, the matrix compliance of the gear ring under bolt loosening fault along the actual meshing line direction is calculated by the displacement data. This principle ensures that the subsequent neural network model can accurately respond to the transient change of the normal load angle.

[0033] (4) The matrix-induced tooth coupling compliance of the tooth ring of the present invention is separately constructed and predicted. In addition, a dual-input feature mechanism based on the measurement point radius of the load-bearing tooth and the force point on the adjacent tooth is specially designed in the neural network model, which enables the present invention to capture the coupling behavior between adjacent teeth under fault conditions with high precision, filling the gap in the existing fault analysis model in coupling compliance calculation. Attached Figure Description

[0034] Figure 1 This is a flowchart of the present invention;

[0035] Figure 2 This is a schematic diagram of the internal meshing of the planetary gear and the gear ring of the present invention;

[0036] Figure 3 This is a schematic diagram showing the forces acting on the loaded tooth and adjacent teeth of the gear ring of the present invention;

[0037] Figure 4 This is a comparison chart of the time-varying meshing stiffness curves obtained by the present invention and the traditional finite element simulation model. Detailed Implementation

[0038] The present invention will be further described in detail below with reference to the accompanying drawings and embodiments.

[0039] To facilitate understanding of the overall technical solution, the various parts of planetary gear 1 and gear ring 2 are labeled and explained, where: 11 represents the base of the planetary gear, 12 represents the tooth of the planetary gear, 21 represents the base of the gear ring, and 22 represents the tooth of the gear ring.

[0040] like Figures 1-3As shown, the method for calculating the time-varying meshing stiffness of an internal meshing spur gear pair, taking into account the effects of bolt loosening faults, includes the following specific steps:

[0041] (1) The tooth compliance of planetary gear 1, the matrix-induced local compliance of planetary gear 1 and the matrix-induced coupled compliance of the tooth, as well as the tooth compliance of gear ring 2, are obtained by calculating the potential energy method.

[0042] (2) Establish a finite element static analysis model for gear ring 2 with bolt loosening fault, specifically as follows:

[0043] (2.1) Establish the finite element static analysis model of gear ring 2;

[0044] (2.2) Set the elastic modulus of the tooth body material of gear ring 2 to 10% of the elastic modulus of the base material of gear ring 2. 5 times;

[0045] (2.3) Apply full degree of freedom constraints to all bolt holes on the base 21 of the gear ring 2 (i.e. tighten all bolts 3) to simulate the assembly boundary conditions under gear ring fastening.

[0046] (2.4) Release the full degree of freedom of at least one bolt hole to obtain the finite element static analysis model of the gear ring 2 with bolt loosening fault;

[0047] (3.1) Apply along the finite element static analysis model of the gear ring 2 with bolt loosening fault established in step (2) to one of the teeth 22. x direction and y Unit load in the direction;

[0048] (3.2) Extracting data using a Python script x direction and y Displacement data of the loaded tooth 221 and adjacent tooth 222 under a unit load applied in the direction are generated into a text file;

[0049] (4) Calculate the matrix-induced local compliance and matrix-induced coupled compliance of gear ring 2 under bolt loosening fault using the displacement data extracted in step (3.2). The relationship between the matrix-induced local compliance of gear ring 2 under bolt loosening fault is as follows:

[0050] ,

[0051] in: Indicates unit composite load f (Right now x direction and y The matrix-induced local compliance of the tooth body when a composite load of unit load in the direction is applied to the loaded tooth 221 of the gear ring 2. This represents the unit composite load acting at the loaded tooth 221. f The angle between its direction and the direction of its tangent on its circumference. and Correspondingly, this indicates that an application is made to the loaded tooth 221 of the gear ring 2. x Under a unit load in the direction, the point of force application of the loaded tooth 221 is at x direction and y Displacement component generated by direction; and Correspondingly, this indicates that an application is made at the same position on the loaded tooth 221 of the gear ring 2. y Under a unit load in the direction, the point of force application of the loaded tooth 221 is at x direction and y Displacement component generated by direction;

[0052] The relationship between the matrix-induced tooth coupling compliance of gear ring 2 under bolt loosening fault is as follows:

[0053] ,

[0054] in: Indicates unit composite load f When the load-bearing tooth 221 of the gear ring 2 is subjected to the load-bearing tooth 221, the adjacent tooth 222 induces tooth body coupling compliance with the matrix of the load-bearing tooth 221. This indicates the angle between the direction of the load on adjacent tooth 222 and the direction of its circumferential tangent; z This indicates the number of teeth on gear 2. and Correspondingly, this indicates that when applying force to the loaded tooth 221 of the gear ring 2... x Under a unit load in the direction, the stress point of adjacent tooth 222 is at x direction and y Displacement component generated by direction; and Correspondingly, this indicates that an application is made at the same position on the loaded tooth 221 of the gear ring 2. y Under a unit load in the direction, the stress point of adjacent tooth 222 is at x direction and y The displacement components generated by the direction; each displacement component here is the displacement data extracted in step (3.2);

[0055] (5) Establish a multilayer perceptron (MLP) neural network model, and measure the radius of the force points on the loaded teeth 221 of the gear ring 2 in step (4). r 1. The measuring radius of the force-bearing point on adjacent tooth 222 r 2. As an input feature, the measuring point radius refers to the distance from the force-bearing point on the tooth body 22 of the gear ring 2 to the center O of the gear ring. The calculated value is then compared with the measuring point radius. r 1. rThe matrix-induced local compliance and matrix-induced coupled compliance of the gear ring 2 under bolt loosening fault at the stress point of 2 are used as output features. Furthermore, a measurement point radius based on the loaded tooth 221 is established for the matrix-induced local compliance of the gear ring 2. r A single-input single-output network is used to establish the matrix-induced tooth coupling compliance of gear ring 2 based on the radius of the dual measurement points of the loaded tooth 221 and the adjacent tooth 222. r 1. r The dual-input single-output network of 2 is used to train the MLP neural network. After training, the MLP neural network is used to predict the matrix-induced local compliance and matrix-induced coupling compliance of the tooth ring 2 when the radius of the measuring point of different stress points of the loaded tooth 221 under bolt loosening fault is used as the input feature.

[0056] (6) Establish a gear bearing contact analysis (LTCA) model, and input the tooth body compliance of planetary gear 1, the base-induced tooth body local compliance and the base-induced tooth body coupling compliance of planetary gear 1, the tooth body compliance of gear ring 2, and the base-induced tooth body local compliance and base-induced tooth body coupling compliance of gear ring 2 under bolt loosening fault into the LTCA model, solve for the time-varying meshing stiffness value of the loaded tooth 221 at different meshing positions in one meshing cycle of planetary gear 1 and gear ring 2, and thus obtain the time-varying meshing stiffness curve.

[0057] The time-varying meshing stiffness curves obtained using this method (i.e., analytical-finite element simulation model AFEM) are compared with those obtained using the traditional finite element simulation model FEM. Figure 4 As shown, the relative error of the time-varying meshing stiffness of the two is no more than 5%, which proves that the method can ensure the accuracy of acquisition while improving the acquisition efficiency of the time-varying meshing stiffness curve, and provides a high-fidelity and fast-response theoretical support for the condition monitoring and dynamic evaluation of planetary gear systems with bolt loosening faults.

[0058] The scope of protection of this invention includes, but is not limited to, the above embodiments. The scope of protection is defined by the claims. Any substitutions, modifications, or improvements to this technology that are easily conceived by those skilled in the art fall within the scope of protection of this invention.

Claims

1. A method for calculating the time-varying meshing stiffness of an internal meshing spur gear pair considering the effects of bolt loosening faults, characterized in that... The specific steps include the following: (1) The tooth body compliance, base compliance and tooth body compliance of the planetary gear are obtained by analytical calculation; (2) Establish a finite element static analysis model of the gear ring with bolt loosening fault. The establishment process is as follows: (2.1) Establish the finite element static analysis model of the gear ring; (2.2) Set the elastic modulus of the gear body material of the gear ring to 10% of the elastic modulus of the base material of the gear ring. 5 times; (2.3) Apply full-degree-of-freedom constraints to all bolt holes on the base of the gear ring to simulate the assembly boundary conditions under gear ring fastening; (2.4) Release the full degree of freedom of at least one bolt hole to obtain the finite element static analysis model of the gear ring with bolt loosening fault; (3) Apply a unit load to the finite element static analysis model of the gear ring and extract the displacement data of the loaded tooth and the adjacent tooth; (4) Calculate the matrix compliance of the gear ring under bolt loosening fault using displacement data; (5) Establish an MLP neural network model, take the measured radius of the load-bearing tooth and the force-bearing point on the adjacent tooth of the gear ring as the input feature, and take the calculated matrix compliance of the gear ring under bolt loosening fault as the output feature. Train the MLP neural network. After training, predict the matrix compliance of the gear ring when the measured radius of the load-bearing point of the load-bearing tooth under bolt loosening fault is taken as the input feature through the MLP neural network. (6) Establish a gear bearing contact analysis model, and input the tooth body flexibility, base flexibility and tooth body flexibility of the planetary gear obtained by analytical method, and the base flexibility of the gear ring under bolt loosening fault into the gear bearing contact analysis model, solve for the time-varying meshing stiffness value of the planetary gear and the gear ring at different meshing positions of the loaded teeth in one meshing cycle, and thus obtain the time-varying meshing stiffness curve.

2. The method for calculating the time-varying meshing stiffness of an internal meshing spur gear pair considering the influence of bolt loosening faults as described in claim 1, characterized in that: In step (1), the analytical method uses the potential energy method.

3. The method for calculating the time-varying meshing stiffness of an internal meshing spur gear pair considering the influence of bolt loosening faults as described in claim 1, characterized in that: In step (1), the matrix compliance of the planetary gear includes matrix-induced tooth body local compliance and matrix-induced tooth body coupling compliance.

4. The method for calculating the time-varying meshing stiffness of an internal meshing spur gear pair considering the influence of bolt loosening faults as described in claim 1, characterized in that: In step (3), the specific extraction process of the displacement data of the loaded tooth and adjacent teeth of the gear ring is as follows: (3.1) Apply unit loads along the x and y directions respectively to one of the teeth of the finite element static analysis model of the gear ring with bolt loosening fault established in step (2). (3.2) Extract the displacement data of the loaded tooth and adjacent teeth under unit load applied in the x and y directions using a Python script and form a text file.

5. The method for calculating the time-varying meshing stiffness of an internal meshing spur gear pair considering the influence of bolt loosening faults as described in claim 1, characterized in that: In step (4), the matrix compliance of the gear ring includes matrix-induced local compliance of the tooth body and matrix-induced coupled compliance of the tooth body. The relationship between the matrix-induced local compliance of the gear ring under bolt loosening fault is as follows: ， in: This represents the matrix-induced local compliance of the tooth structure when a unit composite load f is applied to the loaded tooth of the gear ring. This represents the angle between the direction of the unit combined load f acting on the loaded tooth and the direction of its circumferential tangent. and Correspondingly, it represents the displacement components in the x and y directions generated by the force-bearing point of the gear when a unit load in the x direction is applied to the loaded tooth of the gear ring; and Correspondingly, it represents the displacement components in the x and y directions of the load-bearing tooth when a unit load in the y direction is applied to the same position of the loaded tooth of the gear ring. The relationship between the matrix-induced tooth coupling compliance of the gear ring under bolt loosening fault is as follows: ， in: This represents the tooth body coupling compliance induced by adjacent teeth to the matrix of the loaded tooth when a unit composite load f is applied to the loaded tooth of the gear ring. The angle between the direction of the load on an adjacent tooth and the direction of its circumferential tangent is indicated; z represents the number of teeth on the gear ring. and Correspondingly, it represents the displacement components in the x and y directions generated by the force-bearing point of the adjacent tooth when a unit load in the x direction is applied to the loaded tooth of the gear ring; and Correspondingly, it represents the displacement components in the x and y directions generated by the force-bearing point of adjacent teeth when a unit load in the y direction is applied to the same position of the loaded teeth of the gear ring.

6. The time-varying meshing stiffness calculation method for internal meshing spur gear pairs considering the influence of bolt loosening faults as described in claim 5, characterized in that: In step (5), when training the MLP neural network, a single-input single-output network based on the measurement point radius of the loaded tooth is established for the matrix-induced local flexibility of the tooth body of the gear ring, and a dual-input single-output network based on the dual measurement point radius of the loaded tooth and the adjacent tooth is established for the matrix-induced coupling flexibility of the tooth body of the gear ring.

Images