Published: 15 November 2015

Synthesis of running RMS-induced non-Gaussian random vibration based on Weibull distribution

Fei Xu1
Chuanri Li2
Tongmin Jiang3
Kjell Ahlin4
1Blekinge Tekniska Högskola, 371 79 Karlskrona, Sweden
1, 2, 3School of Reliability and System Engineering, Beihang University, 37 Xue Yuan Road, Hai Dian District, Beijing 100191, China
4Xielalin Consulting, S. Skogrundan 38 Akersberga, SE-18463, Sweden
Corresponding Author:
Fei Xu
Views 213
Reads 95
Downloads 1303

Abstract

Gaussian signal is produced by ordinary random vibration controllers to test the products in the laboratory, while the field data usually is non-Gaussian. To synthesize non-Gaussian random vibration, both the probability density function (PDF) and the damage effects must be considered. A new method is presented in this paper to synthesize non-Gaussian random vibration that is characterized by running RMS (root mean square). The essential idea is to model the non-Gaussian signal by a Gaussian signal multiplied by an amplitude modulation function (AMF). A two-parameter Weibull distribution is used to model the PDF of the running RMS and to create the AMF. The shock response spectrum (SRS) is used to detect significant shocks within the non-Gaussian signal. A case study is presented to show that the synthesized non-Gaussian signal has the same power spectral density (PSD), kurtosis, PDF and fatigue damage spectrum (FDS) as the field data.

1. Introduction

Random vibration testing is usually used to bring a test item to failure to identify weaknesses in the product or to verify if the product can survive a particular random vibration environment. Random vibration controllers accomplish this goal by producing a power spectral density (PSD) that would expose the test item to the type of vibratory environment that the test item would experience in a real-world setting. Since PSD does not contain phase information, vibration controllers always assume the phase of each frequency component follows a 0 to 2π uniform distribution when producing time history using the PSD. This leads to data with Gaussian distribution. Non-Gaussian vibration is usually encountered, especially in the road transportation. W. H. Connon and D. Charles points out that when vehicles are running on irregular roads, the random excitation placed on vehicles is often non-Gaussian. The tail of measured acceleration PDF is wider and the middle is sharper when compared with Gaussian distribution [1, 2]. This difference may lead to totally different accumulated fatigue damage, because most fatigue damage is caused by 2 to 4σ [3]. When peak value exceeds 3σ, products fail quickly. MIL-STD810 points out that we should always check if the field tested data is non-Gaussian and the testing hardware and software is appropriate [4]. The time waveform replication (TWR) is frequently referred to as a methodology for non-Gaussian testing. The basic idea of TWR is to reproduce a sequence of instantaneous values of the vibration process. Such a test may be non-Gaussian, however this is only a replication of one particular measured record, not a simulation of a specified road type. Besides, MIL-STD810F says that TWR is usually used for controlling transient signals or short time lasting random vibration [4]. With these drawbacks of TWR, many researchers presented some other methods to simulate non-Gaussian vibration. Smallwood used three zero-memory nonlinear (ZMNL) functions to transform Gaussian vibration to non-Gaussian vibration [5, 6]. These functions all give similar, but not identical results. Each covers similar but slightly different ranges of skewness and kurtosis. Each will result in slightly different distributions. Initially, this method was recommended in the form where non-Gaussian time histories are first prepared offline before running a test and, then, these time histories are reproduced on the shaker in the Time Waveform Replication mode. It was only later when this method was extended and modified into a method that can be used in the closed-loop frequency domain control mode [8, 9]. This method is easy to implement but it significantly reduces the controller’s dynamic range compared to the traditional Gaussian random testing and may disrupt the PSD shape. Stein wolf developed a special phase selection method [7-9]. The essence of this method is that the amplitudes are still responsible for the PSD, whereas the phases are not all random as they were in the classic IFFT generation. Instead, some of the phases are now used to adjust kurtosis. Thus, both characteristics, the PSD and the kurtosis, are controlled independently of each other. Rouillard et al. presented a novel technique by which non-Gaussian vibrations are synthesized by generating a sequence of random Gaussian processes of varying RMS levels and durations [10-15]. John et al. studied the relationship between kurtosis and fatigue damage spectrum (FDS) [16]. He declares that the fatigue damage spectrum shows the increasing fatigue damage across all frequencies that accompany the increased kurtosis levels that cause the faster product failure.

