A gear dynamic response analysis calculation method, device, equipment and medium

Abstract

The application discloses a gear dynamic response analysis calculation method, device, equipment and medium, the method comprises the steps that S1, the gear inherent frequency obtained from modal calculation, transient dynamics calculation and checking is screened out as the final excitation frequency which does not produce beat vibration phenomenon;S2, according to the coincidence degree of gear, the part of tooth number of gear loading is determined, the final excitation frequency, the determined part of tooth number is used to simulate loading to replace full tooth simulation loading to obtain excitation load curve;S3, according to the material characteristics of aviation gear, the corresponding gear damping ratio is determined;S4, according to the final excitation frequency, excitation load curve and gear damping ratio, dynamic response is solved, and the dynamic response result of target position is obtained.The application improves the calculation precision of the real dynamic response of gear in the calculation of damping state;Through the cyclic loading on a small amount of teeth, the meshing loaded process of gear is simulated, and the simulation efficiency and computer memory consumption are significantly improved.

Description

A method, apparatus, equipment and medium for dynamic response analysis and calculation of gears. Technical Field

[0001] This application relates to the field of aero-engine gear technology, and in particular to a method, apparatus, equipment and medium for dynamic response analysis and calculation of gears. Background Technology

[0002] Aero engines operate at high speeds over a wide range, making it difficult to completely avoid designing gears to resonate at all engine operating speeds. While gear modal analysis can calculate the resonant frequency and speed, it cannot determine the actual dynamic stress at the resonant speed, thus failing to accurately assess the risk of resonant failure at that speed. Currently, to verify whether gears will experience resonant failure, dynamic stress measurement tests are often conducted on the gears during engine testing or in component testing equipment.

[0003] Currently, gear design can only provide modal analysis results and the operating speed at the resonant speed point within the working range. It's difficult to accurately predict the vibration response of the gear at the resonant speed point. The reliability of the gear can only be determined by measuring dynamic stress in its physical state after machining. However, testing through experiments has several drawbacks:

[0004] 1) Testing and experimentation are costly and have a time lag. Once it is found that the dynamic stress of the gear is too high at the resonance speed point, the gear structure needs to be adjusted. Sometimes, major adjustments to the gear structure parameters are required to reduce the dynamic stress, and the adjustment process is often a continuous trial and error process.

[0005] 2) Sufficient space is required to install the testing equipment during the testing process. Therefore, not all gears can be used for dynamic stress measurement. Engines need to be modified for testing, which will also increase the processing cost and delay the research progress. Summary of the Invention

[0006] This application provides a method for analyzing and calculating the dynamic response of gears, which solves the technical problems of high testing and experimentation costs, lag, limited applicability, and delays in scientific research progress in existing technologies.

[0007] This application is achieved through the following solution:

[0008] A method for analyzing and calculating the dynamic response of gears, comprising the following steps:

[0009] S1. Select the excitation frequency that does not produce beat vibration from the gear natural frequencies obtained by modal calculation, transient dynamics calculation and verification as the final excitation frequency;

[0010] S2. Determine the number of teeth of the gear to be loaded based on the overlap ratio of the gears, and use the final excitation frequency and the determined number of teeth of the gear to perform simulated loading instead of full-tooth simulated loading to obtain the excitation load curve.

[0011] S3. Determine the appropriate gear damping ratio based on the material properties of the aircraft gears;

[0012] S4. Solve the dynamic response based on the final excitation frequency, excitation load curve and gear damping ratio to obtain the dynamic response result at the target position.

[0013] Further, step S1 specifically includes the following steps:

[0014] S11. Perform modal calculations on the gears;

[0015] S12. Select the natural frequency of a certain mode within the working range. ;

[0016] S13. Perform transient dynamics calculations, adding the frequency during the calculation. The excitation is a sinusoidal function, and the damping is defined as 0;

[0017] S14. Based on the beat vibration results given by transient dynamics calculation, calculate the excitation frequency at which beat vibration will not occur. :

[0018] ;

[0019] Among them, the shooting frequency Based on the beat response calculation results: ,in The period of the shot can be obtained directly by reading the response result;

