Blade high-cycle fatigue analysis method and system considering crystal orientation angle influence

Abstract

The invention discloses a blade high-cycle fatigue analysis method and system considering crystal orientation angle influence, and the method comprises the steps: calculating the maximum vibration stress at all crystal orientation angles according to the measured maximum vibration stress and the actual crystal orientation angle of a measured blade, and taking the maximum vibration stress as the statistical maximum vibration stress; according to the method, the influence of the crystal angle change on the statistical maximum vibration stress is fully considered, the confidence coefficient of the blade high-cycle fatigue evaluation result obtained through analysis is higher, high-cycle fatigue analysis can be carried out on the blade more accurately, and the risk of blade high-cycle fatigue failure is reduced.

Description

Technical Field

[0001] The present invention relates to the technical field of blade fatigue analysis, and in particular, to a high-cycle fatigue analysis method and system for blades considering the influence of crystal orientation angles, an electronic device, and a computer-readable storage medium. Background Art

[0002] Aeroengine blades usually adopt single-crystal blades or directionally solidified blades. Among them, the performance of single-crystal blades is better than that of directionally solidified blades. In order to ensure the safety of the engine, it is necessary to conduct high-cycle fatigue assessment on the blades. The high-cycle fatigue assessment of aeroengine blades is also an important part of the "Airworthiness Regulations for Aeroengines". The high-cycle fatigue assessment method based on vibration strain test data is the conventional means for current blade high-cycle fatigue assessment. This method directly converts and solves the maximum strain corresponding to each order obtained from blade tests to obtain the vibration stress value, and uses the above vibration stress value as the maximum vibration stress value in the statistical sense to carry out high-cycle fatigue assessment. However, the maximum vibration stress needs to consider the influence of different crystal orientation angles, different working conditions, and different processing errors on the vibration stress. In particular, the influence of crystal orientation angles on the vibration stress is relatively large. When the actual crystal orientation angle of the blade differs greatly from the theoretical crystal orientation angle, the maximum vibration stress obtained by the existing method may be significantly smaller than the actual maximum vibration stress. This means that when the high-cycle fatigue margin of the blade evaluated by the existing method is small, the vibration stress of the actually produced aeroengine blade may exceed the allowable value. Therefore, the existing blade high-cycle fatigue analysis method does not consider the influence of crystal orientation angles, resulting in poor accuracy of blade high-cycle fatigue analysis. Summary of the Invention

[0003] The present invention provides a high-cycle fatigue analysis method and system for blades considering the influence of crystal orientation angles, an electronic device, and a computer-readable storage medium. It fully considers the influence of crystal angle changes on the maximum vibration stress in statistics, and the confidence level of the obtained high-cycle fatigue assessment result of the blade is higher. It can carry out high-cycle fatigue analysis on the blade more precisely and reduce the risk of blade high-cycle fatigue failure.

[0004] According to one aspect of the present invention, a high-cycle fatigue analysis method for blades considering the influence of crystal orientation angles is provided, including the following steps: Obtain the steady-state stress of the blade under various working conditions; Obtain the vibration characteristics of the blade under various working conditions, and conduct vibration characteristic analysis to determine the vibration orders concerned with high-cycle fatigue of the blade; Conduct a dynamic stress measurement test on the blade to obtain the measured maximum vibration stress of the measured blade, and obtain the actual crystal orientation angle of the measured blade; For each vibration order, the maximum vibration stress at all crystal orientation angles is calculated based on the measured maximum vibration stress of the blade under test and the actual crystal orientation angle, and it is taken as the maximum vibration stress statistically. Based on the maximum vibration stress statistically and the steady-state stress, a GoodMan diagram is plotted to conduct high-cycle fatigue analysis of the blade.

[0005] Furthermore, the process of calculating the maximum vibration stress at all crystal orientation angles based on the measured maximum vibration stress of the blade under test and the actual crystal orientation angle includes the following: Conduct simulation analysis according to the actual crystal orientation angle of the blade under test to obtain the sensitivity of the strain gauge patch position corresponding to the measured maximum vibration stress at the actual crystal orientation angle, and convert the measured maximum vibration stress according to the sensitivity to obtain the actual maximum vibration stress of the blade under test. Carry out finite element modal analysis of the blade to obtain the stress coefficients of the blade at different crystal orientation angles, and screen out the maximum stress coefficient. Based on the maximum stress coefficient, the stress coefficient at the actual crystal orientation angle, and the actual maximum vibration stress of the blade under test, calculate the maximum vibration stress of the blade at all crystal orientation angles.

