# On the Shaker Simulation of Wind-Induced Non-Gaussian Random Vibration

Authors: Fei Xu; Chuanri Li; Tongmin Jiang
Journal: Shock and Vibration
 (2016)

---

## Abstract

Gaussian signal is produced by ordinary random vibration controllers to test the products in the laboratory, while the field data is usually non-Gaussian. Two methodologies are presented in this paper for shaker simulation of wind-induced non-Gaussian vibration. The first methodology synthesizes the non-Gaussian signal offline and replicates it on the shaker in the Time Waveform Replication (TWR) mode. A new synthesis method is used to model the non-Gaussian signal as a Gaussian signal multiplied by an amplitude modulation function (AMF). A case study is presented to show that the synthesized non-Gaussian signal has the same power spectral density (PSD), probability density function (PDF), and loading cycle distribution (LCD) as the field data. The second methodology derives a damage equivalent Gaussian signal from the non-Gaussian signal based on the fatigue damage spectrum (FDS) and the extreme response spectrum (ERS) and reproduces it on the shaker in the closed-loop frequency domain control mode. The PSD level and the duration time of the derived Gaussian signal can be manipulated for accelerated testing purpose. A case study is presented to show that the derived PSD matches the damage potential of the non-Gaussian environment for both fatigue and peak response.


---

## Body

## 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. Historically, 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. However, a random process can be characterized by a PSD only when it is stationary and Gaussian distributed. In practice, non-Gaussian vibration is usually encountered, especially in the road transportation [1, 2]. The difference of the probability density function (PDF) between the measured acceleration and Gaussian distribution may lead to totally different accumulated fatigue damage [3]. 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 [4]. 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. Smallwood presented a zero memory nonlinear function (ZMNL) method to synthesize a given non-Gaussian data [5]. Analytical solution was presented later to extend the original offline method to a method that can be used in the closed-loop frequency domain control mode. Iteration method was presented to address the PSD distortion problem introduced by the nonlinear transformation [6]. A. Steinwolf presented a phase selection method to simulate the road vehicle measured data [7]. In particular, selected phase is transformed from random to deterministic in order to obtain a prescribed kurtosis. PSD and PDF are controlled independently. An analytical relation between kurtosis, amplitude, and phase at specific frequencies was presented later to make this method applicable in a closed-loop control implemented by the shake controller [8–10]. These methods mentioned above are applicable when the non-Gaussianity is caused by significant shocks in vehicle vibration data. For wind-induced vibration data, however, the non-Gaussianity may be caused by other factors, for example, by the running root mean square. Rouillard and Sek presented a novel technique to synthesize non-Gaussian vibrations by generating a sequence of random Gaussian processes of varying RMS levels and durations [11].Two methodologies are presented in this paper for shaker simulation of wind-induced vibration. The essential idea of the first methodology is to synthesize the non-Gaussian signals offline and replicate them on the shaker in the Time Waveform Replication mode. The non-Gaussian signal is modeled as a Gaussian signal multiplied by an amplitude modulation function (AMF). A two-parameter Weibull distribution is used to create the AMF. A case study is presented to show that the synthesized non-Gaussian signal has the same power spectral density (PSD), probability density function (PDF), and loading cycle distribution (LCD) as the field data. The rainflow cycle counting (RFCC) method [12] and Dirlik method [13] are used for the LCD calculation.The fatigue damage spectrum (FDS) and the extreme response spectrum (ERS) have been used as a quantitative way to take the fatigue damage and overstress damage into consideration for field data [14–20]. The FDS and ERS are used in the second methodology to derive a damage equivalent Gaussian signal (characterized by a PSD and time duration) from the non-Gaussian signal, so that they can be easily produced on the shaker in the closed-loop frequency domain control mode. Improvement of the derived PSD by taking its spectrum shape into consideration is also presented. The PSD level and the duration time of the derived Gaussian signal can be manipulated for accelerated testing purpose. A case study is presented to show that the derived PSD matches the damage potential of non-Gaussian environment for both fatigue and peak response.


## 2. Theory

### 2.1. Non-Gaussian Random Signal

If a Gaussian process has a zero mean value and is ergodic, then the PDF of the instantaneous acceleration valuesx
(
t
) is given by the Gaussian distribution with zero mean: (1)
p
x
=
1
2
π
σ
x
exp
⁡
-
x
2
2
σ
x
2
,where σ
x is standard deviation.For large durationsT, the variance σ
x
2 is given by(2)
σ
x
2
≈
1
T
∫
0
T
x
2
d
t
=
∫
0
∞
G
x
x
f
d
f
,where G
x
x
f is the single-sided PSD.It shows that random Gaussian processes with zero mean can be completely described by the PSD function. One useful method for establishing how well a random process can be described by the Gaussian distribution is by computing the higher order moments of the process defined as (continuous and discrete form)(3)
m
n
=
∫
-
∞
∞
x
n
p
x
d
x
=
1
N
∑
j
=
1
N
x
j
n
.When the mean value is zero,(4)
σ
x
2
=
1
N
∑
j
=
1
N
x
j
2
=
m
2
,
S
=
1
N
∑
j
=
1
N
x
j
3
σ
x
3
=
m
3
σ
x
3
,
K
=
1
N
∑
j
=
1
N
x
j
4
σ
x
4
=
m
4
σ
x
4
,where S is the skewness and K is the kurtosis.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.


