Next Article in Journal
A Novel Fast Parallel Batch Scheduling Algorithm for Solving the Independent Job Problem
Next Article in Special Issue
Learning Collision Situation to Convolutional Neural Network Using Collision Grid Map Based on Probability Scheme
Previous Article in Journal
The Effect of Shape on Chloride Penetration of Circular Reinforcement Concrete Columns and Its Durability Design
Previous Article in Special Issue
Structural Colors Based on Amorphous Arrays Comprised Solely of Silica Particles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simplified Vibration PSD Synthesis Method for MIL-STD-810

Department of Reliability Assessment, Korea Institute of Machinery and Materials, Daejeon 34103, Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(2), 458; https://doi.org/10.3390/app10020458
Submission received: 30 November 2019 / Revised: 22 December 2019 / Accepted: 6 January 2020 / Published: 8 January 2020
(This article belongs to the Special Issue Selected Papers from the ICMR 2019)

Abstract

:

Featured Application

In MIL-STD-810, a DP (damage potential) formula, a concise version of FDS (fatigue damage spectrum) formula, is used to synthesize vibration conditions. In this paper, a more simplified formula that yields the same synthesis results as the DP formula is proposed.

Abstract

In MIL-STD-810, the environmental engineering consideration and laboratory tests by US military, the DP (damage potential) formula is used to calculate the fatigue damage by vibration environments which is a simplified version of a FDS (fatigue damage spectrum) formula. DP, however, was originally made for comparison between different test standards and is not an optimized formula for the synthesis of vibration environments. This paper presents the GRS (Gaussian Random Synthesis) formula for only vibration test synthesis featuring simplified DP formula and produce the same result as DP. Although the GRS provides insight into the synthesis of vibration tests, it must be used with care because it inherits the constraints of the original FDS formula.

1. Introduction

Products must perform adequately under all environmental conditions associated with its service life. Laboratory environmental tests are used extensively instead of more time-consuming and less cost effective field exposure tests. The MIL-STD-810, Environmental Engineering Considerations and Laboratory Tests, establishes laboratory test methods that replicate the effects of environments on the equipment rather than imitating the environments themselves. Vibration test among dozens of environmental test methods are also performed to verify that materiel will function in and withstand the vibration exposures of a life cycle. The most dynamically accurate method to reproduce the full exposure would be to sequentially vibrate the product to all the individual, uncompressed events representing its full lifecycle. However, such an approach is generally not feasible from both schedule and economic perspectives, and some compromises must be made to realize the benefits of testing in the laboratory. Time compression techniques based on fatigue equivalency are typically employed such that vibration testing can be performed in a timely and economic manner. The objective of the laboratory vibration test is to generate the equivalent fatigue damage accumulated in the products by vibration throughout their lifetime [1].
This fatigue damage by vibration can be quantified by the fatigue damage spectrum (FDS) methods which is a spectral representation of fatigue damage as a function of the SDOF system’s natural frequency [2]. The FDS can be computed in time domain [3] and in frequency domain [4,5,6]. In frequency domain, the formula calculating FDS can be obtained based on some assumptions. In case of a low damping SDOF system and Gaussian random input, in which the peak distribution of the response becomes a Rayleigh one, the FDS equation gets much simpler [7]. By simplifying this further, Henderson and Piersol proposed a descriptor for comparing random test conditions and named it as damage potential (DP) [7]. MIL-STD-810 uses this damage potential formula for synthesizing random vibration conditions. The purpose of DP is not to calculate the accurate fatigue damage, but to compare the relative fatigue damage between test specifications. However, as the synthesized vibration condition using DP is the same as using FDS, it is more efficient to use DP which is a much simpler formula. However, DP is not an optimized index for vibration condition synthesis. In this paper, we propose the slimmest formula in which unnecessary factors are removed from the DP to suit the purpose of vibration condition synthesis. This simplified equation, called GRS (Gaussian Random Synthesis), produce the same results as DP in vibration synthesis application only.
To introduce GRS, we first examine the concept of FDS thoroughly in time domain and how the FDS formula, which is the basis of DP, was derived in frequency domain, and then look at how to synthesize random vibration conditions under Miner’s cumulative damage law in MIL-STD-810. Then the more simplified formula, GRS that is newly presented in this paper is derived. Actually, it reveals a concise expression in the process of composing the formula of itself. Finally, we present an example of synthesis using FDS, DP, and the new GRS equations and discuss the results.
GRS formula aims just to present a more convenient method for those who wish to create vibration test conditions by the method introduced in MIL-STD-810. Therefore, the constraints and problems that the method has inherited must be checked as well.

