# Optimization of Damage Equivalent Accelerated Test Spectrum Derivation Using Multiple Non-Gaussian Vibration Data

Authors: Fei Xu; Kjell Ahlin; Binyi Wang
Journal: Journal of Sensors
 (2021)

---

## Abstract

The response spectra are widely used in the damage assessment of non-Gaussian random vibration environments and the derivation of damage equivalent accelerated test spectrum. The effectiveness of the latter is strongly affected by modal parameter uncertainties, multiple field data processing, and the nonsmooth shape of the derived power spectral density (PSD). Optimization of accelerated test spectrum derivation based on dynamic parameter selection and iterative update of spectrum envelope is presented in this paper. The extreme response spectrum (ERS) envelope of the field data is firstly taken as the limiting spectrum, and the corresponding relationship between damping coefficient, fatigue exponent, and damage equivalent PSD under different test times is constructed to achieve the dynamic selection of uncertain parameters in the response spectrum model. Then, an iterative update model based on the weighted sum of fatigue damage spectrum (FDS) error is presented to reduce the error introduced by the nonsmooth shape of the derived PSD. The case study shows that undertest can be effectively avoided by the dynamic selection of model parameters. The weighted error is reduced from 80.1% to 7.5% after 7 iterations. Particularly, the error is close to 0 within the peak and valley frequency band.


---

## Body

## 1. Introduction