### 2.2. Response Spectra

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 that are 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 Figures1 and 2.Figure 1
SDOF systems for calculation of FDS.Figure 2
Calculation procedure of the FDS.It has been shown that the stress is roughly proportional to pseudovelocity [21]. Thus, it is more reasonable to calculate the pseudovelocity FDS. For a SDOF system with a natural frequency f
n and a damping ratio ξ, the output pseudovelocity x
pv to an input acceleration x
a can be calculated using a ramp invariant digital filter method [16]:(5)
x
pv
=
F
filter
b
,
a
,
x
a
,
a
=
1
,
-
2
C
,
E
2
,
b
=
-
1
w
2
π
f
i
C
-
1
Q
+
q
S
+
w
,
-
2
C
w
+
1
-
E
2
Q
-
2
q
S
,
E
2
w
+
1
Q
-
C
Q
+
q
S
,
k
1
=
1
-
1
4
Q
2
,
k
2
=
π
Q
f
s
,
k
3
=
2
π
f
s
k
1
,
q
=
1
/
2
Q
-
1
k
1
,
w
=
2
π
f
n
f
s
,
Q
=
1
2
ξ
,
A
=
k
2
f
n
,
B
=
k
3
f
n
,
E
=
exp
⁡
-
A
,
C
=
E
cos
⁡
B
,
S
=
E
sin
⁡
B
,where f
s is the sampling frequency and F
filter indicates filtering the input signal with a ramp invariant digital filter.With the outputx
pv, the cumulative damage can be calculated in both time domain and frequency domain. In time domain, RFCC method is usually used for stress cycle counting. With cyclic numbers of different stress level, S-N curve and Miner’s rule are combined to calculate fatigue damage [22] as follows:(6)
N
i
=
c
S
i
-
b
,
D
t
=
∑
i
=
1
p
n
i
N
i
=
∑
i
=
1
p
n
i
S
i
b
c
,where n
i is the number of cycles exposure at S
i (i
=
1,2
,
…
,
p), p is the number of stress levels considered, N
i is fatigue life at stress S
i, b is the fatigue exponent, and D
t is the total damage index calculated in time domain.Assume that the stress is proportional to pseudovelocity as follows:(7)
S
=
k
x
pv
.The total damage can be calculated as(8)
D
t
=
∑
i
=
1
p
n
i
N
i
=
∑
i
=
1
p
n
i
k
b
x
pv
,
i
b
c
=
k
b
c
∑
i
=
1
p
n
i
x
pv
,
i
b
.In frequency domain, Rayleigh distribution of response stress maxima is assumed and used to calculate the stress cycles(9)
p
S
=
S
σ
S
2
e
-
S
2
/
2
σ
S
2
,where S is the stress value of peaks and σ
S is the RMS of the stress time history.The total damage can be calculated as(10)
D
f
=
n
p
+
T
∫
0
+
∞
p
S
N
S
d
S
=
n
p
+
T
c
∫
0
+
∞
p
S
S
b
d
S
=
f
n
T
c
∫
0
+
∞
p
S
S
b
d
S
,where n
p
+ is the number of positive maxima per unit time in the stress history, which can be replaced by the natural frequency f
n for a lightly damped SDOF system response and T is the total time of exposure to the stress environment.Solving (10) with (9) leads to(11)
D
f
=
f
n
T
c
2
σ
S
2
b
/
2
Γ
1
+
b
2
=
f
n
T
c
k
b
/
2
2
σ
pv
2
b
/
2
Γ
1
+
b
2
,where Γ is the Gamma function and σ
pv is the RMS of pseudovelocity.The RMS of pseudovelocity can be calculated using(12)
σ
pv
=
∫
f
low
f
hi
H
f
2
∗
G
x
x
f
∗
d
f
,where H
f is the transmissibility of a SDOF system (pseudovelocity/acceleration).If the input PSD is relatively flat in the half power band width of each SDOF system, then Mile’s equation can be used to calculateσ
pv in a closed form approximately as follows:(13)
σ
pv
=
G
x
x
f
n
Q
8
π
f
n
.With (11) and (13), the PSD at each natural frequency can be derived from the FDS at the corresponding natural frequency:(14)
G
x
x
f
n
=
D
f
n
Γ
1
+
b
/
2
f
n
T
2
/
b
4
π
f
n
Q
.Note that if the flat spectrum assumption does not hold around some natural frequencies, the FDS calculated with (11) and (12) using the derived PSD in (14) will deviate from the FDS calculated with (11) and (13). In such case, an iteration procedure is needed to take the shape of spectrum into account; see Section 2.4.The ERS shows the damage caused by vibration from an overstress viewpoint. Like the FDS, the ERS is essentially a plot that shows the responses of a number of SDOF systems to a base input acceleration time history. The only difference is that the ERS is generated by calculating the maximum response of a SDOF system to the input. The final plot, the ERS, shows the largest response encountered for a particular SDOF system anywhere within the analyzed time. For a SDOF system with a natural frequencyf
n and a damping ratio ξ, the ERS can also be calculated in both time domain and frequency domain [23] as follows: (15)
ERS
t
=
max
⁡
F
filter
b
,
a
,
x
a
,
(16)
ERS
f
=
C
∗
∫
f
low
f
hi
H
f
2
∗
G
x
x
f
∗
d
f
,
(17)
C
=
2
ln
⁡
f
n
∗
T
+
0.577
2
ln
⁡
f
n
∗
T
.In this paper, the FDS and ERS are used to determine the damage equivalent PSD and test duration.Since the essential idea of the new synthesis method is to model the non-Gaussian signal by a Gaussian signal multiplied by an AMF,(18)
x
nonGaussian
=
x
Gaussian
∗
f
AMF
,the field data must be analyzed to see if it contains many shocks. Shocks must be removed if they are the main reason to cause the non-Gaussianity and processed with separate shock tests. To detect if there are 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 ERS computed from the PSD, which is calculated from the whole time history using Welch method, is used to show where the shocks are [24]:(19)
R
n
,
i
=
20
∗
lg
⁡
y
n
,
i
ERS
n
,
n
=
1,2
,
…
,
N
,
i
=
1,2
,
…
,
M
,where y
n
,
i indicates the maximum value calculated for the i
t
h segment and the n
t
h SDOF system and E
R
S
n indicates the SRS calculated form PSD for the n
t
h SDOF system.


