Kurtosis transmission law of amplitude modulated non-Gaussian signals

Abstract

There are many methods published about how to create non-Gaussian signals with given Power Spectral Density (PSD) and kurtosis, to use for vibration testing and/or fatigue life evaluation. However, some signals show significantly different fatigue damage characteristics, expressed as Fatigue Damage Spectrum (FDS), compared to others, highlighting the need for a proper understanding of their properties. A specific non-Gaussian signal generation method is studied in depth in this paper, where a Gaussian signal with defined PSD is multiplied with a slowly varying Amplitude Modulation Signal (AMS). This AMS is created by raising the absolute value of a Gaussian signal to a power p. A new algebraic expression for calculation of AMS parameter from given kurtosis is then presented. To investigate how kurtosis transmits from the excitation to the system response, simulations were conducted on single degree of freedom (SDOF) systems with varying AMS and system properties. A new kurtosis transfer model is then proposed and validated through vibration experiments on a cantilever beam. The results show that with a given input kurtosis, the response kurtosis can be accurately predicted from the SDOF resonance frequency divided by AMS cutoff frequency and SDOF system Q factor, which lays the foundation for fatigue life evaluation of structures under non-Gaussian excitations.

Keywords

non-Gaussian vibration amplitude modulation kurtosis control kurtosis transmission law

Introduction

Random vibration tests are widely used to evaluate fatigue life of structures. Using the Power Spectral Density (PSD) and Inverse Fast Fourier Transform (IFFT), Gaussian signals are generated by modern vibration controllers to simulate random excitation. In practice, however, field measured data are usually non-Gaussian,^1,2 which may result in non-Gaussian response and a shorter fatigue life.³

In recent years, there is an increasing interest in generating non-Gaussian signals with given PSD and kurtosis.^4,5 The generated non-Gaussian signals can be broadly divided into two families: steady non-Gaussian with peaks and burst non-Gaussian. Commonly used methods for steady non-Gaussian signal generation include polynomial transformation,⁶ filtered Poisson process,⁷ phase manipulation,⁸ and combination of shocks and random signals.⁹ A comprehensive review can be found in.¹⁰ For burst non-Gaussian signal generation, the most popular method is the amplitude modulation method.¹¹ A stationary Gaussian signal is generated and then multiplied with a slowly varying Amplitude Modulation Signal (AMS) that is independent of the Gaussian signal. The crux of this method is the modeling of AMS. The kurtosis of a signal generated using this method is governed by the kurtosis of AMS. Different methods have been proposed to model the AMS, for example, Beta distribution,¹¹ Gamma distribution¹² and Weibull distribution.¹³ Among these methods, a wide range of kurtosis was covered by modeling the AMS with the absolute values of the standard normal distribution (called ‘absolute Gaussian method’ in this paper),¹⁴ and only one parameter of the AMS needs to be iterated to reach target kurtosis value. The problem with current amplitude modulation methods is that the parameters of the AMS have not been linked to the target kurtosis of the generated burst non-Gaussian signal such that an efficient control of kurtosis can be realized, that is, the mathematical relationship between target kurtosis and the parameters of AMS has not been treated in much details.

Ideally, an accelerated fatigue test can be conducted using the generated non-Gaussian signals as excitations. However, the kurtosis transmission from excitation to structural response remains a key issue that needs to be addressed. By decomposing the burst non-Gaussian excitation signal by 3.2 times the half-power bandwidth at the damped natural frequency of the structure, a linear model for the transmission of kurtosis values between the decomposed signal of the excitation and the structural response is established in Lei et al.¹⁵ Braccesi et al.¹⁶ used four different PSDs with increasing bandwidth to show that the bandwidth of the input burst non-Gaussian signal does not significantly affect the kurtosis of structural response. Using higher-order spectra (HOS), Trap¹⁷ showed that the non-Gaussianity due to non-stationarity is fully transferred to the response and is independent of the bandwidth of ideal bandpass filters. Kihm et al.¹⁸ used a ‘kurtosis rate’ to describe the input-output kurtosis relationship in the case of burst non-Gaussian excitation with the aid of numerical simulations. The ‘kurtosis rate’ was shown to be related to the time required for the Single Degree of Freedom (SDOF) system IRF to decrease to some fraction of its initial amplitude, and the period of the bursts. The kurtosis of the output signal varies between 3.0 and the kurtosis of the input signal. By weighting and synthesizing the frequency response characteristics and the PSD components of the excitation, Lei et al.¹⁹ established a kurtosis transmission prediction model under complex vibration conditions. Although an analytical expression was derived, the effects of key parameters such as the cutoff frequency of the non-Gaussian excitation and the system quality factor were not considered. Ren et al.²⁰ developed a modified time-domain randomization method for generating and controlling non-stationary non-Gaussian vibration signals in multi-shaker systems, incorporating amplitude modulation to regulate time-varying RMS characteristics. Although the transmission law of kurtosis has been studied by some researchers, problems remain unsolved. The explicit relationship between the kurtosis of the system response and the characteristics of system and AMS needs to be further developed.