[0020] S15. Perform transient dynamic calculation verification. During the calculation verification, add an excitation frequency of... The excitation frequency is determined by applying a sinusoidal function excitation with zero damping. If the response no longer exhibits beat vibration and increases continuously with time, then the excitation frequency is determined accordingly. This is the final excitation frequency.

[0021] Furthermore, step S2 specifically includes the following steps:

[0022] S21. Determine the number of teeth of the gear under load based on the gear overlap ratio ε.

[0023] S22. Simulate the alternating meshing of partially paired gears during operation based on the number of teeth Z0 of the loaded gear and the final excitation frequency to obtain the excitation load curve.

[0024] Preferably, in step S21, determining the number of teeth of the gear under load based on the gear overlap ratio ε specifically involves:

[0025] Z0 = [ε] + 1;

[0026] Where [ε] represents the degree of overlap ε, which is rounded down to the nearest integer and the decimal part is discarded.

[0027] Further, in step S22, when the overlap ratio 1 < ε < 2, the obtained excitation load curve is a trapezoidal wave. The amplitude of the flat top of the trapezoidal wave is determined according to the actual load on the gear, including alternating single-pair meshing areas and double-pair meshing areas, wherein:

[0028] Gear meshing cycle:

[0029] ;

[0030] Duration of engagement zone for a single pair of teeth:

[0031] ;

[0032] Duration of double-pair tooth engagement zone:

[0033] .

[0034] Further, in step S22, when the overlap ratio ε=2, the obtained excitation load curve changes from a trapezoidal wave to a triangular wave. The peak amplitude of the triangular wave is determined according to the actual load on the gear, and the peak time of the excitation load curve is shortened to 0. It includes alternating single-pair meshing areas and double-pair meshing areas, wherein:

[0035] Gear meshing cycle:

[0036] ;

[0037] Duration of engagement zone for a single pair of teeth:

[0038] ;

[0039] Duration of double-pair tooth engagement zone:

[0040] .

[0041] Further, in step S22, when the overlap ratio 2 < ε < 3, the obtained excitation load curve is a trapezoidal wave. The amplitude of the flat top of the trapezoidal wave is determined according to the actual load on the gear, including alternating double-pair meshing areas and triple-pair meshing areas, wherein:

[0042] Gear meshing cycle:

[0043] ;

[0044] Duration of double-pair tooth engagement zone:

[0045] t2=(3-ε)T;

[0046] Duration of the three pairs of teeth in the meshing zone:

[0047] t3=(ε-2)T.

[0048] This application also provides a gear dynamic response analysis and calculation device, including:

[0049] The excitation frequency determination module is used to select the excitation frequency that does not produce beat vibration as the final excitation frequency from the gear natural frequencies obtained by modal calculation, transient dynamics calculation and verification.

[0050] The excitation load curve determination module is used to determine the number of teeth of the gear that are partially loaded based on the gear overlap ratio, and to obtain the excitation load curve by using the final excitation frequency and the determined number of teeth of the gear that are partially loaded to perform simulated loading instead of full-tooth simulated loading.

[0051] The damping ratio determination module is used to determine the corresponding gear damping ratio based on the material properties of the aircraft gears.

[0052] The dynamic response solution module is used to solve the dynamic response based on the final excitation frequency, excitation load curve and gear damping ratio, and obtain the dynamic response result at the target position.

[0053] This application also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and running on the processor, wherein the processor executes the computer program to implement the steps of the gear dynamic response analysis and calculation method.

[0054] This application also provides a storage medium including a stored program that, when the program is executed, controls the device where the storage medium is located to perform the steps of the gear dynamic response analysis and calculation method.

[0055] Compared with the prior art, this application has the following advantages:

[0056] 1. The focus of this application is on the preprocessing of gear dynamic response calculation. Specifically, based on the beat vibration theory, a method for calculating the beat frequency is created by combining gear modal and transient dynamic calculations. This method can accurately calculate the beat frequency of the gear, effectively eliminating the difference between the gear's natural frequency and the gear's excitation frequency, thereby improving the calculation accuracy of the gear's real dynamic response under damped conditions.

[0057] 2. This application creates a gear loading method that simulates the gear meshing loading process by cyclically loading a small number of gear teeth. Compared with the full-tooth cyclic loading according to the real situation, this method can greatly reduce the amount of calculation and significantly improve the simulation efficiency and computer memory consumption.