2. Fatigue Damage Spectrum (FDS)

2.1. FDS in General

Fatigue Damage Spectrum represents the cumulative damage due to the random responses of a series of Single Degree of Freedom (SDOF) systems, with varying natural frequency, excited by random acceleration, both in the time and frequency domain [5]. So, the calculation of FDS is based on the response of SDOF system, which is the virtual specimen, to the input vibration. It has been shown that the velocity response is proportional to stress [7,8]. Once a stress time history is given, a stress histogram can be derived by the various counting methods and then damage index can be calculated for each stress level in case the slope of S-N (stress–number) curve in log–log plot is given and finally the total damage can be simply summed thanks to the cumulative damage law. The FDS of the frequency range of interest is derived by changing the resonant frequency of the virtual specimen and repeating the entire process as shown in Figure 1.
The concept of FDS was explained in time domain above, but a lot of researchers calculate FDS in frequency domain. The paper cited in many FDS literature is that of Miles [9]. This is an article on the structural fatigue of aircraft, which describes the basic concept well, but the input is not vibrations but aerodynamic forces. The most extensive book on vibration fatigue is the Lalanne’s five-volume books [5], which presents generalized formulas for accurate calculations and simplified equation based on some assumptions. McNeil described the application of FDS based on the pseudo-velocity response of the SDOF systems [3]. Halfpenny et al. discussed FDS to represent the cumulative fatigue damage caused by relative displacement response not velocity or pseudo velocity [4]. MIL-STD-810 describes the synthesis process of vibration conditions in detail from ‘G w/Change1’ version published in 2014 [1]. The DP (Damage Potential) formula used in the standard is from the paper by Henderson and Piersol [7]. They simplified the FDS formula and proposed the DP formula to compare the power spectral density and test duration of different excitation test specifications in terms of fatigue damage [7]. In case of the NATO environmental test specification, AECTP-200 (2006) [10] uses Lalanne’s formula and AECTP-240 (2009) [2] uses Henderson and Piersol’s Damage Potential.

2.2. FDS in MIL-STD-810

