Thermoacoustic fatigue equivalent load spectrum induction method and device

Abstract

The invention provides a thermoacoustic fatigue load spectrum induction method and device, and belongs to the technical field of aircraft durability design, and the method comprises the steps: obtaining thermoacoustic fatigue load spectrums under different working conditions, and determining an equivalent value Tse of a structure surface static temperature according to the thermoacoustic fatigue load spectrums; dynamic fatigue S-N curves of the structural material at different temperatures are obtained, and equivalent characteristic parameters be and Ce of the dynamic fatigue S-N curves are determined according to the equivalent value of the static temperature of the surface of the structure; determining characteristic parameters and conversion parameters of the single-degree-of-freedom system; based on a fatigue damage spectrum theory of random vibration, fatigue damage spectrums caused by thermo-acoustic loads under different working conditions are calculated, and fatigue damage under the same inherent frequency under all working conditions is linearly superposed to obtain an accumulated damage spectrum; taking the accumulated damage spectrum as a target, and enabling a fatigue damage spectrum caused by the external force power spectrum density of the damage equivalent sound pressure to be equal to the accumulated damage spectrum by constructing the external force power spectrum density of the damage equivalent sound pressure, so as to obtain an equivalent fatigue load spectrum of the thermoacoustic load.

Description

Technical Field

[0001] This application belongs to the field of aircraft durability design technology, and specifically relates to a method and apparatus for summarizing the equivalent load spectrum of thermoacoustic fatigue. Background Technology

[0002] With the rapid development of aerospace technology, structural fatigue problems under thermoacoustic coupling environments are becoming increasingly prominent, especially for hypersonic vehicles. Due to the complexity of their flight environment, structures such as skin, engine casing, flame tube panels, and turbine blades are affected by high temperatures and strong noise. Fatigue life prediction under thermoacoustic loads is one of the important research directions in the design of such vehicles. Because hypersonic vehicles operate under diverse conditions, the temperature and noise environments experienced by the structure vary under different conditions; each condition corresponds to a thermoacoustic fatigue load spectrum. If this thermoacoustic fatigue load spectrum is directly used in fatigue analysis, it is necessary to calculate the damage caused by the thermoacoustic fatigue load spectrum under each condition and then perform cumulative damage assessment. This will greatly increase the workload of computational analysis and correspondingly increase the cycle and complexity of thermoacoustic fatigue testing.

[0003] To address this issue, engineering practice typically requires summarizing and simplifying the thermoacoustic fatigue load spectrum for each operating condition into a single thermoacoustic load spectrum equivalent to fatigue damage. Current technologies, based on Miner's linear cumulative damage criterion, "compress" the sound pressure power spectral density of different operating condition data segments into an equivalent power spectral density. However, this method does not consider the influence of temperature and cannot be used to summarize the equivalent thermoacoustic fatigue spectrum. Summary of the Invention

[0004] The purpose of this application is to provide a method and apparatus for summarizing the equivalent load spectrum of thermoacoustic fatigue, so as to solve or alleviate at least one problem in the background art.

[0005] On the one hand, the technical solution of this application is: a method for summarizing the equivalent load spectrum of thermoacoustic fatigue, comprising:

[0006] Thermoacoustic fatigue load spectra under different working conditions are obtained, and the equivalent value Ts of the static temperature of the structural surface is determined based on the thermoacoustic fatigue load spectra. e ;

[0007] Obtain the dynamic fatigue SN curves of the structural material at different temperatures, and determine the equivalent characteristic parameter b of the dynamic fatigue SN curve based on the equivalent value of the static temperature of the structural surface. e and C e ;

[0008] Determine the characteristic parameters and transformation parameters of a single-degree-of-freedom system;

[0009] Based on the fatigue damage spectrum theory of random vibration, the fatigue damage spectrum caused by thermoacoustic load under different working conditions is calculated, and the cumulative damage spectrum is obtained by linearly superimposing the fatigue damage under the same natural frequency for all working conditions.

[0010] By taking the cumulative damage spectrum as the target, and constructing the external force power spectral density of the damage equivalent sound pressure, the fatigue damage spectrum caused by the external force power spectral density of the damage equivalent sound pressure is equal to the cumulative damage spectrum, thus obtaining the equivalent fatigue load spectrum of the thermoacoustic load.

[0011] In at least one embodiment of this application, the thermoacoustic fatigue load spectrum includes the action time, the static temperature of the structural surface, and the sound pressure power spectral density, wherein the equivalent value of the static temperature of the structural surface is the static temperature of the structural surface corresponding to the maximum action time in the thermoacoustic fatigue load spectrum.

[0012] In at least one embodiment of this application, the following is employed: Characterizing the SN curve, where S p Let N be the stress amplitude, N be the number of fatigue cycles, and b and C be characteristic parameters. The equivalent characteristic parameter b of the dynamic fatigue SN curve is... e and C e Select the equivalent value Ts of the static temperature of the structure surface. e The parameters b and C of the corresponding dynamic fatigue SN curve.

[0013] In at least one embodiment of this application, the characteristic parameters and conversion parameters of the single-degree-of-freedom system include the system mass m, the system critical damping ratio ζ, the sound pressure area A, and the displacement-to-stress conversion factor K, wherein the system mass m is a positive value within 20, the system critical damping ratio ζ is a positive value within 0.03, the sound pressure area A is a positive value within 1, and the displacement-to-stress conversion factor K is a positive value within 1000.

[0014] In at least one embodiment of this application, the process of calculating the fatigue damage spectrum caused by thermoacoustic load under different working conditions based on the fatigue damage spectrum theory of random vibration is as follows:

[0015] 1) Extract the sound pressure power spectral density Φ under working condition i p,i (f);

[0016] 2) The sound pressure power spectral density Φ p,i (f) Multiply by the square of the area A where the sound pressure is applied to obtain the power spectral density Φ of the dynamic force. F,i (f), that is: Φ F,i (f)=A 2 ·Φ p,i (f);

[0017] 3) Set the natural frequency f of a single-degree-of-freedom system j. n,jThe power spectral density Φ of the dynamic force is calculated based on the displacement response calculation formula of a single-degree-of-freedom system under the action of external force power spectral density. F,i (f) Displacement response caused by:

[0018]

[0019] In the formula: H z,j (f) is the transfer function of the j-response displacement z and dynamic force F of a single-degree-of-freedom system; m j Let f be the mass of the single-degree-of-freedom system j; n,j ζ is the natural frequency of the single-degree-of-freedom system j; j Let be the critical damping ratio of a single-degree-of-freedom system j;

