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

Abstract

The present invention discloses a blade high-cycle fatigue analysis method and system that considers the influence of crystal orientation angle. The blade high-cycle fatigue analysis method of the present invention 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 uses it as the statistically maximum vibration stress. The influence of the crystal angle change on the statistically maximum vibration stress is fully considered. The confidence level of the blade high-cycle fatigue assessment result obtained by analysis is higher, and high-cycle fatigue analysis of the blade can be carried out more accurately, reducing the risk of high-cycle fatigue failure of the blade.

Description

Technical Field

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

[0002] Aircraft engine blades typically utilize single crystal or directionally crystallized blades. Single crystal blades offer superior performance to directionally crystallized blades. To ensure engine safety, high-cycle fatigue (HCF) assessments are necessary for blades. This is also a key component of the Aircraft Engine Airworthiness Regulations. Currently, a conventional method for evaluating blade HCF is based on vibration strain test data. This method directly converts the maximum strain corresponding to each order obtained from the blade test into a vibration stress value, which is then used as the statistically significant maximum vibration stress value for HCF assessment. However, the maximum vibration stress must account for the effects of different crystal orientations, operating conditions, and machining errors. Crystal orientation, in particular, has a significant impact on vibration stress. When the actual crystal orientation angle of the blade differs significantly from the theoretical one, the maximum vibration stress obtained by existing methods may be significantly lower than the actual maximum stress. This means that if the HCF margin obtained using existing methods is small, the vibration stress in actual engine blades may exceed the allowable value. Therefore, existing blade HCF analysis methods fail to consider the influence of crystal orientation angle, resulting in poor accuracy. Summary of the Invention

[0003] The present invention provides a blade high-cycle fatigue analysis method and system, electronic equipment, and computer-readable storage medium that take into account the influence of crystal orientation angle. It fully considers the influence of crystal angle changes on the statistically maximum vibration stress. The confidence level of the blade high-cycle fatigue assessment results obtained by analysis is higher, and high-cycle fatigue analysis of the blade can be performed more accurately, reducing the risk of blade high-cycle fatigue failure.

[0004] According to one aspect of the present invention, a blade high cycle fatigue analysis method considering the influence of crystal orientation angle is provided, comprising the following contents:

[0005] Obtain the steady-state stress of the blade under various working conditions;

[0006] Obtain the vibration characteristics of the blade under various working conditions, and perform vibration characteristic analysis to determine the vibration order of concern for high-cycle fatigue of the blade;

[0007] Perform a dynamic stress measurement test on the blade to obtain the maximum vibration stress of the blade under test and the actual crystal orientation angle of the blade under test;

[0008] For each vibration order, the maximum vibration stress at all crystal orientation angles is calculated based on the measured maximum vibration stress and the actual crystal orientation angle of the blade under test, and is used as the statistically maximum vibration stress;

[0009] 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.

[0010] Furthermore, the process of 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 includes the following steps:

[0011] A simulation analysis is performed based on the actual crystal orientation angle of the blade being measured to obtain the sensitivity of the strain gauge patch position corresponding to the measured maximum vibration stress at the actual crystal orientation angle. The actual maximum vibration stress of the blade being measured is then converted based on the sensitivity and the measured maximum vibration stress.

[0012] 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;

[0013] 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.

[0014] Furthermore, the maximum vibration stress of the blade at all crystal orientation angles is calculated based on the following formula:

[0015] ;

[0016] in, Indicates the actual crystal orientation 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.

[0017] Furthermore, after calculating the maximum vibration stress at all crystal orientation angles, the following is also included:

[0018] Considering the influence of crystal orientation angle on the exciting force, the maximum vibration stress at all crystal orientation angles is corrected.

[0019] Furthermore, the process of correcting the maximum vibration stress at all crystal orientation angles by considering the influence of the crystal orientation angle on the exciting force includes the following steps:

[0020] 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 is obtained by speed interpolation. Then, according to the principle that the aerodynamic excitation force is proportional to the aerodynamic pressure, the maximum vibration stress at all crystal orientation angles is corrected 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.

[0021] Furthermore, the corrected maximum vibration stress is calculated based on the following formula:

[0022] ;

[0023] in, Indicates the actual crystal orientation 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.

[0024] In addition, the present invention also provides a blade high cycle fatigue analysis system that considers the influence of crystal orientation angle, comprising:

[0025] Steady-state stress analysis module, used to obtain the steady-state stress of the blade under various working conditions;