### 2.3. Synthesis of Non-Gaussian Signal

The procedure for synthesizing wind-induced non-Gaussian vibration signal consists of three phases.(
1
) Check the characteristics of measured vibration signal.Calculate the kurtosis of the measured signal to see if it is non-Gaussian signal. Then, calculate the PSD of the measured signal using Welch’s method. Use (19) to detect shocks.(
2
) Synthesize the non-Gaussian random vibration.Create a stationary Gaussian signal using the PSD calculated in phase 1. Calculate the probability distribution functionF
R
M
S of the running RMS of the field data. 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(20)
F
Weibull
x
=
1
-
exp
⁡
-
x
L
.
∧
K
,where x is the random variable, L is the scale parameter, and K is the shape parameter.To find the correctL and K, the probability distribution function of the assumed Weibull distribution F
W
e
i
b
u
l
l and the running RMS F
R
M
S are used to form an objective function using least squares method:(21)
F
obj
=
std
F
RMS
-
F
Weibull
=
std
F
RMS
-
1
-
exp
⁡
-
x
L
i
.
∗
K
i
,where x are the values of the running RMS, depending on how many bins are used to calculate the probability distribution function of the running RMS, i indicates the i
t
h iteration, and s
t
d indicates standard deviation. Note that both F
R
M
S and F
W
e
i
b
u
l
l are vectors.The MATLAB functionfminsearch is used to find the right L and K, so that the standard deviation of error between these two probability distribution functions is minimized.The initial AMF is essentially a vector withN elements following the estimated Weibull distribution:(22)
N
=
T
total
T
c
,where N is the total number of elements in the initial AMF, T
t
o
t
a
l is the total time of the field data, and T
c (correlation time) is the time interval (in second) that is used to construct the AMF.The AMF is then resampled to the same length of the Gaussian signal created from the PSD. The non-Gaussian signal is synthesized by multiplying the AMF with the Gaussian signal. Using the kurtosis of the field data as an objective function, update the AMF from the Weibull distribution until the kurtosis of the field data is reached by the synthesized non-Gaussian signal. The synthesized non-Gaussian signal is scaled to have the same RMS as the field data.(
3
) Compare PSD, PDF, and LCD of the synthesized non-Gaussian signal and the field data.


### 2.4. Derivation of Damage Equivalent Gaussian Signal

To avoid overtesting or undertesting the specimen, both FDS and ERS of the damage equivalent Gaussian signal should be updated so that they match the FDS and ERS of the field data. The procedure of deriving the damage equivalent Gaussian signal from the non-Gaussian signal can be implemented as follows.(
1
) Calculate the FDS from the field measured non-Gaussian signal in time domain, and resample the field data if needed.As pointed out by Lalanne [25], Q does not affect the PSD of derived equivalent Gaussian signal, and a smaller value of b leads to more conservative result. As recommended in MIL-STD 810G [4], Q
=
10 and b
=
4 are used to calculate the pseudovelocity FDS from the field data.(
2
) Calculate the ERS from the time history in time domain.(
3
) Predetermine a test time duration and use FDS from step (
1
) and (14) to derive an initial PSD.(
4
) Calculate the FDS with derived PSD in frequency domain using (11) and (12).(
5
) Use an iteration procedure to update the derived PSD in step (
3
) if the error between two FDS in step (
1
) and step (
4
) is significant(23)
G
m
+
1
f
n
=
G
m
f
n
∗
D
env
f
n
D
m
,
num
f
n
,where G
f
n indicates the PSD, m indicates the iteration number, D
e
n
v
f
n indicates the FDS envelope of field data, and D
m
,
n
u
m
f
n indicates the FDS calculated from the PSD in iteration number m using numerical integration as shown in (11) and (12).(
6
) Calculate the ERS with (updated) PSD in frequency domain using (16) and compare it with that calculated in step (
2
).(
7
) Manipulate the level of PSD and the test time duration to match both FDS and ERS between derived Gaussian signal and measured non-Gaussian signal.