[0020] 4) Calculate the mean displacement response per unit time and the number of positive crossings. and its root mean square value z rms,j :

[0021]

[0022]

[0023] In the formula: f c For displacement response Φ z,ij (f) The corresponding cutoff frequency;

[0024] 5) Since the response of a single-degree-of-freedom system is a narrow-band process, its peak probability density function tends to a Rayleigh distribution. Therefore, fatigue damage D j Represented as:

[0025] In the formula: K is the displacement to stress conversion factor; T i Γ is the duration of the thermoacoustic load; Γ is the gamma function; C i Let b be the material fatigue performance constant under the i-th working condition; i It is the negative reciprocal of the slope of the material fatigue SN curve under the i-th working condition;

[0026] 6) Select the natural frequency f of the next single-degree-of-freedom system. n,j+1 The corresponding sound pressure power spectral density Φ p,i (f) For the next center frequency in the frequency series, repeat steps 3 to 5 above to obtain the fatigue damage D at that natural frequency. j+1 ;

[0027] 7) Repeat step 6 until the natural frequency reaches the cutoff frequency f. c Thus, a set (f) is obtained. n -D) i Data, namely the fatigue damage spectrum of working condition i;

[0028] 8) Based on the calculation process of the fatigue damage spectrum under the above single working condition, the fatigue damage spectrum (f) under different working conditions is obtained. n -D) i+1 、(f n -D) i+2 ... (f n -D) num , where num is the total number of operating conditions.

[0029] In at least one embodiment of this application, the process of constructing an external force power spectral density of the damage equivalent sound pressure, such that the fatigue damage spectrum caused by the external force power spectral density of the damage equivalent sound pressure is equal to the cumulative damage spectrum, is as follows:

[0030] 1) Construct a white noise spectrum of dynamic pressure, with the same frequency range as the fatigue cumulative damage spectrum, and an amplitude of 1 Pa. 2 / Hz;

[0031] 2) Set the natural frequency f of a single-degree-of-freedom system. n,j Calculate the fatigue damage D caused by white noise spectrum. j In the calculation of fatigue damage caused by white noise spectrum, the duration of thermoacoustic load is taken as the total duration of thermoacoustic load under all working conditions, i.e., the duration of equivalent fatigue load spectrum T. e ;

[0032] 3) Extract the natural frequency f n,j The sum of the corresponding cumulative damage spectrum amplitude and fatigue loss (∑D) i ) j Calculate the sum of fatigue losses using the following formula (∑D) i ) j With fatigue damage D j The ratio Am j ;

[0033]

[0034] In the formula: b e The negative reciprocal of the slope of the equivalent material fatigue SN curve, i.e., the equivalent characteristic parameter b. e ;

[0035] 4) The ratio Am j As the natural frequency f n,j The amplitude of the corresponding damage equivalent sound pressure external force power spectral density;

[0036] 5) By changing the natural frequency of the single-degree-of-freedom system, repeat steps 2 to 4 to obtain the amplitude of the external force power spectral density of the damage equivalent sound pressure at different natural frequencies, i.e., f. n -Am.

[0037] In at least one embodiment of this application, the method further includes: verifying the validity of the equivalent thermoacoustic fatigue equivalent load spectrum summarization results, the process of which includes:

[0038] Using the equivalent fatigue load spectrum of the thermoacoustic load as the thermoacoustic load, the fatigue damage spectrum caused by the thermoacoustic load is calculated according to the calculation process of the fatigue damage spectrum caused by the thermoacoustic load. The fatigue damage spectrum and the cumulative damage spectrum are plotted on one graph to check the degree of agreement between the two curves.

[0039] The relative error of fatigue damage at the same natural frequency is calculated based on the data from the two curves, and the validity of the induction results of the equivalent thermoacoustic fatigue equivalent load spectrum is verified based on the relative error.

[0040] On the other hand, the technical solution provided in this application is: a thermoacoustic fatigue equivalent load spectrum induction device, comprising:

[0041] The load spectrum acquisition module is used to acquire thermoacoustic fatigue load spectra under different working conditions, and determine the equivalent value Ts of the static temperature of the structural surface based on the thermoacoustic fatigue load spectra. e ;

[0042] The SN curve acquisition module is used to acquire the dynamic fatigue SN curves of structural materials at different temperatures, and to determine the equivalent characteristic parameter b of the dynamic fatigue SN curve based on the equivalent value of the static temperature of the structural surface. e and C e ;

[0043] The parameter determination module determines the characteristic parameters and transformation parameters of a single-degree-of-freedom system.

[0044] The damage spectrum calculation module is used to calculate the fatigue damage spectrum caused by thermoacoustic load under different working conditions based on the fatigue damage spectrum theory of random vibration. The cumulative damage spectrum is obtained by linearly superimposing the fatigue damage at the same natural frequency under all working conditions.

[0045] The induction module takes the cumulative damage spectrum as the target and constructs the external force power spectral density of the damage equivalent sound pressure so that the fatigue damage spectrum caused by the external force power spectral density of the damage equivalent sound pressure is equal to the cumulative damage spectrum, thereby obtaining the equivalent fatigue load spectrum of the thermoacoustic load.

[0046] In at least one embodiment of this application, the thermoacoustic fatigue load spectrum includes the action time, the static temperature of the structural surface, and the sound pressure power spectral density, wherein the equivalent value of the static temperature of the structural surface is the static temperature of the structural surface corresponding to the maximum action time in the thermoacoustic fatigue load spectrum.

[0047] In at least one embodiment of this application, the following is employed: Characterizing the SN curve, where S pLet N be the stress amplitude, N be the number of fatigue cycles, and b and C be characteristic parameters. The equivalent characteristic parameter b of the dynamic fatigue SN curve is... e and C e Select the equivalent value Ts of the static temperature of the structure surface. e The parameters b and C of the corresponding dynamic fatigue SN curve.

[0048] In at least one embodiment of this application, the characteristic parameters and conversion parameters of the single-degree-of-freedom system include the system mass m, the system critical damping ratio ζ, the sound pressure area A, and the displacement-to-stress conversion factor K, wherein the system mass m is a positive value within 20, the system critical damping ratio ζ is a positive value within 0.03, the sound pressure area A is a positive value within 1, and the displacement-to-stress conversion factor K is a positive value within 1000.