[0026] The vibration characteristics analysis module is used to obtain the vibration characteristics of the blade under various working conditions and perform vibration characteristics analysis to determine the vibration order of concern for high-cycle fatigue of the blade;

[0027] Dynamic stress measurement and crystal orientation analysis module, used to perform dynamic stress measurement tests on blades, obtain the maximum vibration stress of the blades under test, and obtain the actual crystal orientation angle of the blades under test;

[0028] A maximum vibration stress calculation module is used to calculate the maximum vibration stress at all crystal orientation angles for each vibration order based on the measured maximum vibration stress and the actual crystal orientation angle of the blade under test, and use it as the statistical maximum vibration stress;

[0029] The blade high-cycle fatigue analysis module is used to draw GoodMan diagrams based on statistically maximum vibration stress and steady-state stress to conduct blade high-cycle fatigue analysis.

[0030] Furthermore, the maximum vibration stress calculation module is further configured to, after calculating the maximum vibration stress at all crystal orientation angles, correct the maximum vibration stress at all crystal orientation angles by taking into account the influence of the crystal orientation angle on the exciting force.

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

[0032] In addition, the present invention also provides a computer-readable storage medium for storing a computer program for performing high-cycle fatigue analysis of blades considering the influence of crystal orientation angles, wherein the computer program executes the steps of the method described above when running on a computer.

[0033] The present invention has the following beneficial effects:

[0034] The blade high-cycle fatigue analysis method considering the influence of crystal orientation angle of the present invention 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 uses it as the statistically maximum vibration stress. It fully considers the influence of the crystal angle change on the statistically maximum vibration stress. The confidence level of the blade high-cycle fatigue assessment result obtained by analysis is higher, and high-cycle fatigue analysis of the blade can be carried out more accurately, reducing the risk of high-cycle fatigue failure of the blade.

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

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

[0037] The accompanying drawings, which constitute part of this application, are intended to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are intended to explain the present invention and do not constitute an undue limitation of the present invention. In the accompanying drawings:

[0038] Figure 1 This is a flow chart of a blade high cycle fatigue analysis method considering the influence of crystal orientation angle according to a preferred embodiment of the present application;

[0039] Figure 2 yes Figure 1 Schematic diagram of the sub-process of step S4;

[0040] Figure 3 Schematic diagram of the crystal orientation angle of a single crystal blade in a preferred embodiment of the present application;

[0041] Figure 4 Schematic diagram of a typical curve showing the change of blade aerodynamic pressure with rotational speed in a preferred embodiment of the present application;

[0042] Figure 5 This is a schematic diagram of the module structure of a blade high-cycle fatigue analysis system that takes into account the influence of crystal orientation angles in another embodiment of the present application. DETAILED DESCRIPTION

[0043] It should be noted that, in the absence of conflict, the embodiments and features of the embodiments in this application can be combined with each other. The present application will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.

[0044] Reference Figure 1 The preferred embodiment of the present application provides a blade high cycle fatigue analysis method considering the influence of crystal orientation angle, including the following contents:

[0045] Step S1: Obtaining the steady-state stress of the blade under various working conditions;

[0046] Step S2: Obtain the vibration characteristics of the blade under various working conditions, and perform vibration characteristic analysis to determine the vibration order of concern for high-cycle fatigue of the blade;

[0047] Step S3: performing 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;

[0048] Step S4: For each vibration order, the maximum vibration stress at all crystal orientation angles is calculated based on the measured maximum vibration stress and the actual crystal orientation angle of the measured blade, and is used as the statistical maximum vibration stress;

[0049] Step S5: Draw a GoodMan diagram based on the statistical maximum vibration stress and steady-state stress to perform high-cycle fatigue analysis of the blade.

[0050] It can be understood that the blade high-cycle fatigue analysis method of this embodiment that considers 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 measured blade, and uses it as the statistically maximum vibration stress, fully considering the influence of the crystal angle change on the statistically maximum vibration stress. The confidence level of the blade high-cycle fatigue assessment results obtained by analysis is higher, and high-cycle fatigue analysis of the blade can be carried out more accurately, reducing the risk of high-cycle fatigue failure of the blade.

[0051] Specifically, in step S1, the steady-state stress of the blade under various actual operating conditions is first analyzed using the finite element method. The specific analysis process of the steady-state stress belongs to the prior art and is not further described here. In addition, the blades referred to in the present invention include single crystal blades and directionally crystallized blades. The following description uses single crystal blades as an example for illustrative purposes.

