A validation of the stress power spectral density method with zero-order moment

Zhuang, N. N.; Chen, H. B.; Wang, Y. Y.

doi:10.21595/vp.2020.21444

Vibroengineering Procedia

Browse Procedia

Published: 07 May 2020

Check for updates

A validation of the stress power spectral density method with zero-order moment

N. N. Zhuang¹

H. B. Chen²

Y. Y. Wang³

^{1, 2, 3}Department of Modern mechanics, University of Science and Technology of China, CAS Key Laboratory of Mechanical Behavior and Design of Materials, Hefei, China

³Chengdu Aircraft Design and Research Institute, Chengdu, China

Corresponding Author:

H. B. Chen

Cite the article Download PDF

Downloads 1099

Abstract

To give a fatigue life prediction for problems with bimodal stress power spectral density which is commonly observed in the dynamic response of structures, the bimodal triangular model is established first and a large sum of bimodal stress power spectral densities with different characteristics are built. Then the error analysis of the fatigue life prediction results of bimodal stress power spectral densities with different shapes are calculated by the stress power spectral density method with zero-order moment (SPSD⁰) and the classical frequency-domain fatigue life prediction method, which is under single frequency band process and broadband process. The simulation results show that, for single frequency band process or broadband process being divided into appropriate segments, the fatigue damage is not sensitive to the real curve shape of bimodal stress power spectral density.

1. Introduction

Fatigue damage to a structure has always been one of the most important damage factors in engineering [1-3]. According to statistics, 50 %-90 % of the damage of mechanical components is fatigue damage [4]. During the service period of engineering structures, random loads are very common. Therefore, the fatigue life prediction and reliability prediction of the structure under random load is particularly important. The bimodal stress power spectral density is a very typical form of the structural dynamic response. It is characterized by two separate vibration modes [5], for example, the response of aircraft’s shell and plate subsystem under the excitation of complex pulsating noise. Its main shape of the stress power spectral density corresponds to the first two natural frequencies of the structure is shown as a double peak curve.

In recent years, through a large number of numerical simulations, Wang [6] found that fatigue analysis using the stress root mean square value of each frequency band usually has high engineering accuracy, based on the proposed SPSD⁰. However, he did not analyze the representative bimodal stress power spectral density. In this work, we selects the triangular model for the bimodal stress power spectral density. The fatigue life was predicted using the SPSD⁰, and compared with the fatigue life prediction results of Dirlik method. Thus, the validity of SPSD⁰ was further verified.

2. Fatigue life prediction method

2.1. Classical frequency-domain fatigue life prediction method (Dirlik method)

The frequency-domain fatigue life prediction method uses the structure’s stress power spectral density function to describe the stress amplitude or peak information in the frequency domain to obtain the probability density function. According to the S-N fatigue life curve of the material and the linear cumulative damage criterion, the fatigue life of a structure can be predicted. The frequency-domain fatigue life prediction method is simple in idea and has high calculation efficiency and is currently the most widely used fatigue life prediction method in engineering.

Dirlik [7] proposed a frequency domain fatigue life prediction method based on a large number of numerical simulation results, which is the most widely used in engineering. This method presents the probability density function of stress amplitude as a sum of an exponential distribution and two Rayleigh distributions [7]:

1

p_{D i r} (S) = \frac{1}{\sqrt{m_{0}}} [\frac{D_{1}}{Q} e^{- \frac{Z}{Q}} + \frac{D_{2} Z}{R^{2}} e^{- \frac{Z^{2}}{2 R^{2}}} + D_{3} Z e^{- \frac{Z^{2}}{2}}],

where:

Z = \frac{S}{\sqrt{m_{0}}}, D_{1} = \frac{2 (x_{m} - α_{2}^{2})}{1 + α_{2}^{2}}, D_{2} = \frac{1 - α_{2}^{2} - D_{1} + D_{1}^{2}}{1 - R}, D_{3} = 1 - D_{1} - D_{2},

Q = \frac{1.25 (α_{2} - D_{3} - D_{2} R)}{D_{1}}, R = \frac{α_{2} - x_{m} - D_{1}^{2}}{1 - α_{2} - D_{1} + D_{1}^{2}},

x_{m} = α_{1} ∙ α_{2} = \frac{m_{1}}{\sqrt{m_{0} ∙ m_{2}}} ∙ \frac{m_{2}}{\sqrt{m_{0} ∙ m_{4}}} .

The fatigue damage is:

2

{E [D]}_{D i r} = E [P] \frac{T}{C} \int_{0}^{\infty} S^{k} p_{D i r} (S) d S,