[0049] In at least one embodiment of this application, the process of calculating the fatigue damage spectrum caused by thermoacoustic load under different working conditions based on the fatigue damage spectrum theory of random vibration is as follows:

[0050] 1) Extract the sound pressure power spectral density Φ under working condition i p,i (f);

[0051] 2) The sound pressure power spectral density Φ p,i (f) Multiply by the square of the area A where the sound pressure is applied to obtain the power spectral density Φ of the dynamic force. F,i (f), that is: Φ F,i (f)=A 2 ·Φ p,i (f);

[0052] 3) Set the natural frequency f of a single-degree-of-freedom system j. n,j The power spectral density Φ of the dynamic force is calculated based on the displacement response calculation formula of a single-degree-of-freedom system under the action of external force power spectral density. F,i (f) Displacement response caused by:

[0053]

[0054] In the formula: H z,j (f) is the transfer function of the j-response displacement z and dynamic force F of a single-degree-of-freedom system; m j Let f be the mass of the single-degree-of-freedom system j; n,j ζ is the natural frequency of the single-degree-of-freedom system j; j Let be the critical damping ratio of a single-degree-of-freedom system j;

[0055] 4) Calculate the mean displacement response per unit time and the number of positive crossings. and its root mean square value z rms,j :

[0056]

[0057]

[0058] In the formula: f c For displacement response Φ z,ij (f) The corresponding cutoff frequency;

[0059] 5) Since the response of a single-degree-of-freedom system is a narrow-band process, its peak probability density function tends to a Rayleigh distribution. Therefore, fatigue damage D j Represented as:

[0060] In the formula: K is the displacement to stress conversion factor; T i Γ is the duration of the thermoacoustic load; Γ is the gamma function; C i Let b be the material fatigue performance constant under the i-th working condition; i It is the negative reciprocal of the slope of the material fatigue SN curve under the i-th working condition;

[0061] 6) Select the natural frequency f of the next single-degree-of-freedom system. n,j+1 The corresponding sound pressure power spectral density Φ p,i (f) For the next center frequency in the frequency series, repeat steps 3 to 5 above to obtain the fatigue damage D at that natural frequency. j+1 ;

[0062] 7) Repeat step 6 until the natural frequency reaches the cutoff frequency f. c Thus, a set (f) is obtained. n -D) i Data, namely the fatigue damage spectrum of working condition i;

[0063] 8) Based on the calculation process of the fatigue damage spectrum under the above single working condition, the fatigue damage spectrum (f) under different working conditions is obtained. n -D) i+1 、(f n -D) i+2 ... (f n -D) num , where num is the total number of operating conditions.

[0064] In at least one embodiment of this application, the process of constructing an external force power spectral density of the damage equivalent sound pressure, such that the fatigue damage spectrum caused by the external force power spectral density of the damage equivalent sound pressure is equal to the cumulative damage spectrum, is as follows:

[0065] 1) Construct a white noise spectrum of dynamic pressure, with the same frequency range as the fatigue cumulative damage spectrum, and an amplitude of 1 Pa. 2 / Hz;

[0066] 2) Set the natural frequency f of a single-degree-of-freedom system. n,j Calculate the fatigue damage D caused by white noise spectrum. j In the calculation of fatigue damage caused by white noise spectrum, the duration of thermoacoustic load is taken as the total duration of thermoacoustic load under all working conditions, i.e., the duration of equivalent fatigue load spectrum T. e ;

[0067] 3) Extract the natural frequency f n,j The sum of the corresponding cumulative damage spectrum amplitude and fatigue loss (∑D) i ) j Calculate the sum of fatigue losses using the following formula (∑D) i ) j With fatigue damage D j The ratio Am j ;

[0068]

[0069] In the formula: b e The negative reciprocal of the slope of the equivalent material fatigue SN curve, i.e., the equivalent characteristic parameter b. e ;

[0070] 4) The ratio Am j As the natural frequency f n,j The amplitude of the corresponding damage equivalent sound pressure external force power spectral density;

[0071] 5) By changing the natural frequency of the single-degree-of-freedom system, repeat steps 2 to 4 to obtain the amplitude of the external force power spectral density of the damage equivalent sound pressure at different natural frequencies, i.e., f. n -Am.

[0072] In at least one embodiment of this application, it further includes: a verification module for verifying the validity of the equivalent thermoacoustic fatigue equivalent load spectrum summarization results, the process of which includes:

[0073] Using the equivalent fatigue load spectrum of the thermoacoustic load as the thermoacoustic load, the fatigue damage spectrum caused by the thermoacoustic load is calculated according to the calculation process of the fatigue damage spectrum caused by the thermoacoustic load. The fatigue damage spectrum and the cumulative damage spectrum are plotted on one graph to check the degree of agreement between the two curves.

[0074] The relative error of fatigue damage at the same natural frequency is calculated based on the data from the two curves, and the validity of the induction results of the equivalent thermoacoustic fatigue equivalent load spectrum is verified based on the relative error.

[0075] Thirdly, this application provides an electronic device, comprising:

[0076] One or more processors;

[0077] Memory;

[0078] One or more applications, which are stored in the memory and configured to be executed by the one or more processors, are configured to implement the thermoacoustic fatigue equivalent load spectrum induction method as described in any of the preceding claims.

[0079] Finally, this application provides a computer-readable storage medium storing at least one instruction, at least one program, code set, or instruction set, wherein the at least one instruction, the at least one program, the code set, or the instruction set is loaded and executed by a processor to implement the thermoacoustic fatigue equivalent load spectrum induction method as described in any of the preceding claims.

[0080] The method and apparatus of this application can reasonably and effectively determine the equivalent fatigue load spectrum of thermoacoustic loads under different working conditions, and help improve the efficiency of analysis and testing, which has important guiding significance for the thermoacoustic fatigue analysis and testing of structures. Attached Figure Description

[0081] To more clearly illustrate the technical solutions provided in this application, the accompanying drawings will be briefly described below. Obviously, the drawings described below are merely some embodiments of this application.

[0082] Figure 1 This is a schematic diagram of the thermoacoustic fatigue equivalent spectrum induction method of this application.

[0083] Figure 2 This is the sound pressure power spectral density of operating condition 1 in an embodiment of this application.

[0084] Figure 3 The sound pressure power spectral density is for operating condition 2 in an embodiment of this application.

[0085] Figure 4 The sound pressure power spectral density is for operating condition 3 in an embodiment of this application.