[0006] Furthermore, the maximum vibration stress of the blade at all crystal orientation angles is calculated based on the following formula: ; where represents the actual crystal orientation angle of the blade under test, represents the crystal orientation angle corresponding to the maximum stress coefficient, represents the stress coefficient at the actual crystal orientation angle, represents the maximum stress coefficient, represents the actual maximum vibration stress of the blade under test, represents the maximum vibration stress of the blade at all crystal orientation angles.

[0007] Furthermore, after calculating the maximum vibration stress at all crystal orientation angles, the following is also included: Considering the influence of the crystal orientation angle on the exciting force, correct the maximum vibration stress at all crystal orientation angles.

[0008] Furthermore, the process of considering the influence of the crystal orientation angle on the exciting force and correcting the maximum vibration stress at all crystal orientation angles includes the following: First, based on the typical curve of the blade aerodynamic pressure varying with the rotational speed, the aerodynamic pressures of the blade at the actual crystal orientation angle and the crystal orientation angle corresponding to the maximum stress coefficient are obtained respectively through rotational speed interpolation. Then, according to the principle that the aerodynamic excitation force is proportional to the aerodynamic pressure, the ratio of the aerodynamic pressures at the actual crystal orientation angle and the crystal orientation angle corresponding to the maximum stress coefficient is used to correct the maximum vibration stress at all crystal orientation angles, and the corrected maximum vibration stress is calculated.

[0009] Further, the corrected maximum vibration stress is calculated based on the following formula: ; Where, represents the actual crystal orientation angle of the measured blade, represents the crystal orientation angle corresponding to the maximum stress coefficient, represents the stress coefficient at the actual crystal orientation angle, represents the maximum stress coefficient, represents the actual maximum vibration stress of the measured blade, represents the maximum vibration stress of the blade at all crystal orientation angles, represents the crystal orientation angle the aerodynamic pressure in the resonance state at, represents the aerodynamic pressure in the resonance state at the actual crystal orientation angle a0.

[0010] In addition, the present invention also provides a blade high-cycle fatigue analysis system considering the influence of the crystal orientation angle, including: A steady-state stress analysis module for obtaining the steady-state stress of the blade under various working conditions; A vibration characteristic analysis module for obtaining the vibration characteristics of the blade under various working conditions and performing vibration characteristic analysis to determine the vibration order concerned about the blade high-cycle fatigue; A dynamic stress measurement and crystal orientation analysis module for conducting a dynamic stress measurement test on the blade to obtain the measured maximum vibration stress of the measured blade and obtaining the actual crystal orientation angle of the measured blade; A maximum vibration stress calculation module for, for each vibration order, calculating the maximum vibration stress at all crystal orientation angles based on the measured maximum vibration stress and the actual crystal orientation angle of the measured blade and taking it as the statistically maximum vibration stress; A blade high-cycle fatigue analysis module for drawing a GoodMan diagram based on the statistically maximum vibration stress and the steady-state stress and conducting blade high-cycle fatigue analysis.

[0011] Further, the maximum vibration stress calculation module is also used to correct the maximum vibration stress at all crystal orientation angles considering the influence of the crystal orientation angle on the excitation force after calculating the maximum vibration stress at all crystal orientation angles.

[0012] In addition, the present invention further provides an electronic device, including a processor and a memory, wherein a computer program is stored in the memory, and the processor is configured to execute the steps of the method described above by invoking the computer program stored in the memory.

[0013] In addition, the present invention further provides a computer-readable storage medium for storing a computer program for blade high-cycle fatigue analysis considering the influence of crystal orientation angle, and the computer program executes the steps of the method described above when running on a computer.

[0014] The present invention has the following beneficial effects: The blade high-cycle fatigue analysis method of the present invention considering the influence of crystal orientation angle calculates the maximum vibration stress at all crystal orientation angles based on the measured maximum vibration stress and the actual crystal orientation angle of the blade to be measured, and uses it as the maximum vibration stress in statistics, fully considering the influence of crystal angle change on the maximum vibration stress in statistics. The confidence level of the obtained blade high-cycle fatigue evaluation result is higher, and the blade high-cycle fatigue analysis can be carried out more accurately, reducing the risk of blade high-cycle fatigue failure.

[0015] In addition, the blade high-cycle fatigue analysis system of the present invention considering the influence of crystal orientation angle also has the above advantages.