[0052] In step S2, the finite element method is used to analyze the vibration characteristics of the single crystal blade under various actual working conditions, and the vibration order of concern for high-cycle fatigue of the blade, that is, the dangerous order, is determined based on the vibration characteristics analysis. The specific analysis process belongs to the existing technology and will not be repeated here.

[0053] In step S3, since the actual blade vibration stress cannot currently be determined through computational analysis, a dynamic stress measurement test is conducted on the single-crystal blade according to the relevant specifications for "Blade Vibration Stress Measurement Test Point Design" and the "Blade Dynamic Stress Measurement Test Technical Requirements." The strain corresponding to each vibration order is measured by multiple strain gauges. The maximum strain for each vibration order is then screened to calculate the maximum vibration stress of the measured blade at each vibration order. Furthermore, prior to the test, the actual crystal orientation angle of the measured blade is obtained through crystal orientation analysis. The specific crystal orientation analysis process is prior art and will not be detailed here.

[0054] It can be understood that the high-cycle fatigue analysis of blades for airworthiness certification is to prove 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 believed that the blade will not fail due to high-cycle fatigue. The maximum vibration stress needs to take into account 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 specific blade tested to the maximum vibration stress value of all blades of the mass-produced engine under different service conditions through conversion or correction, that is, the maximum vibration stress value under all crystal orientation angles, and the statistically maximum vibration stress value can be obtained. Therefore, in the step S4, for each vibration order determined in step S2, the maximum vibration stress under all crystal orientation angles is calculated based on the measured maximum vibration stress and the actual crystal orientation angle of the blade being tested, and it is used as the statistically maximum vibration stress. Among them, if Figure 2 As shown, the process of 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 includes the following:

[0055] Step 41: Perform simulation analysis based on 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 actual maximum vibration stress of the blade under test based on the sensitivity and the measured maximum vibration stress;

[0056] Step 42: Conduct finite element modal analysis of the blade to obtain the stress coefficients of the blade at different crystal orientation angles, and select the maximum stress coefficient;

[0057] Step 43: Calculate the maximum vibration stress of the blade at all crystal orientation angles 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.

[0058] Specifically, due to the inherent properties of the strain gauge, in order to obtain high-precision blade strain values, when conducting dynamic stress measurement tests, the strain gauge patch position needs to be in a flat area with relatively low temperature, and avoid being attached to uneven positions such as fillets, chamfers, and holes. Excessive temperature will cause the strain gauge to fail. The above-mentioned area is generally not at the actual maximum vibration stress position of the blade being measured. Therefore, it is necessary to convert it to the actual maximum vibration stress value of the blade being measured based on the sensitivity, where the sensitivity is the ratio of the modal stress at the patch position to the maximum modal stress value of the blade. It should be noted here that the stress value solved by the blade modal calculation is converted from strain normalization, and its absolute value has no practical meaning. 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 sensitivity of single crystal blades with different crystal orientation angles is different. Therefore, it is necessary to first perform a simulation analysis of the blade vibration stress based on the actual crystal orientation angle of the blade being measured, calculate the sensitivity of the strain gauge patch position corresponding to the measured maximum vibration stress at the actual crystal orientation angle, and then convert the actual maximum vibration stress of the blade being measured based on this sensitivity and the measured maximum vibration stress.

[0059] In addition, the single crystal blade needs to consider the effect of the angle between the crystal coordinate axis and the overall coordinate axis of the model on the vibration stress. The analysis needs to cover all possible crystal orientation angles of the actual workpiece. Generally speaking, the single crystal blade of an aerospace engine generally controls the material crystal principal axis angle in the blade span direction

[001] . θ , while releasing the crystal angle direction of the other two crystal main axes ω 、 κ , the crystal orientation angle of the single crystal blade is as follows Figure 3 Assume that the maximum allowable angle of the

[001] direction is the crystal angle θ max ,but θ 、 ω 、 κ At least consider θ =0°, θ max , ω =0°, 90°, 180°, 270°, κ=0°], wherein, due to the large amount of calculation required for vibration stress analysis and the long calculation cycle, it is impossible to obtain the maximum vibration stress parameters through massive calculations within a reasonable time. Therefore, it is necessary to tailor the crystal orientation angles used in the calculation and analysis. The above-mentioned angle values are 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 at different crystal orientation angles varies, 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, the present invention can obtain the stress coefficients of the blade at different crystal orientation angles by performing finite element modal analysis on the blade, screen out the maximum stress coefficient, and then convert the maximum vibration stress of the blade based on the maximum stress coefficient to obtain the maximum vibration stress of the blade, that is, the maximum value of the vibration stress corresponding to all crystal orientation angles.

