PERFORMANCE ASSESSMENT FOR UNREINFORCED MASONRY BUILDINGS IN LOW SEISMIC HAZARD AREAS

Size: px
Start display at page:

Download "PERFORMANCE ASSESSMENT FOR UNREINFORCED MASONRY BUILDINGS IN LOW SEISMIC HAZARD AREAS"

Transcription

1 13 th World Conference on Earthquake Engneerng Vancouver, B.C., Canada August 1-6, 4 Paper No. 49 PERFORMANCE ASSESSMENT FOR UNREINFORCED MASONRY BUILDINGS IN LOW SEISMIC HAZARD AREAS Rcardo BONETT 1, Alex H. BARBAT, Lus G. PUJADES, Sergo LAGOMARSINO 3 and Andrea PENNA 4 SUMMARY The Example dstrct of Barcelona has an mportant hstorcal, archtectural and cultural value, coverng about 75 hectares of the cty. Most of ts housngs are unrenforced masonry buldngs havng an average age of about 1 years. These buldngs are tall and they have been desgned and bult wthout any earthquake resstant consderaton. Furthermore, they show some partcular features, typcal of the constructve technques at that tme, whch have been dentfed as addtonal potental damage sources. In order to evaluate the expected sesmc performance of these buldngs, a typcal sx-story unrenforced masonry buldng was modeled. The buldng was desgned and constructed n 188 and contans detals whch are typcal of that constructve perod of the Example dstrct. The dynamc behavor was studed by means of a structural analyss procedure, whch uses macro elements to model the masonry panels. Ths model descrbes the nonlnear n-plane mechancal behavor of the panels and assesses the expected damage n masonry buldngs due to earthquakes. Monte Carlo smulaton has been used to take nto account the uncertantes n the mechancal propertes of the materals. In ths way, the mean sesmc capacty curves of the buldng and ther correspondng standard devatons have been obtaned. The sesmc demand has been consdered by usng response spectra proposed by the Cartographc Insttute of Catalona (ICC, Irzarry [1]). The results here obtaned for the sesmc performance of ths type of buldngs, make clear ther hgh vulnerablty and, therefore, t s advsable retrofttng them n order to mprove ther sesmc behavor. INTRODUCTION The emblematc zone of the central dstrct of Barcelona, Span, denomnated Example, was desgned n the mddle of the nneteenth century. Ths urban area has an mportant hstorcal, archtectural and cultural value and covers approxmately 75 hectares of the cty. The most representatve typology of ths dstrct corresponds to unrenforced masonry buldngs (URM), whch are ncorporated nto numerous 1 Professor - Natonal Unversty of Colomba, Medelln Colomba Professor - Techncal Unversty of Catalona, Barcelona Span 3 Professor Unversty of Genoa, Genoa Italy 4 Researcher - European Center for Tranng and Research n Earthquake Engneerng, Pava Italy

2 almost square blocks, denomnated slands. The constructon of these buldngs took place between 186 and 194, wth 5 buldngs n average for each block. They were desgned only to vertcal statc loads, wthout any consderaton of sesmc desgn crtera, because they were bult pror to the frst Spansh sesmc code. All of the exstng URM buldngs n ths area already have exceeded ther lfe perod and only a small part of them are new renforced concrete buldngs. The slabs of these buldngs are wooden, or are made of renforced concrete or steel (accordng to the buldng perod) wth ceramc celng vaults. Due to the great heght of the frst floor of these buldngs, almost all of them have soft frst floors. Moreover, due to the need of bgger commercal areas n the frst floors, cast ron columns were used nstead of masonry walls, reducng even more the stffness of the buldngs. Therefore, we expected hgh vulnerablty for ths buldng typology. Recently, researchers at the ICC, have re-evaluated the sesmc hazard of Barcelona, approachng the problem from two ponts of vew: determnstc and probablstc (Irzarry [1]). So, two types of response spectra are avalable: the frst one corresponds to the bggest hstorcal earthquake n the cty (determnst case) and the second s the 475 year return perod earthquake, namely, the earthquake whose ntensty has a 1% probablty to be exceeded n a 5 years perod (probablstc case). Obtanng these two elastc response spectra has been an mportant contrbuton to the defnton of the sesmc hazard, because we are now able to apply capacty-demand based analyses to evaluate the sesmc performance of the buldngs n the cty. The N method proposed by Fajfar [] s used wth ths am. Startng from the obtaned demand of spectral dsplacement and usng the damage states proposed by Calv [3] for URM structures, the expected damage grade and thus the performance of the buldng can be easly determned. For many years n the past, the mechancal propertes of the materals used to buld the structures of the Example, were determned emprcally. Therefore these propertes may show a wde range of varablty and a hgh uncertanty. In order to keep these uncertantes wthn a reasonable range n the case of our buldng class, we have requested expert opnons and we have used Monte Carlo smulatons to evaluate probablstc capacty spectra. As a result, the man parameters defnng the mechancal characterstcs of the model are defned by random varables whch, startng from smple assumptons about ther probablty densty functon, can be characterzed by a mean value and ts covarance. A number of buldng samples are then generated n such a way that the numercal values for the parameters and propertes nvolved n the model ft well the correspondng probablty densty functon. Ths smulaton process allows descrbng the behavor of a wde group of buldngs showng smlar geometrcal and constructve features. Furthermore, the Monte Carlo technques allow studyng the nfluence of the uncertantes n the structural parameters on the evaluaton of the sesmc performance level. STRUCTURAL TYPOLOGY The typcal Example URM buldng, whch has been studed, has sx stores, brck walls of 3 cm for the façade walls and 15 cm n the other walls. The two frst floors have metallc beams and ceramc celng vaults smply supported on metallc man beams and cast ron columns. Rubble s placed on the upper part of the vaults and above t there s a lme mortar layer and the pavement (see Fgure 1). For the other stores, the slab s made of wooden beams, supportng the ceramc celng vaults, as t can be seen n Fgure. At the basement and ground floors, the masonry bearng walls of the upper part of the structure are supported on metallc man beams, whch n turn, are supported on cast ron columns. The columns are supported on a block, whch s supported on the masonry foundaton, beng ths type of connecton very deformable.

