SEISMIC DRIFT PERFORMANCE-BASED DESIGN OPTIMIZATION OF REINFORCED CONCRETE BUILDINGS

Size: px
Start display at page:

Download "SEISMIC DRIFT PERFORMANCE-BASED DESIGN OPTIMIZATION OF REINFORCED CONCRETE BUILDINGS"

Transcription

1 3 th World Conference on Earthqake Engneerng Vancover, B.C., Canada Agst -6, 2004 Paper No. 223 SEISMIC DRIFT PERFORMANCE-BASED DESIGN OPTIMIZATION OF REINFORCED CONCRETE BILDINGS Xao-Kang ZO and Chn-Man CHAN 2 SMMARY Performance-based desgn s a modern approach to sesmc engneerng, n whch the desgn am s to delver a strctre capable of meetng certan predctable performance objectves nder dfferent levels of earthqake motons. Performance-based desgn sng nonlnear pshover analyss, whch generally nvolves tedos and ntensve comptatonal effort, s a hghly teratve process to meet desgnerspecfed and code reqrements. Ths paper presents an effectve compter-based technqe that ncorporates pshover analyss together wth nmercal optmzaton procedres to atomate the pshover drft performance desgn of renforced concrete (RC) bldngs. Steel renforcement, as compared wth concrete materals, appears to be the more cost-effectve materal that can be effectvely sed to control drft beyond the occrrence of frst yeldng and to provde the reqred dctlty of RC bldng frameworks. In ths stdy, steel renforcement ratos are taken as desgn varables drng the desgn optmzaton process. sng the prncple of vrtal work, the nonlnear nelastc sesmc drft responses generated by the pshover analyss can be explctly expressed n terms of steel renforcement desgn varables. An Optmalty Crtera technqe s presented n ths paper to solve the explct performancebased sesmc desgn optmzaton problem for RC bldngs. One bldng frame example s presented to llstrate the effectveness and practcalty of the proposed optmal desgn method. INTRODCTION The concept of performance-based desgn has become the ftre drecton of sesmc desgn codes [-3]. In the newly developed performance-based sesmc desgn approach, nonlnear analyss procedres become mportant n dentfyng the patterns and levels of damage to assess a strctre s nelastc behavor and to nderstand the modes of falre of the strctre drng severe sesmc events. Pshover analyss s a smplfed, statc, nonlnear procedre n whch a predefned pattern of earthqake loads s appled ncrementally to framework strctres ntl a plastc collapse mechansm s reached. Ths analyss method generally adopts a lmped-plastcty approach that tracks the spread of nelastcty throgh the formaton of nonlnear plastc hnges at the frame element s ends drng the ncremental loadng process. In general, the determnaton of the satsfactory performance response that flflls both the system level Research Assocate. Department of Cvl Engneerng, Hong Kong nversty of Scence and Technology, Hong Kong, P.R. Chna. E-mal: cezxk@st.hk. 2 Assocate Professor. Department of Cvl Engneerng, Hong Kong nversty of Scence and Technology, Hong Kong, P.R. Chna. E-mal: cecmchan@st.hk.

2 response and element level response reqres a hghly teratve tral-and-error desgn procedre even wth the ad of today s engneerng compter software. It has been recognzed that the nterstory drft performance of a mltstory bldng s an mportant measre of strctral and non-strctral damage of the bldng nder varos levels of earthqake moton [4]. In performance-based desgn, nterstory drft performance has become a prncpal desgn consderaton [,2]. The system performance levels of a mltstory bldng are evalated based on the nterstory drft vales along the heght of the bldng nder dfferent levels of earthqake moton [5]. The control of nterstory drft can also be consdered as a means to provde nform dctlty over all stores of the bldng. A large story drft may reslt n the occrrence of a weak story that may case catastrophc bldng collapse n a sesmc event. Therefore, nform story dctlty over all stores for a mltstory bldng s sally desred n sesmc desgn [6]. Althogh lateral drft performance s a prncpal concern n the sesmc desgn of strctres, economcally desgnng elements of bldng strctres for varos levels of elastc and nelastc lateral drft performance nder mltple levels of earthqake load s generally a rather dffclt and challengng task. Lateral drft desgn reqres the consderaton of a proper dstrbton of the stffness of all strctral elements and, n a severe sesmc event, also the occrrence and redstrbton of plastcty n the strctral elements. Strctral engneers are faced wth the problem of effcently proportonng strctral materals throghot the bldng to lmt the nelastc sesmc drft responses of a strctre. De to the lack of an atomated optmzaton technqe, performance-based sesmc drft desgn s sally carred ot by tral-and-error methods based on ntton and experence. Chan [7] developed an effcent compterbased optmzaton technqe for lateral stffness desgn of tall bldngs. Althogh ths research has reslted n actal applcatons to nmeros notable tall bldng projects n Hong Kong, t shold be noted that the research has been prmarly focsed on the elastc wnd drft performance of tall bldngs. Mch effort s stll needed to extend the crrent optmzaton technqe to nelastc sesmc desgn of mlt-story bldngs. Ths paper presents an effectve optmzaton technqe for the nelastc drft performance desgn of RC bldng frames nder pshover loadng. Attempts have been made to atomate the performancebased sesmc desgn of RC bldngs sng an optmzaton procedre. The qanttes of steel renforcement, the only effectve materal that provdes dctlty to RC bldng frameworks, are consdered as desgn varables n the nelastc sesmc drft optmzaton. Wth carefl trackng of the formaton of plastc hnges, the pshover drft can be explctly expressed n terms of the szng varables sng the prncple of vrtal work. The optmzaton methodology for the solton of the nonlnear sesmc drft desgn of bldngs s fndamentally based on an Optmalty Crtera (OC) approach. A tenstory, two-bay planar frame bldng s then presented to llstrate the detals of the OC optmzaton method for nelastc sesmc drft performance-based desgn. OPTIMAL INELASTIC SEISMIC DESIGN PROBLEM Implct Desgn Optmzaton Problem In sesmc desgn, t s commonly assmed that a bldng behaves lnear-elastcally nder mnor earthqakes and may respond nonlnear-nelastcally when sbjected to moderate and severe earthqakes. nder sch an assmpton, the entre desgn optmzaton process can therefore be decomposed nto two phases [8, 9]. In the frst phase, the strctral concrete cost s mnmzed sbject to elastc drft responses nder mnor earthqake loadng sng elastc response spectrm analyss. In ths phase, concrete member szes are consdered as desgn varables snce the concrete materal plays a more domnant role n mprovng the elastc drft performance of the bldng. Once the optmal strctral member szes are determned at the end of the frst phase of the optmzaton, the steel renforcement qanttes can then be consdered as desgn varables n the second phase. In controllng the nelastc drft responses, steel renforcement s the only effectve materal that provdes dctlty to an RC bldng strctre beyond frst yeldng. In ths second desgn phase, the member szes are kept nchanged and the cost of the steel

