Probabilistic assessment of dynamic instability of frame structures under seismic excitations

Size: px
Start display at page:

Download "Probabilistic assessment of dynamic instability of frame structures under seismic excitations"

Transcription

1 Probabilisic assessmen of dynamic insabiliy of frame srucures under seismic exciaions D. Deniz & J. Song Universiy of Illinois a Urbana-Champaign, Urbana, IL, USA J.F. Hajjar & T.H. Nguyen Norheasern Universiy, Boson, MA, USA ABSTRACT: Miigaion of collapses of srucural sysems caused by a srong earhquake shaking is crucial o reduce he poenial casualies, injuries and economic losses. Hence, accurae risk assessmen of srucural collapse under seismic exciaions is criical in effors o promoe hazard-resilience of he sociey. This paper summarizes he auhors recen effors for accurae predicion of srucural collapse wih sysemaic incorporaion of uncerainy. Compuaional simulaion models are developed for collapse es frames in he lieraure and validaed using experimenal daa. Through he validaed compuaional simulaions of collapse, alernaive collapse crieria are proposed in erms of dynamic insabiliy, i.e. he loss of he abiliy o susain he graviy loads. Using he new collapse crieria, key parameers ha govern collapse capaciy and collapse limi sae funcions are idenified for more effecive risk assessmen. A probabilisic framework is also being developed for sysemaic reamen of uncerainies in he ground moion ime hisories and for risk-informed design of frame srucures under earhquake hazards. 1 INTRODUCTION Hazard-resilience of a srucure can be defined as he abiliy o recover he full funcionaliy quickly afer a hazardous even. In effors o enhance he resilience of srucures agains earhquake hazards, collapse prevenion is hus considered one of he mos imporan objecives. Therefore, i has become crucial o undersand he causes and effecs of seismic collapse of srucures in order o develop key documens such as naional building codes, regional emergency response plans, and risk managemen sraegies. While numerous research effors are repored in he lieraure o esimae collapse capaciy of srucures, accurae predicion of srucural collapse wih sysemaic incorporaion of uncerainy sill remains elusive. This is mainly due o he lack of collapse crieria based on dynamic insabiliy of he srucure, which are needed o idenify he impac of uncerainies in earhquake loads on he dynamic insabiliy. The incremenal dynamic analysis (IDA) approach (Vamvasikos & Cornell 2002, 2004, Zareian & Krawinkler 2007, Liel e al. 2009), and a similar approach adoped by a recen projec of he Applied Technology Council (ATC-63, 2009) are considered sae-of-he-ar approaches o accoun for he uncerainy in nonlinear dynamic response, including collapse predicion. However, here have been no horough invesigaions using validaed compuaional simulaions regarding he impacs of a srucural model selecion and he seleced se of he ground moions on he predicions of collapse capaciy, and poenial conribuions of various performance measures or predicive parameers o collapse predicions. The IDA-based approach idenifies he collapse capaciy of a srucural sysem for a given earhquake ground moion based on he behavior of he IDA curve, which is he relaionship beween an engineering demand parameer (EDP) (e.g., maximum iner-sory drif raio) and an inensiy measure (IM) (e.g., specral acceleraion of an earhquake ground moion) idenified by nonlinear dynamic analyses a incremenally increased inensiy levels. The main premise of he approach is ha he srucural sysem collapses when he IDA curve becomes almos fla, i.e. a sligh increase in IM causes an exceedingly large increase in EDP. However, he IDA curve could flaen due o large residual EDPs and hus may no necessarily indicae he inabiliy o susain graviy loads. I is also noed ha mos of he recen research effors based on he IDA-based approach used only one EDP (mosly maximum sory drif raio) while one migh need alernaive EDP or

2 muliple EDPs o predic he collapse more accuraely. I is also noeworhy ha he collapse capaciy of a srucure evaluaed by he IDA-approach may be sensiive o a paricular selecion of ground moions as well as possible chaoic behavior of he IDA curve such as srucural resurrecion. Alhough some deerminisic rules have been proposed o handle such unusual behaviors of IDA (Vamvasikos & Cornell 2002, 2004), i appears ha here is a need for developing collapse crieria based on simulaed collapse phenomena and a sysemaic procedure o idenify impac of uncerainies in ground moions on he collapse capaciy of a srucural sysem. In order o overcome hese challenges, a new probabilisic framework has been developed for accurae assessmen of collapse of frame srucures under sochasic ground moions. Firs, nonlinear dynamic analyses are performed for seleced experimenal case sudies repored in he lieraure (Kanvinde 2003, Rodgers & Mahin 2004, Lignos e al. 2008) by use of OpenSees, an objec-oriened sofware framework developed by Pacific Earhquake Engineering Cener (PEER). Using OpenSees compuaional models validaed by corresponding experimenal resuls, new dynamic-insabiliy-based collapse crieria are developed in erms of energy from he inpu ground moions and he graviy loads. The seleced case sudies are hen used o es he new collapse crieria and o idenify key parameers ha govern he collapse of a srucural sysem. Currenly, procedures are being developed for probabilisic predicion of collapse limi saes and corresponding srucural demands, and for idenificaion of criical parameers in collapse predicions. Using hese procedures, he impac of uncerain ground moion deails on he collapse limi saes and srucural demands is being invesigaed o provide guidelines on selecion of ground moions for IDAbased sudies and designs. This paper presens ongoing research aciviies in he proposed framework. Firs, he deails of compuaional simulaions of collapse es cases in he lieraure are presened. Second, he proposed collapse crieria based on dynamic insabiliy, i.e. he loss of he abiliy o susain he graviy loads are inroduced. Using he seleced experimenal case sudies, he new collapse crieria are compared wih he radiional IDA-based approach. Finally, he paper inroduces ongoing research aciviies for probabilisic assessmen of collapse limi saes, srucural demands a collapse levels and for idenificaion of criical parameers for collapse predicions. 2 VALIDATED COMPUTATIONAL SIMULATION OF COLLAPSE In order o develop he framework described in he previous secion, i is necessary o build compuaional simulaion models of collapse ha are validaed by available experimenal es resuls. This secion briefly describes he simulaion ool and provides deails of he compuaional simulaion models of he seleced collapse-case sudies considered in his sudy. 2.1 Simulaion ool OpenSees The Open Sysem for Earhquake Engineering Simulaion is a finie-elemen program developed by PEER o simulae he seismic behavior of srucural and geoechnical sysems (PEER 2004). OpenSees has been exensively used by many researchers in earhquake engineering for various finie-elemen applicaions because of is advanced capabiliies in consiuive models, elemens and soluion algorihms. Moreover, i is open-source sofware plaform ha allows researchers o conribue o he framework. Therefore, his sudy uses Open- Sees as he simulaion ool o perform nonlinear dynamic collapse analysis for seleced case sudies for which collapse or near-collapse level experimenal resuls are available. Advanced high-fideliy analyical models can accoun for many facors needed for accurae simulaion of srucural collapse process, bu are compuaionally demanding. These models are prone o convergence problems during nonlinear dynamic analysis because of he complexiy of he modeling deails. Therefore, hese advanced models can be impracical for he aforemenioned sochasic framework especially in design process. For his reason, his research employs macro-models available in OpenSees ha correlae well wih experimen resuls of seleced case sudies of collapse o achieve he aforemenioned research objecives and perform large-scale parameric sudies in collapse assessmen of srucures. 2.2 Case sudies on srucural collapse Nonlinear dynamic collapse analyses are performed for seleced experimenal case sudies repored in he lieraure by use of OpenSees. So far, hree case sudies have been explored using advanced capabiliies of OpenSees in consiuive models, elemens and soluion algorihms: Kanvinde (2003), Rodgers & Mahin (2004), and Lignos e al. (2008). In developing he compuaional simulaion models of he seleced case sudies, an emphasis was given on validaion of collapse a he sysem level by considering he maximum and residual sory drif responses as well as a he componen level by considering he momen-roaion response obained a he plasic locaions a he elemen ends. Two of hese case sudies (Kanvinde 2003, and Lignos e al. 2008) are presened in his paper for describing hese research aciviies.

