A method and apparatus suitable for measuring vibrational stresses in turbine cooling vanes

Abstract

The present disclosure relates to the technical field of materials, and specifically provides a vibration stress measurement method, a vibration stress measurement device, an electronic device and a storage medium. The vibration stress measurement method comprises: obtaining a sensitivity parameter of a sensor; determining an elastic modulus parameter of a blade along a measurement direction according to a crystal orientation angle corresponding to the blade; obtaining a strain of an actual measurement position by the sensor when the sensor is located at the actual measurement position and located in the measurement direction; and calculating a stress of a theoretical measurement position by the strain, the elastic modulus parameter and the sensitivity parameter. The present disclosure provides a vibration stress measurement method, and solves the problem of inaccurate vibration stress measurement results.

Description

Technical Field

[0001] This disclosure relates to the field of materials, and more specifically, to a method for measuring the vibration stress of turbine cooling blades, an apparatus for measuring the vibration stress of turbine cooling blades, electronic equipment, and a storage medium. Background Technology

[0002] High-pressure turbine blades in aero-engines are subjected to harsh operating environments, including high temperatures, high pressures, and centrifugal forces generated by high-speed rotation. This makes them susceptible to high-cycle fatigue failure, potentially leading to major aviation accidents. The key to preventing high-cycle fatigue failure in high-pressure turbine cooling blades lies in whether the vibration stress during actual operation exceeds the fatigue strength limit. Therefore, accurately measuring the specific vibration stress values ​​of the blades during aero-engine development and testing is a pressing technical problem that needs to be solved.

[0003] In related technologies, vibration stress is measured by placing resistance strain gauges on the blade surface, measuring the vibration strain generated during operation, and then calculating the magnitude of the vibration stress using material constitutive relations. However, existing engine blades are made of nickel-based single-crystal high-temperature alloys, which exhibit significant crystal orientation effects. This affects the sensitivity of the strain gauges and also leads to large differences in the elastic modulus parameters in different directions at the same temperature. Errors in the elastic modulus parameters directly affect the measurement results of vibration stress. Since related technologies do not consider the crystal orientation effect of the material, the measurement results will have large errors and low accuracy.

[0004] It should be noted that the information disclosed in the background section above is only used to enhance the understanding of the background of this disclosure, and therefore may include information that does not constitute prior art known to those skilled in the art. Summary of the Invention

[0005] The purpose of this disclosure is to provide a method, an apparatus, an electronic device, and a storage medium suitable for measuring the vibration stress of turbine cooling blades, which can improve the accuracy of measurement results and reduce measurement errors.

[0006] According to one aspect of this disclosure, a method for measuring vibration stress of turbine cooling blades is provided, comprising: acquiring a sensitivity parameter of a sensor; wherein the sensitivity parameter is used to indicate the proportional relationship between the stress at a theoretical measurement position and an actual measurement position on the blade under the same external factors, the theoretical measurement position being the position on the blade most affected by vibration, and determining the actual measurement position based on the theoretical measurement position; determining the elastic modulus parameter of the blade along the measurement direction based on the crystal orientation angle corresponding to the blade; wherein the measurement direction is the direction of maximum normal stress at the actual measurement position; acquiring the strain at the actual measurement position through the sensor when the sensor is located at the actual measurement position and in the measurement direction; wherein the strain is used to indicate the relative change in the shape of the blade under external load, the external load being the product of stress and the area of ​​force application, and the stress being the internal force of interaction within the blade; and calculating the stress at the theoretical measurement position using the strain, the elastic modulus parameter, and the sensitivity parameter.

[0007] In an exemplary embodiment of this disclosure, the method further includes: determining a corresponding first vibration stress field based on resonance analysis results; wherein the resonance analysis results are used to indicate the resonance properties of the blade, the first vibration stress field is a relative vibration stress field, and the relative vibration stress field is used to indicate the relative magnitude of the stress at each position of the blade when the blade resonates; calculating a second vibration stress field based on the load-bearing properties of the blade; wherein the second vibration stress field is an allowable vibration stress field, and the allowable vibration stress field is used to indicate the maximum value of vibration stress that each position on the blade can withstand under the condition that high-cycle fatigue failure does not occur; and determining the theoretical measurement position on the blade based on the first vibration stress field and the second vibration stress field.

[0008] In one exemplary embodiment of this disclosure, the method for determining the actual measurement position based on the theoretical measurement position further includes: calculating the equivalent stress at the theoretical measurement position; and determining the actual measurement position based on the equivalent stress.

[0009] In an exemplary embodiment of this disclosure, the method for determining the actual measurement position based on equivalent stress further includes: determining a plurality of candidate measurement positions that satisfy preset conditions based on equivalent stress; wherein the preset conditions are that the ratio of the normal stress to the equivalent stress at the candidate measurement position exceeds a preset threshold; and determining the actual measurement position with the largest normal stress among the plurality of candidate measurement positions.

[0010] In an exemplary embodiment of this disclosure, the method for calculating the first vibration stress field based on resonance analysis results further includes: establishing a Cartesian coordinate system; wherein the first coordinate axis points to the blade axial direction, the second coordinate axis points to the blade circumferential direction, and the third coordinate axis points to the blade radial direction; obtaining the dendrite trunk of the crystal, and determining the first material coordinate axis, the second material coordinate axis, and the third material coordinate axis based on the dendrite trunk, so as to establish a material coordinate system based on the first material coordinate axis, the second material coordinate axis, and the third material coordinate axis; wherein the first material coordinate axis and the second material coordinate axis are the two principal directions of the secondary dendrite trunk, and the third material coordinate axis is the principal direction of the primary dendrite trunk, and the direction of the primary dendrite trunk is... The dendritic structure grows longitudinally along the preferred direction, which is the same as the

[001] crystal direction. The direction of the secondary dendrite trunk is the transverse preferred growth direction of the dendritic structure, which is the same as the

[010] and

[100] crystal directions. The

[001] crystal direction is perpendicular to the plane containing the

[010] and

[100] crystal directions. The angular deviation parameters between the Cartesian coordinate system and the material coordinate system are obtained. The angular deviation parameters are the angular deviation between the first coordinate axis and the first material coordinate axis, the angular deviation between the second coordinate axis and the second material coordinate axis, and the angular deviation between the third coordinate axis and the third material axis. The resonance analysis results are obtained by performing finite element analysis on the blade vibration characteristics and response using the angular deviation parameters.

[0011] In an exemplary embodiment of this disclosure, the method further includes: acquiring stress components at the actual measurement location; wherein the stress components include stress at the actual measurement location along a first coordinate axis, stress along a second coordinate axis, and stress along a third coordinate axis, the first, second, and third coordinate axes being coordinate axes of a Cartesian coordinate system established with the actual measurement location as the origin; determining multiple candidate direction vectors based on the actual measurement location and other locations on the blade besides the actual measurement location; calculating the normal stress on the blade surface along each candidate direction vector based on the stress components, and determining the candidate direction vector with the largest normal stress as the target direction vector; wherein the target direction vector is used to indicate the measurement direction.

[0012] In an exemplary embodiment of this disclosure, the method for determining the elastic modulus parameter of the blade along the measurement direction based on the crystal orientation angle corresponding to the blade further includes: establishing a Cartesian coordinate system on the blade and determining the Cartesian direction vector of the measurement direction in the Cartesian coordinate system; wherein the first coordinate axis points to the blade axial direction, the second coordinate axis points to the blade circumferential direction, the third axis points to the blade radial direction, and the Cartesian direction vector is the direction vector of the measurement direction in the Cartesian coordinate system; obtaining the dendrite trunk of the crystal and determining the first material coordinate axis, the second material coordinate axis, and the third material coordinate axis based on the dendrite trunk, so as to establish a material coordinate system based on the first material coordinate axis, the second material coordinate axis, and the third material coordinate axis; wherein the first material coordinate axis and the second material coordinate axis are the two principal directions of the secondary dendrite trunk, and the third material coordinate axis is the direction of the primary dendrite trunk, the direction of the primary dendrite trunk being the preferred longitudinal growth direction of the dendrite structure. Similar to the

[001] crystal orientation, the direction of the secondary dendrite trunk is the preferred transverse growth direction of the dendrite structure, similar to the

[010] and

[100] crystal orientations. The

[001] crystal orientation is perpendicular to the plane containing the

[010] and

[100] crystal orientations. The current temperature of the actual measurement location is obtained. The first elastic modulus parameter and the second elastic modulus parameter of the actual measurement location are obtained. The first elastic modulus parameter is the elastic modulus parameter along a preset direction when the actual measurement location is at the first temperature, and the second elastic modulus parameter is the elastic modulus parameter along a preset direction when the actual measurement location is at the second temperature. The first temperature, the second temperature, and the current temperature are different from each other. The material direction vector is calculated based on the Cartesian direction vector. The material direction vector is the direction vector of the measurement direction in the material coordinate system. The elastic modulus parameter is calculated based on the current temperature, the first elastic modulus parameter, the second elastic modulus parameter, and the material direction vector.