To solve the aforementioned problems, the absolute Gaussian method is firstly studied in depth, giving algebraic expressions for calculation of AMS parameter from given kurtosis of burst non-Gaussian signal. A mathematical model is then proposed to quantify how kurtosis transmits from the excitation to the SDOF system response, incorporating the effects of AMS cut-off frequency and system Q factor. Experimental validation is conducted to verify the proposed model.

The remainder of this paper is organized as follows. In Section 2, the absolute Gaussian method is reviewed and a new algebraic expression for calculation of AMS parameter from given kurtosis is presented. In Section 3, the response of SDOF is computed using digital filters, and the effects of the cut-off frequency of AMS (U_cut) and Q factor on kurtosis transmission are analyzed. In Section 4, a new kurtosis transfer model is presented. In Section 5, a cantilever beam is used for experimental validation of the proposed model. Finally, the results and discussion are presented in Section 6.

Theoretical background

The n-th central moment M_n of a random signal x(t) is defined as $M_{n} (x) = E {(x - E {x})^{n}}$ , where E{} denotes expected value. Kurtosis is defined as a normalized fourth central moment:

K_{x} = \frac{M_{4} (x)}{{(M_{2} (x))}^{2}}

(1)

For a Gaussian signal, the kurtosis is equal to 3. For signals with kurtosis larger than 3, the probability for high amplitudes is larger than for a Gaussian signal with the same Root Mean Square (RMS) value. The Probability Density Function (PDF) of signals with different kurtosis but same RMS value of 6 m/s² is shown in Figure 1. From Figure 1 we can see that a higher kurtosis value leads to a wider tail. A wider tail means that the probability of large peak values is bigger, which leads to larger fatigue damage and quicker fail.

Figure 1.

PDF for signals with the same RMS but different kurtosis values.

To generate a non-Gaussian signal y(t), a Gaussian signal G(t) with a given PSD is firstly created and then multiplied by a slowly varying amplitude modulation signal u(t):

y (t) = G (t) \cdot u (t)

(2)

Here, u(t) is defined as the absolute value of a Gaussian signal g(t), which is generated to have a cutoff frequency U_cut, raised to a power p.¹⁵

u (t) = {| g (t) |}^{p}

(3)

Since the Gaussian signal G(t) has a mean value of zero, the non-Gaussian signal y(t) also has a mean value of zero. G(t) and u(t) are independent, the kurtosis K_y for y(t) is:

\begin{matrix} K_{y} = \frac{E {y^{4}}}{E {y^{2}} \cdot E {y^{2}}} = \frac{E {{(G \cdot u)}^{4}}}{E {{(G \cdot u)}^{2}} \cdot E {{(G \cdot u)}^{2}}} = \\ \frac{E {G^{4}} \cdot E {u^{4}}}{E {G^{2}} \cdot E {G^{2}} \cdot E {u^{2}} E {u^{2}}} = K_{G} \cdot K_{u} = 3 \cdot K_{u} \end{matrix}

(4)

It has to be noted that the mean of u(t) is not zero, so K_u is called ‘kurtosis with mean’.

In order not to change the PSD of G(t), u(t) is normalized to have a mean square value of one. Equation (4) can be simplified as:

K_{y} = 3 \cdot E {u^{4}}

(5)