[0016] In addition to the purposes, features and advantages described above, the present invention has other purposes, features and advantages. The present invention will be further described in detail below with reference to the drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] The drawings constituting a part of this application are used to provide a further understanding of the present invention. The schematic embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation to the present invention. In the drawings: Figure 1 is a schematic flowchart of a blade high-cycle fatigue analysis method considering the influence of crystal orientation angle according to a preferred embodiment of the present application; Figure 2 is Figure 1 a sub-flowchart of step S4 in Figure 3 is a schematic diagram of the crystal orientation angle of a single crystal blade according to a preferred embodiment of the present application; Figure 4 is a schematic diagram of a typical curve of blade aerodynamic pressure varying with rotational speed according to a preferred embodiment of the present application; Figure 5 is a schematic diagram of the module structure of a blade high-cycle fatigue analysis system considering the influence of crystal orientation angle according to another embodiment of the present application. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0018] It should be noted that, without conflict, the embodiments in the present application and the features in the embodiments may be combined with each other. The following will detail the present application with reference to the accompanying drawings and in combination with the embodiments.

[0019] Referring to Figure 1 , a preferred embodiment of the present application provides a method for high-cycle fatigue analysis of blades considering the influence of crystal orientation angles, including the following: Step S1: Obtain the steady-state stress of the blade under various working conditions; Step S2: Obtain the vibration characteristics of the blade under various working conditions, and conduct vibration characteristic analysis to determine the vibration orders concerned with high-cycle fatigue of the blade; Step S3: Conduct a dynamic stress measurement test on the blade to obtain the measured maximum vibration stress of the measured blade, and obtain the actual crystal orientation angle of the measured blade; Step S4: For each vibration order, calculate the maximum vibration stress at all crystal orientation angles based on the measured maximum vibration stress and the actual crystal orientation angle of the measured blade, and take it as the statistically maximum vibration stress; Step S5: Draw a Goodman diagram based on the statistically maximum vibration stress and the steady-state stress, and conduct high-cycle fatigue analysis of the blade.

[0020] It can be understood that the method for high-cycle fatigue analysis of blades considering the influence of crystal orientation angles in this embodiment calculates the maximum vibration stress at all crystal orientation angles based on the measured maximum vibration stress and the actual crystal orientation angle of the measured blade, and takes it as the statistically maximum vibration stress, fully considering the influence of crystal angle changes on the statistically maximum vibration stress. The confidence level of the obtained high-cycle fatigue assessment result of the blade is higher, and high-cycle fatigue analysis of the blade can be carried out more precisely, reducing the risk of high-cycle fatigue failure of the blade.

[0021] Specifically, in the step S1, first analyze and obtain the steady-state stress of the blade under various actual working conditions by the finite element method. Among them, the specific analysis process of the steady-state stress belongs to the prior art and will not be elaborated here. In addition, the blades referred to in the present invention include single-crystal blades and directionally solidified blades. The following description will take the single-crystal blade as an example for illustrative purposes.

[0022] In the step S2, analyze and obtain the vibration characteristics of the single-crystal blade under various actual working conditions by the finite element method, and determine the vibration orders concerned with high-cycle fatigue of the blade, that is, the dangerous orders, according to the vibration characteristic analysis. Among them, the specific analysis process belongs to the prior art and will not be elaborated here.

[0023] In the step S3, at the current stage, since the actual vibration stress of the blade cannot be obtained through calculation and analysis, it is necessary to conduct a dynamic stress measurement test on the single-crystal blade according to the relevant specifications of "Measurement Point Design for Blade Vibration Stress Measurement Test" and the "Technical Requirements for Blade Dynamic Stress Measurement Test" to obtain the strains corresponding to each vibration order measured by multiple strain gauges, and then screen out the maximum strain of each vibration order, and the maximum vibration stress measured of the blade to be measured at each vibration order can be obtained through conversion. In addition, the actual crystal orientation angle of the blade to be measured is obtained through crystal orientation analysis before the test. Among them, the specific process of crystal orientation analysis belongs to the prior art and will not be elaborated here.