[0013] According to one aspect of this disclosure, an apparatus suitable for measuring the vibration stress of a turbine cooling blade is provided, comprising a first parameter acquisition module for acquiring a sensitivity parameter of a sensor; wherein the sensitivity parameter indicates the ratio of the stress experienced at a theoretical measurement position to that at an actual measurement position on the blade under the same external factors, the theoretical measurement position being the position on the blade most affected by vibration, and the actual measurement position being determined based on the theoretical measurement position; a second parameter acquisition module for determining the elastic modulus parameter of the blade along the measurement direction based on the crystal orientation angle corresponding to the blade; wherein the measurement direction is the direction of the maximum normal stress at the actual measurement position; a third parameter acquisition module for acquiring the strain at the actual measurement position via the sensor when the sensor is located at the actual measurement position and in the measurement direction; wherein the strain indicates the relative change in the shape of the blade under external load, the external load being the product of stress and the area of ​​force application, and the stress being the internal force of interaction within the blade; and a calculation module for calculating the stress at the theoretical measurement position using the strain, the elastic modulus parameter, and the sensitivity parameter.

[0014] According to one aspect of this disclosure, an electronic device is provided, comprising:

[0015] A processor; and a memory for storing executable instructions of the processor; wherein the processor is configured to perform a method suitable for measuring the vibration stress of a turbine cooling blade according to any of the above embodiments by executing the executable instructions.

[0016] According to one aspect of this disclosure, a computer-readable storage medium is provided that, when instructions in the storage medium are executed by an electronic device processor, enables the electronic device to perform a method for measuring the vibration stress of a turbine cooling blade according to any of the above embodiments.

[0017] The exemplary embodiments disclosed herein may have some or all of the following beneficial effects:

[0018] In a method for measuring the vibration stress of turbine cooling blades provided in the exemplary embodiments of this disclosure, the sensitivity parameters of a sensor are obtained, the elastic modulus parameters of the blade along the measurement direction are determined based on the crystal orientation angle corresponding to the blade, and the strain at the actual measurement position is obtained by the sensor when the sensor is located at the actual measurement position and in the measurement direction. The stress at the theoretical measurement position is calculated by using the strain, the elastic modulus parameters, and the sensitivity parameters. On the one hand, considering the influence of the crystal orientation angle on the elastic modulus parameters, the elastic modulus parameters are determined based on the crystal orientation angle and added to the final vibration stress calculation process, reducing the measurement result error caused by the elastic modulus parameter error due to crystal orientation, and improving the measurement accuracy. On the other hand, since the effect of crystal orientation is considered in the calculation of each stress field, the influence of crystal orientation on the sensitivity parameters is considered, thus ensuring the accuracy of the sensitivity. Furthermore, considering the three factors of temperature, steady-state stress, and relative vibration stress, the theoretical measurement position can be accurately identified, providing more practical reference value.

[0019] It should be understood that the above general description and the following detailed description are exemplary and explanatory only, and are not intended to limit this disclosure. Attached Figure Description

[0020] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with this disclosure and, together with the description, serve to explain the principles of this disclosure. It is obvious that the drawings described below are merely some embodiments of this disclosure, and those skilled in the art can obtain other drawings based on these drawings without any inventive effort.

[0021] Figure 1The illustration schematically shows a system architecture for a method of measuring the vibration stress of turbine cooling blades according to one embodiment of the present disclosure.

[0022] Figure 2 The flowchart schematically illustrates a method for measuring the vibration stress of a turbine cooling blade according to one embodiment of the present disclosure.

[0023] Figure 3 A flowchart illustrating the determination of a theoretical measurement location is shown schematically according to one embodiment of the present disclosure.

[0024] Figure 4 The flowchart illustrating the determination of the actual measurement position based on equivalent stress is shown in one embodiment of the present disclosure.

[0025] Figure 5 A flowchart illustrating the determination of the actual measurement location is shown in one embodiment of the present disclosure.

[0026] Figure 6 The flowchart illustrating resonance analysis based on crystal orientation angle is shown schematically in one embodiment of the present disclosure.

[0027] Figure 7 The diagram illustrates a Cartesian coordinate system and a material coordinate system according to one embodiment of the present disclosure.

[0028] Figure 8 The diagram illustrates a top view of a Cartesian coordinate system and a material coordinate system according to one embodiment of the present disclosure.

[0029] Figure 9 A flowchart illustrating the determination of measurement direction according to one embodiment of the present disclosure is shown schematically.

[0030] Figure 10 The flowchart illustrating the determination of the elastic modulus parameter based on the crystal orientation angle is shown schematically in one embodiment of the present disclosure.

[0031] Figure 11 The flowchart schematically illustrates a method for measuring the vibration stress of a turbine cooling blade according to one embodiment of the present disclosure.

[0032] Figure 12 A block diagram schematically illustrates an apparatus suitable for measuring the vibration stress of turbine cooling blades according to one embodiment of the present disclosure.

[0033] Figure 13 A schematic diagram of the structure of a computer system suitable for implementing the embodiments of the present disclosure is shown. Detailed Implementation

[0034] Example embodiments will now be described more fully with reference to the accompanying drawings. However, example embodiments can be implemented in many forms and should not be construed as preceding the examples set forth herein; rather, these embodiments are provided to make this disclosure more comprehensive and complete, and to fully convey the concept of the example embodiments to those skilled in the art. The described features, structures, or characteristics can be combined in any suitable manner in one or more embodiments. In the following description, numerous specific details are provided to give a full understanding of embodiments of this disclosure. However, those skilled in the art will recognize that the technical solutions of this disclosure can be practiced with one or more specific details omitted, or other methods, components, apparatus, steps, etc., can be employed. In other instances, well-known technical solutions are not shown or described in detail to avoid obscuring various aspects of this disclosure.

[0035] The flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of this disclosure. In this regard, each block in a flowchart or block diagram may represent a module, segment, or portion of code containing one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions indicated in the blocks may occur in a different order than those indicated in the drawings. For example, two consecutively indicated blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each block in a block diagram or flowchart, and combinations of blocks in a block diagram or flowchart, may be implemented using a dedicated hardware-based system that performs the specified function or operation, or using a combination of dedicated hardware and computer instructions.

[0036] The units described in the embodiments of this disclosure can be implemented in software or hardware, and the described units can also be located in a processor. The names of these units do not necessarily limit the unit itself.

[0037] Figure 1 A schematic diagram of a system architecture for an exemplary application environment in which a method and apparatus for measuring the vibration stress of turbine cooling blades, according to embodiments of the present disclosure, can be applied.

[0038] like Figure 1As shown, system architecture 100 may include one or more of terminal devices 101, 102, and 103, a network 104, and a server 105. Network 104 serves as the medium for providing communication links between terminal devices 101, 102, and 103 and server 105. Network 104 may include various connection types, such as wired, wireless communication links, or fiber optic cables. Terminal devices 101, 102, and 103 may be various electronic devices with displays, including but not limited to desktop computers, laptops, smartphones, and tablets. It should be understood that... Figure 1 The number of terminal devices, networks, and servers shown is merely illustrative. Depending on implementation needs, there can be any number of terminal devices, networks, and servers. For example, server 105 could be a server cluster composed of multiple servers.

[0039] The method for measuring the vibration stress of turbine cooling blades provided in this embodiment can be executed on server 105. Specifically, the method involves acquiring the sensitivity parameters of a sensor, determining the elastic modulus parameters of the blade along the measurement direction based on the crystal orientation angle corresponding to the blade, and acquiring the strain at the actual measurement location when the sensor is located at the actual measurement position and in the measurement direction. The stress at the theoretical measurement location is then calculated using the strain, elastic modulus parameters, and sensitivity parameters. The method for measuring the vibration stress of turbine cooling blades provided in this embodiment can also be executed by terminal devices 101, 102, and 103. Furthermore, the method can be jointly executed by terminal devices 101, 102, and 103 and server 105. This exemplary embodiment does not impose any special limitations on this approach.

[0040] High-pressure turbine blades in aero-engines are subjected to harsh operating environments, including high temperatures, high pressures, and centrifugal forces generated by high-speed rotation. This makes them susceptible to high-cycle fatigue failure, potentially leading to major aviation accidents. The key to preventing high-cycle fatigue failure in high-pressure turbine cooling blades lies in whether the vibration stress during actual operation exceeds the fatigue strength limit. Therefore, accurately measuring the specific vibration stress values ​​of the blades during aero-engine development and testing is a pressing technical problem that needs to be solved.

[0041] In related technologies, vibration stress is measured by placing resistance strain gauges on the blade surface, measuring the vibration strain generated during operation, and then calculating the magnitude of the vibration stress using materials mechanics. However, existing engine blades are made of orthotropic materials with significant crystal orientation effects, resulting in large differences in elastic modulus parameters in different directions at the same temperature. Errors in elastic modulus parameters directly affect the measurement results of vibration stress. Since related technologies do not consider the crystal orientation effect of the material, the measurement results will have large errors and low accuracy.

[0042] One exemplary embodiment of this disclosure provides a method suitable for measuring the vibration stress of turbine cooling blades, with reference to... Figure 2 As shown, a method for measuring the vibration stress of turbine cooling blades may include the following steps:

[0043] Step S210: Obtain the sensitivity parameters of the sensor; wherein, the sensitivity parameters are used to indicate the ratio of the stress on the theoretical measurement position and the actual measurement position on the blade under the same external factors. The theoretical measurement position is the position on the blade that is most affected by vibration. The actual measurement position is determined based on the theoretical measurement position.

[0044] Step S220: Determine the elastic modulus parameter of the blade along the measurement direction based on the crystal orientation angle corresponding to the blade; wherein, the measurement direction is the direction of maximum normal stress at the actual measurement position;

[0045] Step S230: When the sensor is located at the actual measurement position and in the measurement direction, the strain at the actual measurement position is obtained through the sensor; wherein, the strain is used to indicate the relative change in the shape of the blade under the action of external load, the external load is the product of stress and the area of ​​force application, and the stress is the internal force of the interaction inside the blade.