where $E [P] = \frac{1}{2 π} \sqrt{\frac{m_{4}}{m_{2}}}$ , $m_{i}$ is the $i$ -th moment of the stress power spectral density function $G (f)$ , $m_{i} = \int f^{i} G (f) d f$ .

2.2. The stress power spectral density method with zero-order moment (SPSD⁰) [6]

SPSD⁰divides the total stress power spectral density into $N$ small frequency bands. The upper and lower frequency limits of the $i$ -th small frequency band are $f_{U i}$ (Upper frequency) and $f_{L i}$ (Lower frequency) respectively, and the $n$ -th ( $n =$ 0, 1, 2, 4) moment of the stress power spectral density of the $i$ -th small frequency band is ${(m_{n})}_{i}$ .

When the simplified stress power spectral density is used for calculation, the stress power spectral density in the small frequency band is a fixed value. If the value of the stress power spectral density in the $i$ -th small frequency band is recorded as $h_{i}$ , the following relationship is satisfied: ${(m_{0})}_{i} = \int_{f_{L i}}^{f_{U i}} h_{i} d f = h_{i} (f_{U i} - f_{L i})$ .

Therefore, $h_{i} = \frac{{(m_{0})}_{i}}{f_{U i} - f_{L i}}$ .

In turn, each moment of the entire stress power spectral density can be obtained.

Zero-order moment of stress power spectral density: $m_{0} = \sum_{i = 1}^{N} {(m_{0})}_{i}$ .

First moment of stress power spectral density:

m_{1} = \sum_{i = 1}^{N} \int_{f_{L i}}^{f_{U i}} (h_{i} ∙ f) d f = \sum_{i = 1}^{N} \frac{f_{U i} + f_{L i}}{2} {(m_{0})}_{i} .

Second moment of stress power spectral density:

m_{2} = \sum_{i = 1}^{N} \frac{f_{U i}^{2} + f_{U i} f_{L i} + f_{L i}^{2}}{3} {(m_{0})}_{i} .

Fourth moment of stress power spectral density:

m_{4} = \sum_{i = 1}^{N} \frac{f_{U i}^{4} + f_{U i}^{3} {f_{L i} + f_{U i}^{2} f_{L i}^{2} + f}_{U i} f_{L i}^{3} + f_{L i}^{4}}{5} {(m_{0})}_{i} .

From the fatigue formula of the Dirlik method mentioned above, we know that the fatigue life can be predicted after zero, first, second and fourth order moments of the stress power spectral density are obtained. Here, SPSD⁰ combined with Dirlik method is called the zero-order moment Dirlik method, and the fatigue damage obtained is recorded as ${E [D]}_{D i r 0}$ .

Compared with the classical frequency-domain fatigue life prediction method, SPSD⁰ has its unique advantages: it is only necessary to know the zero-order moment of the stress power spectral density in each frequency band, in other words, without knowing the specific stress power spectral density the fatigue life prediction can be performed.

3. Bimodal stress power spectral density

The bimodal triangular stress power spectral density model (Fig. 1) can be formed by superimposing two triangular single peaks. On the whole, the overall shape depends on the specific characteristics of the two single peaks, specifically, the area of the two triangular spectra ( $A_{1}$ or $A_{2}$ ), center frequency ( $f_{1}$ or $f_{2}$ ), area ratio $R$ (Fig. 1(a)), frequency ratio $b$ (Fig. 1(b)) and bandwidth ratio $r$ (Fig. 1(c)) and Fig. 1(d)): $R = A_{2} / A_{1}$ , $b = f_{2} / f_{1}$ , $r = ∆ f / f_{c}$ .

Fig. 1The triangular model of bimodal stress power spectral density

Fig. 2The bimodal stress power spectral density

In the formula, $A_{1}$ and $A_{2}$ are the areas of the first and second peaks, $f_{1}$ and $f_{2}$ are the center frequencies of the first and second peaks, respectively. By changing the base and height of the peak, single peaks with the same area but different shapes can be obtained.

In this paper, the calculation is mainly carried out by numerical simulation on the premise of ensuring that the zero-order moment of the stress power spectral density in the frequency band is a fixed value. The effect of the shape of bimodal stress power spectral density on the amount of fatigue damage is examined. For the stress power spectral density in the frequency band $f_{0} ~ f_{3}$ , the shape of bimodal stress power spectral density is shown in Fig. 2, and the corresponding stress power spectral density function equation is as follows:

3