To synthesize the non-Gaussian random vibration, both the PDF and the damage effects must be considered. A new method is presented in this paper to synthesize the non-Gaussian random vibration that is characterized by the running RMS (root mean square). The essential idea is to model the non-Gaussian signal by a Gaussian signal multiplied by a slowly varying amplitude modulation function (AMF). A two-parameter Weibull distribution is used to model the PDF of the running RMS and to create the AMF. The shock response spectrum (SRS) is used to detect significant shocks within the non-Gaussian field signal. No sequence of random Gaussian processes of varying RMS levels and durations needs to be generated in this method, which makes it much easier to be implemented than that Rouillard et al. used [10-15]. A case study is presented to show that the synthesized non-Gaussian signal has the same PSD, kurtosis, PDF and FDS as the field data.

2. Theory

2.1. Non-Gaussian random signal

It is well known that the PSD is sufficient for describing the characteristics of stationary random signals only if the process under consideration is Gaussian. To check if a random process is Gaussian, the PDF of the process must be compared with the Gaussian model:

1
p(x)=12πσxexp-x-μx22σx2,

where μx is the mean value, σx is the standard deviation.

Two principal characters describing non-Gaussian PDF features are skewness and kurtosis:

2
S=m3σx3, K=m4σx4,

where mn is the central moments of the PDF:

3
mn=-x-μxnpxdx.

Eq. (3) is the theoretical definitions of the central moments. For a given digitized data record, it can be calculated by time averaging:

4
mn=1Nj=1NxjΔt-μxn, μx=1Nj=1NxjΔt.

For a truly Gaussian process, the skewness is 0 and the kurtosis is 3. Any deviation from these indicates that the process is non-Gaussian. From the perspective of simulating non-Gaussian vibration, kurtosis is a more important parameter than skewness, because it represents the probability of peak values in time history. The PDF of data produced by vibration testing can be more realistic using this parameter. The kurtosis describes the sharpness or flatness in the middle range of PDF and the narrowness or wideness in the tail of PDF. Kurtosis larger than 3 indicating the middle range of PDF of field tested signal is sharper than a Gaussian signal with the same standard deviation, and the tail is wider. Wider tail means that the probability of peak values is larger, which leads to larger fatigue damage and quicker fail.

2.2. Fatigue damage analysis

The Rain Flow Cycle Counting (RFCC) and Dirlik methods are frequently referred to as two methodologies for fatigue damage analysis. The RFCC method was initially proposed by M. Matsuiski [17] and T. Endo [18] to count the cycles or the half-cycles of strain-time signals. For fatigue damage analysis under vibration, the cyclic numbers nRFCCi at ith stress level (i= 1, 2,..., p) is counted directly from the stress time history. The number of stress levels p depends on how many bins are used to perform the RFCC. See the counting algorithm in [19] for detail.

T. Dirlik [20] established empirical expressions of the PDF of the ordinary half-ranges and those counted with the RFCC method using a digital simulation. Let x(t) be the response of a SDOF system to a Gaussian signal, xrms the rms value of x(t), δ a range defined starting from RFCC and u the reduced variable corresponding to a half-cycle:

5
u=δ2xrms.

The PDF of u is expressed as:

6
p(u)=D1Ae-uA+D2uB2e-u22B2+D3ue-u22,

where:

7
D1+D2+D3=1, D1=2C-r21+r2, D2=1-r-D1+D121-B,
8
A=1.25r-D3-D2BD1, B=r-C-D121-r-D1+D12, C=M1M0M2M4,
9
r=M2M0M4, Mn=fnGxxfdf.

Using the given PSD, the mean number of positive maxima per second can be expressed as [21]:

10
np+=M4M2.

With Eqs. (6) and (10), the cyclic numbers nDirliki at ith stress level can be expressed as:

11
nDirliki=Tnp+pudu u=ui=TM4M2pudu u=ui, i=1, 2,...,p,

where T is the time duration of the vibration environment that the product experienced, in seconds.