[0058] In addition to the purposes, features, and advantages described above, this application has other purposes, features, and advantages. A further detailed description of this application will be provided below with reference to the figures. Attached Figure Description

[0059] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with this application and, together with the description, serve to explain the principles of this application.

[0060] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, for those skilled in the art, other drawings can be obtained based on these drawings without creative effort, wherein:

[0061] Figure 1 is a schematic flowchart of the gear dynamic response analysis and calculation method according to a preferred embodiment of this application;

[0062] Figure 2 is a schematic diagram of the beat vibration of a preferred embodiment of this application;

[0063] Figure 3 is a gear loading curve and a schematic diagram of single and double tooth meshing regions of a preferred embodiment of this application (ε=1.7);

[0064] Figure 4 is a schematic diagram of the gear loading curve and the meshing regions of double and triple pairs of teeth in a preferred embodiment of this application (ε=2.3);

[0065] Figure 5 is a schematic diagram of the gear dynamic response analysis and calculation device module of a preferred embodiment of this application;

[0066] Figure 6 is a schematic block diagram of an electronic device according to a preferred embodiment of this application;

[0067] Figure 7 is an internal structural diagram of a computer device according to a preferred embodiment of this application. Detailed Implementation

[0068] It should be understood that the specific embodiments described herein are merely illustrative of the technical solutions of this application and are not intended to limit this application.

[0069] To better understand the technical solution of this application, a detailed description will be provided below in conjunction with the accompanying drawings and specific implementation methods.

[0070] It should be noted that the executing entity in this embodiment can be a computing service device with data processing, network communication, and program execution functions, such as a tablet computer, personal computer, or mobile phone, or a gear dynamic response analysis computing device capable of performing the above functions. The following description uses a gear dynamic response analysis computing device as the executing entity to illustrate this embodiment and the subsequent embodiments.

[0071] As shown in Figure 1, a preferred embodiment of this application provides a method for analyzing and calculating the dynamic response of a gear, including the following steps:

[0072] S1. Select the excitation frequency that does not produce beat vibration from the gear natural frequencies obtained by modal calculation, transient dynamics calculation and verification as the final excitation frequency;

[0073] S2. Determine the number of teeth of the gear to be loaded based on the overlap ratio of the gears, and use the final excitation frequency and the determined number of teeth of the gear to perform simulated loading instead of full-tooth simulated loading to obtain the excitation load curve.

[0074] S3. Determine the corresponding gear damping ratio based on the material properties of the aircraft gear. For aircraft gears, the material is generally steel. Although the damping ratios of different gear structures are different, they are generally in the range of 0.06%-0.12%. Therefore, the value can be selected within this range based on engineering experience.

[0075] S4. Based on the final excitation frequency, excitation load curve and gear damping ratio, perform dynamic response calculation (dynamic response calculation can be performed using currently mature commercial software) to obtain the dynamic response results at the target position.

[0076] Compared with the prior art, this embodiment has the following advantages:

[0077] 1. The focus of this embodiment is on the preprocessing of gear dynamic response calculation. Specifically, based on beat vibration theory, a method for calculating beat frequency is created by combining gear modal and transient dynamic calculations. This method can accurately calculate the beat frequency of the gear, effectively eliminating the difference between the gear's natural frequency and the gear's excitation frequency, thereby improving the calculation accuracy of the gear's real dynamic response under damped conditions.

[0078] 2. This embodiment creates a gear loading method that simulates the gear meshing loading process by cyclically loading a small number of gear teeth. Compared with the full-tooth cyclic loading according to the real situation, this method can greatly reduce the amount of calculation and significantly improve the simulation efficiency and computer memory consumption.

[0079] Preferably, step S1 specifically includes the following steps:

[0080] S11. Perform modal calculations on the gears;

[0081] S12. Select the natural frequency of a certain mode within the working range. ;

[0082] S13. Perform transient dynamics calculations, adding the frequency during the calculation. The excitation is a sinusoidal function, and the damping is defined as 0;

[0083] S14. Based on the beat vibration results given by transient dynamics calculation, calculate the excitation frequency at which beat vibration will not occur. :

[0084] ;

[0085] Among them, the shooting frequency The results are derived from the beat response calculation (see Figure 2): ,in The period of the shot can be obtained directly by reading the response result;