To get a target kurtosis K_y, the only thing needed is the fourth moment of u(t), E{u⁴}.

Using the following equation²¹

\int_{0}^{\infty} x^{b} \cdot e^{- a x^{2}} dx = \frac{Γ (\frac{b + 1}{2})}{2 \sqrt{a^{b + 1}}}

(6)

A general formula for the m-th moment of u(t) as function of m and p is:

M_{m} (p) = π^{(m / 4 - 1 / 2)} \cdot \frac{Γ (m \cdot p / 2 + 1 / 2)}{{Γ (p + 1 / 2)}^{m / 2}}

(7)

where $Γ$ is the Gamma function.

When m = 4, equation (5) can be presented as:

\begin{matrix} K_{y} = 3 \cdot K_{u} = 3 \cdot E {u^{4}} = 3 \cdot \\ M_{4} (p) = 3 \cdot \sqrt{π} \frac{Γ (2 p + 1 / 2)}{{Γ (p + 1 / 2)}^{2}} \end{matrix}

(8)

A check of the validity of equation (8) is given in Figure 2. The fourth moment of u(t) is calculated as function of p using equation (8) and compared with a simulation of E{u⁴} based on equation (3). As we can see from Figure 2, a good agreement is achieved.

Figure 2.

The fourth moment of the AMS u(t).

From equation (8), given a target kurtosis K_y, the parameter p can be obtained by solving the following nonlinear equation:

f (p) = \frac{Γ (2 p + \frac{1}{2})}{{(Γ (p + \frac{1}{2}))}^{2}} - \frac{K_{y}}{3 \sqrt{π}} = 0

(9)

With the calculated p, the AMS can be created using equation (3) and the non-Gaussian signal can be generated using equation (2).

A case study is presented here to demonstrate the aforementioned process. The vibration test specification NAVMAT profile, which refers to the U.S. Navy standard NAVMAT P-9492 random vibration test specification, is used as the reference PSD,²² as shown in Figure 3.

Figure 3.

PSD for the NAVMAT profile used for simulations. RMS = 6.06 g.

The target kurtosis is set from 3.25 to 12 in steps of 0.25. The cut-off frequency U_cut for the AMS is set to 5 Hz. With a low U_cut, the synthesis of the low frequent Gaussian g(t) has a rather big random variation, so the function contains an iteration to fit the required tolerance. A good match is shown in Figure 4.

Figure 4.

Kurtosis for simulated time signal is shown as a function of given kurtosis.

Influence of U_cut and Q-factor on kurtosis transmission

Damage potential of non-Gaussian signals is usually evaluated by Fatigue Damage Spectrum (FDS), which is based on fatigue calculated on responses of a series of SDOF systems, using cycle count, Miner’s rule and Wohler exponent coefficient. The damage is related to the response kurtosis, and how input signal kurtosis is transmitted to the response kurtosis is called kurtosis transmission. To study kurtosis transmission, a non-Gaussian signal with NAVMAT PSD, kurtosis 8 and U_cut 1 Hz is loaded to a SDOF system shown in Figure 5.

Figure 5.

Linear SDOF system.

The pseudo velocity transfer function of the SDOF system can be expressed as²³:

H (s) = \frac{- 2 π f_{n}}{s^{2} + \frac{2 π f_{n} \cdot s}{Q} + {(2 π f_{n})}^{2}}

(10)

Where, $f_{n}$ is the SDOF resonance frequency, $s$ is the Laplace variable, Q-factor is the resonance gain.

The pseudo velocity response of the SDOF system is calculated by the ramp invariant digital filter method. The response kurtosis of a series of SDOF systems are calculated and shown in Figure 6. From Figure 6 we can see that the response kurtosis is close to the input kurtosis of 8 for all $f_{n}$ .

Figure 6.

Kurtosis for SDOF filter responses for non-Gaussian signal with kurtosis 8 and U_cut 1 Hz.

If the same test is done for kurtosis 8 but U_cut 40 Hz, the result is quite different, as shown in Figure 7.

Figure 7.

Kurtosis for SDOF filter responses for non-Gaussian signal with kurtosis 8 and U_cut 40 Hz.