At each stress level, a comparion between nDirliki and nRFCCi for both Gaussian and non-Gaussian vibrations will be made later, to show the limitation of the Dirlik method when dealing with non-Gaussian signal.

2.3. Fatigue damage spectrum and shock response spectrum

All structures experience fatigue as they are repeatedly exposed to adequate levels of stress. The total damage a product experiences in a particular time period can be calculated from field data and plotted for a specific range of frequencies. The FDS is essentially a plot that shows the responses of a number of single degree of freedom (SDOF) systems to a base input acceleration time history. Many SDOF systems tuned to a range of natural frequencies are assessed using the same input. The final plot, the FDS, shows the fatigue damage encountered for a particular SDOF system anywhere within the analyzed time, as shown in Fig. 1 and Fig. 2.

Fig. 1SDOF systems for calculation of FDS

SDOF systems for calculation of FDS

Fig. 2Calculation procedure of the FDS

Calculation procedure of the FDS

For a SDOF system with a natural frequency fn and a damping ratio ξ, the output relative displacement xd to an input acceleration xa can be calculated using a ramp invariant digital filter method [22]. With the output xd, the cumulative damage can be calculated in both time domain and frequency domain. In time domain, RFCC method is usually used to count cyclic numbers at different stress levels. Then S-N curve [23] and Miner’s rule [24] are combined to calculate fatigue damage:

12
Ni=cSi-b,
13
Dt=i=1pniNi=i=1pniSibc,

where ni is the number of cycles exposure at Si (i=1, 2,..., p). p is the number of stress levels considered. Ni is fatigue life at stress Si, b is the fatigue exponent, Dt is the total damage index calculated in time domain.

The stress is assumed to be proportional to relative displacement [25]:

14
S=kxd.

The total damage can be calculated as:

15
Dt=i=1pniNi=i=1pnikbxd,ibc=kbci=1pnixd,ib.

In frequency domain, Rayleigh distribution of response stress maxima is assumed and used to calculate the stress cycles [26]:

16
pS=SσS2e-S22σS2,

where S is the stress value of peaks, σS is the RMS of the stress time history.

The total damage can be calculated as:

17
Df=np+T0+pSNSdS=np+Tc0+pSSbdS=fnTc0+pSSbdS,

where np+ can be replaced by the natural frequency fn for a lightly damped SDOF system response, T is the total time of exposure to the stress environment.

Solving Eq. (17) with Eq. (16) leads to:

18
Df=fnTc2σ S2b2Γ1+b2=fnTckb2σ d2b2Γ1+b2,

where Γ is the Gamma function, σd is the RMS of relative displacement.

The RMS of relative displacement can be calculated using:

19
σd=f=flowf=fhiHf2Gxxfdf,

where Hf is the transmissibility of a SDOF system with a natural frequency of fn and a damping ratio of ξ.

Since different kurtosis leads to different number of peak accelerations in the test, and these peak accelerations are primarily responsible for damage that a product experiences in the field, the FDS should show the difference when calculated for Gaussian and non-Gaussian vibrations.

The effects of the shock loads are generally depicted in a SRS. Like the FDS, the SRS is essentially a plot that shows the responses of a number of SDOF systems to a base input acceleration time history, as shown in Fig. 1. The only difference is that the SRS is generated by calculating the maximum response of a SDOF system to the input. The final plot, the SRS, shows the largest response encountered for a particular SDOF system anywhere within the analyzed time.

Since the essential idea of the new synthesis method is to model the non-Gaussian signal by a Gaussian signal multiplied by an AMF:

20
xnonGaussian=xGaussianfAMF,

significant shocks in the field data must be first identified and separated. To detect if there are significant shocks in the field data, the response of a number of SDOF systems to the field data is firstly calculated in time domain. Then the response is divided into M consecutive segments and the maximum value is calculated for each segment. Finally the ratio between this maximum value and the SRS computed from the PSD, which is calculated from the whole time history using Welch method, is used to show where the significant shock is:

21
Rn,i=20lgyn,iSRSn, n=1, 2,...,N, i=1, 2,...,M,