3 Figure 1. Specimen configuraion (Kanvinde 2003). Figure 2. Displacemen ime hisory resuls by OpenSees for he es case by Kanvinde (2003) under he es ground moion record of he 1994 Norhridge a Obregon Park. sof and radi in he plo refer o he sofening amoun and radius of ransiion from elasic o plasic branches in he springs respecively. Kanvinde (2003) conduced shake able ess on a single-sory specimen configuraion measured 12 by 24 in plan (he longer dimension aligned in he direcion of moion) and 10 in heigh o invesigae he concep of dynamic insabiliy of srucures during earhquakes. The specimen configuraion was in he form of four fla columns conneced o he base plae and a seel mass on op served as a rigid diaphragm as shown in Figure 1. For he proposed asks described in his paper, a srucural model was buil in OpenSees using he 2-D analyical model deails given in Kanvinde (2003). Elasic elemens were assigned o he columns and beam, and he beam was assumed o behave rigidly. Concenraed nodal masses were placed a he ends of op beam. Inelasic single-degree-of-freedom (SDOF) zero-lengh roaional springs were modeled a he plasic locaions a he ends of he columns by assuming Giufré- Menegoo-Pino plasiciy model (Menegoo & Pino 1973) for he spring hysereic response. The co-roaional formulaion was used in order o include he nonlinear geomeric effecs hrough he specimen. Nonlinear dynamic collapse analyses under he es ground moion of Obregon Park were performed o provide he OpenSees resul for he es case no. 11 in Figure 2, which almos coincides wih available experimen daa even a he near-collapse level. Lignos e al. (2008) performed a series of collapse shake-able ess of a 4-sory, 2-bay seel frame wih reduced-beam secions (RBS) in 1/8 scale. Figure 3 shows he seup of he es frame on he NEES mass simulaor a he Universiy a Buffalo, which consiss of elasic members wih plasic hinges a he ends. The mass simulaor is conneced o he es frame by means of axially rigid horizonal links hrough which he simulaor ransfers P-Dela effecs acing as a leaning column on he es frame. An analyical model for he 1/8 scale 4-sory es frame was developed in OpenSees based on he deerioraion parameers and mahemaical model properies given in Lignos e al. (2008). The roaional springs were used o analyically model he plasic hinges in he frame wih a modified Ibarra-Krawinkler deerioraion model (Lignos e al. 2008) available in OpenSees, calibraed based on a seel componen daabase of seel beams wih RBS under cyclic loading. Panel zones were modeled a he connecions considering he shear disorions. Furhermore, offses from he panel zones were applied o ake RBS ino accoun following he mehod used by he researchers. Effecs of he panel zones on he srucural response were explored comparing o hose of a developed clear span model. Time hisory analysis and IDA were performed using he OpenSees model and he resuls were compared wih hose by experimen (See Figure 4), which show good agreemen. Noe ha a se of nonlinear dynamic analyses wih four differen levels of seismic inensiy scales were performed sequenially o simulae he acual loading sequences of he ess. Figure 3. Shake-able-es of a 1/8 scale 4-sory, 2-bay seel frame wih reduced beam secions (Lignos e al. 2008).

4 corresponding experimenal resuls, new dynamicinsabiliy-based collapse crieria have been developed in erms of he energy from he inpu ground moions and he graviy loads. The seleced case sudies are hen used o es he new collapse crieria of a srucural sysem. Figure 4. Comparison of experimenal es resuls and simulaion resuls of laeral displacemen ime hisory a he op of he frame for he es case by Lignos e al. (2008). Noe ha he simulaion model here depends on he clear span model, and was coninuously subjeced o he ground moion record of he 1994 Norhridge earhquake a Canoga Park wih a scale facor of 0.4, 1.0, 1.5,1.9, and 2.2 following he es procedure. 2.3 Virual collapse simulaions In order o develop a new sochasic framework for idenifying collapse limi-sae and imporan parameers in he collapse assessmen of srucures subjeced o seismic loads, exensive IDAs are performed using validaed OpenSees simulaion models o obain a large sample for muliple DMs and corresponding IM and for muliple ground moions. Since he ground moions considered in he mehodology of ATC-63 projec were seleced in such a way ha he mehodology can be generally applied o building srucures a any sie, he Far-Field record se of ATC-63 projec have been chosen in he developmen of he sochasic framework. This record se consiss of weny-wo ground moion pairs (wolaeral componens) recorded a sies locaed wihin 10km of faul rupure. Records were seleced from srong earhquake ground moions wih a magniude changing from 6.5 o 7.9. As a fuure work in he sudy, virual collapse simulaions considering a wide array of geomeric and maerial parameers (Seelman & Hajjar 2009) will also be performed using he validaed analyical models o accoun for he impacs of a srucural model selecion on he collapse predicion of srucures. 3 COLLAPSE CRITERIA BASED ON DYNAMIC INSTABILITY OF FRAME STRUCTURES This secion presens new collapse crieria of frame srucures based on he dynamic insabiliy, i.e. he loss of he abiliy o susain he graviy loads. Using he OpenSees compuaional models validaed by 3.1 Tradiional IDA-based collapse limi saes The flaening of IDA curves may be caused by large residual displacemen, no by he occurrence of collapse necessarily. I is also noed ha unusual behaviors of IDA curves such as non-monoonic behavior and disconinuiies may occur (Vamvasikos & Cornell 2002). As a resul, he collapse capaciy idenified from he IDA-based approach migh no accoun for he impac of he uncerainies in ground moion appropriaely. To idenify collapse capaciies from he behavior of IDA curves, some rules have been proposed (Vamvasikos & Cornell 2002). According o he IM-based rule (Figure 5a), a building reaches is collapse capaciy when he slope of he curve is reduced o 20% of he iniial slope. If an IDA curve does no fulfill he IM-based rule, hen one checks if he drif raio exceeds an assumed global drif capaciy, say 10% (DM-based rule; see Figure 5b). However, hese IDA-based collapse idenificaion rules depend on assumed hreshold values on IM and DM, herefore no sufficien for objecive and physicsbased idenificaion of a srucural collapse based on acual dynamic insabiliy of a srucure. Therefore, new dynamic-insabiliy-based collapse crieria are developed in erms of energies from he inpu ground moions, and he graviy loads. Firs, energy balance of a srucural sysem under seismic exciaion is inroduced in he following secion. Then, deails of he new dynamic-insabiliy-based collapse crieria are presened. (a) (b) Figure 5. a) IM-based and b) DM-based rules.