[0024] It can be understood that the high-cycle fatigue analysis of the blade for airworthiness certification is to demonstrate that the maximum vibration stress under actual service conditions is less than the high-cycle fatigue limit of the material. At this time, it is considered that the blade will not suffer from high-cycle fatigue failure. The maximum vibration stress needs to consider the influence of different crystal orientation angles, different working conditions, and different processing errors on the vibration stress. Since the dynamic stress test can only measure a limited number of specific blades, it is necessary to convert the vibration stress of the measured specific blades to the maximum vibration stress values of all blades of the batch-produced engine under different service conditions through conversion or correction, that is, the maximum vibration stress values at all crystal orientation angles, so as to obtain the maximum vibration stress value in statistics. Therefore, in the step S4, for each vibration order determined in the step S2, the maximum vibration stress at all crystal orientation angles is calculated according to the measured maximum vibration stress and the actual crystal orientation angle of the blade to be measured, and it is used as the maximum vibration stress in statistics. Among them, as Figure 2 shown, the process of calculating the maximum vibration stress at all crystal orientation angles according to the measured maximum vibration stress and the actual crystal orientation angle of the blade to be measured includes the following contents: Step 41: Conduct a simulation analysis according to the actual crystal orientation angle of the blade to be measured to obtain the sensitivity of the strain gauge patch position corresponding to the measured maximum vibration stress at the actual crystal orientation angle, and convert the actual maximum vibration stress of the blade to be measured according to the sensitivity and the measured maximum vibration stress; Step 42: Carry out a finite element modal analysis of the blade to obtain the stress coefficients of the blade at different crystal orientation angles, and screen out the maximum stress coefficient; Step 43: Calculate the maximum vibration stress of the blade at all crystal orientation angles according to the maximum stress coefficient, the stress coefficient at the actual crystal orientation angle, and the actual maximum vibration stress of the blade to be measured.

[0025] Specifically, due to the inherent properties of strain gauges, in order to obtain high-precision blade strain values, during the dynamic stress measurement test, the strain gauge bonding position needs to be in a flat and relatively low-temperature area, avoiding bonding on uneven positions such as rounded corners, chamfers, and holes. Excessive temperature can cause the strain gauge to fail, and the above-mentioned area is generally not at the actual maximum vibration stress position of the measured blade. Therefore, it is necessary to convert it to the actual maximum vibration stress value of the measured blade according to the sensitivity. Here, the sensitivity is the ratio of the modal stress at the bonding position to the maximum modal stress value of the blade. It should be noted here that the stress value obtained by blade modal calculation and solution is converted from strain normalization, and its absolute value has no practical significance. Its stress distribution is the same as the actual resonance state, so its stress is called modal stress. Similarly, its deformation is modal deformation. The sensitivities of single-crystal blades with different crystal orientation angles are different. Therefore, it is necessary to first perform a simulation analysis of the blade vibration stress according to the actual crystal orientation angle of the measured blade, calculate the sensitivity of the strain gauge bonding position corresponding to the measured maximum vibration stress at the actual crystal orientation angle, and then convert the actual maximum vibration stress of the measured blade according to this sensitivity and the measured maximum vibration stress.

[0026] In addition, for single-crystal blades, the influence of the angle between the crystal coordinate axis and the overall coordinate axis of the model on the vibration stress needs to be considered. The analysis needs to cover all possible crystal orientation angles of the actual machined parts. Generally speaking, for aero-engine single-crystal blades, the crystal main axis angle of the blade in the spanwise

[001] direction is generally controlled θ while releasing the crystal angle directions of the other two crystal main axes ω 、 κ . The crystal orientation angles of single-crystal blades are as Figure 3 shown. Assuming that the maximum allowable angle in the

[001] direction is the crystal angle θ max , then θ 、 ω 、 κ at least need to consider θ =0°, θ max , ω =0°, 90°, 180°, 270°, κ= 0°, where, due to the large amount of calculation and long calculation period for vibration stress analysis, it is impossible to obtain the maximum vibration stress parameters through a large number of calculations within a reasonable time. Therefore, it is necessary to cut the crystal orientation angles for calculation and analysis. These several angle values above are the reference angles recommended by the airworthiness authorities. Of course, in other embodiments of the present invention, the angle values can be selected according to actual needs to design different crystal orientation angles. Since the vibration stress magnitudes of different crystal orientation angles are different, the actual maximum vibration stress of the measured blade is generally not the maximum vibration stress of the blade at other crystal orientation angles. Therefore, through the finite element modal analysis of the blade, the present invention can obtain the stress coefficients of the blade at different crystal orientation angles, screen out the maximum stress coefficient, and then convert the maximum vibration stress of the blade according to the maximum stress coefficient, that is, the maximum value among the vibration stresses corresponding to all crystal orientation angles.

[0027] Among them, according to the finite element principle, for single crystal blades and directionally solidified blades with such a linear elastic model, at a certain crystal orientation angle the stress coefficient can be expressed as: , where represents the maximum modal stress value corresponding to the m-th order modal obtained by finite element modal calculation, represents the maximum modal displacement value corresponding to the m-th order modal obtained by finite element modal calculation. Therefore, only by obtaining and through finite element modal calculation, the can be solved.