where yn,i indicates the maximum value calculated for the ith segment and the nth SDOF system, SRSn indicates the SRS calculated form PSD for the nth SDOF system.

3. Synthesis method

A new method is presented to synthesize the non-Gaussian random vibration characterized by the running RMS. The procedure uses three phases to streamline the process:

Phase 1: Check the characteristics of measured vibration signal.

1. Obtain the environment acceleration time history, resample if needed.

2. Calculate the kurtosis of measured signal to see if it is Non-Gaussian or not.

3. Calculate the PSD of measured signal using Welch’s method.

4. Use SRS to detect significant shocks (see Eq. (21)).

Phase 2: synthesize the Non-Gaussian random vibration.

1. Create the stationary Gaussian signal using the PSD calculated in step 3, Phase 1.

2. Calculate the probability distribution function FRMS of the running RMS of the field data.

3. Use a two parameter Weibull distribution to model the probability distribution of the running RMS.

The cumulative distribution function of Weibull distribution can be expressed as:

22
FWeibullx=1-exp-xLK,

where x is the random variable, L is the scale parameter, K is the shape parameter.

To find the correct L and K, the probability distribution function of the assumed Weibull distribution FWeibull and the running RMS FRMS are used to form an objective function using least squares method:

23
Fobj=stdFRMS-FWeibull=stdFRMS-1-exp-xLiKi,

where x are the values of the running RMS, depending on how many bins are used to calculated the probability distribution function of the running RMS, i indicates the ith iteration, std indicates standard deviation. Note that both FRMS and FWeibull are vectors.

The MATLAB function “fminsearch” is used find the right L and K, so that the standard deviation of error between these two probability distribution functions is minimized.

1. Use the determined L and K to create the AMF and to get the required Non-Gaussian signal.

2. To create the appropriate AMF, the kurtosis of the field data is used as an objective value. The AMF is created from the Weibull distribution function and multiplied with the Gaussian signal created in step 1, Phase 2, to get the synthesized signal. Compare the kurtosis of the synthesized signal with the objective value. Do iterations until the objective value is reached. The synthesized non-Gaussian signal is scaled to have the same RMS as the field data.

Phase 3: compare the characteristics of the synthesized Non-Gaussian signal and the field data.

1. Compare the PSD.

2. Compare the PDF.

3. Compare the fatigue, e.g., cyclic numbers and FDS.

4. Case study

In this case study, a new method for synthesizing the non-Gaussian random vibration is developed and examined by the field data. The field data provided by Prof. Kjell Ahlin is from an Ericsson Mast Project and is shown in Fig. 3. The test item and setup are shown in Fig. 4.

Fig. 3Field data

Field data

Fig. 4Test item and setup

Test item and setup
Test item and setup

This is a non-Gaussian signal with kurtosis equals to 9.4262. Eq. (21) is used to detect significant shocks (see Fig. 5). The ratio is shown in the third dimension by the color. As we can see from Fig. 5, there are only some insignificant shocks superimposed on the underlying random vibration signal.

The PSD of the field data is calculated using Welch method with 4096 points hanging window and 50 % overlapping. Gaussian signal is generated using the same but smoothed PSD and is shown in Fig. 6.

Fig. 5Detection of shocks

Detection of shocks

Fig. 6Synthesized Gaussian signal using the same PSD calculated from field data

Synthesized Gaussian signal using the same PSD calculated from field data

The PSD and PDF of both the field data and the synthesized Gaussian signal are shown in Fig. 7 and Fig. 8. As we can see, the PDF of the field data is obviously different from that of the synthesized Gaussian signal, although the PSD are basically the same.

Fig. 7PSD of field data and synthesized Gaussian signal

PSD of field data and synthesized  Gaussian signal

Fig. 8PDF of the field data and the synthesized Gaussian signal

PDF of the field data and the synthesized Gaussian signal

Fig. 9The running RMS

The running RMS

Fig. 10PDF of the running RMS using different number of bins and Rayleigh distribution

PDF of the running RMS using different number of bins and Rayleigh distribution

Fig. 11The iteration to determine L and K

The iteration to determine L and K