[0046] Step S240: Calculate the stress at the theoretical measurement location using strain, elastic modulus parameters, and sensitivity parameters.

[0047] The above steps will now be explained in more detail.

[0048] In one exemplary embodiment of this disclosure, a method suitable for measuring the vibration stress of turbine cooling blades is provided, with reference to... Figure 2 As shown, a method for measuring the vibration stress of turbine cooling blades includes the following steps S210 to S240:

[0049] In step S210, the sensitivity parameters of the sensor are obtained.

[0050] In one example embodiment of this disclosure, the sensitivity parameter of the sensor is obtained. This sensitivity parameter indicates the ratio of stress experienced at a theoretical measurement location to that experienced at an actual measurement location on the blade under the same external factors. The theoretical measurement location is often a stress concentration area on the turbine cooling blade. These locations are often unsuitable for direct measurement due to large stress gradients or curvatures. The actual measurement location is determined based on the theoretical measurement location. For example, the blade can be a nickel-based single-crystal high-temperature alloy blade, the sensor can be a resistance strain gauge, and the theoretical measurement location can be a location unsuitable for sensor placement, such as the film cooling hole or the blade root chamfer. The actual measurement location can be any suitable location for sensor placement other than the theoretical measurement location. The resistance strain gauge measures the magnitude of the vibration strain generated by the blade during operation, and then calculates the stress based on the strain value.

[0051] For example, if the blade is a nickel-based single-crystal high-temperature alloy blade, the theoretical measurement position can be the location of a certain air film hole on the blade. The Von Mises stress at the theoretical measurement position is calculated as the equivalent stress. After determining the actual measurement position based on the theoretical measurement, the normal stress at the actual measurement position is calculated. The normal stress is the stress perpendicular to the stress section. The ratio of the equivalent stress to the normal stress is used as the sensitivity parameter of the sensor.

[0052] It is understood that the equivalent stress at the theoretical measurement location is not limited to Von Mises stress. Other methods can also be used to calculate the equivalent stress at the theoretical measurement location. For example, the principal stress at the theoretical measurement location can also be calculated as the equivalent stress. This example embodiment does not limit this.

[0053] In step S220, the elastic modulus parameter of the blade along the measurement direction is determined based on the crystal orientation angle corresponding to the blade.

[0054] In one exemplary embodiment of this disclosure, the elastic modulus parameter of the blade along the measurement direction is determined based on the crystal orientation angle corresponding to the blade. The measurement direction is the direction of maximum normal stress at the actual measurement location.

[0055] Specifically, since the normal stress at the actual measurement location is different in each direction, the direction with the largest normal stress is selected as the measurement direction. After obtaining the measurement direction, the elastic modulus parameter along the measurement direction at the actual measurement location is determined.

[0056] For example, if the blade is a nickel-based single-crystal superalloy blade, at the actual measurement position on the blade surface, the elastic modulus parameter along the measurement direction with the maximum normal stress is determined according to the crystal orientation angle corresponding to the blade. The crystal orientation angle of the blade is used to indicate the covalent crystal orientation of the nickel-based single crystal.

[0057] The directionality of covalent crystals refers to the formation of covalent bonds in a crystal along a specific direction. According to the quantum theory of covalent bonds, the strength of a covalent bond depends on the degree of overlap of the electron clouds. Due to the asymmetry of the distribution of electrons in non-full shells, bonding always occurs in the direction of maximum electron cloud density.

[0058] In step S230, when the sensor is located at the actual measurement position and in the measurement direction, the strain at the actual measurement position is obtained through the sensor.

[0059] In one exemplary embodiment of this disclosure, when the sensor is located at the actual measurement position and in the measurement direction, the strain at the actual measurement position is acquired by the sensor. The strain is used to indicate the relative change in the shape of the blade under external load, and the stress is the internal force of the interaction within the blade.

[0060] Specifically, a sensor is installed at the actual measurement location of the blade, and the measurement direction of the sensor is consistent with the direction of the maximum normal stress at the actual measurement location. The strain of the blade at the actual measurement location is obtained through the sensor.

[0061] For example, the blade is a nickel-based single-crystal high-temperature alloy blade, and the sensor is a resistance strain gauge. The resistance strain gauge is set at the actual measurement position of the blade, and the strain of the blade along the measurement direction at the actual measurement position is obtained through the resistance strain gauge.

[0062] In step S240, the stress at the theoretical measurement location is calculated using strain, elastic modulus parameters, and sensitivity parameters.

[0063] In one example embodiment of this disclosure, the stress at the theoretical measurement location is calculated using strain, elastic modulus parameters, and sensitivity parameters. The stress at the actual measurement location can be calculated using strain and elastic modulus parameters, and the ratio of the stress at the actual measurement location to the stress at the theoretical measurement location satisfies the sensitivity parameters.

[0064] Specifically, since the ratio of the stress at the actual measurement location to the stress at the theoretical measurement location satisfies the sensitivity parameter, the stress at the actual measurement location can be calculated from the strain measured by the sensor and the elastic modulus parameter. Then, the stress at the theoretical measurement location can be calculated based on the stress at the actual measurement location and the sensitivity parameter.

[0065] In one exemplary embodiment of this disclosure, a corresponding first vibration stress field is confirmed based on the resonance analysis results, a second vibration stress field is calculated based on the load properties of the blade, and a theoretical measurement position is determined on the blade based on the first and second vibration stress fields. Specifically, refer to... Figure 3 As shown, determining the theoretical measurement location may include the following steps S310 to S330:

[0066] In step S310, the first vibration stress field is calculated based on the resonance analysis results.

[0067] In one example embodiment of this disclosure, a corresponding first vibration stress field is determined based on the resonance analysis results. The resonance analysis results are used to indicate the resonance properties of the blade, and the first vibration stress field is a relative vibration stress field, which indicates the relative magnitude of the stress experienced at different positions on the blade when it resonates.

[0068] In step S320, the second vibration stress field is calculated based on the load properties of the blade.

[0069] In one example embodiment of this disclosure, a second vibration stress field is calculated based on the load-bearing properties of the blade. The load-bearing properties indicate the maximum allowable stress the blade can withstand, and the second vibration stress field is an allowable vibration stress field, which indicates the maximum distribution of vibration stress that can be withstood at various locations on the blade without high-cycle fatigue failure.

[0070] In step S330, the theoretical measurement position is determined on the blade based on the first vibration stress field and the second vibration stress field.

[0071] In one exemplary embodiment of this disclosure, a theoretical measurement location is determined on the blade based on a first vibration stress field and a second vibration stress field. This theoretical measurement location is the position with the minimum vibration intensity reserve, which is most susceptible to resonance. For example, the theoretical measurement location could be a blade film vent, a blade root chamfer, or a blade tail slit.

[0072] For example, the first vibration stress field is a dimensionless relative vibration stress field, and the second vibration stress field is an allowable vibration stress field. Based on the allowable vibration stress field and the relative vibration stress field, the position with the minimum vibration intensity reserve can be obtained. If these positions occur at the air film hole position, the air film hole position is taken as the theoretical measurement position.

[0073] In one exemplary embodiment of this disclosure, a first vibration stress field is obtained by establishing a finite element modal analysis model, and a second vibration stress field is obtained by establishing a finite element static strength analysis model. Specifically, a finite element modal analysis model is established and a Campbell diagram is plotted to perform resonance analysis on the blade, determining the vibration stress field, frequency, and mode shape, etc., used to indicate the vibration stress distribution of the blade. The vibration stress field is divided by the maximum vibration stress of the blade to obtain a dimensionless relative vibration stress field, which is used to indicate the relative magnitude of the stress at each position of the blade during resonance. A steady-state stress field for the load-bearing properties of the blade is obtained by establishing a finite element static strength analysis model. The allowable vibration stress field of the blade is obtained by plotting the steady-state stress field using a Goodman curve, which indicates the maximum value of vibration stress that each position on the blade can withstand.

[0074] Through steps S310 to S330 above, the corresponding first vibration stress field is confirmed based on the resonance analysis results, the second vibration stress field is calculated based on the load properties of the blade, and the theoretical measurement position on the blade is determined based on the first and second vibration stress fields. The analysis of the first and second vibration stress fields can accurately identify the locations on the blade where high-cycle fatigue failure is more likely to occur, improving applicability.

[0075] In one example embodiment of this disclosure, the equivalent stress at the theoretical measurement location is calculated, and the actual measurement location is determined based on the equivalent stress. Specifically, refer to... Figure 4 As shown, the flowchart for determining the actual measurement location based on equivalent stress may include the following steps S410 to S420:

[0076] In step S410, the equivalent stress at the theoretical measurement location is calculated.

[0077] In one exemplary embodiment of this disclosure, the equivalent stress at the theoretically measured location is calculated. The equivalent stress is the equivalent stress experienced by the blade at the point of strength failure, calculated from the material constitutive relation. For example, the equivalent stress can be either Von Mises stress or principal stress.

[0078] For example, the theoretical measurement location is the location of the film cooling pores on the blade, and the Von Mises stress at that film cooling pore location is calculated as the equivalent stress.

[0079] It is understood that the equivalent stress at the theoretical measurement location is not limited to Von Mises stress. Other methods can also be used to calculate the equivalent stress at the theoretical measurement location. For example, the principal stress at the theoretical measurement location can also be calculated as the equivalent stress. This example embodiment does not limit this.