## 2.1. Non-Gaussian Random Signal

If a Gaussian process has a zero mean value and is ergodic, then the PDF of the instantaneous acceleration valuesx
(
t
) is given by the Gaussian distribution with zero mean: (1)
p
x
=
1
2
π
σ
x
exp
⁡
-
x
2
2
σ
x
2
,where σ
x is standard deviation.For large durationsT, the variance σ
x
2 is given by(2)
σ
x
2
≈
1
T
∫
0
T
x
2
d
t
=
∫
0
∞
G
x
x
f
d
f
,where G
x
x
f is the single-sided PSD.It shows that random Gaussian processes with zero mean can be completely described by the PSD function. One useful method for establishing how well a random process can be described by the Gaussian distribution is by computing the higher order moments of the process defined as (continuous and discrete form)(3)
m
n
=
∫
-
∞
∞
x
n
p
x
d
x
=
1
N
∑
j
=
1
N
x
j
n
.When the mean value is zero,(4)
σ
x
2
=
1
N
∑
j
=
1
N
x
j
2
=
m
2
,
S
=
1
N
∑
j
=
1
N
x
j
3
σ
x
3
=
m
3
σ
x
3
,
K
=
1
N
∑
j
=
1
N
x
j
4
σ
x
4
=
m
4
σ
x
4
,where S is the skewness and K is the kurtosis.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.