[0086] Figure 5 The sound pressure power spectral density is for operating condition 4 in an embodiment of this application.

[0087] Figure 6 This is a fatigue damage spectrum under different working conditions according to an embodiment of this application.

[0088] Figure 7 This is the cumulative damage spectrum fn-∑Di of an embodiment of this application.

[0089] Figure 8 This is an embodiment of the thermoacoustic fatigue equivalent sound pressure spectrum of this application.

[0090] Figure 9This is a comparison diagram of the fatigue damage spectrum and the cumulative damage spectrum caused by the equivalent sound pressure spectrum according to an embodiment of this application.

[0091] Figure 10 This is a schematic diagram of the relative error curve of fatigue damage according to an embodiment of this application.

[0092] Figure 11 This is a schematic diagram of the thermoacoustic fatigue equivalent spectrum summarization device of this application. Detailed Implementation

[0093] To make the objectives, technical solutions, and advantages of this application clearer, the technical solutions in the embodiments of this application will be described in more detail below with reference to the accompanying drawings.

[0094] To address the problem that existing methods for summarizing equivalent acoustic fatigue load spectra cannot be applied to the equivalent processing of thermoacoustic fatigue load spectra, this application provides a method for summarizing multiple equivalent thermoacoustic fatigue load spectra under different working conditions into a single thermoacoustic fatigue load spectrum, thereby improving the efficiency of structural thermoacoustic fatigue analysis and testing.

[0095] like Figure 1 As shown, the thermoacoustic fatigue equivalent load spectrum induction method provided in this application includes:

[0096] Step S10: Obtain the thermoacoustic fatigue load spectrum under different working conditions.

[0097] In engineering practice, thermal and acoustic fatigue loads on structures of interest (usually important load-bearing structures) under different working conditions can typically be obtained through simulation or experimentation. After data processing, a thermoacoustic fatigue load spectrum can be obtained, as shown in Table 1. The dimensions of the action time T can be seconds, minutes, or hours; the static temperature of the structural surface Ts is typically taken as degrees Celsius (°C); and the sound pressure power spectral density Φ... p (f) Contains two columns of data. The first column is the frequency, taken as the center frequency of 1 / 3 octave band, in Hz. The second column is the amplitude, in Pa. 2 / Hz. Additionally, the equivalent value of the static temperature Ts of the structural surface needs to be determined. e In this application, the static temperature of the structural surface corresponding to the maximum action time max{T1,T2,T3,…} is preferably selected as the equivalent value Ts of the static temperature of the structural surface. e .

[0098] Table 1 Schematic diagram of structural thermoacoustic fatigue load spectrum

[0099] Operating condition number Duration of action static temperature of structural surface Sound pressure power spectral density 1 <![CDATA[T1]]> <![CDATA[Ts1]]> <![CDATA[Φ p,1 (f)]]> 2 <![CDATA[T2]]> <![CDATA[Ts2]]> <![CDATA[Φ p,2 (f)]]> 3 <![CDATA[T3]]> <![CDATA[Ts3]]> <![CDATA[Φ p,3 (f)]]> … … … …

[0100] As shown in Table 2, the structural thermoacoustic fatigue load spectrum provided in this embodiment of the application, wherein the sound pressure power spectral density Φ p(f) The horizontal axis represents the center frequency of 1 / 3 octave band, with a frequency range of 5Hz to 2500Hz, and the dimension is Hz; the vertical axis represents the amplitude, and the dimension is Pa. 2 / Hz. Since operating condition 3 has the longest duration, the corresponding static temperature of the structural surface is taken as the equivalent value of the static temperature of the structural surface, i.e., Ts. e =190℃.

[0101] Table 2 Structural thermoacoustic fatigue load spectrum

[0102] Operating condition number Action time / hour Static temperature of structural surface / °C Sound pressure power spectral density 1 2619 60 See Figure 2 As shown 2 987 170 See Figure 3 As shown 3 3629 190 See Figure 4 As shown 4 765 260 See Figure 5 As shown

[0103] Step 20: Obtain the dynamic fatigue SN curves of the structural material at different temperatures.

[0104] In engineering practice, the dynamic fatigue SN curves (Sf, St, and Ss) of structural materials at different temperatures can typically be obtained through thermal vibration fatigue testing. p (where N is the stress amplitude and N is the number of fatigue cycles) The Basquin formula is preferred in this application. The SN curves are characterized by dynamic fatigue. Through data fitting and interpolation methods, the characteristic parameters b and C corresponding to the dynamic fatigue SN curves at static temperatures for different structural surfaces in Table 1 can be obtained. The specific data format is shown in Table 3. Additionally, the equivalent characteristic parameter b of the dynamic fatigue SN curve needs to be determined. e and C e In this application, the equivalent value Ts of the static temperature of the structural surface is preferably adopted. e The corresponding dynamic fatigue SN curve parameters b and C are used as the equivalent characteristic parameters b of the dynamic fatigue SN curve. e and C e .

[0105] Table 3. Schematic diagram of characteristic parameters of fatigue SN curves at different temperatures.

[0106] Serial Number static temperature of structural surface Parameter b Parameter C 1 <![CDATA[Ts1]]> <![CDATA[b1]]> <![CDATA[C1]]> 2 <![CDATA[Ts2]]> <![CDATA[b2]]> <![CDATA[C2]]> 3 <![CDATA[Ts3]]> <![CDATA[b3]]> <![CDATA[C3]]> … … … …

[0107] Table 4 shows the fatigue SN curve characteristic parameters at different temperatures in this embodiment of the application. Based on the equivalent value Ts of the static temperature of the structural surface in step S10... e =190℃, therefore, the characteristic parameter of the dynamic fatigue SN curve corresponding to row 3 in Table 4 is taken as the equivalent characteristic parameter b of the fatigue SN curve. e and C e That is, b e =5.2327、C e =2.6443×10 16 .

[0108] Table 4 Characteristic parameters of fatigue SN curves at different temperatures

[0109] Serial Number Static temperature of structural surface / °C Parameter b Parameter C 1 60 4.8451 <![CDATA[1.2845×10 16 ]]> 2 170 5.1851 <![CDATA[2.4296×10 16 ]]> 3 190 5.2327 <![CDATA[2.6443×10 16 ]]> 4 260 7.9801 <![CDATA[4.4672×10 20 ]]>

[0110] Step S30: Determine the characteristic parameters and transformation parameters of the single-degree-of-freedom system.