[0080] In step S420, the actual measurement position is determined based on the equivalent stress.

[0081] In one exemplary embodiment of this disclosure, the actual measurement location is determined based on equivalent stress. The theoretical measurement location may be unsuitable for direct measurement due to a large stress gradient or curvature; therefore, another actual measurement location needs to be determined based on the theoretical measurement location. For example, the actual measurement location could be a position on the blade where a sensor can be directly installed to measure the stress.

[0082] Specifically, the actual measurement location is determined based on the equivalent stress at the theoretical measurement location, taking into account the engineering design requirements.

[0083] For example, after calculating the equivalent stress at the theoretical measurement location, the actual measurement location is selected based on engineering design requirements, where the ratio of VonMises stress to equivalent stress is greater than 20% and there are no air film holes or chamfers.

[0084] Through steps S410 to S420 above, the equivalent stress at the theoretical measurement location is calculated, and the actual measurement location is determined based on the equivalent stress. Stress measurement at the actual location is more accurate than at the theoretical location, and it also increases the sensor's survivability during blade testing, reducing the difficulty of stress measurement.

[0085] In one example embodiment of this disclosure, multiple candidate measurement locations that meet preset conditions are determined based on equivalent stress, and the actual measurement location with the maximum normal stress is determined from among the multiple candidate measurement locations. Specifically, refer to... Figure 5 As shown, determining the actual measurement location may include the following steps S510 to S520:

[0086] In step S510, multiple candidate measurement positions that meet preset conditions are determined based on equivalent stress.

[0087] In one exemplary embodiment of this disclosure, multiple candidate measurement locations satisfying preset conditions are determined based on equivalent stress. The preset conditions include a ratio of normal stress to equivalent stress exceeding a preset threshold at the candidate measurement location, and the candidate measurement locations being locations on the blade other than the theoretical measurement locations that satisfy the preset conditions. For example, the preset conditions might include a ratio of equivalent stress at the candidate measurement location to equivalent stress at the theoretical measurement location exceeding 20%, the absence of film pores, and blade root chamfering, among other conditions.

[0088] For example, after determining the candidate measurement locations and their equivalent stresses, all locations on the blade that satisfy the requirement that the ratio of the equivalent stress to the theoretical measurement location is greater than 20%, that there are no film pores, and that there is no blade root chamfer are selected as candidate measurement locations.

[0089] In step S520, the actual measurement location with the maximum normal stress is determined from multiple candidate measurement locations.

[0090] In one exemplary embodiment of this disclosure, the actual measurement location with the highest normal stress is determined from multiple candidate measurement locations. Here, normal stress is tensile or compressive stress in a certain direction, and the location with the highest normal stress among the multiple candidate measurement locations is the most suitable actual measurement location for setting up the sensor.

[0091] Specifically, the normal stress at each candidate measurement location is calculated, and the candidate measurement location with the highest normal stress is taken as the actual measurement location.

[0092] Through steps S510-S520 above, multiple candidate measurement locations that meet preset conditions are determined based on equivalent stress. The actual measurement location with the highest normal stress is then selected from these candidate locations. This selection of the most suitable actual measurement location for sensor placement further reduces the difficulty of stress measurement.

[0093] In one exemplary embodiment of this disclosure, a Cartesian coordinate system is established, the dendrite trunk of the crystal is obtained, and a first material coordinate axis, a second material coordinate axis, and a third material coordinate axis are determined based on the dendrite trunk. This allows a material coordinate system to be established based on the first, second, and third material coordinate axes. An angular deviation parameter between the Cartesian coordinate system and the material coordinate system is obtained, and the vibration characteristics of the blade are analyzed using this angular deviation parameter to obtain resonance analysis results. Specifically, refer to... Figure 6 As shown, vibrational characteristic analysis based on crystal orientation angle may include the following steps S610 to S640:

[0094] In step S610, a Cartesian coordinate system is established.

[0095] In one exemplary embodiment of this disclosure, a Cartesian coordinate system is established. The first coordinate axis points along the blade axial direction, the second axis points along the blade circumferential direction, and the third axis points along the blade radial direction.

[0096] Specifically, a Cartesian coordinate system is established on the blade in a virtual environment. The first coordinate axis of the Cartesian coordinate system points to the blade axis, the second coordinate axis points to the blade circumferential direction, and the third coordinate axis points to the blade height direction.

[0097] In step S620, the dendrite trunk of the crystal is obtained, and the first material coordinate axis, the second material coordinate axis and the third material coordinate axis are determined according to the dendrite trunk, so that a material coordinate system is established according to the first material coordinate axis, the second material coordinate axis and the third material coordinate axis.

[0098] In one example embodiment of this disclosure, the dendritic trunk of a crystal is obtained, and a first material coordinate axis, a second material coordinate axis, and a third material coordinate axis are determined based on the dendritic trunk, so that a material coordinate system is established based on the first material coordinate axis, the second material coordinate axis, and the third material coordinate axis. Wherein, the first material coordinate axis and the second material coordinate axis are the two principal directions of the secondary dendritic trunk, respectively, and the third material coordinate axis is the direction of the primary dendritic trunk. The direction of the primary dendritic trunk is the longitudinal preferred growth direction of the dendritic structure, which is the same as the

[001] crystal direction. The direction of the secondary dendritic trunk is the transverse preferred growth direction of the dendritic structure, which is the same as the

[010] and

[100] crystal directions. The

[001] crystal direction is perpendicular to the plane containing the

[010] and

[100] crystal directions. For example, the blade is made of nickel-based single crystal high-temperature alloy. The primary dendrite trunk direction is the preferred growth direction of the primary dendrite trunk of the nickel-based single crystal, and the secondary dendrite trunk direction is the two main preferred growth directions of the secondary dendrite trunk of the nickel-based single crystal. The first material coordinate axis points to the

[100] crystal direction, the second material coordinate axis points to the

[010] crystal direction, and the third material coordinate axis points to the

[001] crystal direction.

[0099] [l1l2l3] is a way of indicating crystal orientation in crystal structure. Crystal orientation is the orientation of a crystal row, which is a straight line determined by any two lattice points in the crystal. l1, l2, and l3 are coprime constants, which are generally called crystal row indices.

[0100] Specifically, based on the crystal orientation characteristics of the blade material, the crystal axis is determined, and a material coordinate system is established based on the crystal axis.

[0101] For example, refer to Figure 7 As shown, the blade is a nickel-based single-crystal high-temperature alloy blade. In the virtual environment, a Cartesian coordinate system is established on the blade. The X-axis (first coordinate axis) points to the blade axis, the Y-axis (second coordinate axis) points to the blade circumferential direction, and the Z-axis (third coordinate axis) points to the blade height direction. Then, based on the crystal orientation characteristics of the nickel-based single crystal, the crystal axis is determined and a material coordinate system is established based on the crystal axis. The

[100] axis (first material coordinate axis) is the coordinate axis pointing to the

[100] crystal direction, the

[010] axis (second material coordinate axis) is the coordinate axis pointing to the

[010] crystal direction, and the

[001] axis (third material coordinate axis) is the coordinate axis pointing to the

[001] crystal direction.

[0102] In step S630, the angular deviation parameters between the Cartesian coordinate system and the material coordinate system are obtained.

[0103] In one exemplary embodiment of this disclosure, angular deviation parameters between the Cartesian coordinate system and the material coordinate system are obtained. These angular deviation parameters are the angular deviations between the first coordinate axis and the first material coordinate axis, the second coordinate axis and the second material coordinate axis, and the third coordinate axis and the third material coordinate axis. For example, the angular deviation parameters can be obtained by rotating the material coordinate system in ZYZ Euler angle sequence, using the Cartesian coordinate system as a reference system.

[0104] For example, using the Cartesian coordinate system as the reference frame, the material coordinate system is rotated in the ZYZ Euler angle sequence, with the reference frame as the reference. Figure 7 As shown, the first rotation is a rotation around the Z-axis that makes the first coordinate axis coincide with the first material coordinate axis, thus obtaining the angular deviation between the first coordinate axis and the first material coordinate axis. The second rotation is a rotation around the Y-axis after the first rotation that makes the second coordinate axis coincide with the second material coordinate axis, thus obtaining the angular deviation between the second coordinate axis and the second material coordinate axis. The third rotation is a rotation around the Z-axis after the second rotation that makes the third coordinate axis coincide with the third material coordinate axis, thus obtaining the angular deviation between the third coordinate axis and the third material coordinate axis.

[0105] A top-down view of the Cartesian coordinate system and the material coordinate system in a virtual environment, as shown below. Figure 8 As shown, the X-axis is the first coordinate axis, the Y-axis is the second coordinate axis, the

[100] axis is the first material coordinate axis, and the

[010] axis is the second material coordinate axis.

[0106] In step S640, the resonance analysis results are obtained by performing finite element analysis on the blade vibration characteristics using the angle deviation parameter.

[0107] In one exemplary embodiment of this disclosure, vibration finite element analysis of the blade is performed using angular deviation parameters to obtain resonance analysis results. The angular deviation parameters are input parameters for the finite element analysis.

[0108] For example, a finite element modal analysis model is constructed and a Campbell diagram is plotted. The angular deviation parameter is input into the finite element vibration analysis model to analyze the vibration characteristics and response of the blade, and the resonance analysis results used to indicate the resonance properties of the blade are obtained.