From Figure 7 we can see that the response kurtosis changes from about 3 to 8 over the range of used $f_{n}$ . The kurtosis behavior in Figure 7 needs an explanation. It is well understood¹⁶ and more, that the meaning of ‘slow enough’ for u(t) is coupled to the length of the SDOF Impulse Response Function (IRF). The IRF with envelope is shown in Figure 8. As the length of the SDOF IRF is depending on $f_{n}$ /Q, the response kurtosis for different Q is calculated and shown in Figure 9. Each Q value is tested with a separate simulation. In the figure, a polynomial fit for each curve is added.

Figure 8.

IRF and envelope for SDOF 200 Hz, Q = 10, pseudo velocity.

Figure 9.

SDOF response kurtosis for different Q values.

As the length of SDOF IRF is depending on $f_{n} / Q$ , a new plot is made with that as x-axis and shown in Figure 10. The alignment of the kurtosis curves for different Q values confirms the assumption that response kurtosis is a function of SDOF IRF length. The result for simulations with different U_cut, keeping Q = 10, is shown in Figure 11. As for the test with different Q values, each U_cut value is tested with a separate simulation.

Figure 10.

SDOF response kurtosis for different Q values against $f_{n}$ /Q.

Figure 11.

SDOF response kurtosis for different U_cut values.

As 1/U_cut determines some correlation time length for the signal, a new plot is made with $f_{n}$ /Q/U_cut as x-axis and shown in Figure 12. The result confirms that the kurtosis may be written as a function of $f_{n}$ /Q/U_cut.

Figure 12.

SDOF response kurtosis for different U_cut values against $f_{n} / Q / U_{cut}$ .

A new kurtosis transmission law

A simulation is made with 64 combinations of input kurtosis $K_{inp}$ , Q-values, and U_cut values ( $K_{inp}$ 5, 6, 7, 8, U_cut3, 7, 11, 15, and Q 10, 15, 20, 25) and compared with the calculated response kurtosis from equation (13). The result is shown in Figure 13.

Figure 13.

SDOF response kurtosis for 64 combinations simulations.

By assuming that there is a function Φ independent of input kurtosis, an algebraic expression for the response kurtosis is given by the equation

\begin{matrix} K_{res} = 3 + Φ (K_{inp} - 3) \\ Φ = \frac{K_{res} - 3}{K_{inp} - 3} \end{matrix}

(11)

That is validated in Figure 14, where the result from all 64 simulations is recalculated to function Φ.

Figure 14.

Function Φ.

A simple least curve fit to the function Φ in Figure 14 is made:

Φ = \frac{1.4899 x + 0.5955 x^{2}}{1 + 1.3705 x + 0.6008 x^{2}}

(12)

Where,

x = f_{n} / Q / U_{cut}

(13)

The fit is added to Figure 14 giving Figure 15 and to Figure 13 to give Figure 16.

Figure 15.

Function Φ and added curve fit.

Figure 16.

Comparison of SDOF response kurtosis between 64 simulations and new kurtosis transmission law.

With the fit given by equation (11), the estimated SDOF response kurtosis may be found as a function of input kurtosis and U_cut, defining the non-Gaussian signal, $f_{n}$ and Q, defining the SDOF system. Many simulations with different PSD and kurtosis specifications seem to show that the general kurtosis shape as in Figure 16 and equation (11) is not influenced by the specific PSD shape and/or kurtosis.

To verify whether AMS statistics influence kurtosis shape, we simulated four signals with kurtosis 9 for comparison with a commercial signal (Ho 1 in Figure 17) showing measured kurtosis 9.06, whose generation method remains unknown, but there is a ‘transition frequency’ parameter set to 10 Hz. Three distinct AMS signals were used in simulation: a simple low-pass Gaussian signal multiplied by a NAVMAT Gaussian signal (Gauss*Gauss in Figure 17) yielding simulated kurtosis 9.03. The second AMS had Weibull distribution (Weibull in Figure 17), total signal kurtosis 9.03 and the third signal was synthesized from equation (9) (psdkurtsyntag in Figure 17), kurtosis 9.04. All those three had a nominal U_cut set to 10 Hz. Figure 17 displays the PSDs of these four signals.

Figure 17.

PSD for four different high kurtosis simulations and nominal low frequency cut-off for AMS.

