CHAPTER 7 MODAL ANALYSIS

Size: px
Start display at page:

Download "CHAPTER 7 MODAL ANALYSIS"

Transcription

1 CHAPER 7 MODAL ANALYSIS 7. Itroducto I Chater we systeatcally defed the equatos of oto for a ult dof (dof) syste ad trasfored to the s doa usg the Lalace trasfor. Chater 3 dscussed frequecy resoses ad udaed ode shaes. Chater 5 dscussed the state sace for of equatos of oto wth arbtrary dag. It also covered the subject of colex odes. Heavly daed structures or structures wth exlct dag eleets, such as dashots, result colex odes ad requre state sace soluto techques usg the orgal couled equatos of oto. Lghtly daed structures are tycally aalyed wth the oral ode ethod, whch s the subject of ths chater. he ablty to th about vbratg systes ters of odal roertes s a very owerful techque that serves oe well both erforg aalyss ad uderstadg test data. he ey to oral ode aalyss s to develo tools whch allow oe to recostruct the overall resose of the syste as a suerosto of the resoses of the dfferet odes of the syste. I aalyss, the odal ethod allows oe to relace the -couled dfferetal equatos wth -ucouled equatos, where each ucouled equato reresets the oto of the syste for that ode of vbrato. If atural frequeces ad ode shaes are avalable for the syste, the t s easy to vsuale the oto of the syste each ode, whch s the frst ste beg able to uderstad how to odfy the syste to chage ts characterstcs. Suarg the odal aalyss ethod of aalyg lear echacal systes ad the beefts derved: ) Solve the udaed egevalue roble, whch detfes the resoat frequeces ad ode shaes (egevalues ad egevectors), useful theselves for uderstadg basc otos of the syste. ) Use the egevectors to ucoule or dagoale the orgal set of couled equatos, allowg the soluto of -ucouled sdof robles stead of solvg a set of -couled equatos.

2 3) Calculate the cotrbuto of each ode to the overall resose. hs also allows oe to reduce the se of the roble by elatg odes that caot be excted ad/or odes that have o oututs at the desred dof s. Also, hgh frequecy odes that have lttle cotrbuto to the syste at lower frequeces ca be elated or aroxately accouted for, further reducg the se of the syste to be aalyed. 4) Wrte the syste atrx, A, by secto. Asseble the ut ad outut atrces, B ad C, usg arorate egevector ters. Frequecy doa ad forced traset resose robles ca be solved at ths ot. If colete egevectors are avalable, tal codto traset robles ca also be solved. For lghtly daed systes, roortoal dag ca be added, whle stll allowg the equatos to be ucouled. 7. Egevalue Proble 7.. Equatos of Moto We wll start by wrtg the udaed hoogeeous (uforced) equatos of oto for the odel Fgure 7.. he we wll defe ad solve the egevalue roble. F F 3 F 3 3 Fgure 7.: Udaed tdof odel. Fro (.5) wth ad c c 0: && + 0 (7.)

3 0 0 && && && (7.) 7.. Prcal (Noral) Mode Defto Sce the syste s coservatve (t has o dag), oral odes of vbrato wll exst. Havg oral odes eas that at certa frequeces all ots the syste wll vbrate at the sae frequecy ad hase,.e., all ots the syste wll reach ther u ad axu dslaceets at the sae ot te. Havg oral odes ca be exressed as (Weaver 990): Where: j t ( ) ω +φ s ω t +φ I(e ) (7.3) vector of dslaceets for all dof s at the th frequecy the th egevector, the ode shae for the th resoat frequecy ω the th egevalue, th resoat frequecy φ a arbtrary tal hase agle For our tdof syste, for the th frequecy, the equato would aear as: s ω t+φ ( ) 3 3, (7.4) where the dces the ter rereset the th dof ad the th ode of the odal atrx Egevalues / Characterstc Equato Sce the equato of oto && + 0 (7.5)

4 ad the for of the oto ( t ) s ω +φ (7.6) are ow, ca be dfferetated twce ad substtuted to the equato of oto: ( t ) && ω s ω +φ (7.7) ( ) ( ) ω s ω t + φ + s ω t +φ 0 (7.8) Cacelg the se ters: ω + (7.9) 0 ω (7.0) Equato (7.0) s the egevalue roble ostadard for, where the stadard for s (Strag 998): A λ (7.) he soluto of the sultaeous equatos whch ae u the stadard for egevalue roble s a vector such that whe s ultled by A, the roduct s a scalar ultle of tself. he ostadard roble s ostadard because the ass atrx falls o the rght-had sde. he for of the atrx resets o roble for had calculatos, but for couter calculatos t s best trasfored to stadard for. Rewrtg the ostadard for egevalue roble as a hoogeeous equato: ( ) ω (7.) 0 A trval soluto, 0, exsts but s of o cosequece. he oly ossblty for a otrval soluto s f the deterat of the coeffcet atrx s ero (Strag 998). Exadg the atrx etres:

5 ω 0 (7.3) Perforg the atrx subtracto: ω 0 ω 0 (7.4) 0 ω Settg the deterat of the coeffcet atrx equal to ero: ω 0 ω 0 0 ω (7.5) he deterat results a olyoal ω, the characterstc equato, where the roots of the olyoal are the egevalues, oles, or resoat frequeces of the syste. ω + ω ω ( ) 3 4 ω ω + ω (7.6a,b) wo of the roots are at the org: Solvg for ω as a quadratc (7.6b) above: ω ω 0 (7.7) ( ) 4 ± ± 3

6 6, 3, (7.8) ω ± ω ± 3 3 (7.9) For each of the three egevalue ars, there exsts a egevector, whch gves the ode shae of the vbrato at that frequecy Egevectors o obta the egevectors of the syste, ay oe of the degrees of freedo, say, s selected as a referece. he, all but oe of the equatos of oto s wrtte wth that value o the rght-had sde: ( ) ω (7.0) 0 ( ω ) 0 ( ω) ( ) ω (7.) Exadg the frst ad secod equatos, drog the subscrts ad : ( ) ( ) ω 0 + ω 0 3 (7.a,b) Rewrtg wth the rato fro (7.a): ter o the rght-had sde ad solvg for the ( / ) ( ) ω (7.3)

7 Solvg for the ( 3 / ) rato fro (7.b): ω (7.4) ( ) ω (7.5) 3 ω ( ) 3 ω 3 ω ( ) (7.6) (7.7) 3 ( ω)( ω) (7.8) ω 3ω + (7.9) 4 3 We ow have the geeral equatos for the egevector values. If a value s chose for, say.0, the the two ratos above ca be solved for corresodg values of ad 3 for each of the three egevalues. Sce at each egevalue there are (+) uows ( ω, ) for a syste wth equatos of oto, the egevectors are oly ow as ratos of dslaceets, ot as absolute agtudes. For the frst ode of our tdof syste the uows are ω,, ad 3 ad we have oly three equatos of oto. Substtutg values for the three egevalues to the geeral egevector rato equatos above, assug, : For ode, ω 0 (7.30)

8 (7.3) 3 ( )( ) (7.3) (7.33) 3 Arbtrarly assgg : (7.34) Rgd-Body Mode, 0 rad/sec Fgure 7.: Mode shae lot for rgd body ode, where all asses ove together wth o stress the coectg srgs. For ode, ω 0 (7.35) 0 (7.36) 3 ( )( 0) (7.37)

9 (7.38) 3 0 (7.39) - Secod Mode, Mddle Mass Statoary, rad/sec Fgure 7.3: Mode shae lot for secod ode, ddle ass statoary ad the two ed asses ove out of hase wth each other wth equal altude. For ode 3, ω (7.40) (7.4) ( )( ) (7.4) (7.43) 3 3 (7.44)

10 - hrd Mode,.73 rad/sec Fgure 7.4: Mode shae lot for thrd ode, wth two ed asses ovg hase wth each other ad out of hase wth the ddle ass, whch s ovg wth twce the altude of the ed asses Iterretg Egevectors For the frst ode, f all the asses start wth ether ero or the sae tal velocty ad wth tal dslaceets of soe scalar ultle of [ ], where s the trasose, the syste wll ether rea at rest or wll cotue ovg at that velocty wth o relatve oto betwee the asses. For the secod ad thrd odes, f the syste s released wth ero tal veloctes but wth tal dslaceets of soe scalar ultle of that egevector, the the syste wll vbrate oly that ode wth all the asses reachg ther u ad axu ots at the sae ot te. Ay other cobato of tal dslaceets wll result a oto whch s a cobato of the three egevectors Modal Matrx Now that the three egevectors have bee defed, the odal atrx wll be troduced. he odal atrx s a (x) atrx wth colus corresodg to the syste egevectors, startg wth the frst ode the frst colu ad so o:

11 ode: 3 3 DOF 3 DOF (7.45) DOF 3 3 For our tdof roble: 0 (7.46) 7.3 Ucoulg the Equatos of Moto At ths ot the syste s well defed ters of atural frequeces ad odes of vbrato. If ay further forato such as traset or frequecy resose s desred, solvg for t would be laborous because the syste equatos are stll couled. For traset resose, the equatos would have to be solved sultaeously usg a uercal tegrato schee uless the roble were sle eough to allow a closed for soluto. o calculate the daed frequecy resose, a colex equato solvg route would have to be used to vert the colex coeffcet atrx at each frequecy. I order to facltate solvg for the traset or frequecy resoses, t s useful to trasfor the -couled secod order dfferetal equatos to - ucouled secod order dfferetal equatos by trasforg fro the hyscal coordate syste to a rcal coordate syste. I lear algebra ters, the trasforato fro hyscal to rcal coordates s ow as a chage of bass. here are ay otos for chage of bass, but we wll show that whe egevectors are used for the trasforato the rcal coordate syste has a hyscal eag; each of the ucouled sdof systes reresets the oto of a secfc ode of vbrato. he -ucouled equatos the rcal coordate syste ca the be solved for the resoses the rcal coordate syste usg well-ow solutos for sgle degree of freedo systes. he -resoses the rcal coordate syste ca the be trasfored bac to the hyscal coordate syste to rovde the actual

12 resose hyscal coordates. hs rocedure s show scheatcally Fgure 7.5. PHYSICAL COORDINAES Couled Equatos of Moto Ital Codtos Forcg Fuctos rasfor PRINCIPAL COORDINAES Ucouled Equatos of Moto Ital Codtos Forcg Fuctos Soluto Bac-rasfor PHYSICAL COORDINAES Soluto Fgure 7.5: Roada for Modal Soluto he rocedure above s aalogous to usg Lalace trasfors for solvg dfferetal equatos, where the dfferetal equato s trasfored to a algebrac equato, solved algebracally, ad bac trasfored to get the soluto of the orgal roble. We ow eed a eas of dagoalg the ass ad stffess atrces, whch wll yeld a set of ucouled equatos. he codto to guaratee dagoalato s the exstece of -learly deedet egevectors, whch s always the case f the ass ad stffess atrces are both syetrc or f there are -dfferet (oreeated) egevalues (Strag 998). Gog bac to the orgal hoogeeous equato of oto: && + 0 (7.47)

13 Havg oral odes eas that at frequecy : s ( ω t +φ ) (7.48) Dfferetatg twce to get accelerato: && ω s ( ω t +φ ) (7.49) Substtutg bac to the equato of oto: { } { } ω s( ω t + φ ) + s( ω t +φ ) 0 (7.50) Cacelg se ters: ω + (7.5) 0 Rearragg ad wrtg the above equato for both the th ad j th odes: ω (7.5) ω (7.53) j j j ad j are the th ad j th egevectors, the th ad j th colus of the odal atrx. Preultlyg (7.5) by the trasose of, : j j ω (7.54) j j ag the trasose of (7.53), where the trasose of a roduct s the roduct of the dvdual trasoses tae reverse order,.e., [ ] AB B A : j ωj j, (7.55) sce ad are syetrcal,, ad : ω (7.56) j j j

14 Postultlyg (7.56) by ω (7.57) j j j Now, subtractg (7.57) fro (7.54): ω j j j j j ( ω ) 0 ( ω ω ) j j (7.58) Whe / j, the ter ( ω ω j ) caot be equal to ero, eag that the ter j ust be equal to ero. j 0 (7.59) Loog at the ses of the atrces ultled: j x x x (7.60) Equato (7.59) ca be rewrtte: ( x) x( x) x( x) ( x) scalar (7.6) 0, (7.6) j j where j s a off-dagoal ter the ass atrx of the rcal coordate syste. he two egevectors j ad are sad to be orthogoal wth resect to, where orthogoalty s defed as the roerty that causes all the off-dagoal ters the rcal ass atrx to be ero. j, ω ω 0. hus the roduct ca be set equal to ay arbtrary costat,, a dagoal ter the rcal ass atrx. Returg to (7.6), for ( j )

15 (7.63) hs s where varous oralato techques for egevectors coe to lay, dscussed the ext secto. he stffess atrx,, s oraled the sae aer. I ractce, stead of dagoalg the ass ad stffess atrces ter by ter by re- ad ostultlyg by dvdual egevectors, the etre odal atrx s used to dagoale oe oerato usg two atrx ultlcatos: (7.64) (7.65) 7.4 Noralg Egevectors Because egevectors are oly ow as ratos of dslaceets, ot as absolute agtudes, we ca choose how to orale the. U to ow, whe calculatg egevectors we have arbtrarly set the altude of the frst dof to. We wll ow dscuss two of the ost cooly used egevector oralato techques. Dfferet oralg techques result dfferet fors of the resultg ucouled dfferetal equatos Noralg wth Resect to Uty Oe ethod s to orale wth resect to uty, ag the largest eleet each egevector equal to uty by dvdg each colu by ts largest value. We ow add the otato, where the refers to a oraled odal atrx (7.66) 0.5 Usg the uty oraled odal atrx to trasfor the ass atrx two atrx ultlcatos:

16 (7.67) (7.68) Slarly trasforg the stffess atrx: (7.69) (7.70) Note that the orgal flled stffess atrx s ow dagoal. Also ote that f the dagoal eleets of the stffess atrx (7.70) are dvded by the corresodg dagoal eleets of the ass atrx (7.69), the three ters are the squares of the resectve egevalues Noralg wth Resect to Mass Aother ethod s to orale wth resect to ass usg the equato:, (7.7).0 ag each dagoal ass ter equal.0. hs s the ethod used by default ANSYS. Oce aga, ote that odal atrx subscrt sgfes the oraled th egevector. Each oraled egevector s defed as follows: q (7.7)

17 Where q s defed as: q j j j (7.73) For a dagoal ass atrx, q ca be slfed sce all the ero: j ters are q (7.74) hus, by oeratg o by,the ass atrx should be trasfored to the detty atrx. Startg wth ad the q values fro above: 0 (7.75) () () () q () ( ) ( ) q () ( ) () q (7.76a,b,c) he odal atrx oraled wth resect to ass becoes: (7.77) Usg to trasfor the ass atrx:

18 (7.78) (7.79) (7.80) he orgal ass atrx has bee trasfored to the detty atrx. Slarly trasforg the stffess atrx:

19 (7.8) (7.8) (7.83)

20 (7.84) Note that the oraled stffess atrx s ow dagoal ad that the dagoal ters are the squares of the corresodg three egevalues. he oraled stffess atrx s also ow as the sectral atrx (Weaver 990). Because oralg wth resect to ass results a detty rcal ass atrx ad squares of the egevalues o the dagoal the rcal stffess atrx, we wll use oly ths oralato the future. Sce we ow the for of the rcal atrces whe oralg wth resect to ass, o ultlyg of odal atrces s actually requred: the hoogeeous rcal equatos of oto ca be wrtte by secto owg oly the egevalues. 7.5 Revewg Equatos of Moto Prcal Coordates Mass Noralato 7.5. Equatos of Moto Physcal Coordate Syste 0 0 && && && [ ] (7.85) Egevalues: ω 0 (7.86) ω ± ω ± 3 3 (7.87a,b) Egevectors, oraled wth resect to ass:

21 3 6 0 (7.88) Equatos of Moto Prcal Coordate Syste && 0 0 && && [ ] (7.89) Exadg Matrx Equatos of Moto Both Coordate Systes Physcal Coordates && + 0 && && Prcal Coordates && 0 && && hese equatos are couled ad have to be solved sultaeously. hese hoogeeous equatos are ucouled ad ca be solved deedetly. able 7.: Suary of equatos of oto hyscal ad rcal coordates.

22 7.6 rasforg Ital Codtos ad Forces Now that we ow how to costruct the hoogeeous ucouled equatos of oto for the syste, we eed to ow how to trasfor tal codtos ad forces to the rcal coordate syste. We ca the solve for traset ad forced resoses the rcal coordate syste usg the ucouled equatos. Startg wth the orgal o-hoogeeous equatos of oto hyscal coordates: Preultlyg by && + F (7.90), the trasose of the odal atrx: && + F (7.9) Isertg the detfy atrx, I : Rewrtg ad regroug ters: && { + { F (7.9) I I && { + {{ F {, (7.93) 3 && F where ad were show to dagoale the ass ad stffess atrces the revous secto. Defg ters: (x) dagoal rcal ass atrx (x) dagoal rcal stffess atrx && && accelerato vector rcal coordates dslaceet vector rcal coordates

23 F F force vector rcal coordates I the revous secto, the deftos for acceleratos ad dslaceets hyscal ad rcal coordates were show to be: && && (7.94) he sae relatoshs hold for tal codtos of dslaceet ad velocty: & o o & o o (7.95) I (7.95), o ad & o are vectors of tal dslaceets ad veloctes, resectvely, the rcal coordate syste, ad o ad & o are vectors of tal dslaceets ad veloctes, resectvely, the hyscal coordate syste. ag the verse of the odal atrx to covert tal codtos requres that the odal atrx be square, wth as ay egevectors as uber of degrees of freedo. We wll see future chaters that there are staces where ot all egevectors are avalable. I oe case, we ay choose to oly calculate egevalues ad egevectors u to a certa frequecy order to save calculato te or because the roble oly requres owledge of resose a certa frequecy rage. I aother case, we ay buld a reduced odel where oly the ost sgfcat odes are retaed. Fortuately, a large ajorty of real lfe robles volve ero tal codtos. 7.7 Suarg Equatos of Moto Both Coordate Systes he two sets of equatos, hyscal ad rcal coordates, are show able 7.:

24 Physcal Coordates Prcal Coordates && + F && + F 3 && + F IC's:,,,&,&,& 3 3 && F && + F 3 && F3 IC's:,,,&,&,& 3 3 able 7.: Suary of equatos of oto hyscal ad rcal coordates. he varables hyscal coordates are the ostos ad veloctes of the asses. he varables rcal coordates are the dslaceets ad veloctes of each ode of vbrato. he equatos rcal coordates ca be easly solved, sce the equatos are ucouled, yeldg the dslaceets. We ow eed to bac trasfor the results the rcal coordate syste to the hyscal coordate syste to get the fal aswer. 7.8 Bac-rasforg fro Prcal to Physcal Coordates We showed revously that the relatosh betwee hyscal ad rcal coordates s: (7.96) Preultlyg by : ( ) 3 I (7.97) (7.98) hus, the dslaceet vector hyscal coordates s obtaed by reultlyg the vector of dslaceets rcal coordates by the oraled odal atrx.

25 Slarly for velocty: & & (7.99) 7.9 Reducg the Model Se Whe Oly Selected Degrees of Freedo are Requred So far we have hted at the fact that oly ortos of the egevector atrx are eeded f selected dof s have forces aled ad other (or the sae) dof s are eeded for outut. hs secto wll show how the reducto dof s occurs. hs reducto s oe of the ey stes to be used later the boo whe we cover how to reduce the se of odels derved fro large fte eleet sulatos. Revewg the stes the odal soluto, startg wth the equatos of oto ad tal codtos hyscal coordates: && + F && + F 3 && + F (7.00) Ital Codtos:,,,&,&,& Solve for egevalues: ω, ω, ω 3 Solve for egevectors, orale wth resect to ass ad for the odal atrx fro colus of egevectors: 3 3 (7.0) rasfor forces fro hyscal to rcal coordates: F F (7.0) Wrte the equatos of oto rcal coordates:

26 && F && +ω F && +ω F IC's:,,,&,&,& (7.03a,b,c,d) Solve the equatos rcal coordates ether te or frequecy doa ad the bac trasfor to hyscal coordates: & & (7.04) Note that the two crtcal trasforatos (assug ero tal codtos) volve reultlyg by the trasose of the odal atrx ( F F ) (7.0) or the odal atrx ( ) (7.04). Let us frst exae the force trasforato by exadg the equatos: F F (7.05) 3 F 3 F 3 3 F F F F F F + F + 3F3 F F 3F F + 3F + 33F3 (7.06) Note that the ultlers of F the frst colu are the eleets of the frst row of the odal atrx, the ultlers of F the secod colu are the eleets of the secod row of the odal atrx ad the ultlers of F 3 the thrd colu are the eleets of the thrd row of the odal atrx. Suose that force s to be aled at oly ass, F, the oly the frst row of the odal atrx s requred to trasfor the force hyscal coordates to the force rcal coordates.

27 Now let us exae the dslaceet trasforato by exadg the equatos: (7.07) (7.08) Note that the coeffcets of the rcal dslaceets the frst row above are the eleets of the frst row of the odal atrx. Slarly, coeffcets of the secod ad thrd rows are the eleets of the secod ad thrd rows of the odal atrx. Suose that the oly hyscal dslaceet we are terested s that of ass,, the oly the secod row of the odal atrx s requred to trasfor the three dslaceets,, 3 rcal coordates to. hs leads to the followg cocluso about reducg the se of the odel: Oly the rows of the odal atrx that corresod to degrees of freedo to whch forces are aled ad/or for whch dslaceets are desred are requred to colete the odel. For ths tdof odel, reducg the se of the roble s ot requred; however, we wll see later that a realstc fte eleet odel, wth hudreds of thousads of degrees of freedo, resets a etrely dfferet roble. Havg the ablty to reduce the roble se s crtcal order to use the detaled results of a colcated fte eleet odel to rovde accurate results a lower order MALAB odel. 7.0 Dag Systes wth Prcal Modes 7.0. Overvew Dag colex bult-u echacal systes s ossble to redct wth the reset state of the art. We wll dscuss ths secto the codtos whch detere f a dag atrx ca be dagoaled, ad the crtero to eable the daed equatos to be dagoaled. I geeral, a arbtrary dag atrx caot be dagoaled by the udaed egevectors, as the

28 ass ad stffess atrces ca. hs leads to usg what s called roortoal dag ost fte eleet sulatos. If a echacal syste s desged wth a secfc vscous dag eleet, for exale a dashot, that doates the sall aout of heret structural dag reset, the that eleet ca be added to the syste as a vscous daer. he resultg syste s lear, but robably does ot exhbt oral odes as dscussed Secto 7... I geeral ths leads to the ablty to dagoale ad ucoule the equatos of oto, requrg a state sace soluto of the orgal, couled equatos of oto. Vscoelastc dag treatets (dag elastoers) have bee used for years ds drves, ost tycally as costraed layer daers o the th sheet etal susesos whch suort the read/wrte head. he effect of ths vscoelastc dag ca be aroxated at a secfc teerature ad frequecy as roortoal dag by usg the odal stra eergy techque assocato wth a fte eleet structural odel (Johso 98). Igorg secfc vscous, coulob, ad vscoelastc dag eleets, dag tycal structures arses fro hysteress losses the aterals as they are straed, soe cases fro vscous losses due to structure/flud teracto but ore ortatly fro relatve oto at the terfaces ad boudares where dfferet arts are attached or grouded. Uless a secfc dag eleet s used a structural desg, ost structures have dag whch vares fro ode to ode ad wll be the rage of 0.05% to % of crtcal dag. he odes ths chater are all real or oral odes as defed earler. Oce aga, havg oral odes eas that at certa frequeces all ots the syste wll vbrate at the sae frequecy ad hase,.e., all ots the syste wll reach ther u ad axu dslaceets at the sae ot te. Chater 5 dscussed colex odes, odes whch all ots the syste do ot reach ther u ad axu dslaceets at the sae ot te Codtos Necessary for Exstece of Prcal Modes Daed Syste Wth a coservatve (o dag) syste, oral odes of vbrato wll exst. I order to have oral odes a daed syste, the ode shaes ust be the sae as for the udaed case, ad the varous arts of the syste ust ass through ther u ad axu ostos at the sae stat te, exressed as:

29 cos ( t ) th ω +φ for the ode (7.09) A suffcet codto for the exstece of daed oral odes s that the dag atrx be a lear cobato of the ass ad stffess atrces. We ow that ad are dagoaled by oeratg o the wth the odal atrx. Whe c s a lear cobato of ad, the the dag atrx c s also ucouled (dagoaled) by the sae re- ad ostultlcato oeratos by the odal atrx as wth the ad atrces (Weaver 990, Crag 98). he daed equatos of oto the becoe: where the dag atrx s a lear cobato of && + c& + F, (7.0) ad : c a+ b (7.) c c, (7.) ad where s the oraled (wth resect to ass) odal atrx. Wrtg out the colete equato: && + c& + F (7.3) { + + {{ {{ { && c & F (7.4) 3 I && c & F Loog at the c to c coverso where c a+ b : c a+ b (7.5) c a + b ai+ b, (7.6) where s a dagoal atrx whose eleets are the squares of the egevalues.

30 he equato for the th ode s: Rewrtg, defg ( ) && + a+ bω & +ω F (7.7) c, the ( a b ) + ω ter, usg otato: c a+ bω ζ ω (7.8) Where ζ s the ercetage of crtcal dag for the th ode, defed as: c c c ζ c ω cr (7.9) he: a+ bω ζ ω (7.0) Rewrtg the equato rcal coordates: && + ζ ω & +ω F (7.) hs tye of dag s ow as roortoal dag, where the dag for each ode (they ca all be dfferet) s roortoal to the crtcal dag for that ode. Sce the dag s also roortoal to velocty, t s of a vscous ature. If the sae dag value s used for all odes, t wll be referred to as ufor dag. Dag whch the dag value for each ode ca be set dvdually wll be referred to as o-ufor dag Dfferet yes of Dag Sle Proortoal Dag Vscous dag each ode s tae to be a arbtrary ercetage, ζ, of crtcal dag:

31 && + ζω & +ω F + ζ & + && F (7.) hs s aalogous to the falar otato used for a sgle degree of freedo syste: && + c& + F c F && + & + (7.3) Defe crtcal dag c cr ad defe the ter ultlyg velocty to be: c ζω c c cr c c (7.4) Rewrtg: F && + ζω & +ω (7.5) Proortoal to Stffess Matrx Relatve Dag Recogg that the hgher odes of vbrato da out qucly, relatve dag yelds dag roorto to frequeces oral odes, bascally lettg the a ter for ζ go to ero: a+ bω bω ζ ω a 0 7.6) If a value of ζ, for the frst ode, s assued, a value ca be defed for b :

32 b ζ, (7.7) ω ad the value for ay other ode s: ω ζ ζ (7.8) ω Proortoal to Mass Matrx Absolute Dag Absolute dag s based o ag b equal to ero, whch case the ercetage of crtcal dag s versely roortoal to the atural frequecy of each ode. hs wll gve decreasg dag for odes as ther frequeces crease. a+ bω a ζ ω ω b 0 (7.9) If a value of ζ, for the frst ode, s assued, a value ca be defed for a : a ζ ω, (7.30) ad the value for ay other ode s: ω ζ ζ (7.3) ω

33 7.0.4 Defg Dag Matrx Whe Proortoal Dag s Assued Fgure 7.6: wo degree of freedo for dag exale. F F c c c 3 A terestg questo to as s what the eleets of the dag atrx should be the two degree of freedo (dof) roble show Fgure 7.6 order to be able to dagoale the equatos of oto. We wll use the egevectors fro the udaed case to orale the dag atrx. he we wll solve for the secfc values of the dvdual daers whch wll allow the dagoalato. We wll see how o-tutve the values of c,c adc 3 are order to be able to dagoale. (See Crag [98] for a geeral exresso to calculate the hyscal dag atrx whe gve roortoal dag values, the orgal ass atrx, the dagoaled ass atrx ad the egevalues ad egevectors.) Solvg for Dag Values Startg wth the udaed egevalues ad egevectors: 0 c + c c 0 c c c + c (7.3) Solve for the dagoaled dag atrx, assug roortoal dag, ad owg that the dagoaled stffess atrx eleets are squares of the egevalues:

34 ω 0 c c ζ ζ 0 ω (7.33) Preultlyg by ( ) ad ostultlyg by ( ) : ( ) c ( ) ζ( ) ( ) 443 I 443 I (7.34) ( ) ( ) c ζ (7.35) Solvg for the verses above, otg that for ths dof syste, the erforg the oeratos o :, ad he verse of a x atrx ca be foud by:. Iterchagg the two dagoal eleets.. Chagg the sgs of the two off-dagoal eleets. 3. Dvdg by the deterat of the orgal atrx. d b a b c a c d a b c d able 7.: Iverse of x atrx. ( ) ( ) (7.36)

35 (7.37) (7.38) c ζ ζ (7.39) Now we ca solve for the secfc values for the three daers: ( ) c ζ 3 ( ) c ζ 3 (7.40) ζ ( ).73

36 3 ( ) c + c c + c ζ + 3 ( ) c c ζ + 3 c 3 ( ) ( ) ζ (7.4) ζ ( ) ζ Suarg: c c3 ζ (7.4) ( ) c ζ.73 (7.43) Note that the values for the three daers are ot at all tutve ad would have bee very dffcult f ossble to guess to be able to costruct a dagoalable dag atrx. If defg the dagoalable dag atrx for ths x roble s dffcult, age tryg to defe t for a real lfe fte eleet roble wth thousads of degrees of freedo. Also, t s hghly robable that the bac-calculated dag values hyscal coordates would atch the actual dag the structure Checg Raylegh For of Dag Matrx We have ow defed the values of the c,c adc 3, daers whch allow dagoalg the equatos of oto. Aother terestg questo s whether the Raylegh for has bee satsfed: Is c a lear cobato of ad? + 3 3? 0 c ζ a + b (7.44) + We have two uows, a ad b, ad essetally two equatos, sce the two dagoal eleets are the sae ad the two off dagoal eleets are the sae. Frst, let us loo at the two off dagoal ters, equatg ters o the two sdes above:

37 ( ) ( ) ( ) ζ 3 a 0 + b ( ) 3 b ζ ζ 3 ( ) (7.45) Now, equatg the dagoal ters: ( ) ζ + 3 a + b a + ζ ( 3 ) (7.46) ( ) a + ζ 3 ( ) ( ) a ζ + 3 ζ 3 ζ (7.47) ξ 3 3 a ζ 3 3 (7.48) Checg the two values for a ad b by substtutg bac to (7.46).

38 ζ ζ ( ) ( ) +ζ 3 ( )( ) ζ + 0 ( 3 3 ) ζ (7.49) So c s a lear cobato of ad ad the Raylegh crtero holds. Probles Note: All the robles refer to the two dof syste show Fgure P.. P7. Set, ad solve for the egevalues ad egevectors of the udaed syste. Norale the egevectors to uty, wrte out the odal atrx ad had lot the ode shaes P7. Norale the egevectors P7. wth resect to ass ad dagoale the ass ad stffess atrces. Idetfy the ters the oraled ass ad stffess atrces. Wrte the hoogeeous equatos of oto hyscal ad rcal coordates. P7.3 Covert the followg ste forcg fucto ad tal codtos hyscal coordates to rcal coordates: a) F, F 3 0,&,,& b) P7.4 Usg the results of P7. ad P7.3, wrte the colete equatos of oto hyscal ad rcal coordates assug roortoal dag.

PRACTICAL CONSIDERATIONS IN HUMAN-INDUCED VIBRATION

PRACTICAL CONSIDERATIONS IN HUMAN-INDUCED VIBRATION PRACTICAL CONSIDERATIONS IN HUMAN-INDUCED VIBRATION Bars Erkus, 4 March 007 Itroducto Ths docuet provdes a revew of fudaetal cocepts structural dyacs ad soe applcatos hua-duced vbrato aalyss ad tgato of

More information

Ahmed Elgamal. MDOF Systems & Modal Analysis

Ahmed Elgamal. MDOF Systems & Modal Analysis DOF Systems & odal Aalyss odal Aalyss (hese otes cover sectos from Ch. 0, Dyamcs of Structures, Al Chopra, Pretce Hall, 995). Refereces Dyamcs of Structures, Al K. Chopra, Pretce Hall, New Jersey, ISBN

More information

Chapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II

Chapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II CEE49b Chapter - Free Vbrato of Mult-Degree-of-Freedom Systems - II We ca obta a approxmate soluto to the fudametal atural frequecy through a approxmate formula developed usg eergy prcples by Lord Raylegh

More information

7.0 Equality Contraints: Lagrange Multipliers

7.0 Equality Contraints: Lagrange Multipliers Systes Optzato 7.0 Equalty Cotrats: Lagrage Multplers Cosder the zato of a o-lear fucto subject to equalty costrats: g f() R ( ) 0 ( ) (7.) where the g ( ) are possbly also olear fuctos, ad < otherwse

More information

for each of its columns. A quick calculation will verify that: thus m < dim(v). Then a basis of V with respect to which T has the form: A

for each of its columns. A quick calculation will verify that: thus m < dim(v). Then a basis of V with respect to which T has the form: A Desty of dagoalzable square atrces Studet: Dael Cervoe; Metor: Saravaa Thyagaraa Uversty of Chcago VIGRE REU, Suer 7. For ths etre aer, we wll refer to V as a vector sace over ad L(V) as the set of lear

More information

2/20/2013. Topics. Power Flow Part 1 Text: Power Transmission. Power Transmission. Power Transmission. Power Transmission

2/20/2013. Topics. Power Flow Part 1 Text: Power Transmission. Power Transmission. Power Transmission. Power Transmission /0/0 Topcs Power Flow Part Text: 0-0. Power Trassso Revsted Power Flow Equatos Power Flow Proble Stateet ECEGR 45 Power Systes Power Trassso Power Trassso Recall that for a short trassso le, the power

More information

Non-degenerate Perturbation Theory

Non-degenerate Perturbation Theory No-degeerate Perturbato Theory Proble : H E ca't solve exactly. But wth H H H' H" L H E Uperturbed egevalue proble. Ca solve exactly. E Therefore, kow ad. H ' H" called perturbatos Copyrght Mchael D. Fayer,

More information

Some Different Perspectives on Linear Least Squares

Some Different Perspectives on Linear Least Squares Soe Dfferet Perspectves o Lear Least Squares A stadard proble statstcs s to easure a respose or depedet varable, y, at fed values of oe or ore depedet varables. Soetes there ests a deterstc odel y f (,,

More information

KURODA S METHOD FOR CONSTRUCTING CONSISTENT INPUT-OUTPUT DATA SETS. Peter J. Wilcoxen. Impact Research Centre, University of Melbourne.

KURODA S METHOD FOR CONSTRUCTING CONSISTENT INPUT-OUTPUT DATA SETS. Peter J. Wilcoxen. Impact Research Centre, University of Melbourne. KURODA S METHOD FOR CONSTRUCTING CONSISTENT INPUT-OUTPUT DATA SETS by Peter J. Wlcoxe Ipact Research Cetre, Uversty of Melboure Aprl 1989 Ths paper descrbes a ethod that ca be used to resolve cossteces

More information

( q Modal Analysis. Eigenvectors = Mode Shapes? Eigenproblem (cont) = x x 2 u 2. u 1. x 1 (4.55) vector and M and K are matrices.

( q Modal Analysis. Eigenvectors = Mode Shapes? Eigenproblem (cont) = x x 2 u 2. u 1. x 1 (4.55) vector and M and K are matrices. 4.3 - Modal Aalyss Physcal coordates are ot always the easest to work Egevectors provde a coveet trasformato to modal coordates Modal coordates are lear combato of physcal coordates Say we have physcal

More information

MOLECULAR VIBRATIONS

MOLECULAR VIBRATIONS MOLECULAR VIBRATIONS Here we wsh to vestgate molecular vbratos ad draw a smlarty betwee the theory of molecular vbratos ad Hückel theory. 1. Smple Harmoc Oscllator Recall that the eergy of a oe-dmesoal

More information

2006 Jamie Trahan, Autar Kaw, Kevin Martin University of South Florida United States of America

2006 Jamie Trahan, Autar Kaw, Kevin Martin University of South Florida United States of America SOLUTION OF SYSTEMS OF SIMULTANEOUS LINEAR EQUATIONS Gauss-Sedel Method 006 Jame Traha, Autar Kaw, Kev Mart Uversty of South Florda Uted States of Amerca kaw@eg.usf.edu Itroducto Ths worksheet demostrates

More information

Department of Mathematics UNIVERSITY OF OSLO. FORMULAS FOR STK4040 (version 1, September 12th, 2011) A - Vectors and matrices

Department of Mathematics UNIVERSITY OF OSLO. FORMULAS FOR STK4040 (version 1, September 12th, 2011) A - Vectors and matrices Deartet of Matheatcs UNIVERSITY OF OSLO FORMULAS FOR STK4040 (verso Seteber th 0) A - Vectors ad atrces A) For a x atrx A ad a x atrx B we have ( AB) BA A) For osgular square atrces A ad B we have ( )

More information

Chapter 9 Jordan Block Matrices

Chapter 9 Jordan Block Matrices Chapter 9 Jorda Block atrces I ths chapter we wll solve the followg problem. Gve a lear operator T fd a bass R of F such that the matrx R (T) s as smple as possble. f course smple s a matter of taste.

More information

A New Method for Solving Fuzzy Linear. Programming by Solving Linear Programming

A New Method for Solving Fuzzy Linear. Programming by Solving Linear Programming ppled Matheatcal Sceces Vol 008 o 50 7-80 New Method for Solvg Fuzzy Lear Prograg by Solvg Lear Prograg S H Nasser a Departet of Matheatcs Faculty of Basc Sceces Mazadara Uversty Babolsar Ira b The Research

More information

2SLS Estimates ECON In this case, begin with the assumption that E[ i

2SLS Estimates ECON In this case, begin with the assumption that E[ i SLS Estmates ECON 3033 Bll Evas Fall 05 Two-Stage Least Squares (SLS Cosder a stadard lear bvarate regresso model y 0 x. I ths case, beg wth the assumto that E[ x] 0 whch meas that OLS estmates of wll

More information

Lecture 07: Poles and Zeros

Lecture 07: Poles and Zeros Lecture 07: Poles ad Zeros Defto of poles ad zeros The trasfer fucto provdes a bass for determg mportat system respose characterstcs wthout solvg the complete dfferetal equato. As defed, the trasfer fucto

More information

On the characteristics of partial differential equations

On the characteristics of partial differential equations Sur les caractérstques des équatos au dérvées artelles Bull Soc Math Frace 5 (897) 8- O the characterstcs of artal dfferetal equatos By JULES BEUDON Traslated by D H Delhech I a ote that was reseted to

More information

2. Independence and Bernoulli Trials

2. Independence and Bernoulli Trials . Ideedece ad Beroull Trals Ideedece: Evets ad B are deedet f B B. - It s easy to show that, B deedet mles, B;, B are all deedet ars. For examle, ad so that B or B B B B B φ,.e., ad B are deedet evets.,

More information

Some results and conjectures about recurrence relations for certain sequences of binomial sums.

Some results and conjectures about recurrence relations for certain sequences of binomial sums. Soe results ad coectures about recurrece relatos for certa sequeces of boal sus Joha Cgler Faultät für Matheat Uverstät We A-9 We Nordbergstraße 5 Joha Cgler@uveacat Abstract I a prevous paper [] I have

More information

Random Variables. ECE 313 Probability with Engineering Applications Lecture 8 Professor Ravi K. Iyer University of Illinois

Random Variables. ECE 313 Probability with Engineering Applications Lecture 8 Professor Ravi K. Iyer University of Illinois Radom Varables ECE 313 Probablty wth Egeerg Alcatos Lecture 8 Professor Rav K. Iyer Uversty of Illos Iyer - Lecture 8 ECE 313 Fall 013 Today s Tocs Revew o Radom Varables Cumulatve Dstrbuto Fucto (CDF

More information

1 Onto functions and bijections Applications to Counting

1 Onto functions and bijections Applications to Counting 1 Oto fuctos ad bectos Applcatos to Coutg Now we move o to a ew topc. Defto 1.1 (Surecto. A fucto f : A B s sad to be surectve or oto f for each b B there s some a A so that f(a B. What are examples of

More information

A Conventional Approach for the Solution of the Fifth Order Boundary Value Problems Using Sixth Degree Spline Functions

A Conventional Approach for the Solution of the Fifth Order Boundary Value Problems Using Sixth Degree Spline Functions Appled Matheatcs, 1, 4, 8-88 http://d.do.org/1.4/a.1.448 Publshed Ole Aprl 1 (http://www.scrp.org/joural/a) A Covetoal Approach for the Soluto of the Ffth Order Boudary Value Probles Usg Sth Degree Sple

More information

MATH 247/Winter Notes on the adjoint and on normal operators.

MATH 247/Winter Notes on the adjoint and on normal operators. MATH 47/Wter 00 Notes o the adjot ad o ormal operators I these otes, V s a fte dmesoal er product space over, wth gve er * product uv, T, S, T, are lear operators o V U, W are subspaces of V Whe we say

More information

Long blade vibration model for turbine-generator shafts torsional vibration analysis

Long blade vibration model for turbine-generator shafts torsional vibration analysis Avalable ole www.ocpr.co Joural of Checal ad Pharaceutcal Research, 05, 7(3):39-333 Research Artcle ISSN : 0975-7384 CODEN(USA) : JCPRC5 Log blade vbrato odel for turbe-geerator shafts torsoal vbrato aalyss

More information

d dt d d dt dt Also recall that by Taylor series, / 2 (enables use of sin instead of cos-see p.27 of A&F) dsin

d dt d d dt dt Also recall that by Taylor series, / 2 (enables use of sin instead of cos-see p.27 of A&F) dsin Learzato of the Swg Equato We wll cover sectos.5.-.6 ad begg of Secto 3.3 these otes. 1. Sgle mache-fte bus case Cosder a sgle mache coected to a fte bus, as show Fg. 1 below. E y1 V=1./_ Fg. 1 The admttace

More information

Algorithms behind the Correlation Setting Window

Algorithms behind the Correlation Setting Window Algorths behd the Correlato Settg Wdow Itroducto I ths report detaled forato about the correlato settg pop up wdow s gve. See Fgure. Ths wdow s obtaed b clckg o the rado butto labelled Kow dep the a scree

More information

Channel Models with Memory. Channel Models with Memory. Channel Models with Memory. Channel Models with Memory

Channel Models with Memory. Channel Models with Memory. Channel Models with Memory. Channel Models with Memory Chael Models wth Memory Chael Models wth Memory Hayder radha Electrcal ad Comuter Egeerg Mchga State Uversty I may ractcal etworkg scearos (cludg the Iteret ad wreless etworks), the uderlyg chaels are

More information

Stationary states of atoms and molecules

Stationary states of atoms and molecules Statoary states of atos ad olecules I followg wees the geeral aspects of the eergy level structure of atos ad olecules that are essetal for the terpretato ad the aalyss of spectral postos the rotatoal

More information

The Number of the Two Dimensional Run Length Constrained Arrays

The Number of the Two Dimensional Run Length Constrained Arrays 2009 Iteratoal Coferece o Mache Learg ad Coutg IPCSIT vol.3 (20) (20) IACSIT Press Sgaore The Nuber of the Two Desoal Ru Legth Costraed Arrays Tal Ataa Naohsa Otsua 2 Xuerog Yog 3 School of Scece ad Egeerg

More information

Mu Sequences/Series Solutions National Convention 2014

Mu Sequences/Series Solutions National Convention 2014 Mu Sequeces/Seres Solutos Natoal Coveto 04 C 6 E A 6C A 6 B B 7 A D 7 D C 7 A B 8 A B 8 A C 8 E 4 B 9 B 4 E 9 B 4 C 9 E C 0 A A 0 D B 0 C C Usg basc propertes of arthmetc sequeces, we fd a ad bm m We eed

More information

ELEMENTS OF NUMBER THEORY. In the following we will use mainly integers and positive integers. - the set of integers - the set of positive integers

ELEMENTS OF NUMBER THEORY. In the following we will use mainly integers and positive integers. - the set of integers - the set of positive integers ELEMENTS OF NUMBER THEORY I the followg we wll use aly tegers a ostve tegers Ζ = { ± ± ± K} - the set of tegers Ν = { K} - the set of ostve tegers Oeratos o tegers: Ato Each two tegers (ostve tegers) ay

More information

3D Reconstruction from Image Pairs. Reconstruction from Multiple Views. Computing Scene Point from Two Matching Image Points

3D Reconstruction from Image Pairs. Reconstruction from Multiple Views. Computing Scene Point from Two Matching Image Points D Recostructo fro Iage ars Recostructo fro ultple Ves Dael Deetho Fd terest pots atch terest pots Copute fudaetal atr F Copute caera atrces ad fro F For each atchg age pots ad copute pot scee Coputg Scee

More information

Derivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations

Derivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations Dervato of -Pot Block Method Formula for Solvg Frst Order Stff Ordary Dfferetal Equatos Kharul Hamd Kharul Auar, Kharl Iskadar Othma, Zara Bb Ibrahm Abstract Dervato of pot block method formula wth costat

More information

SEMI-TIED FULL-COVARIANCE MATRICES FOR HMMS

SEMI-TIED FULL-COVARIANCE MATRICES FOR HMMS SEMI-TIED FULL-COVARIANCE MATRICES FOR HMMS M.J.F. Gales fg@eg.ca.ac.uk Deceber 9, 997 Cotets Bass. Block Dagoal Matrces : : : : : : : : : : : : : : : : : : : : : : : : : : : :. Cooly used atrx dervatve

More information

Standard Deviation for PDG Mass Data

Standard Deviation for PDG Mass Data 4 Dec 06 Stadard Devato for PDG Mass Data M. J. Gerusa Retred, 47 Clfde Road, Worghall, HP8 9JR, UK. gerusa@aol.co, phoe: +(44) 844 339754 Abstract Ths paper aalyses the data for the asses of eleetary

More information

Debabrata Dey and Atanu Lahiri

Debabrata Dey and Atanu Lahiri RESEARCH ARTICLE QUALITY COMPETITION AND MARKET SEGMENTATION IN THE SECURITY SOFTWARE MARKET Debabrata Dey ad Atau Lahr Mchael G. Foster School of Busess, Uersty of Washgto, Seattle, Seattle, WA 9895 U.S.A.

More information

C-1: Aerodynamics of Airfoils 1 C-2: Aerodynamics of Airfoils 2 C-3: Panel Methods C-4: Thin Airfoil Theory

C-1: Aerodynamics of Airfoils 1 C-2: Aerodynamics of Airfoils 2 C-3: Panel Methods C-4: Thin Airfoil Theory ROAD MAP... AE301 Aerodyamcs I UNIT C: 2-D Arfols C-1: Aerodyamcs of Arfols 1 C-2: Aerodyamcs of Arfols 2 C-3: Pael Methods C-4: Th Arfol Theory AE301 Aerodyamcs I Ut C-3: Lst of Subects Problem Solutos?

More information

A Family of Non-Self Maps Satisfying i -Contractive Condition and Having Unique Common Fixed Point in Metrically Convex Spaces *

A Family of Non-Self Maps Satisfying i -Contractive Condition and Having Unique Common Fixed Point in Metrically Convex Spaces * Advaces Pure Matheatcs 0 80-84 htt://dxdoorg/0436/a04036 Publshed Ole July 0 (htt://wwwscrporg/oural/a) A Faly of No-Self Mas Satsfyg -Cotractve Codto ad Havg Uque Coo Fxed Pot Metrcally Covex Saces *

More information

are positive, and the pair A, B is controllable. The uncertainty in is introduced to model control failures.

are positive, and the pair A, B is controllable. The uncertainty in is introduced to model control failures. Lectue 4 8. MRAC Desg fo Affe--Cotol MIMO Systes I ths secto, we cosde MRAC desg fo a class of ult-ut-ult-outut (MIMO) olea systes, whose lat dyacs ae lealy aaetezed, the ucetates satsfy the so-called

More information

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem Joural of Amerca Scece ;6( Cubc Nopolyomal Sple Approach to the Soluto of a Secod Order Two-Pot Boudary Value Problem W.K. Zahra, F.A. Abd El-Salam, A.A. El-Sabbagh ad Z.A. ZAk * Departmet of Egeerg athematcs

More information

Basic Concepts in Numerical Analysis November 6, 2017

Basic Concepts in Numerical Analysis November 6, 2017 Basc Cocepts Nuercal Aalyss Noveber 6, 7 Basc Cocepts Nuercal Aalyss Larry Caretto Mecacal Egeerg 5AB Sear Egeerg Aalyss Noveber 6, 7 Outle Revew last class Mdter Exa Noveber 5 covers ateral o deretal

More information

Unit 9. The Tangent Bundle

Unit 9. The Tangent Bundle Ut 9. The Taget Budle ========================================================================================== ---------- The taget sace of a submafold of R, detfcato of taget vectors wth dervatos at

More information

MATH 371 Homework assignment 1 August 29, 2013

MATH 371 Homework assignment 1 August 29, 2013 MATH 371 Homework assgmet 1 August 29, 2013 1. Prove that f a subset S Z has a smallest elemet the t s uque ( other words, f x s a smallest elemet of S ad y s also a smallest elemet of S the x y). We kow

More information

ESS Line Fitting

ESS Line Fitting ESS 5 014 17. Le Fttg A very commo problem data aalyss s lookg for relatoshpetwee dfferet parameters ad fttg les or surfaces to data. The smplest example s fttg a straght le ad we wll dscuss that here

More information

Solutions to problem set ); (, ) (

Solutions to problem set ); (, ) ( Solutos to proble set.. L = ( yp p ); L = ( p p ); y y L, L = yp p, p p = yp p, + p [, p ] y y y = yp + p = L y Here we use for eaple that yp, p = yp p p yp = yp, p = yp : factors that coute ca be treated

More information

1 Lyapunov Stability Theory

1 Lyapunov Stability Theory Lyapuov Stablty heory I ths secto we cosder proofs of stablty of equlbra of autoomous systems. hs s stadard theory for olear systems, ad oe of the most mportat tools the aalyss of olear systems. It may

More information

Sebastián Martín Ruiz. Applications of Smarandache Function, and Prime and Coprime Functions

Sebastián Martín Ruiz. Applications of Smarandache Function, and Prime and Coprime Functions Sebastá Martí Ruz Alcatos of Saradache Fucto ad Pre ad Core Fuctos 0 C L f L otherwse are core ubers Aerca Research Press Rehoboth 00 Sebastá Martí Ruz Avda. De Regla 43 Choa 550 Cadz Sa Sarada@telele.es

More information

Engineering Vibration 1. Introduction

Engineering Vibration 1. Introduction Egeerg Vbrato. Itroducto he study of the moto of physcal systems resultg from the appled forces s referred to as dyamcs. Oe type of dyamcs of physcal systems s vbrato, whch the system oscllates about certa

More information

CHAPTER 4 RADICAL EXPRESSIONS

CHAPTER 4 RADICAL EXPRESSIONS 6 CHAPTER RADICAL EXPRESSIONS. The th Root of a Real Number A real umber a s called the th root of a real umber b f Thus, for example: s a square root of sce. s also a square root of sce ( ). s a cube

More information

Training Sample Model: Given n observations, [[( Yi, x i the sample model can be expressed as (1) where, zero and variance σ

Training Sample Model: Given n observations, [[( Yi, x i the sample model can be expressed as (1) where, zero and variance σ Stat 74 Estmato for Geeral Lear Model Prof. Goel Broad Outle Geeral Lear Model (GLM): Trag Samle Model: Gve observatos, [[( Y, x ), x = ( x,, xr )], =,,, the samle model ca be exressed as Y = µ ( x, x,,

More information

We have already referred to a certain reaction, which takes place at high temperature after rich combustion.

We have already referred to a certain reaction, which takes place at high temperature after rich combustion. ME 41 Day 13 Topcs Chemcal Equlbrum - Theory Chemcal Equlbrum Example #1 Equlbrum Costats Chemcal Equlbrum Example #2 Chemcal Equlbrum of Hot Bured Gas 1. Chemcal Equlbrum We have already referred to a

More information

TESTS BASED ON MAXIMUM LIKELIHOOD

TESTS BASED ON MAXIMUM LIKELIHOOD ESE 5 Toy E. Smth. The Basc Example. TESTS BASED ON MAXIMUM LIKELIHOOD To llustrate the propertes of maxmum lkelhood estmates ad tests, we cosder the smplest possble case of estmatg the mea of the ormal

More information

L5 Polynomial / Spline Curves

L5 Polynomial / Spline Curves L5 Polyomal / Sple Curves Cotets Coc sectos Polyomal Curves Hermte Curves Bezer Curves B-Sples No-Uform Ratoal B-Sples (NURBS) Mapulato ad Represetato of Curves Types of Curve Equatos Implct: Descrbe a

More information

The theoretical background of

The theoretical background of he theoretcal backgroud of -echologes he theoretcal backgroud of FactSage he followg sldes gve a abrdged overvew of the ajor uderlyg prcples of the calculatoal odules of FactSage. -echologes he bbs Eergy

More information

3.1 Introduction to Multinomial Logit and Probit

3.1 Introduction to Multinomial Logit and Probit ES3008 Ecooetrcs Lecture 3 robt ad Logt - Multoal 3. Itroducto to Multoal Logt ad robt 3. Estato of β 3. Itroducto to Multoal Logt ad robt The ultoal Logt odel s used whe there are several optos (ad therefore

More information

4 Inner Product Spaces

4 Inner Product Spaces 11.MH1 LINEAR ALGEBRA Summary Notes 4 Ier Product Spaces Ier product s the abstracto to geeral vector spaces of the famlar dea of the scalar product of two vectors or 3. I what follows, keep these key

More information

Block-Based Compact Thermal Modeling of Semiconductor Integrated Circuits

Block-Based Compact Thermal Modeling of Semiconductor Integrated Circuits Block-Based Compact hermal Modelg of Semcoductor Itegrated Crcuts Master s hess Defese Caddate: Jg Ba Commttee Members: Dr. Mg-Cheg Cheg Dr. Daqg Hou Dr. Robert Schllg July 27, 2009 Outle Itroducto Backgroud

More information

Assignment 5/MATH 247/Winter Due: Friday, February 19 in class (!) (answers will be posted right after class)

Assignment 5/MATH 247/Winter Due: Friday, February 19 in class (!) (answers will be posted right after class) Assgmet 5/MATH 7/Wter 00 Due: Frday, February 9 class (!) (aswers wll be posted rght after class) As usual, there are peces of text, before the questos [], [], themselves. Recall: For the quadratc form

More information

Functions of Random Variables

Functions of Random Variables Fuctos of Radom Varables Chapter Fve Fuctos of Radom Varables 5. Itroducto A geeral egeerg aalyss model s show Fg. 5.. The model output (respose) cotas the performaces of a system or product, such as weght,

More information

Unique Common Fixed Point of Sequences of Mappings in G-Metric Space M. Akram *, Nosheen

Unique Common Fixed Point of Sequences of Mappings in G-Metric Space M. Akram *, Nosheen Vol No : Joural of Facult of Egeerg & echolog JFE Pages 9- Uque Coo Fed Pot of Sequeces of Mags -Metrc Sace M. Ara * Noshee * Deartet of Matheatcs C Uverst Lahore Pasta. Eal: ara7@ahoo.co Deartet of Matheatcs

More information

Coherent Potential Approximation

Coherent Potential Approximation Coheret Potetal Approxato Noveber 29, 2009 Gree-fucto atrces the TB forals I the tght bdg TB pcture the atrx of a Haltoa H s the for H = { H j}, where H j = δ j ε + γ j. 2 Sgle ad double uderles deote

More information

Singular Value Decomposition. Linear Algebra (3) Singular Value Decomposition. SVD and Eigenvectors. Solving LEs with SVD

Singular Value Decomposition. Linear Algebra (3) Singular Value Decomposition. SVD and Eigenvectors. Solving LEs with SVD Sgular Value Decomosto Lear Algera (3) m Cootes Ay m x matrx wth m ca e decomosed as follows Dagoal matrx A UWV m x x Orthogoal colums U U I w1 0 0 w W M M 0 0 x Orthoormal (Pure rotato) VV V V L 0 L 0

More information

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation PGE 30: Formulato ad Soluto Geosystems Egeerg Dr. Balhoff Iterpolato Numercal Methods wth MATLAB, Recktewald, Chapter 0 ad Numercal Methods for Egeers, Chapra ad Caale, 5 th Ed., Part Fve, Chapter 8 ad

More information

An alternative way to write Fermat s equation and notes on an elementary proof of FLT

An alternative way to write Fermat s equation and notes on an elementary proof of FLT A alteratve way to wrte Ferat s equato ad otes o a eleetary roof of FLT C. Sloae Vctora Uversty, NZ (chrssloae70@gal.co) (Feb 08) Abstract We dscovered a way to wrte the equato x +y -z =0 frst studed by

More information

MULTIOBJECTIVE NONLINEAR FRACTIONAL PROGRAMMING PROBLEMS INVOLVING GENERALIZED d - TYPE-I n -SET FUNCTIONS

MULTIOBJECTIVE NONLINEAR FRACTIONAL PROGRAMMING PROBLEMS INVOLVING GENERALIZED d - TYPE-I n -SET FUNCTIONS THE PUBLIHING HOUE PROCEEDING OF THE ROMANIAN ACADEMY, eres A OF THE ROMANIAN ACADEMY Volue 8, Nuber /27,.- MULTIOBJECTIVE NONLINEAR FRACTIONAL PROGRAMMING PROBLEM INVOLVING GENERALIZED d - TYPE-I -ET

More information

Theory study about quarter-wave-stack dielectric mirrors

Theory study about quarter-wave-stack dielectric mirrors Theor tud about quarter-wave-tack delectrc rror Stratfed edu tratted reflected reflected Stratfed edu tratted cdet cdet T T Frt, coder a wave roagato a tratfed edu. A we kow, a arbtrarl olared lae wave

More information

Third handout: On the Gini Index

Third handout: On the Gini Index Thrd hadout: O the dex Corrado, a tala statstca, proposed (, 9, 96) to measure absolute equalt va the mea dfferece whch s defed as ( / ) where refers to the total umber of dvduals socet. Assume that. The

More information

for periodic signals with period Np, it is obvious that

for periodic signals with period Np, it is obvious that Dgtal Seech Processg Lecture 3 Lear Predctve Codg (LPC)- Itroducto LPC Methods LPC ethods are the ost wdely used seech codg, seech sythess, seech recogto, seaer recogto ad verfcato ad for seech storage

More information

Part 4b Asymptotic Results for MRR2 using PRESS. Recall that the PRESS statistic is a special type of cross validation procedure (see Allen (1971))

Part 4b Asymptotic Results for MRR2 using PRESS. Recall that the PRESS statistic is a special type of cross validation procedure (see Allen (1971)) art 4b Asymptotc Results for MRR usg RESS Recall that the RESS statstc s a specal type of cross valdato procedure (see Alle (97)) partcular to the regresso problem ad volves fdg Y $,, the estmate at the

More information

å 1 13 Practice Final Examination Solutions - = CS109 Dec 5, 2018

å 1 13 Practice Final Examination Solutions - = CS109 Dec 5, 2018 Chrs Pech Fal Practce CS09 Dec 5, 08 Practce Fal Examato Solutos. Aswer: 4/5 8/7. There are multle ways to obta ths aswer; here are two: The frst commo method s to sum over all ossbltes for the rak of

More information

Construction of Composite Indices in Presence of Outliers

Construction of Composite Indices in Presence of Outliers Costructo of Coposte dces Presece of Outlers SK Mshra Dept. of Ecoocs North-Easter Hll Uversty Shllog (da). troducto: Oftetes we requre costructg coposte dces by a lear cobato of a uber of dcator varables.

More information

Lecture 3 Probability review (cont d)

Lecture 3 Probability review (cont d) STATS 00: Itroducto to Statstcal Iferece Autum 06 Lecture 3 Probablty revew (cot d) 3. Jot dstrbutos If radom varables X,..., X k are depedet, the ther dstrbuto may be specfed by specfyg the dvdual dstrbuto

More information

CHAPTER 6. d. With success = observation greater than 10, x = # of successes = 4, and

CHAPTER 6. d. With success = observation greater than 10, x = # of successes = 4, and CHAPTR 6 Secto 6.. a. We use the samle mea, to estmate the oulato mea µ. Σ 9.80 µ 8.407 7 ~ 7. b. We use the samle meda, 7 (the mddle observato whe arraged ascedg order. c. We use the samle stadard devato,

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 Aalyss of Varace ad Desg of Exermets-I MODULE II LECTURE - GENERAL LINEAR HYPOTHESIS AND ANALYSIS OF VARIANCE Dr Shalabh Deartmet of Mathematcs ad Statstcs Ida Isttute of Techology Kaur Tukey s rocedure

More information

Introduction to local (nonparametric) density estimation. methods

Introduction to local (nonparametric) density estimation. methods Itroducto to local (oparametrc) desty estmato methods A slecture by Yu Lu for ECE 66 Sprg 014 1. Itroducto Ths slecture troduces two local desty estmato methods whch are Parze desty estmato ad k-earest

More information

Beam Warming Second-Order Upwind Method

Beam Warming Second-Order Upwind Method Beam Warmg Secod-Order Upwd Method Petr Valeta Jauary 6, 015 Ths documet s a part of the assessmet work for the subject 1DRP Dfferetal Equatos o Computer lectured o FNSPE CTU Prague. Abstract Ths documet

More information

Bounds on the expected entropy and KL-divergence of sampled multinomial distributions. Brandon C. Roy

Bounds on the expected entropy and KL-divergence of sampled multinomial distributions. Brandon C. Roy Bouds o the expected etropy ad KL-dvergece of sampled multomal dstrbutos Brado C. Roy bcroy@meda.mt.edu Orgal: May 18, 2011 Revsed: Jue 6, 2011 Abstract Iformato theoretc quattes calculated from a sampled

More information

UNIT 2 SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS

UNIT 2 SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS Numercal Computg -I UNIT SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS Structure Page Nos..0 Itroducto 6. Objectves 7. Ital Approxmato to a Root 7. Bsecto Method 8.. Error Aalyss 9.4 Regula Fals Method

More information

ANALYSIS ON THE NATURE OF THE BASIC EQUATIONS IN SYNERGETIC INTER-REPRESENTATION NETWORK

ANALYSIS ON THE NATURE OF THE BASIC EQUATIONS IN SYNERGETIC INTER-REPRESENTATION NETWORK Far East Joural of Appled Mathematcs Volume, Number, 2008, Pages Ths paper s avalable ole at http://www.pphm.com 2008 Pushpa Publshg House ANALYSIS ON THE NATURE OF THE ASI EQUATIONS IN SYNERGETI INTER-REPRESENTATION

More information

Laboratory I.10 It All Adds Up

Laboratory I.10 It All Adds Up Laboratory I. It All Adds Up Goals The studet wll work wth Rema sums ad evaluate them usg Derve. The studet wll see applcatos of tegrals as accumulatos of chages. The studet wll revew curve fttg sklls.

More information

Ideal multigrades with trigonometric coefficients

Ideal multigrades with trigonometric coefficients Ideal multgrades wth trgoometrc coeffcets Zarathustra Brady December 13, 010 1 The problem A (, k) multgrade s defed as a par of dstct sets of tegers such that (a 1,..., a ; b 1,..., b ) a j = =1 for all

More information

Queueing Networks. γ 3

Queueing Networks. γ 3 Queueg Networks Systes odeled by queueg etworks ca roughly be grouped to four categores. Ope etworks Custoers arrve fro outsde the syste are served ad the depart. Exaple: acket swtched data etwork. γ µ

More information

Several Theorems for the Trace of Self-conjugate Quaternion Matrix

Several Theorems for the Trace of Self-conjugate Quaternion Matrix Moder Aled Scece Setember, 008 Several Theorems for the Trace of Self-cojugate Quatero Matrx Qglog Hu Deartmet of Egeerg Techology Xchag College Xchag, Schua, 6503, Cha E-mal: shjecho@6com Lm Zou(Corresodg

More information

02/15/04 INTERESTING FINITE AND INFINITE PRODUCTS FROM SIMPLE ALGEBRAIC IDENTITIES

02/15/04 INTERESTING FINITE AND INFINITE PRODUCTS FROM SIMPLE ALGEBRAIC IDENTITIES 0/5/04 ITERESTIG FIITE AD IFIITE PRODUCTS FROM SIMPLE ALGEBRAIC IDETITIES Thomas J Osler Mathematcs Departmet Rowa Uversty Glassboro J 0808 Osler@rowaedu Itroducto The dfferece of two squares, y = + y

More information

C.11 Bang-bang Control

C.11 Bang-bang Control Itroucto to Cotrol heory Iclug Optmal Cotrol Nguye a e -.5 C. Bag-bag Cotrol. Itroucto hs chapter eals wth the cotrol wth restrctos: s boue a mght well be possble to have scotutes. o llustrate some of

More information

Lecture Oct. 10, 2013

Lecture Oct. 10, 2013 CS 229r: Algorths for Bg Data Fall 203 Prof. Jela Nelso Lecture Oct. 0, 203 Scrbe: To Morga Overvew I ths lecture ad the revous oe, we focused o fdg Johso-Ldestrauss atrces that allow for faster ultlcato.

More information

CS286.2 Lecture 4: Dinur s Proof of the PCP Theorem

CS286.2 Lecture 4: Dinur s Proof of the PCP Theorem CS86. Lecture 4: Dur s Proof of the PCP Theorem Scrbe: Thom Bohdaowcz Prevously, we have prove a weak verso of the PCP theorem: NP PCP 1,1/ (r = poly, q = O(1)). Wth ths result we have the desred costat

More information

Chapter 4 (Part 1): Non-Parametric Classification (Sections ) Pattern Classification 4.3) Announcements

Chapter 4 (Part 1): Non-Parametric Classification (Sections ) Pattern Classification 4.3) Announcements Aoucemets No-Parametrc Desty Estmato Techques HW assged Most of ths lecture was o the blacboard. These sldes cover the same materal as preseted DHS Bometrcs CSE 90-a Lecture 7 CSE90a Fall 06 CSE90a Fall

More information

IS 709/809: Computational Methods in IS Research. Simple Markovian Queueing Model

IS 709/809: Computational Methods in IS Research. Simple Markovian Queueing Model IS 79/89: Comutatoal Methods IS Research Smle Marova Queueg Model Nrmalya Roy Deartmet of Iformato Systems Uversty of Marylad Baltmore Couty www.umbc.edu Queueg Theory Software QtsPlus software The software

More information

( ) ( ) ( ( )) ( ) ( ) ( ) ( ) ( ) = ( ) ( ) + ( ) ( ) = ( ( )) ( ) + ( ( )) ( ) Review. Second Derivatives for f : y R. Let A be an m n matrix.

( ) ( ) ( ( )) ( ) ( ) ( ) ( ) ( ) = ( ) ( ) + ( ) ( ) = ( ( )) ( ) + ( ( )) ( ) Review. Second Derivatives for f : y R. Let A be an m n matrix. Revew + v, + y = v, + v, + y, + y, Cato! v, + y, + v, + y geeral Let A be a atr Let f,g : Ω R ( ) ( ) R y R Ω R h( ) f ( ) g ( ) ( ) ( ) ( ( )) ( ) dh = f dg + g df A, y y A Ay = = r= c= =, : Ω R he Proof

More information

THE ROYAL STATISTICAL SOCIETY 2016 EXAMINATIONS SOLUTIONS HIGHER CERTIFICATE MODULE 5

THE ROYAL STATISTICAL SOCIETY 2016 EXAMINATIONS SOLUTIONS HIGHER CERTIFICATE MODULE 5 THE ROYAL STATISTICAL SOCIETY 06 EAMINATIONS SOLUTIONS HIGHER CERTIFICATE MODULE 5 The Socety s provdg these solutos to assst cadtes preparg for the examatos 07. The solutos are teded as learg ads ad should

More information

ENGI 3423 Simple Linear Regression Page 12-01

ENGI 3423 Simple Linear Regression Page 12-01 ENGI 343 mple Lear Regresso Page - mple Lear Regresso ometmes a expermet s set up where the expermeter has cotrol over the values of oe or more varables X ad measures the resultg values of aother varable

More information

COMPUTERISED ALGEBRA USED TO CALCULATE X n COST AND SOME COSTS FROM CONVERSIONS OF P-BASE SYSTEM WITH REFERENCES OF P-ADIC NUMBERS FROM

COMPUTERISED ALGEBRA USED TO CALCULATE X n COST AND SOME COSTS FROM CONVERSIONS OF P-BASE SYSTEM WITH REFERENCES OF P-ADIC NUMBERS FROM U.P.B. Sc. Bull., Seres A, Vol. 68, No. 3, 6 COMPUTERISED ALGEBRA USED TO CALCULATE X COST AND SOME COSTS FROM CONVERSIONS OF P-BASE SYSTEM WITH REFERENCES OF P-ADIC NUMBERS FROM Z AND Q C.A. MURESAN Autorul

More information

Introduction to Matrices and Matrix Approach to Simple Linear Regression

Introduction to Matrices and Matrix Approach to Simple Linear Regression Itroducto to Matrces ad Matrx Approach to Smple Lear Regresso Matrces Defto: A matrx s a rectagular array of umbers or symbolc elemets I may applcatos, the rows of a matrx wll represet dvduals cases (people,

More information

best estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best

best estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best Error Aalyss Preamble Wheever a measuremet s made, the result followg from that measuremet s always subject to ucertaty The ucertaty ca be reduced by makg several measuremets of the same quatty or by mprovg

More information

The internal structure of natural numbers, one method for the definition of large prime numbers, and a factorization test

The internal structure of natural numbers, one method for the definition of large prime numbers, and a factorization test Fal verso The teral structure of atural umbers oe method for the defto of large prme umbers ad a factorzato test Emmaul Maousos APM Isttute for the Advacemet of Physcs ad Mathematcs 3 Poulou str. 53 Athes

More information

Ordinary Least Squares Regression. Simple Regression. Algebra and Assumptions.

Ordinary Least Squares Regression. Simple Regression. Algebra and Assumptions. Ordary Least Squares egresso. Smple egresso. Algebra ad Assumptos. I ths part of the course we are gog to study a techque for aalysg the lear relatoshp betwee two varables Y ad X. We have pars of observatos

More information

MEASURES OF DISPERSION

MEASURES OF DISPERSION MEASURES OF DISPERSION Measure of Cetral Tedecy: Measures of Cetral Tedecy ad Dsperso ) Mathematcal Average: a) Arthmetc mea (A.M.) b) Geometrc mea (G.M.) c) Harmoc mea (H.M.) ) Averages of Posto: a) Meda

More information

Polyphase Filters. Section 12.4 Porat

Polyphase Filters. Section 12.4 Porat Polyphase Flters Secto.4 Porat .4 Polyphase Flters Polyphase s a way of dog saplg-rate coverso that leads to very effcet pleetatos. But ore tha that, t leads to very geeral vewpots that are useful buldg

More information