[0111] To calculate the fatigue damage caused by the thermoacoustic fatigue load spectrum, it is also necessary to determine the characteristic parameters and transformation parameters of the single-degree-of-freedom system, including the system mass m, the system critical damping ratio ζ, the sound pressure area A, and the displacement-to-stress conversion factor K, as detailed in Table 5. In this application, the preferred recommended values ​​are: system mass m = positive values ​​within 20; system critical damping ratio ζ = positive values ​​within 0.03; sound pressure area A = positive values ​​within 1; and displacement-to-stress conversion factor K = positive values ​​within 1000.

[0112] Table 5 Characteristic parameters and transformation parameters of a single-degree-of-freedom system

[0113] Parameter name Parameter symbol Parameter dimensions Parameter values System quality m kg Pending System critical damping ratio ζ Dimensionless Pending sound pressure area A <![CDATA[m 2 ]]> Pending Displacement to stress conversion factor K MPa / m Pending

[0114] For example, in this embodiment of the application, based on the actual usage of the structure and engineering experience, the characteristic parameters and transformation parameters of the single-degree-of-freedom system are determined as shown in Table 7.

[0115] Table 6 Characteristic parameters and transformation parameters of a single-degree-of-freedom system

[0116] Parameter name Parameter symbol Parameter values System quality m 10 kg System critical damping ratio ζ 0.025 sound pressure area A <![CDATA[0.05 m 2 ]]> Displacement to stress conversion factor K 100 MPa / m

[0117] Step S40: Calculate the cumulative damage spectrum of thermoacoustic fatigue load spectrum under different working conditions.

[0118] This application summarizes the equivalent spectrum of thermoacoustic fatigue based on the fatigue damage spectrum (FDS) theory of random vibration. This theory states that the fatigue damage spectrum of random vibration is a series of curves showing the change of fatigue damage caused by random excitation in a single-degree-of-freedom linear system under a given damping ratio with the natural frequency of the single-degree-of-freedom system.

[0119] FDS theory is based on the following assumptions:

[0120] 1) Simplify the structure into a series of single-degree-of-freedom linear systems;

[0121] 2) The Basquin formula is used to characterize the fatigue SN curve of the material, i.e. ;

[0122] 3) It is assumed that the stress amplitude S p With the maximum system displacement z p Proportional, i.e., S p =K·z p ;

[0123] 4) The number of peak response times was counted using the rainflow method;

[0124] 5) The Miner criterion is used to characterize cumulative damage.

[0125] According to FDS theory, the calculation process of fatigue damage spectrum caused by thermoacoustic load under a certain working condition is as follows:

[0126] 1) Extract the power spectral density Φ of the sound pressure under condition i. p,i (f);

[0127] 2) The power spectral density Φ p,i (f) Multiply by the square of the area A where the sound pressure is applied to obtain the power spectral density Φ of the dynamic force. F,i (f), that is: Φ F,i (f)=A 2 ·Φ p,i (f). (1)

[0128] 3) Set the natural frequency f of a single-degree-of-freedom system j. n,j (Usually starting from 5Hz), calculate the power spectral density Φ of the dynamic force based on the displacement response calculation formula of a single-degree-of-freedom system under the action of external force power spectral density (PSD). F,i (f) Displacement response caused by:

[0129] (2)

[0130] In the formula: H z,j (f) is the transfer function of the j-response displacement z and dynamic force F of a single-degree-of-freedom system; m j Let f be the mass of the single-degree-of-freedom system j; n,j ζ is the natural frequency of the single-degree-of-freedom system j; j Let be the critical damping ratio of a single-degree-of-freedom system j.

[0131] 4) Calculate the mean displacement response per unit time and the number of positive crossings. and its root mean square value z rms,j :

[0132] (3)

[0133] (4)

[0134] In the formula: f c For displacement response Φ z,ij (f) The corresponding cutoff frequency, for example, f in this embodiment of the application. c =2500Hz.

[0135] 5) Since the response of a single-degree-of-freedom system is a narrow-band process, its peak probability density function tends to a Rayleigh distribution. Therefore, fatigue damage D jIt can be represented as: (5)

[0136] In the formula: K is the displacement to stress conversion factor; T i Γ is the duration of the thermoacoustic load; Γ is the gamma function; C i Let b be the material fatigue performance constant under the i-th working condition; i It is the negative reciprocal of the slope of the material fatigue SN curve under the i-th working condition.

[0137] 6) Select the natural frequency f of the next single-degree-of-freedom system. n,j+1 , corresponding to Φ p,i (f) Repeat steps 3 to 5 above for the next center frequency in the frequency series (e.g., 6.3Hz after 5Hz, and so on) to obtain the fatigue damage D at that natural frequency. j+1 ;

[0138] 7) Repeat step 6 until the natural frequency reaches the cutoff frequency f. c This gives us a set (f) n -D) i The data refers to the fatigue damage spectrum of working condition i.

[0139] 8) Based on the calculation method and steps for the fatigue damage spectrum under a single working condition described above, the fatigue damage spectrum (f) under different working conditions can be obtained. n -D) i+1 、(f n -D) i+2 ... (f n -D) num , where num is the total number of operating conditions.

[0140] like Figure 6 The figure shows the fatigue damage spectrum (f) under the above four working conditions in this embodiment of the application. n -D)1、(f n -D)2、(f n -D)4、(f n -D)4.

[0141] By linearly superimposing the fatigue damage under the same natural frequency for all operating conditions, the cumulative damage spectrum f can be obtained. n -∑D i ,like Figure 7 As shown.

[0142] Step S50: Summarize the equivalent fatigue load spectrum.

[0143] The cumulative damage spectrum f n -∑D iAs a goal, the external force power spectral density (PSD) of the damage equivalent sound pressure is constructed such that the fatigue damage spectrum induced by it is equal to the cumulative damage spectrum. The specific construction process includes:

[0144] 1) Construct a white noise spectrum of dynamic pressure, with the same frequency range as the fatigue cumulative damage spectrum, and an amplitude of 1 Pa. 2 / Hz;

[0145] 2) Set the natural frequency f of a single-degree-of-freedom system. n,j The fatigue damage D caused by the white noise spectrum is calculated according to equations (1) to (5). j Note that in equation (5) T i Take the total duration of the thermoacoustic load under all operating conditions, i.e., the duration of the equivalent fatigue load spectrum, T. e For example, in this embodiment, the effective fatigue load spectrum's duration T e =2619+987+3629+765=8000 hours;