There are various formulas calculating FDS, but here we thoroughly review the derivation process of the DP formula. As a special case, if the input is a random waveform featuring (1) Gaussian distribution and (2) constant PSD, and also (3) the damping coefficient of the SDOF system is low (ζ < 0.1), then the response of the system is also normally distributed random waveform and the peak distribution shows a Rayleigh one [11] as depicted in Figure 2.
This peak distribution corresponds to the histogram e in Figure 1, and is represented by the following equation.
n ( s ) = f n T s σ 2 exp ( s 2 2 σ 2 ) = f n T · p ( s )
p(s) is the probability density function of the Rayleigh distribution, where σ is the scale parameter, which is equal to the standard deviation from the response waveform. Assuming an SDOF system with a low damping coefficient, the response waveform oscillates based on its natural frequency fn, so the total number of peaks during time T can be approximated by the product of fn and T. Multiplying the Rayleigh function by the total number of peaks yields the distribution function of peaks over stress, n(s). The standard deviation, σ, of the normally distributed output random waveform is squared root of the area of the stress response PSD, Gs(f), which can be obtained by multiplying the input PSD with the squared amplitude of the system’s frequency response function (FRF), H(f), in the frequency domain as shown in Equation (2). G is set to be a constant to enable the integration and k is conversion factor from velocity to stress. H is the FRF of pseudo velocity output to the acceleration input. These FRF relationships are well documented in ISO 18431-4 [12].
σ 2 = 0 G s ( f ) d f = 0 k 2 G · | H ( f ) | 2 d f = k 2 G 0 | H ( f ) | 2 d f = ( k 2 π f n ) 2 G π f n 4 ς
| H ( f ) | 2 = | p v a | 2 = ( 1 / 2 π f n ) 2 ( 1 ( f / f n ) 2 ) 2 + ( 2 ζ f / f n ) 2
Since n(s) has been completely defined, we calculate the damage index by stress and add them all up to get the cumulative damage index D with Equation (4). This corresponds to step f in Figure 1. Here, the Basquin’s Equation (5) is used to represent N(s) which is the number to failure.
D ( f n ) = 0 n ( s ) N ( s ) d s = 0 f n T p ( s ) C s b d s = f n T C 0 p ( s ) s b d s
N ( s ) = C s b
The last integral term in Equation (4) is the b-th moment of the Rayleigh function, p(s), and can be expressed using the Gamma function [13] as below.
0 p ( s ) s b d s = σ 2 2 b / 2 Γ [ 1 + b / 2 ]
Applying the results of Equations (2) and (6) to Equation (4) yields Equation (7).
D ( f n ) = f n T C [ G ( f n ) k 2 8 π f n ζ ] b / 2 Γ [ 1 + b / 2 ]
The equation above is a formula for calculating the fatigue damage index accumulated in a lightly damped SDOF system excited by Gaussian random vibration, and it can be called the FDS equation.
The constant terms are collected in Equation (7) as follows,
D ( f n ) = k b C ( 8 π ) b / 2 Γ [ 1 + b / 2 ] f n T [ G ( f n ) f n ζ ] b / 2 = c o n s t f n T [ G ( f n ) f n ζ ] b / 2
DP is extracted from Equation (8) and defined as follows.
D P ( f n ) = f n T [ G ( f n ) f n ζ ] b / 2
Since the DP evolves from a proportional relationship with fatigue damage index, the acceleration PSD can be in the more common units of g2/Hz. And, the value is not bounded by unity, as is the actual damage D estimated by Equation (7) [7].

3. Vibration PSD Synthesis

3.1. Vibration PSD Synthesis in MIL-STD-810

MIL-STD-810 states that since fatigue damage is based on a cumulative effect of various environments or conditions, a cumulative fatigue damage index can be calculated as the sum of the fatigue damage spectra for individual environments [1]. Thus
D P t o t a l ( f n ) = D P 1 + D P 2 + + D P N = i = 1 N D P i ( f n )
The individual damage indexes, D P i , obtained by Equation (9) are cumulated by Equation (10) to generate the total damage index D P t o t a l .
Then the PSD is derived by Equation (11) which is the inverse form of Equation (9).
G t e s t ( f n ) = f n ζ [ D P t o t a l ( f n ) f n T t e s t ] 2 / b
Equation (11) reads that in order to apply damage, DPtotal, in Ttest seconds to the virtual specimen, the PSD must be Gtest.
In addition to Equation (9), the damage index can be calculated by various methods as discussed in Section 2. The important thing is that the formula calculating FDS and the inverse formula calculating PSD must be consistent. If the Equation (7) is used instead to calculate FDS, the PSD for testing shall be obtained using the following formula. The ‘const’ is the same as used in Equation (8).
G t e s t ( f n ) = 8 π f n ζ k 2 [ D t o t a l ( f n ) · C f n T t e s t · Γ [ 1 + b / 2 ] ] b / 2 = f n ζ [ D ( f n ) c o n s t · f n T ] 2 / b
In order to create a test specification, in summary, by synthesizing vibration environment conditions according to the MIL-STD-810, calculate each damage index with the Equation (9) and add them to calculate the cumulative damage index, and then substitute it into the Equation (11) along with the test time T. T and G will be adjusted to prevent the G from becoming too large or the test time from being unrealistically extended.