## 2.2. Response Spectra

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 that are 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 Figures1 and 2.Figure 1
SDOF systems for calculation of FDS.Figure 2
Calculation procedure of the FDS.It has been shown that the stress is roughly proportional to pseudovelocity [21]. Thus, it is more reasonable to calculate the pseudovelocity FDS. For a SDOF system with a natural frequency f
n and a damping ratio ξ, the output pseudovelocity x
pv to an input acceleration x
a can be calculated using a ramp invariant digital filter method [16]:(5)
x
pv
=
F
filter
b
,
a
,
x
a
,
a
=
1
,
-
2
C
,
E
2
,
b
=
-
1
w
2
π
f
i
C
-
1
Q
+
q
S
+
w
,
-
2
C
w
+
1
-
E
2
Q
-
2
q
S
,
E
2
w
+
1
Q
-
C
Q
+
q
S
,
k
1
=
1
-
1
4
Q
2
,
k
2
=
π
Q
f
s
,
k
3
=
2
π
f
s
k
1
,
q
=
1
/
2
Q
-
1
k
1
,
w
=
2
π
f
n
f
s
,
Q
=
1
2
ξ
,
A
=
k
2
f
n
,
B
=
k
3
f
n
,
E
=
exp
⁡
-
A
,
C
=
E
cos
⁡
B
,
S
=
E
sin
⁡
B
,where f
s is the sampling frequency and F
filter indicates filtering the input signal with a ramp invariant digital filter.With the outputx
pv, the cumulative damage can be calculated in both time domain and frequency domain. In time domain, RFCC method is usually used for stress cycle counting. With cyclic numbers of different stress level, S-N curve and Miner’s rule are combined to calculate fatigue damage [22] as follows:(6)
N
i
=
c
S
i
-
b
,
D
t
=
∑
i
=
1
p
n
i
N
i
=
∑
i
=
1
p
n
i
S
i
b
c
,where n
i is the number of cycles exposure at S
i (i
=
1,2
,
…
,
p), p is the number of stress levels considered, N
i is fatigue life at stress S
i, b is the fatigue exponent, and D
t is the total damage index calculated in time domain.Assume that the stress is proportional to pseudovelocity as follows:(7)
S
=
k
x
pv
.The total damage can be calculated as(8)
D
t
=
∑
i
=
1
p
n
i
N
i
=
∑
i
=
1
p
n
i
k
b
x
pv
,
i
b
c
=
k
b
c
∑
i
=
1
p
n
i
x
pv
,
i
b
.In frequency domain, Rayleigh distribution of response stress maxima is assumed and used to calculate the stress cycles(9)
p
S
=
S
σ
S
2
e
-
S
2
/
2
σ
S
2
,where S is the stress value of peaks and σ
S is the RMS of the stress time history.The total damage can be calculated as(10)
D
f
=
n
p
+
T
∫
0
+
∞
p
S
N
S
d
S
=
n
p
+
T
c
∫
0
+
∞
p
S
S
b
d
S
=
f
n
T
c
∫
0
+
∞
p
S
S
b
d
S
,where n
p
+ is the number of positive maxima per unit time in the stress history, which can be replaced by the natural frequency f
n for a lightly damped SDOF system response and T is the total time of exposure to the stress environment.Solving (10) with (9) leads to(11)
D
f
=
f
n
T
c
2
σ
S
2
b
/
2
Γ
1
+
b
2
=
f
n
T
c
k
b
/
2
2
σ
pv
2
b
/
2
Γ
1
+
b
2
,where Γ is the Gamma function and σ
pv is the RMS of pseudovelocity.The RMS of pseudovelocity can be calculated using(12)
σ
pv
=
∫
f
low
f
hi
H
f
2
∗
G
x
x
f
∗
d
f
,where H
f is the transmissibility of a SDOF system (pseudovelocity/acceleration).If the input PSD is relatively flat in the half power band width of each SDOF system, then Mile’s equation can be used to calculateσ
pv in a closed form approximately as follows:(13)
σ
pv
=
G
x
x
f
n
Q
8
π
f
n
.With (11) and (13), the PSD at each natural frequency can be derived from the FDS at the corresponding natural frequency:(14)
G
x
x
f
n
=
D
f
n
Γ
1
+
b
/
2
f
n
T
2
/
b
4
π
f
n
Q
.Note that if the flat spectrum assumption does not hold around some natural frequencies, the FDS calculated with (11) and (12) using the derived PSD in (14) will deviate from the FDS calculated with (11) and (13). In such case, an iteration procedure is needed to take the shape of spectrum into account; see Section 2.4.The ERS shows the damage caused by vibration from an overstress viewpoint. Like the FDS, the ERS is essentially a plot that shows the responses of a number of SDOF systems to a base input acceleration time history. The only difference is that the ERS is generated by calculating the maximum response of a SDOF system to the input. The final plot, the ERS, shows the largest response encountered for a particular SDOF system anywhere within the analyzed time. For a SDOF system with a natural frequencyf
n and a damping ratio ξ, the ERS can also be calculated in both time domain and frequency domain [23] as follows: (15)
ERS
t
=
max
⁡
F
filter
b
,
a
,
x
a
,
(16)
ERS
f
=
C
∗
∫
f
low
f
hi
H
f
2
∗
G
x
x
f
∗
d
f
,
(17)
C
=
2
ln
⁡
f
n
∗
T
+
0.577
2
ln
⁡
f
n
∗
T
.In this paper, the FDS and ERS are used to determine the damage equivalent PSD and test duration.Since the essential idea of the new synthesis method is to model the non-Gaussian signal by a Gaussian signal multiplied by an AMF,(18)
x
nonGaussian
=
x
Gaussian
∗
f
AMF
,the field data must be analyzed to see if it contains many shocks. Shocks must be removed if they are the main reason to cause the non-Gaussianity and processed with separate shock tests. To detect if there are 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 ERS computed from the PSD, which is calculated from the whole time history using Welch method, is used to show where the shocks are [24]:(19)
R
n
,
i
=
20
∗
lg
⁡
y
n
,
i
ERS
n
,
n
=
1,2
,
…
,
N
,
i
=
1,2
,
…
,
M
,where y
n
,
i indicates the maximum value calculated for the i
t
h segment and the n
t
h SDOF system and E
R
S
n indicates the SRS calculated form PSD for the n
t
h SDOF system.


## 2.3. Synthesis of Non-Gaussian Signal

The procedure for synthesizing wind-induced non-Gaussian vibration signal consists of three phases.(
1
) Check the characteristics of measured vibration signal.Calculate the kurtosis of the measured signal to see if it is non-Gaussian signal. Then, calculate the PSD of the measured signal using Welch’s method. Use (19) to detect shocks.(
2
) Synthesize the non-Gaussian random vibration.Create a stationary Gaussian signal using the PSD calculated in phase 1. Calculate the probability distribution functionF
R
M
S of the running RMS of the field data. 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(20)
F
Weibull
x
=
1
-
exp
⁡
-
x
L
.
∧
K
,where x is the random variable, L is the scale parameter, and K is the shape parameter.To find the correctL and K, the probability distribution function of the assumed Weibull distribution F
W
e
i
b
u
l
l and the running RMS F
R
M
S are used to form an objective function using least squares method:(21)
F
obj
=
std
F
RMS
-
F
Weibull
=
std
F
RMS
-
1
-
exp
⁡
-
x
L
i
.
∗
K
i
,where x are the values of the running RMS, depending on how many bins are used to calculate the probability distribution function of the running RMS, i indicates the i
t
h iteration, and s
t
d indicates standard deviation. Note that both F
R
M
S and F
W
e
i
b
u
l
l are vectors.The MATLAB functionfminsearch is used to find the right L and K, so that the standard deviation of error between these two probability distribution functions is minimized.The initial AMF is essentially a vector withN elements following the estimated Weibull distribution:(22)
N
=
T
total
T
c
,where N is the total number of elements in the initial AMF, T
t
o
t
a
l is the total time of the field data, and T
c (correlation time) is the time interval (in second) that is used to construct the AMF.The AMF is then resampled to the same length of the Gaussian signal created from the PSD. The non-Gaussian signal is synthesized by multiplying the AMF with the Gaussian signal. Using the kurtosis of the field data as an objective function, update the AMF from the Weibull distribution until the kurtosis of the field data is reached by the synthesized non-Gaussian signal. The synthesized non-Gaussian signal is scaled to have the same RMS as the field data.(
3
) Compare PSD, PDF, and LCD of the synthesized non-Gaussian signal and the field data.