3 renforcement s mnmzed sbject to desgn constrants on nelastc nterstory drft prodced by the nonlnear pshover analyss. The emphass of ths paper s on the second phase of the desgn optmzaton, the nelastc sesmc drft desgn optmzaton. The detals of the frst phase elastc sesmc drft desgn optmzaton can be fond n the work of Zo [9]. For an RC bldng havng =, 2,, N members and 2N plastc hnges (assmng one hnge at each end of a member), the tenson steel renforcement rato, ρ, and the compresson steel renforcement rato, ρ, for a rectanglar cross secton are taken as desgn varables n the desgn optmzaton, whereas the member szes, B (wdth) and D (depth), are fxed. If the topology of a bldng s strctral system s predefned, the steel renforcement cost of the RC framework s mnmzed as N Mnmze: steel cost = ws ( Ls ρ + Ls ρ ) () = where w s s the cost coeffcent for steel renforcements; and L s and L s are respectvely the lengths of the tenson and compresson steel renforcements for member. Only the longtdnal flexral renforcement of member sectons s consdered as desgn varables n ths stdy, whle the transverse shear renforcement s consdered nvarant nder the assmpton that adeqate shear capacty strength s provded for each member. In the performance-based desgn, t s necessary to check the capacty of a strctre aganst the demand of an earthqake at the performance pont whch s the ntersecton of the pshover capacty and demand spectrm crves. In ths stdy, the nterstory drft responses of a bldng, generated by a specfed earthqake demand, are checked aganst approprate lmts correspondng to a gven performance level. Namely, for a mltstory bldng strctre, the nterstory drft rato cased by pshover loadng shold comply wth the followng reqrement: j j j = d j (2) h j h j where j s the nterstory drft of the jth story; j and j are the respectve story dsplacement of two adjacent j and j- floor levels; h j s the jth story heght; and d j s the specfed nterstory drft rato lmt representng the damage threshold for the jth story. Besdes the consderatons of the nterstory drft responses, local element responses sch as the sectonal plastc rotaton and member strength at the performance pont mst also be checked aganst certan acceptablty lmts. The plastc rotaton, θ, at the hth end secton of a member (where sbscrpt ph h represents one end of a member and h=, 2) shold be checked as θ ph θ p (3) where θ s the plastc rotaton lmt correspondng to a specfc performance level. Once the desgner p determnes the performance levels of the strctre (e.g., Immedate Occpancy, Lfe Safety, Collapse Preventon), the lmtng vales of θ for all members s then determned. nless specfc desgn crtera p are stated, otherwse FEMA-273 [3] and ATC-40 [2] provde gdelnes for estmatng the lmtng vales of plastc rotaton of a flexral member for varos levels of performance crtera of an RC frame. In practcal mltstory bldng strctres, excessve nmber of desgn constrants may create enormos ncreases n comptatonal effort. In order to redce the practcal bldng desgn problem to a manageable sze, the strength desgn of each member s not consdered explctly as a desgn constrant; rather, the strength-based steel renforcement ratos n accordance wth code specfcatons are frst calclated and these vales are then taken as the lower sze bond for each member n the nelastc sesmc drft desgn optmzaton.

4 In addton to the desgn performance constrants on the system-level story drft and element-level sectonal plastc rotaton, the steel renforcement varables are mposed wthn the mnmm and maxmm steel renforcement ratos as L L ρ ρ ρ ; ρ ρ ρ (4a, b) where the sperscrpts L and denote the mnmm and maxmm lmts of the desgn varables, ρ and ρ. In order to facltate a nmercal solton of the drft desgn problem, t s necessary that the mplct story drft constrant (Eq. 2) and the plastc rotaton constrant (Eq. 3) be expressed explctly n terms of the desgn varables, ρ and ρ. Explct Drft Formlaton Based on the nternal element forces and moments of the strctre obtaned from the pshover analyss at the performance pont, the prncple of vrtal work can be employed to express the pshover dsplacement. The pshover story dsplacement, j, at the performance pont ncldes the vrtal work, j, memb, prodced by the strctral members and the vrtal work, j, hnge, generated by the plastc hnges. That s, j = j, memb + j, hnge (5) n whch N C 0j C j C2 j = j, memb( B, D ) = B D B D B D (6) N 2 0 j, hnge = mpjhθ ph (7) = h= In Eq. (6), the dsplacement, j, memb, s expressed n terms of wdth ( B ) and depth ( D ) [7, 9]. Drng the nelastc drft desgn optmzaton process, j, memb s kept nchanged snce B and D of each member secton are fxed. The emphass here s on the dsplacement,, cased by the formaton of j, hnge 0 the plastc hnges. In Eq. (7), m pjh s the vrtal end moment at the locaton of the hth hnge of a member; θ ph s the actal plastc rotaton experenced by the hth plastc hnge, whch s eqal to zero when no plastc hnge s fond. As shown n Fg., the behavor of a plastc hnge s modeled as a blnear crve: the elastc segment, AB, and the hardenng segment, BC, where Pont A corresponds to the nloaded condton, Pont B s the frst yeld moment pont, Pont C s the ltmate moment capacty, whch generally corresponds to the strctral stablty performance level n ATC-40 [2]. Based on the lne segments A-B-C, the plastc rotaton, θ, can be gven as follows p M M y θ p = θ p θ p (8) M M y where θ p s the ltmate plastc rotaton whch can be establshed based on expermental tests or can be obtaned drectly from desgn gdelnes sch as the ATC-40 [2]; M s the appled moment at the locaton of the plastc hnge; M y s the bendng moment at the frst yeldng of the tensle steel; and M s the ltmate moment resstance. Gven the qantty of the steel renforcement sed n a concrete secton, the vales of M and M can then be determned. For smplcty, M can be approxmately related to M y y