[0146] 3) Extract the natural frequency f n,j The sum of the corresponding cumulative damage spectrum amplitude and fatigue loss (∑D) i ) j Calculate the sum of fatigue losses (∑D) according to formula (6). i ) j With fatigue damage D j The ratio Am j ;

[0147] (6)

[0148] In the formula: b e The negative reciprocal of the slope of the equivalent material fatigue SN curve (i.e., the equivalent characteristic parameter b) e ).

[0149] 4) Am j As the natural frequency f n,j The amplitude of the corresponding damage equivalent sound pressure level (PSD);

[0150] 5) By changing the natural frequency of the single-degree-of-freedom system and repeating steps 2 to 4, the amplitude of the damage equivalent sound pressure (PSD) at different natural frequencies can be obtained, i.e., f. n -Am.

[0151] After the above processing, the equivalent fatigue load spectrum of the thermoacoustic load is finally obtained as shown in Table 7. The equivalent fatigue load spectrum of the thermoacoustic load obtained in the embodiments of this application is shown in Table 8.

[0152] Table 7. Schematic diagram of equivalent fatigue load spectrum for thermoacoustic load.

[0153] Duration of action static temperature of structural surface Sound pressure power spectral density <![CDATA[T e ]]> <![CDATA[Ts e ]]> <![CDATA[f n -Am]]>

[0154] Table 8. Schematic diagram of equivalent thermoacoustic fatigue load spectrum

[0155] Action time / hour Static temperature of structural surface / °C Sound pressure power spectral density 8000 190 See Figure 8

[0156] Step S60: Verify the validity of the induction results of the equivalent thermoacoustic fatigue equivalent load spectrum.

[0157] Will Figure 8 The thermoacoustic fatigue equivalent sound pressure spectrum shown is used as a thermoacoustic load. The fatigue damage spectrum caused by this load is calculated according to steps 1 to 7 in step S40. This fatigue damage spectrum is then compared with... Figure 7 The cumulative damage spectrum shown is plotted on a single graph, as follows: Figure 9 As shown, you can check the degree of agreement between the two curves. Additionally, according to... Figure 9 The data from the two curves can be used to calculate the relative error of fatigue damage at the same natural frequency. (See...) Figure 10 As can be seen, the absolute value of the relative error of fatigue damage across the entire frequency band is within 8%, indicating that the method of this application is feasible and effective.

[0158] like Figure 11 As shown, based on this, this application also provides a thermoacoustic fatigue equivalent load spectrum induction device, the device 100 comprising:

[0159] The load spectrum acquisition module 101 is used to acquire the thermoacoustic fatigue load spectrum under different working conditions, and determine the equivalent value Ts of the static temperature of the structural surface based on the thermoacoustic fatigue load spectrum. e ;

[0160] The SN curve acquisition module 102 is used to acquire the dynamic fatigue SN curves of structural materials at different temperatures, and to determine the equivalent characteristic parameter b of the dynamic fatigue SN curve based on the equivalent value of the static temperature of the structural surface. e and C e ;

[0161] Parameter determination module 103 determines the characteristic parameters and transformation parameters of a single-degree-of-freedom system;

[0162] The damage spectrum calculation module 104 is used to calculate the fatigue damage spectrum caused by thermoacoustic load under different working conditions based on the fatigue damage spectrum theory of random vibration, and to obtain the cumulative damage spectrum by linearly superimposing the fatigue damage under the same natural frequency of all working conditions.

[0163] The summarization module 105 takes the cumulative damage spectrum as the target and constructs the external force power spectral density of the damage equivalent sound pressure so that the fatigue damage spectrum caused by the external force power spectral density of the damage equivalent sound pressure is equal to the cumulative damage spectrum, thereby obtaining the equivalent fatigue load spectrum of the thermoacoustic load.

[0164] In addition, the device 100 also includes a verification module 106 for verifying the validity of the equivalent thermoacoustic fatigue equivalent load spectrum summarization results.

[0165] The processing procedures of each module in the device 100 of this application can be referred to the above-mentioned thermoacoustic fatigue equivalent load spectrum induction method, and will not be repeated here.

[0166] In addition, this application also provides an electronic device comprising: one or more processors;

[0167] Memory;

[0168] One or more applications, which are stored in the memory and configured to be executed by the one or more processors, are configured to implement the thermoacoustic fatigue equivalent load spectrum induction method as described in any of the preceding claims.

[0169] Finally, this application also provides a computer-readable storage medium, characterized in that the computer-readable storage medium stores at least one instruction, at least one program, code set, or instruction set, wherein the at least one instruction, the at least one program, the code set, or the instruction set is loaded and executed by a processor to implement the thermoacoustic fatigue equivalent load spectrum induction method described above.

[0170] The above description is merely a specific embodiment of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method of inducing a thermal-mechanical fatigue equivalent load spectrum, characterized in that, The method comprises the following steps: Obtain a thermal-acoustic fatigue load spectrum under different working conditions, and determine an equivalent value Ts of a structure surface static temperature according to the thermal-acoustic fatigue load spectrum e ; The dynamic fatigue S-N curve of the structural material at different temperatures is obtained, and the equivalent characteristic parameter b of the dynamic fatigue S-N curve is determined according to the equivalent value of the static temperature of the structure surface e and C e ; determining characteristic parameters and conversion parameters of a single degree of freedom system; calculating fatigue damage spectra caused by thermoacoustic loads under different working conditions based on fatigue damage spectrum theory of random vibration, linearly superimposing fatigue damages under the same natural frequency of all working conditions to obtain a cumulative damage spectrum; taking the cumulative damage spectrum as a target, constructing an external force power spectrum density of damage equivalent sound pressure, making fatigue damage spectrum caused by the external force power spectrum density of damage equivalent sound pressure equal to the cumulative damage spectrum, and thus obtaining an equivalent fatigue load spectrum of the thermoacoustic load.

2. The thermo-mechanical fatigue equivalent load spectrum condensation method of claim 1, wherein, The thermoacoustic fatigue load spectrum comprises an action time, a structure surface static temperature and a sound pressure power spectrum density, and an equivalent value of the structure surface static temperature is a structure surface static temperature corresponding to a maximum value of the action time in the thermoacoustic fatigue load spectrum.

3. The thermo-mechanical fatigue equivalent load spectrum condensation method of claim 2, wherein, use Characterizing the SN curve, where S p Let N be the stress amplitude, N be the number of fatigue cycles, and b and C be characteristic parameters. The equivalent characteristic parameter b of the dynamic fatigue SN curve is... e and C e Select the equivalent value Ts of the static temperature of the structure surface. e The parameters b and C of the corresponding dynamic fatigue SN curve.