3 Fgure 1. Detal of the slab wth steel beam and ceramcs celng vaults (taken from the orgnal archtectonc plans of the buldng). Fgure. Detal of the slab wth wood beam and ceramcs celng vaults (taken from the orgnal archtectonc plans of the buldng). The Example dstrct has about 9 housngs and about 7% of them correspond to the URM typology here descrbed. Therefore, ths sx-story URM buldng constructed n 188, wth detals typcal of that constructve perod n the Example dstrct, has been chosen for a detaled study of ths type of buldngs. The man purpose of ths work has been evaluatng the dynamc behavor and sesmc performance of the buldngs of Barcelona whch are well represented by ths typology. The dstrbuton n plant of the buldng s almost rectangular (18.9 m 4.5 m) and the buldng has a central and two lateral squared patos. In elevaton, the buldng shows certan rregulartes, such as: cast ron columns at the ground floor, masonry bearng walls drectly supported on metallc man beams, whch, n turn, are supported on the mentoned columns. Therefore, there s a consderable varaton of the stffness wth the heght of the structure, reducng ts sesmc capacty n such a way that we may expect the typcal collapse mechansm produced by the presence of a soft floor. SEISMIC DEMAND Barcelona, cty located n the northeast of Span, has a moderate sesmc hazard and low tectonc actvty. Startng from 1998 a detaled analyss of the mcrozonaton of the cty has been undertaken. Ths analyss allowed classfyng the sol of the cty n four types correspondng to 4 homogeneous areas (Cd [4]). The

4 sesmc hazard, consderng the sze of the acton n terms of the ntensty and spectral acceleratons for perods of,.3,.6, 1. and. s, has been recently reevaluated. The problem has been analyzed startng both from the determnstc and probablstc ponts of vew. Fnally, the sesmc demand was defned by means of the elastc response spectra for the four zones of the cty (Irzarry [1]). The Acceleraton Dsplacement Response Spectra (ADRS), correspondng to the determnstc and probablstc hazard scenaros, can be seen n Fgure 3..6 S a (g).5 ICC probablstc case -spectrum.4.3. ICC determnstc case - spectrum S d (cm) Fgure 3. Response spectra for determnstc and probablstc scenaros. STRUCTURAL CAPACITY Buldng Structural model TREMURI program has been used for modelng the consdered buldng type. Ths program was developed by Galasco [5]. A non-lnear macro-element model, able to reproduce earthquake damage to masonry buldngs and falure modes observed n expermental testng, s mplemented n the program: t allows 3-dmensonal modellng and several sesmc analyss procedures. The 3-dmensonal modellng of whole URM buldngs starts from some hypotheses on ther structural and sesmc behavour: the bearng structure, both referrng to vertcal and horzontal loads, s dentfed, nsde the constructon, wth walls and floors (or vaults); the walls are the bearng elements, whle the floors, n addton to share vertcal loads to the walls, are consdered as planar stffenng elements (orthotropc 3-4 nodes membrane elements), on whch the horzontal actons dstrbuton between the walls depend; the local flexural behavour of the floors and the walls out-of-plane response are not computed because they are consdered neglgble wth respect to the global buldng response, whch s governed by ther n-plane behavour (a global sesmc response s possble only f vertcal and horzontal elements are properly connected). A frame-type representaton of the n-plane behavour of masonry walls s adopted: each wall of the buldng s subdvded nto pers and lntels ( nodes macro-elements) connected by rgd areas (nodes). Earthquake damage observaton shows, n fact, that only rarely (very rregular geometry or very small openngs) cracks appear n these areas of the wall: for ths reason these regons deformaton s assumed to be neglgble, relatvely to the macro-elements non-lnear deformatons governng the sesmc response. The presence of strngcourses (beam elements), te-rods (noncompressve spar elements), prevous damage, heterogeneous masonry portons, gaps and rregulartes can be easly ncluded n the structural model.

5 The non-lnear macro-element model, representatve of a whole masonry panel, proposed by Gambarotta [7], permts, wth a lmted number of degrees of freedom (8), to represent the two man masonry falure modes, bendng-rockng and shear-sldng (wth frcton) mechansms, on the bass of mechancal assumptons. Ths model consders, by means of nternal varables, the shear-sldng damage evoluton, whch controls the strength deteroraton (softenng) and the stffness degradaton. Fgure 4 shows the three substructures, whch a macro element s dvded: Two layers, nferor1 and superor3, n whch s concentrated the bendng and axal effects. Fnally a central part, ths one suffers shear deformatons and presents no evdence of axal or bendng deformatons. A complete cnematc model should take nto account the three degrees of freedom for each node and j on the extremtes: axal dsplacement w, horzontal dsplacement u and rotaton ϕ. There are two degrees of freedom for the central zone: axal dsplacement δ and rotaton φ (Fgure 4). h 3 ϕ j w j u j j ϕ w δ φ u M N T n m M j 3 N j j T j T M N (a) 1 ϕ 1 1 w b w 1 u 1 u ϕ s 1 T 1 M 1 N 1 (b) M 1 1 T 1 1 N 1 T M N Fgure 4. Cnematc model for the macro element [6]. Thus, cnematc s descrbed by an eght degree freedom vector, a T = {u w ϕ u j w j ϕ j δ φ}, whch s obtaned for each macro element. It s assumed, for ths hypothess, that the extremtes have an nfntesmal wdth ( ). The overturnng mechansm, whch happens because the materal does not resst tracton stress s modeled by a mono lateral elastc contact between 1 and 3 nterfaces. The consttutve equatons between the cnematc varables w, ϕ and the correspondent statc quanttes n and m are uncoupled to the lmt m b condton, when the secton s smaller than the entre compresson zone. n 6 For substructure 1 the followng equatons are obtaned: ( w ) N N = ka δ + (1) 1 M = kab ϕ φ + M 1 ( ) ()

6 Where A = s b, corresponds to the transversal secton of the panel. The nelastc contrbuton N and M are obtaned from the unlateral condton of perfect elastc contact: N 1 [ ϕ φ b + ( δ w )] H e, k A = b 8ϕ φ 6 (3) 1 [ ] H e k A (4) M = [( ϕ φ) b ( δ w )] ϕ φ b + ( δ w ) b 4( ϕ φ) ϕ φ 6 Where H () s the Heavsde s functon. The panel s shear response s expressed consderng a unform shear deformaton dstrbuton u u j γ = + φ n the central part and mposng a relatonshp between the cnematc quanttes u, h u j andφ, and the shear stress T = Tj. The crackng damage s usually located on the dagonals, where the dsplacement take place along the jonts and s represented by an nelastc deformaton component, whch s actvated when the Coulomb s lmt frcton condton s reached. From the effectve shear deformaton correspondng to module and ndcatng the elastc shear module as G, the consttutve equatons can be expressed as: ( u u + h) T GA T = j φ + h (5) GA cα h (6) T = u u j + φh + f h 1+ cα GA Where the nelastc component T ncludes the frcton stress f effect, opposed to the sldng mechansm, and nvolves a damage parameter α and an un-dmensonal coeffcent whch controls the nelastc deformaton c. In ths model, the frcton plays the role of an ntern varable defned by the followng lmt condton [6]: Φ S = f µ N (7) Where µ corresponds to the frcton coeffcent. These consttutve equatons can represent the panel s resstance varaton due to changes on axal stresses N j = N. The damage and ts effects upon panel s mechancal characterstcs are descrbed by the damage varable α whch grows accordng to falure crtera [5]: Φ d = Y ( S ) R( α ), (8) Where Y = 1 cq s the rate of energy lberaton by damage; R s the resstance functon and T S = { t n m } s the nternal stress vector. Assumng R as a growng functon of α to the crtcal value α = 1and decreasng for hgher values; the model can represent the stffness degradaton, the resstance C degradaton and pnchng effect.

7 The complete consttutve model, for the macro element, can be expressed n the followng fnte form: { T N M T N M N M } = Ka Q (9) Q + Q = j j j contans the nonlnear terms evaluated by the evoluton equatons for the damage varable α and the frcton, f. Fnally K s the elastc stffness matrx: GA / h GA / h K = GA ka ka kab kab /1 /1 GA / h GA / h GA ka ka kab kab /1 /1 ka ka ka GA kab /1 GA kab /1 GAh + kab / 6 (1) The nonlnear terms N and M are defned though the followng equaton: N = N N ; M = M M + T h (11) j The macro element shear model s a smplfcaton of a more complex contnuous model (see Gambarotta [7]) whose parameters are drectly correlated wth the masonry elements mechancals propertes. The macro model parameters should be consdered as a representatve of an average behavor. In addton to ts geometrcal characterstcs, the macro element s defned from sx parameters: The shear module G, the axal stffness K, the shear resstance of the masonry f vq, the un-dmensonal coeffcent that controls the nelastc deformaton c, the global frcton coeffcent f and β factor whch controls the softenng. The last factor s defned by pllar, as well as lntels. The macro-element used n the program to assemble the wall model keeps also nto account the effect (especally n bendng-rockng mechansms) of the lmted compressve strength of masonry (Penna [13]). Toe crushng effect s modelled by means of phenomenologcal non-lnear consttutve law wth stffness degrade n compresson: the effect of ths modellzaton on the cyclc vertcal dsplacement-rotaton nteracton s represented n Fgure 5. In order to perform non-lnear sesmc analyses of URM buldngs a set of analyss procedures has been mplemented: ncremental statc (Newton-Raphson) wth force or dsplacement control, 3D pushover analyss wth fxed load pattern and 3D tme-hstory dynamc analyss (Newmark ntegraton method; Raylegh vscous dampng). The pushover procedure, wth an effectve algorthm, transforms the problem of pushng a structure mantanng constant ratos between the appled forces nto an equvalent ncremental statc analyss wth one d.o.f. dsplacement control. Addtonal nformaton and further descrptons of the non-lnear macro-element modellng and analyss of URM buldngs can be found n Galasco [14] j

8 w ϕ Fgure 5. Cyclc vertcal dsplacement-rotaton nteracton wth (red lne) and w/o toe crushng (blue dots) [13] Macro element model for the studed buldng of the Example Fgure 6 shows a three-dmensonal vew and n plant of the model used for the representatve buldng of the Example. The model s defned by 8 walls n the x drecton (walls M1 to M8) and 6 walls n the y drecton (walls M9 to M14). Each wall has been modeled as an assemblage of pers, lntels and frame elements (n some cases) connected to the nodes of the model by means of rgd jonts. All the nodes have 5 degrees of freedom (3 dsplacement components and rotaton components correspondng to the axes x and y) except the base nodes of the model. The slabs have been modeled as an orthotropc fnte element daphragm, defned by 3 or 4 nodes connected to the three-dmensonal nodes of each level. A man analyss drecton s dentfed, whch s characterzed by a Young s modulus E 1 and the drecton perpendcular to ths one s characterzed by a Young s modulus E. Fgure 7 shows the macro element model correspondng to walls 1 and. M8 M7 M6 M14 M1 3 M1 M5 M4 M11 M1 M9 M3 M M1 Fgure 6. Three-dmensonal model of the analyzed typcal buldng of the Example. In order to analyze the constructve system of the URM buldngs of the Example, t s necessary to have a good knowledge on the materals used for ther man elements. Brcks are the basc materal of these buldngs, beng used wdely n walls, stars and slabs. The typcal dmensons of the used brcks are 3 15 cm and ther thckness vares between 3 and 11 cm. Ths knd of man-made brcks were used untl the begnnng of the XX th century. Later, mechancal systems were used, consderably mprovng ther qualty and compactness. Lme mortar was used n the constructve process of the buldngs of the

9 Example. The wde use of ths materal s assocated to constructve tradton, to consumpton habts and, apparently, to ts strength whch was consdered to be adequate at that perod. N n5 417 n6 418 n7 419 n8 4 N N n1 41 n 413 n3 414 n4 415 N N67 46 n17 47 n18 48 n19 49 n 41 N N41 41 n13 4 n14 43 n15 44 n16 45 N N15 n9 n1 n11 n1 N N n5 39 n6 393 n7 394 n8 395 N19 N n N3 438 N n48 44 N N n45 43 N N n N N69 46 n43 47 N7 N n44 43 N N43 41 n41 4 N44 N n4 45 N N17 n37 n38 N18 N19 n39 n 4 N N191 n33 n34 N19 N193 n35 n 36 N N163 n1 n n3 n4 N N165 n9 n3 N166 N167 n31 n 3 N168 Fgure 7. Macro-element model. Walls 1 and. As sad before, n ths work, probablty densty functons, pdf, are used to defne the most mportant parameters of the model. These functons are characterzed by a mean value and a covarance. The defnton of the mean value of each parameter has been defned usng the opnon of experts, who provded suffcent nformaton for defnng a model. Nevertheless, due to the subjectve character of ths nformaton, the man parameters have been consdered as random varables wth ther uncertantes. The most mportant mechancal propertes of the materals used n the analyss of the buldng of the Example are descrbed below. Masonry Young s modulus of the wall E = N/m Shear modulus G = N/m Shear strength τ = N/m Softenng factor for the pers β p =.5 Softenng factor for the lntels β d =.5 Cast ron columns Young s modulus E s = N/m Specfc weght γ s = 785 kg/m 3 Concrete columns Young s modulus E h = N/m Specfc weght γ h = 5 kg/m 3 Slabs Young s modulus n the man drecton E 1 = N/m Young s modulus n the orthogonal drecton E = N/m Shear modulus G = N/m Among all these characterstcs, those shown n Table 1 have been defned as random varables because they have an mportant nfluence on the structural response of ths type of buldngs. The normal probablty dstrbuton functon has been used for the three varables, where the mean value of each

10 parameter corresponds to the values proposed by experts. The covarance has been defned n such a way to cover a reasonable varaton range for each parameter. Table 1. Probablty dstrbuton functons for random varables. Mean value and covarance. Parameter fdp Mean Covarance Young s modulus E Normal.11 9 N/m.3 Shear strength τ Normal N/m.3 Softenng factor β p Normal.5.3 Capacty curve The capacty curve generally corresponds to the frst mode of vbraton of the structure, based on the assumpton that the fundamental mode of vbraton contans the predomnant response of the structure. In the case of the analyzed buldng, a force dstrbuton was establshed correspondng to the bendng modal shape orented along the y axs. Therefore, the walls 9, 1, 11, 1, 13 and 14 are nvolved n the analyss (see Fgure 7). However, for the sake of smplcty, the loads only wll be appled to the walls 9 and 14, whch really provde the greater stffness n that drecton. Fgure 8. Mean, mean + 1σ and mean 1σ capacty spectra. The capacty curve s obtaned by performng a pushover analyss wth ths load pattern. Ths curve descrbes the relatonshp between base shear and the roof dsplacement of an equvalent sngle degree of freedom model, characterzed by the perod and the modal mass of the thrd mode of vbraton. The response of the model of the typcal URM buldng s defned by means of the capacty curves obtaned by means of the Monte Carlo smulaton technque. Thus, 1 samples for each varable were generated and a structural model was defned for each sample group. One hundred capacty curves were thus obtaned. The advanced computatonal tool STAC [8] has been used n the smulaton process. Fgure 8 shows the mean capacty spectra together wth ther standard devatons. Ths type of the representaton shows the senstvty of these methods to the uncertantes n the structural parameters. The blnear representaton s obtaned for these three spectra usng the values of the spectral dsplacement and acceleraton for the yeldng pont ( D y, S ) ay and for the pont of the ultmate capacty ( D u ), Sau. Table shows these values for the mean capacty spectrum.

11 Table. Blnear representaton parameters of the capacty spectrum. Capacty Spectrum D y ( cm) ( g) D ( cm ) u ( g) S ay S au x DAMAGE STATE LIMITS In order to obtan the damage state lmts or the performance levels of the URM buldng of the Example, there are nether laboratory tests nor avalable values calbrated from observed damage durng earthquakes. Addtonally, the values of the mechancal propertes of the materals used n ths structural typology are not completely known. Takng nto account all these aspects, the thresholds of the spectral dsplacement for the dscrete damage states are defned based on the blnear representaton of the capacty spectrum. Table 3 shows the expressons proposed by Lagomarsno [9] to defne the varaton ntervals of the spectral dsplacement for the fve damage states here consdered: no damage, slght, moderate, severe and complete. Table 3. Spectral dsplacement for the damage states [9, 1]. Damage state Spectral dsplacement, S d No damage d <.7D Slght.7D y Moderate Dy Extensve D.5( D D ) y u y S y < d S D y S D +.5( D D ) < d + < S d Complete S d > y u y Startng from the expressons of Table 3 and usng the values of D y and D u obtaned for the sx-story buldng (see Table 4), the thresholds of the spectral dsplacement are obtaned for the fve damage states. D u D u Table 4. Thresholds for the spectral dsplacement. Damage state Threshold, S d (cm) No damage Sd <.48 Slght.48 < S d.69 Moderate.69 < S d 1.17 Extensve 1.17 < S d.61 Complete Sd >.61 SEISMIC PERFORMANCE In order to evaluate the sesmc performance of the typcal URM buldng of the Example, the N method proposed by Fajfar [] was used. Startng from a frst verson, publshed n 1987, the method has been

12 revsed and updated to the present verson, n whch the Acceleraton-Dsplacement format s used. Nowadays, the method combnes the vsual representaton advantages of the capacty spectrum method [1] wth the physcal bass of the nelastc demand spectrum [11]. The basc characterstcs of the method are: use of two dfferent mathematcal models, applcaton of the response spectrum, nonlnear statc analyss (pushover analyss) and the selecton of a model, whch takes nto account the cumulatve damage. Ths last aspect s very mportant for exstng buldngs, whch frequently have not been desgned to resst many hysteretc cycles wthn nelastc ranges []. In ths case, two response spectra (one determnstc and one probablstc) have been used to descrbe the sesmc demand. For each of them, the spectral dsplacement demand s obtaned and the performance pont s evaluated (see Table 5). Fgures 9 and 1 show the graphcal representatons of the performance pont correspondng to the determnstc and probablstc cases of the sesmc demand, respectvely. Table 5. Damage state and performance levels. Demand spectrum S d (cm) Damage state Performance levels Determnstc.67 Slght Operatonal Probablstc 1.13 Moderate Lfe-safe.4 S a (g).3 ICC determnstc case spectrum - PGA =.141 g Blnear representaton. Perfomance pont S dp =.67 cm S ap =.1 g S d (cm) Fgure 9. Sesmc performance pont (determnstc case). FRAGILITY CURVES Fraglty curves have been generated startng from the assumpton that the cumulatve probablty of reachng or exceedng a partcular damage state follows a lognormal dstrbuton. Therefore, for a gven spectral dsplacement and damage state, ths probablty can be obtaned by means of the followng equaton:

13 (1) 1 [ ] S d P DS DS / Sd = Φ ln β DS Sd, DS S d, DS s the mean value of the spectral dsplacement at whch the buldng reaches the damage state threshold DS, β ED s the standard devaton of the natural logarthm of ths spectral dsplacement and Φ s the cumulatve standard normal dstrbuton functon. The subscrpt stays for the damage state: slght ( = 1), moderate ( = ), extensve ( = 3) and complete ( = 4). In order to calculate the probabltes startng from the dstrbuton functon Φ[ ] (equaton 1), t s necessary to defne S, and for each damage state. Table 6 and Fgure 11 show the parameters and the fraglty curves obtaned for the studed URM buldng. d ED β ED.5 S a (g).4.3. ICC probablstc case spectrum - PGA =.194 g Blneal representaton Demand spectrum - Ductlty = 1.7 Performance pont S dp = 1.13 cm S ap =.14 g S d (m) Fgure 1. Sesmc performance pont (probablstc case). In order to estmate the expected damage for each sesmc hazard scenaro (determnstc and probablstc), we use the spectral dsplacements n Table 5 and the fraglty curves n Fgure 11 to obtan the probabltes of each damage state. We can see n Fgure 1 how the Slght damage state s the most probable damage state n the determnstc case, whle t s the Severe damage state n the probablstc case. Table 6. Parameters of the lognormal dstrbuton functon. Damage S d, DS β DS state Slght Moderate Extensve Complete.61.65

14 FD = Prob (DS > ds / Sd = Sd ) NO DAMAGE SLIGHT MODERATE EXTENSIVE COMPLETE S d (cm) Fgure 11. Fraglty curves for sx-story URM buldng of the Example. Fgure 1. Damage probabltes for the determnstc and probablstc sesmc scenaros. DISCUSSION AND CONCLUSIONS The sesmc performance of a typcal unrenforced masonry buldng of the Example dstrct n Barcelona, Span, has been analyzed. The capacty of the buldng was studed by usng a structural model, whch uses macro elements for the masonry panels. The expected demand has been defned by two response spectra proposed by the Cartographc Insttute of Catalona. The frst one corresponds to the bggest hstorcal earthquake n the cty (determnstc case) whle the second corresponds to a 475 years return perod earthquake (probablstc case). The mechancal propertes of the materals used for the constructon of the URM buldngs n Barcelona, show a hgh varablty and Monte Carlo smulaton has been performed to take nto account the uncertantes. In ths way we have obtaned mean sesmc capacty curves together wth ther correspondng standard devatons. The results show an mportant dsperson, whch can also be observed n the expected damage. The performance pont of the URM buldng of the Example for the determnstc case remans wthn the elastc range and the most probable damage state s the slght. Nevertheless, when the probablstc case s analyzed, the most probable damage state s the

15 severe or pre-collapse state. Ths stuaton s typcal of areas wth low to moderate sesmc hazard. Any way, n both cases, probablstc and determnstc, we found sgnfcant probabltes for the severe and collapse damage states, ndcatng the hgh vulnerablty of most of the buldngs n the Example dstrct. Ths hgh vulnerablty and expected damage s due to the neglect of any sesmc consderaton n the cty. Ths fact s ncreased because the low sesmc requrements planned n the Spansh sesmc codes for the cty and would result n consderable damage n the case of a relatvely low earthquake. Therefore, an mportant concluson of ths work s that t s very convenent to serously consder retrofttng and upgradng the sesmc performance of the buldngs of the cty, partcularly those whose functon s mportant n the post-earthquake emergency, as for example, hosptals. ACKNOWLEDGMENTS Ths work has been partally sponsored by the Spansh Mnstry of Scence and Technology and wth FEDER funds (projects: REN--174-C5-1/RIES, REN C4-1 y REN-3365/- RIES), by the European Commsson (RISK-UE Project, contract EVK4-CT--14) and by the Cvl Engneerng School of Barcelona (UPC). REFERENCES 1. Irzarry J., Goula X., Roca A. and T. Susagna. Earthquake rsk scenaros for monuments of Barcelona. 13 th World Conference on Earthquake Engneerng. Vancouver B.C. 4. Paper N pp.. Fajfar, P. and Gaspersc, P. The N method for the sesmc damage analyss of RC buldngs. Earthquake Engneerng and Structural Dynamcs 1996, 5: Calv, M. C. A dsplacement-based approach for vulnerablty evaluaton of classes of buldngs, Journal of Earthquake Engneerng 1999, 3 (3): Cd, J. Zonfcacón sísmca de la cudad de Barcelona basada en métodos de smulacón numérca de efectos locales. Tess Doctoral, Unversdad Poltécnca de Cataluña, Barcelona, Galasco, M., Lagomarsno, S. y Penna, A. TREMURI Program: Sesmc Analyser of 3D Masonry Buldngs, Unversdad de Genoa,. 6. Brencch, A. and Lagomarsno, S. A macroelement dynamc model for masonry shear walls, Computer Methods n Structural Masonry-4, edted by G.N.Pande, J. Mddleton and B. Kralj, Gambarotta, L. and Lagomarsno, S. A mcrocrack damage model for brttle materals. Internatonal Journal Solds and Structures 1993, 3: STAC program. Stochastc análss computatonal, Internatonal Center for Numercal Methods n Engneerng (CIMNE), Barcelona,. 9. Lagomarsno, S. and Penna, A. Gudelnes for the mplementaton of the II level vulnerablty methodology. WP4: Vulnerablty assessment of current buldngs. RISK-UE project: An advanced approach to earthquake rsk scenaros wth applcaton to dfferent European towns, Freeman, S. A. Development and use of capacty spectrum method. Proceedngs of the 6 th U. S. Natonal Conference on Earthquake Engneerng, Seattle, EERI. Oakland. 11. Fajfar, P. Capacty spectrum method based on nelastc demand spectra. Earthquake Engneerng and Structural Dynamcs 1999, 8: RISK-UE project. An advanced approach to earthquake rsk scenaros wth applcaton to dfferent European towns, European Commsson, contract EVK4-CT--14, Penna, A A macro-element procedure for the dynamc non-lnear analyss of masonry buldngs, Ph.D. Dssertaton,, Poltecnco d Mlano (n talan). 14. Galasco, A., Lagomarsno S., Penna, A., Resemn, S. Non-lnear sesmc analyss of masonry structures, Proc. 13th WCEE, 4, Vancouver, paper Nº. 843.

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

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

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES Manuel J. C. Mnhoto Polytechnc Insttute of Bragança, Bragança, Portugal E-mal: mnhoto@pb.pt Paulo A. A. Perera and Jorge

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

INCREMENTAL DYNAMIC ANALYSIS AND PUSHOVER ANALYSIS. A PROBABILISTIC COMPARISON

INCREMENTAL DYNAMIC ANALYSIS AND PUSHOVER ANALYSIS. A PROBABILISTIC COMPARISON COMPDYN 2 III ECCOMAS Thematc Conference on Computatonal Metho n Structural Dynamcs and Earthquake Engneerng M. Papadrakaks, M. Fragadaks, V. Plevrs (e.) Corfu, Greece, 25 28 May 2 INCREMENTAL DYNAMIC

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

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

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

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

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

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 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

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

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

NUMERICAL RESULTS QUALITY IN DEPENDENCE ON ABAQUS PLANE STRESS ELEMENTS TYPE IN BIG DISPLACEMENTS COMPRESSION TEST

NUMERICAL RESULTS QUALITY IN DEPENDENCE ON ABAQUS PLANE STRESS ELEMENTS TYPE IN BIG DISPLACEMENTS COMPRESSION TEST Appled Computer Scence, vol. 13, no. 4, pp. 56 64 do: 10.23743/acs-2017-29 Submtted: 2017-10-30 Revsed: 2017-11-15 Accepted: 2017-12-06 Abaqus Fnte Elements, Plane Stress, Orthotropc Materal Bartosz KAWECKI

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

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

RELIABILITY ASSESSMENT

RELIABILITY ASSESSMENT CHAPTER Rsk Analyss n Engneerng and Economcs RELIABILITY ASSESSMENT A. J. Clark School of Engneerng Department of Cvl and Envronmental Engneerng 4a CHAPMAN HALL/CRC Rsk Analyss for Engneerng Department

More information

Lecture 8 Modal Analysis

Lecture 8 Modal Analysis Lecture 8 Modal Analyss 16.0 Release Introducton to ANSYS Mechancal 1 2015 ANSYS, Inc. February 27, 2015 Chapter Overvew In ths chapter free vbraton as well as pre-stressed vbraton analyses n Mechancal

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

Modal Strain Energy Decomposition Method for Damage Detection of an Offshore Structure Using Modal Testing Information

Modal Strain Energy Decomposition Method for Damage Detection of an Offshore Structure Using Modal Testing Information Thrd Chnese-German Jont Symposum on Coastal and Ocean Engneerng Natonal Cheng Kung Unversty, Tanan November 8-16, 2006 Modal Stran Energy Decomposton Method for Damage Detecton of an Offshore Structure

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

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth

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

Lifetime prediction of EP and NBR rubber seal by thermos-viscoelastic model

Lifetime prediction of EP and NBR rubber seal by thermos-viscoelastic model ECCMR, Prague, Czech Republc; September 3 th, 2015 Lfetme predcton of EP and NBR rubber seal by thermos-vscoelastc model Kotaro KOBAYASHI, Takahro ISOZAKI, Akhro MATSUDA Unversty of Tsukuba, Japan Yoshnobu

More information

Chapter 13: Multiple Regression

Chapter 13: Multiple Regression Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to

More information

A Robust Method for Calculating the Correlation Coefficient

A Robust Method for Calculating the Correlation Coefficient A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal

More information

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

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

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

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

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

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

COMPONENT-LEVEL AND SYSTEM-LEVEL SENSITIVITY STUDY FOR EARTHQUAKE LOSS ESTIMATION

COMPONENT-LEVEL AND SYSTEM-LEVEL SENSITIVITY STUDY FOR EARTHQUAKE LOSS ESTIMATION 13 th World Conference on Earthquae Engneerng Vancouver, B.C., Canada August 1-6, 2004 Paper No. 1070 COMPONEN-LEVEL AND SYSEM-LEVEL SENSIIVIY SUDY FOR EARHQUAKE LOSS ESIMAION Hesameddn ASLANI 1, Eduardo

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

Composite Hypotheses testing

Composite Hypotheses testing Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter

More information

Week 9 Chapter 10 Section 1-5

Week 9 Chapter 10 Section 1-5 Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,

More information

Constitutive Modelling of Superplastic AA-5083

Constitutive Modelling of Superplastic AA-5083 TECHNISCHE MECHANIK, 3, -5, (01, 1-6 submtted: September 19, 011 Consttutve Modellng of Superplastc AA-5083 G. Gulano In ths study a fast procedure for determnng the constants of superplastc 5083 Al alloy

More information

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed (2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected

More information

Statistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA

Statistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA 14 th Internatonal Users Conference Sesson: ALE-FSI Statstcal Energy Analyss for Hgh Frequency Acoustc Analyss wth Zhe Cu 1, Yun Huang 1, Mhamed Soul 2, Tayeb Zeguar 3 1 Lvermore Software Technology Corporaton

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

One Dimensional Axial Deformations

One Dimensional Axial Deformations One Dmensonal al Deformatons In ths secton, a specfc smple geometr s consdered, that of a long and thn straght component loaded n such a wa that t deforms n the aal drecton onl. The -as s taken as the

More information

Operating conditions of a mine fan under conditions of variable resistance

Operating conditions of a mine fan under conditions of variable resistance Paper No. 11 ISMS 216 Operatng condtons of a mne fan under condtons of varable resstance Zhang Ynghua a, Chen L a, b, Huang Zhan a, *, Gao Yukun a a State Key Laboratory of Hgh-Effcent Mnng and Safety

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

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method

More information

THE EFFECT OF BEAM TO COLUMN CONNECTION IN ARC PORTAL FRAME

THE EFFECT OF BEAM TO COLUMN CONNECTION IN ARC PORTAL FRAME THE EFFECT OF BEAM TO COLUMN CONNECTON N ARC PORTAL FRAME Asko Keronen Rakenteden Mekankka, Vol. 26 No 2 1993, ss. 35-5 SUMMARY A full scale rc (renforced concrete) portal frame has been bult n order to

More information

Appendix B: Resampling Algorithms

Appendix B: Resampling Algorithms 407 Appendx B: Resamplng Algorthms A common problem of all partcle flters s the degeneracy of weghts, whch conssts of the unbounded ncrease of the varance of the mportance weghts ω [ ] of the partcles

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

STUDY OF A THREE-AXIS PIEZORESISTIVE ACCELEROMETER WITH UNIFORM AXIAL SENSITIVITIES

STUDY OF A THREE-AXIS PIEZORESISTIVE ACCELEROMETER WITH UNIFORM AXIAL SENSITIVITIES STUDY OF A THREE-AXIS PIEZORESISTIVE ACCELEROMETER WITH UNIFORM AXIAL SENSITIVITIES Abdelkader Benchou, PhD Canddate Nasreddne Benmoussa, PhD Kherreddne Ghaffour, PhD Unversty of Tlemcen/Unt of Materals

More information

Turbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH

Turbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH Turbulence classfcaton of load data by the frequency and severty of wnd gusts Introducton Oscar Moñux, DEWI GmbH Kevn Blebler, DEWI GmbH Durng the wnd turbne developng process, one of the most mportant

More information

Here is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y)

Here is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y) Secton 1.5 Correlaton In the prevous sectons, we looked at regresson and the value r was a measurement of how much of the varaton n y can be attrbuted to the lnear relatonshp between y and x. In ths secton,

More information

x = , so that calculated

x = , so that calculated Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to

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

LIMIT ANALYSIS IN LARGE DISPLACEMENTS OF MASONRY ARCHES SUBJECTED TO VERTICAL AND HORIZONTAL LOADS

LIMIT ANALYSIS IN LARGE DISPLACEMENTS OF MASONRY ARCHES SUBJECTED TO VERTICAL AND HORIZONTAL LOADS SAC2014 9 th Internatonal Conference on Structural Analyss of storcal Constructons F. Peña & M. Chávez (eds.) Mexco Cty, Mexco, 14 17 October 2014 LIMIT ANALYSIS IN LARGE DISPLACEMENTS OF MASONRY ARCES

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

Uncertainty in measurements of power and energy on power networks

Uncertainty in measurements of power and energy on power networks Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:

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

Simulation and Probability Distribution

Simulation and Probability Distribution CHAPTER Probablty, Statstcs, and Relablty for Engneers and Scentsts Second Edton PROBABILIT DISTRIBUTION FOR CONTINUOUS RANDOM VARIABLES A. J. Clark School of Engneerng Department of Cvl and Envronmental

More information

Numerical analysis of buried pipes subjected to lateral soil movements

Numerical analysis of buried pipes subjected to lateral soil movements Numercal analss of bured ppes subjected to lateral sol movements P. Casamchele, M. Mauger & E. Motta Department of Cvl and Envronmental Engneerng, Catana Unverst, Ital Abstract he response of bured ppes

More information

Effect of loading frequency on the settlement of granular layer

Effect of loading frequency on the settlement of granular layer Effect of loadng frequency on the settlement of granular layer Akko KONO Ralway Techncal Research Insttute, Japan Takash Matsushma Tsukuba Unversty, Japan ABSTRACT: Cyclc loadng tests were performed both

More information

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient 58:080 Expermental Engneerng 1 OBJECTIVE Lab 2e Thermal System Response and Effectve Heat Transfer Coeffcent Warnng: though the experment has educatonal objectves (to learn about bolng heat transfer, etc.),

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

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

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

DETERMINATION OF UNCERTAINTY ASSOCIATED WITH QUANTIZATION ERRORS USING THE BAYESIAN APPROACH

DETERMINATION OF UNCERTAINTY ASSOCIATED WITH QUANTIZATION ERRORS USING THE BAYESIAN APPROACH Proceedngs, XVII IMEKO World Congress, June 7, 3, Dubrovn, Croata Proceedngs, XVII IMEKO World Congress, June 7, 3, Dubrovn, Croata TC XVII IMEKO World Congress Metrology n the 3rd Mllennum June 7, 3,

More information

General displacement arch-cantilever element method for stress analysis of arch dam

General displacement arch-cantilever element method for stress analysis of arch dam Water Scence and Engneerng, 009, (): 3-4 do:0.388/j.ssn.674-370.009.0.004 http://kkb.hhu.edu.cn e-mal: wse@hhu.edu.cn General dsplacement arch-cantlever element method for stress analyss of arch dam Hao

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

Statistics II Final Exam 26/6/18

Statistics II Final Exam 26/6/18 Statstcs II Fnal Exam 26/6/18 Academc Year 2017/18 Solutons Exam duraton: 2 h 30 mn 1. (3 ponts) A town hall s conductng a study to determne the amount of leftover food produced by the restaurants n the

More information

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI 2017 2nd Internatonal Conference on Electrcal and Electroncs: echnques and Applcatons (EEA 2017) ISBN: 978-1-60595-416-5 Study on Actve Mcro-vbraton Isolaton System wth Lnear Motor Actuator Gong-yu PAN,

More information

Analytical simplified procedure for the evaluation of the RC buildings.

Analytical simplified procedure for the evaluation of the RC buildings. Analytcal smplfed procedure for the evaluaton of the RC buldngs. N. Gattesco, R. Franceschns, F. Zorzn, Unversty of Treste, Italy SUMMARY: Numerous exstng buldngs are located n sesmc prone areas and many

More information

/ n ) are compared. The logic is: if the two

/ n ) are compared. The logic is: if the two STAT C141, Sprng 2005 Lecture 13 Two sample tests One sample tests: examples of goodness of ft tests, where we are testng whether our data supports predctons. Two sample tests: called as tests of ndependence

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

PRATICAL STATIC CALCULATION METHOD FOR ESTIMATING ELASTO-PLASTIC DYNAMIC RESPONSES OF SPACE FRAMES

PRATICAL STATIC CALCULATION METHOD FOR ESTIMATING ELASTO-PLASTIC DYNAMIC RESPONSES OF SPACE FRAMES th World Congress on Computatonal Mechancs (WCCM XI) 5th European Conference on Computatonal Mechancs (ECCM V) 6th European Conference on Computatonal lud ynamcs (EC VI) E. Oñate, J. Olver and A. Huerta

More information

Experiment 1 Mass, volume and density

Experiment 1 Mass, volume and density Experment 1 Mass, volume and densty Purpose 1. Famlarze wth basc measurement tools such as verner calper, mcrometer, and laboratory balance. 2. Learn how to use the concepts of sgnfcant fgures, expermental

More information

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown

More information

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY. Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR

More information

Credit Card Pricing and Impact of Adverse Selection

Credit Card Pricing and Impact of Adverse Selection Credt Card Prcng and Impact of Adverse Selecton Bo Huang and Lyn C. Thomas Unversty of Southampton Contents Background Aucton model of credt card solctaton - Errors n probablty of beng Good - Errors n

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

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

Modeling of Dynamic Systems

Modeling of Dynamic Systems Modelng of Dynamc Systems Ref: Control System Engneerng Norman Nse : Chapters & 3 Chapter objectves : Revew the Laplace transform Learn how to fnd a mathematcal model, called a transfer functon Learn how

More information

PARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY

PARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY POZNAN UNIVE RSITY OF TE CHNOLOGY ACADE MIC JOURNALS No 86 Electrcal Engneerng 6 Volodymyr KONOVAL* Roman PRYTULA** PARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY Ths paper provdes a

More information

Determination of Dose Factors for External Gamma Radiation in Dwellings

Determination of Dose Factors for External Gamma Radiation in Dwellings Determnaton of Dose Factors for External Gamma daton n Dwellngs M.F. Maduar and G. Hromoto Insttuto de Pesqusas Energétcas e Nucleares, CP 049, 05422-970 São Paulo, Brazl INTRODUCTION e largest contrbuton

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

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

STATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION

STATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION STATIC ANALYSIS OF TWO-LERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION Ákos József Lengyel István Ecsed Assstant Lecturer Emertus Professor Insttute of Appled Mechancs Unversty of Mskolc Mskolc-Egyetemváros

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

risk and uncertainty assessment

risk and uncertainty assessment Optmal forecastng of atmospherc qualty n ndustral regons: rsk and uncertanty assessment Vladmr Penenko Insttute of Computatonal Mathematcs and Mathematcal Geophyscs SD RAS Goal Development of theoretcal

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

More information

Supplementary Notes for Chapter 9 Mixture Thermodynamics

Supplementary Notes for Chapter 9 Mixture Thermodynamics Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects

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

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

Statistics Chapter 4

Statistics Chapter 4 Statstcs Chapter 4 "There are three knds of les: les, damned les, and statstcs." Benjamn Dsrael, 1895 (Brtsh statesman) Gaussan Dstrbuton, 4-1 If a measurement s repeated many tmes a statstcal treatment

More information

A Particle Filter Algorithm based on Mixing of Prior probability density and UKF as Generate Importance Function

A Particle Filter Algorithm based on Mixing of Prior probability density and UKF as Generate Importance Function Advanced Scence and Technology Letters, pp.83-87 http://dx.do.org/10.14257/astl.2014.53.20 A Partcle Flter Algorthm based on Mxng of Pror probablty densty and UKF as Generate Importance Functon Lu Lu 1,1,

More information

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur Analyss of Varance and Desgn of Experment-I MODULE VII LECTURE - 3 ANALYSIS OF COVARIANCE Dr Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Any scentfc experment s performed

More information

Polynomial Regression Models

Polynomial Regression Models LINEAR REGRESSION ANALYSIS MODULE XII Lecture - 6 Polynomal Regresson Models Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Test of sgnfcance To test the sgnfcance

More information

Response Analysis of Multi-Storey RC Buildings under Equivalent Static and Dynamic Loads According to Egyptian Code

Response Analysis of Multi-Storey RC Buildings under Equivalent Static and Dynamic Loads According to Egyptian Code nternatonal Journal of Cvl and tructural Engneerng Research N 2348-7607 (Onlne) Vol. 2, ssue 1, pp: (79-88), Month: Aprl 2014 - eptember 2014, Avalable at: www.researchpublsh.com Response Analyss of Mult-torey

More information

Study on Non-Linear Dynamic Characteristic of Vehicle. Suspension Rubber Component

Study on Non-Linear Dynamic Characteristic of Vehicle. Suspension Rubber Component Study on Non-Lnear Dynamc Characterstc of Vehcle Suspenson Rubber Component Zhan Wenzhang Ln Y Sh GuobaoJln Unversty of TechnologyChangchun, Chna Wang Lgong (MDI, Chna [Abstract] The dynamc characterstc

More information

CS-433: Simulation and Modeling Modeling and Probability Review

CS-433: Simulation and Modeling Modeling and Probability Review CS-433: Smulaton and Modelng Modelng and Probablty Revew Exercse 1. (Probablty of Smple Events) Exercse 1.1 The owner of a camera shop receves a shpment of fve cameras from a camera manufacturer. Unknown

More information