5 as M =.M y [2]. M M M M y B C A θ y θ p θ p Fg.. Moment-rotaton crve For the explct problem formlaton, t s necessary that the plastc rotaton, θ p, be accrately expressed n terms of the desgn varables (.e., ρ and ρ ' ). Frthermore, a good formlaton shold reflect accrately the change n the plastc rotaton, θ, de to a change n the desgn varables drng the p optmzaton reszng process. In other words, any change n the desgn varables, ρ and ρ, drng the nelastc optmzaton process reqres a correspondng pdate on the vales of M and M y. θ D d ρ' ρ B d' d' d kd f y d' f c f s Netral axs f s Bdρ' 0.5 f c Bkd f y Bdρ M y Secton Stress Resltant force Fg. 2. Doble renforced member secton at frst yeld In pshover analyss, moment hnges are assmed and are generally assgned to the two ends of each beam or colmn. By the force eqlbrm shown n Fg. 2, where f c s the stress at the extreme compresson concrete fber, f s s the stress n the compresson steel, f y s the yeld strength of the tenson steel, and d s the effectve depth, whch s eqal to the dstance from the extreme compresson fber to the centrod of the tenson steel, M y for a moment hnge can be expressed n terms of desgn varables, ρ and ρ, as kd M y = 0.5 f c Bkd( d') + f y Bd( d d')ρ (9) 3 where k s the netral axs depth factor at the frst yeld and t s gven as 2 2 d' k = ( ρ + ρ ) nsc + 2( ρ + ρ ) nsc ( ρ + ρ ) nsc (0) d

6 E n whch n s sc = where E c and E s are the modl of elastcty of the concrete and of the steel, Ec respectvely. To take nto accont the change n θ p de to a change n ρ and ρ whle mantanng an nstantaneosly fxed vale of M, a second-order Taylor seres approxmaton for evalatng the vale of θ s gven as p 2 θ θ θ ( ρ) = p p p θ p ( ρ ρ ) + 0 ( ρ ρ ) () ρ= ρ ρ= ρ 2 ρ ρ ρ 2 ρ = where the tenson steel rato, ρ, s consdered as the major desgn varable; for smplcty, the compresson steel rato, ρ, s assmed to be lnearly related to ρ for beams and to be the same as ρ for colmns. Gven the explct expresson of M y as a fncton of ρ from Eq. (9), the gradent, 2 θ p θ, and the second-order term, p, can be analytcally calclated from Eq. (8). 2 ρ ρ By sbstttng the explct plastc rotaton, (ρ), gven n Eq. () nto Eq. (7), the pshover dsplacement, θ p j, n Eq. (5) can also be explctly expressed n terms of the desgn varable, Plastc Rotaton Constrant and Szng Constrant In ths desgn optmzaton, when the plastc rotaton, θ p, s to be modfed wth changes n the desgn varable, ρ, t s necessary to make sre that θ p does not exceed the specfed threshold of plastc rotaton, θ p, for each specfed performance level. Moreover, n order to prevent drastc changes n the nternal element force and moment redstrbton de to the changes n the desgn varables resltng n flctaton of solton convergence drng the pshover reanalyss and desgn optmzaton processes, t s necessary that each plastc hnge reman plastc once t appears drng the reszng teraton of the desgn varables. Frthermore, to mantan the accracy of the Taylor approxmaton of the pshover dsplacement n Eq. (), t s necessary to ensre that the varaton of ρ for the members wth plastc hnges be restrcted wthn a relatvely small range. It s fond from Fg. that, n order to mantan the relatonshp of 0 θ, the nternal moment, p θ p M, leadng to the occrrence of a plastc hnge mst satsfy the followng condton: M M (2) y M As a reslt, based on Eq. (2), the lower and pper bonds of ρ for each plastc hnge can be nstantaneosly establshed drng the OC teratve reszng process. It shold be noted that the proper establshment of the lower and pper bonds of ρ not only lmts the changes n the steel renforcement desgn varables, bt also t satsfes the local performance-based constrants on the control of the local plastc rotaton at the ends of members. Explct Desgn Problem Formlaton pon establshng the explct nelastc drft formlaton, the optmzaton problem of mnmzng the steel constrcton cost of a mltstory RC bldng can be explctly wrtten n terms of the desgn varable, ρ, as Mnmze: N = = ρ. F( ρ ) w s ρ (3)

7 sbject to: where N N g j ( ρ ) = j 0 + α + d j = ( ρ ρ ) α 2 ( ρ ρ ) ρ ρ h j = 2 = j =, 2,..., N ) (4) ( j ρ ρ ρ =, 2,..., N ) (5) L = j = 2 0 ρ = ρ ρ h= ( θ 0 ph α m jh 0 (6a) ρ = ρ ρ 2 j ph = = 2 2 θ 0 α 2 0 m jh 0 (6b) 2 ρ = ρ 2 ρ = ρ ρ h= ρ In Eq. (3), w s s the cost coeffcent for the steel renforcement, ρ. Eq. (4) defnes the set of sesmc nterstory drft performance constrants nder specfed earthqake grond motons. Eq. (5) defnes the L szng constrants for the steel renforcement, where ρ and ρ correspond to the lower and pper sze bonds specfed for the tensle steel renforcement varable, ρ, and they shold be pdated after each nonlnear pshover analyss. Once the desgn optmzaton problem s explctly expressed n terms of desgn varables, the next task s to adopt a stable method for solvng the problem. The OC approach s proposed to solve the explct nelastc drft optmzaton problem Eqs. (3)-(5). When sng the OC technqe, a set of necessary optmalty condtons for the desgn are derved and then a recrsve algorthm s appled to resze the strctre to satsfy the optmalty condtons and ths ndrectly optmze the strctre. Frther detals of the OC technqe can be fond n the reference by Zo [9]. Intal Preprocessor Althogh the OC method does not mpose any restrctons on the ntal vales of the desgn varables, the rate of convergence of the OC process depends on the ntal desgn vales. To speed p the convergence of the OC process, t s mportant to begn wth a reasonably good startng desgn. One effectve scalng approach was proposed by Chan [0]. In sch an approach, the desgn optmzaton wth a sngle drft constrant s frst consdered and a smple closed form solton for the problem s derved analytcally and exploted as an ntal preprocessor for the teratve OC process. The advantage of ths approach s that a reasonable ntal desgn can be qckly establshed based on a representatve sngle drft constrant. Experence ndcates that the ntal preprocessor can generally lead to steady and rapd solton convergence of the mltple nelastc drft optmzaton problem. OVERALL DESIGN OPTIMIZATION PROCEDRE The overall desgn optmzaton procedre for lmtng lateral elastc and nelastc drfts of a renforced concrete bldng strctre s lsted as follows. () Establsh an ntal desgn wth optmal member szes, whch can be obtaned from the elastc sesmc desgn optmzaton by mnmzng the concrete cost of an RC strctre sbjected to a mnor earthqake loadng sng the elastc response spectrm analyss method [9]. (2) Determne the desgn spectrm, correspondng to a severe earthqake event, that wll be sed n the nonlnear pshover analyss. (3) Condct a statc vrtal load analyss to obtan the member nternal forces that wll be sed n formlatng nelastc drft responses by employng the prncple of vrtal work.

8 (4) Based on the optmal member sze, determne the mnmm and maxmm sze bonds of the steel renforcement ratos, ρ and ρ, n accordance wth the strength-based code reqrements. (5) Apply the ntal preprocessor based on a representatve sngle drft constrant to establsh a reasonable startng set of steel renforcement desgn varables for the mltple drft constraned optmzaton. (6) Carry ot the nonlnear pshover analyss sng commercally avalable software to determne the performance pont of the strctre and the assocated nelastc drft responses of the strctre at the performance pont. (7) Track down the locatons of the plastc hnges, establsh the nstantaneos lower and pper bond move lmts of ρ for those members wth plastc hnges based on Eq. (2) and determne the vales of the frst-order and second-order dervatves of the drft responses sng Eqs. (6a) and (6b). (8) Establsh the explct nterstory drft constrants sng a second-order Taylor seres approxmaton and formlate the explct desgn problem, Eqs. (3)-(5). (9) Apply the recrsve OC optmzaton algorthm to resze all steel renforcement desgn varables and to dentfy the actve nelastc drft constrants. (0) Check convergence of the steel cost and the nelastc drft performance of the strctre. Termnate wth the optmm desgn f the solton convergence s fond; otherwse, retrn to Step 6. ILLSTRATIVE EXAMPLE A ten-story, two-bay planar frame s sed to llstrate the proposed optmal desgn method. The geometry of the example s gven n Fg. 3. Concrete wth the cylnder strength of 20MPa and steel renforcement wth the yeld strength, f y, of 335MPa are sed for all members. The loads consdered n the pshover analyss are lateral sesmc loads and vertcal gravty loads. Whle the lateral loads are ncrementally appled, the gravty loads are mantaned to be nchanged drng the nonlnear pshover analyss process. A nformly dstrbted gravty load of 30kN/m s to be appled to the beams of each story. Intal member szes of the framework sed to commence the nelastc desgn optmzaton are shown n Table 3. Intal steel renforcement ratos are frst calclated based on the strength reqrements of the RC members n accordance wth the Chnese sesmc desgn code []. Sch strength-based renforcement ratos are taken ntally as the lower bonds for the nelastc desgn optmzaton. The pper sze bonds of the steel renforcement ratos are assmed to be 6.0% for colmns and 4.0% for beams. For smplcty, symmetrcal arrangement of the steel renforcement of each member s assmed sch that ρ = ρ. Flexral moment hnges are assgned to the end locatons of the beams and colmns and the ltmate plastc hnge rotaton, θ p, s assmed to be 0.02 radan The 5% damped desgn spectrm wth an ntal peak acceleraton of.4g accordng to the Chnese sesmc desgn code [] s modfed by the spectral redcton method n the pshover analyss of ths example. A typcal nt constrcton cost of the steel renforcement (ncldng the costs of the steel materal and the labor) of S$950/tonne s assmed. Interstory drft constrants are consdered wth an assmed allowable nterstory drft rato lmt of /00. The ntal preprocessor, n whch analytcal optmzaton wth the top dsplacement constrant alone s consdered, s appled before the mltple nelastc nterstory drft constraned optmzaton s nvoked. The desgn process s deemed to converge when the dfference n the strctre costs for two sccessve desgn cycles s wthn 0.5% and when the dfference between the actve nterstory drft vale and ts allowable lmt at the performance pont s wthn 0.5%.

9 q=30kn/m 3.0*0=30m C B C2 5m B2 C3 5m Fg. 3 A Ten-Story, Two-Bay Frame The optmal desgn hstory of the example s presented n Fg. 4. The ntal preprocessng s frst completed and the mltple nelastc nterstory drft constraned optmzaton s then commenced. It s fond that there s a relatvely large ncrease of 4% n the steel cost from the ntal S$693 to S$93 after the ntal preprocessng wth only the top dsplacement constrant. However, the optmal desgn process wth the mltple drft constrant converges slowly bt steadly wthn desgn cycles, wth only a slght dfference of 2% n the steel cost from S$93 to the fnal S$978. Relatvely slow, bt steady, convergence s fond de to the need for mantanng a small change n the steel renforcement ratos drng the nelastc desgn optmzaton process Steel renforcement cost (S$) Desgn hstory Fg. 4 Desgn Hstory of Steel Renforcement Costs

10 Table Intal and Fnal Steel Renforcement Ratos Element Story Member Intal member szes Steel ratos type level grop Wdth Depth Intal Optmal (mm) (mm) (%) (%) 9th~0th C,C C th C,C C Colmn 7th C,C C th C,C C th C,C C th C,C C rd C,C C nd C,C C st C,C C th~0th B,B th B,B Beam 7th B,B th B,B th B,B th B,B rd B,B nd B,B st B,B Table presents the ntal and optmal steel renforcement ratos. Intally, the startng desgn wth strength-based steel renforcement s fond to be nfeasble n terms of the assmed allowable nterstory drft lmt. After the optmzaton, the steel renforcement ratos of the beams greatly ncrease partclarly n the lower levels of the strctre, whle those of colmns are fond wth lttle changes. Fg. 5 presents the ntal and fnal performance ponts respectvely. The performance pont P of the ntal strctre has a spectral acceleraton capacty of 0.068g and a spectral dsplacement capacty of 0.262m. The optmzed strctre correspondng to the fnal performance pont P2 has a spectral acceleraton capacty of 0.086g and a spectral dsplacement of 0.2m. Sch a reslt of shftng the spectral dsplacement from 0.262m to 0.2m ndcates that, for the optmzed nelastc frame, the nelastc lateral load resstance has been enhanced throgh optmal reszng of the steel renforcement by the OC procedre. Also, shftng the ltmate spectral acceleraton capacty from 0.068g to 0.085g ndcates that the optmzed strctre attracts an ncrease n the sesmc loadng acton and therefore, reqres the strctre to be stffened. The OC procedre developed s fond to be able to atomatcally drve from any ntal performance pont to the fnal performance pont resltng n the mnmm cost desgn.

11 Spectral Acceleraton, Sa (g) Fnal demand Intal demand Fnal capacty Intal capacty Fnal performance pont (0.2, 0.086g) P2 P Intal performance pont (0.262, 0.068g) Spectral Dsplacement, Sd (m) Fg. 5 Comparson of Performance Ponts 0 8 Drft rato lmt=/00 Intal drft rato Fnal drft rato Story level Interstory drft rato Fg. 6 Intal and Fnal Interstory Drft Ratos The ntal and fnal nterstory drft ratos are shown n Fg. 6. The ntal nterstory drft constrants at the second throgh the eghth floors are fond to volate sbstantally the allowable nterstory drft rato lmt of /00, resltng n the occrrence of the weak stores on the these floor levels of the bldng. However, these pshover nterstory drft constrants are fond to be close to and wthn the allowable vales after the optmzaton, ndcatng that a rather nform nterstory drft dstrbton over the heght of the bldng has been acheved and the occrrence of weak story has been prevented at the optmm performance pont.

12 Nmber of plastc hnges Total nmber B-IO IO-LS LS-CP of plastc hnges Intal desgn Fnal desgn (a) Nmber of Plastc Hnges A M B IO C CP LS B-IO: IO-LS: LS-CP: θ (b) Intal Desgn (c) Fnal Desgn Fg. 7 Intal and Fnal Plastc Hnge Dstrbton Fg. 7(a) ncldes a table showng the nmber of plastc hnges at three dfferent performance states. Fgs. 7(b)-7(c) show the ntal and fnal plastc hnge dstrbtons nder the pshover loadng at the performance pont of the strctre. No plastc hnge rotaton s fond to exceed the specfed threshold of plastc rotaton. As shown n Fg. 7(b), the rotatons of twenty plastc hnges of the ntal desgn are fond to be located between the LS-CP state. However, after the optmzaton, most of the plastc hnges are fond to be n the B-IO and IO-LS states and only one hnge s n the LS-CP state, as can be observed from the optmzed framework n Fg. 7(c). Frthermore, the nterstory drfts along the heght of the bldng are also fond to be almost all flly constraned at the optmm, resltng n a rather lnear deflected profle of the nelastc desgn. Sch a reslt frther ndcates that the optmzaton method developed can atomatcally resze the steel renforcements of all members to attan a nform dctlty demand along the heght of the mltstory bldng.

13 CONCLSIONS It has been demonstrated that steel renforcement plays a sgnfcant role n controllng the lateral drft beyond frst yeldng and n provdng dctlty to an RC bldng framework. sng the prncple of vrtal work and the Taylor seres approxmaton, the nelastc performance-based sesmc desgn problem has been explctly expressed n terms of the steel renforcement desgn varables. It s demonstrated that the OC desgn method s able to mprove atomatcally and gradally a performance-based nterstory drft desgn to attan optmal performance. Also, ths OC desgn method developed s able to atomatcally shft any ntal performance pont to acheve the fnal optmal performance pont. However, the restrctve move lmt mposed on the steel renforcement desgn varables s necessary to ensre a smooth and steady convergence of the nelastc drft desgn process. At optmm, a nform lateral drft or dctlty demand over all stores of the bldng wth the mnmm cost s acheved. It s also beleved that ths optmzaton methodology provdes a powerfl compter-based technqe for performance-based desgn of mltstory RC bldng strctres. The proposed optmzaton methodology provdes a good bass for more comprehensve performance-based optmzaton of strctres as more accrate nonlnear pshover procedres takng nto the hgher mode effects are developed and mltple levels of performance crtera and desgn objectves are to be smltaneosly consdered. ACKNOWLEDGMENTS The athors are gratefl for the fnancal spport provded by the Research Grants Concl of Hong Kong nder Project No. HKST6249/00E. REFERENCES. SEAOC Vson 2000 Commttee. Performance Based Sesmc Engneerng of Bldngs, Part 2: Conceptal Framework. Sacramento, Calforna: Strctral Engneers Assocaton of Calforna, ATC-40. Sesmc Evalaton and Retroft of Concrete Bldngs. Volme, ATC-40 Report. Redwood Cty, Calforna: Appled Technology Concl, FEMA. NEHRP Gdelnes for the Sesmc Rehabltaton of Bldngs. Washngton, D.C., SA: Bldng Sesmc Safety Concl for the Federal Emergency Management Agency (Report No. FEMA 273), Moehle JP, Mahn SA. Observatons on the Behavor of Renforced Concrete Bldngs Drng Earthqakes. Ghosh SK, Edtor. Earthqake-Resstant Concrete Strctres Inelastc Response and Desgn. Amercan Concrete Insttte SP-27, Ghobarah A, Aly NM, El-Attar M. Performance Level Crtera and Evalaton. Fajfar P, Krawnkler H, Edtors. Proceedngs of the nternatonal workshop on Sesmc Desgn Methodologes for the Next Generaton of Codes, Slovena, Chopra AK. Dynamcs of Strctres: theory and applcatons to earthqake engneerng. Englewood Clffs, New Jersey, Chan CM. Optmal Lateral Stffness Desgn of Tall Bldngs of Mxed Steel and Concrete Constrcton. Jornal of Strctral Desgn of Tall Bldngs 200; 0: Zo XK, Chan CM. Optmal Drft Performance Desgn for Nonlnear Pshover Response of Concrete Strctres. WCSMO-4: Proceedngs of the Forth Wold Congress of Strctral and Mltdscplnary Optmzaton. Jne 4-8, Dalan, Chna, Zo XK. Optmal Sesmc Performance-Based Desgn of Renforced Concrete Bldngs. Ph.D. Dssertaton, Hong Kong nversty of Scence and Technology, 2002.

14 0. Chan CM. How to Optmze Tall Steel Bldng Frameworks. Gde to Strctral Optmzaton, ASCE Manals and Reports on Engneerng Practce No.90, Amercan Socety of Cvl Engneers, 997; Natonal Standard of the People s Repblc of Chna. Chnese Code for Sesmc Desgn Bldngs (GBJ-89). Bejng, Chna: New World Press, 994.

Research Note NONLINEAR ANALYSIS OF SEMI-RIGID FRAMES WITH RIGID END SECTIONS * S. SEREN AKAVCI **

Research Note NONLINEAR ANALYSIS OF SEMI-RIGID FRAMES WITH RIGID END SECTIONS * S. SEREN AKAVCI ** Iranan Jornal of Scence & Technology, Transacton, Engneerng, Vol., No., pp 7-7 rnted n The Islamc Repblc of Iran, 7 Shraz Unversty Research Note NONLINER NLYSIS OF SEMI-RIGID FRMES WIT RIGID END SECTIONS

More information

AE/ME 339. K. M. Isaac. 8/31/2004 topic4: Implicit method, Stability, ADI method. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.

AE/ME 339. K. M. Isaac. 8/31/2004 topic4: Implicit method, Stability, ADI method. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept. AE/ME 339 Comptatonal Fld Dynamcs (CFD) Comptatonal Fld Dynamcs (AE/ME 339) Implct form of dfference eqaton In the prevos explct method, the solton at tme level n,,n, depended only on the known vales of,

More information

MEMBRANE ELEMENT WITH NORMAL ROTATIONS

MEMBRANE ELEMENT WITH NORMAL ROTATIONS 9. MEMBRANE ELEMENT WITH NORMAL ROTATIONS Rotatons Mst Be Compatble Between Beam, Membrane and Shell Elements 9. INTRODUCTION { XE "Membrane Element" }The comple natre of most bldngs and other cvl engneerng

More information

Second Order Analysis

Second Order Analysis Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

NON LINEAR ANALYSIS OF STRUCTURES ACCORDING TO NEW EUROPEAN DESIGN CODE

NON LINEAR ANALYSIS OF STRUCTURES ACCORDING TO NEW EUROPEAN DESIGN CODE October 1-17, 008, Bejng, Chna NON LINEAR ANALYSIS OF SRUCURES ACCORDING O NEW EUROPEAN DESIGN CODE D. Mestrovc 1, D. Czmar and M. Pende 3 1 Professor, Dept. of Structural Engneerng, Faculty of Cvl Engneerng,

More information

Structural Dynamics and Earthquake Engineering

Structural Dynamics and Earthquake Engineering Structural Dynamcs and Earthuake Engneerng Course 9 Sesmc-resstant desgn of structures (1) Sesmc acton Methods of elastc analyss Course notes are avalable for download at http://www.ct.upt.ro/users/aurelstratan/

More information

SE Story Shear Frame. Final Project. 2 Story Bending Beam. m 2. u 2. m 1. u 1. m 3. u 3 L 3. Given: L 1 L 2. EI ω 1 ω 2 Solve for m 2.

SE Story Shear Frame. Final Project. 2 Story Bending Beam. m 2. u 2. m 1. u 1. m 3. u 3 L 3. Given: L 1 L 2. EI ω 1 ω 2 Solve for m 2. SE 8 Fnal Project Story Sear Frame Gven: EI ω ω Solve for Story Bendng Beam Gven: EI ω ω 3 Story Sear Frame Gven: L 3 EI ω ω ω 3 3 m 3 L 3 Solve for Solve for m 3 3 4 3 Story Bendng Beam Part : Determnng

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

Connection and Tension Member Design

Connection and Tension Member Design ACH 631 Note Set 17.1 F2015abn Connecton and Tenson Member Desgn Notaton: A = area (net = wth holes, bearng = n contact, etc...) A e = effectve net area fond from the prodct of the net area A n by the

More information

DUE: WEDS FEB 21ST 2018

DUE: WEDS FEB 21ST 2018 HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant

More information

Optimum Design of Steel Frames Considering Uncertainty of Parameters

Optimum Design of Steel Frames Considering Uncertainty of Parameters 9 th World Congress on Structural and Multdscplnary Optmzaton June 13-17, 211, Shzuoka, Japan Optmum Desgn of Steel Frames Consderng ncertanty of Parameters Masahko Katsura 1, Makoto Ohsak 2 1 Hroshma

More information

BAR & TRUSS FINITE ELEMENT. Direct Stiffness Method

BAR & TRUSS FINITE ELEMENT. Direct Stiffness Method BAR & TRUSS FINITE ELEMENT Drect Stness Method FINITE ELEMENT ANALYSIS AND APPLICATIONS INTRODUCTION TO FINITE ELEMENT METHOD What s the nte element method (FEM)? A technqe or obtanng approxmate soltons

More information

CLOSED-FORM CHARACTERIZATION OF THE CHANNEL CAPACITY OF MULTI-BRANCH MAXIMAL RATIO COMBINING OVER CORRELATED NAKAGAMI FADING CHANNELS

CLOSED-FORM CHARACTERIZATION OF THE CHANNEL CAPACITY OF MULTI-BRANCH MAXIMAL RATIO COMBINING OVER CORRELATED NAKAGAMI FADING CHANNELS CLOSED-FORM CHARACTERIZATION OF THE CHANNEL CAPACITY OF MULTI-BRANCH MAXIMAL RATIO COMBINING OVER CORRELATED NAKAGAMI FADING CHANNELS Yawgeng A. Cha and Karl Yng-Ta Hang Department of Commncaton Engneerng,

More information

Finite Element Modelling of truss/cable structures

Finite Element Modelling of truss/cable structures Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures

More information

CS 3750 Machine Learning Lecture 6. Monte Carlo methods. CS 3750 Advanced Machine Learning. Markov chain Monte Carlo

CS 3750 Machine Learning Lecture 6. Monte Carlo methods. CS 3750 Advanced Machine Learning. Markov chain Monte Carlo CS 3750 Machne Learnng Lectre 6 Monte Carlo methods Mlos Haskrecht mlos@cs.ptt.ed 5329 Sennott Sqare Markov chan Monte Carlo Importance samplng: samples are generated accordng to Q and every sample from

More information

FUZZY FINITE ELEMENT METHOD

FUZZY FINITE ELEMENT METHOD FUZZY FINITE ELEMENT METHOD RELIABILITY TRUCTURE ANALYI UING PROBABILITY 3.. Maxmum Normal tress Internal force s the shear force, V has a magntude equal to the load P and bendng moment, M. Bendng moments

More information

Indeterminate pin-jointed frames (trusses)

Indeterminate pin-jointed frames (trusses) Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all

More information

A New Method for Marshaling of Freight Train with Generalization Capability Based on the Processing Time

A New Method for Marshaling of Freight Train with Generalization Capability Based on the Processing Time Proceedngs of the Internatonal MltConference of Engneers and Compter Scentsts 205 Vol I, A New Method for Marshalng of Freght Tran wth Generalzaton Capablty Based on the Processng Tme Yoch Hrashma Abstract

More information

GEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE

GEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE GEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE Prof. J. N. Mandal Department of cvl engneerng, IIT Bombay, Powa, Mumba 400076, Inda. Tel.022-25767328 emal: cejnm@cvl.tb.ac.n Module - 9 LECTURE - 48

More information

Composite model 131: Ply types 1 and 7 calibration

Composite model 131: Ply types 1 and 7 calibration Composte model 3: Ply types and 7 calbraton Model calbraton Composte Global Ply model 3 for elastc, damage and falre PAM-CRASH materal model 3 s for mlt-layered composte shell elements. Wthn ths model

More information

9.2 Seismic Loads Using ASCE Standard 7-93

9.2 Seismic Loads Using ASCE Standard 7-93 CHAPER 9: Wnd and Sesmc Loads on Buldngs 9.2 Sesmc Loads Usng ASCE Standard 7-93 Descrpton A major porton of the Unted States s beleved to be subject to sesmc actvty suffcent to cause sgnfcant structural

More information

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,

More information

Torsion Stiffness of Thin-walled Steel Beams with Web Holes

Torsion Stiffness of Thin-walled Steel Beams with Web Holes Torson Stffness of Thn-walled Steel Beams wth Web Holes MARTN HORÁČEK, JNDŘCH MELCHER Department of Metal and Tmber Structures Brno Unversty of Technology, Faculty of Cvl Engneerng Veveří 331/95, 62 Brno

More information

Suppression of Low-frequency Lateral Vibration in Tilting Vehicle Controlled by Pneumatic Power

Suppression of Low-frequency Lateral Vibration in Tilting Vehicle Controlled by Pneumatic Power Challenge D: A world of servces for passengers Sppresson of Low-freqency Lateral Vbraton n ltng Vehcle Controlled by Pnematc Power A. Kazato, S.Kamoshta Ralway echncal Research Insttte, okyo, Japan. Introdcton

More information

Proceedings of the 11th WSEAS International Conference on Sustainability in Science Engineering

Proceedings of the 11th WSEAS International Conference on Sustainability in Science Engineering 4,50m 4,00m 4,00m 4,00m 4,00m 4,00m 4,00m 4,00m 4,00m 4,00m 40,50m 7,50m 7,50m 7,50m Proceedngs of the 11th WSEAS Internatonal Conference on Sustanablty n Scence Engneerng OERSTREGTH DEMADS I MULTISTOREY

More information

Introduction to elastic wave equation. Salam Alnabulsi University of Calgary Department of Mathematics and Statistics October 15,2012

Introduction to elastic wave equation. Salam Alnabulsi University of Calgary Department of Mathematics and Statistics October 15,2012 Introdcton to elastc wave eqaton Salam Alnabls Unversty of Calgary Department of Mathematcs and Statstcs October 15,01 Otlne Motvaton Elastc wave eqaton Eqaton of moton, Defntons and The lnear Stress-

More information

Preliminary Design of Moment-Resisting Frames

Preliminary Design of Moment-Resisting Frames Prelmnary Desgn of Moment-Resstng Frames Preprnt Aamer Haque Abstract A smple method s developed for prelmnary desgn of moment-resstng frames. Preprnt submtted to Elsever August 27, 2017 1. Introducton

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion Assessment of Ste Amplfcaton Effect from Input Energy Spectra of Strong Ground Moton M.S. Gong & L.L Xe Key Laboratory of Earthquake Engneerng and Engneerng Vbraton,Insttute of Engneerng Mechancs, CEA,

More information

COMBINED HIGH AND LOW CYCLE FATIGUE MODEL FOR PREDICTION OF STEEL BRIDGE LIVES

COMBINED HIGH AND LOW CYCLE FATIGUE MODEL FOR PREDICTION OF STEEL BRIDGE LIVES 505 COMBINED HIGH AND LOW CYCLE FATIGUE MODEL FOR PREDICTION OF STEEL BRIDGE LIVES M. OHGA *, K. KARUNANANDA, P.B.R. DISSANAYAKE 2 and S.A.S.C. SIRIWARDANE 3 Department of Cvl and Envronmental Engneerng,

More information

Reading Assignment. Panel Data Cross-Sectional Time-Series Data. Chapter 16. Kennedy: Chapter 18. AREC-ECON 535 Lec H 1

Reading Assignment. Panel Data Cross-Sectional Time-Series Data. Chapter 16. Kennedy: Chapter 18. AREC-ECON 535 Lec H 1 Readng Assgnment Panel Data Cross-Sectonal me-seres Data Chapter 6 Kennedy: Chapter 8 AREC-ECO 535 Lec H Generally, a mxtre of cross-sectonal and tme seres data y t = β + β x t + β x t + + β k x kt + e

More information

Lecture Notes on Linear Regression

Lecture Notes on Linear Regression Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume

More information

Experimental Study on Ultimate Strength of Flexural-Failure-Type RC Beams under Impact Loading

Experimental Study on Ultimate Strength of Flexural-Failure-Type RC Beams under Impact Loading xpermental Study on Ultmate Strength of Flexural-Falure-Type RC Beams under Impact Loadng N. Ksh 1), O. Nakano 2~, K. G. Matsuoka 1), and T. Ando 1~ 1) Dept. of Cvl ngneerng, Muroran Insttute of Technology,