G (f) = \{\begin{array}{l} 0, & f_{0} \leq f \leq f_{1} - d_{1}, \\ \frac{h_{1}}{d_{1}} f + \frac{{d_{1} - f}_{1}}{d_{1}} h_{1}, & f_{1} - d_{1} \leq f \leq f_{1}, \\ - \frac{h_{1}}{d_{1}} f + \frac{{d_{1} + f}_{1}}{d_{1}} h_{1}, & f_{1} \leq f \leq f_{1} + d_{1}, \\ 0, & f_{1} + d_{1} \leq f \leq f_{2} - d_{2}, \\ \frac{h_{2}}{d_{2}} f + \frac{{d_{2} - f}_{2}}{d_{2}} h_{2}, & f_{2} - d_{2} \leq f \leq f_{2}, \\ - \frac{h_{2}}{d_{2}} f + \frac{{d_{2} + f}_{2}}{d_{2}} h_{2}, & f_{2} \leq f \leq f_{2} + d_{2}, \\ 0, & f_{2} + d_{2} \leq f \leq f_{3}, \end{array}

where:

d_{1} = \frac{∆ f_{1}}{2} = \frac{r_{1} f_{1}}{2}, d_{2} = \frac{∆ f_{2}}{2} = \frac{r_{2} f_{2}}{2} = \frac{r_{2} {b f}_{1}}{2},

and ( $r_{1}$ and $r_{2}$ are the bandwidth ratios of the two peaks):

h_{1} = \frac{2 A_{1}}{r_{1} f_{1}}, h_{2} = \frac{2 A_{2}}{r_{2} f_{2}} = \frac{2 {R A}_{1}}{r_{2} {b f}_{1}} .

4. Fatigue life prediction analysis

4.1. Sensitivity analysis of fatigue damage in single frequency band process

In the frequency band $f_{0} ~ f_{3}$ , the zero-order moment of the stress power spectral density is a fixed value, that is $A_{1} + A_{2} = (1 + R) A_{1} = m_{0}$ , and $A_{1} = \frac{1}{1 + R} m_{0}$ can be obtained. Combining with the equation of the double-peak stress power spectral density in Eq. (3), we can see that the specific shape of the stress power spectral density can be determined by the parameters $f_{1}$ , $b$ , $R$ , $r_{1}$ and $r_{2}$ .

Considering the influence of randomness to construct the bimodal stress power spectral density, when the parameters take different values, the specific stress power spectral density constructed also varies. The relative variables are chosen as the following rules: the frequency band $f_{0} ~ f_{3}$ is divided into $N$ equal parts with the same frequency resolution as $1 / 200 f_{0}$ , and the position of $f_{1}$ locates at each $N$ equal point. Once $f_{1}$ determined, divide the frequency band $f_{1} ~ f_{3}$ into six equal parts, and the position of $f_{2}$ locates at each six equal point, then, the frequency ratio $b$ picks 5 values in the interval $({1, f}_{3} / f_{1})$ . In order to ensure that each single peak is a narrow-band spectrum, the bandwidth ratios $r_{1}$ and $r_{2}$ take the values of 0.1, 0.2, 0.3, 0.4 and 0.5 respectively; the area ratio $R$ is defined as 0.01, 0.05, 0.25, 0.5, 1, 2, 5, 10, respectively. By traversing these parameters, a large number of bimodal stress power spectral densities with different characteristics can be obtained, and the error analysis of fatigue damage can be carried out through numerical simulation for different fatigue life prediction methods.

We define the relative error of fatigue damage to represent the difference between the fatigue damage calculated by the constructed stress power spectral density and the simplified stress power spectral density (the stress power spectral density in each subsection is constant). The relative error of fatigue damage is:

4

ε_{D i r 0} = |\frac{{E [D]}_{D i r} - {E [D]}_{D i r 0}}{{E [D]}_{D i r 0}}| \times 100 %,

where ${E [D]}_{D i r}$ is the fatigue damage calculated by the Dirlik method, and ${E [D]}_{D i r 0}$ is the fatigue damage calculated by the zero-order moment Dirlik method. For different stress power spectral densities with the same $m_{0}$ value, their ${E [D]}_{D i r 0}$ value is the same, but the ${E [D]}_{D i r}$ value is different, so here the ${E [D]}_{D i r 0}$ value is used as a reference value.

Finally, take out the maximum relative error of the fatigue damage of all stress power spectral densities constructed from multiple cycles:

5

ε_{D i r 0}^{m} = M a x (ε_{D i r 0}),

where $ε_{D i r 0}^{m}$ represents the maximum relative error of the fatigue damage calculated by the Dirlik method between the bimodal stress power spectral density and the simplified stress power spectral density that are respectively constructed when the value of $m_{0}$ is fixed.

The variation curve of the maximum relative error of the fatigue damage calculated by the zero-order moment Dirlik method and the Dirlik method with the frequency ratio $f_{3} / f_{0}$ is shown in Fig. 3.

The abscissa represents the frequency ratio of the frequency band, and the relationship satisfies ${f_{3} / f_{0} = 10}^{1 / i}$ . For different frequency ratios, the variables $f_{1}$ , $b$ , $R$ , $r_{1}$ and $r_{2}$ are traversed to determine the number of simulations. The calculation results of ${E [D]}_{D i r 0}$ and ${E [D]}_{D i r}$ are related to the value of the S-N curve material parameter $k$ . For metal materials commonly used in engineering, the value of k is generally 2 to 6.

It can be seen from Fig. 3 that as $i$ increases (i.e., the bandwidth of the frequency band decreases gradually), the maximum relative error decreases gradually, and the change in $k$ value has little effect on the result. When $i =$ 10 (the frequency ratio is about 1.26), the maximum relative error of ${E [D]}_{D i r 0}$ is about 9 %. Such a result is an acceptable error in fatigue life prediction. Therefore, we can conclude that for a single frequency-band process with a frequency ratio of 1.26, if its $m_{0}$ is kept unchanged, the fatigue damage is not sensitive to the specific shape of the bimodal stress power spectral density.

Fig. 3The maximum relative error of E[D]Dir0 for a single frequency-band

Fig. 4The maximum relative error of E[D]Dir0 for a broad frequency-band

4.2. Sensitivity analysis of fatigue damage in broadband process

Take the broadband bimodal stress power spectral density of $f_{3} / f_{0} =$ 10 as an example for calculation, divide the broadband process into $n$ segments, the frequency ratio of each segment is $10^{1 / n}$ . Set the zero-order moment value of each small section of the simplified stress power spectral density and the bimodal stress power spectral density in the part of the frequency band to be equal, then the zero-order moment of the entire broadband process remains unchanged. The classical frequency-domain fatigue life prediction process is the same as that predicted by the single frequency band, while the prediction process of SPSD⁰ becomes a little more complicated than in the single frequency band. The maximum relative error of fatigue damage calculated by the zero-order moment Dirlik method is shown in Fig. 4.

It can be seen from Fig. 4 that for the broadband process with a frequency ratio of 10, keep its $m_{0}$ unchanged, when it is divided into 10 sub-sections, the frequency ratio of each small frequency band is about 1.26. The maximum relative error of ${E [D]}_{D i r 0}$ is about 10 % at this time, and such result is still an acceptable error in fatigue life prediction. It shows that when the broadband process is divided into enough number of segments, the fatigue damage is not sensitive to the specific shape of the bimodal stress power spectral density.

5. Conclusions

In this work, the common bimodal stress power spectral density in engineering is used to calculate its fatigue life using the SPSD⁰and Dirlik method. It is proved that for narrow-band processes (such as the ratio of the upper frequency limit to the lower frequency limit not exceeding 1.26) or broadband processes that are divided into appropriate segments (such as the ratio of the upper frequency limit to the lower frequency limit of each sub-section not exceeding 1.26), it is only necessary to know the zero-order moment of the stress power spectral density to predict the fatigue life.

SPSD⁰ has a simple calculation process and little information is required for fatigue life prediction, so it has strong application potential in the future. In this work, the bimodal stress power spectral density is used to systematically examine its calculation accuracy, which proves that SPSD⁰ has a high calculation accuracy. Especially for the broadband process, when the entire frequency band is divided into a large number of frequency bands (such as 10), the zero-order moment Dirlk method can even achieve the accuracy equivalent to the Dirlk method.

References

Miles J. W. On structural fatigue under random loading. Journal of the Acoustical Society of America, Vol. 29, Issue 1, 1957, p. 176-176.

Publisher
Carmona E., Paneque M., Poveda M. L., et al. Damage and failure of composite materials. AIAA Journal, Vol. 52, Issue 5, 2012, p. 317-323.

Search CrossRef
Zhang H. W., He Y. T., Fan C. H., et al. Fatigue life prediction method for aircraft metal material under alternative corrosion /fatigue process. Acta Aeronautica et Astronautica Sinica, Vol. 34, Issue 5, 2013, p. 1114-1121, (in Chinese).

Search CrossRef
Zhao S. B., Wang B. Z. Anti-Fatigue Design. Mechanical Industry Press, Beijing, 1997, (in Chinese).

Search CrossRef
Jiao G., Moan T. Probabilistic analysis of fatigue due to Gaussian load processes. Probabilistic Engineering Mechanics, Vol. 5, Issue 90, 1990, p. 76-83.

Publisher
Wang Y. Y. Research on High-Frequency Structural-Acoustic Coupling Analysis and Fatigue Life Prediction Methods. University of Science and Technology of China, 2016, (in Chinese).

Search CrossRef
Dirlik T. Application of Computers in Fatigue Analysis. University of Warwick, 1985.

Search CrossRef

About this article

Received

30 April 2020

Accepted

05 May 2020

Published

07 May 2020

SUBJECTS

Fault diagnosis based on vibration signal analysis

DOI

https://doi.org/10.21595/vp.2020.21444

Keywords

bimodal stress power spectral density

fatigue life prediction

the stress power spectral density method with zero-order moment

frequency domain fatigue life prediction method

Acknowledgements

This work was financially supported by the National Natural Science Foundation of China (NSFC) under Grant No. 11772322, and the Strategic Priority Research Program of the Chinese. Academy of Sciences under Grant No. XDB22040502.

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Previous article in issue Previous Next article in issue Next

Research article

2017 05 30

Fatigue life estimation algorithm of structures with high frequency random vibration based on hybrid energy finite element method

J. Zeng, H. B. Chen, Y. Y. Wang

Research article

2016 12 08

Fatigue analysis of simplified automobile model based on SEA and hybrid FE-SEA method

K. J. Li, H. B. Chen, Y. Y. Wang, C. Wang

Research article

2016 08 15

Time-frequency techniques for modal parameters identification of civil structures from acquired dynamic signals

Carlos Andres Perez-Ramirez, Juan Pablo Amezquita-Sanchez, Hojjat Adeli, Martin Valtierra-Rodriguez, Rene de Jesus Romero-Troncoso, Aurelio Dominguez-Gonzalez, Roque Alfredo Osornio-Rios

Research article

2015 11 15

A vibration fatigue analysis model considering interaction effects

Shan Jiang, Wei Zhang, Zili Wang

N. N. Zhuang, H. B. Chen, and Y. Y. Wang, “A validation of the stress power spectral density method with zero-order moment,” Vibroengineering PROCEDIA, Vol. 31, pp. 27–32, May 2020, https://doi.org/10.21595/vp.2020.21444

Copy Extrica

Copied to clipboard!

TY  - JOUR
DO  - 10.21595/vp.2020.21444
UR  - https://doi.org/10.21595/vp.2020.21444
TI  - A validation of the stress power spectral density method with zero-order moment
T2  - Vibroengineering PROCEDIA
AU  - Chen, H. B.
AU  - Zhuang, N. N.
AU  - Wang, Y. Y.
PY  - 2020
DA  - 2020/05/07
PB  - JVE International Ltd.
SP  - 27-32
VL  - 31
SN  - 2345-0533
SN  - 2538-8479
ER  - 

Copy Ris

Copied to clipboard!

@article{Chen_2020,
	doi = {10.21595/vp.2020.21444},
	url = {https://doi.org/10.21595/vp.2020.21444},
	year = 2020,
	month = {may},
	publisher = {{JVE} International Ltd.},
	volume = {31},
	pages = {27--32},
	author = {H. B. Chen and N. N. Zhuang and Y. Y. Wang},
	title = {A validation of the stress power spectral density method with zero-order moment},
	journal = {Vibroengineering {PROCEDIA}}
}

Copy Bibtex

Copied to clipboard!

[1]H. B. Chen, N. N. Zhuang, and Y. Y. Wang, “A validation of the stress power spectral density method with zero-order moment,” Vibroengineering PROCEDIA, vol. 31, pp. 27–32, May 2020, doi: 10.21595/vp.2020.21444.

Copy IEEE

Copied to clipboard!

Chen, H. B., N. N. Zhuang, and Y. Y. Wang. “A Validation of the Stress Power Spectral Density Method with Zero-Order Moment.” Vibroengineering PROCEDIA 31 (May 7, 2020): 27–32. https://doi.org/10.21595/vp.2020.21444.

Copy Chicago

Copied to clipboard!

A validation of the stress power spectral density method with zero-order moment

Abstract

1. Introduction

2. Fatigue life prediction method

2.1. Classical frequency-domain fatigue life prediction method (Dirlik method)

2.2. The stress power spectral density method with zero-order moment (SPSD0) [6]

3. Bimodal stress power spectral density

4. Fatigue life prediction analysis

4.1. Sensitivity analysis of fatigue damage in single frequency band process

4.2. Sensitivity analysis of fatigue damage in broadband process

5. Conclusions

References

About this article

Related Articles

2.2. The stress power spectral density method with zero-order moment (SPSD⁰) [6]