5 3.2 Energy balance of a srucural sysem subjeced o seismic forces The equaion of moion a ime for a muli-degreeof-freedom (MDOF) srucure under horizonal earhquake loads and graviy loads is given as M u() + C u() + F S () = M u F () (1) where u is he relaive nodal displacemen vecor, u is he relaive nodal velociy vecor, u is he relaive nodal acceleraion vecor, M is he srucural mass marix, C is he srucural damping marix, F S is he srucural resoring force vecor, and u F is he acceleraion vecor of he applied loads. An insigh ino he dynamic insabiliy of srucures can be gained by considering energy balance of a srucural sysem under dynamic and graviy forces. If each erm in Equaion 1 is inegraed wih respec o u, he energy balance of he srucural sysem can be derived as (Uang & Berero 1990): M u() du + C u() 0 du + F 0 S () du = M u 0 F () du (2) 0 The inegrals in Equaion 2 give he following energy componens of a srucural sysem respecively, i.e. EK + ED + ES = EI (3) where EK is he relaive kineic energy, ED is he damping energy, ES is he srain energy, and EI is he relaive dynamic inpu energy. Using du = u()d, he inpu energy componen can be decomposed as follows: E I = M u F () u() d (4) 0 If he srucure is subjeced o he horizonal earhquake exciaion and he graviy loads only, hen he inpu energy can be separaed ino dynamic inpu energy due o horizonal seismic acions, EEQ, and graviy energy due o applied graviy loads on he srucure, EG, i.e. E EQ = M u () u EQ x () d (5) 0 E G = M g u 0 y () d = Mgu y () (6) Where u EQ is he acceleraion of he horizonal earhquake exciaion, u x is he relaive nodal velociy vecor in x-direcion, u y is he relaive nodal displacemen vecor in y-direcion, and g denoes he graviy acceleraion, which is consan. The earhquake energy applied on he srucure is dissipaed hrough he work done by he damping and hysereic forces. Therefore, he damping and hyseric energies are irrecoverable while he elasic and kineic energies are recoverable vibraional energy. If all he individual energy componens are gahered ogeher, he energy balance of a srucure in Equaion 3 is now defined as: EK + ED + ES = EG + EEQ (7) Akiyama (2002) saed ha he graviy energy can be considered as a release of he poenial energy as a resul of he P-Dela effecs, and akes a par in he oal resisance of a srucure agains a seismic exciaion. Therefore, he graviy energy can be alernaively shown on he lef side of he energy balance: EK + ED + ES EG = EEQ (8) 3.3 Idenificaion of dynamic insabiliy by srucural graviy energy Dynamic insabiliy is a complex phenomenon, which canno be effecively prediced by a ground inensiy measure and/or an engineering parameer roughly represening srucural damage. The mos commonly used crierion for idenificaion of simulaed collapse under dynamic loads is boundless drifs owards collapse. However, his approach requires checking displacemen demands a each degree-of-freedom (DOF), bu mos sudies consider only he roof or sory drifs o check he sabiliy of he global srucural behavior. Therefore, he dynamic insabiliy of srucural sysems is invesigaed from he viewpoin of he energy balance because energy is an overall indicaor of he sysems. Progressive accumulaion of permanen laeral drifs during a srong ground shaking may render graviy forces he dominan forces and make he srucure collapse under significan P-dela effecs due o governing graviy forces (Jennings & Husid 1968). Therefore, new collapse crieria have been developed based on comparison beween he amoun of dynamic energy released by he earhquake o he srucure and he amoun of graviaional work done by he verical saic loads during he dynamic analysis. The sudden increase of he graviy energy, which evenually causes he graviy energy o exceed he dynamic energy, can be considered as an indicaor of he dominaion of graviy loads over dynamic loads. Figure 6 presens he inpu-energy-ime hisories for he validaed single-sory model of Kanvinde (2003) under he ground moion record of he 1994 Norhridge earhquake a Obregon Park, Los Angeles. If he inensiy of he ground moion is no srong enough o rigger he large geomeric effecs in he frame (e.g. he non-collapse case a he scale of 0.8), he srucure obains a seady sae in erms of he graviaional energy, which is found insignifican comparing o he quaniy of dynamic inpu energy coming from he ground moion (Figure 6a). A he inensiy scale of 1.0 (Figure 6b), geomeric non-

6 lineariies in he srucure become significan near he collapse, causing he frame o show very large displacemen in verical direcions and hus resul in a sudden increase in he graviaional energy as he srucure ges close o collapse. Since his approach employs sysem-level measures, i.e. graviy and dynamic inpu energies, one does no need o check each degree-of-freedom of he srucure o check he dynamic insabiliy. Moreover, he approach may faciliae developing a mahemaical descripion of dynamic insabiliy, which can be paricularly useful as limi-sae funcions during srucural reliabiliy analysis of collapse. (a) level of ground moion a which he srucure loses he dynamic sabiliy as he collapse capaciy. However, he srucural collapse capaciy should be evaluaed based on he maximum inensiy level before occurrence of dynamic insabiliy (Krawinkler e al. 2009, Haselon e al. 2009). The capaciy a his inensiy level is acual represenaion of he larges srucural resisance agains dynamic collapse. The validaed model for he case sudy by Lignos e al. (2008) is uilized here o perform nonlinear dynamic analyses using he ATC-63 far field se. Figure 7 shows he IDA curves of peak ground acceleraion (PGA) o op laeral displacemen obained from he validaed OpenSees model. Tradiional IDA-based rules are compared o he new crieria called energy rule based on he maximum inensiy level observed before he dynamic insabiliy occurs, i.e., graviy energy exceeds dynamic energy. Much variabiliy is observed in collapse capaciy level for all rules due o he effec of randomness in he seleced ground moions on srucural collapse. I is also noed ha he radiional IDA-based approach and he proposed crieria resul in significan differences in predicions of collapse capaciy (in PGA) and he srucural demand a he collapse level (in roof displacemen). (b) Figure 6. Inpu energy componens near collapse under he es earhquake of 1994 Norhridge earhquake a Obregon Park a) non-collapse case a he ground moion scale of 0.8 and b) collapse case a he ground moion scale of Comparison of new collapse crieria wih radiional IDA-based rules Tradiional IDA-based rules are compared o he new crieria called energy rule in erms of he maximum inensiy level observed before he dynamic insabiliy occurs, i.e., graviy energy exceeds dynamic energy. Mos recen research effors based on he IDA-based approach assume he inensiy Figure 7. IDA curves and comparison of collapse crieria for he es case of Lignos e al. (2008) subjeced o ATC-63 far field record se (44 ground moions). 4 PROBABILISTIC ASSESSMENT OF COLLAPSE USING NEW COLLAPSE CRITERIA Collapse assessmen of srucures based on IDA curves depends on he selecion of IM and DM used o consruc hese IDA curves as well as variabiliy in he se of ground moions considered in he analysis (Villaverde 2007). IDA curves usually reach a fla plaeau as an indicaion of collapse (i.e., a large increase in he srucural response corresponding o a small increase in he ground moion inensiy), bu

