"Non UBC"@en .
"DSpace"@en .
"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."@en .
"International Conference on Applications of Statistics and Probability (12th : 2015 : Vancouver, B.C.)"@en .
"Spence, Seymour M. J."@en .
"2015-05-22T16:57:50Z"@en .
"2015-07"@en .
"Recently a number of efficient reliability-based design optimization methodologies have\nbeen proposed for optimizing uncertain dynamic systems subject to stochastic excitation. While these\nmethods are capable of handling a large number of uncertain parameters, they are generally applicable to\nproblems characterized by small design variable vectors. This paper focuses on the development of a new\nreliability-based design optimization methodology for uncertain dynamic systems subject to stationary\nstochastic wind excitation that is capable of handling large design variable vectors, a characteristic of\nmany practical design problems, while considering system level performance constraints."@en .
"https://circle.library.ubc.ca/rest/handle/2429/53385?expand=metadata"@en .
"12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12Vancouver, Canada, July 12-15, 2015The SL-AVV Approach to System Level Reliability-Based DesignOptimization of Large Uncertain and Stochastic Dynamic SystemsSeymour M.J. SpenceAssistant Professor, Department of Civil Engineering and Environmental Engineering,University of Michigan, Ann Arbor, MI 48109, USAABSTRACT: Recently a number of efficient reliability-based design optimization methodologies havebeen proposed for optimizing uncertain dynamic systems subject to stochastic excitation. While thesemethods are capable of handling a large number of uncertain parameters, they are generally applicable toproblems characterized by small design variable vectors. This paper focuses on the development of a newreliability-based design optimization methodology for uncertain dynamic systems subject to stationarystochastic wind excitation that is capable of handling large design variable vectors, a characteristic ofmany practical design problems, while considering system level performance constraints.1. INTRODUCTIONReliability based design optimization (RBDO) isa powerful tool for obtaining structural systemsthat satisfy a number of probabilistic constraintsposed with the aim of ensuring a satisfactory per-formance of the system. The advantages of suchan approach over more classic deterministic op-timization strategies are well known. However,this approach is far more computationally involvedcompared to its deterministic counterpart. Thishas hindered the widespread adoption of RBDO,especially for large and dynamic structural sys-tems. This computational burden led to many ofthe early approaches to RBDO being based onclassic analytical approximations of the reliabil-ity integrals of the probabilistic constraints. Re-cent advances in computational capabilities as wellas achievements in the field of simulation-basedreliability assessment have spawned a new gen-eration of simulation-based RBDO strategies thatare capable of handling problems characterizedby large numbers of random variables, a typicalproperty of uncertain structural systems subject tostochastic excitation (Schu\u00C3\u00ABller and Jensen, 2008;Valdebenito and Schu\u00C3\u00ABller, 2010). This paper fo-cuses on the development of a novel simulation-based RBDO strategy aimed at solving problemsthat are not only characterized by a high number ofrandom variables, but also by a high-dimensionaldesign variable vector which, significantly com-plicates the situation (Valdebenito and Schu\u00C3\u00ABller,2010). This class of problems is often encoun-tered in structural design optimization where alarge number of members are to be designed(Spence and Kareem, 2014), or in applications oftopology optimization where a high number ofdesign variables are required in order to providean adequate discretization of the design domain(Bobby et al., 2014). In particular, the method pro-posed in this work is based on the generalizationof a recently proposed component-level simulation-based RBDO strategy (Spence and Kareem, 2014;Spence and Gioffr\u00C3\u00A8, 2012) to wind excited systemscharacterized by system level constraints.2. PROBLEM DEFINITIONThe RBDO problems of interest to this work my becast in the following form:Find x = {x1, ...,xm}T (1)to minimize W = f (x) (2)s. t. Pf (x) \u00E2\u0089\u00A4 P0 (3)xk \u00E2\u0088\u0088 Xk k = 1, ...,m (4)where x is a high-dimensional design variable vec-tor containing the parameters that fully define the112th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12Vancouver, Canada, July 12-15, 2015state of the system, e.g. section sizes, W is the ma-terial weight of the structural system, Pf is the sys-tem level failure probability, P0 is the target systemlevel failure probability while Xi is the discrete setto which the kth design variable must belong.The RBDO problem outlined above is character-ized, for the problems that are of interest to thisstudy, by a high-dimensional design variable vectorx, a high-dimensional uncertain vectorU describingthe model uncertainty as well as a multi-variate sta-tionary stochastic process, F, describing the windexcitation. In order to describe damage, the follow-ing system level demand to capacity ratio will beconsidered:d(u, r\u00CB\u0086,x) = maxi=1,...,Nr{maxt\u00E2\u0088\u0088[0,T ]|ri(t;u,x)|ci}(5)where Nr is the total number of components defin-ing the system response, T is the event duration,ri(t) are the component response processes, r\u00CB\u0086 isthe vector collecting the largest values of ri to oc-cur during a given event of duration T while ci isthe capacity of the system. Obviously ci is charac-terized by a significant amount of uncertainty andwill therefore be modeled as a random variable inthis work. Under these circumstances, a predefineddamage state will occur if d is larger than 1. There-fore the limit state function of interest is the follow-ing:g(u, r\u00CB\u0086,x) = 1\u00E2\u0088\u0092d(u, r\u00CB\u0086,x) (6)while the failure probability of interest is:Pf (x) =P(g(u, r\u00CB\u0086,x) \u00E2\u0089\u00A4 0)=\u00E2\u0088\u00AB \u00E2\u0088\u00ABg(u,r\u00CB\u0086,x)60p(r\u00CB\u0086|u)p(u)dr\u00CB\u0086du (7)where p indicates the conditional and non-conditional joint probability density functions ofthe uncertain vectors U and R\u00CB\u0086. As indicated inEq. (7), the random Nr-dimensional vector R\u00CB\u0086 willin general depend on U.3. MECHANICAL MODELING3.1. ExcitationIt is common to describe the intensity of windstorms through the maximum wind speed, vH tooccur during the event at a height of interest H(e.g. building height) averaged over a fixed intervalT (e.g. an hour) while considering a site specificroughness length z0. Generally, wind speed data, v,is only available averaged over a period \u00CF\u0084 (often 3 s)and collected at a meteorological height Hmet at re-gional airports characterized by a roughness lengthz01. A probabilistic model for transforming this in-formation into site specific data is the following:vH(T,z0) = e7e3(\u00CF\u0084,T )(e5z0e6z01)e4\u00CE\u00B4ln[H/(e5z0)]ln[Hmet/(e6z01)]e2e1v(\u00CF\u0084,Hmet ,z01)(8)where e1 and e2 are random variables modeling ob-servational and sampling errors in v; e4, e5, ande6 are random variables modeling the uncertaintieswith respect to the actual values of the empiricalconstant \u00CE\u00B4 = 0.0706 and of the roughness lengthsz0 and z01, respectively; e3(\u00CF\u0084,T ) is the conversionfactor that accounts for the uncertainty in convert-ing between wind speed averaging times; while e7is a model uncertainty to be used in the case of hur-ricanes and tornadoes.The stochastic wind loads can be estimated di-rectly from wind tunnel tests carried out on rigidscale models. In particular, each realization of themulti-variate stationary stochastic process definingthe wind loads may be related to a scaled (lengthand time scales) realization of its wind tunnel coun-terpart, fw(t), through the expression:f(t;u) = w1w2w3(vHvHm)2fw(t) (9)where vHm is the simulated hazard intensity used inthe wind tunnel tests while w1, w2 and w3 are com-ponents of U and model the uncertainties associ-ated with the estimation of building aerodynamicsthrough the use of wind tunnels.3.2. Component ResponseIn this work it is assumed that the response process,ri(t), associated with the ith failure mode of the sys-tem may be written in the following form:r(t;u) = s1\u00CE\u0093Tr[K\u00CE\u00A6nqrn(t)+ f(t)](10)212th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12Vancouver, Canada, July 12-15, 2015where \u00CE\u0093r is a vector of influence functions givingthe response in r due to a unit load acting at eachdegree of freedom of the system, \u00CE\u00A6n = [\u00CF\u00861, ...,\u00CF\u0086n]is the mass normalized mode shape matrix of ordern, qrn(t) = {qr1(t), ...,qrn(t)}T is the vector of reso-nant modal displacement responses, K is the nom-inal (mean) stiffness matrix while s1 is a randomvariable modeling the epistemic uncertainty in theload effect model of Eq. (10).In Eq. (10) each component of qrn(t) is given bythe solution of the following uncertain modal equa-tion:q\u00C2\u00A8 j(t)+2s3 j\u00CE\u00B6 js2 j\u00CF\u0089 jq\u00CB\u0099 j(t)+(s2 j\u00CF\u0089 j)2 q j(t) = \u00CF\u0086Tj f(t)(11)where q j, q\u00CB\u0099 j and q\u00C2\u00A8 j are the jth generalized displace-ment, velocity and acceleration response, \u00CF\u0089 j is thenominal value of the jth circular frequency, s2 j isan uncertain parameter modeling the variability inthe estimate of \u00CF\u0089 j, while \u00CE\u00B6 j is the nominal valueof the generalized damping ratio with s3 j an uncer-tain parameter modeling the uncertainty that existsin the nominal value of \u00CE\u00B6 j.The jth component of qrn(t) is simply givenby qr j(t) = q j(t)\u00E2\u0088\u0092 qb j(t), where the backgroundmodal displacement qb j is given by:qb j(t) =1(s2 j\u00CF\u0089 j)2\u00CF\u0086Tj f(t) (12)4. RELIABILITY PROBLEMAs outlined in Eq. (7), in order to calculate the fail-ure probability of the system, the conditional distri-butions of the largest values of the response func-tions ri(t) are needed. In particular, to this end, itis convenient to consider the following random re-duced variate:\u00CE\u00A8\u00CB\u0086i(u) =R\u00CB\u0086i(u)\u00E2\u0088\u0092\u00C2\u00B5ri(u)\u00CF\u0083ri(u)(13)where \u00C2\u00B5ri and \u00CF\u0083ri are the mean and standard devia-tion of the stationary response process ri(t) condi-tioned on u. For Gaussian systems (a typical char-acteristic of the stochastic response of multistorybuildings), the conditional distribution of \u00CF\u0088\u00CB\u0086i can beestimated from classic results of time-variant relia-bility:P(\u00CE\u00A8\u00CB\u0086i \u00E2\u0089\u00A4 \u00CF\u0088\u00CB\u0086i|u) = [1\u00E2\u0088\u0092G0(\u00CF\u0088\u00CB\u0086i)]exp[\u00E2\u0088\u0092\u00CF\u0085+(\u00CF\u0088\u00CB\u0086i,u)T1\u00E2\u0088\u0092G0(\u00CF\u0088\u00CB\u0086i)](14)where \u00CF\u0085+(\u00CF\u0088\u00CB\u0086i,u) is the up-crossing rate of thenormalized threshold \u00CF\u0088\u00CB\u0086i conditional on u whileG0(\u00CF\u0088\u00CB\u0086i,u) is the probability that the system responseis above the threshold \u00CF\u0088\u00CB\u0086i at time equal to zero andis given by:G0(\u00CF\u0088\u00CB\u0086i) = exp[\u00E2\u0088\u0092\u00CF\u0088\u00CB\u0086i2](15)The conditional crossing rate may be estimated as:\u00CF\u0085+(\u00CF\u0088\u00CB\u0086i,u) = \u00CE\u00BA(\u00CF\u0088\u00CB\u0086i,u)r+(\u00CF\u0088\u00CB\u0086i,u) (16)where r+ is given by the conditioned Rice formulaas:r+(\u00CF\u0088\u00CB\u0086i,u) =\u00CF\u0083r\u00CB\u0099i(u)2pi\u00CF\u0083ri(u)G0(\u00CF\u0088\u00CB\u0086i) (17)with \u00CF\u0083r\u00CB\u0099i the standard deviation of the derivative ofthe response process ri(t) while \u00CE\u00BA is the correc-tion factor accounting for any dependency betweensuccessive crossings of the normalized threshold \u00CF\u0088\u00CB\u0086iand can be modeled as:\u00CE\u00BA(\u00CF\u0088\u00CB\u0086i,u) = 1\u00E2\u0088\u0092 exp[\u00E2\u0088\u0092(k(u))1.2(2pi)0.1\u00CF\u0088\u00CB\u0086i](18)where k is the spectral shape factor given by:k(u) =\u00E2\u0088\u009A2pi(1\u00E2\u0088\u0092\u00CE\u00B321 (u)\u00CE\u00B30(u)\u00CE\u00B32(u))(19)where \u00CE\u00B3p for p = 0,1,2 are the spectral moments ofri given by:\u00CE\u00B3p(u) =\u00E2\u0088\u00AB\u00E2\u0088\u009E0\u00CF\u0089pSri(\u00CF\u0089;u)d\u00CF\u0089 (20)where Sri is the one-sided spectrum of ri while \u00CF\u0089 isthe circular frequency.312th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12Vancouver, Canada, July 12-15, 20154.1. Conditional Response StatisticsIn defining the distributions of the reduced vari-ates \u00CF\u0088\u00CB\u0086i, the conditional (on u) response statisticsare needed. To this end, the mean conditioned re-sponse, \u00C2\u00B5ri , is simply given by the expected valueof Eq. (10) given u and therefore by:\u00C2\u00B5ri(u) = s1\u00CE\u0093Tri f\u00C2\u00AF (21)where f\u00C2\u00AF is the expected value of f. The second or-der conditioned response statistic, \u00CF\u00832ri , may be esti-mated in the frequency domain as:\u00CF\u00832ri(u) = \u00CF\u00832rib(u)+\u00E2\u0088\u00AB\u00E2\u0088\u009E0Srrn (\u00CF\u0089 ;u)d\u00CF\u0089 (22)where \u00CF\u00832rib is the background response variancewhile Srrn is the one-sided resonant response spec-trum estimated considering the participation of nvibration modes. The estimate of the backgroundcontribution to \u00CF\u00832r is given by:\u00CF\u00832rib(u) = s21\u00CE\u0093TriCf(u)\u00CE\u0093ri (23)where Cf(u) is the covariance matrix of the excita-tion f conditioned on u. The resonant response con-tribution to \u00CF\u00832ri can be efficiently estimated througha double modal spectral proper orthogonal de-composition (POD) of f(t) (Carassale et al., 2001;Spence and Kareem, 2013). Following this frame-work, Srrn is estimated as:Srrn (\u00CF\u0089;u) =l\u00E2\u0088\u0091i=1F\u00CB\u009Ci(\u00CF\u0089;u)F\u00CB\u009C\u00E2\u0088\u0097i (\u00CF\u0089;u) (24)where the symbol \u00E2\u0088\u0097 indicates the transposed com-plex conjugate, l is the number of spectral load-ing modes used in the estimation of Srrn while F\u00CB\u009Ciis given by:F\u00CB\u009Ci(\u00CF\u0089 ;u) = s1\u00CE\u0093TriK\u00CE\u00A6nHrn(\u00CF\u0089 ;u)\u00CE\u00A6Tn \u00CF\u0087 i(\u00CF\u0089)\u00E2\u0088\u009A\u00CE\u009Bi(\u00CF\u0089)(25)where \u00CF\u0087 i is the ith frequency dependent load-ing eigenvector, \u00CE\u009Bi is the corresponding ith fre-quency dependent loading eigenvalue while Hrn isthe diagonal resonant mechanical transfer function(Spence and Kareem, 2013).The other spectral moments, i.e. \u00CE\u00B31 and \u00CE\u00B32 =\u00CF\u00832r\u00CB\u0099i , may be efficiently estimated through Eq. (20)where Sri is assessed in a strictly modal setting bysimply substituting into F\u00CB\u009Ci (Eq. 25) the classic me-chanical transfer function instead of Hrn .4.2. Solution StrategyIn defining an efficient solution strategy for the reli-ability integral Eq. (7), it is convenient to write theintegral in terms of the vector of reduced variates,\u00CE\u00A8\u00CB\u0086 = {\u00CF\u0088\u00CB\u00861, ..., \u00CF\u0088\u00CB\u0086Nr}T :Pf (x) =\u00E2\u0088\u00AB \u00E2\u0088\u00ABg(u,\u00CF\u0088\u00CB\u0086,x)60p(\u00CF\u0088\u00CB\u0086|u)p(u)d\u00CF\u0088\u00CB\u0086du (26)For the problems of interest to this work, the aboveintegral will be of high dimensions therefore rul-ing out the use of classical reliability methods(Schu\u00C3\u00ABller et al., 2003). In particular, in this workMonte Carlo simulation is adopted. In generatingconditional samples of the random vector \u00CE\u00A8\u00CB\u0086, thecomponents are considered independent and there-fore fully described by their marginal distributionsgiven in Eq. (14). The grounds for making thisassumption lies in the fact that the dependency be-tween the components of \u00CE\u00A8\u00CB\u0086models the dependencybetween the peaks of the normalized response pro-cesses, which are only weakly dependent on uthrough the crossing rate. Therefore \u00CE\u00A8\u00CB\u0086 will be onlyweakly dependent on parameters such as the meanwind speed vH defining the intensity of the windevent, which significantly contributes to the depen-dency of the peaks of the various response functionsri(t).5. THE SYSTEM LEVEL RBDO APPROACHThis section presents a new system level RBDOalgorithm, SL-AVV, that leverages the recentlyintroduced concept of Auxiliary Variable Vector(Spence and Kareem, 2014).5.1. Problem DefinitionThe failure probability of Eq. (7) is a system fail-ure probability defined as the union of the followingcomponent level failure events:Fi ={|r\u00CB\u0086i|ci> 1}i = 1, ...,Nr (27)412th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12Vancouver, Canada, July 12-15, 2015Here the failure probability of Eq. (7) is first writtenin terms of the component failure modes as:Pf = P{Nr\u00E2\u008B\u0083i=1Fi}=Nr\u00E2\u0088\u0091i=1P{Fi}\u00E2\u0088\u0092\u00CE\u00B1(Nr\u00E2\u0088\u0091i=1P{Fi})(28)where \u00CE\u00B1 is simply given by:\u00CE\u00B1 =(Nr\u00E2\u0088\u00921\u00E2\u0088\u0091i=1Nr\u00E2\u0088\u0091j>iP{Fi\u00E2\u0088\u00A9Fj}\u00E2\u0088\u0092Nr\u00E2\u0088\u00922\u00E2\u0088\u0091i=1Nr\u00E2\u0088\u00921\u00E2\u0088\u0091j>iNr\u00E2\u0088\u0091k> jP{Fi\u00E2\u0088\u00A9Fj \u00E2\u0088\u00A9Fk}+ . . .\u00E2\u0088\u0092 (\u00E2\u0088\u00921)Nr\u00E2\u0088\u00921P{Fi\u00E2\u0088\u00A9 ...\u00E2\u0088\u00A9FNr})/(Nr\u00E2\u0088\u0091i=1P{Fi})(29)It is then assumed that the component failure prob-abilities can be approximated as exponential dis-tributions. This allows the failure probability tobe written as (where the dependency on the designvariables is now indicated):Pf (x) = (1\u00E2\u0088\u0092\u00CE\u00B1(x))(Nr\u00E2\u0088\u0091i=1exp[\u00E2\u0088\u0092 1\u00C2\u00B5Di(x)])+\u00E2\u0088\u0086(x)(30)where \u00C2\u00B5Di is the expected value of the random com-ponent demand to capacity ratio given by R\u00CB\u0086i/Ciwhile \u00E2\u0088\u0086 is the error term introduced to account forthe assumption made on the component failure dis-tributions. If it is assumed that \u00CE\u00B1 and \u00E2\u0088\u0086 of Eq.(30) are independent of the design variable vectorx, then the dependency of the system level fail-ure probability on x is exclusively in terms of themean values of the random component demand tocapacity ratios Di. The problem therefore becomesthe description of the dependency of \u00C2\u00B5Di on x. Tothis end the concept of AVV (Spence and Kareem,2014) can be leveraged.Before continuing, it should be observed that theassumption of independence of \u00CE\u00B1 from x is equiv-alent to assuming that the change in the probabilityof the joint occurrence of the various failure modesdue to a change in x, will in general follow the trendof the sum of the probabilities of the individual fail-ure events.5.2. The AVVsIn order to derive the aforementioned AVVs, it isfirst convenient to consider the following variablesdefined for each realization of U and \u00CE\u00A8\u00CB\u0086 and for thecurrent design variable vector x0:\u00CF\u0092i(u, \u00CF\u0088\u00CB\u0086i,x0) =\u00CF\u0088\u00CB\u0086i(u,x0)CL(u,x0)\u00CE\u0093ri(x0)\u00CF\u0083ri(u,x0)(31)where CL is the covariance matrix of the followingvector:L(t;u,x0) = s1[K\u00CE\u00A6nqrn(t)+ f(t)](32)Dividing by the value of the capacity contained inu, the following static relationship is determined forthe damage ratio di associated with the realizationsof u and \u00CF\u0088\u00CB\u0086 :di(u,x0) =1ci\u00CE\u0093Tri(x0)\u00CF\u0092i(u,x0) = \u00CE\u0093Tri(x0)\u00CF\u0092di(u,x0)(33)where \u00CF\u0092di is simply given by the ratio between \u00CF\u0092iand ci. The realizations of the vector \u00CF\u0092d j , generatedduring the simulation process used to estimate thereliability of the system, may be used to define thefollowing AVV:\u00CF\u0092\u00CB\u009Ci(x0) = \u00CF\u0092\u00C2\u00AFdi(x0) (34)where \u00CF\u0092\u00C2\u00AFdi is the expected value of \u00CF\u0092di .The significance of the AVV, and so of \u00CF\u0092\u00CB\u009Ci, is thatif it is statically applied to the nominal structure itwill cause a response in ri that is equal to the ex-pected value of the damage ratio Di. In other wordsthe following relationship holds:\u00C2\u00B5Di(x0) = \u00CE\u0093Tri(x0)\u00CF\u0092\u00CB\u009Ci(x0) (35)This relationship is particularly useful as it allowsthe failure probability of Eq. (30) to be written interms of what may be considered, for all intents andpurposes, a series of static load distributions.5.3. The Approximate SubproblemThe definition \u00CF\u0092\u00CB\u009Ci together with the formulation ofthe failure probability given in Eq. (30) can be usedto define a approximate optimization subproblem,the sequential definition and solution of which will512th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12Vancouver, Canada, July 12-15, 2015lead to a final optimal solution. To this end, if theassumption of independence of \u00CE\u00B1 and \u00E2\u0088\u0086 from thedesign variable vector is extended to the AVV \u00CF\u0092\u00CB\u009Ci,then the following optimization problem can be castthat is completely defined from the results of a sin-gle reliability analysis carried out in the current de-sign point x0:Find x = {x1, ...,xm}T (36)to minimize W = f (x) (37)s. t.(1\u00E2\u0088\u0092\u00CE\u00B1)(Nr\u00E2\u0088\u0091i=1exp[\u00E2\u0088\u0092 1\u00CE\u0093Tri(x)\u00CF\u0092\u00CB\u009Ci])+\u00E2\u0088\u00866 P0(38)xk \u00E2\u0088\u0088 Xk k = 1, ...,m (39)Not only does the subproblem outlined above de-couple the reliability problem from the optimiza-tion loop, but it also takes on an extremely conve-nient form as it can be easily made explicit in termsof x by simply defining an explicit expression of \u00CE\u0093riin terms of x. In defining this explicit relationship, anumber of classic results can be used that were de-veloped for the optimization of deterministic struc-tures under static loads. In particular, here the solu-tion proposed by Chan et al. (1995) was adopted.Once the subproblem has been made explicit inx, any gradient-based optimization strategy can beused to find solutions. Here the pseudo-discrete op-timality criteria (Chan et al., 1995) is used, there-fore allowing the discrete nature of the design spaceto be fully considered.Once a solution is found to a subproblem definedin x0, a new sub-problem must be formulated in thedesign point found from the solution of the previoussub-problem. This updating procedure defines a de-sign cycle and will continue until the design pointof two successive design cycles are identical.6. CASE STUDY6.1. DescriptionThe case study consists of a 45-story rectangu-lar building with an offset core (Fig. 1a). Thecolumns consist of steel box sections and aregrouped in plan as indicated in Figure 1b (C1 toC18). In particular, it is required that the mid-line diameter of the box sections, Di, belong to theHV YX(a)180 mXY(b)BracingB1B2B2B3B3 B4B5B6C1C1C2C2C3C4C5C6C7C8C7C6C5C4C3C8C9 C10C11C12C13C12C11C10C9C17C18C17C15C14C14C15C16C1672 m 36 mDesignColumnFigure 1: (a) 3D view of the 45-story building; (b)structural layout of the building showing beam andcolumn assignments and the critical column.discrete set {0.2 m,0.25 m,0.3 m,0.35 m, ...,2 m}while the flange thickness is fixed at Di/20. Thecolumns are grouped over three consecutive floors.The beams, B1 to B6, are grouped in plan as in-dicated in Figure 1b and are required to belong tothe family of AISC W24 steel profiles. The beamsare also grouped three floors at a time (for a total of6\u00C3\u0097 15 groups). The diagonal bracings are also re-quired to belong to the AISC W24 steel profiles andare grouped as pairs over the height of the building.The aforementioned grouping results in 375 inde-pendent design variables. Initially the structure isdesigned with all columns having a mid-line diame-ter equal to 0.6 m while all beams and diagonals areset to W24\u00C3\u0097176 profiles. The first three structuralmodes, with mean circular frequencies \u00CF\u00891 = 1.023rad/s, \u00CF\u00892 = 1.118 rad/s and \u00CF\u00893 = 1.847 rad/s, areconsidered sufficient for describing the resonant re-sponse. The mean modal damping ratios were takenas 1.5 %. The uncertain parameters S1, S2 j and S3 jwith j = 1, ...,3 were modeled as independent log-normal random variables with coefficients of varia-tion 0.025, 0.015, and 0.3 respectively.Non-structural damage is to be controlled forwind blowing down the X direction (Fig. 1). In par-ticular, it is assumed that damage is associated withthe X and Y interstory drift response of the criti-cal column line illustrated in Fig. 1b. The capaci-ties,Ci, are taken as independent lognormal random612th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12Vancouver, Canada, July 12-15, 2015variables with mean 1/400 of the story height andcoefficient of variation 0.25.The building is considered to be located inthe Miami area of Florida, USA. The extremewind hazard is given by the risk of hurricanes. Inparticular, the site of the building is characterizedby a roughness length z0 = 2 m while an averagingtime of T = 3600 s was considered. The set of windspeeds used to defined the distribution of v (Eq.(8)) was obtained from the simulated hurricanedatabase of the National Institute of Standardsand Technology while considering milepost 1450where a roughness length of z0 = 0.05 m wasassumed together with a meteorological height ofHmet = 10 m. The averaging time for this datawas \u00CF\u0084 = 60 s. The distributions and parametervalues assumed for E1 to E7 can be found in Table2 of Spence and Kareem (2014) while consideringa mean for E3 of 0.8065. The stochastic loadsacting on the building were derived from the TokyoPolytechnic University Wind Pressure Database.In particular, the mean wind speed during thetests was vHm = 11 m/s. Through integrationand appropriate scaling, realizations of fw wereobtained. The distributions and parameter valuesassumed for W1 to W3 can be found in Table 2 ofSpence and Kareem (2014) while considering acoefficient of variation of 0.05 for W1. For the casestudy under investigation the reliability integralneeds to be calculated over the space of U ={C1, ...,C90,S1,ST2n,ST3n ,W1, ...,W3,E1, ....,E7,V}T ,with n = 3, and \u00CE\u00A8\u00CB\u0086, resulting in a total dimensionof 198 which constitutes a large scale reliabilityproblem. Considering the size of x, m = 375,the problem truly represents a large scale RBDOoptimization problem.6.2. ResultsTwo cases are presented in this section with thefirst considering the components of random vec-tor \u00CE\u00A8\u00CB\u0086 independent and the second considering thecomponents of \u00CE\u00A8\u00CB\u0086 fully correlated with the aim ofquantifying the assumption of Sec. 4.2. P0 wasset at 2\u00C3\u0097 10\u00E2\u0088\u00922. Figure 2 shows the convergencehistory of the objective function which here coin-cides with the weight of the material composing thestructural system while Fig. 3 reports the history ofDesign Cycle0 2 4 6 8 10W[kN]\u00C3\u009710511.21.41.61.8UncorrelatedCorrelatedFigure 2: Design history of the objective function.Design Cycle0 2 4 6 8 10P f10-210-1100UncorrelatedCorrelatedFigure 3: Failure probability history.the system level failure probability. From Figs.2 and 3 it is immediately evident the strong con-vergence properties of the proposed SL-AVV algo-rithm. Indeed, the problem practically convergesafter only 6 design cycles which means that thesimulation-based reliability analysis only had to beinvoked 6 times before solutions are found. It isalso evident that the assumption that the randomvector \u00CE\u00A8\u00CB\u0086 has independent components would, atleast for this example, seem to be perfectly reason-able with little difference between the two extremecases. Figure 4 reports the variation of \u00CE\u00B1 (Eq. (29))during the optimization process. While \u00CE\u00B1 is seento change, its variation quickly drops off after de-sign cycle 4 therefore allowing for the quick over-all convergence seen in Figs 2 and 3. Fig. 5 showsthe initial and final failure distributions for the twocases. Once again the effectiveness of the proposed712th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12Vancouver, Canada, July 12-15, 2015Design Cycle0 2 4 6 8 10\u00CE\u00B100.20.40.60.81UncorrelatedCorrelatedFigure 4: History of the parameter \u00CE\u00B1 (Eq. (29)).D10-1 100 101P f10-210-1100Initial uncorr.Initial corr.Final uncorr.Final corr.Figure 5: Initial and final failure probability distribu-tions.method is evident. Finally, it should be observedthat the failure probability of the converged designis exact as for any design cycle in which x does notchange (e.g. at convergence), all approximationsused in the proposed approach are exact.7. CONCLUSIONSThis paper presented a new simulation-based sys-tem level reliability based design optimization al-gorithm, SL-AVV, for large scale uncertain anddynamic systems excited by stationary stochasticwind excitation. The method was based on theconcept of decoupling the reliability analysis fromthe optimization loop through the definition of asequence of high quality subproblems defined interms of a number of Auxiliary Variable Vectors(AVVs) that allow the system failure probability tobe written in terms of a sum of exponential func-tions modeling the failure modes of the system. Thestructure of the AVVs allows the subproblem to beeasily made explicit in the design variable vector,and therefore any gradient-based optimization al-gorithm may be used for its resolution. The effec-tiveness of the proposed approach lies in how eachsubproblem can be defined exclusively in terms of asingle simulation-based reliability analysis carriedout in the current design point. The practicalityand strong convergence properties of the proposedapproach were illustrated on a full scale buildingexample characterized by a high-dimensional de-sign variable vector as well as a high-dimensionalsystem-level reliability integral.8. REFERENCESBobby, S., Spence, S. M. J., Bernardini, E., and Kareem,A. (2014). \u00E2\u0080\u009CPerformance-based topology optimiza-tion for wind-excited tall buildings: A framework.\u00E2\u0080\u009DEngineering Structures, 74, 242\u00E2\u0080\u0093255.Carassale, L., Piccardo, G., and Solari, G. (2001). \u00E2\u0080\u009CDou-ble modal transformation and wind engineering appli-cations.\u00E2\u0080\u009D J. Eng. Mech., 127(5), 432\u00E2\u0080\u0093438.Chan, C. M., Grierson, D. E., and Sherbourne, A. N.(1995). \u00E2\u0080\u009CAutomatic optimal design of tall steel build-ing frameworks.\u00E2\u0080\u009D J. Struct. Eng., 121(5), 838\u00E2\u0080\u0093847.Schu\u00C3\u00ABller, G. I. and Jensen, H. A. (2008). \u00E2\u0080\u009CComputa-tional methods in optimization considering uncertain-ties - an overview.\u00E2\u0080\u009DComput. Methods Appl. Mech. En-grg., 198(1), 2\u00E2\u0080\u009313.Schu\u00C3\u00ABller, G. I., Pradlwarter, H. J., and Koutsoure-lakis, P. S. (2003). \u00E2\u0080\u009CA comparative study of reliabil-ity estimation procedures for high dimensions.\u00E2\u0080\u009D 16thASCE Engineering Mechanics Conference, July 16-18, Seattle, Washington, USA.Spence, S. M. J. and Gioffr\u00C3\u00A8, M. (2012). \u00E2\u0080\u009CLarge scalereliability-based design optimization of wind excitedtall buildings.\u00E2\u0080\u009D Prob. Eng. Mech., 28, 206\u00E2\u0080\u0093215.Spence, S. M. J. and Kareem, A. (2013). \u00E2\u0080\u009CData-enableddesign and optimization (DEDOpt): Tall steel build-ings frameworks.\u00E2\u0080\u009D Comput. Struct., 134(12), 134\u00E2\u0080\u0093147.Spence, S. M. J. and Kareem, A. (2014). \u00E2\u0080\u009CPerformance-based design and optimization of uncertain wind-excited dynamic building systems.\u00E2\u0080\u009D Eng. Struct., 78,133\u00E2\u0080\u0093144.Valdebenito, M. A. and Schu\u00C3\u00ABller, G. I. (2010). \u00E2\u0080\u009CA sur-vey on approaches for reliability-based optimization.\u00E2\u0080\u009DStruct. Multidisc. Optim., 42, 645\u00E2\u0080\u0093663.8"@en .
"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."@en .
"Conference Paper"@en .
"10.14288/1.0076161"@en .
"eng"@en .
"Unreviewed"@en .
"Vancouver : University of British Columbia Library"@en .
"Attribution-NonCommercial-NoDerivs 2.5 Canada"@en .
"http://creativecommons.org/licenses/by-nc-nd/2.5/ca/"@en .
"Faculty"@en .
"The SL-AVV approach to system level reliability-based design optimization of large uncertain and stochastic dynamic systems"@en .
"Text"@en .
"http://hdl.handle.net/2429/53385"@en .