[0028] Among them, the specific finite element modal calculation process is as follows: First, establish a finite element model of the blade. During the model establishment process, it is necessary to make the z-direction of the global Cartesian coordinate system coincide with the blade span direction, so that the established model material coordinate system (i.e., the crystal coordinate system) completely coincides with the global coordinate system. At this time, the crystal orientation angle is =[0,0,0]; then, define and establish a local coordinate system through the command stream (CLOCAL, cnum, 0, 0, 0, 0, thxy, thyz, thzx). cnum is the number of the local coordinate system, and the relationship between the coordinate rotation angle and the crystal orientation angle is one-to-one as [thxy, thyz, thzx] → . This command stream can realize the local coordinate system x’y’z’ of the single crystal micro-element crystal from the initial coincidence with the global Cartesian coordinate system xyz , through rotation around the y’ axis by → rotation around the z’ axis by → rotation around the x’ axis by , the local coordinate system can be rotated to the given crystal orientation angle; then through the following command stream: CSYS, cnum; ESYS, ALL, cnum, where CSYS is to activate the local coordinate system and ESYS is to set the element coordinate system, rotate the element coordinate system of the blade finite element model (by default, the element coordinate system is parallel to the global Cartesian coordinate system, and the single crystal main axis direction is determined by the element coordinate system) to the local coordinate system. Since the model does not rotate with the element coordinate system, the included angles between the overall coordinate system of the model and the main axes of the single crystal are the same as the given crystal orientation direction; finally, conduct the blade vibration characteristic analysis on the rotated finite element model, and read the maximum modal displacement of the model results at a certain crystal orientation angle and the maximum modal stress , the stress coefficient of the blade at the crystal orientation angle can be calculated .

[0029] After calculating the stress coefficients of the blade at different crystal orientation angles , screen out the maximum stress coefficients at all crystal orientation angles , where represents the crystal orientation angle corresponding to the maximum stress coefficient. Since the vibration stress is proportional to the stress coefficient, the maximum stress coefficient corresponds to the maximum vibration stress at all crystal orientation angles. Then, based on the maximum stress coefficient , the stress coefficient at the actual crystal orientation angle and the actual maximum vibration stress of the measured blade , the maximum vibration stress of the blade at different crystal orientation angles can be calculated. Among them, the maximum vibration stress of the blade at all crystal orientation angles is specifically calculated based on the following formula: ; where represents the actual crystal orientation angle of the measured blade, represents the crystal orientation angle corresponding to the maximum stress coefficient, represents the stress coefficient at the actual crystal orientation angle, represents the maximum stress coefficient, represents the actual maximum vibration stress of the measured blade, represents the maximum vibration stress of the blade at all crystal orientation angles. After calculating the maximum vibration stress of the blade at all crystal orientation angles, it is taken as the maximum vibration stress in statistics.

[0030] It can be understood that the present invention constructs a conversion path between the vibration stress of the measured blade and the vibration stress at different crystal orientation angles through the stress coefficient formula and the maximum vibration stress conversion formula, which can more precisely consider the change in vibration stress caused by the change in the crystal orientation angle of the blade. Therefore, high-cycle fatigue analysis can be more accurately carried out on single-crystal blades and directionally solidified blades, and the confidence level of the high-cycle fatigue analysis results of the blades is higher, reducing the risk of high-cycle fatigue failure.

[0031] In addition, in the step S5, based on the maximum vibration stress in statistics and the steady-state stress value obtained in the step S1, for each vibration order, the corresponding GoodMan diagram can be drawn, and high-cycle fatigue analysis of the blade can be carried out. Among them, the drawing principle of the GoodMan diagram and the process of high-cycle fatigue analysis of the blade belong to the prior art and will not be elaborated here.

[0032] Optionally, after calculating the maximum vibration stress at all crystal orientation angles in the step S4, the following content is further included: Considering the influence of the crystal orientation angle on the exciting force, the maximum vibration stress at all crystal orientation angles is corrected.

[0033] Specifically, the process of considering the influence of the crystal orientation angle on the exciting force and correcting the maximum vibration stress at all crystal orientation angles includes the following content: First, based on the typical curve of the blade aerodynamic pressure varying with the rotational speed, the aerodynamic pressures of the blade at the actual crystal orientation angle and the crystal orientation angle corresponding to the maximum stress coefficient are respectively obtained through rotational speed interpolation. Then, according to the principle that the aerodynamic exciting force is proportional to the aerodynamic pressure, the maximum vibration stress at all crystal orientation angles is corrected by using the ratio of the aerodynamic pressures at the actual crystal orientation angle and the crystal orientation angle corresponding to the maximum stress coefficient, and the corrected maximum vibration stress is calculated.

