A widely used approach for the first crossing reliability evaluation of structures subject to nonstationary Gaussian random input is represented by the direct extension to the nonstationary case of the solution based on the qualified envelope, originally proposed for stationary cases. The most convenient way to approach this evaluation relies on working in the time domain, where a common assumption used is to adopt the modulation of stationary envelope process instead of the envelope of modulated stationary one, by utilizing the so-called “preenvelope” process. The described assumption is demonstrated in this work, also showing that such assumption can induce some errors in the envelope mean crossing rate.
1. Introduction
One of the most important indexes for engineers and designers is the accurate quantification of structural safety. Its correct evaluation is connected with loads’ nature and failure type. Moreover, a wide class of engineering problems deals with reliability evaluation of structures subject to loads characterized by an intrinsic probabilistic nature and whose correct description could be obtained by means of random process. Besides, there are situations where suitable modelling of these loads requires taking their nonstationary characteristics into account [1–4]. Examples of such conditions include ground acceleration during an earthquake and structural vibration induced by gust loading.
In these cases, the probabilistic characterization of structural loads turns the random dynamic analysis to be the most useful method in achieving an accurate and qualified structural reliability evaluation [5–10]. On the other hand, using this method considerably improves the complexity of analysis and the computational costs for reliability estimation. For instance, exact solutions exist only in few cases. These limitations have substantially reduced the applications of the stochastic methods from design engineering in real structural safety problems. However, many different approximations have been proposed to overcome these practical restrictions.
In particular for Gaussian inputs and linear systems, the first- and second-order moments completely define the statistics of the response. This problem does not encounter serious difficulties in both spectral and time domains relative to a stationary case and with regard to first crossing failure. However, the application of Rice’s reliability formulation presents some problems related to the required knowledge of hazard function, defined as the probability density function (PDF) of collapse at a given instant, under the restriction that crisis has not happened yet. Due to complexity in its evaluations, a simplified approach is based on the independent crossing hypothesis, using then an unconditioned PDF as hazard function. This hypothesis, which assumes threshold crossing distributed as a Poisson process, can be acceptable only for wideband process and high barrier values, where barrier crossing can be really assumed to be independent. Besides these restrictions, Poisson approach is too poor and conservative for a correct reliability evaluation [11, 12].
A more recent approach for the first passage problem has been developed in [13] for linear and nonlinear systems driven by Poissonian and normal white noise input. In addition, in [14], a numerical path integral solution approach is developed for determining the response and first passage probability density functions (PDFs) of nonlinear oscillators subject to evolutionary broad-band stochastic excitations. A strategy for estimating first excursion probabilities for linear dynamical systems involving uncertain structural parameters subject to Gaussian excitation has been developed recently in [15].
Among different methods, an approach, widely used nowadays, is that proposed by Vanmarcke [16] and is based on the envelope process statistics, which have lower dependencies on threshold crossing events. This approach has been described for probabilistic assessments of structural failure using the barrier crossing rates obtained by using the envelope probabilistic information of a number of parameters derived from the PDF generally known as spectral moments. Their first interpretation has described geometrical moments of the one-sided PSD with respect to the central origin [17]. This frequency definition has been extended to time domain first by Di Paola [18] who has found suitable relationships between the correlation function, its Hilbert transform, and the spectral moments. This approach, also known as “nongeometrical formulation” [19], is able to extend the stationary approach to nonstationary cases, so eliminating the significant difficulties derived by using a geometrical interpretation [20]. It is based on the use of the so-called “preenvelope” process defined as the complex process whose real part is the given process and the imaginary part is the Hilbert transform of the real part, A=X2+X⌢2 [21]. It is worth mentioning that different definitions have been advanced for envelope processes such as a process joining the peaks of the response [22]. Using a “nongeometrical” approach by covariances of the preenvelope process, defined in time domain, it is possible to obtain the same quantities needed for reliability evaluation derived by a “geometrical” one [19, 23]. The two approaches coincide in a stationary case; however, this result is lost in a nonstationary condition. The reason can be explained as follows: for physical and mathematical questions (different from zero before beginning of input), the exact nonstationary preenvelope process is replaced by a simplified one, the modulation of stationary preenvelope process, which has a more realistic characterization. In addition, this simplification, commonly used by different authors for reliability in nonstationary conditions [24–26], induces some mistake in the evaluation of the joint probability density function (JPDF) of the preenvelope and its first time derivative process. In the present work, the agreement between the approximate analytical formulation and the numerical simulation is shown for different transient conditions and structural damping. Moreover, this study reveals that this simplified approach can be used to transpose important characteristic of covariances matrices needed for reliability evaluation. By utilizing this result and by adopting a differential approach, it is demonstrated that useful simplifications can be applied to the computation of the mean crossing rate. Finally, a benchmark between simulation and different analytical formulae for mean crossing rate is developed. Also, a simple modified formula is finally proposed to increase the agreement with numerical results.
2. Structural Reliability Evaluation by Classical Poisson and Vanmarcke Approaches
Mechanical safety or reliability R(t) at a fixed time t throughout the life of a structure subject to static or dynamic loads can be defined as the collapse survival probability, where the collapse is a partial or a total static damage in a fixed time interval [t0, t]. Time t0 is typically the initial observation time assumed as the end of the manufacturing process. In this situation, the collapse probability at time t0 is usually assumed equal to 0. It is clear that the definition of structural collapse (or in general, failure) plays an important role in reliability assessment, considering the numerous possible meanings (not only of mechanical kind) of this condition. Two of the most common mechanical undertakings are related to fatigue phenomenon and first threshold crossing failure.
In this section, a simple one-degree-of-freedom, viscous-elastic system subject to random dynamic actions is analysed, for which failure is related to displacement barrier crossing. Precisely, the system crisis happens when structural displacement X(t) exceeds an admissible value β for the first time. In many structural situations, typically if the response is symmetric, the most used definition of reliability R(t) is given by the probability that the system will not have a failure associated with a bilateral threshold crossing of level β by process displacement Xt (double barrier problem) in the time interval [t0, t]. Its formal definition can be written as(1)Rt=PXt<β; t∈t0,t.
An alternative method is using its complementary collapse probability (failure) Pf(t) = 1 − R(t), which can be defined as the probability that the double barrier β will be exceeded by the process Xt during the observation interval [11]. Reliability evaluation requires the knowledge of the hazard function h(t), also named as “decay rate function”, which is related to R(t) according to the following equation:(2)Rt+dt=1−htdtRt.
The hazard function h(t) is defined as the quantity whose product with the infinitesimal time interval dt supplies the structural collapse probability for threshold crossing during [t, t + dt], under the condition that at time t, the response process has not yet shown any other threshold crossing, which is formally represented as(3)htdt=PC∣S,where C is the event “threshold crossing” in [t, t + dt] and S is the absence of barrier excursions before time t.
Using (2), it is possible to write(4)R′t=limdt⟶0Rt+dt−Rtdt=limdt⟶01−htdt−1Rtdt=−htRt.
The formal solution of (4) is ∫dR/R=−∫htdt, which gives at last the well-known integral reliability formulation:(5)Rt=R0exp−∫t0thtdt.
The exact hazard function formulation is still an open question due to the fact that it has been solved only in few cases. In the original Rice’s formulation [27], the hazard function has the following form:(6)ht=∫−∞+∞x˙pXX˙β,x˙;t∣Sdx˙,where pXX˙β,x˙;t∣S is the conditioned JPDF of X and its first time derivative X˙t processes. The difficulties related to its evaluation have induced the development of approximate solutions. One of the most common approximations uses the unconditioned JPDF instead of the conditioned one. Using this approximation, the hazard function h(t) is replaced by υt, known as unconditioned threshold crossing rate. In addition, the relative reliability is evaluated in the so-called Poisson hypothesis, because the independent threshold crossings assumption means, according to a Poisson process, that they are rare and without memory events. This approximation has been demonstrated to be suitable for wideband processes and high barrier levels, when the threshold crossing is effectively a rare event, so that the distribution of crossings has an asymptotic Poisson distribution in this situation. On the contrary, this approach could be particularly poor for narrowband processes and for relatively low barrier levels, since it is too conservative. In fact, the barrier crossings in these situations tend to happen in clumps, completely neglecting the independence hypothesis.
A frequently used solution to overcome the limitation discussed above is obtained by using the approach of Vanmarcke [16, 17], which is based on a two-state Markov crossing assumption. This formulation uses the statistical information of “envelope process A(t),” which is introduced to consider the effects of the clumps, and thus reducing the errors that originate from the independent crossing assumption. The envelope process has the general physical meaning of a nonnegative random process that joins the peaks of Xt. Using the Vanmarcke formulation, a two-state Markov process is introduced to describe the two possible conditions (over and below) of the envelope process A(t) with respect to the assigned barrier β and to define the so-called E-crossing event (i.e., an up-crossing of the envelope). The author has developed an analytical formulation for hazard function by using the Poisson assumption for the arrivals of E-crossing and by taking into account only the cases where a D-crossing (up-crossing of Xt) follows it (qualified envelope process). The resulting formulation has a stronger agreement with numerical results as compared to the original simpler one that is based on the Poisson assumption directly applied to the system response process [11, 13].
For a stationary process, the formulation that furnishes the hazard function for a bilateral threshold level β is(7)υVMβ=υXβ1−expυA+β/υXβ1−υXβ/υX0,where υXβ=2υX+β is the mean crossing rate for a single barrier evaluated with the Poisson hypothesis:(8)υXβ=∫−∞+∞x˙pXX˙β,x˙dx˙.
And the single barrier mean crossing rate of the envelope process υA+β is defined using the same approach as(10)υA+β=∫0+∞a˙pAA˙β,a˙da˙.
Its evaluation requires the JPDF pAA˙a,a˙ of At and A˙t, which means the stochastic characterization of these two processes. Using this solution for a nonstationary process, the hazard function for threshold level β is [11](11)υVmβ,t=2υX+β,t1−expυA+β,t/2υX+β,t1−υX+β,t/υX+0,t,where υX+β,t is the mean crossing rate for a single barrier β evaluated under the Poisson hypothesis.
3. Envelope Process Definition
In a stationary condition, the envelope process can be defined as the modulus of a complex process whose real part is the original one and the imaginary part has to be opportunely defined in order to satisfy the following equation:(12)Ast=Ust=Xs2t+iYs2t.
The choice of Yst has to guarantee the physical meaning of the envelope process, which implies a smoothly curve able to connect response peaks. Different approaches have been employed to choose a suitable formulation of Yst [22]. Usually, the preferred ones are those based on the use of the response processs first time derivative Yst=X˙st/ωc [21], or on its Hilbert transform (Hilbert transform of function ft is Hft=f⌢t=1/π∫−∞+∞fτ/t−τdτ) Yst=X^st [21]. The latter is the approach adopted in the current work.
White noise or filtered Gaussian uniform modulated processes is a widely used model for representing nonstationary loads [28]. As shown in (13), it is obtained as the product of a stationary Gaussian white noise process St and a deterministic time modulation amplitude controlling function φt, which must satisfy these conditions: φt≥0 for t>0; φt=0 for t<0; and maxφt=1:(13)Nt=φtSt.
The stationary process is usually represented by a zero-mean white noise process wtRWWτ=2πS0δτ; the covariance of the corresponding modulated process Ft=φtwtμNt=φtμSt,RNNt1,t2=φt1φt2RSSt1,t2,σN2t=φ2tσS2t is given by(14)RFFt1,t2=2πS0φt1φt2δt2−t1,where δt is the Dirac delta function.
A multiplicative factor equal to 1/2 has to be taken into account due to the consideration that, for a generic analytical stochastic process vt=1/2ut+iu⌢t (where ut is real), the covariance is given by(15)Rvτ=Ruτ+iR⌢uτ.
In this case, Ust becomes a complex process whose imaginary part is the Hilbert transform of the real one; that is, it can be expressed as the system response to the Hilbert transform of a real input modelled as a white noise process wst:(16)Ust=12Xst+iX^st,where Xst=∫0tht−τwsτdτ and X^st=∫0tht−τwsτdτ; ht is the impulse response function, which for a SDOF system can be written as ht=1/ω¯e−ξωtsinω¯t, where ω¯=ω1−ξ2.
In a more compact form, Ust can be defined as the analytical process response of a linear system to the analytical input process:(17)Fst=12wst+iw^st.
The covariance matrix of the stationary vector process Ust can be expressed as RUsUsτ=RXsXsτ+iRXsX⌢sτ, due to the following properties, generally true only in the stationary case [28]:(18a)RXsXsτ=RX⌢sX⌢sτ,(18b)RXsX⌢sτ=−RX⌢sXsτ=R⌢XsXsτ.
In a nonstationary condition, the structural response envelope can be defined in a way similar to the stationary case, by introducing a nonstationary complex process Ut. However, a different approach is necessary to define the complex part, due to physical and analytical reasons.
Using a nonstationary modulated process input Ft=φtwt, the direct Hilbert transforms applied to the linear system response(19)X^t=∫0tht−τHφτwτdτ,is different from zero in the time interval −∞,0, and the corresponding envelope too, conflicting with the physical meaning of an event that actually begins at time t0=0. For this reason, the modulated Hilbert transform input process F⌢t=φtw⌢t, first proposed by Di Paola and Muscolino [18, 29, 30], is used to define the so-called analytical preenvelope input process:(20)Ft=12φtwt+iw⌢t.which corresponds to the analytical preenvelope response:(21)Ut=12Xt+iX⌢t,where the real and imaginary parts are(22a)Xt=∫0tht−τφτwτdτ,(22b)X˜t=∫0tht−τφτw^τdτ.
In this way, it can be written in [18]:(23)RUUt1,t2=RXXt1,t2+iRXX⌢t1,t2.
By extending the following properties to a nonstationary case:(24a)RXXt1,t2=RX⌢X⌢t1,t2,(24b)RXX⌢t1,t2=−RX⌢Xt1,t2.which are in general not true using the direct Hilbert transformation when applied to the modulated input process.
4. Standard Space State Covariance Response
As a main direct consequence of the above assumptions, working with the two-state space vectors is expressed in [6]:(25a)Zt=XtX˙t,(25b)Z˜t=X˜tX˙˜t.
The correlation matrix is given by(26)Z¯tZ¯t=RZZt=σX2tγXX˙tγX˙XtσX˙2t.and is symmetric as it is well known, being XtX˙t=γX˙Xt=γXX˙t.
Similarly the matrix(27)Z¯tZ¯⌢t=RZZ⌢t=0γXX˙⌢tγX˙X⌢t0.is asymmetric, being X˙tX⌢t=γX˙X⌢t=−γXX˙⌢t.
Then, it is possible to write the differential equations in the space state as(28a)Z˙¯t=AZ¯t+B¯t,(28b)Z˙⌢¯t=AZ⌢¯t+B⌢¯t,where the state matrix A in the case of a linear SDF system is(29)A=01−ω02−2ξ0ω0.
And the input vectors B¯t and B⌢¯t are, respectively(30a)B¯t=φt0wt,(30b)B⌢¯t=φt0w⌢t.
As shown earlier in (10), the nonstationary envelope process mean rate threshold crossing υAβ,t needs the JPDF of At and A˙t, which are related to both processes Ut and U˙t, being A˙t=U˙t=1/2X˙t+iX˙⌢t. It is then possible to introduce the analytical envelope state vector process as(31)Ψ¯t=UtU˙t=XtX˙t+iX⌢tX˙⌢t.
The covariance matrix is given by(32)RΨ¯Ψ¯t=Ψ¯tΨ¯tT=Z¯tZ⌢¯tZ¯tTZ⌢¯tT=Z¯Z¯TZ¯Z⌢¯TZ⌢¯Z¯TZ⌢¯Z⌢¯T.
It can be rewritten in a more compact form by decomposing the 4 × 4 matrix into four 2 × 2 submatrices [6]:(33)RΨ¯Ψ¯t=RZZtRZZ⌢tRZ⌢ZtRZ⌢Z⌢t.
Using the mean value Hilbert properties previously shown, it is possible to demonstrate that the 4 submatrices are two by two equal (RZZt=RZ⌢Z⌢t) and opposite (RZZ⌢t=−RZ⌢Zt).
Therefore, RΨ¯Ψ¯t can be defined just knowing the two 2 × 2 submatrices RZZt and RZZ⌢t.
The first one can be obtained using the well-known Lyapunov differential matrix covariance equation:(34)R˙ZZ=ARZZ+RZZAT+BZT+BZTT,where the last two matrices, in the case of modulated white noise input, can be expressed as(35)BZT+BZTT=Pt=2πS0φ2t0001.
Due to the symmetry of RZZ, it can be replaced by a vectorial form by the three first-order differential equations(36a)λ˙1=2λ2,(36b)λ˙2=−ω12λ1−2ξωλ2+λ3,(36c)λ˙3=−2ωλ2−4ξωλ3+2πS0φ2t,where λ1=σX2t, λ2=γXX˙t, and λ3=σX˙2t.
5. Hilbert Transform Space State Covariance Response
By applying a similar approach, it is possible to obtain the matrix RZZ˜t=Z¯Z˜¯T.
By making its first-order temporal derivation, it can be written as(37)ddtRZZ⌢t=R˙ZZ⌢t=Y˙tY˜tT+YtY˙⌢tT.
Moreover, by adopting the two-state vector differential equation [6], the left quantities can be expressed as(38a)Z˙Z⌢T=AZZ⌢T+BZ⌢T,(38b)ZZ˙⌢T=ZZ⌢TAT+ZB⌢T.
So finally obtaining the following 1st order differential equation matrix:(39)R˙ZZ⌢t=ARZZ⌢t+RZZ⌢tAT+BY˜T+YB⌢T.
Since the two last matrices are asymmetrical, it results as(40)BY˜T+YB⌢T=P˜t=0bt−bt0,where bt is the convolution integral:(41)bt=2Soφt∫0thuuφt−udu.
Since RZZ⌢t is asymmetrical, the following single first-order differential equation is obtained:(42)γ˙XX˙⌢t=−2ξωγXX˙⌢t+bt.whose solution depends on the instantaneous value of bt. Its evaluation can be quite expensive from a computational point of view, due to the involvement of a convolution integral; in fact, it has to be reevaluated at each instant t of the considered time interval. This problem can be solved if the modulation function φt is a separable one. This means that ϕt−u can be expressed as a summation of finite functions just depending on t or u:(43)φt−u=∑k=1mBktCku.
In this case, bt assumes the following form:(44)bt=2Soφt∑r=1mNrt=2Soφt∑r=1mBrtIrt,where Nrt=BrtIrt while Irt=∫0thρ/ρCrτdρ4 is a standard integral.
Equation (42) can be so rewritten as:(45)γ˙XX˙⌢t=−2ξωγXX˙⌢t+2Soφt∑r=1mBrtIrt
The covariance matrix RΨ¯Ψ¯t is then completely evaluated, by taking into account only the four evolutionary parameters λ1=σX2t, λ2=γXX˙t, λ3=σX˙2t, and λ4=γXX˙⌢t, and by using the following four first-order differential equations:(46a)λ˙1=2λ2,(46b)λ˙2=−ω12λ1−2ξωλ2+λ3,(46c)λ˙3=−2ωλ2−4ξωλ3+2πS0φ2t,(46d)λ˙4=−2ξωλ4+2Soφt∑r=1pArtIrt,where only the first three are coupled.
6. Envelope JPDF Evaluation
The covariance matrix RΨ¯Ψ¯t is equal to the matrix RUUt, corresponding to the vector U¯t=Xt,X˙t,X⌢t,X˙⌢tT,(47)RUU=U¯U¯T=σX2tγXX˙t0γXX˙⌢tγXX˙tσX˙2t−γXX˙⌢t00−γXX˙⌢tσX2tγXX˙tγXX˙⌢t0γXX˙tσX˙2t.whose JPDF is:(48)pU¯x,x˙,x⌢,x˙⌢;t=14π2RΨ¯Ψ¯t1/2exp−12U¯TtRΨ¯Ψ¯−1tU¯t.
As demonstrated before, a central point of Vanmarcke formulation for hazard function evaluation is the knowledge of the preenvelope process mean threshold crossing rate υA+β,t. Its determination requires the knowledge of the JPDF paa˙a,a˙;t, as follows:(49)υA+β,t=∫0+∞a˙pAA˙β,a˙;tda˙.
By applying the definition of the preenvelope process through the original process Xt and its nonstationary Hilbert transform X⌢t, it is useful to introduce the 4th order evolutionary zero mean random vector U¯t, whose joint PDF can be written as(50)pU¯x,x˙,x⌢,x˙⌢=14π2RUUexp−12U¯TRUU−1U¯.
By introducing the following vector:(51)V¯t=AtA˙tθtθ˙tT,where θ=tan−1X⌢t/Xt, it is possible to obtain the JPDF pV¯v¯,t [23] by using the well-known rule pU¯a,a˙,θ,θ˙;t=JpU¯x,x˙,X⌢,x˙⌢;t, J=∂U¯/∂V¯=a2 being the norm of the transformation. The joint density of V¯t is thus given by(52)pV¯a,a˙,θ,θ˙;t=aσxσx˙24π2ΔA2exp−aσxσx˙22ΔA2·a2σx˙2+a2θ˙2σx2+a˙2σx2−2aa˙γxx˙−2a˙2θ˙γxx˙⌢,where ΔAt=1−ρXX˙2t−ρXX˙⌢2t. The phase θ does not appear in (52), and then having a range over 0,2π, it is independent from the other three vector components, characterized by a uniform marginal distribution over the same interval. By applying the saturation with regard to both θ and θ˙, the required second-order JPDF can be written as(53)pAA˙a,a˙;t=∫02π∫−∞+∞pV¯a,a˙,θ,θ˙;tdθdθ˙=aBatexp−12a2Cat−2aa˙Dat+a˙2Eat,where Ba, Ca, Da, and Fa are reported in Appendix A as obtained in [25]. It is hence straightforward to obtain the marginal distribution of a and a˙ making saturation with respect to the other variable at each time, as expressed by(54a)pAa;t=aσx2texp−a22σx2t,(54b)pA˙a˙;t=a˙2πqt+ρxx˙2tσx˙texp−a˙2qtσx˙2.
By using the preenvelope density function pAA˙a,a˙;t, it is possible to establish the analytical form of the mean threshold crossing rate by solving the following integral:(55)vA+β,t=ηt2ΔAtσX˙tσXtexp−12η2tχdAt,where the adimensional threshold ηt is defined as the evolutionary function ηt=β/σXt. The function χy is given by χy=exp−y2/2+2πyΦy, where Φy is the normal cumulative function Φy=1/2π∫−∞yexp−u2/2du, and dAt=1/2ηtρXX˙⌢t/ΔAt.
By extending the original approach proposed by Vanmarcke [16, 17] and by using the nonstationary mean threshold crossing rate for the Gaussian process Xt in the Poisson hypothesis, it results:(56)vX+β,t=1πΔXtσX˙tσXtexp−12η2tχdXt,where ΔXt=1−ρXX˙2t and dXt=1/2ηtρXX˙t/ΔXt. With reference to the formulation from Michelov et al. [23](57)vVLM+ξ,t=1π1−ρXX˙2σX˙tσXt1−expvA+ξ,t/2vX+ξ,texpη2t/2−1.
It can be rewritten as:(58)vVLM+ξ,t=1πΔXtσX˙tσXtχdXt1−exp−π/2ΔAt/ΔXtχdAt/χdXtexpη2t/2−χdXt,which is the same relation proposed in [23] except for χdXt whose complete expression is(59)χϑ,t=e−ϑ2t+πϑtΦϑt2,where(60)ϑt=ηtρXX˙t1−ρXX˙2t,
Expression in (60) is function of time through both the dimensionless threshold value ηt and the ratio ρXX˙t/1−ρXX˙2t. Moreover, since ϑt can assume just positive values (in stationary case, it is exactly the dimensionless threshold η), its values are close to the unit, that is, ϑt is smaller than 1, growing up for greater values (Figure 1). In particular, in practical cases, it is quite close to one, as shown in the numerical applications.
Factor χ.
This formulation can be improved by introducing the parameter k, originally proposed by Vanmarcke [16, 17] for the stationary case, to enhance the agreement between the numerical simulation and the analytical formulation. This parameter is an empirical exponent of the “equivalent bandwidth factor” qXs defined for the stationary case. Its value has been evaluated to be equal to 1.2 and has been subsequently used in different studies.
In the nonstationary case, the parameter k has to be taken into account as the exponent of(61)qt=ΔAt/ΔXt=1−γXX˙⌢2t1−γXX˙2t.which coincides with the stationary value when ρXX˙2t=0.
The mean threshold crossing rate vVLM+ξ,t can be finally rewritten as(62)vVLM+ξ,t=1πΔXtσX˙tσXtχdXt1−exp−π/2qktχdAt/χdXtexpη2t/2−e−ϑ2t+πϑtΦϑt/2.
7. Numerical Results
Different modulation functions have been proposed in the scientific literature to model the behavior of nonstationary loads in a suitable way [5]. These functions typically depend on two or three control parameters. Aiming at achieving an accurate reliability evaluation, the present work uses a single parameter parabolic modulation function, whose equation is(63)φpt=−ttm2+2ttm.for a time interval equal to 0,2tm. The control of the amplitude modulation variation velocity is done by setting the time tm, which is the time at which the modulation reaches its maximum value, to be equal to one, as shown in Figure 2.
Parabolic modulation function for different values of tm.
Regarding the specific modulation function in (61), the nonstationary structural responses for an SDOF system have been attained by integrating the stationary white noise simulations, which are generated and then modulated.
Structural characteristics differ for damping ratio (ξ0=0.01−0.05−0.10), while the modulation control parameter is expressed as the dimensionless ratio between the time of maximum modulation and the natural structural system period T0:(64)τm=tmT0.with τ=t/T0. The values of τm used in this study are 2.5, 5, 10, and 15.
In addition to evaluating the structural response in the “classical” space state Z¯t (25) by the integration of the generated φtwt, the determination of this response is carried out in the “exact” and “approximate” space state Z^¯t and Z˜¯t (26), by integrating Hφtwt and φtHwt, respectively. The Hilbert transforms of the generated stationary white noise and the modulated one are accomplished by MATLAB standard algorithms. Figures 3 and 4 show displacement results and their envelopes evaluated by “exact” and “approximate” ways in the cases of damping ratio equal to 0.01 and 0.10. It is quite evident that there is an acceptable agreement between the two different envelopes, except for a few instants at the beginning of the excitation. This is an effect exactly of the difference of the proposed methodology that has a zero value at the starting time. This consideration, which has originally been suggested by Borino et al. [31], is not exactly correct when referring to the first time derivative of the envelope, as illustrated in Figures 5 and 6. In these two figures, it can be observed that some discrepancies exist between the approximate and exact nonstationary solutions, and they increase as the damping ratio diminishes. In this case, it is evident that at the end of the process, there are discrepancies, especially for case 2 (Tm = 10T0) and for low damping. In this specific case, there is a separation between the two results at the end when the approximate solution has a strongly oscillatory shape. This result is inferred more clearly in Figures 7–10 (each for different values of damping ratio: 0.01, 0.05, and 0.10), where the JPDF of At and A˙t is shown for different values of τ/τm=0.5−1.0−1.5.
Structural displacement Xt and its envelope At evaluated by “exact” and “approximate” ways, for damping ratio equal to 0.01.
Structural displacement Xt and its envelope At evaluated by “exact” and “approximate” ways, for damping ratio equal to 0.1.
Approximate and exact nonstationary solutions for Tmax = 2.5T0.
Approximate and exact nonstationary solutions for Tmax = 10T0.
JPDF of At and A˙t for ξ = 0.1 and Tmax = 2.5T0.
JPDF of At and A˙t for ξ = 0.01 and Tmax = 2.5T0.
JPDF of At and A˙t for ξ = 0.1 and Tmax = 10T0.
JPDF of At and A˙t for ξ = 0.01 and Tmax = 10T0.
Furthermore, the good accordance between the approximate and exact numerical evaluations of pAa;t is obvious in all cases. However, this agreement does not hold for the PDF of A˙t, where the two cases tend to diverge especially for high values of the envelope velocity. This can be attributed reasonably to numerical reasons. This condition (differences between pA˙a˙;t evaluated by the two different approaches) can lead to a nonperfect agreement between the approximate and exact approaches for the mean crossing rate of the envelope process, which depends on the JPDF of At and A˙t. Additionally, and as stated above, the adoption of the approximate solution in (22b) instead of those in (19) leads to a direct and serious reduction of the analytical and numerical efforts in covariance “preenvelope” evaluations and thus in reliability.
For a better evaluation of differences for threshold crossing problem deriving by the various approaches, in Figures 11 and 12, the numerical results for conditioned crossing rates are shown for different values of the barrier. The barriers are in a range between 1 and 4 in Figure 11 (short phenomena) and 1 and 3 in Figure 12 (long phenomena). In both cases, there is a reasonable agreement with the proposed method, which is anyway quite close to the Michelov et al. [19] formulation, if an exponent k equal to one is adopted. Using a value equal to 1.2 gives some results just a little bit smaller, but the difference is usually negligible. It should also be noticed that the Michelov and the proposed formulations are in a good agreement with numerical results in the first part of the phenomenon, during the growing phase, but generally underestimate the real threshold crossing rate in the decay part.
Mean threshold crossing rate for ξ = 0.01 and Tmax = 5T0.
Mean threshold crossing rate for ξ = 0. 1 and Tmax = 10T0.
Finally, in order to evaluate the empirical correction exponent k in (62), numerical and analytical results are compared for different values of k in Figure 13 (it is commonly assumed as 1.2 in stationary case, as originally proposed by Vanmarcke).
Numerical and analytical results in function of k.
8. Conclusions
This study dealt with the first crossing reliability evaluation of structures subject to nonstationary Gaussian random input. The solution based on the qualified envelope, originally derived for stationary cases, has been extended to the nonstationary case. At this aim, the modulation of stationary envelope has been adopted by utilizing the so-called “preenvelope” process. A modified reliability formulation has also been proposed, based on the extension of the formulation of the empirical bandwidth factor exponent developed under the stationary hypothesis.
To evaluate the accuracy of the reliability evaluation achieved by using the previous assumption, the mean crossing rate and the envelope JPDF obtained by numerical simulations for a SDOF system subject to Gaussian white noise have been compared to the ones obtained in analytical way. Different damping ratios and velocities in modulation amplitude variation have been taken into account. Results show that only a partial agreement is obtained, based not directly on the final reliability, but also on a suitable measure of accuracy of hazard function. These analyses are carried out on a wide number of cases taking into account different damping factors and nonstationary levels.
Aiming at performing a parametric analysis on nonstationary input characterization, in the case of load conditions such as earthquakes [32, 33] a modulation function controlled by a single parameter is used as follows:(B.1)φt=−ttm2+2ttm,where tm is the maximum intensity time for which ϕtm=1. This function is symmetric, and then its total duration is 2tm. In this case, it results:(B.2)ϕt−τ=−t¯2+2t¯−τ¯¯2−2τ¯¯+2t¯τ¯,where t¯=t/tm;τ¯=τ/tm.
Therefore, it is possible to rewrite bt in the following form:(B.3)bt=2So−ttm2+2ttm∫0thττ−t¯2+2t¯−τ¯¯2−2τ¯¯1+τ¯¯dτ=I1t+I2t+I3t,where(B.4)I1t=2So−ttm2+2ttm2∫0te−ξωuusinω¯udu=2Sotm−t2+2ttmΨt,ω¯,−ξω,I2t=2Sotm2−ttm2+2ttm∫0tτe−ξωτsinω¯τdτ=2Sotm4−t2+2ttm1−2ξ2+e−ξωτ−ξωt−2ξ2+1cosω¯t+ωt−2ξ1−ξ2sinω¯t,I3t=2Sotm−ttm2+2ttm1+ttm∫0te−ξωτsinω¯τdτ=2Sotm4−t2+2ttmtm−11ω1−ξ2−e−ξωτξsinω¯t+1−ξ2cosω¯t.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
AtesS.Seismic behaviour of isolated multi-span continuous bridge to nonstationary random seismic excitation201267126328210.1007/s11071-011-9976-72-s2.0-82255164053Hong-YuJ.De-YiZ.Shi-XiongZ.Wei-ChauX.MaheshD. P.Local site effects on a high-pier railway bridge under tridirectional spatial excitations: nonstationary stochastic analysis2013525569JinX. L.HuangZ. L.Nonstationary probability densities of strongly nonlinear single-degree-of-freedom oscillators with time delay2010591-219520610.1007/s11071-009-9532-x2-s2.0-77949424568JunH. E.An approximation of the first passage probability of systems under nonstationary random excitation200930225526210.1007/s10483-009-0213-y2-s2.0-59449092534GrecoR.FioreA.MaranoG. C.The role of modulation function in nonstationary stochastic earthquake model201485145001510.1142/s17934311145001582-s2.0-84927722302MaranoG. C.AccianiG.FioreA.AbresciaA.Integration algorithm for covariance non-stationary dynamic analysis of SDOF systems using equivalent stochastic linearization2015152145004410.1142/s02194554145004482-s2.0-84914132256GrecoR.FioreA.BriseghellaB.Influence of soil type on damping reduction factor: a stochastic analysis based on peak theory201810436536810.1016/j.soildyn.2017.10.0202-s2.0-85034111831MaranoG. C.TrentadueF.GrecoR.Optimum design criteria for elastic structures subject to random dynamic loads200638785387110.1080/030521506009130282-s2.0-33749456022GrecoR.MaranoG. C.Robust optimization of base isolation devices under uncertain parameters201622385386810.1177/10775463145326702-s2.0-84954287764MaranoG. C.TrentadueF.GrecoR.Stochastic optimum design criterion of added viscous dampers for buildings seismic protection2007251213710.12989/sem.2007.25.1.0212-s2.0-33845432691NigamN. C.1983Cambridge, MA, USAMIT PressLutesL. D.SarkaniS.2004Oxford, UKElsevier Butterworth-HeinemannBucherC.Di PaolaM.Efficient solution of the first passage problem by Path Integration for normal and Poissonian white noise20154112112810.1016/j.probengmech.2015.06.0072-s2.0-84936950216KougioumtzoglouI. A.SpanosP. T. D.Response and first-passage statistics of nonlinear oscillators via a numerical path integral approach201813991207121710.1061/(asce)em.1943-7889.0000564ValdebenitoM. A.JensenH. A.LabarcaA. A.Estimation of first excursion probabilities for uncertain stochastic linear systems subject to Gaussian load2014138364810.1016/j.compstruc.2014.02.0102-s2.0-84897086004VanmarckeE. H.On the distribution of the first passage time for normal stationary random process197542121522010.1115/1.34235212-s2.0-0016483905VanmarckeE. H.Properties of spectral moments with applications to random vibrations1972982425446Di PaolaM.Transient spectral moments of linear systems198510225243MichelovG.SarkanyS.LutesL. D.Spectral characteristics of nonstationary random process-a critical review199921322324410.1016/s0167-4730(99)00022-32-s2.0-0033356639CorotisR. B.VanmarckeE. H.CornellA. C.First passage of nonstationary random processes1972982401414DugundjiJ.Envelope and pre-envelope of real waveforms195841535710.1109/tit.1958.10574352-s2.0-0002674160LangleyR. S.On various definition of the envelope of a random process1986105350351210.1016/0022-460x(86)90175-62-s2.0-0022677458MichelovG.SarkanyS.LutesL. D.Spectral characteristics of nonstationary random process-response of a simple oscillator199921324526710.1016/s0167-4730(99)00019-32-s2.0-0033328515LangleyR. S.Structural response to non-stationary non-white stochastic ground motion198614690992410.1002/eqe.42901406072-s2.0-0022811128MichelovG.LutesL. D.SarkanyS.Extreme value of response to nonstationary excitation20011272135236310.1061/(asce)0733-9399(2001)127:4(352)2-s2.0-0035310726YangJ. N.Non stationary envelope process and first excursion probability19721223124810.1080/036012172089053412-s2.0-84963482911CrandallS. H.MarkW. D.1963New York, NY, USAAcademic PressPapoulisA.1984New York, NY, USAMcGraw-HillDi PaolaM.MuscolinoG.On the convergent parts of spectral moments1986110223324510.1016/s0022-460x(86)80207-32-s2.0-0022791595MuscolinoG.Nonstationary envelope in random vibration theory198811481396141310.1061/(asce)0733-9399(1988)114:8(1396)2-s2.0-0024067351BorinoG.Di PaolaM.MuscolinoG.Non stationary spectral moments of base excited MDOF systems198816574575610.1002/eqe.42901605092-s2.0-0024053529NutiC.RasuloA.VanziI.Seismic assessment of utility systems: Application to water, electric power and transportation networks3Proceedings of the Joint ESREL and SRA-Europe Conference-Safety, Reliability and Risk analysis: Theory, Methods and Applications200925192529LavoratoD.VanziI.NutiC.MontiG.Generation of non-synchronous earthquake signals2017201716919810.1007/978-3-319-52425-2_82-s2.0-85014075385