More information

A Reinforcement Learning for Freight Train with Augmented Collective Motions in a Marshaling

A Reinforcement Learning for Freight Train with Augmented Collective Motions in a Marshaling Proceedngs of the Internatonal MltConference of Engneers and Compter Scentsts 207 Vol I, IMECS 207, March 5-7, 207, Hong Kong A Renforcement earnng for Freght Tran wth Agmented Collectve Motons n a Marshalng

More information

CONTINUUM SOLUTION OF LATERAL LOADING OF LARGE PILE GROUPS CUED/D-SOILS/TR 334 (July 2004) By A. Klar, A.D. Spasojevic and K. Soga

CONTINUUM SOLUTION OF LATERAL LOADING OF LARGE PILE GROUPS CUED/D-SOILS/TR 334 (July 2004) By A. Klar, A.D. Spasojevic and K. Soga CONTINUUM SOLUTION OF LATERAL LOADIN OF LARE PILE ROUPS CUED/D-SOILS/TR 334 (Jl 4) B A. Klar, A.D. Spasoevc and K. Soga Klar, A. Spasoevc, A.D. and Soga, K. Contnm Solton of Lateral Loadng of Large Ple

More information

The Impact of Instrument Transformer Accuracy Class on the Accuracy of Hybrid State Estimation