[0060] Among them, according to the finite element principle, for the linear elastic model of single crystal blades and directional crystallization blades, at a certain crystal orientation angle Stress coefficient under It can be expressed as: ,in, Indicates the maximum modal stress value corresponding to the mth order mode obtained by finite element modal calculation, Indicates the maximum modal displacement value corresponding to the mth order mode obtained by finite element modal calculation. Therefore, it is only necessary to obtain and , we can solve it and get .

[0061] The specific finite element modal calculation process is as follows: first establish the blade finite element model. During the model establishment process, the z direction of the global Cartesian coordinate system needs to coincide with the blade span direction, so that the established model material coordinate system (i.e., crystal coordinate system) completely coincides with the global coordinate system. At this time, the crystal direction angle is [ ]=[0,0,0]; Then, define and establish the 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 one-to-one correspondence between the coordinate rotation angle and the crystal orientation angle is [thxy,thyz,thzx]→[ ], this command stream can realize the local coordinate system of single crystal microelement crystal x'y'z' From the initial and global Cartesian coordinate systems xyz Overlap, after winding y’ Axis rotation →Wrap z’ Axis rotation →Wrap x’ Axis rotation , the local coordinate system can be rotated to the given crystal orientation angle; then the following command flow is used: CSYS,cnum; ESYS,ALL,cnum, where CSYS is used to activate the local coordinate system and ESYS is used to set the unit coordinate system. The unit coordinate system of the blade finite element model (under default conditions, the unit coordinate system is parallel to the global Cartesian coordinate system, and the direction of the main axis of the single crystal is determined by the unit coordinate system) is rotated to the local coordinate system. Since the model will not rotate with the unit coordinate system, the angle between the overall coordinate system of the model and the main axis of the single crystal is the same as the given crystal orientation direction; finally, the blade vibration characteristics analysis is carried out on the rotated finite element model, and the maximum modal displacement of the model result of a certain crystal orientation angle is read. and maximum modal stress , you can calculate the blade's crystal orientation angle Stress coefficient under .

[0062] The stress coefficient of the blade at different crystal orientation angles is calculated Finally, the maximum stress coefficient at all crystal orientation angles is screened out. ,in, Indicates the crystal orientation angle corresponding to the maximum stress coefficient. Since the vibration stress is proportional to the stress coefficient, the maximum stress coefficient That is, it corresponds to the maximum vibration stress at all crystal orientation angles, which can be calculated based on the maximum stress coefficient , stress coefficient at actual crystal orientation angle and the actual maximum vibration stress of the blade under test , the maximum vibration stress of the blade under different crystal orientation angles can be calculated The maximum vibration stress of the blade at all crystal orientation angles is calculated based on the following formula:

[0063] ;

[0064] in, Indicates the actual crystal orientation 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, Indicates the maximum vibration stress of the blade at all crystal orientation angles. Then, it is taken as the statistical maximum vibration stress.

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

[0066] Furthermore, in step S5, based on the statistically determined maximum vibration stress and the steady-state stress value obtained in step S1, a corresponding GoodMan diagram can be drawn for each vibration order, and high-cycle fatigue analysis of the blade can be performed. The principles for drawing the GoodMan diagram and the blade high-cycle fatigue analysis process are prior art and will not be further elaborated here.

[0067] Optionally, after calculating the maximum vibration stress at all crystal orientation angles, step S4 further includes the following:

[0068] Considering the influence of crystal orientation angle on the exciting force, the maximum vibration stress at all crystal orientation angles is corrected.

[0069] Specifically, the process of correcting the maximum vibration stress at all crystal orientation angles by considering the influence of the crystal orientation angle on the exciting force includes the following steps:

[0070] 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 is obtained by speed interpolation. Then, according to the principle that the aerodynamic excitation force is proportional to the aerodynamic pressure, the maximum vibration stress at all crystal orientation angles is corrected 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.

[0071] It is understandable that since the maximum vibration stress corresponds to the blade resonance state, the present invention also considers the change in the excitation force when the blade is in the resonance state due to the change in the resonant frequency corresponding to different crystal orientation angles. According to engineering experience, the aerodynamic excitation force of the same blade is directly proportional to the aerodynamic pressure of the blade. Therefore, the present invention is based on the typical curve of the blade aerodynamic pressure changing with the rotational speed, such as Figure 4 As shown, a0 and a1 can be obtained by speed interpolation. Specifically, the crystal orientation angle corresponding to the maximum stress coefficient is obtained through finite element analysis. The resonant frequency , converted to the corresponding resonant speed according to the excitation frequency N (i.e. the engine operating speed when resonance can be aroused), the conversion formula is: , and then the function relationship P(n) is obtained by fitting the typical curve of blade aerodynamic pressure changing with speed. Substituting into P(n), the crystal orientation angle can be calculated Aerodynamic pressure in the lower resonance state Similarly, the aerodynamic pressure in the resonance state under the actual crystal orientation angle a0 can also be calculated .

[0072] Then, according to the principle that the pneumatic excitation force is proportional to the pneumatic pressure, the crystal orientation angle can be converted to The ratio of the aerodynamic excitation force under the crystal orientation angle a0 is calculated as follows: ,in, Indicates the crystal orientation angle The aerodynamic excitation force in the lower resonance state, It represents the aerodynamic excitation force in the resonant state at the crystal orientation angle a0.

[0073] Since the response under the same order resonance state is proportional to the aerodynamic excitation force, the maximum vibration stress under all crystal orientation angles can be corrected based on the following formula, and the corrected maximum vibration stress is then used as the statistical maximum vibration stress:

[0074] .

[0075] It can be understood that the present invention also takes into account for the first time the changes in exciting forces caused by different crystal orientation angles and thus the changes in the maximum vibration stress of the blade. It can more accurately consider the maximum vibration stress in actual engineering applications, guide the high-cycle fatigue design of the blade to reduce the risk of high-cycle fatigue failure of the blade.

[0076] In addition, if Figure 5 As shown, another embodiment of the present invention further provides a blade high cycle fatigue analysis system considering the influence of crystal orientation angle, preferably using the blade high cycle fatigue analysis method considering the influence of crystal orientation angle as described above, including:

[0077] Steady-state stress analysis module, used to obtain the steady-state stress of the blade under various working conditions;

[0078] The vibration characteristics analysis module is used to obtain the vibration characteristics of the blade under various working conditions and perform vibration characteristics analysis to determine the vibration order of concern for high-cycle fatigue of the blade;

[0079] Dynamic stress measurement and crystal orientation analysis module, used to perform dynamic stress measurement tests on blades, obtain the maximum vibration stress of the blades under test, and obtain the actual crystal orientation angle of the blades under test;

[0080] A maximum vibration stress calculation module is used to calculate the maximum vibration stress at all crystal orientation angles for each vibration order based on the measured maximum vibration stress and the actual crystal orientation angle of the blade under test, and use it as the statistical maximum vibration stress;

[0081] The blade high-cycle fatigue analysis module is used to draw GoodMan diagrams based on statistically maximum vibration stress and steady-state stress to conduct blade high-cycle fatigue analysis.

[0082] It can be understood that the blade high-cycle fatigue analysis system of this embodiment, which takes into account 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 measured blade, and uses it as the statistically maximum vibration stress. It fully considers the influence of crystal angle changes on the statistically maximum vibration stress. The confidence level of the blade high-cycle fatigue assessment results obtained by analysis is higher, and high-cycle fatigue analysis of the blade can be carried out more accurately, reducing the risk of high-cycle fatigue failure of the blade.

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

[0084] It can be understood that the various modules of the system embodiment correspond to the various steps of the above method embodiment, so the specific working principles of each module will not be repeated here, and the corresponding references can be made to the various steps of the above method embodiment.

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

[0086] In addition, another embodiment of the present invention further provides a computer-readable storage medium for storing a computer program for performing high-cycle fatigue analysis of blades considering the influence of crystal orientation angles, wherein the computer program executes the steps of the method described above when running on a computer.

[0087] Common forms of computer-readable storage media include floppy disks, flexible disks, hard disks, magnetic tape, any other magnetic media, CD-ROMs, any other optical media, punch cards, paper tape, any other physical medium with a pattern of holes, random access memory (RAM), programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), flash-EPROM, any other memory chip or cartridge, or any other medium that can be read by a computer. Instructions can further be transmitted or received via a transmission medium. The term transmission medium may include any tangible or intangible medium that can be used to store, encode, or carry instructions for execution by a machine, including digital or analog communication signals or other intangible media that facilitate communication of such instructions. Transmission media include coaxial cables, copper wire, and fiber optics, including the wires of a bus used to transmit computer data signals.

[0088] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may take the form of a fully hardware embodiment, a fully software embodiment, or an embodiment combining software and hardware. Furthermore, the present application may 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 the embodiments of the present application may be implemented in various computer languages, such as the object-oriented programming language Java and the interpreted scripting language JavaScript.

[0089] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0090] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 a process or multiple processes and / or boxes Figure 1 The function specified in one or more boxes.

[0091] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0092] Although the preferred embodiments of the present application have been described, those skilled in the art may make additional changes and modifications to these embodiments once they have learned the basic creative concept. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments and all changes and modifications that fall within the scope of the present application.

[0093] Obviously, those skilled in the art may make various changes and modifications to this application without departing from the spirit and scope of this application. Thus, if these modifications and variations of this application fall within the scope of the claims of this application and their equivalents, this application is intended to include these modifications and variations.

[0094] The foregoing description is merely a preferred embodiment of the present invention and is not intended to limit the present invention. Those skilled in the art will readily appreciate that various modifications and variations of the present invention are possible. Any modifications, equivalent substitutions, or improvements made within the spirit and principles of the present invention are intended to be within the scope of protection 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 blade under various working conditions, and perform vibration characteristic analysis to determine the vibration order of concern for high-cycle fatigue of the blade; Perform a dynamic stress measurement test on the blade to obtain the maximum vibration stress of the blade under test and the actual crystal orientation angle of the blade under test; For each vibration order, the maximum vibration stress at all crystal orientation angles is calculated based on the measured maximum vibration stress and the actual crystal orientation angle of the blade under test, and is used as the statistically maximum vibration stress; Based on the statistical maximum vibration stress and steady-state stress, a GoodMan diagram is drawn to conduct high-cycle fatigue analysis of the blade; The process of 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 includes the following steps: A simulation analysis is performed based on the actual crystal orientation angle of the blade being measured to obtain the sensitivity of the strain gauge patch position corresponding to the measured maximum vibration stress at the actual crystal orientation angle. The actual maximum vibration stress of the blade being measured is then converted 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.

2. The blade high cycle fatigue analysis method considering the influence of crystal orientation angle according to claim 1 is characterized in that: The maximum vibration stress of the blade at all crystal orientation angles is calculated based on the following formula: Among them, a0 represents the actual crystal orientation angle of the blade being measured, a max 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.

3. The blade high cycle fatigue analysis method considering the influence of crystal orientation angle according to claim 1 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.

4. The blade high cycle fatigue analysis method considering the influence of crystal orientation angle according to claim 3 is characterized in that: The process of correcting the maximum vibration stress at all crystal orientation angles by considering the influence of the crystal orientation angle on the exciting force 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 is obtained by speed interpolation. Then, according to the principle that the aerodynamic excitation force is proportional to the aerodynamic pressure, the maximum vibration stress at all crystal orientation angles is corrected 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.

5. The blade high cycle fatigue analysis method considering the influence of crystal orientation angle according to claim 4 is characterized in that: The corrected maximum vibration stress is calculated based on the following formula: Among them, a0 represents the actual crystal orientation angle of the blade being measured, a max 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 a max The aerodynamic pressure in the lower resonance state, Indicates the aerodynamic pressure in the resonant state at the actual crystal orientation angle a0.

6. A blade high cycle fatigue analysis system considering the influence of crystal orientation angle, using the blade high cycle fatigue analysis method considering the influence of crystal orientation angle according to any one of claims 1 to 5, 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 perform vibration characteristics analysis to determine the vibration order of concern for high-cycle fatigue of the blade; Dynamic stress measurement and crystal orientation analysis module, used to perform dynamic stress measurement tests on blades, obtain the maximum vibration stress of the blades under test, and obtain the actual crystal orientation angle of the blades under test; A maximum vibration stress calculation module is used to calculate the maximum vibration stress at all crystal orientation angles for each vibration order based on the measured maximum vibration stress and the actual crystal orientation angle of the blade under test, and use it as the statistical maximum vibration stress; The blade high-cycle fatigue analysis module is used to draw GoodMan diagrams based on statistically maximum vibration stress and steady-state stress to conduct blade high-cycle fatigue analysis.

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

8. An electronic device, characterized in that: The method comprises 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 according to any one of claims 1 to 5 by calling the computer program stored in the memory.

9. A computer-readable storage medium for storing a computer program for performing high-cycle fatigue analysis of blades considering the influence of crystal orientation angles, characterized in that: When the computer program is run on a computer, the steps of the method according to any one of claims 1 to 5 are executed.

Images