7 his plaeau may occur a several differen inensiy levels of ground moions. Large dispersion observed in collapse-causing inensiies of ground moions and damage hresholds has iniiaed search for alernaive performance measures o assess collapse capaciy of srucures. However, mos of recen research on he IDA-based approach sill uses only one DM and one IM (mosly maximum sory drif raio and elasic specral acceleraion). This secion describes ongoing research effors o obain opimal selecion and/or combinaion of muliple performance measures ha describe he limi-sae mos effecively and he benefi of having more han one DM or IM for collapse capaciy predicion. 4.1 Saisical analysis on collapse capaciy by new crieria The collapse capaciy daa obained for he case sudy by Lignos e al. (2008) using ATC-63 far field se was invesigaed here o perform saisical analyses on he levels of IMs (denoed by IMcol) and he condiional disribuions of he corresponding DMs (denoed by DMcol) Empirical cumulaive disribuion funcions (CDF) for IMcol idenified by he energy rule are presened in Figure 8. CDFs of differen kinds of IMcol (normalized by heir means) such as peak ground acceleraion-velociy-displacemen (PGA- PGV-PGD), specral acceleraion for 2% damping (Sa (2%)), undamped inensiy of Arias (IA (0%), Arias 1970), roo of inegral of square of ground acceleraion-ime hisory (ars), average cycle of crossings over zero in he srong par of accelerogram (Tv, srong, where srong duraion is based on definiion by Trifunac & Brady 1975), and earhquake inpu energy (EEQ) were compared. These CDFs provide probabilisic esimaion of he collapse capaciies (in IM). In order o make a numerical comparison, heir normalized sandard deviaion, which is coefficien of variaion (cov), were obained in Table 1. IM wih he smalles cov is found o be ars from which collapse can be prediced wih more confidence comparing o oher candidaes. Currenly, he following linear combinaions of IMs and/or DMs in Equaion 9 or nonlinear combinaions are also being invesigaed o obain new performance measures (PM) ha can predic he collapse capaciies wih smaller covs: n PM = a i IM i + b i DM i (9) i=1 where a i and b i are he coefficiens of he IM and DM ha resul in smalles covs, found by saisical analysis. Condiional disribuion of srucural demand (DMcol) given IMcol is also being obained. Parial descripors such as condiional mean and variance of DMcol given IMcol can be obained hrough linear/nonlinear regression analysis. Such regression models provide probabilisic esimaion of he srucural demand a near-collapse level. 4.2 Effec of variabiliy in ground moions Effecs of uncerain characerisics of ground moions used for IDA on he probabilisic models described in secion 4.1 are being invesigaed o develop guidelines regarding selecion of a suie of ground moions o be used in nonlinear collapse analyses. Moreover, a new analysis framework is being developed for probabilisic evaluaion of collapse using arificial models represening possible ground moions a a given sie insead of a suie of seleced ground moions. This approach allows us o perform sie-specific probabilisic collapse assessmen. The approach will perform random vibraion analysis by employing a discree represenaion of a ground moion daabase (Rezaeian & Der Kiureghian 2008, 2010). The developed probabilisic mehod will be inegraed wih performance-based earhquake engineering (PBEE) framework, which will provide collapse fragiliy models hrough sysemaic reamen of uncerainies in seismic capaciy, demand and models. Figure 8. Cumulaive disribuion funcion of IM-collapse level for he es case of Lignos e al. (2008) subjeced o ATC-63 far field record se (44 ground moions). Table 1. Collapse-iniiaing inensiy measure and coefficien of variaions (cov) IM col cov ars 0.31 T v, srong 0.33 PGV 0.41 PGA 0.42 E EQ 0.48 Sa (2%) 0.56 IA 0.68 PGD 0.87

8 5 CONCLUSIONS AND FUTURE WORK New collapse crieria were developed based on dynamic insabiliy of srucures near collapse hrough validaed compuaional simulaions of collapse es frames repored in he lieraure. Using he crieria, various quaniaive crieria are being explored o faciliae performing srucural reliabiliy analysis and idenifying criical measures. New collapse crieria based on compuaional simulaions of dynamic insabiliy give rise o new research opporuniies o gain beer undersanding of complex collapse of srucural sysems, idenify key parameers requiring more comprehensive and accurae measuremens during experimens, achieve more accurae and sysemaic predicion of collapse, and incorporae uncerainies ino collapse predicion. The sudy described here is expeced o have poenial impac across several srucural engineering research and pracice consiuencies seeking o improve building code provisions for prevening disproporionae collapse; regional emergency response plans and risk managemen sraegies ha rely on accurae assessmen of collapse wihin fragiliy analysis; and collapse assessmen of new srucural sysems. In addiion, hrough his work, life safey will be enhanced, as avoiding srucural collapse due o exreme loads is a criical componen o ensuring a safe infrasrucure and is resilience agains earhquake hazards. 6 ACKNOWLDEGEMENT The auhors would like o hank Naional Science Foundaion (NSF) for heir funding in his research under gran number CMMI Any opinions, findings and conclusions or recommendaions expressed in his maerial are hose of he auhors and do no necessarily reflec he views of NSF. REFERENCES Akiyama, H Collapse modes of srucures under srong moions of earhquake. Annals of Geophysics 45(6): Applied Technology Council, Quanificaion of building seismic performance facors. ATC-63 Projec Repor FE- MA. Arias, A A measure of earhquake inensiy. In R. J. Hansen (ed), Seismic design for nuclear power plans: Cambridge: M. I. T. Press. Haselon, C.B, Liel, A.B., & Deierlein, G.G Simulaing srucural collapse due o earhquakes: Model idealizaion, model Calibraion, and numerical soluion algorihms. 2nd Inernaional Conference on Compuaional Mehods in Srucural Dynamics and Earhquake Engineering, Rhodes, Greece, June. Jennings, P. C. & Husid, R Collapse of yielding srucures during earhquakes. Journal of Engineering Mechanics Division 94(5): Kanvinde, M Mehods o evaluae he dynamic sabiliy of srucures shake able ess and nonlinear dynamic analyses. EERI Annual Suden Paper Compeiion, Earhquake Engineering Research Insiue, Oakland, Calif. Krawinkler, H., Zareian, F., Lignos, D., & Ibarra, L. F Predicion of collapse of srucures under earhquake exciaions. 2nd Inernaional Conference on Compuaional Mehods in Srucural Dynamics and Earhquake Engineering, Rhodes, Greece, June. Liel, A.B., Haselon, C.B., Deierlein, G.G., & Baker, J.W Incorporaing modeling uncerainies in he assessmen of seismic collapse risk of buildings. Srucural Safey 31: Lignos, D.G., Krawinkler, H., & Whiaker, A.S Shaking able collapse ess of wo scale models of a 4-sory momen resising seel frame. Proceedings of he 14h World Conference on Earhquake Engineering, Beijing, China, Ocober. Menegoo, M., and Pino, E Mehod of analysis for cyclically loaded reinforced concree plane frames including changes in geomery and non-elasic behavior of elemens under combined normal force and bending. Proceedings of IABSE Symposium, Lisbon, Porugal. Pacific Earhquake Engineering Research Cener Open Sysem for Earhquake Engineering Simulaion. Universiy of California, Berkeley, CA. Rezaeian, S. & Der Kiureghian, A A sochasic ground moion model wih separable emporal and specral nonsaionariies. Earhquake Engineering and Srucural Dynamics 37: Rezaeian, S. & Der Kiureghian, A Simulaion of synheic ground moions for specified earhquake and sie characerisics. Earhquake Engineering and Srucural Dynamics 39: Rodgers, J.E., & Mahin, S.A Effecs of connecion hysereic degradaion on he seismic behavior of seel momen-resising frame. Pacfic Earhquake Engineering Research Cener Repor, PEER 2003/13. Seelman, J. S. & Hajjar, J. F Influence of inelasic seismic response modeling on regional loss esimaion. Engineering Srucures 31(12): Trifunac, M.D., & Brady, A.G On he correlaion of seismic inensiy scales wih he peaks of recorded ground moion. Bullein of he Seismological Sociey of America 65(1): Uang, C. M. & Berero, V. V Evaluaion of seismic energy in srucures. Earhquake Engineering and Srucural Dynamics 19: Vamvasikos, D. & Cornell, C.A Incremenal dynamic analysis. Earhquake Engineering and Srucural Dynamics 31 (3): Vamvasikos, D. & Cornell, C.A Applied incremenal dynamic analysis. Earhquake Specra 20(2): Villaverde, R Mehods o assess he seismic collapse capaciy of building srucures: Sae of he ar. Journal of Srucural Engineering 33(1): Zareian, F. & Krawinkler, H Assessmen of probabiliy of collapse and design for collapse safey. Earhquake Engineering and Srucural Dynamics 36(13):