[0086] S15. Perform transient dynamic calculation verification. During the calculation verification, add an excitation frequency of... The excitation frequency is determined by applying a sinusoidal function excitation with zero damping. If the response no longer exhibits beat vibration and increases continuously with time, then the excitation frequency is determined accordingly. This is the final excitation frequency.

[0087] Simulation analysis of gear dynamic stress typically employs transient dynamics methods. However, when calculating the dynamic response, if the excitation frequency is very close to the gear's resonant frequency, beat vibration can occur, causing periodic fluctuations in the gear's dynamic stress. This effect is particularly pronounced in low-damped systems (damping ratio ≤ 0.2%). Consequently, the obtained dynamic stress results may differ significantly from the actual results. Therefore, it is crucial that the excitation frequency closely matches the gear's natural frequency.

[0088] When solving for the transient response of gears, the direct integration method is often used. This method often leads to a difference between the natural frequencies of the model and those obtained from modal calculations. If the natural frequencies obtained from modal calculations are used directly, beat vibrations will occur when calculating the gear response, resulting in significant deviations in the calculation results. Therefore, this embodiment calculates the gear's natural frequencies without beat vibrations before calculating the transient dynamics. By accurately calculating the gear's beat frequencies, the difference between the gear's natural frequencies and the gear's excitation frequencies is effectively eliminated, thereby improving the calculation accuracy of the gear's true dynamic response under damped conditions.

[0089] Preferably, step S2 specifically includes the following steps:

[0090] S21. Determine the number of teeth in the loaded portion of the gear based on the gear overlap ratio ε:

[0091] Z0 = [ε] + 1;

[0092] Where [ε] represents the degree of overlap ε, which is rounded down to the nearest integer and the decimal part is discarded;

[0093] S22. Simulate the alternating meshing of partially paired gears during operation based on the number of teeth Z0 of the loaded gear and the final excitation frequency to obtain the excitation load curve.

[0094] During gear meshing, the number of teeth involved in meshing changes over time, and the meshing position also changes. It is impossible to accurately simulate the actual meshing situation. Therefore, it is necessary to perform equivalence. In this embodiment, the number of teeth loaded on the gear is determined according to the overlap ratio ε of the gear, Z0 = [ε] + 1, where [ε] represents the overlap ratio ε. The integer part is taken directly, and the decimal part is discarded. This embodiment only uses a portion of the teeth to simulate the alternating meshing of different pairs of gears during the working process, which will greatly improve the simulation calculation efficiency.

[0095] Preferably, in step S22, when the overlap ratio 1 < ε < 2, taking ε = 1.7 as an example, the obtained excitation load curve is a trapezoidal wave. The amplitude of the flat top of the trapezoidal wave is determined according to the actual load on the gear, including alternating single-pair meshing areas and double-pair meshing areas, wherein:

[0096] Gear meshing cycle:

[0097] ;

[0098] Duration of engagement zone for a single pair of teeth:

[0099] ;

[0100] Duration of double-pair tooth engagement zone:

[0101] .

[0102] This embodiment takes an overlap ratio of 1.7 as an example, and loads two adjacent teeth. The loading curves of the two teeth are shown in Figure 3 (the amplitude in the figure represents the magnitude of the gear pair meshing force). This embodiment uses only two teeth to simulate the alternating single and double meshing during gear operation, which greatly improves the simulation calculation efficiency. If a gear has 31 teeth, and the load is calculated for 10 rotations, in real working conditions, each gear is loaded. During simulation, even unloaded gears would be loaded, albeit with a load of 0, requiring calculation and significantly increasing the computational load and computer memory consumption. This embodiment uses only two teeth to achieve this full-tooth loading substitution, thus greatly reducing the computational load. The gear load is assumed to be a trapezoidal wave (a trapezoidal wave is a relatively ideal load curve, which can be further optimized based on the actual meshing load, but has little impact on the response result, and the optimized curve is also applicable to the method of this embodiment). The peak amplitude of the trapezoidal wave is determined based on the actual gear load.

[0103] Preferably, in step S22, when the overlap ratio ε=2, similar to the case of 1<ε<2, the obtained excitation load curve changes from a trapezoidal wave to a triangular wave. The peak amplitude of the triangular wave is determined according to the actual load on the gear, and the peak time of the excitation load curve is shortened to 0. It includes alternating single-pair meshing areas and double-pair meshing areas, wherein:

[0104] Gear meshing cycle:

[0105] ;

[0106] Duration of engagement zone for a single pair of teeth:

[0107] ;

[0108] Duration of double-pair tooth engagement zone:

[0109] .

[0110] Preferably, in step S22, when the overlap ratio 2 < ε < 3, similar to the case of 1 < ε < 2, the obtained excitation load curve is a trapezoidal wave. The amplitude of the flat top of the trapezoidal wave is determined according to the actual load on the gear, including alternating double-pair meshing areas and triple-pair meshing areas, wherein:

[0111] Gear meshing cycle:

[0112] ;

[0113] Duration of double-pair tooth engagement zone:

[0114] t2=(3-ε)T;

[0115] Duration of the three pairs of teeth in the meshing zone:

[0116] t3=(ε-2)T.

[0117] Taking ε=2.3 as an example, the loading curve is shown in Figure 4. In the figure, (1,2), (1,3), and (2,3) represent the gears under load as (tooth 1, tooth 2), (tooth 1, tooth 3), and (tooth 2, tooth 3), respectively. The case of ε≥3 can be deduced by analogy.

[0118] As shown in Figure 5, another preferred embodiment of this application also provides a gear dynamic response analysis and calculation device, including:

[0119] The excitation frequency determination module is used to select the excitation frequency that does not produce beat vibration as the final excitation frequency from the gear natural frequencies obtained by modal calculation, transient dynamics calculation and verification.

[0120] The excitation load curve determination module is used to determine the number of teeth of the gear that are partially loaded based on the gear overlap ratio, and to obtain the excitation load curve by using the final excitation frequency and the determined number of teeth of the gear that are partially loaded to perform simulated loading instead of full-tooth simulated loading.

[0121] The damping ratio determination module is used to determine the corresponding gear damping ratio based on the material properties of the aircraft gears.

[0122] The dynamic response solution module is used to solve the dynamic response based on the final excitation frequency, excitation load curve and gear damping ratio, and obtain the dynamic response result at the target position.

[0123] The gear dynamic response analysis and calculation device provided in this embodiment adopts the gear dynamic response analysis and calculation method in the above embodiments, solving the technical problems of high testing and experimentation costs, lag, limited applicability, and delays in scientific research progress in existing technologies. Compared with the prior art, the beneficial effects of the gear dynamic response analysis and calculation device provided in this embodiment are the same as those of the gear dynamic response analysis and calculation method provided in the above embodiments, and other technical features in the gear dynamic response analysis and calculation device are the same as those disclosed in the methods of the above embodiments, and will not be repeated here.

[0124] As shown in Figure 6, a preferred embodiment of this embodiment also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the steps of the gear dynamic response analysis and calculation method in the above embodiment.

[0125] This embodiment provides an electronic device that employs the gear dynamic response analysis and calculation method described in the above embodiments, solving the technical problems of high testing and experimental costs, lag, limited applicability, and delays in scientific research progress in existing technologies. Compared with the prior art, the beneficial effects of the electronic device provided in this embodiment are the same as those of the gear dynamic response analysis and calculation method provided in the above embodiments, and other technical features of the electronic device are the same as those disclosed in the methods of the above embodiments, and will not be repeated here.

[0126] As shown in Figure 7, a preferred embodiment of this invention also provides a computer device, which can be a terminal or a liveness detection server, and its internal structure is shown in Figure 7. The computer device includes a processor, memory, and a network interface connected via a system bus. The processor provides computing and control capabilities. The memory includes a non-volatile storage medium and internal memory. The non-volatile storage medium stores an operating system and computer programs. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage medium. The network interface is used to communicate with other external computer devices via a network connection. When the computer program is executed by the processor, it implements the steps of the aforementioned gear dynamic response analysis and calculation method.

[0127] Those skilled in the art will understand that the structure shown in FIG7 is merely a block diagram of a portion of the structure related to the solution of this embodiment, and does not constitute a limitation on the computer device to which the solution of this embodiment is applied. The specific computer device may include more or fewer components than shown in the figure, or combine certain components, or have different component arrangements.