4. The thermo-mechanical fatigue equivalent load spectrum condensation method of claim 3, wherein, The characteristic parameters and conversion parameters of the single degree of freedom system comprise a system mass m, a system critical damping ratio ζ, a sound pressure action area A and a displacement-to-stress conversion factor K, wherein the system mass m takes a positive value within 20, the system critical damping ratio ζ takes a positive value within 0.03, the sound pressure action area A takes a positive value within 1, and the displacement-to-stress conversion factor K takes a positive value within 1000.

5. The thermo-mechanical fatigue equivalent load spectrum condensation method of claim 4, wherein, The process of calculating fatigue damage spectra caused by thermoacoustic loads under different working conditions based on fatigue damage spectrum theory of random vibration is as follows: 1) Extract the sound pressure power spectral density Φ under the working condition i p,i (f); 2) multiplying the sound pressure power spectral density Φ p,i (f) by the square of the sound pressure action area A to obtain the dynamic force power spectral density Φ F,i (f), i.e. Φ F,i (f) = A 2 • Φ p,i (f) ; 3) Set the natural frequency f of a single degree of freedom system j n,j , the power spectral density of the dynamic force Φ F,i (f) induced displacement response: ； where: H z,j (f) is the transfer function of the response displacement z to the dynamic force F for the single degree of freedom system j; m j is the mass of the single degree of freedom system j; f n,j is the natural frequency of the single degree of freedom system j; ζ j is the critical damping ratio of the single degree of freedom system j; 4) Calculate the number of mean positive crossings of the displacement response per unit of time and its root mean square value z rms,j : ； ； wherein: f c is the displacement response Φ z,ij (f) the corresponding cut-off frequency; 5) Since the response of a single degree of freedom system belongs to a narrow band process, the probability density function of its peak tends to a Rayleigh distribution, therefore the fatigue damage D j is expressed as: ; where: K is a conversion factor from displacement to stress; T i is the time of action of the thermoacoustic load; Γ is the gamma function; C i is the fatigue performance constant of the material under the i th working condition; b i is the negative reciprocal of the slope of the fatigue S-N curve of the material under the i th working condition 6) Select the natural frequency f of the next single degree of freedom system n,j+1 , corresponding to the sound pressure power spectral density Φ p,i (f) The next center frequency of the frequency column, repeat steps 3-5 above to get the fatigue damage D at this natural frequency j+1 ; 7) Repeat step 6 until the natural frequency is taken to the cutoff frequency f c Thus, a set of (f n -D) i Data, i.e. fatigue damage spectrum for the working condition i; 8) According to the calculation process of fatigue damage spectrum under the single working condition, the fatigue damage spectrum (f n -D) i+1 , (f n -D) i+2 , …, (f n -D) num , where num is the total number of working conditions.

6. The thermo-mechanical fatigue equivalent load spectrum condensation method of claim 5, wherein, The process of making fatigue damage spectrum caused by the external force power spectrum density of damage equivalent sound pressure equal to the cumulative damage spectrum is as follows: 1) Construct a white noise spectrum of dynamic pressure with a frequency range identical to that of the fatigue cumulative damage spectrum and an amplitude of 1 Pa 2 / Hz; 2) Set the natural frequency f of a single degree of freedom system n,j , calculate the fatigue damage D caused by white noise spectrum j , wherein the action time of thermoacoustic load in the fatigue damage calculation process caused by white noise spectrum is the total action time of all working conditions of thermoacoustic load, that is, the action time T of equivalent fatigue load spectrum e ; 3) Extracting the natural frequency f n,j The sum of the amplitude of the corresponding cumulative damage spectrum and the fatigue loss (∑D i ) j The sum of the fatigue loss (∑D i ) j The ratio Am j of the fatigue damage D j ​ ； wherein: b e is the negative inverse of the slope of the equivalent material fatigue S-N curve, i.e. the equivalent characteristic parameter b e ; 4) the ratio Am j as the natural frequency f n,j corresponding to the amplitude of the damage equivalent sound pressure external force power spectral density; 5) Repeat steps 2-4 by changing the natural frequency of the single degree of freedom system, and obtain the amplitude of the external force power spectral density of the equivalent sound pressure of the damage at different natural frequencies, i.e. f n -Am.

7. The method of thermomechanical fatigue equivalent load spectrum condensation according to any one of claims 1 to 6, characterized in that The method further comprises the following steps: checking effectiveness of the equivalent thermoacoustic fatigue equivalent load spectrum induction result, the process comprising: taking the equivalent fatigue load spectrum of the thermoacoustic load as a thermoacoustic load, calculating fatigue damage spectrum caused by the thermoacoustic load according to the calculation process of the fatigue damage spectrum, and checking coincidence of two curves by plotting the fatigue damage spectrum and the cumulative damage spectrum on a same graph; calculating a relative error of fatigue damage under the same natural frequency according to data of the two curves, and checking the effectiveness of the equivalent thermoacoustic fatigue equivalent load spectrum induction result according to the relative error.

8. A device for inducing a thermo-acoustic fatigue equivalent load spectrum, characterized in that, The method comprises the following steps: The load spectrum acquisition module is configured to acquire a thermal-acoustic fatigue load spectrum under different working conditions, and determine an equivalent value Ts of the structure surface static temperature according to the thermal-acoustic fatigue load spectrum e ; The SN curve acquisition module is configured to acquire dynamic fatigue S-N curves of the structural material at different temperatures, and determine an equivalent characteristic parameter b of the dynamic fatigue S-N curve according to an equivalent value of the static temperature of the structure surface e and C e ; a parameter determination module for determining characteristic parameters and conversion parameters of a single degree of freedom system; a damage spectrum calculation module for calculating fatigue damage spectra caused by thermoacoustic loads under different working conditions based on fatigue damage spectrum theory of random vibration, linearly superimposing fatigue damages under the same natural frequency of all working conditions to obtain a cumulative damage spectrum; an induction module for taking the cumulative damage spectrum as a target, constructing an external force power spectrum density of damage equivalent sound pressure, making fatigue damage spectrum caused by the external force power spectrum density of damage equivalent sound pressure equal to the cumulative damage spectrum, and thus obtaining an equivalent fatigue load spectrum of the thermoacoustic load.