The Impact of Instrument Transformer Accuracy Class on the Accuracy of Hybrid State Estimation 1 anel Sesson: Addressng Uncertanty, Data alty and Accracy n State Estmaton The mpact of nstrment Transformer Accracy Class on the Accracy of Hybrd State Estmaton Elas Kyrakdes and Markos Aspro KOS Research

More information

DESIGN OF STEEL PLATE SHEAR WALLS CONSIDERING BOUNDARY FRAME MOMENT RESISTING ACTION. B. Qu 1 and M.Bruneau 2 ABSTRACT

DESIGN OF STEEL PLATE SHEAR WALLS CONSIDERING BOUNDARY FRAME MOMENT RESISTING ACTION. B. Qu 1 and M.Bruneau 2 ABSTRACT Proceedngs of the 9th U.S. Natonal and 0th Canadan Conference on Earthquake Engneerng Compte Rendu de la 9ème Conférence Natonale Amércane et 0ème Conférence Canadenne de Géne Parassmque July 5-9, 00,

More information

KULLBACK-LEIBER DIVERGENCE MEASURE IN CORRELATION OF GRAY-SCALE OBJECTS

KULLBACK-LEIBER DIVERGENCE MEASURE IN CORRELATION OF GRAY-SCALE OBJECTS The Second Internatonal Conference on Innovatons n Informaton Technology (IIT 05) KULLBACK-LEIBER DIVERGENCE MEASURE IN CORRELATION OF GRAY-SCALE OBJECTS M. Sohal Khald Natonal Unversty of Scences and