[0109] Through the above steps S610 to S640, a Cartesian coordinate system is established, the dendrite trunk of the crystal is obtained, and the first material coordinate axis, the second material coordinate axis, and the third material coordinate axis are determined based on the dendrite trunk, so that a material coordinate system is established based on the first material coordinate axis, the second material coordinate axis, and the third material coordinate axis. The angular deviation parameter between the Cartesian coordinate system and the material coordinate system is obtained, and the resonance analysis result is obtained by performing finite element analysis on the vibration characteristics and response of the blade through the angular deviation parameter.

[0110] In one example embodiment of this disclosure, the stress components at the actual measurement location are obtained. Multiple candidate direction vectors are determined based on the actual measurement location and other locations on the blade besides the actual measurement location. The normal stress on the blade surface along each candidate direction vector is calculated based on the stress components. The candidate direction vector with the largest normal stress is determined as the target direction vector. Specifically, refer to... Figure 9 As shown, determining the measurement direction may include the following steps S910 to S930:

[0111] In step S910, the stress components at the actual measurement location are obtained.

[0112] In one example embodiment of this disclosure, stress components at the actual measurement location are obtained. These stress components include stress along a first coordinate axis, stress along a second coordinate axis, and stress along a third coordinate axis at the actual measurement location. The first, second, and third coordinate axes are coordinate axes of a Cartesian coordinate system established with the actual measurement location as the origin.

[0113] Specifically, in the Cartesian coordinate system, the stress components along the first coordinate axis, the second coordinate axis, and the third coordinate axis are obtained at the actual measurement position.

[0114] In step S920, multiple candidate direction vectors are determined based on the actual measurement position and other positions on the blade besides the actual measurement position.

[0115] In one exemplary embodiment of this disclosure, multiple candidate direction vectors are determined based on the actual measurement position and other positions on the blade besides the actual measurement position. These candidate direction vectors include all direction vectors on the blade that pass through the actual measurement position.

[0116] Specifically, all direction vectors emanating outward from the actual measurement position on the blade are obtained as candidate direction vectors.

[0117] For example, the coordinates of the actual measurement position and other positions on the blade besides the actual measurement position in the Cartesian coordinate system are read, and the coordinates of the actual measurement point are subtracted from the coordinates of the other positions besides the actual measurement position to obtain multiple candidate direction vectors.

[0118] In step S930, the normal stress on the blade surface along each candidate direction vector is calculated based on the stress components, and the candidate direction vector with the largest normal stress is determined as the target direction vector.

[0119] In one example embodiment of this disclosure, the normal stress on the blade surface along each candidate direction vector is calculated based on the stress components, and the candidate direction vector with the largest normal stress is determined as the target direction vector. The target direction vector is used to indicate the measurement direction.

[0120] Specifically, based on the constitutive theory of materials, the normal stress in each direction is calculated based on the candidate direction vector and stress components, and the candidate direction vector with the largest normal stress is selected as the target direction vector.

[0121] For example, according to the static equilibrium equation, each stress component is multiplied by the cosine of the candidate direction vector in the direction of each stress component to obtain each candidate stress component. The projection of each candidate stress component onto the candidate direction vector is calculated and summed to obtain the normal stress in the candidate direction. The candidate direction vector with the largest normal stress is selected as the target direction vector.

[0122] Through steps S910 to S930, the stress components at the actual measurement location are obtained. Multiple candidate direction vectors are determined based on the actual measurement location and other locations on the blade besides the actual measurement location. The normal stress along each candidate direction vector on the blade surface is calculated based on the stress components. The candidate direction vector with the highest normal stress is determined as the target direction vector. This provides the optimal measurement direction at the actual measurement location, further improving measurement accuracy.

[0123] In one exemplary embodiment of this disclosure, a Cartesian coordinate system is established on the blade, and the Cartesian direction vector of the measurement direction in the Cartesian coordinate system is determined. The dendrite trunk of the crystal is obtained, and a first material coordinate axis, a second material coordinate axis, and a third material coordinate axis are determined based on the dendrite trunk. This allows a material coordinate system to be established based on the first, second, and third material coordinate axes. The current temperature at the actual measurement location is obtained, as are the first and second elastic modulus parameters at the actual measurement location. The material direction vector is calculated based on the Cartesian direction vector, and the elastic modulus parameter is calculated based on the current temperature, the first elastic modulus parameter, the second elastic modulus parameter, and the material direction vector. Specifically, refer to... Figure 10 As shown, determining the elastic modulus parameter based on the crystal orientation angle may include the following steps S1010 to S1060:

[0124] In step S1010, a Cartesian coordinate system is established on the blade, and the Cartesian direction vector of the measurement direction in the Cartesian coordinate system is determined.

[0125] In one exemplary embodiment of this disclosure, a Cartesian coordinate system is established on the blade, and the Cartesian direction vector of the measurement direction in the Cartesian coordinate system is determined. The first coordinate axis points to the blade axial direction, the second coordinate axis points to the blade circumferential direction, the third axis points to the blade radial direction, and the Cartesian direction vector is the direction vector of the measurement direction in the Cartesian coordinate system.

[0126] Specifically, a Cartesian coordinate system is established on the blade in a virtual environment. The first coordinate axis of the Cartesian coordinate system points to the blade axial direction, the second coordinate axis points to the blade circumferential direction, and the third coordinate axis points to the blade radial direction. The direction vector of the measurement direction in the Cartesian coordinate system is also determined.

[0127] In step S1020, the dendrite trunk of the crystal is obtained, and the first material coordinate axis, the second material coordinate axis and the third material coordinate axis are determined according to the dendrite trunk, so that a material coordinate system is established according to the first material coordinate axis, the second material coordinate axis and the third material coordinate axis, and the direction vector of the measurement direction in the Cartesian coordinate system is determined.

[0128] In one example embodiment of this disclosure, the dendritic trunk of a crystal is obtained, and a first material coordinate axis, a second material coordinate axis, and a third material coordinate axis are determined based on the dendritic trunk, so that a material coordinate system is established based on the first material coordinate axis, the second material coordinate axis, and the third material coordinate axis. Wherein, the first material coordinate axis and the second material coordinate axis are the two principal directions of the secondary dendritic trunk, respectively, and the third material coordinate axis is the direction of the primary dendritic trunk. The direction of the primary dendritic trunk is the longitudinal preferred growth direction of the dendritic structure, which is the same as the

[001] crystal direction. The direction of the secondary dendritic trunk is the transverse preferred growth direction of the dendritic structure, which is the same as the

[010] and

[100] crystal directions. The

[001] crystal direction is perpendicular to the plane containing the

[010] and

[100] crystal directions.

[0129] Specifically, based on the crystal orientation characteristics of the blade material, the crystal axis is determined, and a material coordinate system is established based on the crystal axis.

[0130] For example, refer to Figure 7 As shown, the blade is a nickel-based single-crystal high-temperature alloy blade. In the virtual environment, a Cartesian coordinate system is established on the blade. The X-axis (first coordinate axis) points to the blade axis, the Y-axis (second coordinate axis) points to the blade circumferential direction, and the Z-axis (third coordinate axis) points to the blade height direction. Based on the crystal orientation characteristics of the nickel-based single crystal, the crystal axis is determined and a material coordinate system is established based on the crystal axis. The

[100] axis (first material coordinate axis) is the coordinate axis pointing to the

[100] crystal direction, the

[010] axis (second material coordinate axis) is the coordinate axis pointing to the

[010] crystal direction, and the

[001] axis (third material coordinate axis) is the coordinate axis pointing to the

[001] crystal direction. The origin of the Cartesian coordinate system and the origin of the material coordinate system are the same origin.

[0131] In step S1030, the current temperature at the actual measurement location is obtained.

[0132] In step S1040, the first elastic modulus parameter and the second elastic modulus parameter at the actual measurement location are obtained.

[0133] In one example embodiment of this disclosure, a first elastic modulus parameter and a second elastic modulus parameter are obtained at the actual measurement location. The first elastic modulus parameter is the elastic modulus parameter along a preset direction when the actual measurement location is at a first temperature, and the second elastic modulus parameter is the elastic modulus parameter along a preset direction when the actual measurement location is at a second temperature. The first temperature, the second temperature, and the current temperature are different from each other. For example, the preset direction may be the

[111] crystal orientation.

[0134] Specifically, the first elastic modulus parameter of the blade along a preset direction is obtained when the blade is at a first temperature, and the second elastic modulus parameter of the blade along a preset direction is obtained when the blade is at a second temperature.

[0135] For example, with the preset direction being the

[111] crystal orientation, the first elastic modulus parameter along the

[111] crystal orientation is obtained when the blade is at a first temperature, and the second elastic modulus parameter along the

[111] crystal orientation is obtained when the blade is at a second temperature.

[0136] In step S1050, the material orientation vector is calculated based on the Cartesian orientation vector.

[0137] In one exemplary embodiment of this disclosure, the material direction vector is calculated based on the Cartesian direction vector. The material direction vector is the direction vector of the measurement direction in the material coordinate system.

[0138] Specifically, the material direction vector in the material coordinate system is calculated based on the Cartesian direction vector in the Cartesian coordinate system and the angular deviation parameter between the Cartesian coordinate system and the material coordinate system.

[0139] For example, using the Cartesian coordinate system as a reference, the material coordinate system is rotated sequentially by the ZYZ Euler angles to obtain the angle deviation parameters. Then, based on the Cartesian direction vector in the Cartesian coordinate system and the angle deviation parameters, the material direction vector of the measurement direction in the material coordinate system is calculated.