9. The apparatus for induction of thermo-mechanical fatigue equivalent load spectrum of claim 8, wherein, The thermoacoustic fatigue load spectrum comprises an action time, a structure surface static temperature and a sound pressure power spectrum density, and an equivalent value of the structure surface static temperature is a structure surface static temperature corresponding to a maximum value of the action time in the thermoacoustic fatigue load spectrum.

10. The apparatus for induction of a thermo-mechanical fatigue equivalent load spectrum according to claim 9, characterized in that use Characterizing the SN curve, where S p Let N be the stress amplitude, N be the number of fatigue cycles, and b and C be characteristic parameters. The equivalent characteristic parameter b of the dynamic fatigue SN curve is... e and C e Select the equivalent value Ts of the static temperature of the structure surface. e The parameters b and C of the corresponding dynamic fatigue SN curve.

11. The apparatus for induction of a thermo-mechanical fatigue equivalent load spectrum according to claim 10, characterized in that The characteristic parameters and conversion parameters of the single degree of freedom system include system mass m, system critical damping ratio ζ, sound pressure acting area A and displacement to stress conversion factor K, wherein the system mass m takes a positive value within 20, the system critical damping ratio ζ takes a positive value within 0.03, the sound pressure acting area A takes a positive value within 1, and the displacement to stress conversion factor K takes a positive value within 1000.

12. The apparatus for induction of a thermo-mechanical fatigue equivalent load spectrum according to claim 11, characterized in that Based on the fatigue damage spectrum theory of random vibration, the process of calculating the fatigue damage spectrum caused by the thermoacoustic load under different working conditions is as follows: 1 ) the sound pressure power spectral density Φ in the extraction condition i p,i (f); 2) multiplying the sound pressure power spectral density Φ p,i (f) by the square of the sound pressure action area A to obtain the dynamic force power spectral density Φ F,i (f), i.e. : Φ F,i (f) = A 2 · Φ p,i (f) ; 3) Set the natural frequency f of a single degree of freedom system j n,j , the power spectral density of the dynamic force Φ F,i (f) induced displacement response: ； where: H z,j (f) is the transfer function of the response displacement z to the dynamic force F for the single degree of freedom system j; m j is the mass of the single degree of freedom system j; f n,j is the natural frequency of the single degree of freedom system j; ζ j is the critical damping ratio of the single degree of freedom system j; 4) the number of mean crossings of the displacement response per unit time and its root mean square value z rms,j : ； ； wherein: f c is the displacement response Φ z,ij (f) the corresponding cut-off frequency; 5) Since the response of a single degree of freedom system belongs to a narrow band process, the probability density function of its peak tends to a Rayleigh distribution, therefore the fatigue damage D j is expressed as: ; where: K is a conversion factor from displacement to stress; T i is the duration of the thermoacoustic load; Γ is the gamma function; C i is the fatigue performance constant of the material under the i th working condition; b i is the negative reciprocal of the slope of the fatigue S-N curve of the material under the i th working condition 6) Select the natural frequency f of the next single degree of freedom system n,j+1 , corresponding to the sound pressure power spectral density Φ p,i (f) The next center frequency of the frequency column, repeat steps 3-5 above to get the fatigue damage D at this natural frequency j+1 ; 7) Repeat step 6 until the natural frequency is taken to the cutoff frequency f c Thus, a set of (f n -D) i Data, i.e. fatigue damage spectrum for the working condition i; 8) According to the calculation process of the fatigue damage spectrum under the single working condition, the fatigue damage spectrum (f n -D) i+1 , (f n -D) i+2 , …, (f n -D) num , where num is the total number of working conditions.

13. The apparatus for induction of thermo-mechanical fatigue equivalent load spectrum of claim 12, wherein, By constructing the external force power spectrum density of the damage equivalent sound pressure, the fatigue damage spectrum caused by the external force power spectrum density of the damage equivalent sound pressure is equal to the cumulative damage spectrum, and the process is as follows: 1) Construct a white noise spectrum of dynamic pressure with a frequency range identical to that of the fatigue cumulative damage spectrum and an amplitude of 1 Pa 2 / Hz; 2) Set the natural frequency f of a single degree of freedom system n,j , calculate the fatigue damage D caused by white noise spectrum j , wherein the action time of thermoacoustic load in the fatigue damage calculation process caused by white noise spectrum is the total action time of all working conditions of thermoacoustic load, that is, the action time T of equivalent fatigue load spectrum e ; 3) Extracting the natural frequency f n,j The sum of the corresponding cumulative damage spectrum amplitude and fatigue loss (∑D i ) j , the sum of fatigue loss (∑D i ) j The ratio Am j of the fatigue damage D j ; ； wherein: b e is the negative inverse of the slope of the equivalent material fatigue S-N curve, i.e. the equivalent characteristic parameter b e ; 4) the ratio Am j as the natural frequency f n,j corresponding to the amplitude of the damage equivalent sound pressure external force power spectral density; 5) Repeat steps 2-4 by changing the natural frequency of the single degree of freedom system, and obtain the amplitude of the external force power spectral density of the equivalent sound pressure of the damage at different natural frequencies, i.e. f n -Am.

14. The apparatus for inducing a spectrum of equivalent thermal-mechanical fatigue loads according to any one of claims 8 to 13, wherein Also includes: The test module is used for testing the validity of the equivalent thermoacoustic fatigue equivalent load spectrum induction result, and the process includes: Taking the equivalent fatigue load spectrum of the thermoacoustic load as the thermoacoustic load, calculating the fatigue damage spectrum caused by the thermoacoustic load according to the calculation process of the fatigue damage spectrum caused by the thermoacoustic load, and comparing the fatigue damage spectrum with the cumulative damage spectrum on a graph to check the coincidence degree of the two curves; According to the data of the two curves, the relative error of the fatigue damage at the same natural frequency is calculated, and the validity of the equivalent thermoacoustic fatigue equivalent load spectrum induction result is tested according to the relative error.

15. An electronic device, comprising: Including: One or more processors; Memory; One or more application programs stored in the memory and configured to be executed by the one or more processors, the one or more application programs being configured for implementing the thermoacoustic fatigue equivalent load spectrum induction method according to any one of claims 1-7.

16. A computer-readable storage medium, characterized in that, The computer readable storage medium stores at least one instruction, at least one program, a code set or an instruction set, which is loaded and executed by the processor to implement the thermoacoustic fatigue equivalent load spectrum induction method according to any one of claims 1-7.

Images