The running RMS and its PDF using different number of bins are calculated and shown in Fig. 9 and Fig. 10. The running RMS is calculated using 10 s frames with 50 % overlapping. From Fig. 10 we can see that, despite of the numbers of bins, the PDF of the running RMS is similar to the PDF of a given Rayleigh distribution. To be more general, a two-parameter Weibull distribution is used to model the probability distribution function of the running RMS Eq. (22).

Iteration method described in Section 3 (step 3, phase 2) is used to determine L and K. The iteration process is shown in Fig. 11. The probability distribution function of resulted Weibull distribution is shown in Fig. 12, together with the probability distribution function of the running RMS.

As the parameters of the Weibull distribution are determined, the AMF can be created and the non-Gaussian signal can be synthesized using the iteration method presented in Section 3 (step 4, phase 2). The synthesized non-Gaussian signal is shown in Fig. 13, together with the field data (the time duration of the synthesized non-Gaussian signal is adjustable). Note that some care must be taken to have the time scale for the AMF selected properly, since inappropriate selection of time scale may lead to inaccurate PDF and/or PSD.

The PSD and PDF of the field data and the synthesized non-Gaussian signal are very similar (see Fig. 14 and 15), which proves the effectiveness of the new method.

For comparison, the cycle counting with the RFCC and Dirlik method can be performed directly on the field data (acceleration signal), the synthesized Gaussian signal and non-Gaussian signal. The results are shown in Fig. 16. As we can see, the Dirlik method is equivalent to the RFCC for Gaussian signal. Big error can be introduced by the Dirlik method for non-Gaussian signal. The RFCC results for the field data and synthesized non-Gaussian signal are very similar.

Fig. 12The probability distribution function of the running RMS and Weibull distribution

The probability distribution function of the running RMS and Weibull distribution

Fig. 13Field data and synthesized non-Gaussian signal

Field data and synthesized non-Gaussian signal

Fig. 14Mast vibration signal power spectral density. The PSD of the field data and the synthesized non-Gaussian signal

Mast vibration signal power spectral density. The PSD of the field data  and the synthesized non-Gaussian signal

Fig. 17 shows the FDS calculated for the field data, the synthesized Gaussian and non-Gaussian signal in the time domain using Eq. (15). For comparison, k and c are normalized to 1. From Fig. 17 we can see that the FDS of the field data and the synthesized non-Gaussian signal are very similar. This again proves the effectiveness of the new method. Besides, we can also see that larger kurtosis leads to larger fatigue damage.

Fig. 15Amplitude probability density. PDF for field data and synthesized Gaussian and non-Gaussian signal

Amplitude probability density. PDF for field data and synthesized Gaussian  and non-Gaussian signal

Fig. 16Number of cycles versus data ranges for RFCC and Dirlik

Number of cycles versus data ranges for RFCC and Dirlik

Fig. 17Fatigue damage spectrum. The FDS of the field data, synthesized Gaussian and non-Gaussian signal

Fatigue damage spectrum. The FDS of the field data, synthesized Gaussian  and non-Gaussian signal

5. Conclusion

A new method for synthesizing the running RMS-induced non-Gaussian Random Vibration is developed and examined by field data. Conclusions are as below:

1. This method is limited to the case of the running RMS-induced non-Gaussian signals, which contain no significant shocks.

2. The synthesized non-Gaussian signal has the same PSD, kurtosis, PDF and FDS with the field data, which proves the effectiveness of the new method.

3. The Dirlik methods are equivalent to the RFCC method for Gaussian signal and can introduce big error for non-Gaussian signal.

4. Larger kurtosis leads to larger fatigue damage.

5. Weibull distribution can be used to simulate the statistical distribution of the running RMS and to synthesize the running RMS-induced non-Gaussian signal.