[0140] In step S1060, the elastic modulus parameter is calculated based on the current temperature, the first elastic modulus parameter, the second elastic modulus parameter, and the material direction vector.

[0141] In one example embodiment of this disclosure, the elastic modulus parameter is calculated based on the current temperature, a first elastic modulus parameter, a second elastic modulus parameter, and a material direction vector. The elastic modulus parameter is the elastic modulus parameter along the measurement direction at the actual measurement location.

[0142] For example, firstly, the elastic modulus parameter of the actual measurement position along the preset direction is calculated based on the current temperature, the first elastic modulus parameter, and the second elastic modulus parameter. Then, the direction vector of the Cartesian direction vector used to indicate the measurement direction in the Cartesian coordinate system is calculated in the material coordinate system. Finally, the elastic modulus parameter of the actual measurement position along the measurement direction is calculated based on the elastic modulus parameter of the actual measurement position along the preset direction and the direction vector in the material coordinate system.

[0143] Through steps S1010 to S1060, a Cartesian coordinate system is established on the blade, and the Cartesian direction vector of the measurement direction in the Cartesian coordinate system is determined. The dendrite trunk of the crystal is obtained, and the first, second, and third material coordinate axes are determined based on the dendrite trunk. This allows the establishment of a material coordinate system based on the first, second, and third material coordinate axes. The current temperature at the actual measurement location is obtained, as are the first and second elastic modulus parameters at the actual measurement location. The material direction vector is calculated based on the Cartesian direction vector, and the elastic modulus parameter is calculated based on the current temperature, the first elastic modulus parameter, the second elastic modulus parameter, and the material direction vector. By comprehensively considering the temperature and the magnitude of the elastic modulus parameter along the patch direction when there is a crystal orientation deviation, the dynamic stress measurement becomes more accurate.

[0144] Figure 11 A flowchart of a method for measuring the vibration stress of turbine cooling blades, which is proposed as an exemplary embodiment of the present disclosure, includes the following steps S1101 to S1114.

[0145] Step S1101: Measure the crystal orientation angle of the blade. Specifically, measure the crystal orientation angle of the blade to which the strain gauge (sensor) needs to be installed, and establish the physical coordinate system and material coordinate system after the blade is assembled.

[0146] Step S1102: Establish a finite element modal analysis model. Specifically, obtain the angular deviations α (angle deviation between the first coordinate axis and the first material coordinate axis), β (angle deviation between the second coordinate axis and the second material coordinate axis), and γ (angle deviation between the third coordinate axis and the third material coordinate axis) between the physical coordinate system and the material coordinate system. Based on α, β, and γ, establish a finite element modal analysis model of the blade to obtain the vibration stress field and natural frequency under different working conditions and at different orders.

[0147] Step S1103: Draw the Campbell diagram. Specifically, based on the natural frequencies of the turbine blades under different operating conditions and orders, a Campbell diagram is drawn for resonance analysis. The analysis yields the corresponding operating conditions and orders that do not meet the resonance margin design requirements (resonance analysis results).

[0148] Step S1104: Calculate the dimensionless relative vibration stress field. Specifically, divide the relative vibration stress field of the corresponding working condition that does not meet the resonance margin design requirements and the corresponding order by the maximum relative vibration stress of that order to obtain the dimensionless relative vibration stress field (first vibration stress field) of that order.

[0149] Step S1105: Establish a finite element static strength model. Specifically, establish a finite element static strength analysis model of the turbine blade under this working condition, and calculate the vibration load properties of the blade when resonance occurs through the finite element static strength analysis model.

[0150] Step S1106: Calculate the allowable vibration stress field. Specifically, based on the finite element static strength analysis results, the allowable vibration stress field for this blade is obtained from the Goodman curve.

[0151] Step S1107: Determine the danger point. Specifically, divide the allowable stress field (second vibration stress field) by the dimensionless relative vibration stress field (first vibration stress field) to obtain the location with the minimum vibration intensity reserve, which is the danger point of vibration (actual measurement location).

[0152] Step S1108 involves blade shape analysis. Specifically, since the blade has numerous film pores, and the critical points (actual measurement locations) are mostly distributed in areas under multiaxial stress, such as the film pores and the blade root chamfer, Von Mises stress (equivalent stress) is selected as the assessment stress at these locations and denoted as σ. max,model .

[0153] Step S1109: Determine the measuring point. Specifically, the sub-steps for determining the measuring point (actual measurement location) are as follows:

[0154] 1) Determine the patchable area (multiple candidate measurement locations). The requirements for the patchable area are: the Von Mises equivalent stress in the area is relatively large, the ratio of the Von Mises equivalent stress to the critical point (target measurement location) should be greater than 20%, the stress gradient is small, the equivalent stress variation range in the selected area is about 5%, and there are no air film holes or chamfers in the area.

[0155] 2) Based on the results of the finite element modal analysis model, the blade vibration modes are divided into the following three types: pure bending (including one bending and two bending), pure torsion (including one torsion and two torsion), and bending-torsion combination.

[0156] 3) The intersection of the possible strain gauge placement areas (candidate measurement locations) under pure bending and combined bending-torsional vibration modes on the blade, or the intersection of the possible strain gauge placement areas (candidate measurement locations) under pure torsion and combined bending-torsional vibration modes, is taken as the possible strain gauge placement area. This allows a single strain gauge (sensor) to measure the vibration stress magnitude under various operating conditions. Finally, the node with the larger Von Mises (equivalent stress) is selected as the measurement point (actual measurement location) within the possible strain gauge placement area.

[0157] Step S1110: Determine the patch orientation. Specifically, the sub-steps for determining the patch orientation (measurement orientation) are as follows:

[0158] 1) Establish a Cartesian coordinate system with the measuring point (actual measurement position) as the origin, where the leaf height direction is taken as the Z-axis, the axial direction is taken as the X-axis, and the circumferential direction is taken as the Y-axis;

[0159] 2) Read the coordinates of the measuring point (actual measurement position) and the nodes around the measuring point on the blade surface (other positions besides the actual measurement position) in the Cartesian coordinate system, and subtract the coordinates of the measuring point from the coordinates of the surrounding points to establish the direction vectors [hkl] (candidate direction vectors) in different directions;

[0160] Calculate the cosines of the angles between the [hkl] direction (candidate direction vector) and the three coordinate axes: cosx, cosy, cosz, and let the three cosine values ​​be a, b, and c, respectively.

[0161] Substituting the stress components at the measuring point (actual measurement location) into the formula for calculating the normal stress of an inclined section in elasticity, the normal stress along different [hkl] directions (candidate direction vectors) is calculated, where σ n For normal stress, σ x Let σ be the stress component along the x-axis. y Let σ be the stress component along the y-axis. z Let σ be the stress component along the z-axis. xy Let σ be the stress component along the bisectors of the x and y axes. yz The stress components along the bisectors of the y-axis and z-axis are σ. xz These are the stress components along the bisecting lines of the x-axis and z-axis.

[0162] σ n =σ x a 2 +σ y b 2 +σ z c 2 +2σ xy ab+2σ yz bc+2σ xz ac

[0163] Choose the direction of maximum normal stress as the patch direction (measurement direction).

[0164] Step S1111: Calculate the sensor's sensitivity parameters. The specific formula is as follows, where σ n The normal stress at the measuring point (actual measurement point), σ max,model θ represents the equivalent stress at the danger point (theoretical measurement location), and θ is the sensitivity parameter.

[0165]

[0166] Step S1112: Determine whether the measuring point (actual measurement location) meets the preset conditions. If not, return to step S1109. If it does, perform strain gauge patching (sensor setup) to measure the vibration strain ε and continue to step S1113. The preset condition is whether the strain gauge sensitivity at the measuring point meets the condition θ≥20%.

[0167] Step S1113, calculate the elastic modulus parameter in the patch direction (measurement direction), the sub-steps of which are as follows:

[0168] 1) Read the temperature value T (current temperature) of the measuring point (actual measurement location). When the temperature value of the measuring point is between T1 and T2 (first temperature and second temperature), perform piecewise linear temperature interpolation on the elastic modulus parameters in the

[111] and

[001] directions respectively to obtain the elastic modulus parameters E of the measuring point in the two directions at that temperature.

[111] and E

[001] The interpolation formula is as follows, where E 1

[111] Let E be the elastic modulus parameter along the

[111] direction at temperature T1. 2

[111] Let E be the elastic modulus parameter along the

[111] direction at temperature T2. 1

[001] E is the elastic modulus parameter along the

[001] direction at temperature T1. 2

[001] Let be the elastic modulus parameter along the

[001] direction at temperature T2.

[0169]

[0170]

[0171] Calculate the direction vector [HKL] (material direction vector) of [hkl] (Cartesian direction vector) in the material coordinate system, where α, β and γ are the angular deviations between the Cartesian coordinate system and the material coordinate system.

[0172]

[0173] The elastic modulus parameter is calculated using the following formula. Where, E [hkl]Here, [HKL] is the elastic modulus parameter, and [HKL] is the direction vector of the measurement direction in the material coordinate system.

[111] E is the elastic modulus parameter along the

[111] direction.

[001] Here is the elastic modulus parameter along the

[001] direction.

[0174]

[0175] Step S1114: Calculate and output the vibration stress at the critical point (theoretical measurement location), using the following formula. Where σ is the vibration stress at the theoretical measurement location, ε is the strain at the actual measurement location, and E... [hkl] Let θ be the elastic modulus parameter along the measurement direction, and θ be the sensor's sensitivity parameter.

[0176]

