Issue
Mechanics & Industry
Volume 18, Number 8, 2017
Experimental Vibration Analysis
Article Number 809
Number of page(s) 12
DOI https://doi.org/10.1051/meca/2017057
Published online 05 September 2018

© AFM, EDP Sciences 2018

1 Introduction

In the context of the design of structures to vibration fatigue, it is common to characterize the stationary random excitation processes by their acceleration power spectral density (PSD), within a range of relevant frequencies, generally covering the LF [1–50 Hz] and MF domains [50–500 Hz], for the case of wheeled vehicles.

Then relying on the conventional assumptions of Gaussian signals, characterization methods of the average fatigue damage generated by these random acceleration processes may be used. This is the case of characterization techniques by the response spectra (RS) method, widely used by manufacturers in the civil and military sectors, as part of the development of design vibration specifications.

The RS technique consists in characterizing the severity of the excitation process through a simple mechanical model consisting of a linear system with one degree of freedom (dof) whose natural frequency is varied f0 within a sufficiently wide domain covering the natural frequencies of the structure to size, and limited by the traditional criterion of the effective masses. In many industrial applications, the criterion of effective masses consists in retaining all the modes in which the sum of the effective masses is 90% of the total mass of the structure to size [1]. At this stage, it is a generalization of the Biot model, historically used in the calculation of the shock response spectrum (SRS) for transient excitations and that is adapted to the case of random vibration [2] referred to as the fatigue damage spectrum (FDS). Thus, for a process of random vibration evaluated generally on a limited portion Tref of the total stress time Tt, it is possible to characterize its severity by applying either a spectral domain approach (stochastic approach) or a temporal domain approach (deterministic approach), each of which has disadvantages in the case of non-Gaussian random processes, as discussed below.

2 Selected standard system

When a vibrational excitation process x¨(t) is applied to the base of a mechanical system with one dof (mass − spring − damper) the standard system responds by generating a relative movement z(t) between the mass and the base of the standard system (see Fig. 1). This relative movement is then dependent on the natural frequency f0 of the standard and is proportional to the uniaxial stress state σ(t) generated in the structure, provided that the latter is considered linear (an assumption commonly used in design). The proportionality constant K (Stress − Relative displacement) is imposed by the designer and is kept constant regardless of the value of the natural frequency f0.

Therefore we obtain: σ(t)=K.z(t),(1) where σ(t): amplitude of the vibration stress (N/m2), K: proportionality constant (N/m3) between stress and relative displacement (K is a constant if it is assumed that the structure to be designed is linear), z(t): amplitude of the relative displacement (m).

In which z(t)=y(t)x(t),(2) where y(t): absolute displacement of the mass m (m), x(t): absolute displacement of the mobile base (m).

The standard system for mass m (in Kg), stiffness k (in N/m) and viscous damping c (in N.s/m) is then characterized by its modal parameters, namely: ξ=c2.k.mandQ=1/2.ξ.(3) and f0=12.π.kmandω0=2.π.f0,(4)where ξ: critical or modal damping,Q: modal overtension factor, f0: natural or modal frequency (Hz), ω0: natural or modal pulsation (rd/s).

The uniaxial stress state of the standard system is characterized through the standard system by assessing its response, in terms of relative displacement which is estimated either in frequency or temporal form, depending on whether a frequency or temporal domain approach is used to evaluate the fatigue damage associated with standard system in question.

thumbnail Fig. 1

Standard system with one dof.

2.1 Frequency domain approach

In the case of a frequency domain approach, the standard system is characterized by its transfer function Hx¨z(ω=2.π.f), defined between the acceleration of x¨(t) the vibrational excitation and the relative displacement z(t) of the response, which is then characterized by its noted unilateral Power Spectral Density PSDzz(ω = 2.π . f).

The expression of the standard system transfer function is then defined as follows: H x¨z(ω)=Z(ω) X¨(ω)=|H x¨z(ω)|.ej.H x¨z(ω),(5) where Z(ω)=TF[z(t)]andX¨(ω)=TF[x¨(t)]|Hx¨z(ω)|: modulusoftheTransfertFunction Hx¨z(ω):phaseoftheTransfertFunction TF[.]:       FourierTransformoperator. And:|Hx¨z(ω)|=1ω02.1(1h2)2+4.ξ2.h2 whereh=ωω0=ff0,ω0=km=2.π.f0andQ=1/2.ξ.

Due to these technical considerations, it is possible to simply assess the pace of unilateral PSD of the relative movement PSDzz(ω) in m2/Hz, knowing the unilateral PSD of the excitation PSDx¨x¨(ω) in (m/s2)2/Hz, based on the classic expressions of the dynamics of linear structures [3] namely:     PSDzz(ω)=|H x¨z(ω)|2.PSD x¨ x¨(ω) whereω[ωmin,ωmax].(8)

2.2 Temporal domain approach