More information

A Mechanics-Based Approach for Determining Deflections of Stacked Multi-Storey Wood-Based Shear Walls

A Mechanics-Based Approach for Determining Deflections of Stacked Multi-Storey Wood-Based Shear Walls A Mechancs-Based Approach for Determnng Deflectons of Stacked Mult-Storey Wood-Based Shear Walls FPINNOVATIONS Acknowledgements Ths publcaton was developed by FPInnovatons and the Canadan Wood Councl based

More information

CHAPTER 9 CONCLUSIONS

CHAPTER 9 CONCLUSIONS 78 CHAPTER 9 CONCLUSIONS uctlty and structural ntegrty are essentally requred for structures subjected to suddenly appled dynamc loads such as shock loads. Renforced Concrete (RC), the most wdely used

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

Fracture analysis of FRP composites using a meshless finite point collocation method

Fracture analysis of FRP composites using a meshless finite point collocation method Forth Internatonal Conference on FRP Compostes n Cvl Engneerng (CICE008) -4Jly 008, Zrch, Swtzerland Fractre analyss of FRP compostes sng a meshless fnte pont collocaton method M. Shahverd & S. Mohammad

More information

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project SE 8 Fnal Project Story Shear Frame u m Gven: u m L L m L L EI ω ω Solve for m Story Bendng Beam u u m L m L Gven: m L L EI ω ω Solve for m 3 3 Story Shear Frame u 3 m 3 Gven: L 3 m m L L L 3 EI ω ω ω