References

  • Connon W. H. Comments on kurtosis of military vehicle vibration data. Journal of the IES, Vol. 34, Issue 6, 1991, p. 38-41.
  • Charles D. Derivation of environment descriptions and test severities from measured road transportation data. Journal of the IES, Vol. 36, Issue 1, 1993, p. 37-42.
  • Lambert R. G. Fatigue life prediction for various random stress peak distributions. Shock and Vibration Inform. Center The Shock and Vibration Bulletin, Vol. 52, Issue 4, 1982, p. 1-10.
  • MIL-STD-810F. Department of Defense Test Method Standard for Environmental Engineering Considerations and Laboratory Tests, USA, 2000.
  • Smallwood D. O. Generation of stationary non-Gaussian time histories with a specified cross-spectral density. Shock and Vibration, Vol. 4, Issue 5, 1997, p. 361-377.
  • Smallwood D. O. Generating non-Gaussian vibration for testing purposes. Sound and Vibration, Vol. 39, Issue 10, 2005, p. 18-23.
  • Steinwolf A. Shaker simulation of random vibration with a high kurtosis value. Journal of the Institute of Environmental Sciences, Vol. 40, Issue 3, 1997, p. 33-43.
  • Steinwolf A. Closed-loop shaker simulation of non-Gaussian random vibrations. Part 1. Discussion and methods. Test Engineering and Management, Vol. 68, Issue 3, 2006, p. 10-13. Part 2. Numerical and experimental results. Test Engineering and Management, Vol. 68, Issue 5, 2006, p. 14-19.
  • Steinwolf A. Random vibration testing beyond PSD limitations. Sound and Vibration, Vol. 40, Issue 9, 2006, p. 12-21.
  • Rouillard V. On the Synthesis of Non-Gaussian Road Vehicle Vibrations. Monash University, 2007.
  • Rouillard V., Lamb M. On the effects of sampling parameters when surveying distribution vibrations. Packaging Technology and Science, Vol. 431, Issue 8, 2008, p. 467-477.
  • Rouillard V., Sek M. A. Synthesizing non-stationary, non-Gaussian random vibrations. Packaging Technology and Science, Vol. 23, Issue 8, 2010, p. 423-439.
  • Garcia-Romeu M. A., Rouillard V. On the statistical distribution of road vehicle vibrations. Journal of Packaging Technology and Science, Vol. 24, Issue 8, 2011, p. 451-467.
  • Rouillard V., Sek M. Creating transport vibration simulation profiles from vehicle and road characteristics. Packaging Technology and Science, Vol. 26, Issue 2, 2013, p. 82-95.
  • Rouillard V. Quantifying the non-stationarity of vehicle vibrations with the run test. Packaging Technology and Science, Vol. 27, Issue 3, 2014, p. 203-219.
  • Van Baren J., Van Baren P., Jenison M. I. The fatigue damage spectrum and kurtosis control. Sound and Vibration, Vol. 46, Issue 10, 2012, p. 10-13.
  • Matsuishi M., Endo T. Fatigue of metals subjected to varying stress. Japan Society of Mechanical Engineers, Fukuoka, Japan, 1968, p. 37-40.
  • Endo T. Damage evaluation of metals for random on varying loading-three aspects of rain flow method. Proceedings of Symposium on MBM, Vol. 1, 1974, p. 374.
  • Lalanne C. Mechanical Vibration and Shock, Fatigue Damage. Wiley, 2010, p. 111-120.
  • Lalanne C. Mechanical Vibration and Shock, Fatigue Damage. Wiley, 2010, p. 172-174.
  • Lalanne C. Mechanical Vibration and Shock, Random Vibration. Wiley, 2010, p. 281.
  • Ahlin K. Comparison of test specifications and measured field data. Sound and Vibration, Vol. 40, Issue 9, 2006, p. 22-25.
  • Lalanne C. Mechanical Vibration and Shock, Fatigue Damage. Wiley, 2010, p. 10-13.
  • Lalanne C. Mechanical Vibration and Shock, Fatigue Damage. Wiley, 2010, p. 49-50.
  • Lalanne C. Mechanical Vibration and Shock, Specification Development. Wiley, 2010, p. 76.
  • Papoulis A. Narrow-band systems and Gaussianity. IEEE Transactions on Information Theory, Vol. 18, Issue 1, 1972, p. 20-27.

About this article

Received
15 July 2015
Accepted
06 October 2015
Published
15 November 2015
SUBJECTS
Vibration generation and control
Keywords
non-Gaussian
random vibration
running RMS
Weibull distribution
fatigue damage spectrum (FDS)
shock response spectrum (SRS)