[0177] In a method for measuring the vibration stress of turbine cooling blades provided in the exemplary embodiments of this disclosure, the sensitivity parameters of a sensor are obtained, the elastic modulus parameters of the blade along the measurement direction are determined based on the crystal orientation angle corresponding to the blade, and the strain at the actual measurement position is obtained by the sensor when the sensor is located at the actual measurement position and in the measurement direction. The stress at the theoretical measurement position is calculated by using the strain, the elastic modulus parameters, and the sensitivity parameters. On the one hand, considering the influence of the crystal orientation angle on the elastic modulus parameters, the elastic modulus parameters are determined based on the crystal orientation angle and added to the final vibration stress calculation process, reducing the measurement result error caused by the elastic modulus parameter error due to crystal orientation, and improving the measurement accuracy. On the other hand, since the effect of crystal orientation is considered in the calculation of each stress field, the influence of crystal orientation on the sensitivity parameters is considered, thus ensuring the accuracy of the sensitivity. Furthermore, considering the three factors of temperature, steady-state stress, and relative vibration stress, the theoretical measurement position can be accurately identified, providing more practical reference value.

[0178] Figure 12 This is a block diagram illustrating an apparatus suitable for measuring the vibration stress of turbine cooling blades according to an exemplary embodiment. (Refer to...) Figure 12 The device 1200 for measuring the vibration stress of turbine cooling blades includes a first parameter acquisition module 1210, a second parameter acquisition module 1220, a third parameter acquisition module 1230, and a calculation module 1240. Wherein:

[0179] The first parameter acquisition module 1210 is used to acquire the sensitivity parameter of the sensor; the second parameter acquisition module 1220 is used to determine the elastic modulus parameter of the blade along the measurement direction based on the crystal orientation angle corresponding to the blade; the third parameter acquisition module 1230 is used to acquire the strain at the actual measurement position through the sensor when the sensor is located at the actual measurement position and in the measurement direction; the calculation module x40 is used to calculate the stress at the theoretical measurement position through the strain, elastic modulus parameter and sensitivity parameter.

[0180] In an exemplary embodiment of this disclosure, based on the foregoing scheme, the apparatus further includes: a first vibration stress field calculation unit, used to calculate a first vibration stress field based on the resonance analysis results; a second vibration stress field calculation unit, used to calculate a second vibration stress field based on the load properties of the blade; and a theoretical measurement position determination unit, used to determine a theoretical measurement position on the blade based on the first vibration stress field and the second vibration stress field.

[0181] In an exemplary embodiment of this disclosure, based on the aforementioned scheme, the actual measurement position is determined according to the theoretical measurement position. The apparatus further includes: an equivalent stress calculation unit for calculating the equivalent stress at the theoretical measurement position; and a first actual measurement position determination unit for determining the actual measurement position based on the equivalent stress.

[0182] In an exemplary embodiment of this disclosure, based on the aforementioned scheme, the actual measurement position is determined according to the equivalent stress. The apparatus further includes: a candidate measurement position determination unit, configured to determine a plurality of candidate measurement positions that meet preset conditions according to the equivalent stress; and a second actual measurement position determination unit, configured to determine the actual measurement position with the largest normal stress among the plurality of candidate measurement positions.

[0183] In an exemplary embodiment of this disclosure, based on the aforementioned scheme, the device calculates a first vibration stress field according to the resonance analysis results. The device further includes: a first coordinate system establishment unit for establishing a Cartesian coordinate system; a second coordinate system establishment unit for acquiring the dendrite trunk of the crystal and determining a first material coordinate axis, a second material coordinate axis, and a third material coordinate axis based on the dendrite trunk, so that a material coordinate system is established based on the first material coordinate axis, the second material coordinate axis, and the third material coordinate axis; an angle deviation parameter acquisition unit for acquiring the angle deviation parameter between the Cartesian coordinate system and the material coordinate system; and a resonance analysis unit for performing resonance analysis on the blade using the angle deviation parameter to obtain the resonance analysis results.

[0184] In an exemplary embodiment of this disclosure, based on the foregoing scheme, the apparatus further includes: a stress component measurement unit for acquiring stress components at the actual measurement location; a candidate direction vector acquisition unit for determining multiple candidate direction vectors based on the actual measurement location and other locations on the blade besides the actual measurement location; and a target direction vector determination unit for calculating the normal stress on the blade surface along each candidate direction vector based on the stress components, and determining the candidate direction vector with the largest normal stress as the target direction vector.

[0185] In an exemplary embodiment of this disclosure, based on the aforementioned scheme, the elastic modulus parameter of the blade along the measurement direction is determined according to the crystal orientation angle corresponding to the blade. The device further includes: a first coordinate system establishment unit, used to establish a Cartesian coordinate system on the blade and determine the Cartesian direction vector of the measurement direction in the Cartesian coordinate system; a second coordinate system establishment unit, used to obtain the dendrite trunk of the crystal and determine a first material coordinate axis, a second material coordinate axis, and a third material coordinate axis according to the dendrite trunk, so that a material coordinate system is established according to the first material coordinate axis, the second material coordinate axis, and the third material coordinate axis; a temperature measurement unit, used to obtain the current temperature at the actual measurement position; an elastic modulus parameter acquisition unit, used to obtain the first elastic modulus parameter and the second elastic modulus parameter at the actual measurement position; a vector calculation unit, used to calculate the material direction vector according to the Cartesian direction vector; and an elastic modulus parameter calculation unit, used to calculate the elastic modulus parameter according to the current temperature, the first elastic modulus parameter, the second elastic modulus parameter, and the material direction vector.

[0186] The apparatus provided in this exemplary embodiment for measuring the vibration stress of turbine cooling blades acquires the sensitivity parameters of a sensor, determines the elastic modulus parameters of the blade along the measurement direction based on the crystal orientation angle corresponding to the blade, and acquires the strain at the actual measurement location when the sensor is located at the actual measurement position and in the measurement direction. The stress at the theoretical measurement location is then calculated using the strain, elastic modulus parameters, and sensitivity parameters. On one hand, considering the influence of crystal orientation angle on the elastic modulus parameters, the elastic modulus parameters are determined based on the crystal orientation angle and incorporated into the final vibration stress calculation, reducing measurement errors caused by elastic modulus parameter errors due to crystal orientation and improving measurement accuracy. On the other hand, since the effect of crystal orientation is considered in the calculation of each stress field, the influence of crystal orientation on the sensitivity parameters is considered, thus ensuring the accuracy of the sensitivity. Furthermore, considering the three factors of temperature, steady-state stress, and relative vibration stress, the theoretical measurement location can be accurately identified, providing greater practical reference value.

[0187] Since the functional modules of the vibration stress measurement and generation device in the example embodiments of this disclosure correspond to the steps of the example embodiments of the vibration stress measurement and generation method described above, for details not disclosed in the device embodiments of this disclosure, please refer to the embodiments of the vibration stress measurement and generation method described above.

[0188] It should be noted that although several modules or units for the device used to perform actions have been mentioned in the detailed description above, this division is not mandatory. In fact, according to embodiments of this disclosure, the features and functions of two or more modules or units described above can be embodied in one module or unit. Conversely, the features and functions of one module or unit described above can be further divided and embodied by multiple modules or units.

[0189] Figure 13 A schematic diagram of the structure of a computer system suitable for implementing the embodiments of the present disclosure is shown.

[0190] It should be noted that, Figure 13 The computer system 1300 of the electronic device shown is merely an example and should not impose any limitation on the functionality and scope of use of the embodiments disclosed herein.

[0191] like Figure 13 As shown, the computer system 1300 includes a central processing unit (CPU) 1301, which can perform various appropriate actions and processes based on programs stored in read-only memory (ROM) 1302 or programs loaded from storage section 1308 into random access memory (RAM) 1303. The RAM 1303 also stores various programs and data required for system operation. The CPU 1301, ROM 1302, and RAM 1303 are interconnected via a bus 1304. An I / O interface 1305 is also connected to the bus 1304.

[0192] The following components are connected to I / O interface 1305: an input section 1306 including a keyboard, mouse, etc.; an output section 1307 including a cathode ray tube (CRT), liquid crystal display (LCD), etc., and speakers, etc.; a storage section 1308 including a hard disk, etc.; and a communication section 1309 including a network interface card such as a LAN card, modem, etc. The communication section 1309 performs communication processing via a network such as the Internet. A drive 1310 is also connected to I / O interface 1305 as needed. Removable media 1311, such as a disk, optical disk, magneto-optical disk, semiconductor memory, etc., are installed on drive 1310 as needed so that computer programs read from them can be installed into storage section 1308 as needed.

[0193] In particular, according to embodiments of this disclosure, the processes described below with reference to the flowcharts can be implemented as computer software programs. For example, embodiments of this disclosure include a computer program product comprising a computer program carried on a computer-readable medium, the computer program containing program code for performing the methods shown in the flowcharts. In such embodiments, the computer program can be downloaded and installed from a network via communication section 1309, and / or installed from removable medium 1311. When the computer program is executed by central processing unit (CPU) 1301, it performs various functions defined in the methods and apparatus of this application. For example, it can perform functions such as... Figure 2 In step S210, the sensitivity parameter of the sensor is obtained; wherein, the sensitivity parameter is used to indicate the ratio of the stress on the theoretical measurement position and the actual measurement position on the blade under the same external factors, the theoretical measurement position is the position on the blade most affected by vibration, and the actual measurement position is determined based on the theoretical measurement position; in step S220, the elastic modulus parameter of the blade along the measurement direction is determined based on the crystal orientation angle corresponding to the blade; wherein, the measurement direction is the direction of the maximum normal stress at the actual measurement position; in step S230, when the sensor is located at the actual measurement position and in the measurement direction, the strain at the actual measurement position is obtained through the sensor; wherein, the strain is used to indicate the relative change in the shape of the blade under the action of external load, the external load is the product of stress and the area of ​​force application, and the stress is the internal force of the interaction within the blade; in step S240, the stress at the theoretical measurement position is calculated through the strain, elastic modulus parameter and sensitivity parameter.