[0128] The computer equipment provided in this application employs the gear dynamic response analysis and calculation method described in the above embodiments, solving the technical problems of high testing and experimental costs, lag, limited applicability, and delays in scientific research progress in existing technologies. Compared with the prior art, the beneficial effects of the computer equipment provided in this embodiment are the same as those of the gear dynamic response analysis and calculation method provided in the above embodiments, and other technical features in the electronic device are the same as those disclosed in the methods of the above embodiments, and will not be repeated here.

[0129] A preferred embodiment of this example also provides a storage medium, which includes a stored program that, when the program is executed, controls the device where the storage medium is located to perform the steps of the gear dynamic response analysis and calculation method described in the above embodiment.

[0130] This application has been successfully used in the dynamic simulation design of gears for a certain aero-engine, and the calculation results are reliable and can meet the engine's requirements for accessory transmission components.

[0131] It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although a logical order is shown in the flowchart, in some cases the steps shown or described may be executed in a different order than that shown here.

[0132] If the functions described in this embodiment are implemented as software functional units and sold or used as independent products, they can be stored in one or more computing device-readable storage media. Based on this understanding, the parts of this embodiment that contribute to the prior art or the technical solution can be embodied in the form of a software product. This software product is stored in a storage medium and includes several instructions to cause a computing device (which may be a personal computer, server, mobile computing device, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this embodiment. The aforementioned storage media include: USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, optical disks, and other media capable of storing program code.

[0133] Those skilled in the art will understand that embodiments of this example can be provided as methods, systems, or computer program products. Therefore, this example can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this example can take the form of a computer program product implemented 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. The solutions in this example can be implemented using various computer languages, such as the object-oriented programming language C++ and the embedded programming language C.

[0134] This embodiment is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to this embodiment. It should 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, create means for implementing the functions specified in one or more blocks of the flowchart illustrations and / or one or more blocks of the block diagrams.

[0135] 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 that implement the functions specified in one or more flowcharts and / or one or more block diagrams.

[0136] These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process, such that the instructions, which execute on the computer or other programmable apparatus, provide steps for implementing the functions specified in one or more flowcharts and / or one or more block diagrams.

[0137] This embodiment also provides a computer program product, including a computer program that, when executed by a processor, implements the steps of the gear dynamic response analysis and calculation method described above.

[0138] The computer program product provided in this embodiment solves the technical problems of high testing and experimentation costs, delays, limited applicability, and slowdowns in scientific research in existing technologies. Compared with the prior art, the beneficial effects of the computer program product provided in this embodiment are the same as those of the gear dynamic response analysis and calculation method provided in the above embodiments, and will not be repeated here.

[0139] Although preferred embodiments of this embodiment have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments as well as all changes and modifications falling within the scope of this embodiment.

[0140] Obviously, those skilled in the art can make various modifications and variations to this embodiment without departing from the spirit and scope of this embodiment. Therefore, if these modifications and variations of this embodiment fall within the scope of the claims of this embodiment and their equivalents, this embodiment is also intended to include these modifications and variations.

Claims

1. A method for analyzing and calculating the dynamic response of gears, characterized in that, The process includes the following steps: S1. Selecting the excitation frequency that does not produce beat vibration from the gear's natural frequencies obtained through modal calculation, transient dynamics calculation, and verification as the final excitation frequency. Specifically, this includes: S11. Performing modal calculations on the gear; S12. Selecting the natural frequency f1 of a certain mode within the working range; S13. Performing transient dynamics calculations, adding a sinusoidal excitation with frequency f1 and defining damping as 0; S14. Calculating the excitation frequency f that will not produce beat vibration based on the beat vibration results given by the transient dynamics calculations. Among them, beat frequency Based on the beat response calculation results: ,in The period of the beat is directly obtained by reading the response result; S15, perform transient dynamic calculation verification. During the calculation verification, add a sinusoidal function excitation with an excitation frequency of f and define the damping as 0. If the response result no longer shows the beat phenomenon and the response increases continuously with time, then the excitation frequency f is determined as the final excitation frequency; S2, determine the number of teeth loaded on the gear according to the gear overlap ratio. Use the final excitation frequency and the determined number of teeth loaded on the gear to perform simulated loading instead of full-tooth simulated loading to obtain the excitation load curve. The specific steps include: S21, determine the number of teeth loaded on the gear according to the gear overlap ratio ε; S22, simulate the situation of partial pair meshing and alternation during the gear operation according to the number of teeth loaded on the gear Z0 and the final excitation frequency to obtain the excitation load curve; S3, determine the corresponding gear damping ratio according to the material properties of the aerospace gear; S4, solve the dynamic response according to the final excitation frequency, excitation load curve and gear damping ratio to obtain the dynamic response result of the target position.

2. The gear dynamic response analysis and calculation method according to claim 1, characterized in that, In step S21, the number of teeth loaded on the gear is determined according to the gear overlap ε as follows: Z0 = [ε] + 1; where [ε] represents the overlap ε being rounded down, taking the integer part directly and discarding the decimal part.

3. The gear dynamic response analysis and calculation method according to claim 2, characterized in that, In step S22, when the overlap ratio 1 < ε < 2, the obtained excitation load curve is a trapezoidal wave. The amplitude of the flat top of the trapezoidal wave is determined according to the actual load on the gear, including alternating single-pair meshing areas and double-pair meshing areas, wherein: gear meshing cycle: Duration of engagement zone for a single pair of teeth: Duration of double-pair tooth engagement zone: 。 4. The gear dynamic response analysis and calculation method according to claim 2, characterized in that, In step S22, when the overlap ratio ε=2, the obtained excitation load curve changes from a trapezoidal wave to a triangular wave. The amplitude of the flat top of the triangular wave is determined according to the actual load on the gear. The flat top time of the excitation load curve is shortened to 0, including alternating single-pair meshing areas and double-pair meshing areas, wherein: gear meshing cycle: Duration of engagement zone for a single pair of teeth: Duration of double-pair tooth engagement zone: 。 5. The gear dynamic response analysis and calculation method according to claim 2, characterized in that, In step S22, when the overlap ratio 2 < ε < 3, the obtained excitation load curve is a trapezoidal wave. The amplitude of the flat top of the trapezoidal wave is determined according to the actual load on the gear, including alternating double-pair meshing areas and triple-pair meshing areas, wherein: gear meshing cycle: Duration of double-pair tooth meshing zone: t2=(3-ε)T; Duration of triple-pair tooth meshing zone: t3=(ε-2)T.

6. A gear dynamic response analysis and calculation device, characterized in that, include: The excitation frequency determination module is used to select the excitation frequency that will not produce beat vibration from the gear's natural frequencies obtained through modal calculation, transient dynamics calculation, and verification as the final excitation frequency. Specifically, it is used for: performing modal calculations on the gear; selecting the natural frequency f1 of a certain mode within the working range; performing transient dynamics calculations, adding a sinusoidal excitation with frequency f1 and defining damping as 0 during the calculation; and calculating the excitation frequency f that will not produce beat vibration based on the beat vibration results given by the transient dynamics calculation. Among them, beat frequency Based on the beat response calculation results: ,in The period of the beat is directly obtained by reading the response result. Transient dynamic calculation verification is performed. During the verification, a sinusoidal excitation with an excitation frequency of f is added, and the damping is defined as 0. If the response result no longer shows beat vibration and the response continuously increases with time, then the excitation frequency f is determined as the final excitation frequency. The excitation load curve determination module is used to determine the number of partially loaded teeth of the gear based on the gear contact ratio. The final excitation frequency and the determined number of partially loaded teeth are used to simulate loading instead of full-tooth simulation to obtain the excitation load curve. Specifically, it is used to: determine the number of partially loaded teeth of the gear based on the gear contact ratio ε; and simulate the alternating pair meshing of the gears during operation based on the number of partially loaded teeth Z0 and the final excitation frequency to obtain the excitation load curve. The damping ratio determination module is used to determine the corresponding gear damping ratio based on the material properties of the aerospace gears. The dynamic response solution module is used to solve the dynamic response based on the final excitation frequency, the excitation load curve, and the gear damping ratio to obtain the dynamic response result at the target position.

7. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the steps of the gear dynamic response analysis and calculation method as described in any one of claims 1 to 5.

8. A storage medium comprising a stored program, characterized in that, When the program is running, it controls the device containing the storage medium to perform the steps of the gear dynamic response analysis and calculation method as described in any one of claims 1 to 5.

Images