[0034] It can be understood that since the maximum vibration stress corresponds to the resonance state of the blade, the present invention also considers the change in the exciting force during the resonance state of the blade caused by the change in the resonance frequency corresponding to different crystal orientation angles. According to engineering experience, the aerodynamic exciting force of the same blade is in a proportional relationship with the blade aerodynamic pressure. Therefore, the present invention first based on the typical curve of the blade aerodynamic pressure varying with the rotational speed, such as Figure 4 shown, the aerodynamic pressures corresponding to a0 and can be respectively obtained by rotational speed interpolation. Specifically, the resonance frequency at the crystal orientation angle corresponding to the maximum stress coefficient is obtained through finite element analysis, and is converted to the corresponding resonance rotational speed (i.e., the engine operating rotational speed when resonance can be excited), and the conversion formula is: , and then the functional relationship P(n) is obtained by fitting the typical curve of the blade aerodynamic pressure varying with the rotational speed, and Substitute into P(n), and the crystal orientation angle can be calculated. The pneumatic pressure in the lower resonance state ; Similarly, the pneumatic pressure in the lower resonance state at the actual crystal orientation angle a0 can also be calculated. .

[0035] Then, according to the principle that the pneumatic excitation force is proportional to the pneumatic pressure, the ratio of the pneumatic excitation forces at the crystal orientation angle and the crystal orientation angle a0 can be converted, and the conversion formula is: , where represents the pneumatic excitation force in the lower resonance state at the crystal orientation angle , represents the pneumatic excitation force in the lower resonance state at the crystal orientation angle a0.

[0036] Since the response in the same-order resonance state is proportional to the pneumatic excitation force, therefore, the maximum vibration stress at all crystal orientation angles can be corrected based on the following formula, and the corrected maximum vibration stress is used as the maximum vibration stress in statistics: .

[0037] It can be understood that the present invention also first considers the change in the excitation force caused by different crystal orientation angles, which in turn causes the change in the maximum vibration stress of the blade, and can more accurately consider the maximum vibration stress in actual engineering applications, guiding the high-cycle fatigue design of the blade to reduce the risk of high-cycle fatigue failure of the blade.

[0038] In addition, as Figure 5 shown, another embodiment of the present invention also provides a blade high-cycle fatigue analysis system considering the influence of crystal orientation angles, preferably adopting the blade high-cycle fatigue analysis method considering the influence of crystal orientation angles as described above, including: A steady-state stress analysis module for obtaining the steady-state stress of the blade in various working states; A vibration characteristic analysis module for obtaining the vibration characteristics of the blade in various working states and performing vibration characteristic analysis to determine the vibration order concerned with the high-cycle fatigue of the blade; A dynamic stress measurement and crystal orientation analysis module for performing a dynamic stress measurement test on the blade to obtain the measured maximum vibration stress of the measured blade and obtaining the actual crystal orientation angle of the measured blade; A maximum vibration stress calculation module for calculating the maximum vibration stress at all crystal orientation angles based on the measured maximum vibration stress and the actual crystal orientation angle of the measured blade for each vibration order, and using it as the maximum vibration stress in statistics; A blade high-cycle fatigue analysis module for drawing a GoodMan diagram based on the maximum vibration stress in statistics and the steady-state stress to conduct blade high-cycle fatigue analysis.

[0039] It can be understood that the blade high-cycle fatigue analysis system considering the influence of crystal orientation angle in this embodiment calculates the maximum vibration stress at all crystal orientation angles based on the measured maximum vibration stress and the actual crystal orientation angle of the blade to be measured, and takes it as the maximum vibration stress in statistics, fully considering the influence of the crystal angle change on the maximum vibration stress in statistics. The confidence level of the obtained blade high-cycle fatigue assessment result is higher, and the high-cycle fatigue analysis of the blade can be carried out more precisely, reducing the risk of blade high-cycle fatigue failure.

[0040] In addition, the maximum vibration stress calculation module is further configured to correct the maximum vibration stress at all crystal orientation angles after calculating the maximum vibration stress at all crystal orientation angles, considering the influence of the crystal orientation angle on the exciting force.

[0041] It can be understood that each module in the system embodiment of the present invention corresponds to each step in the above method embodiment. Therefore, the specific working principle of each module will not be elaborated here, and reference can be made to each step in the above method embodiment correspondingly.

[0042] In addition, another embodiment of the present invention further provides an electronic device, including a processor and a memory. A computer program is stored in the memory, and the processor is configured to execute the steps of the method as described above by calling the computer program stored in the memory.