More information

Some modelling aspects for the Matlab implementation of MMA

Some modelling aspects for the Matlab implementation of MMA Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton

More information

Report on Image warping

Report on Image warping Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.

More information

Modeling Mood Variation and Covariation among Adolescent Smokers: Application of a Bivariate Location-Scale Mixed-Effects Model

Modeling Mood Variation and Covariation among Adolescent Smokers: Application of a Bivariate Location-Scale Mixed-Effects Model Modelng Mood Varaton and Covaraton among Adolescent Smokers: Applcaton of a Bvarate Locaton-Scale Mxed-Effects Model Oksana Pgach, PhD, Donald Hedeker, PhD, Robn Mermelsten, PhD Insttte for Health Research

More information

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS

More information

Chapter 3. Estimation of Earthquake Load Effects

Chapter 3. Estimation of Earthquake Load Effects Chapter 3. Estmaton of Earthquake Load Effects 3.1 Introducton Sesmc acton on chmneys forms an addtonal source of natural loads on the chmney. Sesmc acton or the earthquake s a short and strong upheaval

More information

DESIGN OPTIMIZATION OF CFRP RECTANGULAR BOX SUBJECTED TO ARBITRARY LOADINGS

DESIGN OPTIMIZATION OF CFRP RECTANGULAR BOX SUBJECTED TO ARBITRARY LOADINGS Munch, Germany, 26-30 th June 2016 1 DESIGN OPTIMIZATION OF CFRP RECTANGULAR BOX SUBJECTED TO ARBITRARY LOADINGS Q.T. Guo 1*, Z.Y. L 1, T. Ohor 1 and J. Takahash 1 1 Department of Systems Innovaton, School

More information

EEL 6266 Power System Operation and Control. Chapter 3 Economic Dispatch Using Dynamic Programming

EEL 6266 Power System Operation and Control. Chapter 3 Economic Dispatch Using Dynamic Programming EEL 6266 Power System Operaton and Control Chapter 3 Economc Dspatch Usng Dynamc Programmng Pecewse Lnear Cost Functons Common practce many utltes prefer to represent ther generator cost functons as sngle-

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

An Interactive Optimisation Tool for Allocation Problems

An Interactive Optimisation Tool for Allocation Problems An Interactve Optmsaton ool for Allocaton Problems Fredr Bonäs, Joam Westerlund and apo Westerlund Process Desgn Laboratory, Faculty of echnology, Åbo Aadem Unversty, uru 20500, Fnland hs paper presents

More information

Nonlinear Programming Formulations for On-line Applications

Nonlinear Programming Formulations for On-line Applications Nonlnear Programmng Formlatons for On-lne Applcatons L. T. Begler Carnege Mellon Unversty Janary 2007 Jont work wth Vctor Zavala and Carl Lard NLP for On-lne Process Control Nonlnear Model Predctve Control

More information

STUDY ON SEISMIC BEHAVIOR OF RC COMPOSITE CORE WALLS WITH CONCEALED STEEL TRUSS SUBJECTED TO COMBINED ACTION

STUDY ON SEISMIC BEHAVIOR OF RC COMPOSITE CORE WALLS WITH CONCEALED STEEL TRUSS SUBJECTED TO COMBINED ACTION STUDY ON SEISMIC BEHAVIOR OF RC COMPOSITE CORE WALLS WITH CONCEALED STEEL TRUSS SUBJECTED TO COMBINED ACTION CAO Wanln 1, CHANG Wehua 2, ZHANG Janwe 1 1 College of archtecture and Cvl Engneerng, Bejng

More information

Identification of Instantaneous Modal Parameters of A Nonlinear Structure Via Amplitude-Dependent ARX Model

Identification of Instantaneous Modal Parameters of A Nonlinear Structure Via Amplitude-Dependent ARX Model Identfcaton of Instantaneous Modal Parameters of A Nonlnear Structure Va Ampltude-Dependent ARX Model We Chh Su(NCHC), Chung Shann Huang(NCU), Chng Yu Lu(NCU) Outlne INRODUCION MEHODOLOGY NUMERICAL VERIFICAION

More information

Annexes. EC.1. Cycle-base move illustration. EC.2. Problem Instances

Annexes. EC.1. Cycle-base move illustration. EC.2. Problem Instances ec Annexes Ths Annex frst llustrates a cycle-based move n the dynamc-block generaton tabu search. It then dsplays the characterstcs of the nstance sets, followed by detaled results of the parametercalbraton

More information

APPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam

APPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam APPENDIX F A DISPACEMENT-BASED BEAM EEMENT WITH SHEAR DEFORMATIONS Never use a Cubc Functon Approxmaton for a Non-Prsmatc Beam F. INTRODUCTION { XE "Shearng Deformatons" }In ths appendx a unque development

More information

SIMPLIFIED APPROACH TO THE NON-LINEAR BEHAVIOR OF RC MEMBERS

SIMPLIFIED APPROACH TO THE NON-LINEAR BEHAVIOR OF RC MEMBERS SIMPLIFIED APPROACH TO THE NON-LINEAR BEHAVIOR OF RC MEMBERS Shahd NASIR 1, Supratc GUPTA 2 And Hdetaka UMEHARA 3 SUMMARY In ths paper, a smplfed one-dmensonal analytcal tool based on fnte dfference technque

More information

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009 College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:

More information

In this section is given an overview of the common elasticity models.

In this section is given an overview of the common elasticity models. Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

EMISSION MEASUREMENTS IN DUAL FUELED INTERNAL COMBUSTION ENGINE TESTS

EMISSION MEASUREMENTS IN DUAL FUELED INTERNAL COMBUSTION ENGINE TESTS XVIII IEKO WORLD NGRESS etrology for a Sstanable Development September, 17, 006, Ro de Janero, Brazl EISSION EASUREENTS IN DUAL FUELED INTERNAL BUSTION ENGINE TESTS A.F.Orlando 1, E.Santos, L.G.do Val

More information

Delhi School of Economics Course 001: Microeconomic Theory Solutions to Problem Set 2.

Delhi School of Economics Course 001: Microeconomic Theory Solutions to Problem Set 2. Delh School of Economcs Corse 00: Mcroeconomc Theor Soltons to Problem Set.. Propertes of % extend to and. (a) Assme x x x. Ths mples: x % x % x ) x % x. Now sppose x % x. Combned wth x % x and transtvt

