"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 .
"De, Subhayan"@en .
"Wojtkiewicz, Steven F."@en .
"Johnson, Erik A."@en .
"2015-05-25T18:07:06Z"@en .
"2015-07"@en .
"This paper proposes a computationally efficient framework for design optimization under\nuncertainty for structures with local nonlinearities. To reduce the high computational cost of Monte\nCarlo simulation of such problems, an exact model reduction to a low-order Volterra integral equation\nis used to accelerate each simulation, and variance-reduced sampling is used to reduce the number of\nsimulations required for the uncertainty quantification. This optimization framework is applied to a\nbenchmark cable-stayed bridge problem, designing one pair of passive tuned mass dampers given a pair\nof uncertain passive power law dampers, providing significant gains in computational efficiency, two\norders of magnitude, compared to traditional approaches."@en .
"https://circle.library.ubc.ca/rest/handle/2429/53453?expand=metadata"@en .
"12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12Vancouver, Canada, July 12-15, 2015Efficient Optimal Design-Under-Uncertainty of Passive StructuralControl DevicesSubhayan DeGraduate Research Assistant, Sonny Astani Department of Civil and EnvironmentalEngineering, University of Southern California, Los Angeles, California, USASteven F. WojtkiewiczAssociate Professor, Department of Civil and Environmental Engineering,Clarkson University, Potsdam, New York, USAErik A. JohnsonProfessor, Sonny Astani Department of Civil and Environmental Engineering,University of Southern California, Los Angeles, California, USAABSTRACT: This paper proposes a computationally efficient framework for design optimization underuncertainty for structures with local nonlinearities. To reduce the high computational cost of MonteCarlo simulation of such problems, an exact model reduction to a low-order Volterra integral equationis used to accelerate each simulation, and variance-reduced sampling is used to reduce the number ofsimulations required for the uncertainty quantification. This optimization framework is applied to abenchmark cable-stayed bridge problem, designing one pair of passive tuned mass dampers given a pairof uncertain passive power law dampers, providing significant gains in computational efficiency, twoorders of magnitude, compared to traditional approaches.1. INTRODUCTIONRobust design in presence of uncertainties in ma-terial, loading or topological characteristics of astructure has been investigated using convex op-timization, neural network and evolutionary algo-rithms (Enevoldsen and S\u00C3\u00B8rensen, 1994; Sandgrenand Cameron, 2002; Papadrakakis and Lagaros,2002; Zang et al., 2005; Calafiore and Dabbene,2008). Over the past few decades, reliability-based design optimization of these uncertain struc-tures has been investigated using single and multi-objective optimization both approaches (Frangopol,1985; Gasser and Schu\u00C3\u00ABller, 1997; Tu et al., 1999;Adeli, 2002; Beck et al., 1996, 1999).While many design-level earthquakes generatesuperstructure responses that are elastic and linear,passive structural control devices embedded in thestructure often have nonlinear characteristics (e.g.,power-law, bilinear and/or hysteretic behavior), in-troducing local nonlinearities into an otherwise lin-ear model and, therefore, requiring either simplify-ing approximations or a full nonlinear simulationto evaluate response characteristics. If some partsof such a model are also uncertain (e.g., character-istics of the same, or other, structural control ele-ments or other localized components), then MonteCarlo sampling can be used for uncertainty charac-terization; this, in turn, further increases the compu-tational requirements for analysis of such systems.The optimal design of passive nonlinear control de-vices requires repeated solutions of the nonlinearsystem. Multiplying these three costs \u00E2\u0080\u0094 nonlin-ear simulation, Monte Carlo sampling, and designparameter iterations \u00E2\u0080\u0094 often creates a significantcomputational burden, so there is a clear advantagein the design-under-uncertainty problem for a com-putationally efficient approach to solve for the re-sponse of locally nonlinear systems.A conventional nonlinear solver (e.g., ode45 inMATLAB) cannot exploit the localized nature of the112th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12Vancouver, Canada, July 12-15, 2015nonlinearities. However, a method recently pro-posed by the last two authors (Gaurav et al., 2011)develops an exact nonlinear model reduction forsystems of this type, resulting in a low-order non-linear Volterra integral equation, providing a signif-icant computational speedup compared to ode45.Further, the authors recently proposed (Kamalzareet al., 2015) capitalizing on the local nature of de-sign variables to make design optimization simula-tions more computationally efficient.This paper extends the scope of that work by ad-ditionally incorporating localized uncertainties intothe model to arrive at a computationally-efficientdesign-under-uncertainty framework for systemsthat are mostly linear and deterministic but thathave localized nonlinear and uncertain elements.The proposed framework is illustrated using a fi-nite element model of the Bill Emerson MemorialBridge in Cape Girardeau, Missouri, consisting of579 nodes, 162 beam elements, 128 cable elements,420 rigid links and 134 nodal masses, resultingin a 419 degree-of-freedom (DOF) model after re-moving dependent or boundary DOFs (Dyke et al.,2003). The optimal design of parameters of lin-ear and nonlinear structural control devices is per-formed using the proposed approach while key el-ements of the structure are uncertain (e.g., a previ-ously installed passive damping device). The pro-posed method is shown to provide significant re-duction in required computation times, relative toode45, while delivering the same level of accuracy.2. METHODOLOGY2.1. Efficient response with local uncertaintiesFollowing Gaurav et al. (2011), let the nonlinearstructure model be given in state space byX\u00CB\u0099(t) = AX(t)+Bw(t)+Lugu(Xu(t);\u00CE\u00B4 )+Ldgd(Xd(t);\u00CE\u00B8 )= AX(t)+Bw(t)+Lg(X(t);\u00CE\u00B4 ,\u00CE\u00B8 )Y(t) = CX(t)+Dw(t)+Eugu(Xu(t);\u00CE\u00B4 )+Edgd(Xd(t);\u00CE\u00B8 )= CX(t)+Dw(t)+Eg(X(t);\u00CE\u00B4 ,\u00CE\u00B8 )X(0) = x0(1)where X(t) is the n\u00C3\u00971 state vector; A is then\u00C3\u0097n state matrix, w is an m\u00C3\u00971 external excita-tion; B is the n\u00C3\u0097m influence matrix; gu(\u00C2\u00B7; \u00C2\u00B7) isan ng,u \u00C3\u0097 1 function of a subset Xu(t) = GuX(t)of states and uncertain parameters \u00CE\u00B4 ; gd(\u00C2\u00B7; \u00C2\u00B7) is anng,d \u00C3\u0097 1 function of a subset Xd(t) = GdX(t) ofstates and design parameters \u00CE\u00B8 ; Lu is an n\u00C3\u0097ng,uinfluence matrix mapping to all states from theforce vector gu arising due to uncertain parame-ters; Ld is an n\u00C3\u0097ng,d influence matrix mappingto all states from the force vector gd(\u00C2\u00B7; \u00C2\u00B7) arisingdue to design parameters; D is an ny\u00C3\u0097m influencematrix; Eu and Ed are ny\u00C3\u0097ng,u and ny\u00C3\u0097ng,d in-fluence matrices, respectively; and x0 is the ini-tial condition. State subsets Xu(t) and Xd(t) areno,u\u00C3\u00971 and no,d\u00C3\u00971 vectors, respectively, withno,u,no,d \u001C n. Output Y(t) is an ny\u00C3\u00971 vec-tor. Combining the forces from the uncertaintiesand nonlinearities, L = [LTu LTd ]T is n\u00C3\u0097no andg(X(t);\u00CE\u00B4 ,\u00CE\u00B8 ) = [gTu (Xu(t);\u00CE\u00B4 ) gTd (Xd(t);\u00CE\u00B8 )]T isno\u00C3\u00971, where no \u00E2\u0089\u00A4 (ng,u + ng,d). Similarly, X(t) =[XTu (t) XTd (t)]T = GX(t) for G = [GTu GTd ]T.The nominal linear system corresponding to thenonlinear system in (1) isx\u00CB\u0099(t) = Ax(t)+Bw(t), x(0) = x0 (2)Using the principle of superposition, the responseof the system can be divided into two parts: the so-lution x(t) of the nominal linear system in (2) andthe contribution of the functions g(X) which can bewrittenx(t) =eAtx0 +\u00E2\u0088\u00AB t0HB(t\u00E2\u0088\u0092 s)w(s)ds,x(nl)(t) =\u00E2\u0088\u00AB t0HL(t\u00E2\u0088\u0092 s)g(X(s);\u00CE\u00B4 ,\u00CE\u00B8 )ds(3)where the impulse responses are given by HB(t) =eAtB and HL(t) = eAtL. The total response is thesuperposition of these two responses; i.e., X = x+x(nl). For any value of \u00CE\u00B4 and \u00CE\u00B8 , the solution of (3)can be computed efficiently with the following:p(t) = g(X(t);\u00CE\u00B4 ,\u00CE\u00B8 );X(t) = x(t)+\u00E2\u0088\u00AB t0HL(t\u00E2\u0088\u0092 s)p(s)ds(4)where x(t) = Gx(t) and HL(t) = GHL(t). The setof equations (4) can be combined into:p(t)\u00E2\u0088\u0092g(x+\u00E2\u0088\u00AB t0HL(t\u00E2\u0088\u0092 s)p(s)ds;\u00CE\u00B4 ,\u00CE\u00B8)= 0 (5)212th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12Vancouver, Canada, July 12-15, 2015Equation (5) is a nonlinear vector Volterra inte-gral equation (NVIE) written in nonstandard form,which can be solved by a Newton-Gregory integra-tion scheme (Linz, 1985). In the algorithm pro-posed in the authors\u00E2\u0080\u0099 earlier works, the Newton-Gregory integration is coupled with a recursive FFTformulation to drastically reduce the cost of compu-tation of the required convolution.Of course, any redundant columns in L (e.g., iftwo devices are collocated) can be eliminated andthe dimension of g(X(t);\u00CE\u00B4 ,\u00CE\u00B8 ) can be reduced, sim-plifying the solution of NVIE (5). Further, any re-dundant columns of G can be eliminated and thedimensions of x and X can be reduced.2.2. Design optimization under uncertaintyTwo design-under-uncertainty objectives are em-ployed herein: worst-case design and average de-sign (Calafiore and Dabbene, 2008). A brief de-scription of these two procedures follows.2.2.1. Worst case designIn this design method, the structure or control de-vice is designed for the case when the cost functionis maximized over the domain of uncertainty withconstraints satisfied, known as worst-case design.This design optimization problem can be defined asmin\u00CE\u00B8\u00E2\u0088\u0088\u00CE\u0098max\u00CE\u00B4\u00E2\u0088\u0088\u00E2\u0088\u0086J(Y(t);\u00CE\u00B8 ,\u00CE\u00B4 )subject to h(Y(t);\u00CE\u00B8 ,\u00CE\u00B4 max(\u00CE\u00B8 ))\u00E2\u0089\u00A4 0 a.s.1(6)where J(\u00C2\u00B7) and h(\u00C2\u00B7) may be functionals of the en-tire trajectory of Y; the set of all possible values ofdesign parameter \u00CE\u00B8 is denoted by \u00CE\u0098; \u00E2\u0088\u0086 representsa probability space {\u00E2\u0084\u00A6,P,F} with sample space\u00E2\u0084\u00A6, probability measure P, and \u00CF\u0083 -algebra F corre-sponding to the uncertainty defined for the problem;\u00CE\u00B4 max(\u00CE\u00B8 ) \u00E2\u0088\u0088 \u00E2\u0088\u0086 is the uncertainty realization which,for a particular design \u00CE\u00B8 , maximizes J(Y(t);\u00CE\u00B8 ,\u00CE\u00B4 )subject to the constraint h; and J(\u00C2\u00B7) is assumed con-cave in \u00CE\u00B4 . The analytical solution of (6) may not al-ways be possible. However, with samples {\u00CE\u00B4 i}N\u00CE\u00B4i=1from \u00E2\u0088\u0086, one may approximate the problem asmin\u00CE\u00B8\u00E2\u0088\u0088\u00CE\u0098maxi=1,...,N\u00CE\u00B4J(Y(t);\u00CE\u00B8 ,\u00CE\u00B4 i)subject to h(Y(t);\u00CE\u00B8 ,\u00CE\u00B4 i,max(\u00CE\u00B8 ))\u00E2\u0089\u00A4 0(7)1Almost sure (a.s.) event happens with probability 1.where \u00CE\u00B4 i,max(\u00CE\u00B8 ) corresponds to the sample which,for a particular design \u00CE\u00B8 , maximizes J(Y(t);\u00CE\u00B8 ,\u00CE\u00B4 )subject to the constraint h.2.2.2. Average designIn the second method of design considered, the ex-pected value of cost is minimized while keeping theconstraints satisfied. This optimization problem isdefined asmin\u00CE\u00B8\u00E2\u0088\u0088\u00CE\u0098E\u00CE\u00B4 [J(Y(t);\u00CE\u00B8 ,\u00CE\u00B4 )]subject to h(Y(t);\u00CE\u00B8 ,\u00CE\u00B4 )\u00E2\u0089\u00A4 0 a.s.(8)With samples {\u00CE\u00B4 i}N\u00CE\u00B4i=1 from \u00E2\u0088\u0086,min\u00CE\u00B8\u00E2\u0088\u0088\u00CE\u00981N\u00CE\u00B4N\u00CE\u00B4\u00E2\u0088\u0091i=1J(Y(t);\u00CE\u00B8 ,\u00CE\u00B4 i)subject to h(Y(t);\u00CE\u00B8 ,\u00CE\u00B4 i)\u00E2\u0089\u00A4 0(9)3. NUMERICAL EXAMPLE3.1. Cable-stayed bridge modelThe example used to demonstrate the framework isa numerical model of the Bill Emerson MemorialBridge, a cable stayed bridge, built in 2003 acrossthe Mississippi river between Cape Girardeau, Mis-souri, and East Cape Girardeau, Illinois. A finiteelement model of the bridge superstructure, devel-oped in Dyke et al. (2003) and shown in Fig. 1, con-sists of 579 nodes, 128 cable elements, 162 beamelements, 420 rigid links and 134 nodal masses.(The version of the model used herein has no con-nection between deck and tower except throughthe cables so as to allow energy dissipator devicesplaced between the deck and a tower.) The initial3474 degree-of-freedom (DOF) model, which de-scribes the superstructure\u00E2\u0080\u0099s linear motion about thestatic equilibrium, is reduced to 909 DOFs whenthe boundary conditions are imposed and the slaveDOFs removed. Static condensation is then ap-plied to eliminate DOFs with small contribution tothe global response, resulting in the final 419 DOFmodel, which is available publicly (Dyke et al.,2003) and which is used here. The equation of mo-tion of this final bridge model is given byMsu\u00C2\u00A8s(t)+Csu\u00CB\u0099s(t)+Ksus(t) =\u00E2\u0088\u0092Msru\u00C2\u00A8g(t) (10)312th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12Vancouver, Canada, July 12-15, 2015Figure 1: Finite element model of the bridge (dimensions in m); adapted from Dyke et al. (2003).where Ms, Cs, and Ks are the mass, damping andstiffness matrices, respectively, of the bridge super-structure for the active DOFs; r is the influence vec-tor for the ground acceleration (containing ones inentries corresponding to active displacements in thedirection of the 1-D horizontal excitation and zeroselsewhere); us(t) is the generalized displacementvector; and u\u00C2\u00A8g(t) is the ground acceleration in thelongitudinal direction. The reader is directed to thebenchmark definition paper (Dyke et al., 2003) forfurther details of this model.3.2. Passive damping devicesTwo types of passive devices are investigated: qunonlinear viscous dampers and qd nonlinear tunedmass dampers (TMDs). The damping forces in theviscous dampers follow the power-law relationf ui = cui |\u00E2\u0088\u0086u\u00CB\u0099i|\u00CE\u00B2 ui sgn(\u00E2\u0088\u0086u\u00CB\u0099i), i = 1, . . . ,qu (11)where \u00E2\u0088\u0086u\u00CB\u0099i is the velocity across the ith damper. Thegoverning differential equations of the TMD massmotions are given by,mdi v\u00C2\u00A8i + cdi |\u00E2\u0088\u0086\u00CB\u0099vi|\u00CE\u00B2 di sgn(\u00E2\u0088\u0086\u00CB\u0099vi)+ kdi \u00E2\u0088\u0086vi=\u00E2\u0088\u0092mdi rdi u\u00C2\u00A8g, i = 1, . . . ,qd(12)where vi is the displacement of the TMD relativeto the ground, \u00E2\u0088\u0086vi is the TMD displacement rela-tive to its attachment point on the bridge, and rdi isin [\u00E2\u0088\u00921,1] depending on the TMD\u00E2\u0080\u0099s orientation rel-ative to the earthquake ground motion direction.To be consistent with the symmetric nature ofthe bridge model, the passive devices are placedin identical pairs symmetrically located aboutthe bridge deck centerline; i.e., let cu2i\u00E2\u0088\u00921 = cu2iand \u00CE\u00B2 u2i\u00E2\u0088\u00921 = \u00CE\u00B2 u2i for i = 1, . . . ,qu (dampers) andcd2i\u00E2\u0088\u00921 = cd2i, \u00CE\u00B2 d2i\u00E2\u0088\u00921 = \u00CE\u00B2 d2i, md2i\u00E2\u0088\u00921 = md2i and rd2i\u00E2\u0088\u00921 = rd2ifor i= 1, . . . ,qd (TMDs). Further, this example con-siders each pair of nonlinear viscous dampers as un-certain \u00E2\u0080\u0094 \u00CE\u00B4 =[cu2 \u00CE\u00B2 u2 | cu4 \u00CE\u00B2 u4 | . . . | cuqu \u00CE\u00B2uqu]T\u00E2\u0080\u0094 and the stiffness and damping of each pairof TMDs as the design parameters \u00E2\u0080\u0094 \u00CE\u00B8 =[cd2 \u00CE\u00B2 d2 kd2 | cd4 \u00CE\u00B2 d4 kd4 | . . . | cdqd \u00CE\u00B2dqd kdqd]T.The ratio of TMD masses relative to bridge mass isfixed at a certain value.3.3. FormulationThe bridge equations of motion can be written inform of (1) where,X(t) ={u(t)u\u00CB\u0099(t)}, A =[0 I\u00E2\u0088\u0092M\u00E2\u0088\u00921K \u00E2\u0088\u0092M\u00E2\u0088\u00921C],B ={0r\u00C2\u00AF}, Lu =[0\u00E2\u0088\u0092M\u00E2\u0088\u00921Ru], Ld =[0\u00E2\u0088\u0092M\u00E2\u0088\u00921Rd],u(t) ={us(t)v(t)}, M=[Ms 00 Md], K=[Ks 00 0],and C =[Cs 00 0]; Md is a diagonal matrix ofTMD masses; r\u00C2\u00AF = [rT r1 . . .rqd]T is the influence412th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12Vancouver, Canada, July 12-15, 2015vector for the ground acceleration; Ru and Rd arethe influence matrices for the uncertain damperforces and the TMD forces (from design parame-ters), respectively. Each of the qu columns of Ruand qd columns of Rd transform the pseudoforce ofone passive device to a global force vector, withgu(Xu(t);\u00CE\u00B4 )=\u00EF\u00A3\u00AE\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00B0\u00E2\u0088\u0092cu2|\u00E2\u0088\u0086u\u00CB\u00991|\u00CE\u00B2 u2 sgn(\u00E2\u0088\u0086u\u00CB\u00991)\u00E2\u0088\u0092cu2|\u00E2\u0088\u0086u\u00CB\u00992|\u00CE\u00B2 u2 sgn(\u00E2\u0088\u0086u\u00CB\u00992)...\u00E2\u0088\u0092 cuqu|\u00E2\u0088\u0086u\u00CB\u0099qu\u00E2\u0088\u00921|\u00CE\u00B2 uqu sgn(\u00E2\u0088\u0086u\u00CB\u0099qu\u00E2\u0088\u00921)\u00E2\u0088\u0092cuqu|\u00E2\u0088\u0086u\u00CB\u0099qu|\u00CE\u00B2 uqu sgn(\u00E2\u0088\u0086u\u00CB\u0099qu)\u00EF\u00A3\u00B9\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BB,where Xu(t)=[\u00E2\u0088\u0086u\u00CB\u00991 \u00E2\u0088\u0086u\u00CB\u00992 | . . . | \u00E2\u0088\u0086u\u00CB\u0099qu\u00E2\u0088\u00921 \u00E2\u0088\u0086u\u00CB\u0099qu]T=GuX(t) =[0 RTu]X(t) andgd(Xd(t);\u00CE\u00B8 ) =\u00EF\u00A3\u00AE\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00AF\u00EF\u00A3\u00B0\u00E2\u0088\u0092cd2|\u00E2\u0088\u0086v\u00CB\u00991|\u00CE\u00B2 d2 sgn(\u00E2\u0088\u0086v\u00CB\u00991)\u00E2\u0088\u0092kd1\u00E2\u0088\u0086v\u00CB\u00991\u00E2\u0088\u0092cd2|\u00E2\u0088\u0086v\u00CB\u00992|\u00CE\u00B2 d2 sgn(\u00E2\u0088\u0086v\u00CB\u00992)\u00E2\u0088\u0092kd2\u00E2\u0088\u0086v\u00CB\u00992...\u00E2\u0088\u0092 cdqd|\u00E2\u0088\u0086v\u00CB\u0099qd\u00E2\u0088\u00921|\u00CE\u00B2 dqd sgn(\u00E2\u0088\u0086v\u00CB\u0099qd\u00E2\u0088\u00921)\u00E2\u0088\u0092kdqd\u00E2\u0088\u0086v\u00CB\u0099qd\u00E2\u0088\u00921\u00E2\u0088\u0092cdqd|\u00E2\u0088\u0086v\u00CB\u0099qd|\u00CE\u00B2 dqd sgn(\u00E2\u0088\u0086v\u00CB\u0099qd)\u00E2\u0088\u0092kdqd\u00E2\u0088\u0086v\u00CB\u0099qd\u00EF\u00A3\u00B9\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BA\u00EF\u00A3\u00BBwhere Xd(t)= [\u00E2\u0088\u0086v1 \u00E2\u0088\u0086v\u00CB\u00991 \u00E2\u0088\u0086v2 \u00E2\u0088\u0086v\u00CB\u00992 | . . . | \u00E2\u0088\u0086vqd\u00E2\u0088\u00921\u00E2\u0088\u0086v\u00CB\u0099qd\u00E2\u0088\u00921 \u00E2\u0088\u0086vqd \u00E2\u0088\u0086v\u00CB\u0099qd]T = GdX(t) = [0 RTd ]X(t).3.4. Objective function and constraintsThe goal of this example is to optimize the param-eters of the TMDs to improve the bridge\u00E2\u0080\u0099s perfor-mance given uncertainties in the nonlinear viscousdampers installed elsewhere in the bridge. Differentperformance metrics can be chosen; a set of met-rics, normalized with respect to the uncontrolledand connected deck-tower case suggested by Dykeet al. (2003), are used here subject to the 1940 ElCentro earthquake excitation. As in Dyke et al.(2003): Fbi(t) and Mbi(t) are the base shear andoverturning moment, respectively, at the ith towerat time t; Fdi(t) and Mdi(t) are the correspondingdeck-level shear and overturning moment; || \u00C2\u00B7 (t)||denotes the root mean square (i.e., two-norm overtime) response; (\u00C2\u00B7)max0(\u00C2\u00B7) denotes the maximum, overboth time and tower number, uncontrolled base ordeck-level shear or overturning moment at the baseor deck level; and ||(\u00C2\u00B7)0(\u00C2\u00B7)(t)|| denotes the maxi-mum, over tower number, time-normed uncontrol-lable shear or overturning moment at the base ordeck level. Then, eight of the metrics used to formthe objective function and some constraints are,J1 =maxi,t|Fbi(t)|Fmax0b, J2 =maxi,t|Fdi(t)|Fmax0d,J7 =maxi||Fbi(t)||||F0b(t)||, J8 =maxi||Fdi(t)||||F0d(t)||J3 =maxi|Mbi(t)|Mmax0b, J4 =maxi|Mdi(t)|Mmax0d,J9 =maxi,t||Mbi(t)||||M0b(t)||, J10 =maxi,t||Mdi(t)||||M0d(t)||(13)Additional metrics used in the constraints areJ6 = maxi,t\u00E2\u0088\u00A3\u00E2\u0088\u00A3\u00E2\u0088\u00A3\u00E2\u0088\u00A3xbi(t)x0b\u00E2\u0088\u00A3\u00E2\u0088\u00A3\u00E2\u0088\u00A3\u00E2\u0088\u00A3, J12 = maxi,tfi(t)W, (14)where xbi(t) is the displacement of the deck at thetwo ends (Bent 1 and Pier 4) in the finite elementmodel at time t; x0b is the maximum, over time andend locations, of the uncontrolled displacement ofthe deck; fi(t) is the amount of force exerted bythe ith device, assumed here to be the TMD; andW = 510 MN is the weight of the bridge super-structure. The normed responses are calculated us-ing 200 s of response as suggested by Dyke et al.(2003). (Note: performance metrics J5, J11 and J13through J18, defined in the benchmark (Dyke et al.,2003), are not used herein.)Using these performance metrics, the determin-istic optimization problem is formulated asmin\u00CE\u00B8\u00E2\u0088\u0088\u00CE\u0098J1(\u00CE\u00B8 ) (15)subject to Jk(\u00CE\u00B8 )\u00E2\u0089\u00A4 \u00CE\u00B1Jk,0 \u00E2\u0088\u0080k \u00E2\u0088\u0088 {2\u00E2\u0080\u00934,7\u00E2\u0080\u009310}where \u00CE\u0098 = {\u00CE\u00B8 : klb \u00E2\u0089\u00A4 \u00CE\u00B83 j\u00E2\u0088\u00922 \u00E2\u0089\u00A4 kub,clb \u00E2\u0089\u00A4 \u00CE\u00B83 j\u00E2\u0088\u00921 \u00E2\u0089\u00A4cub,\u00CE\u00B2 lb \u00E2\u0089\u00A4 \u00CE\u00B83 j \u00E2\u0089\u00A4 \u00CE\u00B2 ub, j = 1, . . . ,qd}; Jk,0 cor-responds to the kth performance metric withoutdampers (i.e. cdi = 0, kdi = 0); (\u00C2\u00B7)lb and (\u00C2\u00B7)ub arelower and upper bounds on the design parameters;and \u00CE\u00B1 is set to 1.25 (De et al., 2015). Each TMD isdesigned with a mass that is 2% of the bridge mass;i.e., mdi = 0.02W/g for i = 1, . . . ,qd.512th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12Vancouver, Canada, July 12-15, 2015DamperlocationsTMD locationsFigure 2: Finite element model (side view) of thebridge; adapted from Dyke et al. (2003).For simplicity, consider the case of a single pairof TMDs attached on top of the 2nd tower as shownin Figure 2. The initial values k0 and c0 of thedesign parameters for the optimization problemare taken from the guidelines specified in Hoanget al. (2008) for the first mode of vibration with\u00CE\u00B20 = 1.0 (clearly, any optimization should evalu-ate the effect of different initial values; here, withdifferent initial design points, the objective func-tion minimum was found to be relatively insensitivethough the optimum design point may change). Thelower and upper bounds for stiffness and dampingcoefficients are taken as [50 kN/m, 5 MN/m] and[100 kN\u00C2\u00B7(s/m)\u00CE\u00B2 , 6 MN\u00C2\u00B7(s/m)\u00CE\u00B2 ], respectively, afterpreliminary studies indicated that the optimum pa-rameters are expected to lie in these ranges. Thelower and upper limits for \u00CE\u00B2 are taken as 0.2 and1.8 as suggested by Main and Jones (2002).This optimization problem is then modified toinclude the design under uncertainty formulation.The uncertain parameter vector corresponds to asingle pair of nonlinear dampers (i.e., qu = 2), oneplaced connecting nodes 319 and 186 in the finiteelement model and the other symmetrically con-necting nodes 324 and 119 (see Figure 2). Theuncertain parameters of the passive dampers fol-low the distributions given in Table 1, where themean values were obtained from a deterministicoptimization similar to (15) performed over a setTable 1: Uncertain parameter descriptionVariable Distribution Mean COVDamping Log- 20.5440 0.05coeff.(cu2) normal MN\u00C2\u00B7(s/m)\u00CE\u00B2 u2Exponent Log- 0.9777 0.1(\u00CE\u00B2 u2 ) normalNote: COV = coefficient of variationof possible passive viscous damper locations andparameters with c \u00E2\u0088\u0088 [0.5,30] MN\u00C2\u00B7(s/m)\u00CE\u00B2 and \u00CE\u00B2 \u00E2\u0088\u0088[0.2,1.8].3.5. Optimization ResultsThe optimization is performed using MATLAB\u00E2\u0080\u0099sfmincon with default values for tolerances but lim-ited to 50 function evaluations (actual results mayslightly exceed 50 since each iteration may performmultiple function evaluations). An active-set algo-rithm (sequential quadratic programming) is usedinside fmincon. While fmincon can be providedwith gradient information, a finite difference ap-proximation is used here to evaluate the gradients.A preliminary study of this problem showed thatN\u00CE\u00B4 = 100 samples, using Latin Hypercube sam-pling method with antithetic variates variance re-duction (Ayyub and Lai, 1991), were sufficient togive converged results for the average design op-timization (though the number of samples clearlydepends on the application and the optimization ob-jective and constraints). The results of the designoptimization are shown in Tables 2 and 3. The op-timal values for kd2 are similar for both average andworst case design; however the optimal values forTable 2: Worst-case design of nonlinear TMDReduction of # fcn. evals.cost: max J1 (in %)26.4922 51kd2 cd2 \u00CE\u00B2 d2[MN/m] [MN \u00C2\u00B7 (s/m)\u00CE\u00B2d2 ]0.9086 3.8165 0.2000Table 3: Average design of nonlinear TMDReduction of # fcn. evals.cost: E [J1] (in %)7.3690 53kd2 cd2 \u00CE\u00B2 d2[MN/m] [MN \u00C2\u00B7 (s/m)\u00CE\u00B2d2 ]0.8453 1.9186 1.0826612th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12Vancouver, Canada, July 12-15, 20150 10 20 30 40 50 60\u00E2\u0088\u00920.2\u00E2\u0088\u00920.15\u00E2\u0088\u00920.1\u00E2\u0088\u00920.0500.050.10.150.2Time (s)Displacement (m) ode45 NVIE, \u00E2\u0088\u0086t = 0.46 ms(a) Accuracy in a typical generalized displacement0 20 40 60 80 100 120 140 160 180 20005101520253035Time (s)Cumulative sum of squared displacement (b) Cumulative sum of a typical squared responseFigure 3: Accuracy of the NVIE approach for first 60 s compared to MATLAB\u00E2\u0080\u0099s ode45cd2 and \u00CE\u00B2 d2 differ significantly. In the worst case de-sign, the value of device force performance met-ric J12 increases to 3.19 times that of the averagedesign case. A deterministic optimization for theTMD parameters with mean values for the passivedampers gives kopt = 0.9503 MN/m, copt = 0.3896MN\u00C2\u00B7(s/m)\u00CE\u00B2opt , \u00CE\u00B2opt = 0.7568, which gives a lowerpeak device force and 0.29 times J12 than averagedesign case. Hence, due to the assumed uncertaintyin the structure, the robust optimization requires aTMD capable of exerting greater force.3.6. Computational EfficiencyThe nonlinear Volterra integral equation approachinvolves a one-time cost and repeated costs for eachiteration of the design under uncertainty frame-work. The one-time computational cost involvessolution of x(t) in (3) and calculation of HL. Therepeated cost involves solution of (5) and for theresulting outputs Y(t) for different uncertainty anddesign parameters. To evaluate the computationalcost of the proposed method, it is compared witha fixed time step 4th order Runge-Kutta method(RK-4). The RK-4 method diverges for \u00E2\u0088\u0086t = 1 msor higher because of the numerical stiffness of thesystem of differential equations of motion. Hence,\u00E2\u0088\u0086t = 0.5 ms is used for RK-4, giving a relative ac-curacy of O(10\u00E2\u0088\u00923) computed relative to the ode45with relative tolerance of 10\u00E2\u0088\u00926 and absolute tol-erance of 10\u00E2\u0088\u00928 (for just the first 60 s of responsewhich provides RMS close to the long term asshown in Figure 3b; only 60 s is used because thefull ode45 solution takes too much time). The pro-posed NVIE approach gives relative accuracy inRMS of states in the same order with \u00E2\u0088\u0086t = 0.46 ms.For N\u00CE\u00B4 = 100, and for 10 function evaluations inthe optimization procedure, x(nl)(t) must be evalu-ated a total of 1000 times. Computation times areevaluated using MATLAB\u00E2\u0080\u0099s cputime function on acomputer with a 2.3 GHz Core i7-4850HQ proces-sor, 16 GB RAM, Mac OS X, and running MAT-LAB 2013a. The proposed NVIE approach takes13.53 cpu-min. to compute p(t), 2.16 cpu-min. tocompute Y(t) in each repeated calculation and atotal 10.92 cpu-days for 1000 simulations. How-ever RK-4 requires the full solution a total 1000times, where each simulation takes 77.66 cpu-hrs.,which projects to a total of 8.87 cpu-years for thefull optimization. Hence, the proposed design-under-uncertainty framework provides a computa-tional speedup of 296.33 with comparable accuracy.While only 3 design variables are considered here,comparable or increased gains in computational ef-ficiency, relative to other approaches for simulatingthe system responses for the function evaluations,are expected when more design variables are used.4. CONCLUSIONSStructural design under uncertainties providesmany computational challenges. The proposed712th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12Vancouver, Canada, July 12-15, 2015method shows significant computational advan-tages for such problems. In future work, the reli-ability based design optimization will be incorpo-rated in this framework; further, design optimiza-tion under uncertainty of structural systems, withpassive control devices, subjected to random exci-tation will be considered in subsequent studies.ACKNOWLEDGMENTSThe authors gratefully acknowledge the partial sup-port of this work by the National Science Foun-dation through awards CMMI 13-44937 and 14-36018. Any opinions, findings, and conclusionsor recommendations expressed in this material arethose of the authors and do not necessarily reflectthe views of the National Science Foundation. Theauthors also acknowledge the support of the first au-thor by the Provost\u00E2\u0080\u0099s Ph.D. Fellowship at the Uni-versity of Southern California.5. REFERENCESAdeli, H. (2002). Advances in design optimization. CRCPress.Ayyub, B. M. and Lai, K.-L. (1991). \u00E2\u0080\u009CSelective samplingin simulation-based reliability assessment.\u00E2\u0080\u009D Interna-tional Journal of Pressure Vessels and Piping, 46(2),229\u00E2\u0080\u0093249.Beck, J., Papadimitriou, C., Chan, E., and Irfanoglu, A.(1996). \u00E2\u0080\u009CReliability-based optimal design decisionsin the presence of seismic risk.\u00E2\u0080\u009D Proceedings 11thWorld Conference on Earthquake Engineering, Perg-amon.Beck, J. L., Chan, E., Irfanoglu, A., and Papadim-itriou, C. (1999). \u00E2\u0080\u009CMulti-criteria optimal structuraldesign under uncertainty.\u00E2\u0080\u009D Earthquake Engineeringand Structural Dynamics, 28(7), 741\u00E2\u0080\u0093762.Calafiore, G. C. and Dabbene, F. (2008). \u00E2\u0080\u009COptimiza-tion under uncertainty with applications to designof truss structures.\u00E2\u0080\u009D Structural and MultidisciplinaryOptimization, 35(3), 189\u00E2\u0080\u0093200.De, S., Wojtkiewicz, S. F., and Johnson, E. A. (2015).\u00E2\u0080\u009CComputationally efficient optimal design of passivecontrol devices for a benchmark cable-stayed bridge.\u00E2\u0080\u009DStructural Control and Health Monitoring, in prepa-ration.Dyke, S. J., Caicedo, J. M., Turan, G., Bergman, L., andHague, S. (2003). \u00E2\u0080\u009CPhase I benchmark control prob-lem for seismic response of cable-stayed bridges.\u00E2\u0080\u009DJournal of Structural Engineering, 129(7), 857\u00E2\u0080\u0093872.Enevoldsen, I. and S\u00C3\u00B8rensen, J. D. (1994). \u00E2\u0080\u009CReliability-based optimization in structural engineering.\u00E2\u0080\u009D Struc-tural Safety, 15(3), 169\u00E2\u0080\u0093196.Frangopol, D. M. (1985). \u00E2\u0080\u009CMulticriteria reliability-basedstructural optimization.\u00E2\u0080\u009D Structural Safety, 3(1), 23\u00E2\u0080\u009328.Gasser, M. and Schu\u00C3\u00ABller, G. I. (1997). \u00E2\u0080\u009CReliability-based optimization of structural systems.\u00E2\u0080\u009D Mathemat-ical Methods of Operations Research, 46(3), 287\u00E2\u0080\u0093307.Gaurav, Wojtkiewicz, S. F., and Johnson, E. A. (2011).\u00E2\u0080\u009CEfficient uncertainty quantification of dynamicalsystems with local nonlinearities and uncertainties.\u00E2\u0080\u009DProb. Eng. Mech., 26(4), 561\u00E2\u0080\u0093569.Hoang, N., Fujino, Y., and Warnitchai, P. (2008). \u00E2\u0080\u009COp-timal tuned mass damper for seismic applications andpractical design formulas.\u00E2\u0080\u009D Engineering Structures,30(3), 707\u00E2\u0080\u0093715.Kamalzare, M., Johnson, E. A., and Wojtkiewicz, S. F.(2015). \u00E2\u0080\u009CEfficient optimal design of passive structuralcontrol applied to isolator design.\u00E2\u0080\u009D Smart Structuresand Systems, in press.Linz, P. (1985). Analytical and numerical methods forVolterra equations, Vol. 7. SIAM.Main, J. and Jones, N. (2002). \u00E2\u0080\u009CFree vibrations of tautcable with attached damper. II: Nonlinear damper.\u00E2\u0080\u009DJournal of Engineering Mechanics, 128(10), 1072\u00E2\u0080\u00931081.Papadrakakis, M. and Lagaros, N. D. (2002).\u00E2\u0080\u009CReliability-based structural optimization usingneural networks and Monte Carlo simulation.\u00E2\u0080\u009DComputer Methods in Applied Mechanics andEngineering, 191(32), 3491\u00E2\u0080\u00933507.Sandgren, E. and Cameron, T. (2002). \u00E2\u0080\u009CRobust de-sign optimization of structures through considerationof variation.\u00E2\u0080\u009D Computers & Structures, 80(20), 1605\u00E2\u0080\u00931613.Tu, J., Choi, K. K., and Park, Y. H. (1999). \u00E2\u0080\u009CA new studyon reliability-based design optimization.\u00E2\u0080\u009D Journal ofMechanical Design, 121(4), 557\u00E2\u0080\u0093564.Zang, C., Friswell, M., and Mottershead, J. (2005). \u00E2\u0080\u009CAreview of robust optimal design and its application indynamics.\u00E2\u0080\u009D Computers & Structures, 83(4), 315\u00E2\u0080\u0093326.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.0076287"@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 .
"Graduate"@en .
"Efficient optimal design-under-uncertainty of passive structural control devices"@en .
"Text"@en .
"http://hdl.handle.net/2429/53453"@en .