[0043] In addition, another embodiment of the present invention further provides a computer-readable storage medium for storing a computer program for blade high-cycle fatigue analysis considering the influence of crystal orientation angle. When the computer program runs on a computer, it executes the steps of the method as described above.

[0044] The forms of computer-readable storage media generally include: floppy disks, flexible disks, hard disks, magnetic tapes, any other magnetic media, CD-ROMs, any other optical media, punch cards, paper tapes, any other physical media with a pattern of holes, random access memories (RAMs), programmable read-only memories (PROMs), erasable programmable read-only memories (EPROMs), flash erasable programmable read-only memories (FLASH-EPROMs), any other memory chips or cartridges, or any other media readable by a computer. Instructions can further be transmitted or received by a transmission medium. The term transmission medium can include any tangible or intangible medium that can be used to store, encode, or carry instructions for execution by a machine, and includes digital or analog communication signals or the intangible media that facilitate the communication of the above instructions. Transmission media includes coaxial cables, copper wires, and optical fibers, which include the wires of a bus used to transmit a computer data signal.

[0045] Those skilled in the art will understand that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk memories, CD-ROMs, optical memories, etc.) containing computer-usable program code. The solutions in the embodiments of the present application can be implemented in various computer languages. For example, object-oriented programming languages such as Java and interpreted scripting languages such as JavaScript.

[0046] The present application is described with reference to the flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each flow and / or block in the flowcharts and / or block diagrams, and the combination of flows and / or blocks in the flowcharts and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to the processors of general-purpose computers, special-purpose computers, embedded processors, or other programmable data processing devices to generate a machine, such that the instructions executed by the processors of the computer or other programmable data processing devices generate means for implementing the functions specified in Figure 1 one flow or multiple flows and / or blocks Figure 1 one block or multiple blocks.

[0047] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to operate in a particular manner, such that the instructions stored in the computer-readable memory produce a manufacture including an instruction device that implements the functions specified in one or more of the procedures Figure 1 one or more of the procedures and / or blocks Figure 1 specified in one or more of the blocks or blocks.

[0048] These computer program instructions can also be loaded onto a computer or other programmable data processing apparatus, such that a series of operational steps are performed on the computer or other programmable apparatus to produce a computer-implemented process, whereby the instructions executed on the computer or other programmable apparatus provide steps for implementing the functions specified in one or more of the procedures Figure 1 one or more of the procedures and / or blocks Figure 1 specified in one or more of the blocks or blocks.

[0049] Although the preferred embodiments of the present application have been described, those skilled in the art can make additional changes and modifications once they learn the basic creative concept. Therefore, the appended claims are intended to be construed to include the preferred embodiments as well as all changes and modifications that fall within the scope of the present application.

[0050] Obviously, those skilled in the art can make various changes and modifications to the present application without departing from the spirit and scope of the present application. Thus, if these modifications and variations of the present application fall within the scope of the claims of the present application and their equivalent technologies, the present application is also intended to include these modifications and variations.

[0051] The above are only the preferred embodiments of the present invention and are not intended to limit the present invention. For those skilled in the art, the present invention can have various changes and variations. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present invention shall be included within the protection scope of the present invention.

Claims

1. A blade high cycle fatigue analysis method considering the influence of crystal orientation angle, characterized in that: Includes the following: Obtain the steady-state stress of the blade under various working conditions; Obtain the vibration characteristics of the blades under various working conditions, and perform vibration characteristic analysis to determine the vibration order of concern for high-cycle fatigue of the blades; Conduct a dynamic stress measurement test on the blade to obtain the maximum vibration stress of the measured blade and the actual crystal orientation angle of the measured blade; For each vibration order, the maximum vibration stress at all crystal orientation angles is calculated based on the measured maximum vibration stress of the measured blade and the actual crystal orientation angle, and is used as the statistical maximum vibration stress; Based on the statistical maximum vibration stress and steady-state stress, a GoodMan diagram is drawn to carry out high-cycle fatigue analysis of the blade.

2. The blade high cycle fatigue analysis method considering the influence of crystal orientation angle according to claim 1 is characterized in that: The process of calculating the maximum vibration stress at all crystal orientation angles according to the measured maximum vibration stress of the measured blade and the actual crystal orientation angle comprises the following contents: According to the actual crystal orientation angle of the measured blade, simulation analysis is performed to obtain the sensitivity of the strain gauge patch position corresponding to the measured maximum vibration stress under the actual crystal orientation angle, and the actual maximum vibration stress of the measured blade is obtained by conversion based on the sensitivity and the measured maximum vibration stress; Conduct finite element modal analysis of the blade to obtain the stress coefficient of the blade at different crystal orientation angles and screen out the maximum stress coefficient; The maximum vibration stress of the blade at all crystal orientation angles is calculated based on the maximum stress coefficient, the stress coefficient at the actual crystal orientation angle and the actual maximum vibration stress of the measured blade.