Random vibration tests are widely used to identify defects in product design and to evaluate the reliability and fatigue life of products. To simulate the mechanical environment experienced by products during their life cycle in the lab, the test spectrum is usually derived from the field measured data, which is collected using different sensors (e.g., ultrasonic transducer and accelerometer) [1]. With multiple and complicated field data, the accuracy of signal processing algorithms and data analysis methods are key factors when understanding the mechanical environment and facilitating the test spectrum derivation.Power spectral density (PSD) is used by typical digital shaker controllers to describe the random vibration environment experienced by products, with an implicit assumption that the vibration data follows a Gaussian distribution. However, field recorded time histories, e.g., wind-induced vibration data [2] and road roughness-induced vehicle vibration data [3], usually present non-Gaussian characteristics. Since non-Gaussian (especially super-Gaussian) random vibration tends to shorten the fatigue life of products [4, 5], different vibration control methods were studied to simulate the non-Gaussian vibration environment in the lab. Steinwolf presented a phase manipulation method to simulate the non-Gaussian data [6]. In particular, the selected phase is transformed from random to deterministic in order to obtain a prescribed kurtosis. PSD and probability density function (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 [7]. From the perspective of time-varying PSD and PDF, a non-stationary non-Gaussian stochastic process simulation method based on the zero-memory nonlinear translation relationship between non-Gaussian and Gaussian stochastic processes is proposed by Cui et al. [8]. Fei et al. [9] presented a method 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. Zheng et al. [10] presented a control method for multi-input multioutput stationary non-Gaussian random vibration test using time-domain randomization. The frequency-domain inverse system method is utilized to obtain the desired drive signals for dynamic inputs, and the skewness, kurtosis, and PSD of response signals are controlled efficiently and simultaneously.The major drawback of these non-Gaussian vibration control methods, when used for fatigue life evaluation of high reliability and long-life products, is the long test time and high cost. To address this problem, the response spectrum method was utilized to evaluate the potential damage introduced by the vibration environment. An accelerated test spectrum was then synthesized from the response spectra based on the damage equivalence criterion [11–13]. To perform life cycle potential damage assessment and test tailoring, Lalanne [14] derived an equivalent PSD from multiple task profiles based on the frequency domain fatigue damage spectrum (FDS). Wijker [15] and Decker et al. [16] derived the accelerated test spectrum based on the extreme response spectrum (ERS). The effect of parameter uncertainties (damping coefficient and fatigue exponent) on the root mean square (RMS) value of the derived PSD was studied. Steinwolf and Wolfsteiner [17] indicated that to avoid changes in the failure mechanism and dynamic properties of products, the ERS of field data should be used to limit the accelerated PSD level. Xu et al. [18] showed that the FDS of non-Gaussian field data increases as the kurtosis increases. A significant error would be introduced if FDS and ERS were calculated in the frequency domain directly using PSD. Ahlin [19] calculated FDS and ERS in time domain using the ramp invariant digital filter method. Due to the low computational efficiency of FDS in the time domain, Wolfsteiner [20] decomposed the non-Gaussian load into a combination of Gaussian loads, so that the frequency domain method can be used to improve the computational efficiency. Assuming that the synthesized equivalent PSD is flat within the half-power bandwidth, Cianetti et al. [21] presented the equivalent PSD synthesis method for non-Gaussian data using FDS and Mile’s Equation. These existing methods have demonstrated how to derive the equivalent accelerated test spectrum from a single set of non-Gaussian data and analyzed the influence of parameter uncertainties on the RMS value of test spectrum. However, how to select the correct model parameters using multiple non-Gaussian vibration data, so as to derive the equivalent accelerated test spectrum which can effectively avoid undertest, has not been studied. In addition, how to construct the objective function, so as to rapidly reduce the error introduced by nonsmooth shape of the synthesized PSD, has not been demonstrated. To address these problems, an optimization design procedure of test spectrum derivation based on dynamic parameter selection and iterative spectrum envelope update technique is presented in this paper.The remainder of the paper is organized as follows. Section2 introduces the definition of Gaussian random vibration, non-Gaussian random vibration, FDS, and ERS. The optimized accelerated test spectrum derivation procedure is also presented. In Section 3, validation with field measured non-Gaussian vibration data is presented. Finally, the discussion and conclusions are summarized in Section 4.


## 2. Materials and Methods

### 2.1. Gaussian and Non-Gaussian Random Vibration

For the field measured discrete random vibration signalxt, the higher-order moments can be used to determine if it follows the Gaussian distribution:(1)mn=∫−∞∞xnpxdx=1N∑j=1Nxjn,where mn is the nth order moments, px is the probability density function (PDF), and N is the number of signals.When the mean value is zero,(2)S=1N∑j=1Nxj3σx3=m3σx3,K=1N∑j=1Nxj4σx4=m4σx4,where S is the skewness, K is the kurtosis, and σx is the root mean square (RMS) of xt.For a Gaussian distributed signal, the skewness is 0 and the kurtosis is 3. Deviation from these values indicates that the signal follows a non-Gaussian distribution. From the perspective of response spectrum calculation of a non-Gaussian signal, kurtosis is more important than skewness, since it represents the probability of peak values in time history [22].


### 2.2. Fatigue Damage Spectrum (FDS) and Extreme Response Spectrum (ERS)

The FDS and ERS are widely used response spectra for evaluating potential damage under random vibration environment. The FDS is essentially a plot that shows the response of a series of single-degree-of-freedom (SDOF) systems to the base input acceleration time history. Many SDOF systems tuned to a range of natural frequencies are assessed using the same input. The FDS shows the fatigue damage encountered for a particular SDOF system anywhere within the analyzed time. It has been shown that the stress is roughly proportional to pseudovelocity [23]. For a SDOF system with a natural frequency fn and a damping ratio ξ, the output pseudovelocity xpv to an input acceleration xa can be calculated using a ramp invariant digital filter method [19]:(3)xpv=Ffilterb,a,xa,a=1,−2C,E2,b=−1w2πfnC−1Q+qS+w,−2Cw+1−E2Q−2qS,E2w+1Q−CQ+qS,w=2πfnfs,Q=12ζ,k1=1−14Q2,k2=πQfs,k3=2πfsk1,q=1/2/Q/Q−1k1,A=k2fn,B=k3fn,E=exp−A,C=EcosB,S=EsinB,where fs is the sampling frequency and Ffilter indicates filtering the input signal with a ramp invariant digital filter.With the outputxpv, the cumulative damage can be calculated in both time domain and frequency domain. In time domain, the rain flow cycle counting method is usually used for stress cycle counting. With cyclic numbers of different stress levels, the S-N curve and Miner’s rule are combined to calculate the fatigue damage as follows [24]:(4)Ni=cSi−b,Dt=∑i=1pni/Ni=∑i=1pnikbxpv,ib/c=kb/c∗∑i=1pnixpv,ib,where Ni is the fatigue life at stress Si (i=1,2,⋯,p), p is the number of stress levels considered, ni is the number of cycle exposure at Si, c is a constant, b is the fatigue exponent, k is the proportional coefficient of stress to pseudovelocity, Dt is the total damage index calculated in time domain, and xpv,i is the output pseudovelocity at ith cycle.In the frequency domain, Rayleigh distribution of response stress maxima is assumed and used to calculate the stress cycles:(5)pS=SσS2e−S2/2σS2,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(6)Df=fnTc∫0+∞pSSbdS,where T is the total time of exposure to the stress environment.Substituting Equation (5) into Equation (6) leads to(7)Df=fnTckb2σpv2b/2Γ1+b2,where Γ is the gamma function and σpv is the RMS of pseudovelocity.The RMS of pseudovelocity can be calculated using(8)σpv=∑f=flowf=fhiHf2∗Gxxf∗df,where Hf is the transmissibility of a SDOF system (pseudovelocity/acceleration), Gxxf is the input PSD, flow is the lower limit, and fhi is the higher limit.If the input PSD is relatively flat in the half-power bandwidth of each SDOF system, then Mile’s Equation can be used to calculateσpv in a closed form approximately as follows:(9)σpv=GxxfnQ8πfn.With Equations (7) and (9), the PSD at each natural frequency can be derived from the FDS at the corresponding natural frequency:(10)Gxxfn=DfnΓ1+b/2fnT2/b4πfnQ.Note that if the flat PSD assumption does not hold around some natural frequencies, the FDS calculated with Equations (7) and (8) using the derived PSD in Equation (10) will deviate from the FDS calculated with Equations (7) and (9). In such case, an iterative update procedure is needed to take the shape of PSD into account (see Section 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 frequencyfn and a damping ratio ξ, the ERS can also be calculated in both time domain and frequency domain as follows [14]:(11)ERSt=maxFfilterb,a,xa,(12)ERSf=2lnfn∗T+0.5772lnfn∗T∗∑f=flowf=fhiHf^2∗Gxxf∗df.


### 2.3. Derivation of Damage Equivalent Accelerated Test Spectrum

Since the accuracy of the damage equivalent accelerated test spectrum is strongly affected by uncertain parameters, dispersion of multiple non-Gaussian samples, and the smoothness of the derived PSD, an optimum design procedure is presented as follows:(1)
Calculate the FDS of multiple samples in time domain and take the envelope using different damping ratios (denoted by quality factorQ) and fatigue exponent b(2)
Predetermine the accelerated test time, and synthesize the initial PSD using the FDS envelope and Equation (10)(3)
ChooseQ and b in a dynamic way, so that the most conservative PSD is derived(4)
Calculate the FDS in the frequency domain using the derived PSD, Equation (7), and Equation (8)(5)
Compare the calculated FDS in Step (1) and Step (4), and define the weighted sum of FDS errors:(13)Werror=∑fi=flowestfhiestDnumfn−DenvfnDenvfn∗100∗Denvfn∑fi=flowestfhiestDenvfn,where flowest is the lowest natural frequency, fhiest is the highest natural frequency, Denvfn is the FDS envelope calculated in Step (1), and Dnumfn is the FDS calculated in Step (4)(6)
Update the PSD in Step (3) until a certain error is reached:(14)Gm+1fn=Gmfn∗DenvfnDnumfn2/b,where Gmfn is the damage equivalent PSD after m iterations(7)
Calculate the ERS in frequency domain using the updated PSD and Equation (12) and compare with the ERS envelope calculated in time domain(8)
Update the predetermined test time and PSD level to match both the FDS and ERS envelope of the field data


## 2.1. Gaussian and Non-Gaussian Random Vibration

For the field measured discrete random vibration signalxt, the higher-order moments can be used to determine if it follows the Gaussian distribution:(1)mn=∫−∞∞xnpxdx=1N∑j=1Nxjn,where mn is the nth order moments, px is the probability density function (PDF), and N is the number of signals.When the mean value is zero,(2)S=1N∑j=1Nxj3σx3=m3σx3,K=1N∑j=1Nxj4σx4=m4σx4,where S is the skewness, K is the kurtosis, and σx is the root mean square (RMS) of xt.For a Gaussian distributed signal, the skewness is 0 and the kurtosis is 3. Deviation from these values indicates that the signal follows a non-Gaussian distribution. From the perspective of response spectrum calculation of a non-Gaussian signal, kurtosis is more important than skewness, since it represents the probability of peak values in time history [22].


## 2.2. Fatigue Damage Spectrum (FDS) and Extreme Response Spectrum (ERS)

The FDS and ERS are widely used response spectra for evaluating potential damage under random vibration environment. The FDS is essentially a plot that shows the response of a series of single-degree-of-freedom (SDOF) systems to the base input acceleration time history. Many SDOF systems tuned to a range of natural frequencies are assessed using the same input. The FDS shows the fatigue damage encountered for a particular SDOF system anywhere within the analyzed time. It has been shown that the stress is roughly proportional to pseudovelocity [23]. For a SDOF system with a natural frequency fn and a damping ratio ξ, the output pseudovelocity xpv to an input acceleration xa can be calculated using a ramp invariant digital filter method [19]:(3)xpv=Ffilterb,a,xa,a=1,−2C,E2,b=−1w2πfnC−1Q+qS+w,−2Cw+1−E2Q−2qS,E2w+1Q−CQ+qS,w=2πfnfs,Q=12ζ,k1=1−14Q2,k2=πQfs,k3=2πfsk1,q=1/2/Q/Q−1k1,A=k2fn,B=k3fn,E=exp−A,C=EcosB,S=EsinB,where fs is the sampling frequency and Ffilter indicates filtering the input signal with a ramp invariant digital filter.With the outputxpv, the cumulative damage can be calculated in both time domain and frequency domain. In time domain, the rain flow cycle counting method is usually used for stress cycle counting. With cyclic numbers of different stress levels, the S-N curve and Miner’s rule are combined to calculate the fatigue damage as follows [24]:(4)Ni=cSi−b,Dt=∑i=1pni/Ni=∑i=1pnikbxpv,ib/c=kb/c∗∑i=1pnixpv,ib,where Ni is the fatigue life at stress Si (i=1,2,⋯,p), p is the number of stress levels considered, ni is the number of cycle exposure at Si, c is a constant, b is the fatigue exponent, k is the proportional coefficient of stress to pseudovelocity, Dt is the total damage index calculated in time domain, and xpv,i is the output pseudovelocity at ith cycle.In the frequency domain, Rayleigh distribution of response stress maxima is assumed and used to calculate the stress cycles:(5)pS=SσS2e−S2/2σS2,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(6)Df=fnTc∫0+∞pSSbdS,where T is the total time of exposure to the stress environment.Substituting Equation (5) into Equation (6) leads to(7)Df=fnTckb2σpv2b/2Γ1+b2,where Γ is the gamma function and σpv is the RMS of pseudovelocity.The RMS of pseudovelocity can be calculated using(8)σpv=∑f=flowf=fhiHf2∗Gxxf∗df,where Hf is the transmissibility of a SDOF system (pseudovelocity/acceleration), Gxxf is the input PSD, flow is the lower limit, and fhi is the higher limit.If the input PSD is relatively flat in the half-power bandwidth of each SDOF system, then Mile’s Equation can be used to calculateσpv in a closed form approximately as follows:(9)σpv=GxxfnQ8πfn.With Equations (7) and (9), the PSD at each natural frequency can be derived from the FDS at the corresponding natural frequency:(10)Gxxfn=DfnΓ1+b/2fnT2/b4πfnQ.Note that if the flat PSD assumption does not hold around some natural frequencies, the FDS calculated with Equations (7) and (8) using the derived PSD in Equation (10) will deviate from the FDS calculated with Equations (7) and (9). In such case, an iterative update procedure is needed to take the shape of PSD into account (see Section 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 frequencyfn and a damping ratio ξ, the ERS can also be calculated in both time domain and frequency domain as follows [14]:(11)ERSt=maxFfilterb,a,xa,(12)ERSf=2lnfn∗T+0.5772lnfn∗T∗∑f=flowf=fhiHf^2∗Gxxf∗df.


## 2.3. Derivation of Damage Equivalent Accelerated Test Spectrum

Since the accuracy of the damage equivalent accelerated test spectrum is strongly affected by uncertain parameters, dispersion of multiple non-Gaussian samples, and the smoothness of the derived PSD, an optimum design procedure is presented as follows:(1)
Calculate the FDS of multiple samples in time domain and take the envelope using different damping ratios (denoted by quality factorQ) and fatigue exponent b(2)
Predetermine the accelerated test time, and synthesize the initial PSD using the FDS envelope and Equation (10)(3)
ChooseQ and b in a dynamic way, so that the most conservative PSD is derived(4)
Calculate the FDS in the frequency domain using the derived PSD, Equation (7), and Equation (8)(5)
Compare the calculated FDS in Step (1) and Step (4), and define the weighted sum of FDS errors:(13)Werror=∑fi=flowestfhiestDnumfn−DenvfnDenvfn∗100∗Denvfn∑fi=flowestfhiestDenvfn,where flowest is the lowest natural frequency, fhiest is the highest natural frequency, Denvfn is the FDS envelope calculated in Step (1), and Dnumfn is the FDS calculated in Step (4)(6)
Update the PSD in Step (3) until a certain error is reached:(14)Gm+1fn=Gmfn∗DenvfnDnumfn2/b,where Gmfn is the damage equivalent PSD after m iterations(7)
Calculate the ERS in frequency domain using the updated PSD and Equation (12) and compare with the ERS envelope calculated in time domain(8)
Update the predetermined test time and PSD level to match both the FDS and ERS envelope of the field data


## 3. Results and Discussion

### 3.1. Case Study

Field measured non-Gaussian data are used to demonstrate the procedure in Section4. The field data is from an Ericsson Mast Project. The random acceleration signal induced by wind is taken from a test item installed on a mast, as shown in Figure 1. The test item is a standard Ericsson Micro Radio Base equipment, manufactured in cast lightweight alloy with integrated heat sinks, equipped with B&K WB0179 triaxial accelerometer sets and a wind speed measuring device. The dimensions are 530∗400∗185 mm, and the weight is about 21 kg. The back side of the box was mounted to the mast at a height of 50 meters.Figure 1
Test item.Acceleration signals were collected using 8 channels, as shown in Figure2. From Figure 2, we can see that the kurtosis of each signal is greater than 3, indicating the non-Gaussianity. To show the error introduced by direct calculation of PSD from the non-Gaussian signal, the data collected from channel 5 is used here as an example. The PSD is calculated, and the Gaussian signal is synthesized, as shown in Figure 3. The PSD and the PDF are compared between the field data and the synthesized Gaussian signal, as shown in Figures 4 and 5, respectively. From Figures 4 and 5, we can see that the PDF of the field data is clearly different from that of the synthesized Gaussian signal, although the PSD are basically the same. The pseudovelocity FDS of the measured non-Gaussian signal is significantly higher than that of the synthesized Gaussian load, as shown in Figure 6. From Figures 3–6, we can see that the direct calculation of PSD from a non-Gaussian signal will lead to an obvious undertest.Figure 2
Field measured non-Gaussian signal induced by wind.Figure 3
Comparison of time history between the measured non-Gaussian signal and the synthesized Gaussian signal.Figure 4
Comparison of PSD between the measured non-Gaussian signal and the synthesized Gaussian signal.Figure 5
Comparison of PDF between the measured non-Gaussian signal and the synthesized Gaussian signal.Figure 6
Comparison of FDS between the measured non-Gaussian signal and the synthesized Gaussian signal.To determine the effects of damping ratio and fatigue exponent on the calculation of response spectra and PSD synthesis, the pseudovelocity FDS of field data were calculated with different values ofQ (10, 25, 50) and b (4, 8, 12), as shown in Figure 7. From Figure 7, we can see that the FDS decreases and the dispersion increases, as the value of b increases. The value of Q has little effects on FDS (compared with the effect of b).Figure 7
FDS of measured wind load under differentQ and b values: (a) Q=10, b=4; (b) Q=25, b=4; (c) Q=50, b=4; (d) Q=10, b=8; (e) Q=25, b=8; (f) Q=50, b=8; (g) Q=10, b=12; (h) Q=25, b=12; (i) Q=50, b=12.
(a)(b)(c)(d)(e)(f)(g)(h)(i)Set the value ofQ to 50, calculate FDS with different values of b, and take the envelope. The equivalent PSD under different test times (T=600s,100s,10s) are derived using the FDS envelope and Equation (10), as shown in Figure 8. From Figure 8, we can see that the equivalent PSD level increases as the test time decreases. The effects of the value of b on the equivalent PSD depend on the predetermined accelerated test time. With T equals to 600 s (not accelerated), the equivalent PSD increases as the value of b increases. With T equals to 100 s, the value of b has little effect on the equivalent PSD. With T equals to 10 s, however, the equivalent PSD decreases as the value of b increases.Figure 8
Equivalent PSD under different test times with different values ofb: (a) test time 600 s; (b) test time 100 s; (c) test time 10s.
(a)(b)(c)The FDS of equivalent PSD for different test times remain the same, as shown in Figure9. A comparison is made between the FDS of equivalent PSD and non-Gaussian field data. Due to the nonsmooth shape of the equivalent PSD, the mismatch is clearly seen, especially within the peak and valley frequency band. To account for the nonsmooth shape, the PSD is updated using Equation (14). A better match between the FDS of the field data and the equivalent PSD is reached after the first iteration, as shown in Figure 10.Figure 9
Comparison of FDS between field data and equivalent PSD before correction.Figure 10
Comparison of FDS between field data and equivalent PSD before and after the first correction.The weighted sum of FDS error (see Equation (13)) calculated from the field data and equivalent PSD is reduced from 80.1% to 7.5% after 7 iterations. Particularly, the error is close to 0 within the peak and valley frequency band, as shown in Figures 11 and 12. The deviation of PSD between 150 Hz and 200 Hz is trivial due to the low level.Figure 11
FDS error analysis before and after correction.Figure 12
Comparison of FDS between field data and equivalent PSD before and after the 7th correction.A good match between the ERS calculated from the equivalent PSD and that calculated from the field data is shown in Figure13. Note that although the PSD level can be further increased to reduce the test time, the test item may fail due to the peak response which will not be experienced in the field.Figure 13
Comparison of ERS between equivalent PSD and field data.


## 3.1. Case Study

Field measured non-Gaussian data are used to demonstrate the procedure in Section4. The field data is from an Ericsson Mast Project. The random acceleration signal induced by wind is taken from a test item installed on a mast, as shown in Figure 1. The test item is a standard Ericsson Micro Radio Base equipment, manufactured in cast lightweight alloy with integrated heat sinks, equipped with B&K WB0179 triaxial accelerometer sets and a wind speed measuring device. The dimensions are 530∗400∗185 mm, and the weight is about 21 kg. The back side of the box was mounted to the mast at a height of 50 meters.Figure 1
Test item.Acceleration signals were collected using 8 channels, as shown in Figure2. From Figure 2, we can see that the kurtosis of each signal is greater than 3, indicating the non-Gaussianity. To show the error introduced by direct calculation of PSD from the non-Gaussian signal, the data collected from channel 5 is used here as an example. The PSD is calculated, and the Gaussian signal is synthesized, as shown in Figure 3. The PSD and the PDF are compared between the field data and the synthesized Gaussian signal, as shown in Figures 4 and 5, respectively. From Figures 4 and 5, we can see that the PDF of the field data is clearly different from that of the synthesized Gaussian signal, although the PSD are basically the same. The pseudovelocity FDS of the measured non-Gaussian signal is significantly higher than that of the synthesized Gaussian load, as shown in Figure 6. From Figures 3–6, we can see that the direct calculation of PSD from a non-Gaussian signal will lead to an obvious undertest.Figure 2
Field measured non-Gaussian signal induced by wind.Figure 3
Comparison of time history between the measured non-Gaussian signal and the synthesized Gaussian signal.Figure 4
Comparison of PSD between the measured non-Gaussian signal and the synthesized Gaussian signal.Figure 5
Comparison of PDF between the measured non-Gaussian signal and the synthesized Gaussian signal.Figure 6
Comparison of FDS between the measured non-Gaussian signal and the synthesized Gaussian signal.To determine the effects of damping ratio and fatigue exponent on the calculation of response spectra and PSD synthesis, the pseudovelocity FDS of field data were calculated with different values ofQ (10, 25, 50) and b (4, 8, 12), as shown in Figure 7. From Figure 7, we can see that the FDS decreases and the dispersion increases, as the value of b increases. The value of Q has little effects on FDS (compared with the effect of b).Figure 7
FDS of measured wind load under differentQ and b values: (a) Q=10, b=4; (b) Q=25, b=4; (c) Q=50, b=4; (d) Q=10, b=8; (e) Q=25, b=8; (f) Q=50, b=8; (g) Q=10, b=12; (h) Q=25, b=12; (i) Q=50, b=12.
(a)(b)(c)(d)(e)(f)(g)(h)(i)Set the value ofQ to 50, calculate FDS with different values of b, and take the envelope. The equivalent PSD under different test times (T=600s,100s,10s) are derived using the FDS envelope and Equation (10), as shown in Figure 8. From Figure 8, we can see that the equivalent PSD level increases as the test time decreases. The effects of the value of b on the equivalent PSD depend on the predetermined accelerated test time. With T equals to 600 s (not accelerated), the equivalent PSD increases as the value of b increases. With T equals to 100 s, the value of b has little effect on the equivalent PSD. With T equals to 10 s, however, the equivalent PSD decreases as the value of b increases.Figure 8
Equivalent PSD under different test times with different values ofb: (a) test time 600 s; (b) test time 100 s; (c) test time 10s.
(a)(b)(c)The FDS of equivalent PSD for different test times remain the same, as shown in Figure9. A comparison is made between the FDS of equivalent PSD and non-Gaussian field data. Due to the nonsmooth shape of the equivalent PSD, the mismatch is clearly seen, especially within the peak and valley frequency band. To account for the nonsmooth shape, the PSD is updated using Equation (14). A better match between the FDS of the field data and the equivalent PSD is reached after the first iteration, as shown in Figure 10.Figure 9
Comparison of FDS between field data and equivalent PSD before correction.Figure 10
Comparison of FDS between field data and equivalent PSD before and after the first correction.The weighted sum of FDS error (see Equation (13)) calculated from the field data and equivalent PSD is reduced from 80.1% to 7.5% after 7 iterations. Particularly, the error is close to 0 within the peak and valley frequency band, as shown in Figures 11 and 12. The deviation of PSD between 150 Hz and 200 Hz is trivial due to the low level.Figure 11
FDS error analysis before and after correction.Figure 12
Comparison of FDS between field data and equivalent PSD before and after the 7th correction.A good match between the ERS calculated from the equivalent PSD and that calculated from the field data is shown in Figure13. Note that although the PSD level can be further increased to reduce the test time, the test item may fail due to the peak response which will not be experienced in the field.Figure 13
Comparison of ERS between equivalent PSD and field data.


## 4. Conclusions

Considering the effects of dispersion of field data, parameter uncertainties, and the nonsmooth shape of the PSD, on the accuracy of the derived damage equivalent accelerated test spectrum, an optimized procedure for PSD synthesis using non-Gaussian data is presented in this paper. Conclusions are as follows:(1)
Direct calculation of PSD from non-Gaussian field data leads to undertest. The envelope of FDS and ERS can be used to derive the damage equivalent PSD, considering the dispersion of multiple non-Gaussian field data(2)
Model parameters must be selected dynamically if a priori information is unknown, i.e., the selection of model parameters depends on the predetermined accelerated test time, so as to avoid undertest(3)
The damage equivalent PSD must be updated to account for its nonsmooth shape. The weighted sum of FDS error between field data and derived PSD can be reduced dramatically after only a few iterations. Particularly, the error can be reduced close to 0 within the peak and valley frequency band



---

*Source: 3668726-2021-09-06.xml*