Structural Dynamics and Earthquake Engineering

Structural Dynamics and Earthquake Engineering Srucural Dynamics and Earhquae Engineering Course 1 Inroducion. Single degree of freedom sysems: Equaions of moion, problem saemen, soluion mehods. Course noes are available for download a hp://www.c.up.ro/users/aurelsraan/

More information

Cumulative Damage Evaluation based on Energy Balance Equation

Cumulative Damage Evaluation based on Energy Balance Equation Cumulaive Damage Evaluaion based on Energy Balance Equaion K. Minagawa Saiama Insiue of Technology, Saiama S. Fujia Tokyo Denki Universiy, Tokyo! SUMMARY: This paper describes an evaluaion mehod for cumulaive

More information

Estimation of Kinetic Friction Coefficient for Sliding Rigid Block Nonstructural Components

Estimation of Kinetic Friction Coefficient for Sliding Rigid Block Nonstructural Components 7 Esimaion of Kineic Fricion Coefficien for Sliding Rigid Block Nonsrucural Componens Cagdas Kafali Ph.D. Candidae, School of Civil and Environmenal Engineering, Cornell Universiy Research Supervisor:

More information

Mechanical Fatigue and Load-Induced Aging of Loudspeaker Suspension. Wolfgang Klippel,

Mechanical Fatigue and Load-Induced Aging of Loudspeaker Suspension. Wolfgang Klippel, Mechanical Faigue and Load-Induced Aging of Loudspeaker Suspension Wolfgang Klippel, Insiue of Acousics and Speech Communicaion Dresden Universiy of Technology presened a he ALMA Symposium 2012, Las Vegas

More information

ANALYSIS OF REINFORCED CONCRETE BUILDINGS IN FIRE

ANALYSIS OF REINFORCED CONCRETE BUILDINGS IN FIRE ANALYSIS OF REINFORCED CONCRETE BUILDINGS IN FIRE Dr Zhaohui Huang Universiy of Sheffield 6 May 2005 1 VULCAN layered slab elemens: connecion o beam elemens Plae Elemen Slab nodes y x Reference Plane h

More information

ON THE BEAT PHENOMENON IN COUPLED SYSTEMS

ON THE BEAT PHENOMENON IN COUPLED SYSTEMS 8 h ASCE Specialy Conference on Probabilisic Mechanics and Srucural Reliabiliy PMC-38 ON THE BEAT PHENOMENON IN COUPLED SYSTEMS S. K. Yalla, Suden Member ASCE and A. Kareem, M. ASCE NaHaz Modeling Laboraory,

More information

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec

More information

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix

More information

CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK

CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 175 CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 10.1 INTRODUCTION Amongs he research work performed, he bes resuls of experimenal work are validaed wih Arificial Neural Nework. From he

More information

Combined Bending with Induced or Applied Torsion of FRP I-Section Beams

Combined Bending with Induced or Applied Torsion of FRP I-Section Beams Combined Bending wih Induced or Applied Torsion of FRP I-Secion Beams MOJTABA B. SIRJANI School of Science and Technology Norfolk Sae Universiy Norfolk, Virginia 34504 USA sirjani@nsu.edu STEA B. BONDI

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

A Specification Test for Linear Dynamic Stochastic General Equilibrium Models

A Specification Test for Linear Dynamic Stochastic General Equilibrium Models Journal of Saisical and Economeric Mehods, vol.1, no.2, 2012, 65-70 ISSN: 2241-0384 (prin), 2241-0376 (online) Scienpress Ld, 2012 A Specificaion Tes for Linear Dynamic Sochasic General Equilibrium Models

More information

20. Applications of the Genetic-Drift Model

20. Applications of the Genetic-Drift Model 0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0

More information

At the end of this lesson, the students should be able to understand

At the end of this lesson, the students should be able to understand Insrucional Objecives A he end of his lesson, he sudens should be able o undersand Sress concenraion and he facors responsible. Deerminaion of sress concenraion facor; experimenal and heoreical mehods.

More information

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

The motions of the celt on a horizontal plane with viscous friction

The motions of the celt on a horizontal plane with viscous friction The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and

More information

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

More information

( ) g. Research on energy spectrum of elastic SDOF

( ) g. Research on energy spectrum of elastic SDOF The 14 h World Conference on Earhquake Enineerin Ocober 12-17, 28, Beijin, China Research on enery specrum of elasic SDOF TENG Jun DONG Zhijun Harbin Insiue Universiy of Technoloy Shenzhen School, Shenzhen

More information

Notes on Kalman Filtering

Notes on Kalman Filtering Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren

More information

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology

More information

Curling Stress Equation for Transverse Joint Edge of a Concrete Pavement Slab Based on Finite-Element Method Analysis

Curling Stress Equation for Transverse Joint Edge of a Concrete Pavement Slab Based on Finite-Element Method Analysis TRANSPORTATION RESEARCH RECORD 155 35 Curling Sress Equaion for Transverse Join Edge of a Concree Pavemen Slab Based on Finie-Elemen Mehod Analysis TATSUO NISHIZAWA, TADASHI FUKUDA, SABURO MATSUNO, AND

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

More information

Program: RFEM 5, RSTAB 8, RF-DYNAM Pro, DYNAM Pro. Category: Isotropic Linear Elasticity, Dynamics, Member

Program: RFEM 5, RSTAB 8, RF-DYNAM Pro, DYNAM Pro. Category: Isotropic Linear Elasticity, Dynamics, Member Verificaion Example Program: RFEM 5, RSTAB 8, RF-DYNAM Pro, DYNAM Pro Caegory: Isoropic Linear Elasiciy, Dynamics, Member Verificaion Example: 0104 Canilever Beam wih Periodic Exciaion 0104 Canilever Beam

More information

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,

More information

Numerical Evaluation of an Add-On Vehicle Protection System

Numerical Evaluation of an Add-On Vehicle Protection System Numerical Evaluaion of an Add-On Vehicle Proecion Sysem Geneviève Toussain, Amal Bouamoul, Rober Durocher, Jacob Bélanger*, Benoî S-Jean Defence Research and Developmen Canada Valcarier 2459 Bravoure Road,

More information