3.2. Vibration PSD Synthesis Made Simpler

Originally, the concept of damage potential (DP) was proposed to describe the potential damage delivered to the virtual test items, which are the SDOF systems with different natural frequencies, by random vibration test [7]. So, DP itself is meaningful for comparison purpose. In the combination process, however, DP is an intermediate process, not a final output. Therefore, the Equation (9) that produce the individual DPs and the Equation (10) that combines them is substituted into the Equation (11) that calculates the test PSD, which can be summarized as follows.
G t e s t ( f n ) = f n ζ [ D P t o t a l ( f n ) f n T t e s t ] 2 b = f n ζ [ i = 1 N T i T t e s t ( G i ( f n ) f n ζ ) b / 2 ] 2 / b = ( i = 1 N T i ( G i ( f n ) ) b / 2 T t e s t ) 2 / b
This formula can be rewritten as follows to look simple,
for   each   f n , G t e s t = ( i = 1 N T i G i b / 2 T t e s t ) 2 / b = ( T 1 G 1 b / 2 + T 2 G 2 b / 2 + + T N G N b / 2 T t e s t ) 2 / b
This can be rearranged as follows in a symmetric form which shows the ultimately simple relationship between field conditions and test condition in fatigue damage perspective when input PSDs are Gaussian waveform.
T t e s t G t e s t b / 2 = T 1 G 1 b / 2 + T 2 G 2 b / 2 + + T N G N b / 2
Equations (15) shows that we do not need natural frequency and damping ratio in combining the PSDs. Now we can remove the unnecessary parameters and define a new index for only synthesizing process as follows and named it as Gaussian random synthesis (GRS) to emphasize its inherent constraints and usage at the same time,
G R S ( f n ) = T · G ( f n ) b / 2
The inverse form is
G ( f n ) = ( G R S ( f n ) T ) 2 / b
The same results as Equations (14) and (15) are obtained from other FDS equation. The nCode, fatigue and durability analysis software, uses the FDS formula under the assumption that the relative displacement response instead of the velocity response of the virtual specimen is proportional to stress [4]. The equation below shows the process that the constant terms are collected in the FDS as in Equation (8) when DP is obtained.
D ( f n ) = f n T k b C [ Q · G ( f n ) 2 ( 2 π f n ) 3 ] b / 2 Γ [ 1 + b 2 ] = c o n s t · f n T [ Q · G ( f n ) f n 3 ] b / 2
Here, Q = 1 / 2 ζ , The inverse of Equation (18) calculating the test PSD is
G t e s t ( f n ) = f n 3 Q [ D P t o t a l ( f n ) c o n s t · f n T t e s t ] 2 b = f n 3 Q [ i = 1 N T i T t e s t ( Q · G i ( f n ) f n 3 ) 2 b ] 2 b = ( i = 1 N T i ( G i ( f n ) ) 2 b T t e s t ) 2 b
The last expression is same as Equation (14) again which means we would get the same synthesized PSD using both FDS formula based on different FRF.
Table 1 is a simple example synthesizing two PSDs by using FDS based on pseudo velocity, DP, and GRS. The frequency is fixed at 300 Hz, and b = 4, ζ = 0.05, C = 1, and k = 1 is given then the intermediate calculation results are different but the final test PSDs to generate the equivalent fatigue damage index for 10 min test time are exactly the same, 30 (m/s2)2/Hz.
It should be noted that the intermediate calculation values have no physical meanings and that the unit of the test PSD follows the unit of the input PSDs. Equation (15) is the best representation showing these characteristics in synthesis.

4. Conclusions