## 2.4. Derivation of Damage Equivalent Gaussian Signal

To avoid overtesting or undertesting the specimen, both FDS and ERS of the damage equivalent Gaussian signal should be updated so that they match the FDS and ERS of the field data. The procedure of deriving the damage equivalent Gaussian signal from the non-Gaussian signal can be implemented as follows.(
1
) Calculate the FDS from the field measured non-Gaussian signal in time domain, and resample the field data if needed.As pointed out by Lalanne [25], Q does not affect the PSD of derived equivalent Gaussian signal, and a smaller value of b leads to more conservative result. As recommended in MIL-STD 810G [4], Q
=
10 and b
=
4 are used to calculate the pseudovelocity FDS from the field data.(
2
) Calculate the ERS from the time history in time domain.(
3
) Predetermine a test time duration and use FDS from step (
1
) and (14) to derive an initial PSD.(
4
) Calculate the FDS with derived PSD in frequency domain using (11) and (12).(
5
) Use an iteration procedure to update the derived PSD in step (
3
) if the error between two FDS in step (
1
) and step (
4
) is significant(23)
G
m
+
1
f
n
=
G
m
f
n
∗
D
env
f
n
D
m
,
num
f
n
,where G
f
n indicates the PSD, m indicates the iteration number, D
e
n
v
f
n indicates the FDS envelope of field data, and D
m
,
n
u
m
f
n indicates the FDS calculated from the PSD in iteration number m using numerical integration as shown in (11) and (12).(
6
) Calculate the ERS with (updated) PSD in frequency domain using (16) and compare it with that calculated in step (
2
).(
7
) Manipulate the level of PSD and the test time duration to match both FDS and ERS between derived Gaussian signal and measured non-Gaussian signal.


## 3. Case Study

### 3.1. Synthesis of Non-Gaussian Signal

In this case study, the field measured data is used to validate the new method for synthesizing the non-Gaussian vibration signal. The random acceleration signal induced by wind is taken from a test item installed on a mast, as shown in Figure3, along with the running RMS. The test item is a dummy box containing electronics. The dummy box was equipped with B&K WB0179 triaxial accelerometer set and a wind speed measuring device. The box was a standard Ericsson Micro Radio Base equipment, manufactured in cast lightweight alloy with integrated heat sinks. The dimensions are about 530
∗
400
∗
185 mm and the weight was about 21 kg. The back side of the dummy box was mounted to the mast at a height of 50 meters. The running RMS is calculated using 10 s frames with 50% overlapping.Figure 3
Field data and the running RMS.This is a non-Gaussian signal with kurtosis equaling 9.4262. Equation (19) is used to detect shocks (see Figure 4). The ratio is shown in the third dimension by the color. As we can see from Figure 4, there are only few shocks superimposed on the underlying random vibration signal, which means that the non-Gaussianity is caused by the running RMS.Figure 4
Significant shock detection.The PSD of the field data is calculated using Welch method with 4096-point Hanning window and 50% overlapping. A Gaussian signal is generated using the same PSD. The PDF of both the field data and the synthesized Gaussian signal are shown in Figure9. As we can see, the PDF of the field data is obviously different from that of the synthesized Gaussian signal. The PDF of running RMS is calculated using different number of bins and is shown in Figure 5. From Figure 5, we can see that despite of the numbers of bins used 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.Figure 5
PDF of the running RMS using different number of bins and Rayleigh distribution.Iteration method described in Section2.3 (phase 2) is used to determine the L and K, as shown in Figure 6.Figure 6
The iteration to determineL and K.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 Section2.3 (phase 2). The synthesized non-Gaussian signal is shown in Figure 7, together with the field data (the time duration of the synthesized non-Gaussian signal is adjustable).Figure 7
Field data and synthesized non-Gaussian signal.The PSD and PDF of the field data and the synthesized non-Gaussian signal are very similar (see Figures8 and 9), which proves the effectiveness of the new method.Figure 8
The PSD of the field data and the synthesized non-Gaussian signal.Figure 9
The PDF of the field data and synthesized Gaussian and non-Gaussian signal.For comparison, the cycle counting is performed on the field data, the synthesized Gaussian, and non-Gaussian signal, using the RFCC and Dirlik method. The results are shown in Figure10. 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.Figure 10
Number of cycles versus data ranges for RFCC and Dirlik.Care must be taken when creating the AMF. The AMF can be generated by using different correlation time and being resampled. The AMFs using correlation time of 0.1 s, 10 s, and 100 s are shown in Figure11. The corresponding synthesized time signals are shown in Figure 12. From Figure 12, we can see that inappropriate choice of the correlation time (0.1 s and 100 s) leads to dissimilar time histories to the field data.Figure 11
Time history of field data and synthesized AMF.Figure 12
Time history of field data and synthesized data.The PSD are calculated for the field data and the synthesized non-Gaussian signals using different correlation time (see Figure13). From Figure 13, we can see that the PSD of the synthesized non-Gaussian signal with a correlation time of 0.1 s differs from others at low frequencies (between 1 Hz and 10 Hz). This is caused by new significant frequency components introduced by the AMF when it is created. This indicates that small correlation time may lead to inaccurate PSD. The PDF are calculated for the field data, synthesized Gaussian, and non-Gaussian signals using different correlation time (see Figure 14). From Figure 14, we can see that the PDF of the synthesized non-Gaussian signal using 100 s correlation time shows different characteristics from other non-Gaussian signals. The middle of this PDF is shaper and the tails are narrower. This indicates that large correlation time may lead to inaccurate PDF.Figure 13
The PSD of field signal and synthesized signal.Figure 14
PDF for field data and synthesized data.