3. The blade high cycle fatigue analysis method considering the influence of crystal orientation angle according to claim 2 is characterized in that: The maximum vibration stress of the blade at all crystal orientation angles is calculated based on the following formula: in, Indicates the actual crystal angle of the blade being measured. represents the crystal orientation angle corresponding to the maximum stress coefficient, represents the stress coefficient at the actual crystal orientation angle, represents the maximum stress coefficient, Indicates the actual maximum vibration stress of the blade being tested, It represents the maximum vibration stress of the blade at all crystal orientation angles.

4. The blade high cycle fatigue analysis method considering the influence of crystal orientation angle according to claim 2 is characterized in that: After calculating the maximum vibration stress at all crystal orientation angles, the following is also included: Considering the influence of crystal orientation angle on the exciting force, the maximum vibration stress at all crystal orientation angles is corrected.

5. The blade high cycle fatigue analysis method considering the influence of crystal orientation angle according to claim 4 is characterized in that: The process of considering the influence of the crystal orientation angle on the exciting force and correcting the maximum vibration stress at all crystal orientation angles includes the following: First, based on the typical curve of the blade aerodynamic pressure changing with the rotational speed, the aerodynamic pressure of the blade at the actual crystal orientation angle and the crystal orientation angle corresponding to the maximum stress coefficient are obtained by rotational speed interpolation. Then, according to the principle that the aerodynamic exciting force is proportional to the aerodynamic pressure, the maximum vibration stress at all crystal orientation angles is corrected by using the ratio of the aerodynamic pressure at the actual crystal orientation angle and the crystal orientation angle corresponding to the maximum stress coefficient, and the corrected maximum vibration stress is calculated.

6. The blade high cycle fatigue analysis method considering the influence of crystal orientation angle according to claim 5, characterized in that: The corrected maximum vibration stress is calculated based on the following formula: in, Indicates the actual crystal angle of the blade being measured. represents the crystal orientation angle corresponding to the maximum stress coefficient, represents the stress coefficient at the actual crystal orientation angle, represents the maximum stress coefficient, Indicates the actual maximum vibration stress of the blade being tested, represents the maximum vibration stress of the blade at all crystal orientation angles, Indicates the crystal orientation angle The aerodynamic pressure in the lower resonance state, Indicates the aerodynamic pressure in the resonant state at the actual crystal orientation angle a0.

7. A blade high cycle fatigue analysis system considering the influence of crystal orientation angle, characterized in that: include: Steady-state stress analysis module, used to obtain the steady-state stress of the blade under various working conditions; The vibration characteristics analysis module is used to obtain the vibration characteristics of the blade under various working conditions, and to perform vibration characteristics analysis to determine the vibration order of concern for high-cycle fatigue of the blade; The dynamic stress measurement and crystal orientation analysis module is used to perform a dynamic stress measurement test on the blade, obtain the maximum vibration stress of the measured blade, and obtain the actual crystal orientation angle of the measured blade; A maximum vibration stress calculation module is used to calculate the maximum vibration stress under all crystal direction angles for each vibration order according to the measured maximum vibration stress of the measured blade and the actual crystal direction angle, and use it as the statistical maximum vibration stress; The blade high cycle fatigue analysis module is used to draw the GoodMan diagram based on the statistical maximum vibration stress and steady-state stress to carry out blade high cycle fatigue analysis.

8. The blade high cycle fatigue analysis system considering the influence of crystal orientation angle according to claim 7, characterized in that: The maximum vibration stress calculation module is also used to correct the maximum vibration stress at all crystal orientation angles after calculating the maximum vibration stress at all crystal orientation angles, taking into account the influence of the crystal orientation angle on the exciting force.

9. An electronic device, characterized in that: The method comprises a processor and a memory, wherein the memory stores a computer program, and the processor executes the steps of the method according to any one of claims 1 to 6 by calling the computer program stored in the memory.

10. A computer-readable storage medium for storing a computer program for performing high-cycle fatigue analysis of blades taking into account the influence of crystal orientation angles, characterized in that: When the computer program is run on a computer, the computer program executes the steps of the method according to any one of claims 1 to 6.

Images