More information

SEISMIC ANALYSIS AND DESIGN OF BUILDING STRUCTURES WITH SUPPLEMENTAL LEAD DAMPERS

SEISMIC ANALYSIS AND DESIGN OF BUILDING STRUCTURES WITH SUPPLEMENTAL LEAD DAMPERS SEISMIC ANALYSIS AND DESIGN OF BUILDING STUCTUES WITH SUPPLEMENTAL LEAD DAMPES X LIN 1, Peter J MOSS And Athol J CA 3 SUMMAY The response behavour of buldng structures ncorporatng a new type of dampng

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

, + &. / / 0 # 1) (5

, + &. / / 0 # 1) (5 ! # % # & & + + +, + &. / / 0 # 23 4 55 006 2 54223 7 / / 2 8 Hajrasoulha I & Doostan A 200 A smplfed model for sesmc response predcton of concentrcally braced frames. Advances n Engneerng Software, 43,

More information

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS) Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998

More information

Rehabilitation of existing structures by optimal placement of viscous dampers

Rehabilitation of existing structures by optimal placement of viscous dampers Lfe-Cycle and Sustanablty of Cvl Infrastructure Systems Strauss, Frangopol & Bergmester (Eds) 3 Taylor & Francs Group, London, ISBN 978--45-66-7 Rehabltaton of exstng structures by optmal placement of

More information

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.

More information

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

The Minimum Universal Cost Flow in an Infeasible Flow Network

The Minimum Universal Cost Flow in an Infeasible Flow Network Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran

More information

Robust Design Optimization of Building Frames using Order Statistics and Local Search

Robust Design Optimization of Building Frames using Order Statistics and Local Search The Eghth Chna-Japan-Korea Jont Symposum on Optmzaton of Structural and Mechancal Systems May 5-9, 04, Gyeongju, Korea aper No 0059 Robust Desgn Optmzaton of Buldng Frames usng Order Statstcs and Local

More information

OFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES

OFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES ICAMS 204 5 th Internatonal Conference on Advanced Materals and Systems OFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES VLAD LUPĂŞTEANU, NICOLAE ŢĂRANU, RALUCA HOHAN, PAUL CIOBANU Gh. Asach Techncal Unversty

More information

ORIGIN 1. PTC_CE_BSD_3.2_us_mp.mcdx. Mathcad Enabled Content 2011 Knovel Corp.

ORIGIN 1. PTC_CE_BSD_3.2_us_mp.mcdx. Mathcad Enabled Content 2011 Knovel Corp. Clck to Vew Mathcad Document 2011 Knovel Corp. Buldng Structural Desgn. homas P. Magner, P.E. 2011 Parametrc echnology Corp. Chapter 3: Renforced Concrete Slabs and Beams 3.2 Renforced Concrete Beams -

More information

EQUATION CHAPTER 1 SECTION 1STRAIN IN A CONTINUOUS MEDIUM

EQUATION CHAPTER 1 SECTION 1STRAIN IN A CONTINUOUS MEDIUM EQUTION HPTER SETION STRIN IN ONTINUOUS MEIUM ontent Introdcton One dmensonal stran Two-dmensonal stran Three-dmensonal stran ondtons for homogenety n two-dmensons n eample of deformaton of a lne Infntesmal

More information

A comprehensive study: Boundary conditions for representative volume elements (RVE) of composites

A comprehensive study: Boundary conditions for representative volume elements (RVE) of composites Insttute of Structural Mechancs A comprehensve study: Boundary condtons for representatve volume elements (RVE) of compostes Srhar Kurukur A techncal report on homogenzaton technques A comprehensve study:

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

C PLANE ELASTICITY PROBLEM FORMULATIONS

C PLANE ELASTICITY PROBLEM FORMULATIONS C M.. Tamn, CSMLab, UTM Corse Content: A ITRODUCTIO AD OVERVIEW mercal method and Compter-Aded Engneerng; Phscal problems; Mathematcal models; Fnte element method. B REVIEW OF -D FORMULATIOS Elements and

More information

C PLANE ELASTICITY PROBLEM FORMULATIONS

C PLANE ELASTICITY PROBLEM FORMULATIONS C M.. Tamn, CSMLab, UTM Corse Content: A ITRODUCTIO AD OVERVIEW mercal method and Compter-Aded Engneerng; Phscal problems; Mathematcal models; Fnte element method. B REVIEW OF -D FORMULATIOS Elements and

More information

Geometric Misfitting in Structures An Interval-Based Approach

Geometric Misfitting in Structures An Interval-Based Approach Geometrc Msfttng n Strctres An Interval-Based Approach M.V.Rama Rao, R.L.Mhanna and R.L.Mllen 3 Vasav College of Engneerng, Hyderabad-3, INDIA dr.mvrr@gmal.com School of Cvl and Envronmental Engneerng,

More information

Computer Graphics. Curves and Surfaces. Hermite/Bezier Curves, (B-)Splines, and NURBS. By Ulf Assarsson

Computer Graphics. Curves and Surfaces. Hermite/Bezier Curves, (B-)Splines, and NURBS. By Ulf Assarsson Compter Graphcs Crves and Srfaces Hermte/Bezer Crves, (B-)Splnes, and NURBS By Ulf Assarsson Most of the materal s orgnally made by Edward Angel and s adapted to ths corse by Ulf Assarsson. Some materal

More information

829. An adaptive method for inertia force identification in cantilever under moving mass

829. An adaptive method for inertia force identification in cantilever under moving mass 89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,

More information

Numerical Transient Heat Conduction Experiment

Numerical Transient Heat Conduction Experiment Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

Module 3: Element Properties Lecture 1: Natural Coordinates

Module 3: Element Properties Lecture 1: Natural Coordinates Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers

More information

The Bellman Equation

The Bellman Equation The Bellman Eqaton Reza Shadmehr In ths docment I wll rovde an elanaton of the Bellman eqaton, whch s a method for otmzng a cost fncton and arrvng at a control olcy.. Eamle of a game Sose that or states

More information

HOMOGENEOUS LEAST SQUARES PROBLEM

HOMOGENEOUS LEAST SQUARES PROBLEM the Photogrammetrc Jornal of nl, Vol. 9, No., 005 HOMOGENEOU LE QURE PROLEM Kejo Inklä Helsnk Unversty of echnology Laboratory of Photogrammetry Remote ensng P.O.ox 00, IN-005 KK kejo.nkla@tkk.f RC he

More information

DESIGN SPECTRUM-BASED SCALING OF STRENGTH REDUCTION FACTORS

DESIGN SPECTRUM-BASED SCALING OF STRENGTH REDUCTION FACTORS 13 th World Conference on Earthquake Engneerng Vancouver, B.C., Canada August 1-6, 2004 Paper No. 539 DESIGN SPECTRUM-BASED SCALING OF STRENGTH REDUCTION FACTORS Arndam CHAKRABORTI 1 and Vnay K. GUPTA

More information

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

Chapter 11: Simple Linear Regression and Correlation

Chapter 11: Simple Linear Regression and Correlation Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests

More information

Advanced Mechanical Elements

Advanced Mechanical Elements May 3, 08 Advanced Mechancal Elements (Lecture 7) Knematc analyss and moton control of underactuated mechansms wth elastc elements - Moton control of underactuated mechansms constraned by elastc elements

More information

Thin-Walled Structures Group

Thin-Walled Structures Group Thn-Walled Structures Group JOHNS HOPKINS UNIVERSITY RESEARCH REPORT Towards optmzaton of CFS beam-column ndustry sectons TWG-RR02-12 Y. Shfferaw July 2012 1 Ths report was prepared ndependently, but was

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information