SOFTWARE RELIABILITY GROWTH MODEL WITH LOGISTIC TESTING-EFFORT FUNCTION CONSIDERING LOG-LOGISTIC TESTING-EFFORT AND IMPERFECT DEBUGGING

SOFTWARE RELIABILITY GROWTH MODEL WITH LOGISTIC TESTING-EFFORT FUNCTION CONSIDERING LOG-LOGISTIC TESTING-EFFORT AND IMPERFECT DEBUGGING Inernaional Journal of Compuer Science and Communicaion Vol. 2, No. 2, July-December 2011, pp. 605-609 SOFTWARE RELIABILITY GROWTH MODEL WITH LOGISTIC TESTING-EFFORT FUNCTION CONSIDERING LOG-LOGISTIC TESTING-EFFORT

More information

Probabilistic Models for Reliability Analysis of a System with Three Consecutive Stages of Deterioration

Probabilistic Models for Reliability Analysis of a System with Three Consecutive Stages of Deterioration Yusuf I., Gaawa R.I. Volume, December 206 Probabilisic Models for Reliabiliy Analysis of a Sysem wih Three Consecuive Sages of Deerioraion Ibrahim Yusuf Deparmen of Mahemaical Sciences, Bayero Universiy,

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

G.F. Dargush, Y. Hu, S. Dogruel, G. Apostolakis University at Buffalo MCEER Annual Meeting Washington, DC June 29, 2006

G.F. Dargush, Y. Hu, S. Dogruel, G. Apostolakis University at Buffalo MCEER Annual Meeting Washington, DC June 29, 2006 Simulaion-based Muli-hazard Decision Suppor G.F. Dargush, Y. Hu, S. Dogruel, G. Aposolakis Universiy a Buffalo MCEER Annual Meeing Washingon, DC June 29, 2006 Simulaion-based Muli-hazard Decision Suppor

More information

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time. Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike

More information

Calculus and Seismic Analysis November Ken Mark Bechtel Project Manager

Calculus and Seismic Analysis November Ken Mark Bechtel Project Manager Calculus and Seismic Analysis November 2007 Ken Mark Bechel Projec Manager BART Today 5 Lines 104 Miles Underground 37 Miles Aerial 23 Miles Surface 44 Miles 43 Saions 330,000 Daily Ridership 150,000 Daily

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

Testing for a Single Factor Model in the Multivariate State Space Framework

Testing for a Single Factor Model in the Multivariate State Space Framework esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics

More information

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP). Appendix A: Conservaion of Mechanical Energy = Conservaion of Linear Momenum Consider he moion of a nd order mechanical sysem comprised of he fundamenal mechanical elemens: ineria or mass (M), siffness

More information

An recursive analytical technique to estimate time dependent physical parameters in the presence of noise processes

An recursive analytical technique to estimate time dependent physical parameters in the presence of noise processes WHAT IS A KALMAN FILTER An recursive analyical echnique o esimae ime dependen physical parameers in he presence of noise processes Example of a ime and frequency applicaion: Offse beween wo clocks PREDICTORS,

More information

IB Physics Kinematics Worksheet

IB Physics Kinematics Worksheet IB Physics Kinemaics Workshee Wrie full soluions and noes for muliple choice answers. Do no use a calculaor for muliple choice answers. 1. Which of he following is a correc definiion of average acceleraion?

More information

1. VELOCITY AND ACCELERATION

1. VELOCITY AND ACCELERATION 1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under

More information

Numerical Dispersion

Numerical Dispersion eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal

More information

PROGRESSIVE COLLAPSE ANALYSIS THROUGH STRENGTH DEGRADATION AND FRACTURE IN THE MIXED LAGRANGIAN FRAMEWORK

PROGRESSIVE COLLAPSE ANALYSIS THROUGH STRENGTH DEGRADATION AND FRACTURE IN THE MIXED LAGRANGIAN FRAMEWORK Ocober 1-17, 8, Beijing, China PROGRESSIVE COLLAPSE ANALYSIS THROUGH STRENGTH DEGRADATION AND RACTURE IN THE MIXED LAGRANGIAN RAMEWORK O. Lavan 1, M.V. Sivaselvan, A.M. Reinhorn 3 and G.. Dargush 4 1 Senior

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION DOI: 0.038/NCLIMATE893 Temporal resoluion and DICE * Supplemenal Informaion Alex L. Maren and Sephen C. Newbold Naional Cener for Environmenal Economics, US Environmenal Proecion

More information

OBJECTIVES OF TIME SERIES ANALYSIS

OBJECTIVES OF TIME SERIES ANALYSIS OBJECTIVES OF TIME SERIES ANALYSIS Undersanding he dynamic or imedependen srucure of he observaions of a single series (univariae analysis) Forecasing of fuure observaions Asceraining he leading, lagging

More information

Modeling the Dynamics of an Ice Tank Carriage

Modeling the Dynamics of an Ice Tank Carriage Modeling he Dynamics of an Ice Tank Carriage The challenge: To model he dynamics of an Ice Tank Carriage and idenify a mechanism o alleviae he backlash inheren in he design of he gearbox. Maplesof, a division

More information

v A Since the axial rigidity k ij is defined by P/v A, we obtain Pa 3

v A Since the axial rigidity k ij is defined by P/v A, we obtain Pa 3 The The rd rd Inernaional Conference on on Design Engineering and Science, ICDES 14 Pilsen, Czech Pilsen, Republic, Czech Augus Republic, 1 Sepember 1-, 14 In-plane and Ou-of-plane Deflecion of J-shaped

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

On Multicomponent System Reliability with Microshocks - Microdamages Type of Components Interaction

On Multicomponent System Reliability with Microshocks - Microdamages Type of Components Interaction On Mulicomponen Sysem Reliabiliy wih Microshocks - Microdamages Type of Componens Ineracion Jerzy K. Filus, and Lidia Z. Filus Absrac Consider a wo componen parallel sysem. The defined new sochasic dependences

More information

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8) I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression

More information

STUDY ON A METHOD TO IDENTIFY EXTERNAL FORCES AND DAMPING PARAMETERS OF STRUCTURE FROM RANDOM RESPONSE VALUES

STUDY ON A METHOD TO IDENTIFY EXTERNAL FORCES AND DAMPING PARAMETERS OF STRUCTURE FROM RANDOM RESPONSE VALUES STUDY ON A METHOD TO IDENTIFY ETERNAL FORCES AND DAMPING PARAMETERS OF STRUCTURE FROM RANDOM RESPONSE VALUES 114 M KANDA 1, E MARUTA, D INABA 3 And T MATSUYAMA 4 SUMMARY Proposed in his paper is a mehod

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance

More information

CHAPTER 2 Signals And Spectra

CHAPTER 2 Signals And Spectra CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par

More information

An introduction to the theory of SDDP algorithm

An introduction to the theory of SDDP algorithm An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking

More information

Sub Module 2.6. Measurement of transient temperature

Sub Module 2.6. Measurement of transient temperature Mechanical Measuremens Prof. S.P.Venkaeshan Sub Module 2.6 Measuremen of ransien emperaure Many processes of engineering relevance involve variaions wih respec o ime. The sysem properies like emperaure,

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

MEI STRUCTURED MATHEMATICS 4758

MEI STRUCTURED MATHEMATICS 4758 OXFORD CAMBRIDGE AND RSA EXAMINATIONS Advanced Subsidiary General Cerificae of Educaion Advanced General Cerificae of Educaion MEI STRUCTURED MATHEMATICS 4758 Differenial Equaions Thursday 5 JUNE 006 Afernoon

More information

Computation of the Effect of Space Harmonics on Starting Process of Induction Motors Using TSFEM

Computation of the Effect of Space Harmonics on Starting Process of Induction Motors Using TSFEM Journal of elecrical sysems Special Issue N 01 : November 2009 pp: 48-52 Compuaion of he Effec of Space Harmonics on Saring Process of Inducion Moors Using TSFEM Youcef Ouazir USTHB Laboraoire des sysèmes

More information

ACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H.

ACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H. ACE 564 Spring 2006 Lecure 7 Exensions of The Muliple Regression Model: Dumm Independen Variables b Professor Sco H. Irwin Readings: Griffihs, Hill and Judge. "Dumm Variables and Varing Coefficien Models

More information

The field of mathematics has made tremendous impact on the study of

The field of mathematics has made tremendous impact on the study of A Populaion Firing Rae Model of Reverberaory Aciviy in Neuronal Neworks Zofia Koscielniak Carnegie Mellon Universiy Menor: Dr. G. Bard Ermenrou Universiy of Pisburgh Inroducion: The field of mahemaics

More information

Finite element method for structural dynamic and stability analyses

Finite element method for structural dynamic and stability analyses Finie elemen mehod for srucural dynamic and sabiliy analyses Module- Nonlinear FE Models Lecure-39 Toal and updaed Lagrangian formulaions Prof C Manohar Deparmen of Civil Engineering IIc, Bangalore 56

More information

th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 3256

th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 3256 11111 1 h World Conference on Earhquake Engineering Vancouver, B.C., Canada Augus 1-6, 2004 Paper No. 256 TOUCHING ANALYSIS OF TWO BUILDINGS USING FINITE ELEMENT METHOD Mircea IEREMIA 1, Silviu GINJU 1,

More information

Optimal Path Planning for Flexible Redundant Robot Manipulators

Optimal Path Planning for Flexible Redundant Robot Manipulators 25 WSEAS In. Conf. on DYNAMICAL SYSEMS and CONROL, Venice, Ialy, November 2-4, 25 (pp363-368) Opimal Pah Planning for Flexible Redundan Robo Manipulaors H. HOMAEI, M. KESHMIRI Deparmen of Mechanical Engineering

More information

Inventory Control of Perishable Items in a Two-Echelon Supply Chain

Inventory Control of Perishable Items in a Two-Echelon Supply Chain Journal of Indusrial Engineering, Universiy of ehran, Special Issue,, PP. 69-77 69 Invenory Conrol of Perishable Iems in a wo-echelon Supply Chain Fariborz Jolai *, Elmira Gheisariha and Farnaz Nojavan

More information

Inequality measures for intersecting Lorenz curves: an alternative weak ordering

Inequality measures for intersecting Lorenz curves: an alternative weak ordering h Inernaional Scienific Conference Financial managemen of Firms and Financial Insiuions Osrava VŠB-TU of Osrava, Faculy of Economics, Deparmen of Finance 7 h 8 h Sepember 25 Absrac Inequaliy measures for

More information

Stable block Toeplitz matrix for the processing of multichannel seismic data

Stable block Toeplitz matrix for the processing of multichannel seismic data Indian Journal of Marine Sciences Vol. 33(3), Sepember 2004, pp. 215-219 Sable block Toepliz marix for he processing of mulichannel seismic daa Kiri Srivasava* & V P Dimri Naional Geophysical Research

More information

Flow-Induced Vibration Analysis of Supported Pipes with a Crack

Flow-Induced Vibration Analysis of Supported Pipes with a Crack Flow-Induced Vibraion Analsis of Suppored Pipes wih a Crack Jin-Huk Lee, Samer Masoud Al-Said Deparmen of Mechanical Engineering American Universi of Sharjah, UAE Ouline Inroducion and Moivaion Aeroacousicall

More information

MECHANICS OF MATERIALS Poisson s Ratio

MECHANICS OF MATERIALS Poisson s Ratio Poisson s Raio For a slender bar subjeced o axial loading: ε x x y 0 The elongaion in he x-direcion i is accompanied by a conracion in he oher direcions. Assuming ha he maerial is isoropic (no direcional

More information

Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source

Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source Muli-scale D acousic full waveform inversion wih high frequency impulsive source Vladimir N Zubov*, Universiy of Calgary, Calgary AB vzubov@ucalgaryca and Michael P Lamoureux, Universiy of Calgary, Calgary

More information

Dynamic Analysis of Loads Moving Over Structures

Dynamic Analysis of Loads Moving Over Structures h Inernaional ongress of roaian ociey of echanics epember, 18-, 3 Bizovac, roaia ynamic nalysis of Loads oving Over rucures Ivica Kožar, Ivana Šimac Keywords: moving load, direc acceleraion mehod 1. Inroducion

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

Stability and Bifurcation in a Neural Network Model with Two Delays Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy

More information

Article from. Predictive Analytics and Futurism. July 2016 Issue 13

Article from. Predictive Analytics and Futurism. July 2016 Issue 13 Aricle from Predicive Analyics and Fuurism July 6 Issue An Inroducion o Incremenal Learning By Qiang Wu and Dave Snell Machine learning provides useful ools for predicive analyics The ypical machine learning

More information

Morning Time: 1 hour 30 minutes Additional materials (enclosed):

Morning Time: 1 hour 30 minutes Additional materials (enclosed): ADVANCED GCE 78/0 MATHEMATICS (MEI) Differenial Equaions THURSDAY JANUARY 008 Morning Time: hour 30 minues Addiional maerials (enclosed): None Addiional maerials (required): Answer Bookle (8 pages) Graph

More information

Lab #2: Kinematics in 1-Dimension

Lab #2: Kinematics in 1-Dimension Reading Assignmen: Chaper 2, Secions 2-1 hrough 2-8 Lab #2: Kinemaics in 1-Dimension Inroducion: The sudy of moion is broken ino wo main areas of sudy kinemaics and dynamics. Kinemaics is he descripion

More information

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),

More information

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum MEE Engineering Mechanics II Lecure 4 Lecure 4 Kineics of a paricle Par 3: Impulse and Momenum Linear impulse and momenum Saring from he equaion of moion for a paricle of mass m which is subjeced o an

More information

Energy Storage Benchmark Problems

Energy Storage Benchmark Problems Energy Sorage Benchmark Problems Daniel F. Salas 1,3, Warren B. Powell 2,3 1 Deparmen of Chemical & Biological Engineering 2 Deparmen of Operaions Research & Financial Engineering 3 Princeon Laboraory

More information

1.1. Example: Polynomial Curve Fitting 4 1. INTRODUCTION

1.1. Example: Polynomial Curve Fitting 4 1. INTRODUCTION 4. INTRODUCTION Figure.2 Plo of a raining daa se of N = poins, shown as blue circles, each comprising an observaion of he inpu variable along wih he corresponding arge variable. The green curve shows he

More information

Seismic response analysis of the tunnel with accumulated damage and crack effect in shaking table test

Seismic response analysis of the tunnel with accumulated damage and crack effect in shaking table test Seismic response analysis of he unn wih accumulaed damage and crack effec in shaking able es ABSTRACT : Z.Z. Wang 1, B. Gao 1, T.C. Sun 1, 2, Y.S. Shen 1 1 School of Civil Engineering, Souhwes Jiaoong

More information

Basilio Bona ROBOTICA 03CFIOR 1

Basilio Bona ROBOTICA 03CFIOR 1 Indusrial Robos Kinemaics 1 Kinemaics and kinemaic funcions Kinemaics deals wih he sudy of four funcions (called kinemaic funcions or KFs) ha mahemaically ransform join variables ino caresian variables

More information

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based

More information

5 The fitting methods used in the normalization of DSD

5 The fitting methods used in the normalization of DSD The fiing mehods used in he normalizaion of DSD.1 Inroducion Sempere-Torres e al. 1994 presened a general formulaion for he DSD ha was able o reproduce and inerpre all previous sudies of DSD. The mehodology

More information

MATHEMATICAL DESCRIPTION OF THEORETICAL METHODS OF RESERVE ECONOMY OF CONSIGNMENT STORES

MATHEMATICAL DESCRIPTION OF THEORETICAL METHODS OF RESERVE ECONOMY OF CONSIGNMENT STORES MAHEMAICAL DESCIPION OF HEOEICAL MEHODS OF ESEVE ECONOMY OF CONSIGNMEN SOES Péer elek, József Cselényi, György Demeer Universiy of Miskolc, Deparmen of Maerials Handling and Logisics Absrac: Opimizaion

More information

Transverse and Longitudinal Damped Vibrations of Hydraulic Cylinder in a Mining Prop

Transverse and Longitudinal Damped Vibrations of Hydraulic Cylinder in a Mining Prop Vibraions in Physical Sysems Vol. 7 (06) Transverse and Longiudinal Damped Vibraions of Hydraulic Cylinder in a Mining Prop Wojciech SOCHACK nsiue of Mechanics and Fundamenals of Machinery Design Universiy

More information

Q.1 Define work and its unit?

Q.1 Define work and its unit? CHP # 6 ORK AND ENERGY Q.1 Define work and is uni? A. ORK I can be define as when we applied a force on a body and he body covers a disance in he direcion of force, hen we say ha work is done. I is a scalar

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

Finite Element Analysis of Structures

Finite Element Analysis of Structures KAIT OE5 Finie Elemen Analysis of rucures Mid-erm Exam, Fall 9 (p) m. As shown in Fig., we model a russ srucure of uniform area (lengh, Area Am ) subjeced o a uniform body force ( f B e x N / m ) using

More information

Mathematical Theory and Modeling ISSN (Paper) ISSN (Online) Vol 3, No.3, 2013

Mathematical Theory and Modeling ISSN (Paper) ISSN (Online) Vol 3, No.3, 2013 Mahemaical Theory and Modeling ISSN -580 (Paper) ISSN 5-05 (Online) Vol, No., 0 www.iise.org The ffec of Inverse Transformaion on he Uni Mean and Consan Variance Assumpions of a Muliplicaive rror Model

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

Air Traffic Forecast Empirical Research Based on the MCMC Method

Air Traffic Forecast Empirical Research Based on the MCMC Method Compuer and Informaion Science; Vol. 5, No. 5; 0 ISSN 93-8989 E-ISSN 93-8997 Published by Canadian Cener of Science and Educaion Air Traffic Forecas Empirical Research Based on he MCMC Mehod Jian-bo Wang,

More information

Math Week 15: Section 7.4, mass-spring systems. These are notes for Monday. There will also be course review notes for Tuesday, posted later.

Math Week 15: Section 7.4, mass-spring systems. These are notes for Monday. There will also be course review notes for Tuesday, posted later. Mah 50-004 Week 5: Secion 7.4, mass-spring sysems. These are noes for Monday. There will also be course review noes for Tuesday, posed laer. Mon Apr 3 7.4 mass-spring sysems. Announcemens: Warm up exercise:

More information

Dynamic Analysis of Damped Driven Pendulum using Laplace Transform Method

Dynamic Analysis of Damped Driven Pendulum using Laplace Transform Method , ISSN 0974-570X (Online), ISSN 0974-578 (Prin), Vol. 6; Issue No. 3; Year 05, Copyrigh 05 by CESER PUBLICATIONS Dynamic Analysis of Damped Driven Pendulum using Laplace Transform Mehod M.C. Agarana and

More information

V AK (t) I T (t) I TRM. V AK( full area) (t) t t 1 Axial turn-on. Switching losses for Phase Control and Bi- Directionally Controlled Thyristors

V AK (t) I T (t) I TRM. V AK( full area) (t) t t 1 Axial turn-on. Switching losses for Phase Control and Bi- Directionally Controlled Thyristors Applicaion Noe Swiching losses for Phase Conrol and Bi- Direcionally Conrolled Thyrisors V AK () I T () Causing W on I TRM V AK( full area) () 1 Axial urn-on Plasma spread 2 Swiching losses for Phase Conrol

More information

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion

More information

Acceleration. Part I. Uniformly Accelerated Motion: Kinematics & Geometry

Acceleration. Part I. Uniformly Accelerated Motion: Kinematics & Geometry Acceleraion Team: Par I. Uniformly Acceleraed Moion: Kinemaics & Geomery Acceleraion is he rae of change of velociy wih respec o ime: a dv/d. In his experimen, you will sudy a very imporan class of moion

More information

THE EVALUATION OF RISK-TARGETED SAFETY FACTOR AND BEHAVIOUR FACTOR FOR SELECTED STEEL STRUCTURES

THE EVALUATION OF RISK-TARGETED SAFETY FACTOR AND BEHAVIOUR FACTOR FOR SELECTED STEEL STRUCTURES THE EVALUATION OF RISK-TARGETED SAFETY FACTOR AND BEHAVIOUR FACTOR FOR SELECTED STEEL STRUCTURES Francesca CELANO 1, Jure ŽIŽMOND 2, Majaž DOLŠEK 3 ABSTRACT Seel srucures can ofen be essenial componens

More information

Failure of the work-hamiltonian connection for free energy calculations. Abstract

Failure of the work-hamiltonian connection for free energy calculations. Abstract Failure of he work-hamilonian connecion for free energy calculaions Jose M. G. Vilar 1 and J. Miguel Rubi 1 Compuaional Biology Program, Memorial Sloan-Keering Cancer Cener, 175 York Avenue, New York,

More information

Keywords: thermal stress; thermal fatigue; inverse analysis; heat conduction; regularization

Keywords: thermal stress; thermal fatigue; inverse analysis; heat conduction; regularization Proceedings Inverse Analysis for Esimaing Temperaure and Residual Sress Disribuions in a Pipe from Ouer Surface Temperaure Measuremen and Is Regularizaion Shiro Kubo * and Shoki Taguwa Deparmen of Mechanical

More information

Shiva Akhtarian MSc Student, Department of Computer Engineering and Information Technology, Payame Noor University, Iran

Shiva Akhtarian MSc Student, Department of Computer Engineering and Information Technology, Payame Noor University, Iran Curren Trends in Technology and Science ISSN : 79-055 8hSASTech 04 Symposium on Advances in Science & Technology-Commission-IV Mashhad, Iran A New for Sofware Reliabiliy Evaluaion Based on NHPP wih Imperfec

More information

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H. ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple

More information