In this paper, the DP formula used to synthesize vibration test conditions in MIL-STD-810 was further simplified and GRS formula which is the most concise form for vibration synthesis is proposed. GRS is used under the assumption that vibration conditions are Gaussian random, constant PSD, and the response system has low damping ratio (less than 0.1).
DP is presented as a descriptor for comparing random test conditions in their fatigue damage potential perspective, which MIL-STD-810 uses to synthesize a Gaussian random vibration environment. The results of combining PSDs with DP are no different from the results by DP’s original FDS formula. But in the same context, the same synthesis results can be obtained with the simplified formula, GRS. The numerical verification was done through simple example.
However, both DP and GRS methods are subject to the same constraints as of the original FDS as mentioned. The assumption that the damping coefficient is less than 0.1 when synthesizing vibration using DP or GRS in MIL-STD-810 is valid if the material of the component considering fatigue damage is metal. However, the effectiveness of this method should be reconsidered for synthetic resins with high damping coefficients for example. In addition, it is not valid when the input vibration includes a lot of shocks and is not Gaussian distribution. Otherwise, if the DP can be used, it is efficient to use the simplified GRS formula instead for synthesis of random vibration condition. Lastly, it is emphasized that it would be helpful to have an intuitive understanding of the relationship between PSDs through Equation (15).

Author Contributions

Formal analysis and writing, J.J.; Project administration, review & editing, J.-W.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by “Virtual Engineering Platform Construction for Future Transportation Vehicle Components Manufacturing Project” supported by Ministry of Trade, Industry and Energy.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

n(s)peak distributionaacceleration
fnnatural frequencyN(s)Number of stress cycles to failure
Tdurationbfatigue parameter
sstressCconstant of proportionality
σstandard deviationΓgamma function
p(s)Rayleigh distribution functionDcumulative damage index
GsPSD of stressDPdamage potential
GPSD of input accelerationTi, Giexposure time and PSD of event i
kproportional constantTtestexposure time for test
Hfrequency response functionGtestPSD for test
ζdamping ratioQamplification factor
pvpseudo-velocityGRSGaussian random synthesis

References

  1. MIL-STD-810G w/Change 1, Environmental Engineering Considerations and Laboratory Tests; U.S. Department of Defense: Washington, DC, USA, 2014.
  2. AECTP-240, EDITION 1, Allied Environmental Conditions and Test Publication, Mechanical Conditions; NATO Standardization Agency: Brussels, Belgium, 2009.
  3. McNeill, S. Implementing the Fatigue Damage Spectrum and Fatigue Damage Equivalent Vibration Testing. In Proceedings of the 79th Shock and Vibration Symposium, Orlando, FL, USA, 26–30 October 2008. [Google Scholar]
  4. Halfpenny, A. Accelerated Vibration Testing Based on Fatigue Damage Spectra; nCode International: Sheffield, UK, 2006. [Google Scholar]
  5. Lalanne, C. Mechanical Vibration and Shock, Random Vibration, 2nd ed.; John Willy & Sons, Inc.: Hoboken, NJ, USA, 2009. [Google Scholar]
  6. Plaskitt, R.; Halfpenny, A. The Fatigue Damage Spectrum: Explained. In Proceedings of the 28th ICAF Symposium (ICAF 2015), Helsinki, Finland, 1–5 June 2015; pp. 3–5. [Google Scholar]
  7. Henderson, G.R.; Piersol, A.G. Fatigue Damage Related Descriptor for Random Vibration Test Environments. Sound Vib. 1995, 29, 20–24. [Google Scholar]
  8. Hunt, F.V. Stress and Strain Limits on the Attainable Velocity in Mechanical Vibration. J. Acoust. Soc. Am. 1960, 32, 1123–1128. [Google Scholar] [CrossRef]
  9. Miles, J.W. On structural fatigue under random loading. J. Aeronaut. Sci. 1954, 21, 753–762. [Google Scholar] [CrossRef]
  10. AECTP-200-3, Environmental Conditions; NATO Standardization Agency: Brussels, Belgium, 2006.
  11. Cartwright, D.E.; Longuet-Higgins, M.S. The Statistical Distribution of the Maxima of a Random Function. Proc. R. Soc. Lond. Ser. A Math. Phys. Sci. 1956, 237, 212–232. [Google Scholar]
  12. ISO 18431-4, Mechanical Vibration and Shock—Signal Processing—Part4: Shock-Response Spectrum Analysis; International Standard Organization: London, UK, 2007.
  13. Davis, P.J. Leonard Euler’s Integral: A Historical Profile of the Gamma Function; National Bureau of Standard: Washington, DC, USA, 1959. [Google Scholar]