In the case of a temporal domain approach, the standard system is characterized by its impulse function hx¨z(t), defined between the acceleration x¨(t) of the vibrational excitation and the relative displacement z(t) of the response, which is then characterized by its z-transform,Hx¨z(z) within the meaning of the first order hold (FOH) method, used for the mathematical reconstruction of the response z(t) sampled [4]. The expression of the z-transform of the impulse function of the standard system is then defined as follows [5]: Hx¨z(z)=(z1)2Δt.z.TZ {TL1[Hx¨z(p)p²]}=b0+b1.z1+b2.z21+a1.z1+a2.z2 wherez=exp(p.Δt)andFe=1Δt, where Fe: Sampling frequency of the random excitation process x¨(t) and the response z(t) in Hz.

The advantage of using the z-transform instead of the impulse function basically stems from in the speed of the calculations of the z(t) sampled response, which, as part of the process of characterizing the random non-Gaussian stress processes, is an important issue, as we shall see below with the notion of temporal XFS, explained below in Section 5.2.

Due to these technical considerations, it is then possible to simply evaluate the expression of the sampled response of the relative displacement z(i . Δt) = z(i) over the reference duration Tref corresponding to the acquisition time of the excitation processes x¨(t) to be characterized in terms of vibration fatigue, namely: z(i)=a1.z(i1)a2.z(i2)+b0.x¨(i)+b1.x¨(i1)+b2.x¨(i2) wherei=[1,n]andn=Tréf.Fe. where Tref: acquisition time of the excitation processes x¨(t)

And: a1=2.exp(ξ.ω0.Δt).cos(ωd.Δt)anda2=exp(2.ξ.ω0.Δt)whereωd=ω0.1ξ2andQ=1/(2.ξ)andFe=1/Δt b0=1ω03.Δt.{2.ξ.[exp(ξ.ω0.Δt).cos(ωd.Δt)1]+exp(ξ.ω0.Δt).[ω0ωd.(2.ξ21).sin(ωd.Δt)]+ω0.Δt} b1=1ω03.Δt.{2.ω0.Δt.exp(ξ.ω0.Δt).cos(ωd.Δt)2.ω0ωd.(2.ξ21).exp(ξ.ω0.Δt).sin(ωd.Δt)+2.ξ.[1exp(2.ξ.ω0.Δt)]} b2=1ω03.Δt.{(2.ξ+ω0.Δt).exp(2.ξ.ω0.Δt)+exp(ξ.ω0.Δt).[ω0ωd.(2.ξ21).sin(ωd.Δt)2.ξ.cos(ωd.Δt)]}.

3 Approaches for calculating FDS

As presented in the introduction, there are two possible approaches for calculating the FDS of a stationary random excitation process, depending on whether we wish:

  • to adopt the case of a stochastic approach to estimate the statistics of damage fatigue levels generated by the excitation process x¨(t) that we wish to characterize favorably in terms of mean and (1−α) quantile value, or order quantile (1−α), for the total stress time during which Tt the random excitation process is applied. This stochastic approach is especially important in that it can be used to extrapolate statistically fatigue damage levels over a Tt longer (Tt > Tref) period than that associated with the reference time Tref for which the excitation process x¨(t) is measured or estimated by modeling (as shown in Sect. 3.1)

  • or adopt the case of a deterministic approach allowing only calculate the damage fatigue levels generated by the excitation process x¨(t) acquired on the reference time Tref by measurement or numerical modeling, without being able to extrapolate favorably on the total duration Tt other than by a deterministic approach to proportionality, as presented in Section 3.2. Therefore, it is clear that this method of temporal calculation can be non-conservative and can therefore lead to significant mechanical undersizing for the structures being studied. This is especially true if the duration Tref is short compared with the total duration Tt of the stress in question.

Beyond the identified risk of undersizing in the case of temporal domain approaches, the latter lead to prohibitive computing time, unlike the spectral methods. Indeed, the advantage of spectral methods lies in the fact that the calculation of the FDS is analytical, but the associated up-crossing risk α should be made clear. For example, if one takes for simplicity the case of random stationary and Gaussian processes, it is possible to show that this notion of FDS is implicitly associated with an up-crossing risk which is not constant and equal to 50%, as suggested by many designers and specifiers. Indeed the work of the author in 2007 [6] showed that the level of accumulated damage could favorably be described by a 2-parameter Weibull distribution (β˜,σ˜) and not by a normal distribution conventionally used by some authors. Due to these technical considerations, it is possible to show that the concept of spectral FDS is associated with an up-crossing risk, depending on the shape parameter β˜ of this Weibull distribution, which is then based on the calculation parameters of the Spectral FDS:

  • f0 and Q: modal parameters of the standard system;

  • Tt: total duration of the vibratory stress x¨(t) in question;

  • b: basquin slope of the fatigue behavior of the structure in question, described by equation (12) and the following Figure 2.

N.[σa]b=C=Ab,(12) where b and C = Ab: material parameters of the structure in question, A: alternating failure stress of the material (N/m2), N: average number of failure cycles (cycles), σa: sinusoidal alternating stress, corresponding to local extrema (maxima and/or minima) of the stress process σ(t) = K . z(t) defined by (1) in N/m2.

The fluctuation of the exceedance threshold associated with the conventional expression of the spectral FDS, depending on the frequency f0, is not very reassuring as a factor for designers who are in charge of complex, high-performance structures. It is also for these reasons that the model of the conventional spectral FDS, widely used for establishing vibration specifications [7], has been generalized by the spectral XFS approach, which is the counterpart of the spectral XRS approach also proposed in 2007 by the author [6].

thumbnail Fig. 2

Basquin model and convention of writing.

3.1 Calculation of the spectral FDS

There are different possible approaches to the calculation of spectral FDS according to the counting method (Peak-Valley and/or Rainflow Counting [8]). In the present case, the one proposed here is based on the extrema counting method (Peak-Valley) and is associated with the following factors:

  • the excitation process x¨(t) of duration Tt is stationary, ergodic, Gaussian and has a zero mean (all the spectral moments are therefore constant and calculable from the PSD of the relative movement of the standard system);

  • the standard system adopted is linear, the relative displacement z(t) of the model is itself stationary, ergodic, Gaussian and zero mean. In addition the standard system being a quality factor Q ≥ 3, the random stress process has a narrow band centered on f0 in which the local maxima are characterized by a Rayleigh distribution;

  • the fatigue behavior of the structure is described by its Basquin's law, defined by its parameters b and C = Ab (see relation 12 above);

  • the accumulated damage per cycle is based on Miner's conventional linear accumulation rule [9];

  • only the positive local maxima and negative local minima of the stress process σ(t) = K . z(t) are assumed to be damaging.

The spectral FDS is then calculated from the unilateral PSD of the random stress process σ(t) which is assumed to be Gaussian and therefore from the random process associated with the relative displacement z(t) of the standard system, insofar as the latter is considered here to be linear. For each natural frequency f0, the value of the average cumulative damage E [Dc(ξ)] is calculated from the damage produced by different cycles extracts over the total duration of Tt the excitation process x¨(t) in question and accumulated linearly through the fatigue behavior of the structure in question. Based on the above assumptions, it can be easily shown that the analytical expression of the spectral FDS of a Gaussian random process x¨(t) can be calculated simply by using the Rayleigh distribution fZmax(z) associated with the random process of the local maxima of the relative displacement z(t), i.e.: fZmax(z)=zm0.exp(z22.m0),wherez[0,+[ andm0=zeff2(f0)=f=0+PSDzz(f).dfPSDzz(f):unilateralPSDof the random processz(t).

Hence the conservative expression of the (narrow-band) spectral FDS is: FDSspBE(f0) FDSspBE(f0)={CKb.E[Dc(ξ)]=n.z=0+zb.fZmax(z).dz(f0.Tt).[2.zeff2(f0)]b/2.Γ(1+b2)wheren=f0.Tt.

This simple expression of the narrow-band spectral FDS is conservative and can be generalized to the case of broad-band Gaussian processes, relying on the Rice Theory of stationary random Gaussian processes [10,11]. However, the disadvantage of these spectral methods is that they remain subject to a restrictive assumption of gaussianity, which now needs to be abandoned, in order to positively expand the scope of the RS method.

3.2 Calculation of the temporal (deterministic) FDS

In the case of the temporal FDS, the relative displacement z(t) is obtained by the FOH method presented in section 2.2 above and this for t ∈ [0, Tref], as the excitation x¨(t) is known only for the reference period Tref. The absolute values of the peaks (positive maximum) and valleys (negative minimum) of the relative displacement z(t) evaluated over time Tref, are determined using the conventional counting method of extrema which associates with a maximum or minimum a half cycle of damage to stay in the same calculation assumptions as the spectral FDS. This “Peak-Valley” method can then be used to retrieve from the relative displacement z(t) of duration Tref:

  • the nmax damage half-cycles associated with the positive maximum amplitude Zmax (i) for i=1anmax;

  • the nmin damage half-cycles associated with the negative minimum amplitude Zmin(i)  for i=1anmin.

The half cycles thus retrieved from the random stress process σ(t) = K . z(t), defined on the reference duration Tref are then accumulated linearly in terms of partial damage Dp using Basquin's law defined by (12) and Miner's rule defined as follows: Dp=j=1nmaxdmax(j)+j=1nmindmin(j)avec{dmax(j)=Kb2.C.[Zmax(j)]bdmin(j)=Kb2.C.[|Zmin(j)|]b.(15)

The extrapolation in the temporal domain of this partial damage over the total stress duration Tt is performed in a deterministic manner, using the following rule of proportionality; and can therefore be used to calculate the cumulative damage Dc, respecting the same assumptions used for the calculation of the average cumulative damage E [Dc(ξ)]: Dc=TtTref.Dp=TtTref.Kb2.C.[j=1nmax[Zmax(j)]b+j=1nmin[|Zmin(j)|]b].(16)

This calculation procedure is then applied to all of the systems with one dof in question for each of the natural frequencies f0, which leads to the following expression of the temporal FDS, noted FDStpD(f0), of a deterministic nature: FDStpD(f0)=CKb.Dc=Tt2.Tref.[j=1nmax[Zmax(j)]b+j=1nmin[|Zmin(j)|]b].(17)

Note: As in the case of the SRS [12] and the temporal ERS, it is recommended here to sample the random excitation process x¨(t) at a frequency at least 10 times greater than the maximum frequency of the analysis of the RS in question, in order to correctly evaluate the temporal FDS.

This technique of calculating the temporal FDS of a deterministic nature, although very easy to develop from a numerical point of view, is very often inadequate for the secure sizing of structures, which necessarily requires taking into account the statistics of the cumulative damage of stress processes of a stochastic nature. This therefore necessarily requires the use of probabilistic rather than deterministic methods, that are capable of extrapolating the values of the (1−α) quantiles of cumulative damage for a given total stress duration Tt.

4 Securing structures and design risk control

To secure the fatigue design of mechanical structures subjected to stationary random processes of prolonged stress Tt we should now focus on the statistical characteristics of the accumulated fatigue damage levels generated by these processes over time Tt, applying a risk of exceeding them α, which is determined according to the criticality level of the structure to be sized.

In practice, the three levels of risk, commonly used in the mechanical design of military land vehicles, are:

  • α = 10%: Risk associated with the development of equipment and/or mechanical structures “common”, not part of a functional chain of important performance for the end customer. This is particularly true of the components that are associated with secondary functions of the System;

  • α = 1%: Risk associated with the development of equipment and/or mechanical structures “standards”, part of a functional chain of important performance (sensitivity) for the end customer, but the level of customer demand remains moderate. This is the case of the components that are associated with the main functions of the system, for which the functional performance levels specified by the end customer are achievable, safe, including the vibratory stresses taken into account in the development;

  • α = 0, 1%: Risk associated with the development of “secure” mechanical equipment and/or structures as part of a functional chain of performance which is important for the end user, but for which the level of customer demand is high. This is the case of the components that are associated with the main functions of the system, for which the high functional performance levels specified by the end customer can only be achievable, without risk, if the stresses taken into account in the development are fully controlled.

Due to these technical considerations, it was therefore considered important in the years 2005–2010 [6,13] to propose major changes, to the RS techniques, to better understand the risks of exceeding vibration fatigue damage levels in the case of the excitation processes generated by military land vehicles. This work was first written in the context of random Gaussian processes, which helped to establish the notion of spectral XFS with an up-crossing risk α, which replaces the case of spectral FDS, by definition associated with the accrued level of damage.

In section 5 below, we therefore recall the main results of this notion of spectral XFS, which is subject to the case of Gaussian processes. It will be completed positively by innovative approaches that can be used to propose concepts of temporal XFS, which are now capable of taking into account the non-Gaussian nature of the random processes generated by Land vehicles. These techniques have been favorably used in the ARAVIS, VBCI and CAESAR vehicles, which are products developed by Nexter Systems on wheeled mobility platforms (see Fig. 3).

thumbnail Fig. 3

ARAVIS (4 × 4 of 12,5 Tons), VBCI (8 × 8 of 29 Tons) and CAESAR (Artillery of 155 mm, with carriers 6 × 6 RTD or UNIMOG).

5 XFS computing approach

As presented above, it is essential now to estimate the statistical model of cumulative damage levels associated with random stress processes over the entire duration Tt of the excitation process, regardless of the statistical nature of the excitation process (Gaussian or non-Gaussian). This results in no longer promoting the use of deterministic temporal FDS techniques, which seemed to be widespread among designers, since they did not have probabilistic methods suitable for the target technical issue. Indeed, the objective hitherto sought is to assess the threshold level of cumulative damage [(1−α)-quantile of the random variable Dc] over the total excitation period Tt with a constant up-crossing risk α for each natural frequency f0 imposed by the designer and/or specifier; whereas in the case of spectral FDS (see Sect. 5.1), this risk α is a function of the frequency values f0 retained.

These up-crossing risk RS α were designated XFS in order to differentiate them from the term FDS and now cover the case of Gaussian and/or non-Gaussian stationary random processes, as presented below.

5.1 Calculation of spectral XFS

Based on the work conducted in 2007 by the author on a conservative approach associated with the case of stationary, zero mean and narrow-band Gaussian random processes, an analytical approach of the spectral XFS was established on the basis of the following assumptions, which are complementary to those defined for calculating the spectral FDS:

  • the excitation process x¨(t) of duration Tt is stationary, ergodic, Gaussian, zero mean and broadband of the white noise type (a highly restrictive assumption and not consistent with the physical reality of Defense sector land vehicles);

  • the damping level ξ the standard system retained is weak (ξ ≤ 0, 05) corresponding to a quality factor Q ≥ 10;

  • the numbers of stress cycles associated with the various standard systems are high and must respect the following order relation, namely: (ξ . f0 . Tt ≥ 1);

  • the cumulative damage distribution law Dc(ξ) is assumed to be represented by a 2-parameter Weibull distribution (β˜,σ˜), in which β˜ and σ˜ respectively represent the shape parameter and the scale parameter of Weibull distribution in question.

Indeed, insofar as the cumulative damage Dc(ξ) is a positive support random variable, it is clear that the statistical model to retain cannot be Gaussian as claimed by some authors. The choice of a Weibull model has been preferred, given its exponential decay, which is the same as that associated with the Rayleigh model characterizing the local maxima of the narrow-band random stress process in question. Thus, based on these assumptions, the coefficient of variation of the cumulative damage Dc(ξ) can be calculated analytically and leads to the following approximate result, namely [14]: CV[Dc(ξ)][2.Q.f1(b)n]1/2wheren=f0.Ttwheref1(b)=exp[c2.b2+c1.b+c0]forb[3,20]andc2=2,51419E03,c1=5,04776E01,c0=2,38747E+00.

The shape of the function f1(b) is described in Figure 4 and is an approximate model of the infinite series from which it derives [6].

Due to these technical considerations, it appears that the 2-parameter Weibull model (β˜,σ˜)for the cumulative damage is now adjustable by the method of moments, by endeavoring to readjust for each standard system, the average value of cumulative damage provided by equation (14) and its coefficient of variation provided by the previous equation (18). Therefore, it is shown that the spectral XFS with up-crossing risk α can be calculated from the classical spectral FDS, using the following simple expression: XFSspBE(f0)FDSspBE(f0)=[ln(1/α)]λ˜Γ(1+λ˜)=g4(α,b,Q,n)whereλ˜issuchash(λ˜)=0 whereh(λ)=Γ(1+2.λ)Γ2(1+λ)μandμ=2.Q.f1(b)n+1wheren=f0.Tt.

The coefficient λ˜ obtained by dichotomy according to equation (19) is directly connected to the shape parameter β˜ of the Weibull distribution, selected to describe the statistical model of cumulative damage Dc(ξ) over the period Tt of stress during the random excitation process in question, since one can write: β˜=1/λ˜.(20) Thus in the case of electronic equipment (b = 4) and mechanical equipment (b = 8), we find that the coefficient λ˜ depends on the number n of cycles associated with the stress response process, for a given value of Q (usually taken to be equal to 10 in the case of customization approaches [7,15]). Graphically, we find that the value of this coefficient λ˜ decreases with the number of cycles n, for a given parameter b. In addition, for a given value of numbers of cycles, the greater the value of b (representative of the slope of the Basquin model) increases, the higher the coefficient λ˜, as shown in the following Figure 5.

By now injecting the value of the coefficient λ˜ into the analytic expression of the XFS provided by equation (19), it is possible to plot the relationship between the spectral XFS and the spectral FDS, conventionally used in customizing trials for a α given up-crossing risk. We thus obtain the following graphs.

At this stage, it should be noted that the spectral XFS is particularly relevant for taking into account the fact that the number of stress cycles n in question is low, the risk α selected by the specifier and/or designer is low and that the Basquin slope for the material b of the structure is high, as shown by the curves of Figure 6, taken from work in 2007 [6].

Based on the spectral XFS model presented above and associated with the case of narrow-band stationary Gaussian random processes, it is now possible to assess the analytical expression of the value of the risk α0 that is implicitly associated with the case of expression of the spectral FDS. Indeed, it seems important to correctly identify the up-crossing risk with which the spectral FDS concept is associated, which is very widely used in the sector for the experimental validation of products, as indeed it is today in the field of product design, which increasingly uses this technical contour to substantiate the fatigue behavior of structures subjected to a vibratory environment, in terms of reliability. Due to these technical considerations, according to equations (19) and (20) we can state: α0=exp[[Γ(1+1β˜)]β˜].(21)

Graphically for given values of b and Q, the shape of the up-crossing risk α0 depending on the number n of cycles shown in Figure 7 below, namely:

This therefore shows that the spectral FDS of a stationary Gaussian random excitation process is associated with an up-crossing risk α0 which is not constant and which is not equal to 50% as some authors tend to say. This risk α0 depends on the b and Q parameters, but especially the number of stress cycles, as seen by the structure in terms of stress response, as shown in Figure 7 above. The greater the number of cycles increases, the greater the up-crossing risk α0 increases, reaching an asymptotic value of the order of 57%, in the case of structures characterized by current values of b ranging from 8 to 12, for example. We also note that for a value of given number of cycles n, the higher the value of b, the lower the risk α0, showing the need to clearly identify the nature of the critical areas of failure, the reliability of which needs to be evaluated, while controlling the associated risk.

thumbnail Fig. 4

Analytical expression of the f1(b) function, associated with calculation of spectral XFS.

thumbnail Fig. 5

Coefficient λ˜, associated with spectral XFS for b = 8 and b = 12.

thumbnail Fig. 6

Relationship between the spectral XFS and classical spectral FDS for b = 8 and b = 12.

thumbnail Fig. 7

Up-crossing risk, associated with the classical spectral FDS for b = 8 and b = 12.

5.2 Calculation of the temporal XFS

In the case of a non-stationary Gaussian process of duration Tt characterized by a single stress time sample z(t) of limited duration Tref (Tref<Tt) it is no longer possible to rely on spectral approaches to assess its XFS. Although these non-Gaussian processes will always be characterized in terms of frequency by their PSD, the probabilistic structure of partial damage is only accessible by analyzing their first statistical moments, evaluated over the period of time [0,Tref]. To do so, the temporal approach [15] (called the DBM Method, for Disjoint Block Method) implemented in the case of a stationary non-Gaussian random process consists in now analyzing the local fatigue damage of the stress process by blocks of time Tb.

To simplify, here we are placed in the case of non-Gaussian random processes, which we presume to be order-2 stationary (by restricting the stationarity in terms of rms values). The non-stationarity of the processes, induced for example by the vehicle speed does not constitute a difficulty per se, given that in order to overcome it, it is useful to classify the blocks of time data Tb, retrieved from the x¨(t) non-Gaussian and non-stationary excitation process, in terms of rms values x¨eff.

Due to these technical considerations, and after calculating using the FOH technique on each standard system the relative displacement z(t) over the reference duration Tref of stationary non-Gaussian excitation x¨(t), the responses z(t) of the standard systems are as if x¨(t) cut into blocks of time Tb, without any overlapping. The number of analysis blocks N is then equal to the ratio (Tref/Tb) and the level of damage by partial fatigue d(i) of each block  (i = 1 to N) is estimated. The damage by partial fatigue d(i) of each block of time Tb is then evaluated in the temporal domain by equation (15) above, under Basquin's law and Miner's rule.

For each response z(t) of non-Gaussian type, the sample of N damage values {d (i) ,  i = 1 to N} thus evaluated over time Tb is then an N-random sample that can be used to characterize the first statistical moments of the damage by partial fatigue D(ξ) of z(t), produced over time Tb. The first two statistical moments of the random variable D(ξ) are then easy to calculate and have as their expressions the following relationships in terms of average D¯and variance σˆD2, which here are two non-biased statistical estimators, namely: D¯=1N.i=1Nd(i)   where    N=TrefTb   and    σˆD2=1N1.i=1N[d(i)D¯]2.(22)

The damage by accumulated fatigue Dc(ξ) associated with the (order-2) stationary non-Gaussian and zero mean random stress process over the entire duration Tt is then defined by the equation (23) below. To do so, once again we use the Miner principle, considering an extrapolation factor over time (denoted M), which is associated with the excitation process x¨(t) measured or estimated over Tréf, given that the objective is now to estimate the accumulated fatigue damage over the total period (Tt  >  Tref) by including its statistics. Dc(ξ)=i=1MDi(ξ)avecM=Tt/Tb M:extrapolationfactorovertime.(23)

The M random variables representative of the damage by partial fatigue Di(ξ) are assumed here to be independent and identically distributed (i.i.d), which therefore allows us to apply the central limit theorem (CLT) to define the law of asymptotic probability of cumulative damage Dc(ξ). This hypothesis, related to the use of CLT is only acceptable if the extrapolation factor over time M is sufficiently large, which is a statistical validation criterion to be satisfied first of all, namely:

  • if the random variables laws Di(ξ) are close to a normal distribution, then (M > 4) is sufficient;

  • if the random variables laws Di(ξ) are moderately close to a normal distribution, then (M > 12) is sufficient (for symmetrical distributions around their mean);

  • if the random variables laws Di(ξ) are not close to a normal distribution, then (M ≥ 100) is to be considered (case of highly asymmetric distributions, which corresponds here to our case study, since the cumulative damage by fatigue is by definition positive support.

In practice, one can be satisfied with an extrapolation factor with a minimum duration of 30, but it would be unwise to go below, to avoid estimating irrelevant and unrepresentative fatigue damage (1−α) quantiles of the excitation process x¨(t) in question. Therefore the asymptotic distribution of fDc(y) accumulated fatigue damage is comparable to a normal distribution (under the assumption that the random variables Di(ξ) are i.i.d) centered on D¯c and standard deviation σˆDc. It is thus written as follows: fDc(y)Mhigh1σˆDc.φ(yD¯cσˆDc)whereφ(y)=12.πexp(y22) D¯c=M.D¯andσˆDc=M.σˆDandwhereM:extrapolationfactorovertime, where ϕ(y) represents the reduced centered normal distribution (zero mean and unit variance), often denoted N(0, 1) in statistics.

Due to these technical considerations, it is possible to evaluate the expression of temporal XFS, associated with a given up-crossing risk α, based on the definition of (1−α) quantile of the distribution law fDc(y) defined in equation (24) above. We then obtain the analytical expression of (1−α) quantile, denoted Dc,α, which is homogeneous in the quantity Dc of equation (16) of the deterministic temporal approach. Dc,α(f0)=D¯c.[1+2.CV[Dc(ξ)].erf1(12.α)]whereCV[Dc(ξ)]=σˆDc/D¯canderf(x)=2π.λ=0xexp(λ2).dλ.

Based on the unity of an FDS which corresponds to that of a fatigue damage (dimensionless value), in relation to the standard coefficient (Kb/C), it is then possible to write the analytical expression of the temporal XFS as follows: XFStp(f0)=CKb.Dc,α=M. D¯.[1+2M.CV[D(ξ)].erf1(12.α)] whereCV[D(ξ)]=σˆD/ D¯andwhere(D¯,σˆD)aredefinedby(22), where CV [D(ξ)] is the experimental coefficient of variation of partial damage, associated with the random stress response process, evaluated through on standard system (f0, Q), for a measuring time equal to Tref.

By definition, the temporal FDS noted FDStp(f0), being associated with the up-crossing risk equal to 50%, we see that we can immediately write: XFStp(f0)FDStp(f0)=1+2M.CV[D(ξ)].erf1(12.α)whereCV[D(ξ)]:experimentalcoefficientofvariationofpartialdomage.

This shows that the ratio between the temporal XFS and the temporal FDS is a function of α, of M and of CV [D(ξ)] which can be represented graphically as follows for an extrapolation factor over time M of 10, 100 or 1000 (plotted respectively in Figs. 810  below).

For an extrapolation factor M of given time, we see that the ratio between the temporal XFS and temporal FDS increases with the increase in the experimental coefficient of variation of partial damage CV [D(ξ)] as well as with the reduction in the self-imposed up-crossing risk α.

By drawing on the results of previous graphs, we show that for a given stationary non-Gaussian random stress process (with a given coefficient of variation of partial damage) and a given level of risk, the longer the stress period (the greater the increase in M), the lower the ratio between the temporal XFS and the temporal FDS.

These temporal FDS and/or temporal XFS calculation techniques are conditioned by the choice of the length Tb of the analysis block of partial damage, associated with the non-Gaussian random process in question. The scientific bases developed previously, show that it is important to respect some basic criteria to define the optimal value of this analysis time Tb, namely:

  • reduce Tb as much as possible, in order to obtain a number N (see Eq. (22)) sufficiently high to evaluate, with a good level of accuracy the first two statistical moments of partial damage D(ξ). This also means increasing the extrapolation factor over time M of the process (see Eqs. (23) and (24)) and to substantiate the use of the CLT to evaluate the statistical model of cumulative damage Dc(ξ);

  • increase Tb as much as possible, to ensure the existence of a sufficient number of cycles to evaluate the partial damage that is not strongly correlated between two successive blocks, in order to satisfy the independence assumption of the previous equation (23).

These two criteria being contradictory by nature, the choice of the period Tb of the so-called DBM Method in practice depends on the reference measurement time Tref and the minimum value of f0 required.

thumbnail Fig. 8

XFStp(f0)/FDStp(f0) pour M = 10.

thumbnail Fig. 9

XFStp(f0)/FDStp(f0) pour M = 100.

thumbnail Fig. 10

XFStp(f0)/FDStp(f0) pour M = 1000.

6 Conclusions

The RS techniques, historically developed in terms of temporal FDS for shock phenomena, and spectral FDS for random vibration phenomena, do not make it possible to properly characterize non-Gaussian random processes in terms of fatigue damage.

To do so, temporal XFS techniques were developed and tested under the Land Weapons programs (VBCI, ARAVIS and CAESAR), highlighting their characteristics of robustness and speed of calculation.

Faced with the increased use of RS for the structural design and drafting of qualification test specifications for the equipment, these techniques of temporal XFS calculations have recently been proposed as part of AFNOR standardization work. Taking into account the lessons learned acquired by Nexter Systems on the use of these innovative approaches developed by the author in an industrial context, the members of the AFNOR Group have adopted them as part of NF X50-144-3 [16], for their ability to positively expand RS techniques to very widespread processes in industry, and so far poorly characterizable by conventional temporal FDS and/or spectral FDS approaches.

As illustrated by this standard, these tools (temporal XFS) for characterizing fatigue damage, with an imposed up-crossing risk, thus make it possible to greatly expand the scope of the calculation methods of RS, by now covering the case of non-Gaussian processes (whether narrow-band or broadband).

The non-Gaussian nature of the processes now being taken into by these temporal approaches of a stochastic nature, it is no longer necessary to use the deterministic temporal approaches, described in section 3.2, in particular since they are non-conservative.

References

  1. J.F. Imbert, Analyse des Structures par Eléments Finis- Cépadues Editions [Google Scholar]
  2. C. Lalanne, Vibrations et chocs mécaniques- Hermes Science publications, copyright Lavoisier, 1999: Tome 5-Élaboration des spécifications [Google Scholar]
  3. J.S. Bendat, Principles and applications of random noise theory, John Wiley & Sons, New York, 1958 [Google Scholar]
  4. R.W. Hamming, Numerical methods for scientists and engineers, 2nd edn, ISBN 0486 652416 [Google Scholar]
  5. D.O. Smallwood, An improved recursive formula for calculating shock response spectra, shock and vibration Bulletin: Proceedings of the 51st Symposium on Shock and Vibration, San Diego, California, 1980 [Google Scholar]
  6. B. Colin, La probabilisation des SDF: une réponse à la comparaison des sévérités d'environnement déterministe et stochastique, en termes d'endommagement par fatigue, 25–27 septembre 2007, ASTELAB 2007 [Google Scholar]
  7. NORMDEF 0101, Guide d'application de la démarche de personnalisation en environnement mécanique, DGA 2009, Edition 01, Julliet 2009. [Google Scholar]
  8. B. Colin, Approche stochastique de la Fatigue vibratoire: comparaison entre les méthodes Peak-Valley et Rainflow Counting, Revue de Métallurgie, EDP Sciences, 2012, DOI:10.1051/metal/2012024 [Google Scholar]
  9. M.A. Miner, Cumulative damage in Fatigue, J. Appl. Mech. 67 (1945) A159–A164 [Google Scholar]
  10. S.O. Rice, Mathematical analysis of random noise, Bell Labs Tech. J. 23 (1944) 282–332 [Google Scholar]
  11. A. Preumont, Random vibration and spectral analysis, Kluwer Academics Publishers, 1994 [CrossRef] [Google Scholar]
  12. C. Lalanne, Vibrations et chocs mécaniques- Hermes Science publications, copyright Lavoisier, 1999: Tome 2-Chocs mécaniques [Google Scholar]
  13. B. Colin, La probabilisation des SRE: une réponse à la comparaison des sévérités d'environnement déterministe et stochastique, en termes de contraintes extrêmes, 25–27 Septembre 2007, ASTELAB 2007 [Google Scholar]
  14. Stephen H. Crandall, W.D. Mark, Random vibration in mechanical systems, Hardcover, 1963 [Google Scholar]
  15. B. Colin, Characterization of random Gaussian and non-Gaussian stress processes in terms of extreme responses, in: Matec Web of Conferences, EDP Sciences, DOI:10.1051/matecconf/20152006001. [Google Scholar]
  16. AFNOR, Norme NF X50-144-3, Démonstration de la tenue aux environnements-Conception et réalisation des essais en environnement-Partie 3: Application de la démarche de personnalisation en environnement mécanique, 2014 [Google Scholar]

Cite this article as: B. Colin, Characterization of random Gaussian and non-Gaussian stress processes in terms of vibration fatigue, Mechanics & Industry 18, 809 (2017)

All Figures

thumbnail Fig. 1

Standard system with one dof.

In the text
thumbnail Fig. 2

Basquin model and convention of writing.

In the text
thumbnail Fig. 3

ARAVIS (4 × 4 of 12,5 Tons), VBCI (8 × 8 of 29 Tons) and CAESAR (Artillery of 155 mm, with carriers 6 × 6 RTD or UNIMOG).

In the text
thumbnail Fig. 4

Analytical expression of the f1(b) function, associated with calculation of spectral XFS.

In the text
thumbnail Fig. 5

Coefficient λ˜, associated with spectral XFS for b = 8 and b = 12.

In the text
thumbnail Fig. 6

Relationship between the spectral XFS and classical spectral FDS for b = 8 and b = 12.

In the text
thumbnail Fig. 7

Up-crossing risk, associated with the classical spectral FDS for b = 8 and b = 12.

In the text
thumbnail Fig. 8

XFStp(f0)/FDStp(f0) pour M = 10.

In the text
thumbnail Fig. 9

XFStp(f0)/FDStp(f0) pour M = 100.

In the text
thumbnail Fig. 10

XFStp(f0)/FDStp(f0) pour M = 1000.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.