[0194] In exemplary embodiments of this disclosure, a computer-readable storage medium is also provided, on which a program product capable of implementing the methods described above is stored. In some possible implementations, various aspects of this disclosure may also be implemented as a program product, including program code that, when the program product is run on a terminal, causes the terminal device to perform the various exemplary embodiments and steps according to this disclosure described in the "Exemplary Methods" section of this specification.

[0195] The program product for implementing the above-described method according to embodiments of the present disclosure may employ a portable compact disc read-only memory (CD-ROM) and include program code, and may run on a terminal device, such as a personal computer. However, the program product of the present disclosure is not limited thereto; the readable storage medium may be any tangible medium containing or storing a program that may be used by or in conjunction with an instruction execution system, apparatus, or device.

[0196] The program product may employ any combination of one or more readable media. A readable medium may be a readable signal medium or a readable storage medium. A readable medium may be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples (a non-exhaustive list) of readable storage media include: electrical connections having one or more wires, portable disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof.

[0197] Computer-readable signal media may include data signals propagated in baseband or as part of a carrier wave, carrying readable program code. Such propagated data signals may take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. A readable signal medium may also be any readable medium other than a readable storage medium, which can send, propagate, or transmit programs for use by or in conjunction with an instruction execution system, apparatus, or device.

[0198] The program code contained in the readable medium can be used for transmission on any suitable medium, including but not limited to wireless, wired, fiber optic, RF, or any suitable combination thereof.

[0199] Furthermore, the above figures are merely illustrative of the processes included in the method according to exemplary embodiments of this disclosure and are not intended to be limiting. It is readily understood that the processes shown in the above figures do not indicate or limit the temporal order of these processes. Additionally, it is readily understood that these processes may be executed synchronously or asynchronously, for example, in multiple modules.

[0200] Other embodiments of this disclosure will readily occur to those skilled in the art upon consideration of the specification and practice of the invention disclosed herein. This application is intended to cover any variations, uses, or adaptations of this disclosure that follow the general principles of this disclosure and include common knowledge or customary techniques in the art not disclosed herein. The specification and embodiments are to be considered exemplary only, and the true scope and spirit of this disclosure are indicated by the claims.

Claims

1. A method for measuring the vibration stress of turbine cooling blades, characterized in that, The method for measuring the stress on the blade using sensors includes: Acquire the sensor's sensitivity parameters; wherein the sensitivity parameters are used to indicate the ratio of stress experienced by a theoretical measurement position to that experienced by an actual measurement position on the blade under the same external factors, the theoretical measurement position being the position on the blade most affected by vibration, and the actual measurement position being determined based on the theoretical measurement position; wherein determining the actual measurement position includes: calculating the equivalent stress at the theoretical measurement position; determining multiple candidate measurement positions that meet preset conditions based on the equivalent stress; wherein the preset conditions are that the ratio of the normal stress to the equivalent stress at the candidate measurement position exceeds a preset threshold; and determining the actual measurement position with the highest normal stress among the multiple candidate measurement positions; The elastic modulus parameter of the blade along the measurement direction is determined based on the crystal orientation angle corresponding to the blade; the measurement direction is the direction of maximum normal stress at the actual measurement position; wherein, determining the elastic modulus parameter of the blade along the measurement direction includes: A Cartesian coordinate system is established on the blade, and the Cartesian direction vector of the measurement direction in the Cartesian coordinate system is determined; wherein, the first coordinate axis points to the blade axial direction, the second coordinate axis points to the blade circumferential direction, the third coordinate axis points to the blade radial direction, and the Cartesian direction vector is the direction vector of the measurement direction in the Cartesian coordinate system; Obtain the dendrite trunk of the crystal, and determine the first material coordinate axis, the second material coordinate axis and the third material coordinate axis according to the dendrite trunk, so as to establish a material coordinate system according to the first material coordinate axis, the second material coordinate axis and the third material coordinate axis; wherein, the first material coordinate axis and the second material coordinate axis are the two main directions of the secondary dendrite trunk, the third material coordinate axis is the main direction of the primary dendrite trunk, the direction of the primary dendrite trunk is the longitudinal preferred growth direction of the dendrite structure, which is the same as the [001] crystal direction, the direction of the secondary dendrite trunk is the transverse preferred growth direction of the dendrite structure, which is the same as the [010] and [100] crystal directions, and the [001] crystal direction is perpendicular to the plane where the [010] and [100] crystal directions are located; Obtain the current temperature at the actual measurement location; Obtain a first elastic modulus parameter and a second elastic modulus parameter at the actual measurement location; wherein, the first elastic modulus parameter is the elastic modulus parameter along a preset direction when the actual measurement location is at a first temperature, and the second elastic modulus parameter is the elastic modulus parameter along the preset direction when the actual measurement location is at a second temperature, and the first temperature, the second temperature and the current temperature are different from each other; The material direction vector is calculated based on the Cartesian direction vector; wherein, the material direction vector is the direction vector of the measurement direction in the material coordinate system; The elastic modulus parameter is calculated based on the current temperature, the first elastic modulus parameter, the second elastic modulus parameter, and the material direction vector; When the sensor is located at the actual measurement position and in the measurement direction, the strain at the actual measurement position is obtained through the sensor; wherein, the strain is used to indicate the relative change in the shape of the blade under the action of external load, the external load is the product of stress and the area of ​​force application, and the stress is the internal force of the interaction within the blade; The stress at the theoretical measurement location is calculated using the strain, the elastic modulus parameter, and the sensitivity parameter. The method further includes: The first vibration stress field is determined based on the resonance analysis results; wherein, the resonance analysis results are used to indicate the resonance properties of the blade, and the first vibration stress field is a relative vibration stress field, which is used to indicate the relative magnitude of the stress at each position of the blade when the blade resonates; The second vibration stress field is calculated based on the load-bearing properties of the blade; wherein, the second vibration stress field is an allowable vibration stress field, which is used to indicate the maximum value of vibration stress that each position on the blade can withstand without high-cycle fatigue failure. The theoretical measurement position is determined on the blade based on the first vibration stress field and the second vibration stress field.

2. The method according to claim 1, characterized in that, The calculation of the first vibration stress field based on the resonance analysis results includes: Obtain the angle deviation parameters between the Cartesian coordinate system and the material coordinate system; wherein, the angle deviation parameters are the angle deviation between the first coordinate axis and the first material coordinate axis, the angle deviation between the second coordinate axis and the second material coordinate axis, and the angle deviation between the third coordinate axis and the third material coordinate axis; The resonance analysis results are obtained by performing finite element analysis on the blade vibration characteristics and response using the aforementioned angle deviation parameters.

3. The method according to claim 1, characterized in that, The method further includes: Obtain the stress components at the actual measurement location; wherein, the stress components include the stress at the actual measurement location along the first coordinate axis, the stress along the second coordinate axis, and the stress along the third coordinate axis, and the first coordinate axis, the second coordinate axis, and the third coordinate axis are coordinate axes of a Cartesian coordinate system established with the actual measurement location as the origin; Multiple candidate direction vectors are determined based on the actual measurement position and other positions on the blade besides the actual measurement position; The normal stress on the blade surface along each of the candidate direction vectors is calculated based on the stress components, and the candidate direction vector with the largest normal stress is determined as the target direction vector; wherein, the target direction vector is used to indicate the measurement direction.

4. An apparatus suitable for measuring the vibration stress of turbine cooling blades, used to implement the method as described in claim 1, characterized in that, include: The first parameter acquisition module is used to acquire the sensitivity parameters of the sensor; wherein, the sensitivity parameters are used to indicate the ratio of the stress on the theoretical measurement position and the actual measurement position on the blade under the same external factors, the theoretical measurement position is the position on the blade that is most affected by vibration, and the actual measurement position is determined based on the theoretical measurement position; The second parameter acquisition module is used to determine the elastic modulus parameter of the blade along the measurement direction based on the crystal orientation angle corresponding to the blade; wherein, the measurement direction is the direction of maximum normal stress at the actual measurement position; The third parameter acquisition module is used to acquire the strain at the actual measurement position through the sensor when the sensor is located at the actual measurement position and in the measurement direction; wherein the strain is used to indicate the relative change in the shape of the blade under the action of external load, the external load is the product of stress and the area of ​​force application, and the stress is the internal force of the interaction inside the blade. The calculation module is used to calculate the stress at the theoretical measurement location using the strain, the elastic modulus parameter, and the sensitivity parameter.

5. An electronic device, characterized in that, include: processor; as well as Memory for storing the executable instructions of the processor; The processor is configured to execute the method for measuring the vibration stress of turbine cooling blades according to any one of claims 1-3 by executing the executable instructions.

6. A non-transitory computer-readable storage medium, characterized in that, When the instructions in the storage medium are executed by the processor of the electronic device, the electronic device is able to perform the method for measuring the vibration stress of turbine cooling blades as described in any one of claims 1-3.

Images