### 3.2. Derivation of Damage Equivalent Gaussian Signal

In this case study, the same non-Gaussian signal from Section3.1 is used to validate the procedure of generating damage equivalent Gaussian signal. Figure 15 shows the pseudovelocity FDS calculated from the field data.Figure 15
The FDS of the field data and the equivalent PSD (T = 1200 s).Using the same time duration as the field data (1200 s), an equivalent Gaussian PSD is derived from the FDS using (14). To account for the shape of the PSD, the FDS is calculated again using (11) and (12) and compared with the FDS that is used to derive the PSD. Using the derived PSD as an initial value, (23) is used to update the derived PSD. The FDS calculated from the updated PSD is shown in Figure 15. As we can see from Figure 15, a better match between the FDS of the field data and the equivalent PSD is reached after iteration. The ERS is calculated from the field data and the equivalent PSD, as shown in Figure 16. From Figure 16, we can see that the initial derived equivalent Gaussian signal is not severe enough to match the ERS of the field data.Figure 16
The ERS of the field data and the equivalent PSD (T = 1200 s).The test time duration is further reduced to obtain an appropriate level of PSD so that the ERS calculated from it matches the ERS calculated from the field data, as shown in Figure17. The PSD of derived equivalent Gaussian signals using different duration time is shown in Figure 18. As we can see, the PSD level increases as the duration time decreases so that the FDS remains the same, as shown in Figure 19. Note that although the PSD level can be further increased to reduce the test time, the test item may fail due to peak response which will not be experienced in field.Figure 17
The ERS of the field data and the equivalent PSD (T = 120 s).Figure 18
PSD of equivalent Gaussian signal with different duration time.Figure 19
The FDS of the field data and the equivalent PSD (T = 120 s).


## 3.1. Synthesis of Non-Gaussian Signal

In this case study, the field measured data is used to validate the new method for synthesizing the non-Gaussian vibration signal. The random acceleration signal induced by wind is taken from a test item installed on a mast, as shown in Figure3, along with the running RMS. The test item is a dummy box containing electronics. The dummy box was equipped with B&K WB0179 triaxial accelerometer set and a wind speed measuring device. The box was a standard Ericsson Micro Radio Base equipment, manufactured in cast lightweight alloy with integrated heat sinks. The dimensions are about 530
∗
400
∗
185 mm and the weight was about 21 kg. The back side of the dummy box was mounted to the mast at a height of 50 meters. The running RMS is calculated using 10 s frames with 50% overlapping.Figure 3
Field data and the running RMS.This is a non-Gaussian signal with kurtosis equaling 9.4262. Equation (19) is used to detect shocks (see Figure 4). The ratio is shown in the third dimension by the color. As we can see from Figure 4, there are only few shocks superimposed on the underlying random vibration signal, which means that the non-Gaussianity is caused by the running RMS.Figure 4
Significant shock detection.The PSD of the field data is calculated using Welch method with 4096-point Hanning window and 50% overlapping. A Gaussian signal is generated using the same PSD. The PDF of both the field data and the synthesized Gaussian signal are shown in Figure9. As we can see, the PDF of the field data is obviously different from that of the synthesized Gaussian signal. The PDF of running RMS is calculated using different number of bins and is shown in Figure 5. From Figure 5, we can see that despite of the numbers of bins used 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.Figure 5
PDF of the running RMS using different number of bins and Rayleigh distribution.Iteration method described in Section2.3 (phase 2) is used to determine the L and K, as shown in Figure 6.Figure 6
The iteration to determineL and K.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 Section2.3 (phase 2). The synthesized non-Gaussian signal is shown in Figure 7, together with the field data (the time duration of the synthesized non-Gaussian signal is adjustable).Figure 7
Field data and synthesized non-Gaussian signal.The PSD and PDF of the field data and the synthesized non-Gaussian signal are very similar (see Figures8 and 9), which proves the effectiveness of the new method.Figure 8
The PSD of the field data and the synthesized non-Gaussian signal.Figure 9
The PDF of the field data and synthesized Gaussian and non-Gaussian signal.For comparison, the cycle counting is performed on the field data, the synthesized Gaussian, and non-Gaussian signal, using the RFCC and Dirlik method. The results are shown in Figure10. 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.Figure 10
Number of cycles versus data ranges for RFCC and Dirlik.Care must be taken when creating the AMF. The AMF can be generated by using different correlation time and being resampled. The AMFs using correlation time of 0.1 s, 10 s, and 100 s are shown in Figure11. The corresponding synthesized time signals are shown in Figure 12. From Figure 12, we can see that inappropriate choice of the correlation time (0.1 s and 100 s) leads to dissimilar time histories to the field data.Figure 11
Time history of field data and synthesized AMF.Figure 12
Time history of field data and synthesized data.The PSD are calculated for the field data and the synthesized non-Gaussian signals using different correlation time (see Figure13). From Figure 13, we can see that the PSD of the synthesized non-Gaussian signal with a correlation time of 0.1 s differs from others at low frequencies (between 1 Hz and 10 Hz). This is caused by new significant frequency components introduced by the AMF when it is created. This indicates that small correlation time may lead to inaccurate PSD. The PDF are calculated for the field data, synthesized Gaussian, and non-Gaussian signals using different correlation time (see Figure 14). From Figure 14, we can see that the PDF of the synthesized non-Gaussian signal using 100 s correlation time shows different characteristics from other non-Gaussian signals. The middle of this PDF is shaper and the tails are narrower. This indicates that large correlation time may lead to inaccurate PDF.Figure 13
The PSD of field signal and synthesized signal.Figure 14
PDF for field data and synthesized data.


## 3.2. Derivation of Damage Equivalent Gaussian Signal

In this case study, the same non-Gaussian signal from Section3.1 is used to validate the procedure of generating damage equivalent Gaussian signal. Figure 15 shows the pseudovelocity FDS calculated from the field data.Figure 15
The FDS of the field data and the equivalent PSD (T = 1200 s).Using the same time duration as the field data (1200 s), an equivalent Gaussian PSD is derived from the FDS using (14). To account for the shape of the PSD, the FDS is calculated again using (11) and (12) and compared with the FDS that is used to derive the PSD. Using the derived PSD as an initial value, (23) is used to update the derived PSD. The FDS calculated from the updated PSD is shown in Figure 15. As we can see from Figure 15, a better match between the FDS of the field data and the equivalent PSD is reached after iteration. The ERS is calculated from the field data and the equivalent PSD, as shown in Figure 16. From Figure 16, we can see that the initial derived equivalent Gaussian signal is not severe enough to match the ERS of the field data.Figure 16
The ERS of the field data and the equivalent PSD (T = 1200 s).The test time duration is further reduced to obtain an appropriate level of PSD so that the ERS calculated from it matches the ERS calculated from the field data, as shown in Figure17. The PSD of derived equivalent Gaussian signals using different duration time is shown in Figure 18. As we can see, the PSD level increases as the duration time decreases so that the FDS remains the same, as shown in Figure 19. Note that although the PSD level can be further increased to reduce the test time, the test item may fail due to peak response which will not be experienced in field.Figure 17
The ERS of the field data and the equivalent PSD (T = 120 s).Figure 18
PSD of equivalent Gaussian signal with different duration time.Figure 19
The FDS of the field data and the equivalent PSD (T = 120 s).


## 4. Conclusions

Two methodologies are presented for shaker simulation of wind-induced non-Gaussian vibration. The first methodology utilizes a new method to produce the non-Gaussian signal offline and reproduce it on the shaker in a Time Waveform Replication mode. The second methodology derives a damage equivalent Gaussian signal from the non-Gaussian signal so that the traditional closed-loop frequency control mode can be easily implemented. Conclusions are as below.(1)
The simulation method presented in the first methodology is limited to the case of the running RMS-induced non-Gaussian signal. Using this method, the synthesized non-Gaussian signal has the same PSD, PDF, and LCD as the field data.
(2)
Pseudovelocity FDS and ERS can be used to assess the severity of the field measured non-Gaussian data and to derive a stationary test specification characterized by only PSD and test duration. The derived PSD can be improved by taking its spectrum shape into consideration using iteration method.
(3)
The first methodology applies when test time reduction is not needed, whereas the second methodology applies when accelerated testing is required.



---

*Source: 5450865-2016-01-10.xml*