12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 1 An Approximate Approach for Assessing the Reliability of a Stochastically Excited Softening Duffing Oscillator Yuanjin Zhang Institute for Risk & Uncertainty, University of Liverpool, Liverpool, UK Ioannis A. Kougioumtzoglou Assistant Professor, Dept. of Civil Engineering and Engineering Mechanics, Columbia University, New York, USA ABSTRACT: An approximate analytical technique for assessing the reliability of a softening Duffing oscillator subject to evolutionary stochastic excitation is developed. Specifically, relying on a stochastic averaging treatment of the problem the oscillator time-varying survival probability is determined in a computationally efficient manner. In comparison with previous techniques that neglect the potential unbounded response behavior of the oscillator when the restoring force acquires negative values, the herein developed technique readily takes this aspect into account by introducing a special form for the oscillator non-stationary response amplitude probability density function (PDF). A significant advantage of the technique relates to the fact that it can readily handle cases of stochastic excitations that exhibit strong variability in both the intensity and the frequency content. An illustrative numerical example of a softening Duffing oscillator subject to earthquake excitation is included. Comparisons with pertinent Monte Carlo simulation data demonstrate the efficiency of the technique. 1. INTRODUCTION It is often desirable for risk assessment applications to estimate the probability (also known as survival probability) that the system response stays within a prescribed domain over a given time interval. Indicatively, advanced Monte Carlo simulation (MCS) methodologies such as importance sampling, subset simulation and line sampling have been developed for reliability assessment applications; see Bucher (2011), and Au and Beck (2001) for some indicative references. Note, however, that there are cases of complex systems where MCS can be a computationally demanding, or even a prohibitive task; thus, there is a need for developing efficient approximate analytical and/or numerical techniques for addressing the problem such as Poisson distribution based approximations (e.g. Vanmarcke (1975)), probability density evolution schemes (e.g. Li and Chen (2009)), stochastic averaging/linearization approaches (e.g. Spanos and Kougioumtzoglou (2014a, 2014b)), as well as Wiener path integral techniques (e.g. Kougioumtzoglou and Spanos (2014a), Zhang and Kougioumtzoglou (2014)) and numerical path integration schemes (e.g. Naess and Johnsen (1993), Di Paola and Santoro (2008)). The softening Duffing oscillator is a nonlinear oscillator possessing a linear-plus-cubic restoring force so that the spring has a softening characteristic. This oscillator has received considerable attention in the literature primarily due to its importance in describing the roll motion of a ship model in beam seas (e.g. Spyrou and Thomson (2000), Belenky and Sevastianov (2007)). Note, however, that the softening Duffing oscillator has found applications in diverse other fields of engineering dynamics such as structural system vibration isolation (e.g. Fu et al. (2014)), energy harvesting (e.g. Vandewater and Moss (2013)) and dynamics of timber structures (e.g. Reynolds et al. (2014)). 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 2 Further, although several research efforts have focused on studying the oscillator response under deterministic excitation (e.g. Nayfeh and Sanchez (1989), Brennan et al. (2008)), limited results exist regarding the response analysis of the oscillator when it is subjected to stochastic excitation (e.g. Roberts and Vasta (2000), Cottone et al. (2010)). Specifically, most of the results are based on rather heuristic approaches which inherently assume stationarity and that the probability the response leaves the stable region is extremely small; thus, neglecting important aspects of the analysis such as the possible unbounded response behavior when the restoring force acquires negative values. Recently, a numerical path integral approach was developed in Kougioumtzoglou and Spanos (2014b) for determining the survival probability of a softening Duffing oscillator subject to stochastic excitation. The unbounded character of the response was rigorously taken into account by introducing a special form for the conditional response PDF, while the solution was propagated by utilizing a discrete version of the C-K equation. In this paper, an efficient approximate analytical technique for determining the survival probability of a softening Duffing oscillator subject to evolutionary stochastic excitation is developed. Specifically, relying on a stochastic averaging treatment of the problem and introducing a special form for the oscillator response PDF, the technique developed in Spanos and Kougioumtzoglou (2014b) is adapted and generalized herein to account for the special case of the softening Duffing oscillator. A significant advantage of the technique is that it can readily handle cases of evolutionary stochastic excitation with arbitrary evolutionary power spectrum (EPS) forms, even of the non-separable kind. 2. MATHEMATICAL FORMULATION 2.1. Softening Duffing oscillator response analysis Consider the softening Duffing oscillator whose motion is governed by the equation ?Μ? + 2π0π0?Μ? + π02π₯ + ππ02π₯3 = π€(π‘), (1) where a dot over a variable denotes differentiation with respect to time π‘ ; π < 0 denotes a negative constant representing the magnitude of the nonlinearity degree; π0 is the ratio of critical damping; π0 is the natural frequency corresponding to the linear oscillator (i.e. π = 0) and π€(π‘) represents a Gaussian, zero-mean non-stationary stochastic process possessing an evolutionary broad-band power spectrum ππ€(π, π‘) . Examining Eq.(1), it can be readily seen that there exist values of the response displacement π₯(π‘) for which the oscillator restoring force πΉ(π₯) = π02π₯ + ππ02π₯3 =π02π₯(1 + ππ₯2) reaches zero, and even negative values. Clearly, this may lead to unbounded system response, and a special treatment is necessary to account for this behavior. Next, bearing this qualitative behavior in mind, and focusing on lightly damped systems (i.e. π0 βͺ 1 ), it can be argued (e.g. Spanos and Lutes (1980)) that for πΉ(π₯) = π02π₯(1 + ππ₯2) β₯ 0 , or equivalently π₯2 β₯ β1/π, the oscillator response exhibits a pseudo-harmonic behavior described by the equations π₯(π‘) = πcos[π(π)π‘ + π(π‘)], (2) and ?Μ?(π‘) = βπ(π) πsin[π(π)π‘ + π(π‘)]. (3) In Eqs.(2-3), π and π represent a slowly varying with time phase and a slowly varying with time response amplitude, respectively. Manipulating Eqs.(2-3) yields an expression for the oscillator response amplitude; that is, π(π‘) = βπ₯2(π‘) +?Μ?2(π‘)π(π). (4) It is primarily the assumption of light damping that allows a combination of deterministic and stochastic averaging to be performed next and to approximate the second-order stochastic differential equation 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 3 (SDE) (Eq.(1)) by a first-order SDE governing the response amplitude process π. A more detailed presentation/discussion of the assumptions involved and the corresponding assumed pseudo-harmonic behavior of the response process π₯(π‘) can be found in references such as Spanos and Lutes (1980) , Roberts and Spanos (1986), and Kougioumtzoglou and Spanos (2009). Next, following a stochastic averaging/linearization approach (e.g. Kougioumtzoglou and Spanos (2009)) a linearized version of Eq.(1) becomes ?Μ? + 2π0π0?Μ? + π2(π)π₯ = π€(π‘), (5) where the equivalent natural frequency π(π) is given by the expression π2(π) =π02ππβ« πππ π((π cosπ2π0+ π(π cosπ)3)ππ= π02 (1 +34ππ2). (6) Examining Eq.(6) it can be readily seen that the stiffness element of the equivalent linear oscillator becomes zero at the critical response amplitude value πππ = ββ43π. In this regard, the requirement π₯2 β₯ β1/π for the oscillator of Eq.(1) to have a bounded response is equivalently expressed in the following by the requirement π < πππ . Bearing this qualitative aspect in mind, a special form for the non-stationary response amplitude PDF π(π, π‘) is introduced next; that is, π(π, π‘) =ππ(π‘)exp (βπ22π(π‘)) ππππ‘(π)+ π(π‘)πΏ(π β πβ), (7) where ππππ‘(π) = π’(π) β π’(π β πππ) , π’(. ) denotes the unit step function, π(π‘) is a time-dependent coefficient to be determined, πΏ(. ) denotes the Dirac delta function, and πβ represents an arbitrary response amplitude value with the property πβ β« π β [0, πππ] . Further, the time-dependent factor π(π‘) can be determined by applying the normalization condition β« π(π, π‘)ππβ0= 1; this yields π(π‘) = 1 β β«ππ(π‘)exp (βπ22π(π‘))πππππ0= exp (βπππ22π(π‘)). (8) Examining the form of the non-stationary response amplitude PDF of Eq.(7), it can be readily seen that it comprises two conceptually different terms. The first one represents a truncated Rayleigh PDF for amplitude values in the range [0, πππ] , whereas the factor π(π‘) in the second term represents the probability at a specific time instant that the response grows unbounded, namely the system response asymptotically approaches infinity. The rationale behind the choice of the truncated time-dependent Rayleigh PDF of Eq.(7) relates to the fact that the linear oscillator stationary response amplitude PDF is a Rayleigh one. In fact, as it was shown in Spanos and Lutes (1980), the non-stationary response amplitude PDF of a linear oscillator subject to Gaussian white noise excitation is a time-dependent Rayleigh PDF of the form π(π, π‘) =ππ(π‘)exp (βπ22π(π‘)) with the property limπ‘ββ π(π, π‘) =ππ2ππ₯π (βπ22π2) ; where π2 represents the linear oscillator stationary response variance. In Kougioumtzoglou and Spanos (2009), it was further shown that the Rayleigh representation is suitable for nonlinear oscillators also and under evolutionary stochastic excitation as well. It is pointed out that a significant difference between adopting a PDF of the form π(π, π‘) =ππ(π‘)exp (βπ22π(π‘)) in Kougioumtzoglou and Spanos (2009) and introducing a PDF form of Eq.(7) in the herein developed technique, is that in the former case π(π‘) accounts for the variance of the non-stationary response process π₯ , whereas in the latter case π(π‘) is simply a time-varying coefficient to be determined. 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 4 Next, relying on Eq.(7), it can be argued that an alternative to Eq.(5) equivalent linear system is given in the form ?Μ? + 2π0π0?Μ? + πππ2 (π‘)π₯ = π€(π‘), (9) where the time-dependent stiffness element πππ2 (π‘) is defined as (see also Spanos and Kougioumtzoglou (2014b), Kougioumtzoglou and Spanos (2009)) πππ2 (π‘) = πΈ[π2(π)] = β« π2(π)π(π, π‘)ππβ0. (10) Note that taking into account the form of the amplitude PDF of Eq.(7), the time-varying equivalent stiffness element of Eq.(10) also has two parts. Specifically, for π β [0, πππ] , πππ2 (π‘) has a bounded part, i.e. πππ,π΅2 (π‘) , whereas for π > πππ the stiffness element πππ2 (π‘) exhibits negative values; thus, yielding negative restoring force values resulting potentially in an unbounded system response behavior. In this regard, utilizing Eq.(7) the bounded part πππ,π΅2 (π‘) is determined as πππ,π΅2 (π‘) = β« π2(π)π(π, π‘)πππππ0. (11) Analytical determination of the integral in Eq.(11) yields πππ,π΅2 (π‘) = π02 (1 +32ππ(π‘)(1 β π(π‘))). (12) Examining Eq.(12) it can be readily seen that the stiffness element πππ,π΅2 (π‘) is bounded between the values 0 and π02 . Specifically, assuming that the oscillator is initially at rest yields limπ‘β0+ π(π, π‘) = πΏ(π0) , or in other words, limπ‘β0+ π(π‘) = 0 , which yields limπ(π‘)β0+ πππ,π΅2 (π‘) = π02 . This means that for the very early part of the oscillation duration the oscillator features an approximately linear restoring force. Further, as time increases and the transient phase progresses, the truncated Rayleigh PDF of Eq.(7) broadens as the oscillator exhibits higher amplitude values π(π‘) . Equivalently, the time-varying coefficient π(π‘) increases with time, whereas the equivalent stiffness part πππ,π΅2 (π‘) decreases with time. Taking into account Eqs.(7) and (12) it can be readily shown that in the extreme case limπ(π‘)ββπππ,π΅2 (π‘) = 0. Thus, the equivalent stiffness part πππ,π΅2 (π‘) is a non-negative and bounded quantity varying with time between the values 0 and π02 . This is in agreement with the fact that πππ,π΅2 (π‘) corresponds to amplitude values π β [0, πππ] where the oscillator response is assumed to behave in a bounded manner. Further, focusing on the case where π β [0, πππ] and based on a stochastic averaging approach Eq.(9) can be cast in a first-order SDE governing the evolution in time of the amplitude π(π‘); see Spanos and Lutes (1980) , Roberts and Spanos (1986), and Kougioumtzoglou and Spanos (2009) for a more detailed presentation. Related to this SDE is the Fokker-Planck (F-P) partial differential equation ππ(π, π‘|π1, π‘1)ππ‘= βπππ[πΎ1(π, π‘)π]+12π2ππ2[πΎ2(π, π‘)π], (13) where πΎ1(π, π‘) = βπ0π0π +ππ(πππ,π΅(π‘), π‘)2ππππ,π΅2 (π‘) , (14) and πΎ2(π, π‘) =ππ(πππ,π΅(π‘), π‘)πππ,π΅2 (π‘). (15) The F-P Eq.(13) governs the evolution in time of the transition PDF π(π, π‘|π1, π‘1) for π β [0, πππ] and π1 β [0, πππ] . Next, a solution of the associated F-P equation π(π, π‘|π1 = 0, π‘1 = 0) = π(π, π‘) is attempted in the form of the truncated Rayleigh PDF of Eq.(7). Specifically, substituting the truncated Rayleigh PDF into the associated F-P equation, assuming that the oscillator is initially at rest (i.e. π(π, π‘ = 0) = πΏ(π)), and manipulating yields the first-order nonlinear differential equation ?Μ?(π‘) = β2π0π0π(π‘) +ππ(πππ,π΅(π‘), π‘)πππ,π΅2 (π‘), (16) 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 5 to be solved numerically for the time-varying coefficient π(π‘) . Obviously, once the time-varying coefficient π(π‘) is determined, the time-dependent coefficient π(π‘) can be evaluated via Eq.(8). Further, equations similar to Eq.(16) can be derived for the case of the response amplitude transition PDF in a straightforward manner. Specifically, following a similar analysis as in Spanos and Solomos (1983), the transition amplitude PDF π(π, π‘|π1, π‘1) is sought in the form π(π, π‘|π1, π‘1)= {ππ‘π(π, π‘|π1, π‘1) + π (π‘, π‘1)πΏ(π β πβ) , π1 β (0, πππ)πΏ(π β πβ), π1 > πππ, (17) where ππ‘π(π, π‘|π1, π‘1) =ππ(π‘, π‘1)exp (βπ2 + β2(π‘, π‘1)2π(π‘, π‘1)) πΌ0 (πβ(π‘, π‘1)π(π‘, π‘1)) ππππ‘(π), (18) and π(π‘, π‘1) and β(π‘, π‘1) are time-varying coefficients to be determined. Further, applying the normalization condition β« π(π, π‘|π1, π‘1)β0ππ = 1 yields the time-varying coefficient π (π‘, π‘1) = 1 ββ« ππ‘π(π, π‘|π1, π‘1)πππππ0, (19) where πΌ0(. ) denotes the modified Bessel function of the first kind and of zero order. In a similar manner as before, under the condition that π β [0, πππ] and π1 β [0, πππ] substituting the bounded part of Eq.(17) into Eq. (13) and manipulating yields the first-order differential equations (see Spanos and Solomos (1983) for a more detailed derivation) ππ(π‘, π‘1)ππ‘+ 2π0π0π(π‘, π‘1) βππ(πππ,π΅(π‘), π‘)πππ,π΅2 (π‘)= 0, (20) and πβ(π‘, π‘1)ππ‘+ π0π0β(π‘, π‘1) = 0. (21) Eqs.(20-21) are subject to the initial condition π(π2, π‘1|π1, π‘1) = πΏ(π2 β π1) which states that no change of state can occur if the transition time is zero. 2.2 Softening Duffing oscillator reliability assessment In this section the approach developed in Spanos and Kougioumtzoglou (2014b) is adapted and generalized herein to account for the special case of the softening Duffing oscillator and to determine the oscillator time-dependent survival probability. This is defined as the probability ππ΅(π‘) that the amplitude π stays below the threshold πππ over a given time interval [π‘0, π] ; that is, ππππ[π(π‘) β€ πππ , ππ£ππ [π‘0, π]|π(π‘0) < πππ]. In the following, adopting the dicretization scheme applied in Spanos and Kougioumtzoglou (2014b), the time domain is divided into intervals of the form [π‘πβ1, π‘π ], π = 1,2β¦ ,π, π‘0 = 0, π‘π = π πππ π‘π = π‘πβ1 + ππ πππ(π‘πβ1), (22) where πππ denotes the equivalent natural period of the oscillator given by πππ(π‘) =2ππππ,π΅(π‘) , (23) and ππ is a constant to be selected with the property ππ β (0,1]. In the ensuing analysis, the survival probability is determined assuming that it is approximately constant over the time interval [π‘πβ1, π‘π ]. Clearly, for ππ = 1 the time interval [π‘πβ1, π‘π] corresponds to the equivalent time-dependent natural period of the oscillator. The choice is justified by the fact that the response amplitude π is assumed to be approximately constant over the interval [π‘πβ1, π‘π] , owing to its slowly varying character with respect to time (see section 2.1). Thus, the survival probability ππ΅(π) is assumed to be constant over [π‘πβ1, π‘π] as well. Of course, if higher accuracy is required a smaller value for ππ can be chosen. Further, taking into account the discretization of Eq.(22), the survival probability ππ΅(π) is given by the equation ππ΅(π = π‘π) =β(1 β πΉπ)ππ=1, (24) 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 6 where πΉπ is defined as the probability that π will cross the barrier πππ in the time interval [π‘πβ1, π‘π] , given that no crossings have occurred prior to time π‘πβ1 . Next, invoking the Markovian property for the process π and utilizing the standard definition of conditional probability yields πΉπ =ππππ[π(π‘π) β₯ πππ βπ(π‘πβ1) β€ πππ]ππππ[π(π‘πβ1) β€ πππ]=ππβ1,ππ»πβ1(25) where π¨πβ1 = β« π(ππβ1, π‘πβ1)πππβ1πππ0, (26) and, by utilizing the relationship π(π1, π‘1; π2, π‘2) = π(π1, π‘1)π(π2, π‘2|π1, π‘1), ππβ1,π = β« (β« π(ππ , π‘π|ππβ1, π‘πβ1)πππ+βπππ)πππ0 π(ππβ1, π‘πβ1)πππβ1. (27) Next, taking into account Eqs.(7) and (17), Eqs.(26-27) become π»πβ1 = 1 β exp (πππ22π(π‘πβ1)), (28) and ππβ1,π= β« (β« (ππ‘π(ππ , π‘π|ππβ1, π‘πβ1)+βππππππ0+ π (π‘π, π‘πβ1)πΏ(π β πβ))πππ) π(ππβ1, π‘πβ1)πππβ1. (29) respectively. Taking into account the properties of the Dirac delta function, Eq.(29) becomes ππβ1,π = β« π (π‘π, π‘πβ1)πππ0π(ππβ1, π‘πβ1)πππβ1, (30) and utilizing Eq.(19) yields ππβ1,π = β« π(ππβ1, π‘πβ1)πππ0πππβ1 ββ« (β« ππ‘π(ππ , π‘π|ππβ1, π‘πβ1)ππππππ0)πππ0 π(ππβ1, π‘πβ1) πππβ1. (31) Next, considering Eqs.(26), Eq.(31) takes the form ππβ1,π = π»πβ1 ββ« β« ππ‘π(ππ , π‘π|ππβ1, π‘πβ1)πππ0πππ0 π(ππβ1, π‘πβ1)ππππππβ1. (32) Relying further on the assumption that πππ,π΅(π‘) follows a slowly varying with time behavior, the following approximation over a small time interval [π‘πβ1, π‘π] is introduced; i.e., πππ,π΅(π‘) = πππ,π΅(π‘πβ1) for π‘ β [π‘πβ1, π‘π] . Next, based on the slowly varying with time behavior of the EPS, ππ€(π, π‘) is also treated as a constant over the interval [π‘πβ1, π‘π] . Further, based on the above assumptions, introducing the variables ππ = π‘π β π‘πβ1 , ππ2 =π(π‘πβ1)π(π‘π)(1 β 2π0π0ππ), and applying a first-order Taylor expansion around point ππ = 0 for various quantities, analytical treatment of the involved double integral of Eq.(32) is possible yielding ππβ1,π = π»πβ1 β (π΄0 +βπ΄πππ=1), (33) where π΄0 = (1 β exp (βπππ22π(π‘π)(1 β ππ2))) (1 βexp (βπππ22π(π‘π)(1 β ππ2)))(1 β ππ2), (34) π΄π =ππ2π(1 β ππ2)(π!)2πΏπ , (35) and πΏπ = (Ξ[1 + π, 0] β Ξ [1 + π,πππ22π(π‘πβ1)(1 β ππ2)]) (Ξ[1 + π, 0] β Ξ [1 + π,πππ22π(π‘π)(1 β ππ2)]) . (36) In Eq.(36) Ξ[πΎ, π§] represents the incomplete Gamma function defined as Ξ[πΎ, π§] =β« π‘πΎβ1πβπ‘ππ‘+βπ§. A more detailed presentation of the derivations in this section can be found in Spanos and Kougioumtzoglou (2014b). 3. NUMERICAL EXAMPLES As noted in the introductory section, although the softening Duffing oscillator has been widely utilized to model the nonlinear ship rolling motion in beam seas (e.g. Spyrou and Thomson (2000), Belenky and Sevastianov (2007)), it has also been used in conjunction with structural dynamics/earthquake engineering applications such as structural system vibration isolation (e.g. Fu et al. (2014)), energy harvesting (e.g. Vandewater 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 7 and Moss (2013)) and dynamics of timber structures (e.g. Reynolds et al. (2014)). In this regard, the non-separable excitation EPS of the form ππ€(π, π‘)= π (π5π)2ππ₯π(β0.2π‘)π‘2ππ₯π (β (π10π)2π‘), (37) is considered in this example. This spectrum comprises some of the main characteristics of seismic shaking, such as decreasing of the dominant frequency with time (e.g. Sabetta and Pugliese (1996)). Further, survival probabilities determined via the herein developed approximate technique are compared with pertinent Monte Carlo simulation data (10,000 realizations). To this aim, realizations compatible with the EPS of Eq.(37) are generated based on a spectral representation approach (e.g. Liang et al. (2007)), while a standard fourth-order Runge-Kutta scheme is employed for solving the nonlinear equation of motion (Eq.(1)). The initial distribution chosen for the response amplitude PDF is the Dirac delta function, i.e., π(π0, π‘0 = 0) = πΏ(π0), assuming the system is initially at rest. In the ensuing analysis the value π = 60 is chosen for the terms to be included in the expansion. Fig.(1). Bounded equivalent time-varying natural frequency πππ,π΅(π‘) for a softening Duffing oscillator (π = 1, π02 = π2, π½0 = 0.0628) In Fig.(1), the bounded equivalent natural frequencies (Eq.(12)) of the oscillators with parameter values (π = 1,π02 = π2, π½0 =0.0628, π = β1), (π = 1,π02 = π2, π½0 =0.0628, π = β2), and (π = 1, π02 = π2, π½0 =0.0628, π = β3) are plotted. In Fig.(2) the survival probabilities determined by Eqs.(24) are plotted for various barrier levels πππ =ββ43π; comparisons with MCS (10000 realizations) demonstrate a quite satisfactory agreement. Fig.(2). Survival probability for a softening Duffing oscillator (π = 1, π02 = π2, π½0 = 0.0628, ππ = 0.125) ; comparisons with MCS (10,000 realizations) 4. CONCLUSION In this paper, an approximate analytical technique based on stochastic averaging has been developed for determining the survival probability of a softening Duffing oscillator subject to evolutionary stochastic excitation. In this regard, by introducing a special form for the non-stationary response amplitude PDF that takes into account the potential unbounded behavior of the oscillator, the time-varying survival probability has been determined in a computationally efficient manner. Approximate technique based results regarding a softening Duffing oscillator under earthquake excitation have been compared with pertinent MCS data demonstrating a satisfactory level of accuracy. 5. REFERENCES Au S.-K., Beck J.L., 2001, Estimation of small failure probabilities in high dimensions by subset simulation, Prob Eng Mech, 16 (4): 263β277. Belenky V. L., Sevastianov N. B., 2007. Stability and safety of ships: Risk of capsizing, The Society of Naval Architects and Marine Engineers, 2nd Edition. Brennan M. J., Kovacic I., Carrella A., Waters T. P., 2008. On the jump-up and jump-down frequencies of the Duffing oscillator, Journal of Sound and Vibration, vol. 318: 1250-1261. Bucher C., 2011. Simulation methods in structural reliability, Marine Technology and Engineering, vol. 2: 1071-1086. 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 8 Cottone G., Di Paola M., Ibrahim R., Pirrotta A., Santoro R., 2010. Stochastic ship roll motion via path integral method, International Journal of Naval Architecture and Ocean Engineering, vol. 2: 119-126. Di Paola, M., and Santoro, R. 2008. Path integral solution for nonlinear system enforced by Poisson white noise. Probab. Eng. Mech., 23(2β3), 164β169. Fu Niu, et al, 2014, Design and analysis of a quasi-zero stiffness isolator using a slotted conical disk spring as negative stiffness structure, Journal of Vibroengineering, 16(4): 1392-8716. Kougioumtzoglou I. A., Spanos P. D., 2009, An approximate approach for nonlinear system response determination under evolutionary stochastic excitation, Current Science, 97:1203-1211. Kougioumtzoglou I. A., Spanos P. D., 2014a, Nonstationary Stochastic Response Determination of Nonlinear Systems: A Wiener Path Integral Formalism, J. Eng. Mech., 04014064-1~ 04014064-14. Kougioumtzoglou I. A., Spanos P. D., 2014b. Stochastic response analysis of the softening Duffing oscillator and ship capsizing probability determination via a path integral approach, Probabilistic Engineering Mechanics, vol. 35: 67-74. Li J., Chen J., 2009. Stochastic dynamics of structures, New York: John Wiley & Sons. Liang J., Chaudhuri S. R., Shinozuka M., 2007. Simulation of non-stationary stochastic processes by spectral representation, Journal of Engineering Mechanics, vol. 133: 616-627. Naess, A., and Johnsen, J. M. 1993, Response statistics of nonlinear, compliant offshore structures by the path integral solution method. Probab. Eng. Mech., 8(2), 91β106. Nayfeh A. H., Sanchez N. E., 1989. Bifurcations in a forced softening Duffing oscillator, Int. J. Non-Linear Mech., 24: 483-497. Reynolds Thomas, Harris Richard, Chang Wen-Shao, 2014, Nonlinear pre-yield modal properties of timber structures with large-diameter steel dowel connections, Engineering Structures 76 :235β244. Roberts, J. B., and Spanos, P. D. 1986. Stochastic averaging: An approximate method of solving random vibration problems. Int. J. Nonlinear Mech., 21(2), 111β134. Roberts J. B., Vasta M., 2000. Markov modeling and stochastic identification for nonlinear ship rolling in random waves, Phil. Trans. R. Soc. Lond. A, vol. 358: 1917-1941. Sabetta, F., and Pugliese, A., 1996, βEstimation of Response Spectra and Simulation of Non-Stationary Earthquake Ground Motions,β Bull. Seismol. Soc. Am., 86, pp. 337β352. Spanos, P. D., and Lutes, L. D. 1980. βProbability of response to evolutionary process.β J. Engrg. Mech. Div., 106(2), 213β224. Spanos P. D., Solomos G. P., 1983. Markov approximation to transient vibration, Journal of Engineering Mechanics, vol. 109: 1134-1150. Spanos P. D., Kougioumtzoglou I. A., 2014a. Galerkin scheme based determination of first-passage probability of nonlinear system response, Structure and Infrastructure Engineering, vol. 10: 1285-1294. Spanos P. D., Kougioumtzoglou I. A., 2014b. Survival probability determination of nonlinear oscillators subject to evolutionary stochastic excitation, ASME Journal of Applied Mechanics, vol. 81, 051016: 1-9. Spyrou K. J., Thomson J. M. T., 2000. The nonlinear dynamics of ship motions: a field overview and some recent developments, Phil. Trans. R. Soc. Lond. A, vol. 358: 1735-1760. Vandewater L. A. and Moss S. D., 2013, Probability-of-existence of vibro-impact regimes in a nonlinear vibration energy harvester, Smart Mater. Struct. 22: 094025. Vanmarcke E. H., 1975. On the distribution of the first-passage time for normal stationary random processes, ASME Journal of Applied Mechanics, vol. 42: 215-220. Zhang Y., Kougioumtzoglou I. A., Nonlinear oscillator stochastic response and survival probability determination via the Wiener path integral, ASCE-ASME Journal of Risk and Uncertainty in Engineering Systems, Part B. Mechanical Engineering (Accepted).
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- International Conference on Applications of Statistics and Probability in Civil Engineering (ICASP) (12th : 2015) /
- An approximate approach for assessing the reliability...
Open Collections
International Conference on Applications of Statistics and Probability in Civil Engineering (ICASP) (12th : 2015)
An approximate approach for assessing the reliability of a stochastically excited softening Duffing oscillator Zhang, Yuanjin; Kougioumtzoglou, Ioannis A. Jul 31, 2015
pdf
Page Metadata
Item Metadata
Title | An approximate approach for assessing the reliability of a stochastically excited softening Duffing oscillator |
Creator |
Zhang, Yuanjin Kougioumtzoglou, Ioannis A. |
Contributor | International Conference on Applications of Statistics and Probability (12th : 2015 : Vancouver, B.C.) |
Date Issued | 2015-07 |
Description | An approximate analytical technique for assessing the reliability of a softening Duffing oscillator subject to evolutionary stochastic excitation is developed. Specifically, relying on a stochastic averaging treatment of the problem the oscillator time-varying survival probability is determined in a computationally efficient manner. In comparison with previous techniques that neglect the potential unbounded response behavior of the oscillator when the restoring force acquires negative values, the herein developed technique readily takes this aspect into account by introducing a special form for the oscillator non-stationary response amplitude probability density function (PDF). A significant advantage of the technique relates to the fact that it can readily handle cases of stochastic excitations that exhibit strong variability in both the intensity and the frequency content. An illustrative numerical example of a softening Duffing oscillator subject to earthquake excitation is included. Comparisons with pertinent Monte Carlo simulation data demonstrate the efficiency of the technique. |
Genre |
Conference Paper |
Type |
Text |
Language | eng |
Notes | This collection contains the proceedings of ICASP12, the 12th International Conference on Applications of Statistics and Probability in Civil Engineering held in Vancouver, Canada on July 12-15, 2015. Abstracts were peer-reviewed and authors of accepted abstracts were invited to submit full papers. Also full papers were peer reviewed. The editor for this collection is Professor Terje Haukaas, Department of Civil Engineering, UBC Vancouver. |
Date Available | 2015-05-13 |
Provider | Vancouver : University of British Columbia Library |
Rights | Attribution-NonCommercial-NoDerivs 2.5 Canada |
DOI | 10.14288/1.0076037 |
URI | http://hdl.handle.net/2429/53148 |
Affiliation |
Non UBC |
Citation | Haukaas, T. (Ed.) (2015). Proceedings of the 12th International Conference on Applications of Statistics and Probability in Civil Engineering (ICASP12), Vancouver, Canada, July 12-15. |
Peer Review Status | Unreviewed |
Scholarly Level | Faculty |
Rights URI | http://creativecommons.org/licenses/by-nc-nd/2.5/ca/ |
Aggregated Source Repository | DSpace |
Download
- Media
- 53032-Paper_128_Zhang.pdf [ 488.24kB ]
- Metadata
- JSON: 53032-1.0076037.json
- JSON-LD: 53032-1.0076037-ld.json
- RDF/XML (Pretty): 53032-1.0076037-rdf.xml
- RDF/JSON: 53032-1.0076037-rdf.json
- Turtle: 53032-1.0076037-turtle.txt
- N-Triples: 53032-1.0076037-rdf-ntriples.txt
- Original Record: 53032-1.0076037-source.json
- Full Text
- 53032-1.0076037-fulltext.txt
- Citation
- 53032-1.0076037.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.53032.1-0076037/manifest