Figure 1. Flowchart of fatigue damage spectrum (FDS) calculations, (a) input acceleration, (b) Single Degree of Freedom (SDOF) system which is the virtual specimen, (c) velocity response, (d) stress waveform which is obtained by multiplying k to the velocity response, (e) stress counting and idealized stress–number (S-N) curve, (f) accumulation of damage index, (g) FDS.
Figure 1. Flowchart of fatigue damage spectrum (FDS) calculations, (a) input acceleration, (b) Single Degree of Freedom (SDOF) system which is the virtual specimen, (c) velocity response, (d) stress waveform which is obtained by multiplying k to the velocity response, (e) stress counting and idealized stress–number (S-N) curve, (f) accumulation of damage index, (g) FDS.
Applsci 10 00458 g001
Figure 2. Velocity response of the SDOF system to Gaussian random acceleration input with constant power spectral density (PSD) shows Gaussian distribution. The peak distribution of the response becomes a Rayleigh function. The standard deviation of the response distribution matches the scale parameter of the Rayleigh function. a(t) is input acceleration and s(t) is stress waveform (above). Stress PSD is calculated by multiplying input PSD and squared stress amplitude of system’s frequency response function (below).
Figure 2. Velocity response of the SDOF system to Gaussian random acceleration input with constant power spectral density (PSD) shows Gaussian distribution. The peak distribution of the response becomes a Rayleigh function. The standard deviation of the response distribution matches the scale parameter of the Rayleigh function. a(t) is input acceleration and s(t) is stress waveform (above). Stress PSD is calculated by multiplying input PSD and squared stress amplitude of system’s frequency response function (below).
Applsci 10 00458 g002
Table 1. PSD synthesis with FDS, damage potential (DP), and GRS 1.
Table 1. PSD synthesis with FDS, damage potential (DP), and GRS 1.
EventsCalculations
EventFrequencyTimePSDFDS Equation (7)DP Equation (9)GRS Equation (16)
1300 Hz2 h5 (m/s2) 2/Hz7.6 × 1022.4 × 1051.8 × 105
2300 Hz1 h10 (m/s2) 2/Hz1.5 × 1039.6 × 1053.6 × 105
SynthesisSUM (Event 1 + 2)2.3 × 10312.0 × 1055.4 × 105
test PSD
for test T = 10 min
Equation (12)Equation (11)Equation (17)
30 (m/s2) 2/Hz30 (m/s2) 2/Hz30 (m/s2) 2/Hz
1 b = 4, ζ = 0.05.

Share and Cite

MDPI and ACS Style

Jang, J.; Park, J.-W. Simplified Vibration PSD Synthesis Method for MIL-STD-810. Appl. Sci. 2020, 10, 458. https://doi.org/10.3390/app10020458

AMA Style

Jang J, Park J-W. Simplified Vibration PSD Synthesis Method for MIL-STD-810. Applied Sciences. 2020; 10(2):458. https://doi.org/10.3390/app10020458

Chicago/Turabian Style

Jang, Jin, and Jong-Won Park. 2020. "Simplified Vibration PSD Synthesis Method for MIL-STD-810" Applied Sciences 10, no. 2: 458. https://doi.org/10.3390/app10020458

APA Style

Jang, J., & Park, J. -W. (2020). Simplified Vibration PSD Synthesis Method for MIL-STD-810. Applied Sciences, 10(2), 458. https://doi.org/10.3390/app10020458

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop