instructions for preparing manuscripts using msword

2 downloads 0 Views 1MB Size Report
Santorini, Greece, June 15-18, 2014. 1. RELIABILITY ASSESSMENT OF NONLINEAR MDOF. SYSTEMS SUBJECT TO EVOLUTIONARY STOCHASTIC.
Computational Stochastic Mechanics โ€“ Proc. of the 7th International Conference (CSM-7) G. Deodatis and P.D. Spanos (eds.) Santorini, Greece, June 15-18, 2014

RELIABILITY ASSESSMENT OF NONLINEAR MDOF SYSTEMS SUBJECT TO EVOLUTIONARY STOCHASTIC EXCITATION IOANNIS P. MITSEAS1*, IOANNIS A. KOUGIOUMTZOGLOU2, POL D. SPANOS3, M. BEER1+ 1

Institute for Risk and Uncertainty, University of Liverpool, Brodie Tower, Brownlow Street, L69 3GQ, Liverpool, United Kingdom. E-mail: *[email protected], [email protected] 2 Department of Civil Engineering and Engineering Mechanics, Columbia University, 500 West 120th Street, New York, NY 10027,USA. E-mail: [email protected] 3 L.B. Ryon Chair in Engineering, Rice University MS 321, P.O. Box 1892, Houston, TX 77251, USA. E-mail: [email protected] An approximate analytical technique for determining the survival probability and first-passage probability density function (PDF) of nonlinear multi-degree-of-freedom (MDOF) structural systems subject to a non-stationary stochastic excitation vector is developed. The proposed technique can be construed as a two-stage approach. First, relying on statistical linearization and utilizing a dimension reduction approach the nonlinear n-degree-of-freedom system is cast/decoupled into (n) effective single-degree-of-freedom (SDOF) linear time-variant (LTV) systems corresponding to each and every degree of freedom of the original MDOF system. Second, an approximate technique based on stochastic averaging is employed in conjunction with the time-varying stiffness and damping elements of the effective SDOF LTV oscillator for determining the survival probability and first-passage PDF of the nonlinear MDOF system. Overall, the technique appears to be efficient and versatile since it can handle readily, at a low computational cost, a wide range of nonlinear/hysteretic behaviors as well as various stochastic excitation forms, even of the fully non-stationary in time and frequency kind. A 3-DOF structural system exhibiting hysteresis following the Bouc-Wen model is included as a numerical example. The reliability of the technique is demonstrated by pertinent Monte Carlo simulations. Keywords: first-passage problem, evolutionary stochastic processes, hysteretic system, statistical linearization, stochastic averaging, reliability assessment

1

Introduction

Structural systems are subject to stochastic excitations that exhibit strong variability in both intensity and frequency content. This fact necessitates the representation of this class of structural loads by non-stationary stochastic

1

processes. Further, structural systems under severe dynamic excitations, such as those due to earthquakes, can exhibit significant nonlinear behavior of the hysteretic kind. Thus, there is considerable interest within the engineering community for determining the response and assessing the reliability of

Reliability Assessment of Nonlinear MDOF Systems Subject to Evolutionary Stochastic Excitation I.P. Mitseas, I.A. Kougioumtzoglou, P. D. Spanos, and M. Beer

nonlinear/hysteretic systems subject to evolutionary stochastic excitations (Roberts and Spanos 2003). In engineering dynamics, the evaluation of the probability that the system response stays within prescribed limits for a specified time interval is advantageous for reliability based system design applications. The determination of the above time-variant probability, known as survival probability, has been coined in the literature as the firstpassage problem and has been a persistent challenge in the field of stochastic dynamics. In this regard, several researchers have focused on developing Monte Carlo simulation (MCS) based techniques (e.g. Schueller et al. 2004) for structural system reliability assessment. Note, however, that there are cases where the computational cost can be prohibitive; thus, rendering the development of alternative efficient approximate analytical/numerical techniques necessary (e.g. Kougioumtzoglou and Spanos 2013). In this paper, an approximate analytical technique for determining the survival probability and first-passage probability density function (PDF) of nonlinear multidegree-of-freedom (MDOF) structural systems subject to a non-stationary stochastic excitation vector is developed. Specifically, first relying on a statistical linearization based dimension reduction approach the original MDOF system is cast into effective singledegree-of-freedom (SDOF) linear time-variant (LTV) systems corresponding to every degree of freedom of the original MDOF system. Second, a stochastic averaging based approximate technique is utilized to derive the nonlinear MDOF system survival probability and first-passage PDF at a low computational cost. A Bouc-Wen 3-DOF structural system subject to an evolutionary stochastic earthquake excitation is included as a numerical example. Pertinent MCS data demonstrate the reliability of the technique.

2

2 2.1

MDOF System Dimension Reduction Statistical linearization treatment

Consider an n-degree-of-freedom nonlinear structural system governed by the equation ๐‘ด๐’’ + ๐‘ช๐’’ + ๐‘ฒ๐’’ + ๐’ˆ ๐’’, ๐’’ = ๐‘ญ ๐‘ก

(1)

where q denotes the response acceleration vector, q is the response velocity vector, q is the response displacement vector; M, C and K denote the n ร— n mass, damping and stiffness matrices, respectively. g q, q is an arbitrary nonlinear n ร— 1 vector function of the variables q and q. F t T = (f1 t , f2 t , โ€ฆ , fn t ) is a n ร— 1 zero mean, non-stationary stochastic excitation vector process defined as F t = ฮณa t , where ฮณT = (ฮณ1 , ฮณ2 , โ€ฆ , ฮณn ) is an arbitrary constant vector of weighting coefficients, and a t is a non-stationary process with an evolutionary power spectrum (EPS) Sa ฯ‰, t . In this regard, F t possesses the EPS matrix ๐‘บ๐‘ญ ๐Ž, ๐’• = ๐›พ1 2 ๐‘†๐‘Ž ๐œ”, ๐‘ก 0 โ‹ฏ 0 2 0 ๐›พ2 ๐‘†๐‘Ž ๐œ”, ๐‘ก โ‹ฏ 0 โ‹ฎ โ‹ฑ โ‹ฏ โ‹ฎ 0 0 โ‹ฏ ๐›พ๐‘› 2 ๐‘†๐‘Ž ๐œ”, ๐‘ก

(2)

Further, the non-stationary stochastic excitation process is regarded to be a filtered stationary stochastic process according to the concept proposed by Priestley (Priestley 1965); see also (Dalhaus 1997, Spanos and Kougioumtzoglou 2012, Kougioumtzoglou 2013). Thus, the excitation EPS matrix Eq.(2) can be given by ๐‘บ๐‘ญ ๐œ”, ๐‘ก = ๐‘จ ๐œ”, ๐‘ก ๐‘บ๐‘ญ ๐œ” ๐œœ(๐œ”, ๐‘ก)๐‘ปโˆ—

(3)

where the superscripts (T) and (*) denote matrix transposition and complex conjugation, respectively; A ฯ‰, t is the modulating matrix which serves as a time-variant filter; and

Computational Stochastic Mechanics โ€“ Proc. of the 7th International Conference (CSM-7) G. Deodatis and P.D. Spanos (eds.) Santorini, Greece, June 15-18, 2014

SF ฯ‰ is the power spectrum matrix corresponding to the stationary stochastic vector process F(t). Note that both separable and non-separable EPS can be defined considering Eq.(3). In this manner, realistic excitations exhibiting variability in both intensity and frequency content can be considered. Focusing next on the frequency domain, the response determination problem is defined as seeking the corresponding system response EPS matrix Sq ฯ‰, t . In the following a statistical linearization approach (Roberts and Spanos 2003) is employed for determining the response EPS matrix Sq ฯ‰, t . In this regard, a linearized version of Eq.(1) is given in the form

๐‘ก

๐’‰ ๐‘ก โˆ’ ๐œ ๐‘จ ๐œ”, ๐œ ๐‘’ โˆ’๐‘–๐œ” (๐‘กโˆ’๐œ) ๐‘‘๐œ

(8)

0

In Eq.(8) h t denotes the impulse response function matrix. Furthermore, the time dependent cross โ€“ variance of the response can be evaluated by the expression โˆž

๐ธ ๐‘ž๐‘– ๐‘ž๐‘— = โˆ’โˆž

๐‘†๐‘ž ๐‘– ๐‘ž ๐‘— ๐œ”, ๐‘ก ๐‘‘๐œ”

(9)

Next, omitting the convolution of the impulse response function matrix with the modulating matrix can lead to substantial reduction of the computational effort, especially for the case of MDOF systems (Kougioumtzoglou and Spanos 2013, Mitseas et al. 2014). In this manner, Eq.(8) takes the form

๐‘ด๐’’ + ๐‘ช + ๐‘ช๐’†๐’’ ๐’’ + ๐‘ฒ + ๐‘ฒ๐’†๐’’ ๐’’ = ๐‘ญ ๐‘ก (4)

๐‘ฏ๐’ˆ๐’†๐’ ๐œ”, ๐‘ก = ๐‘ฏ ๐œ” ๐‘จ ๐œ”, ๐‘ก

Following (Roberts and Spanos 2003) and adopting the standard assumption that the response processes are Gaussian, the timedependent elements of the equivalent linear matrices Ceq and K eq are given by the expressions

where H ฯ‰ is the frequency response function (FRF) matrix defined as

๐‘’๐‘ž

๐‘๐‘–,๐‘— = ๐ธ

๐œ•๐‘”๐‘– ๐œ•๐‘ž๐‘—

(5)

๐‘’๐‘ž ๐‘˜๐‘–,๐‘—

๐‘ฏ ๐œ” = (โˆ’๐œ”2 ๐‘ด + ๐‘–๐œ” ๐‘ช + ๐‘ช๐’†๐’’ + โ‹ฏ โ€ฆ + (๐‘ฒ + ๐‘ฒ๐’†๐’’ ))โˆ’๐Ÿ

๐œ•๐‘”๐‘– =๐ธ ๐œ•๐‘ž๐‘—

(6)

Further, for a linear MDOF system subject to evolutionary stochastic excitation a matrix input-output spectral relationship of the form โˆ—

๐‘บ๐’’ ๐œ”, ๐‘ก = ๐‘ฏ๐’ˆ๐’†๐’ ๐œ”, ๐‘ก ๐‘บ๐‘ญ ๐œ” ๐‘ฏ๐‘ป

๐’ˆ๐’†๐’

๐œ”, ๐‘ก (7)

can be derived (Kougioumtzoglou and Spanos 2009), where ๐‘ฏ๐’ˆ๐’†๐’ ๐œ”, ๐‘ก =

3

(11)

Consequently, taking into account Eq.(3), Eq. (7) becomes ๐‘บ๐’’ ๐œ”, ๐‘ก = ๐‘ฏ ๐œ” ๐‘บ๐‘ญ ๐œ”, ๐‘ก ๐‘ฏ๐‘ปโˆ— ๐œ”

and

(10)

(12)

Note that Eq.(12) constitutes a straightforward generalization of the celebrated spectral relationship based on stationarity and on the Wiener-Khinchin theorem. In this regard, the abovementioned approximation can be construed as a quasi-stationary approach which, in general, yields satisfactory accuracy in cases of relatively stiff systems (Hammond 1973, Jangid and Datta 1999). Note in passing that the spectral input-output relationship of Eq.(12) is exact for the case of stationary processes (Roberts and Spanos 2003). Further, adopting the aforementioned quasi-stationary approach, it can be readily seen that for the i-

Reliability Assessment of Nonlinear MDOF Systems Subject to Evolutionary Stochastic Excitation I.P. Mitseas, I.A. Kougioumtzoglou, P. D. Spanos, and M. Beer

th degree of freedom, using Eq.(2), Eq.(9) and Eq.(12) yields โˆž 2

๐ธ ๐‘ž๐‘– ๐‘ก

=

( ๐ป๐‘–1 ๐œ”

2

๐›พ1 + โ‹ฏ

๐œ”2

โˆ’โˆž

โˆ’โˆž

(13)

and โˆž

=

+ ๐ป๐‘–๐‘› (๐œ”)

โˆ’โˆž 2 2 ๐›พ๐‘› )๐‘†๐‘Ž

๐œ”2 ( ๐ป๐‘–1 ๐œ”

2

๐›พ1 2 + โ‹ฏ

๐œ”, ๐‘ก ๐‘‘๐œ”

(14)

Eqs.(13) and (14) hold true in the approximate quasi-stationary sense delineated earlier. Clearly, Eq.(12) constitutes an approximate formula for determining the MDOF system response EPS matrix at a low computational cost; thus, circumventing computationally intensive Monte Carlo simulations. 2.2 Effective SDOF linear time-variant system Adopting next the system dimension reduction approach developed in (Kougioumtzoglou and Spanos 2013), an auxiliary effective SDOF LTV system corresponding to the i-th degree of freedom can be defined as 2 ๐‘ž๐‘– + ๐›ฝ๐‘’๐‘ž ,๐‘– ๐‘ก ๐‘ž๐‘– + ๐œ”๐‘’๐‘ž ,๐‘– ๐‘ก ๐‘ž๐‘– = ๐‘“๐‘– ๐‘ก

(15)

where the time-varying stiffness and damping elements of the effective LTV system can be determined by equating the variances of the response displacement and velocity expressed utilizing Eq.(15) with the corresponding ones (Eqs.(13-14)) evaluated via the statistical linearization approach delineated in the previous section; that is, E ๐‘ž๐‘– 2 ๐‘ก โˆž

โˆ’โˆž

=

1 2 2 )2 + (๐›ฝ 2 (๐œ”๐‘’๐‘ž ๐‘ก โˆ’ ๐œ” ๐‘’๐‘ž ,๐‘– ๐‘ก ๐œ”) ,๐‘–

ร— ๐›พ๐‘– 2 ๐‘†๐‘Ž ๐œ”, ๐‘ก ๐‘‘๐œ”

E ๐‘ž๐‘– 2 ๐‘ก โˆž

2

+ ๐ป๐‘–๐‘› (๐œ”) 2 ๐›พ๐‘› 2 )๐‘†๐‘Ž ๐œ”, ๐‘ก ๐‘‘๐œ”

๐ธ ๐‘ž๐‘– 2 ๐‘ก

and

(16)

4

= 1 2 2 )2 + (๐›ฝ 2 (๐œ”๐‘’๐‘ž ๐‘ก โˆ’ ๐œ” ๐‘’๐‘ž ,๐‘– ๐‘ก ๐œ”) ,๐‘–

ร— ๐›พ๐‘– 2 ๐‘†๐‘Ž ๐œ”, ๐‘ก ๐‘‘๐œ”

(17)

Clearly, Eqs.(16) and (17) in conjunction with Eqs.(13) and (14) constitute a nonlinear system of two algebraic equations for the evaluation of the LTV system time-varying stiffness and damping elements. Further, considering the case where the LTV system of Eq.(15) is a lightly damped system, and relying on a combination of deterministic and stochastic averaging (e.g. Spanos 1978), it was shown in (Spanos and Lutes 1980, Kougioumtzoglou and Spanos 2009) that the non-stationary LTV system response amplitude PDF is given by the time-dependent Rayleigh distribution of the form ๐‘ ๐‘Ž๐‘– , ๐‘ก =

๐‘Ž๐‘– ๐‘Ž๐‘– 2 exp โˆ’ ๐‘๐‘– (๐‘ก) 2๐‘๐‘– (๐‘ก)

(18)

In Eq.(18) ๐‘Ž๐‘– represents the LTV oscillator response amplitude defined as ๐‘Ž๐‘– 2 ๐‘ก = ๐‘ž๐‘–2 ๐‘ก +

๐‘ž๐‘– ๐‘ก ๐œ”๐‘’๐‘ž ,๐‘– (๐‘ก)

2

(19)

where ci (t) accounts for the non-stationary variance of the LTV system of Eq.(15). As it was shown in (Spanos and Lutes 1980, Kougioumtzoglou and Spanos 2009) the nonstationary LTV system response is given by ๐‘๐‘– ๐‘ก = โˆ’๐›ฝ๐‘’๐‘ž ,๐‘– ๐‘ก ๐‘๐‘– ๐‘ก +

๐œ‹๐‘†๐‘“ ๐œ”๐‘’๐‘ž ,๐‘– ๐‘ก , ๐‘ก (20) 2 ๐œ”๐‘’๐‘ž ,๐‘– ๐‘ก

Eq.(20) constitutes a first-order nonlinear ordinary differential equation (ODE) which can be solved utilizing standard numerical schemes such as the Runge-Kutta; see references (Spanos and Lutes 1980, Roberts and Spanos 2003, Kougioumtzoglou and

Computational Stochastic Mechanics โ€“ Proc. of the 7th International Conference (CSM-7) G. Deodatis and P.D. Spanos (eds.) Santorini, Greece, June 15-18, 2014

Spanos 2013) for a more detailed presentation of the topic.

๐‘š

๐‘ƒ๐‘–๐ต ๐‘‡ =

๐ต 1 โˆ’ ๐น๐‘–,๐‘—

(24)

๐‘— =1

3

MDOF Nonlinear System Reliability Assessment

The survival probability ๐‘ƒ๐‘–๐ต is defined as the probability that the system response amplitude ๐‘Ž๐‘– stays below a prescribed barrier ๐ต over the time interval [0, ๐‘‡], given that ๐‘Ž๐‘– ๐‘ก = 0 < ๐ต. In the ensuing analysis, a recently developed approximate technique (Spanos and Kougioumtzoglou 2014) is employed for determining the survival probability and the first-passage PDF of the LTV system of Eq.(15). The first-passage PDF ๐‘๐‘–๐ต ๐‘‡ and the survival probability ๐‘ƒ๐‘–๐ต are related according to the expression ๐‘๐‘–๐ต ๐‘‡ = โˆ’

๐‘‘๐‘ƒ๐‘–๐ต ๐‘‡ ๐‘‘๐‘‡

(21)

Next, adopting the discretization scheme employed in (Spanos and Solomos 1983) yields intervals of the form

๐‘‡๐‘’๐‘ž ,๐‘– (๐‘ก ๐‘–,๐‘— โˆ’1 ) 2

(22)

where Teq ,i represents the LTV system equivalent natural period given by ๐‘‡๐‘’๐‘ž ,๐‘– ๐‘ก =

2๐œ‹ ๐œ”๐‘’๐‘ž ,๐‘– ๐‘ก

where is defined as the probability that the response amplitude ๐‘Ž๐‘– will exceed the prescribed barrier B over the time interval t i,jโˆ’1 , t i,j , given that no crossings have occurred prior to time t i,jโˆ’1 . Next, invoking the Markovian property of the response amplitude ๐‘Ž๐‘– , one gets ๐ต ๐น๐‘–,๐‘— = ๐‘ƒ๐‘Ÿ๐‘œ๐‘ ๐‘Ž๐‘– ๐‘ก๐‘–,๐‘— โ‰ฅ ๐ต โˆฉ ๐‘Ž๐‘– ๐‘ก๐‘–,๐‘— โˆ’1 < ๐ต

๐‘ƒ๐‘Ÿ๐‘œ๐‘ ๐‘Ž๐‘– ๐‘ก๐‘–,๐‘— โˆ’1 < ๐ต

(23)

=

๐ต ๐ป๐‘–,๐‘— โˆ’1,๐‘—

=

(25)

๐ต ๐ป๐‘–,๐‘— โˆ’1

where โˆฉ denotes the intersection of sets. B Utilizing Eq.(18) Hi,jโˆ’1 can be determined as ๐ต ๐ต ๐ป๐‘–,๐‘— โˆ’1 =

๐‘ ๐‘Ž๐‘–,๐‘— โˆ’1 , ๐‘ก๐‘–,๐‘— โˆ’1 ๐‘‘๐‘Ž๐‘–,๐‘— โˆ’1 = 0

= 1 โˆ’ ๐‘’๐‘ฅ๐‘ โˆ’

๐‘ก๐‘–,๐‘— โˆ’1 , ๐‘ก๐‘–,๐‘— , ๐‘— = 1,2, โ€ฆ , ๐‘š, ๐‘ก๐‘–,0 = 0, ๐‘ก๐‘–,๐‘š = ๐‘‡, ๐‘ก๐‘–,๐‘— โˆ’ ๐‘ก๐‘–,๐‘— โˆ’1 =

B Fi,j

๐ต2 2๐‘๐‘– (๐‘ก๐‘–,๐‘— โˆ’1 )

(26)

B whereas Hi,jโˆ’1,j is defined as a double integral of the form โˆž ๐ต ๐ป๐‘–,๐‘— โˆ’1,๐‘— = ๐ต

๐‘‘๐‘Ž๐‘–,๐‘— ร— ๐ต

๐‘ ๐‘Ž๐‘–,๐‘— โˆ’1 , ๐‘ก๐‘–,๐‘— โˆ’1 ; ๐‘Ž๐‘–,๐‘— , ๐‘ก๐‘–,๐‘— ๐‘‘๐‘Ž๐‘–,๐‘— โˆ’1

(27)

0

Considering the slowly time-varying behavior of the response amplitude ๐‘Ž๐‘– , it is assumed that it is constant over the time interval t i,jโˆ’1 , t i,j . In this regard, the survival probability PiB also has a constant value over the same time interval. Obviously, the survival probability is given by

5

Assuming next that ฮฒeq ,i t i,j = ฮฒeq ,i t i,jโˆ’1 and ฯ‰eq ,i t i,j = ฯ‰eq ,i t i,jโˆ’1 over the time interval t i,jโˆ’1 , t i,j , it has been shown in (Spanos and Solomos 1983, Spanos and Kougioumtzoglou 2014) that the joint response amplitude PDF p ฮฑi,jโˆ’1 , t i,jโˆ’1 ; ฮฑi,j , t i,j is given in the form

Reliability Assessment of Nonlinear MDOF Systems Subject to Evolutionary Stochastic Excitation I.P. Mitseas, I.A. Kougioumtzoglou, P. D. Spanos, and M. Beer

๐ฟ๐‘–,๐‘› ๐‘ 2 ๐‘› =1 2๐‘›

๐‘ ๐›ผ๐‘–,๐‘— โˆ’1 , ๐‘ก๐‘–,๐‘— โˆ’1 ; ๐›ผ๐‘–,๐‘— , ๐‘ก๐‘–,๐‘— = ๐›ผ๐‘–,๐‘— โˆ’1 ๐›ผ๐‘–,๐‘— ร— 2 ๐‘๐‘– ๐‘ก๐‘–,๐‘— โˆ’1 ๐‘๐‘– ๐‘ก๐‘–,๐‘— 1 โˆ’ ๐‘Ÿ๐‘–,๐‘— 2 2 ๐›ผ๐‘–,๐‘— ๐‘ ๐‘ก๐‘–,๐‘— โˆ’1 + ๐›ผ๐‘–,๐‘— โˆ’1 ๐‘๐‘– ๐‘ก๐‘–,๐‘— ๐‘’๐‘ฅ๐‘ โˆ’ 2 2๐‘๐‘– ๐‘ก๐‘–,๐‘— โˆ’1 ๐‘๐‘– ๐‘ก๐‘–,๐‘— 1 โˆ’ ๐‘Ÿ๐‘–,๐‘—

with ร—

๐ฟ๐‘–,๐‘› =

2 1 โˆ’ ๐‘Ÿ๐‘–,๐‘—

๐‘๐‘– ๐‘ก๐‘–,๐‘— โˆ’1 ๐‘๐‘– ๐‘ก๐‘–,๐‘—

๐‘›+1

4๐‘› ๐‘๐‘– ๐‘ก๐‘–,๐‘— โˆ’1

๐›ผ๐‘–,๐‘— ๐›ผ๐‘–,๐‘— โˆ’1 ๐‘Ÿ๐‘–,๐‘—

๐ผ0

(33)

(28)

๐›ค๐‘– 1 + ๐‘›,

๐‘›+2

ร—

2 2๐‘๐‘– ๐‘ก๐‘–,๐‘— โˆ’1 1 โˆ’ ๐‘Ÿ๐‘–,๐‘—

โˆ’ ๐›ค๐‘– 1 + ๐‘›, 0

where 2 ๐‘Ÿ๐‘–,๐‘— =

2 ๐‘๐‘– ๐‘ก๐‘–,๐‘— 1 โˆ’ ๐‘Ÿ๐‘–,๐‘— ๐ต2

๐‘๐‘– ๐‘ก๐‘–,๐‘— โˆ’1 ๐‘๐‘– ๐‘ก๐‘–,๐‘—

1 โˆ’ ๐›ฝ๐‘’๐‘ž ,๐‘– ๐‘ก๐‘–,๐‘— โˆ’1 (๐‘ก๐‘–,๐‘— โˆ’1 โˆ’ ๐‘ก๐‘–,๐‘— )

.

(29)

Next, I0 (. ) denotes the modified Bessel function of the first kind of order zero. I0 (. ) can be further expanded in the form ๐‘ฅ2 ๐‘ฅ4 ๐‘ฅ6 ๐ผ0 ๐‘ฅ = 1 + 2 + 2 2 + 2 2 2 + โ‹ฏ 2 2 4 2 4 6

(30)

Considering Eqs.(27), (28) and (30) and manipulating yields N ๐ต ๐ป๐‘–,๐‘— โˆ’1,๐‘—

= Ai,0 +

Ai,n

(31)

n=1

+

๐ต2 ๐‘๐‘– ๐‘ก๐‘–,๐‘—

โˆ’ ๐›ค๐‘– 1 + ๐‘›,

๐ด๐‘–,0 = ๐‘’๐‘ฅ๐‘ 1 โˆ’ ๐‘’๐‘ฅ๐‘

ii. ร— ร—

2 2๐‘๐‘– ๐‘ก๐‘–,๐‘— โˆ’1 1 โˆ’ ๐‘Ÿ๐‘–,๐‘—

iii.

2 1 โˆ’ ๐‘Ÿ๐‘–,๐‘—

(32) iv.

and ๐ด๐‘–,๐‘›

2๐‘› ๐‘Ÿ๐‘–,๐‘—

= ๐‘๐‘– ๐‘ก๐‘–,๐‘— โˆ’1 ๐‘๐‘– ๐‘ก๐‘–,๐‘—

๐‘›+1

2๐‘›+1

1โˆ’

v.

ร—

2 ๐‘Ÿ๐‘–,๐‘—

6

2๐‘๐‘– ๐‘ก๐‘–,๐‘—

โˆ’๐‘›

๐›ค๐‘– 1 + ๐‘› ๐ต2 1 โˆ’ ๐‘Ÿ๐‘—2

๐ต2๐‘› (34)

Concisely, the developed technique comprises of the following steps: i.

2 2๐‘๐‘– ๐‘ก๐‘–,๐‘— 1 โˆ’ ๐‘Ÿ๐‘–,๐‘— 2 โˆ’๐ต

2 1 โˆ’ ๐‘Ÿ๐‘–,๐‘—

๐›ค๐‘– 1 + ๐‘›

In Eq.(34) ฮ“i z represents the Gamma โˆž function defined as ฮ“i z = 0 t zโˆ’1 eโˆ’t dt; and ฮ“i [ฮณ, z] represents the incomplete Gamma โˆž function defined as ฮ“i [ฮณ, z] = z t ฮณโˆ’1 eโˆ’t dt.

where โˆ’๐ต2

๐‘›

2 1 โˆ’ ๐‘Ÿ๐‘–,๐‘—

โˆ’ ๐‘๐‘– ๐‘ก๐‘–,๐‘—

ร—

Determination of the MDOF system response variances (Eqs.(13-14)). Determination of the time-varying elements ๐›ฝ๐‘’๐‘ž ,๐‘– ๐‘ก and ๐œ”๐‘’๐‘ž ,๐‘– ๐‘ก by solving the 2x2 system of algebraic equations (Eqs.(16-17)). Determination of ๐‘๐‘– ๐‘ก via numerically integrating the ODE Eq.(20). Determination of the equivalent natural period ๐‘‡๐‘’๐‘ž ,๐‘– ๐‘ก (Eq.(23)) and discretization of the time domain via Eq.(22). Determination of the parameters ๐ต ๐ต ๐ป๐‘–,๐‘— โˆ’1 and ๐ป๐‘–,๐‘— โˆ’1,๐‘— via Eqs.(26) and (27).

Computational Stochastic Mechanics โ€“ Proc. of the 7th International Conference (CSM-7) G. Deodatis and P.D. Spanos (eds.) Santorini, Greece, June 15-18, 2014

vi.

4

Determination of the survival probability ๐‘ƒ๐‘–๐ต ๐‘‡ via Eq.(24) and of the corresponding first-passage PDF ๐‘๐‘–๐ต ๐‘‡ via Eq.(21). Numerical Application

In this section, a nonlinear three-degree-offreedom structural system following the BoucWen hysteretic model (Wen 1980, Ikhouane and Rodellar 2007) subject to evolutionary stochastic earthquake excitation is considered to demonstrate the reliability of the technique. The survival probabilities and the firstpassage PDFs obtained via the developed approximate technique are compared with survival probability and first-passage PDF estimates obtained via pertinent Monte Carlo simulations (10,000 realizations). A standard fourth-order Runge-Kutta numerical integration scheme is employed for solving the nonlinear system differential equation of motion (Eq.(1)), whereas the barrier level B is expressed as a fraction ฮป of the maximum value of the non-stationary response displacement standard deviation of the first floor, i.e. B = ฮปmaxโก ( c1 (t)). Further, considering relative displacements, the 3-DOF nonlinear structural system is governed by the equation ๐‘ด๐’’ + ๐‘ช๐’’ + ๐‘ฒ๐’’ + ๐’ˆ ๐’’, ๐’’ = ๐‘ญ ๐‘ก

(35)

where ๐ชT = q1 q2 q3 z1 z2 z3

(36)

and ๐Œ=

๐Œ๐Ÿ๐Ÿ ๐Œ๐Ÿ๐Ÿ

๐Œ๐Ÿ๐Ÿ ๐Œ๐Ÿ๐Ÿ

(37)

with

๐Œ๐Ÿ๐Ÿ

m1 m = 2 m3

0 m2 m3

0 0 m3

and 0 ๐Œ๐Ÿ๐Ÿ = ๐Œ๐Ÿ๐Ÿ = ๐Œ๐Ÿ๐Ÿ = 0 0

0 0 0

0 0 0

(39)

In Eq.(38) m1 = 2.0615 ร— 105 , m2 = 5 2.0559 ร— 10 and m3 = 2.0261 ร— 105 . Further, ๐Š=

๐Š ๐Ÿ๐Ÿ ๐Š ๐Ÿ๐Ÿ

๐Š ๐Ÿ๐Ÿ ๐Š ๐Ÿ๐Ÿ

(40)

with ๐Š ๐Ÿ๐Ÿ =

ฮฑk1 0 0

โˆ’ฮฑk 2 ฮฑk 2 0

๐Š ๐Ÿ๐Ÿ (1 โˆ’ ฮฑ)k1 0 = 0

0 โˆ’ฮฑk 3 , ฮฑk 3

โˆ’(1 โˆ’ ฮฑ)k 2 (1 โˆ’ ฮฑ)k 2 0

(41)

0 โˆ’(1 โˆ’ ฮฑ)k 3 (42) (1 โˆ’ ฮฑ)k 3

and 0 ๐Š ๐Ÿ๐Ÿ = ๐Š ๐Ÿ๐Ÿ = 0 0

0 0 0

0 0 0

(43)

In Eqs.(41) and (42) k1 = 3.9668 ร— 108 , k 2 = 3.5007 ร— 108 and k 3 = 2.6927 ร— 108 whereas ฮฑ stands for the rigidity ratio which can be viewed as a form of post-yield to preyield stiffness ratio. Further, the 6-by-6 damping matrix of the structural system C is assumed to be partly proportional to the stiffness matrix; that is, ๐‚=

๐‚๐Ÿ๐Ÿ ๐‚๐Ÿ๐Ÿ

๐‚๐Ÿ๐Ÿ ๐‚๐Ÿ๐Ÿ

with ๐‚๐Ÿ๐Ÿ = c๐Š ๐Ÿ๐Ÿ and

7

(38)

(44)

Reliability Assessment of Nonlinear MDOF Systems Subject to Evolutionary Stochastic Excitation I.P. Mitseas, I.A. Kougioumtzoglou, P. D. Spanos, and M. Beer

0 ๐‚๐Ÿ๐Ÿ = ๐‚๐Ÿ๐Ÿ = 0 0 1 C22 = 0 0

0 1 0

0 0 0

0 0, 0

and (45) ๐Š ๐ž๐ช =

0 0 1

(46)

where c is taken equal to 2 ร— 10โˆ’3 . For the specific example ฮณ = ๐‘š๐‘– , and the loading vector becomes ๐…(t)T = f t f t f t 0 0 0

(47)

Further, ๐  ๐ช, ๐ช

๐“

= ( 0 0 0 โˆ’g1 q1 , z1 โ€ฆ โˆ’g 2 q2 , z2 โˆ’ g 3 q3 , z3 )

(48)

In the Bouc-Wen model the additional state zi is associated with the displacement qi via the equation z i = g i q i , zi

(49)

๐Š ๐ž๐ช๐Ÿ๐Ÿ ๐Š ๐ž๐ช๐Ÿ๐Ÿ

๐Š ๐ž๐ช๐Ÿ๐Ÿ , ๐Š ๐ž๐ช๐Ÿ๐Ÿ

(54)

0 ๐Š ๐ž๐ช๐Ÿ๐Ÿ = ๐Š ๐ž๐ช๐Ÿ๐Ÿ = ๐Š ๐ž๐ช๐Ÿ๐Ÿ = 0 0 k eq 1 0 0

๐Š ๐ž๐ช๐Ÿ๐Ÿ =

0 k eq 2 0

0 0 0

0 0 0

0 0

(55)

(56)

k eq 3

Further, the elements ceq i and k eq i are given by the expressions

ceq i =

2 E(qi zi ) ฮณ + ฮฒ E(zi 2 ) โˆ’ A ฯ€ 2 E(qi )

(57)

2 E(qi zi ) ฮณ E(qi 2 ) + ฮฒ ฯ€ E(zi 2 )

(58)

and k eq i =

where g i q i , zi = โˆ’ฮณ q i zi zi โ€ฆ + ฮ‘q i

nโˆ’1

โˆ’ ฮฒqi zi

n

โ€ฆ (50)

The parameters ฮณ, ฮฒ, ฮ‘ and n are capable of representing a wide range of hysteresis loops (Ikhouane and Rodellar 2007). In the study herein the following values are considered: a = 0.15, ฮฒ=ฮณ=0,5, n=1 and A=1. The equivalent linear matrices take the form (e.g. Roberts and Spanos 2003) ๐‚๐ž๐ช =

๐‚๐ž๐ช๐Ÿ๐Ÿ ๐‚๐ž๐ช๐Ÿ๐Ÿ

๐‚๐ž๐ช๐Ÿ๐Ÿ , ๐‚๐ž๐ช๐Ÿ๐Ÿ

0 ๐‚๐ž๐ช๐Ÿ๐Ÿ = ๐‚๐ž๐ช๐Ÿ๐Ÿ = ๐‚๐ž๐ช๐Ÿ๐Ÿ = 0 0 ๐‚๐ž๐ช๐Ÿ๐Ÿ

ceq 1 = 0 0

0 ceq 2 0

0 0 0 0 0

0 0 0

respectively. Furthermore, the excitation EPS Sฮฑ ฯ‰, t takes the form Sฮฑ ฯ‰, t = g t

2

SCP ฯ‰

(59)

where SCP ฯ‰ represents the widely used in earthquake engineering Clough-Penzien power spectrum (Clough and Penzien 1993) and g t denotes a time-modulating envelope function given by

(51)

g t = k eโˆ’at โˆ’ eโˆ’bt

(52)

where a=0.1 and b=0.3; and k is a normalization constant so that g t max = 1. The Clough-Penzien spectrum is given by

(53)

ceq 3

8

SCP ฯ‰ =

(60)

Computational Stochastic Mechanics โ€“ Proc. of the 7th International Conference (CSM-7) G. Deodatis and P.D. Spanos (eds.) Santorini, Greece, June 15-18, 2014

S0

(ฯ‰/ฯ‰f )4 1 โˆ’ (ฯ‰/ฯ‰f )2 2 + 4ฮพf 2 (ฯ‰/ฯ‰f )2 ฯ‰g 4 + 4(ฮพg )2 ฯ‰g 2 ฯ‰2

(ฯ‰g 2 โˆ’ ฯ‰2 )2 + 4ฮพg 2 ฯ‰g 2 ฯ‰2

ร— (61)

where S0 is the amplitude of the bedrock excitation spectrum, modeled as a white noise process. The parameters values used are S0 = 20m2 /sec 3 , ฮพg = 0.7, ฯ‰g = 2rad/sec, ฮพf = 0.6 and ฯ‰f = 12.5rad/sec (see Fig.1). The total duration of the excitation is 20 sec.

Figure 1.Nonstationary separable excitation power spectrum Sฮฑ ฯ‰, t

Figure 3: Time-varying equivalent linear damping ฮฒeq ,i t coefficients.

Figure 4: Time-varying survival probability of the systemโ€™s first degree-of-freedom for various values of the parameter ฮป; comparisons with MCS (10,000 realizations).

Figure 2: Time-varying equivalent linear stiffness ฯ‰eq .i t coefficients. Figure 5: First-passage PDF of the systemโ€™s first degree-of-freedom for various values of the parameter ฮป; comparisons with MCS (10,000 realizations).

9

Reliability Assessment of Nonlinear MDOF Systems Subject to Evolutionary Stochastic Excitation I.P. Mitseas, I.A. Kougioumtzoglou, P. D. Spanos, and M. Beer

Figure 6: Time-varying survival probability of the systemโ€™s second degree-of-freedom for various values of the parameter ฮป; comparisons with MCS (10,000 realizations).

Figure 7: First-passage PDF of the systemโ€™s second degree-of-freedom for various values of the parameter ฮป; comparisons with MCS (10,000 realizations).

Figure 8: Time-varying survival probability of the systemโ€™s third degree-of-freedom for various values of the parameter ฮป; comparisons with MCS (10,000 realizations).

10

Figure 9: First-passage PDF of the systemโ€™s third degree-of-freedom for various values of the parameter ฮป; comparisons with MCS (10,000 realizations).

In Figs. (2) and (3) the equivalent timevarying stiffness ฯ‰eq .i t and damping ๐›ฝ๐‘’๐‘ž ,i ๐‘ก elements corresponding to each DOF are plotted. Note that the hysteretic/degrading behavior of the structural system is captured by the decreasing with time trend of the stiffness element, as well as the increasing with time trend of the damping element. Further, in Figs. (4) and (5) the survival probabilities PiB T as well as the corresponding first-passage PDFs pBi T for the first DOF of the hysteretic MDOF structural system are plotted for various barrier levels; comparisons with MCS utilizing a spectral representation approach (Shinozuka and Deodatis 1991) (10,000 realizations) demonstrate a satisfactory degree of agreement. Furthermore, in Figs. (6), (7), (8) and (9) the survival probabilities PiB T as well as the associated first-passage PDFs pBi T corresponding to the second and third DOF are plotted for various barrier levels. Comparisons with MCS demonstrate a satisfactory degree of accuracy in these cases as well. 5

Conclusions

An approximate analytical technique for determining the survival probability and firstpassage PDF of nonlinear/hysteretic MDOF structural systems subject to evolutionary

Computational Stochastic Mechanics โ€“ Proc. of the 7th International Conference (CSM-7) G. Deodatis and P.D. Spanos (eds.) Santorini, Greece, June 15-18, 2014

stochastic excitation has been developed. Specifically, based on an efficient dimension reduction approach and relying on the concepts of stochastic averaging/linearization reliability statistics of the nonlinear MDOF system have been determined at a low computational cost. A significant advantage of the technique relates to the fact that it can readily treat a broad range of nonlinear/hysteretic behaviors as well as stochastic excitations with arbitrary EPS forms. Acknowledgments The first author acknowledges the financial support of the State Scholarships Foundation in Greece (IKY).

References Journal Papers Dahlhaus R., 1997. Fitting time series models to non-stationary processes, Ann Stat; 25: 1โ€“37. Hammond J.K., 1973. Evolutionary spectra in random vibration, J Roy Stat Soc; 35: 167โ€“188. Jangid R.S., Datta T.K, 1999. Evaluation of the methods for response analysis under nonstationary excitation, Shock and Vibration; 6: 285โ€“297. Kougioumtzoglou I. A., 2013. Stochastic joint time-frequency response analysis of nonlinear structural systems, Journal of Sound and Vibration, vol. 332: 7153-7173. Kougioumtzoglou I. A., Spanos P. D., 2009. An approximate approach for nonlinear system response determination under evolutionary stochastic excitation, Current Science, Indian Academy of Sciences, vol. 97: 1203-121. Kougioumtzoglou I. A., Spanos P. D., 2012. An analytical Wiener path integral technique for non-stationary response determination of nonlinear oscillators, Probabilistic Engineering Mechanics, vol. 28: 125-131. Kougioumtzoglou I. A., Spanos P. D., 2013. Nonlinear MDOF system stochastic response determination via a dimension reduction approach, Computers and Structures; 126: 135148.

11

Priestley M.B, 1965. Evolutionary spectra and nonstationary processes. J. Roy Stat Soc; 27: 204โ€“ 237. Schueller G. I., Pradlwarter H. J., Koutsourelakis P. S., 2004. A critical appraisal of reliability estimation procedures for high dimensions, Probabilistic Engineering Mechanics, vol. 19: 463-474. Shinozuka M., Deodatis G., 1991. Simulation of stochastic processes by spectral representation, Applied Mechanics Reviews, vol. 44, no. 4, pp. 191-204. Spanos P.D, 1978. Non-stationary random vibration of a linear structure, International Journal of Solids and Structures, vol 14: 861867. Spanos P. D, Kougioumtzoglou I. A., 2014. Survival probability determination of nonlinear oscillators subject to evolutionary stochastic excitation, Journal of Applied Mechanics, 81(5), 051016. Spanos P. D, Kougioumtzoglou I. A., 2012. Harmonic wavelets based statistical linearization for response evolutionary power spectrum determination, Probabilistic Engineering Mechanics; 27: 57-68. Spanos P. D., Lutes L.D, 1980. Probability of response to evolutionary process, J. Eng Mech. Div. Am. Soc. Civil Eng., 106, 213-224. Spanos P. D., Solomos G. P., 1983. Markov approximation to transient vibration, Journal of Engineering Mechanics, vol. 109: 1134-1150. Wen Y. K., 1980. Equivalent linearization for hysteretic systems under random excitation, Journal of Applied Mechanics, vol. 47: 150154. Books Clough R. W., Penzien J., 1993. Dynamics of structures, McGraw-Hill. Ikhouane F., Rodellar J., 2007. Systems with hysteresis: analysis, identification and control using the Bouc-Wen model, John Wiley and Sons. Roberts J. B., Spanos P. D., 2003. Random Vibration and Statistical Linearization. New York: Dover Publications.

Reliability Assessment of Nonlinear MDOF Systems Subject to Evolutionary Stochastic Excitation I.P. Mitseas, I.A. Kougioumtzoglou, P. D. Spanos, and M. Beer

Edited Volumes Mitseas I.P., Kougioumtzoglou I.A., Beer M., Optimal design of nonlinear structures under evolutionary stochastic earthquake excitations, Proceedings of the Intnl Conf. on Engrg. and

12

Appl. Science Optimiz., OPT-i 2014, Kos, Greece, 2014.