SDOF systems response kurtosis for the four simulations is calculated as before. The plot of kurtosis, Figure 18 is against $f_{n}$ /Q/U_cut. The four curves align well, the commercial signal Ho 1 was rather short, so the random variation is bigger for that one. The shapes of PSD for the AMS differ between the methods, so a small adjustment of the cut-off frequency is done for the plot. Some equivalent band width should be used.

Figure 18.

SDOF response kurtosis for four different simulations.

The result of the simulation is that the statistics of the AMS signal has small, if any, effect on the kurtosis curve shape.

It should be noted that although the application scope of the proposed kurtosis transfer model is restricted to linear SDOF systems, one particular useful application is the easy calculation of FDS, which makes the fatigue life evaluation of structures much simpler.

Experimental verification

To validate the proposed kurtosis transmission model, experimental verification was performed on a cantilever beam. The test setup is shown in Figure 19. The cantilever beam was excited by a vibration shaker. An accelerometer (PCB 356A15) at the base served as the input end and was connected to channel 0, while an accelerometer (PCB 352C33) at the free end served as the response end and was connected to channel 2. The excitation signal was monitored via a T-connector and recorded on channel 1. All signals were acquired using a National Instruments 9234 data acquisition device, and a host computer was employed for data storage. The sampling frequency is 2048 Hz, and the signal length is 2²² samples.

Figure 19.

Experimental test setup.

The experimentally measured natural frequency is 23.125 Hz, as shown in Figure 20.

Figure 20.

Measured response PSD.

Due to the nature of the vibration shaker system, the actual excitation reaching the beam becomes distorted, as shown in Figure 21. To ensure that the excitation PSD applied to the cantilever beam is flat around its resonance frequency, the input PSD was compensated for the frequency-dependent transfer characteristics of the shaker–fixture–beam system. This compensation was achieved by measuring the transfer function between the input signal and the beam response and iteratively adjusting the input PSD until the desired flat response around resonance was obtained. The PSD modification is shown in Figure 22.

Figure 21.

Original input PSD and excitation to the beam.

Figure 22.

Modified input PSD and excitation to the beam.

Theoretically, to validate the kurtosis transmission law, multiple cantilever beams with different natural frequencies are needed. In this paper, however, by adjusting the cutoff frequency of the non-Gaussian excitation, only one beam is needed. The equivalent $f_{n}$ is:

f_{equ} = \frac{f_{e} * U_{s}}{U_{e}}

(14)

where:

$f_{equ}$ is the equivalent $f_{n}$ , $U_{s}$ is the cutoff frequency used in simulations (15 Hz), $f_{e}$ is the natural frequency of the cantilever beam (23.125 Hz), $U_{e}$ is the cut off frequency of excitation signal.

The response kurtosis is calculated against (equivalent $f_{n}$ )/Q/U_cut:

X_{axis} = \frac{f_{equ}}{U_{s} * Q} = \frac{f_{e}}{U_{e} * Q}

(15)

Experiments were conducted for eight different $U_{e}$ values. The corresponding equivalent $f_{n}$ and $X_{axis}$ values are summarized in the table below, a Q value of 40 is used (Table 1).

Table 1.

Equivalent $fn$ and X-axis.

Cut off frequency of excitation	Equivalent $f_{n}$	X-axis
0.0578	6001.29	10.00
0.2	1734.37	2.890
0.431	804.81	1.341
0.929	373.38	0.622
2.0	173.43	0.289
4.31	80.48	0.134
9.29	37.33	0.0622
20	17.34	0.0289

By setting the kurtosis of excitation signals to 6, 7 and 8, a total of 24 measurements were made. The results are shown in Figure 23.

Figure 23.

A comparison between kurtosis from many simulations and kurtosis from experiments.

Figure 23 demonstrates excellent agreement between simulations and experimental results, confirming the validity of the proposed new kurtosis transmission law.

To quantitatively evaluate the accuracy of the proposed kurtosis transfer model, the mean absolute error (MAE) and the root mean square error (RMSE) between the simulated kurtosis values and the experimental measurements was calculated:

MAE = \frac{1}{N} \sum_{i = 1}^{N} | K_{m} - K_{s} |

(16)

RMSE = \sqrt{\frac{1}{N} \sum_{i = 1}^{N} (K_{m} - K_{s})^{2}}

(17)

where $K_{m}$ denotes the experimentally measured kurtosis, $K_{s}$ denotes the simulated kurtosis, as listed in Table 2, N is the number of comparison points, with N = 8.

Table 2.

MAE and RMSE.

Input kurtosis	6		7		8
x-Axis	k _m	k _s	k _m	k _s	k _m	k _s
0.0289	3.1559	3.1268	3.2140	3.1646	3.2617	3.4205
0.0622	3.2646	3.2097	3.3690	3.2794	3.4464	3.5100
0.134	3.5176	3.4660	3.6693	3.6146	3.8914	3.8266
0.289	3.9780	3.9826	4.2935	4.2237	4.6535	4.5955
0.622	4.6229	4.5003	5.2103	5.1115	5.7709	5.6650
1.341	5.3583	5.3697	6.1350	6.1688	6.9244	6.6791
2.890	5.7644	5.8661	6.7772	6.5987	7.7266	7.5179
10	5.9737	5.9850	6.9594	7.0911	7.9559	9.0151
MAE	0.0484		0.0920		0.1205
RMSE	0.0634		0.1016		0.1393

These small errors indicate a good match between the predicted and measured kurtosis, showing the high accuracy and reliability of the proposed kurtosis transfer model under non-Gaussian excitation.

It should be noted that the experimental results are obtained based on a single cantilever beam specimen. The validation of the proposed conclusions using multiple cantilever beams with different physical natural frequencies remains an important topic for future investigation.

Results and discussion

In this work, a specific non-Gaussian signal generation method is studied, where a Gaussian signal with defined PSD is multiplied with a slowly varying AMS. This AMS is created by raising the absolute value of a Gaussian signal to a power p. A new algebraic expression for calculation of AMS parameter from given kurtosis is then presented, which is capable of efficiently and accurately generating specified non-Gaussian signals, enabling precise and efficient kurtosis control of specified non-Gaussian signals.

Since non-stationary non-Gaussian signals induce bigger fatigue damage, accurate prediction of structural response kurtosis under non-Gaussian random excitation is particularly crucial for fatigue life calculation. Based on the response characteristics of SDOF systems, this study investigated the effects of the AMS cutoff frequency and Q-factor on kurtosis transmission, thereby establishing a new kurtosis transmission law that links the response kurtosis to both excitation and structural parameters. A novel kurtosis transmission model was proposed, which quantitatively relates the response kurtosis to the input kurtosis, AMS cutoff frequency and system Q-factor. Experimental verification on a cantilever beam confirmed the model’s effectiveness, showing that predicted response kurtosis agrees well with measured results.

In conclusion, the proposed approach provides a clear and practical framework for predicting structural response kurtosis under non-Gaussian excitation, which lays a foundation for vibration fatigue life evaluation under non-Gaussian random excitation.

Footnotes

Acknowledgements

This study was supported by The National Natural Science Foundation of China (Grant # 52102443).

Handling Editor: Aarthy Esakkiappan

ORCID iDs

Fei Xu

Wenjia Wang

Ethical considerations

Not applicable.

Author contributions

Writing—original draft: Fei Xu, Kjell Ahlin, Wenjia Wang, supervision, funding acquisition: Fei Xu, Methodology: Fei Xu, Kjell Ahlin, Simulation: Fei Xu, Wenjia Wang, Validation: Fei Xu, Kjell Ahlin, Writing—review and editing: Fei Xu, Kjell Ahlin, Wenjia Wang. All authors reviewed the results and approved the final version of the manuscript.

Funding

The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was funded by The National Natural Science Foundation of China (Grant # 52102443).

Declaration of conflicting interests

The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Data availability statement

The data used to support the findings of this study are available from the corresponding author upon request.

References

Hui

Kawai

, et al. Non-stationary and non-Gaussian characteristics of wind speeds. Wind Struct 2017; 24(1): 059–078.

Decker

Savaidis

. Measurement and analysis of wheel loads for design and fatigue evaluation of vehicle chassis components. Fatigue Fract Eng Mater Struct 2002; 25(12): 1103–1119.

Steinwolf

Cornelis

Peeters

, et al. On the use of kurtosis control methods in shaker testing for fatigue damage. J Test Eval 2020; 48(1): 538–556.

Smallwood

. Generating non-Gaussian vibration for testing purposes. Sound Vib 2005; 39(10): 18–23.

Steinwolf

. Vibration testing of vehicle components by random excitations with increased kurtosis. Int J Veh Noise Vib 2015; 11(1): 39–66.

Steinwolf

Peeters

Van der Auweraer

(eds). Two methods of generating random excitations with increased kurtosis for in-house testing of vehicle components. In: Proceedings of the 22nd international congress on sound and vibration, 2015. The International Institute of Acoustics and Vibration.

Jiang

. Generation of non-Gaussian random vibration excitation signal for reliability enhancement test. Chin J Aeronaut 2007; 20: 236–239.

Steinwolf

. Random vibration testing with kurtosis control by IFFT phase manipulation. Mech Syst Signal Process 2012; 28: 561–573.

Zheng

Chen

Wang

, et al. Kurtosis and crest factor simultaneous control for non-Gaussian random vibration test. Mech Syst Signal Process 2025; 222: 111760.

10.

Zheng

Chen

. Stationary non-Gaussian random vibration control: a review. Chin J Aeronaut 2021; 34(1): 350–363.

11.

Smallwood

. Vibration with non-Gaussian noise. J IEST 2009; 52(2): 13–30.

12.

Cui

Geng

Wang

, et al. Study of the amplitude modulation method for kurtosis control purposes. Mech Syst Signal Process 2022; 179: 109399.

13.

Jiang

, et al. Synthesis of running RMS-induced non-Gaussian random vibration based on Weibull distribution. J Vibroengineering 2015; 17(7): 3662–3674.

14.

Trapp

Makua

Wolfsteiner

. Fatigue assessment of amplitude-modulated non-stationary random vibration loading. Procedia Struct Integr 2019; 17: 379–386.

15.

Lei

Jiang

Zeng

, et al. Research on the transmission law of kurtosis of SDOF system under nonstationary and non-Gaussian random excitations. Mech Syst Signal Process 2022; 165: 108292.

16.

Braccesi

Cianetti

Palmieri

, et al. The importance of dynamic behaviour of vibrating systems on the response in case of non-Gaussian random excitations. Procedia Struct Integr 2018; 12: 224–238.

17.

Trapp

Hollweck

Wolfsteiner

. On the transmission of non-Gaussian random loading through linear structures. Procedia Struct Integr 2022; 37: 622–631.

18.

Kihm

Rizzi

Ferguson

, et al. Understanding how kurtosis is transferred from input acceleration to stress response and it’s influence on fatigue life. Proc., 11th International Conference on Recent Advances in Structural Dynamics (RASD 2013), Pisa, Italy, 2013.

19.

Lei

Jiang

Tian

, et al. Kurtosis transfer in MDOF systems under amplitude-modulated non-stationary non-Gaussian excitations: a frequency-domain decomposition approach. Mech Syst Signal Process 2025; 236: 113059.

20.

Ren

Chen

Zheng

. Modified time domain randomization technique for multi-shaker non-stationary non-Gaussian random vibration control. Mech Syst Signal Process 2024; 213: 111311.

21.

Tuma

Walsh

. Engineering mathematics handbook. McGraw-Hill, 1979, p.335.

22.

U.S. Navy Department. NAVMAT P-9492. Navy Manufacturing Screening Program, 1979.

23.

ISO 18431-4:2007. Mechanical vibration and shock - signal processing - part 4: shock-response spectrum analysis, ISO 18431-4:2007(E). International Organization for Standardization, ISO, 2007.

Kurtosis transmission law of amplitude modulated non-Gaussian signals

Abstract

Keywords

Introduction

Theoretical background

Influence of Ucut and Q-factor on kurtosis transmission

A new kurtosis transmission law

Experimental verification

Results and discussion

Footnotes

Acknowledgements

ORCID iDs

Ethical considerations

Author contributions

Funding

Declaration of conflicting interests

Data availability statement

References

Influence of U_cut and Q-factor on kurtosis transmission