Hydrodynamic Fluid Film Bearings and Their Effect on the Stability of Rotating Machinery

Size: px
Start display at page:

Download "Hydrodynamic Fluid Film Bearings and Their Effect on the Stability of Rotating Machinery"

Transcription

1 Their Effec on he Sabiliy of Roaing Machinery Lui San André Ma-Child Tribology Profeor Turbomachinery Laboraory, Texa A&M Univeriy College Saion, T USA Lanandre@mengr.amu.edu ABSTRACT The lecure inroduce he baic principle of hydrodynamic lubricaion and he fundamenal equaion of Claical Lubricaion Theory. The analyi proceed o derive he aic and dynamic performance characeriic of hor lengh cylindrical journal bearing, wih applicaion o he dynamic forced performance of a rigid roor uppored on plain bearing. In a radial bearing, he Sommerfeld number define a relaionhip beween he aic load and he journal eccenriciy wihin he bearing. Thi deign parameer how he aic performance of he bearing a roor peed increae. Roordynamic force coefficien are inroduced and heir effec on he abiliy of a roor-bearing yem horoughly dicued. Cro-coupled force coefficien are olely due o journal roaion, and he magniude (and ign) of he cro-iffne deermine roordynamic abiliy. The whirl frequency raio (WR) relae he whirl frequency of ubynchronou moion o a hrehold peed of inabiliy. The deired WR i null; however, plain cylindrical bearing how a whirl raio of ju 0.50, limiing he operaion of roaing machinery o haf peed below wice he yem fir criical peed. The analyi conclude wih a review of pracical (in ue) journal bearing configuraion wih highligh on heir major advanage and diadvanage, including remedie o reduce or enirely avoid ubynchronou whirl inabiliy problem. 1.0 UNDAMENTS O LUID ILM BEARING ANALSIS igure 1 depic an idealized geomery of a fluid film bearing. The major characeriic of a lubrican film, and which allow a major implificaion of i analyi, i ha he hickne of he film (h) i very mall when compared o i lengh (L) or o i radiu of curvaure (R), i.e. (h/l) or (h/r) <<< 1. y (U,V) urface velociie V Vy Vz U h(x,z, Vx z x Lz y Lx h << Lx,Lz y V h(x,z, U Vz Vy Vx Lx x igure 1: Geomery of low Region in a Thin luid ilm Bearing (h << Lx, Lz). San André, L. (006) Hydrodynamic luid ilm Bearing and. In Deign and Analyi of High Speed Pump (pp ). Educaional Noe RTO-EN-AVT-14, Paper 10. Neuilly-ur-Seine, rance: RTO. Available from: hp:// RTO-EN-AVT

2 A plain cylindrical journal bearing, ee igure, comprie of an inner roaing cylinder (journal) of radiu R J and an ouer cylinder (bearing) of radiu R B (>R J ). The wo cylinder are cloely paced and he annular gap beween he wo cylinder i filled wih ome lubrican. The radial clearance c (R B -R J ) i very mall. In mo fluid film bearing wih incompreible liquid, c/r B 0.001; while for ga film bearing, C/R B , ypically. D J R J D B R B igure : Schemaic View of a Cylindrical Bearing. A a conequence of he mallne in film hickne, he effec of he film curvaure are negligible in he operaion of a journal bearing. The analyi of he flow equaion define he circumferenial flow Reynold number (Re) a ρ U c Re * & Re* Re (c/l * ) (1a) µ baed on he characeriic peed (U * ΩR). Re denoe he raio beween fluid ineria (advecion) force and vicou-hear force. And, ρ ω c Re (1b) µ i he queeze film Reynold number repreening he raio beween emporal fluid ineria force due o ranien moion a a characeriic frequency (ω) and vicou-hear force. luid ineria effec in hin film flow are of imporance only in hoe applicaion where boh Reynold number are larger han ONE 1, i.e. Re *, Re >> 1. The film hickne o characeriic lengh raio (c/l * ) in hin film flow i ypically very mall. Thu, fluid ineria erm are o be reained for flow wih Reynold number of he order []: Re > (L * /c) x 1 > 1, 000 for (c/l * ) () Claical lubricaion i baed on he aumpion ha fluid ineria effec are negligible, i.e. (Re, Re * ) 0, in mo pracical applicaion; and hence, rendering effecively an ineriale fluid. 1 In acualiy, Re > 1 for eady uper laminar flow in hin film bearing, a demonraed in [1]. In hin film flow, raniion o urbulence i due o inabiliy of hear driven parallel flow. A raniion o urbulence iniiaed by he appearance of Taylor vorice generaed by cenrifugal force i more peculiar o configuraion wih large clearance, i.e. no common in hin film bearing. The acceped Reynold number for flow urbulence in journal bearing i Re c ρωrc/µ >,000 [] RTO-EN-AVT-14

3 The hin film laminar flow of an incompreible, ineriale and iovicou fluid i governed by he following equaion coninuiy and momenum ranpor equaion: ( v ) ( v y ) ( v ) x x + y + z z 0 P x v y P z v y x x ; 0 + ; 0 + µ µ () where v x,y,z are he velociy componen and he film preure (P) i uniform acro he film hickne (h). The momenum equaion eablih a quai-aic balance of preure and vicou force. The lubricaion equaion denoe low flow condiion wih ime appearing a a parameer, no an independen variable. Table 1 preen he circumferenial flow Reynold number (Re) for a ypical journal bearing applicaion operaing wih differen fluid. The example bearing i a inch (76 mm) diameer (R J ) journal and he clearance o radiu raio (c/r J ) i 0.001, a ypical value for journal bearing. Two roaional peed of 1,000 and 10,000 rpm (Ω104.7 and 1047 rad/) are noed in he able. The calculaed Reynold number (Re) how ha bearing applicaion wih mineral oil and (even) air do no need o include fluid ineria effec, i.e. Re < 1,000. However, proce fluid applicaion uing waer, R14a refrigeran and cryogenic fluid how large Reynold number a a peed of 10,000 rpm. luid Table 1: Imporance of luid Ineria Effec on Several luid ilm Bearing Applicaion (c/r J )0.001, R J 8.1 mm (1.5 inch) Abolue vicoiy (µ) lbm.f. x 10-5 Kinemaic vicoiy (ν) cenioke Re a 1,000 rpm Re a 10,000 rpm Air Thick oil 1, Ligh oil Waer ,588 Liquid hydrogen ,05 Liquid oxygen ,94 Liquid nirogen ,477 R14 refrigeran ,96 Noe ha curren bearing applicaion uing proce liquid o replace mineral oil may operae a peed well above 10,000 rpm. Incidenally, he operaing peed of cryogenic urbopump i on he order of 0 70 krpm, and fuure applicaion (currenly in he work) will operae a peed cloe o 00 krpm! Incidenally, proce ga and liquid eal, iolaing region of high and low preure in a ypical compreor or pump, have larger radial clearance han load uppor fluid film bearing. or example, in waer neck-ring and inerage eal in pump, R/C ~ 50, and hu fluid ineria effec are of imporance even a relaively low roaional peed (~1,000 rpm and larger). The opic of eal i analyzed in he nex lecure. 1.1 Oher luid Ineria Effec luid ineria wihin hin film flow domain can be afely ignored in mo convenional (oil lubricaed) bearing applicaion. However, fluid ineria effec may alo be of grea imporance a he inle o he film and dicharge from he film ecion in a ypical pad bearing or eal applicaion, ee igure. Depending on he flow condiion upream of a udden conracion or a udden enlargemen, a fracion of he dynamic preure head, ypically given a (½ ρ U ), i lo or recovered. RTO-EN-AVT

4 U U P P ~ ½ ρu P igure : Preure Drop & Rie a Sudden Change in ilm Thickne. Sudden preure loe are ypical a he edge of a pocke in a hydroaic bearing and a he inle plane of annular preure eal. The ame phenomenon alo occur a he leading edge of a bearing pad in high peed iling pad bearing. A udden preure recovery i alo quie ypical a he dicharge ecion of a preurized annular or labyrinh eal. Noe ha he imporance of fluid ineria effec may be rericed only o he inle and dicharge ecion, and may no be relevan wihin he hin film flow domain..0 RENOLDS EQUATION AND KINEMATICS O JOURNAL MOTION Lubricaed cylindrical bearing are low fricion, load bearing uppor in roaing machinery. Thee fluid film bearing alo inroduce vicou damping ha aid in reducing he ampliude of vibraion in operaing machinery. igure 4 how a chemaic view of a cylindrical bearing. The journal pin wih angular peed (Ω) and i cener (O J ), due o dynamic load, alo decribe ranlaional moion wihin he bearing clearance. The bearing or houing i aionary in mo applicaion. Noable excepion are hoe of floaing ring journal bearing and auomoive reciprocaing engine uppor rod bearing RTO-EN-AVT-14

5 xrθ Noe: film gap enlarged for decripion purpoe y A Θ θ h Ω e OB OJ Journal φ Bearing r Θθ+φ Clearance CR B -R J OB e ilm hickne: h C + e co(θ) or h C + e co(θ) + e in(θ) e e OJ φ r igure 4: Schemaic View of a Cylindrical Journal Bearing. ixed Coordinae Syem (θ,z) and Moving Coordinae Syem (Θ,z)..1 Reynold Equaion for Journal Bearing The mallne of hi raio allow for a Careian coordinae (xrθ, y, z) o be locaed on he bearing urface (ee igure 4). Then, in Claical Lubricaion, Reynold equaion decribe he generaion of hydrodynamic preure (P) wihin he bearing. Thi equaion arie from inegraion of he momenum equaion () acro he film hickne and ubiuion ino he coninuiy equaion []: Ω 1 ρ h P ρ h P { ρ h} + { ρ h} + (4) Θ R Θ 1 µ Θ z 1 µ z in he flow domain {0 Θ π, -½ L z ½ L}, where h(θ, z, ) i he film hickne, L i he bearing axial lengh, U ΩR J i he journal urface peed, and (ρ, µ) denoe he lubrican deniy and vicoiy, repecively. RTO-EN-AVT

6 The boundary condiion for he hydrodynamic preure in a plain cylindrical bearing are : a) he preure i coninuou and periodic in he circumferenial direcion, i.e. P(Θ, z,) P(Θ + π, z, ) (5) b) he preure equal he dicharge or amopheric value (P a ) on he bearing ide, i.e. P(Θ, ½ L, ) P(Θ, - ½ L, ) P a (6) A a conrain, he hydrodynamic preure need o be greaer han he liquid caviaion preure everywhere in he flow domain, i.e. P P cav in 0 Θ π, - ½ L z ½ L (7) Here P cav repreen he lubrican auraion preure or he ambien preure needed for releae of diolved gae. In pracice, no diincion i made beween hee wo value ince hydrodynamic film preure could be one o wo order of magniude larger han he ambien value. Conider he journal and bearing o be aligned and he journal cener o have an eccenriciy diplacemen e ( c). The film hickne i h c + eco(θ ) (8) Thi formula i accurae for (c/r) raio a large a The film hickne derived aume rigid bearing and journal urface, uniform axial and azimuhal clearance and no journal mialignmen.. Kinemaic of Journal Moion The journal cener O J i diplaced a diance (e) from he bearing cener O B. Thi diance i known a he journal eccenriciy and may vary wih ime depending upon he impoed exernal load on he bearing. The journal eccenriciy canno exceed he bearing clearance, oherwie olid conac and poenial caarophic failure may occur. The eccenriciy componen in he (, ) fixed coordinae yem are: e e co(φ); e e in(φ) (9) where φ i known a he bearing aiude angle, and Θθ+φ. Then, he film hickne alo equal h c + e coθ + e in Θ e inθ (10) and h Θ e in Θ + e coθ; h e coθ + e in Θ (11) where (. ) denoe differeniaion wih repec o ime, i.e. ( / ). Subiuion of he film hickne gradien ino Reynold equaion (4) give he following lubricaion equaion for an incompreible and iovicou fluid: The imple journal bearing model doe no accoun for feeding hole or axial groove for upply of he lubrican ino he bearing. A more deailed dicuion on lubrican caviaion and i phyical model can be found in [] RTO-EN-AVT-14

7 1 R h P h P + e Θ 1µ Θ z 1µ z + e Ω co Θ + e Ω e in Θ (1) An alernaive form of Reynold equaion arie when uing he angular coordinae (θ). Thi angle ar from he locaion of maximum film hickne. A coordinae yem wih radial and angenial (r, ) axe i convenienly defined wih he uni radial vecor along he line joining he bearing and journal cener. Recall ha e e co(φ); e e in(φ), and e e + e. The journal cener velociie in he (, ) and (r, ) coordinae yem are relaed by he ranformaion: OB e V e e coφ inφ inφ e coφ eφ (1) OJ φ r V r Noe ha V e, V e φ are he radial and angenial r componen of he journal ranlaional velociy. rom he film hickne h c + e coθ, i follow Tranlaional velociie of journal cener h h einθ θ θ e coθ e inθ e coθ + e φ inθ V r coθ + V inθ (14) Thu, Reynold equaion (4) for an incompreible and iovicou fluid i alo expreed a 1 R h P h P Ω + e coθ + e φ inθ θ 1µ θ z 1µ z (15) Equaion (15) i of paricular imporance ince i allow u o realize an imporan phyical phenomenon. Conider he journal cener o decribe circular cenered orbi wih a fixed ampliude or radiu, e. Hence de/d0. urhermore, if he frequency of whirl equal o 50% of he roaional peed; φ Ω / ; hen he righ hand ide of Eqn. (15) i null; and hence he preure i zero, P0 wihin he bearing film land. There i no generaion of hydrodynamic preure, hu reuling in a udden lo of load uppor capabiliy. The 50% ped whirl phenomenon i he bai of roordynamic inabiliy, a explained laer.. Bearing Reacion orce Once he preure field i obained, fluid film force acing on he journal urface, ee igure 5, are calculaed by inegraion of he preure field acing on he journal urface. An equal oppoing force ac on he bearing a well. The bearing reacion force are expreed in he fixed (, ) coordinae yem and moving (r, ) coordinae yem a RTO-EN-AVT

8 10-8 RTO-EN-AVT-14 ( ) dz d R z P L Θ Θ Θ Θ in co,, 0 0 π ; ( ) dz d R z P L r θ θ θ θ π in co,, 0 0 (16) θ P.coθ P.inθ P r θ Θ P journal r igure 5: luid ilm orce Acing on Journal Surface. The relaionhip beween he fluid film force in boh coordinae yem i given by: r φ φ φ φ co in in co (16.b) The fluid film force are generic funcion of he journal roaional peed (Ω) and he journal cener ranlaional velociie, i.e. ( ) Ω Ω,,, φ α α α e e e e ; α, or r, (17) An analyical oluion of Reynold equaion for arbirary geomery cylindrical bearing i no feaible. Mo frequenly, numerical mehod are employed o olve Reynold equaion and o obain he performance characeriic of bearing configuraion of paricular inere. There are analyical oluion o Reynold equaion applicable o wo limiing geomerie of journal bearing. Thee are known a he infiniely long and infiniely hor lengh journal bearing model []. In he LONG BEARING MODEL, ee igure 6, he lengh of he bearing i very large, L/D, and conequenly he axial flow i effecively very mall, i.e. ( P/ z) 0.

9 L bearing D journal Ω Axial preure field igure 6: The Long Bearing Model. or large L/D raio, Reynold equaion reduce o: 1 R h P θ 1 µ θ Ω Θ { h} + { h} (18) Thi bearing model give accurae reul for journal bearing wih lenderne raio (L/D) >. Mo modern bearing in high performance urbomachinery applicaion have a mall L/D raio, rarely exceeding one. Thu, he infiniely long journal bearing model i of limied curren inere. Refer o [] for deail on he analyical oluion of Eqn. (18). Thi i no he cae for queeze film damper (SD), however, ince he long bearing model provide a very good approximaion for ighly ealed damper even for mall L/D raio [4]..0 STATIC LOAD PERORMANCE O SHORT LENGTH BEARINGS In hi mo ueful bearing model, ee igure 7, he bearing lengh i hor, L/D 0, and conequenly he circumferenial flow i effecively mall, i.e. ( P/ θ) 0. or hi limiing bearing configuraion, Reynold equaion reduce o Ω h P { h} + { h} Θ z 1 µ z (19) L bearing D journal Ω L/D << 1 dp/dθ 0 Axial preure field igure 7: Shor Lengh Bearing Model. RTO-EN-AVT

10 The hor lengh bearing model provide (urpriingly) accurae reul for plain cylindrical bearing of lenderne raio L/D 0.50 and for mall o moderae value of he journal eccenriciy, e 0.75 c [4]. The hor lengh bearing model i widely ued for quick eimaion of journal bearing aic and dynamic force performance characeriic. Inegraion of equaion (19) lead o he preure diribuion P( θ, z, ) P a Ω 6 µ e coθ + e φ inθ z C H L (0) wih Hh/c 1 + ε co(θ) a he dimenionle film hickne, and ε e/c i he journal eccenriciy raio; [0 ε 1], ε 0 mean cenered operaion (ypically a condiion of no load uppor), and ε 1.0 evidence olid conac of he journal wih i bearing. No lubrican caviaion will occur if he exi or dicharge preure P a i well above he liquid caviaion preure. However, if P a i low, ypically ambien condiion a 1 bar, i i almo cerain ha he bearing will caviae or how air enrainmen when he oule plenum i no flooded wih lubrican. The caviaion model in he hor lengh bearing imply neglec any prediced negaive preure and equae hem o zero. Thi chop procedure alhough heoreically no well juified grap wih ome accuracy he acual phyic [5]. Hence, if P a 0, and from equaion (1), he preure field P>0, when co(θ+α) < 0. Thu, P>0 in he circumferenial region limied by π π θ + α π π α θ1 θ θ α (1) Tha i, regardle of he ype of journal moion, he region of poiive preure ha an exen of π (180 ); hu hen he infamou π film caviaion model widely ued in he lieraure. luid film reacion force on he journal are evaluaed by inegraion of he preure field acing on he journal urface. Wih P a 0, he radial and angenial force ( r, ) are given by r µ R L c J J 0 11 J J 11 0 e Ω e φ () where he J are inegral defined in analyical form by Booker [6]. Noe ha he fluid film force are proporional o he journal cener ranlaional velociie ( e, e φ ) a well a he journal roaional peed (Ω). The reacion force depend linearly on he fluid vicoiy and he bearing radiu and grow rapidly wih he raio (L/C). Hydrodynamic journal bearing are deigned (and implemened) o uppor a aic load W, hereafer aligned wih he axi for convenience, ee igure 8. A he equilibrium condiion, denoed by a journal cener eccenric diplacemen (e) wih an aiude angle (φ), he hydrodynamic bearing generae a reacion force balancing he applied exernal load a he raed roaional peed (Ω). The equaion of aic equilibrium are W W 0 coφ r r inφ + inφ coφ () RTO-EN-AVT-14

11 bearing Saic load W Journal Roaion Ω e W - r r fluid film φ Roor (journal) φ igure 8: orce Equilibrium for Saically Applied Load. or aic equilibrium, e 0, φ 0,and θ 1 0 o θ π. rom equaion (), he aic radial and angenial film reacion force are µ R L Ω ε µ R L Ω π ε r ; + (4) c ( 1 ε ) c 4( 1 ε ) / igure 9 depic he radial and angenial force for a ypical hor lengh bearing. The force are proporional o he lubrican vicoiy and roor urface peed (ΩR), he lengh (L ), and inverely proporional o he radial clearance (c ). Mo imporanly, he bearing force grow rapidly (non-linearly) wih he journal eccenriciy (εe/c) Saic orce for hor lengh bearing Radial and Tangenial force [N] r journal eccenriciy (e/c) * igure 9: Radial and Tangenial orce for Shor Lengh Bearing. µ0.019 Pa., L0.05 m, c0.1 mm,, 000 rpm, L/D0.5. The exernal load (W) i balanced by he fluid film reacion force. Thu, W ( + ) + π ( 1 ε ) ( 1 ) 1 L ε 16ε r µ Ω R L (5) c 4 ε RTO-EN-AVT

12 and he journal aiude angle φ i obained from angφ r π ( 1 ε ) 4 ε (6) Noe ha a he journal eccenriciy ε 0, φ π /, while a ε 1, φ 0..1 Deign of Hydrodynamic Bearing Selecion of Operaing Eccenriciy In he deign of hydrodynamic journal bearing, he bearing aic performance characeriic are relaed o a unique dimenionle parameer known a he Sommerfeld Number (S) defined a S µ N L D W R c (7) where N (Ω/π) i he roaional peed in revoluion/ec. In pracice, he pecific load or preure i known a he raio of applied load o bearing projeced area, i.e. (W/LD). In hor lengh journal bearing, a modified Sommerfeld number (σ) i defined and relaed o (S) by [5, 7]: µ Ω L R L σ π S ( L D) (8) 4W C Subiuion of Eqn. (8) ino Eqn. (5) relae he modified Sommerfeld number o he equilibrium operaing journal eccenriciy (e), i.e. ( 1 ε ) { 16ε + π ( 1 )} µ Ω L R L σ 4W c ε ε A a raed operaing condiion, σ i known ince he bearing geomery (R, L, c), roaional peed (Ω), fluid vicoiy (µ) and applied load (W) are known. Then, equaion (9) provide a relaionhip o deermine (ieraively) he equilibrium journal eccenriciy raio ε (e/c) required o generae he fluid film reacion force balancing he exernally applied load W. igure 10 and 11 depic he modified Sommerfeld number and aiude angle v. journal eccenriciy, repecively. Large Sommerfeld (σ) number; i.e. denoing mall load, high peed Ω or large lubrican vicoiy, deermine mall operaing journal eccenriciie or nearly cenered operaion, ε 0, φ π/ (90 ). Tha i, he journal eccenriciy vecor i nearly orhogonal o he applied load. (9) 10-1 RTO-EN-AVT-14

13 high peed mall load, large vicoiy Sommerfeld number for hor journal bearing Sommerfeld number Modified Sommerfeld number journal eccenriciy (e/c) low peed large load, low vicoiy igure 10: Modified Sommerfeld (σ) Number veru Journal Eccenriciy high peed mall load, large vicoiy aiude angle journal eccenriciy (e/ c) low peed large load, low vicoiy igure 11: Equilibrium Aiude Angle veru Journal Eccenriciy. Small Sommerfeld (σ) number, i.e. denoing large load, low peed Ω or low lubrican vicoiy, deermine large operaing journal eccenriciie, ε 1.0, φ 0 (0 ). Noe ha he journal eccenriciy vecor i nearly parallel o he applied load. igure 1 how he journal diplacemen wihin he bearing clearance for differen operaing condiion. The journal eccenriciy approache he clearance for large load, low haf peed or ligh lubrican vicoiy, and i i aligned wih he load vecor. or mall load, high peed or large lubrican vicoiie RTO-EN-AVT

14 (large Sommerfeld number), he journal ravel oward he bearing cener and i poiion i orhogonal o he applied load. Thi peculiar behavior i he ource of roordynamic inabiliy a will be hown horly. ey/c Low load, high peed, large vicoiy W load e/c Journal locu Clearance circle ex/c aiude angle load increae, low peed, low vicoiy pin direcion clearance circle igure 1: Locu of Journal Cener for Shor Lengh Bearing. 4.0 DNAMICS O A RIGID ROTOR SUPPORTED ON SHORT LENGTH BEARINGS igure 1 depic a ymmeric rigid roor of ma M, and upporing a aic load ( o W) along he axi. The roor i mouned on wo idenical plain hydrodynamic journal bearing. The equaion of moion of he roaing yem a conan roaional peed (Ω) are given by [5]: M + M u Ω in( Ω) + o M + M u Ω co( Ω) (0) where u i he magniude of he imbalance vecor, () and () are he coordinae of he roor ma cener, and (, )are he fluid film bearing reacion force. Since he roor i rigid, he cener of ma diplacemen are idenical o hoe of he journal bearing, i.e. ( ) e ( ), ( ) e ( ) RTO-EN-AVT-14

15 Saic load Rigid haf Dik M Journal bearing dik o u Ω Clearance circle e igure 1: Rigid Roor Suppored on Journal Bearing. (u) Imbalance, (e) Journal Eccenriciy. We are inereed on he roor dynamic behaviour for mall ampliude moion abou he equilibrium poiion defined by: O, 0, e, e or e, φ (1) o O O O O O where (e o,φ o ) denoe he aic equilibrium journal eccenriciy and aiude angle, repecively. The bearing aic reacion force aify o - ro o φ o r o + O O 0 0 o 0 O O ro ro coφ O O inφ + O O coφ inφ Small ampliude journal moion abou he equilibrium poiion, a repreened in igure 14, are defined a: O O () e e + e ( ), e e + e ( ), or + ( ), + ( ) (a) O O O O or converely, e( ) e + e( ), φ( ) φ + φ( ) (b) O O wih d e d d e d e e.. d e e d d e e d (c) RTO-EN-AVT

16 φ o o e e e o e clearance circle e e e φ r igure 14: Small Ampliude Journal Moion abou an Equilibrium Poiion. The journal dynamic diplacemen in he (r, ) coordinae yem are relaed o hoe in he (, ) yem by he linear ranformaion e e coφo inφo inφo e( ) coφ O eo φ( ) (4) Similar relaionhip hold for he journal cener velociie and acceleraion. Noe ha he aumpion of mall ampliude moion require e, e << c, i.e., he journal dynamic diplacemen are maller han he bearing clearance. Recall ha he fluid film force are general funcion of he journal cener diplacemen and velociie, i.e. [ e ( ), e ( ), e ( ), e ( )], α. Now, expre he bearing reacion force a a α α, Taylor Serie expanion around he aic journal poiion (e o, e o ), i.e. O O (5) 4.1 Definiion of Dynamic orce Coefficien in luid ilm Bearing luid film bearing iffne (K ij ) ij, and damping (C ij ) ij, force coefficien are defined a K ij i j ; C ij i i, j, (6) j or example, K - / correpond o a iffne produced by a fluid force in he direcion due o a journal aic diplacemen in he direcion. By definiion, hi coefficien i evaluaed a he equilibrium RTO-EN-AVT-14

17 poiion wih oher journal cener diplacemen and velociie e o zero. The negaive ign in he definiion enure ha a poiive magniude iffne coefficien correpond o a reoraive force. The force coefficien (K, K ) are known a he direc iffne erm, while (K, K ) are referred a cro-coupled. igure 15 provide a picorial repreenaion of he bearing force coefficien a mechanical parameer. Kxx, Cxx Kyy Cyy Kxy, Cxy K ij - i / j journal C ij - i /( j /δ) Kyx, Cyx bearing igure 15: The Phyical Repreenaion of Siffne and Damping Coefficien in Lubricaed Bearing. Ineria or added ma coefficien {M ij } ij, can alo be defined a i M ij ;, are j journal cener acceleraion. Ineria coefficien are of paricular imporance in uper laminar and urbulen flow fluid film bearing and annular eal. The ineria force coefficien or apparen mae have a ound phyical inerpreaion. Thee coefficien are alway preen in a fluid film bearing. Ineria coefficien can be of large magniude, in paricular for dene liquid. However, he effec of ineria force on he dynamic repone of roor-bearing yem i only of imporance a large exciaion frequencie, i.e. high queeze film Reynold number. (Thi fac alo hold for mo mechanical yem ubjeced o fa ranien moion). Wih he given definiion, he bearing reacion force are repreened a i,j, where { } ( ) ( ) O O K K K K C C C C (7) where o o ½W and o 0. Noe ha he defined force coefficien allow he repreenaion of he dynamic fluid film bearing (or eal) force in erm of he fundamenal mechanical parameer {K, C, and M}. However, hi doe no mean ha hee force coefficien mu be accordance wih acceped phyic grounded knowledge. or example, he vicou damping coefficien may be negaive, i.e. nondiipaive, or he iffne coefficien non-reoraive. rom Eqn. (0), he linear equaion for mall ampliude moion of he roor-bearing yem are M O O C + M C C C K + K K K M u Ω coω in Ω (8) RTO-EN-AVT

18 The lieraure preen he force coefficien in dimenionle form according o he definiion: k ij c cω Kij ; cij Cij i,j, (9) 0 0 where o i he aic load applied on each bearing (along he direcion). [Noe ha he oal aic load W o i hared by he wo bearing in a ymmeric roor moun]. Lund [8] derived fir he analyical formula for he hor bearing force coefficien. igure 16 and 17 depic he dimenionle force coefficien, iffne and damping, a funcion of he journal eccenriciy and of he modified Sommerfeld number (σ), repecively. In he figure, boh repreenaion are neceary ince a ime he journal eccenriciy i known a priori; while mo ofen, he deign parameer, i.e. he Sommerfeld number, i known in advance. In general, he phyical magniude of he iffne and damping coefficien increae rapidly (nonlinearly) a he journal eccenriciy increae RTO-EN-AVT-14

19 100 Siffne coefficien x (c/o) Siffne K (c/o) kxx kyy kyx kxy keq journal eccenriciy (e/c) high peed low peed mall load < > large load large vicoiy mall vicoiy 100 Damping coefficien x (c Omega/o) Damping: C (c Omega/o) cxx cyy cyx cxy Modified Sommerfeld number igure 16: Dimenionle Siffne and Damping Coefficien v. Journal Eccenriciy (ε) for Shor Journal Bearing. RTO-EN-AVT

20 100 Siffne coefficien x (C/o) Siffne K (C/o) kxx kyy kyx kxy keq Modified Sommerfeld number low peed high peed large load < > mall load low vicoiy large vicoiy 100 Damping coefficien x (c Omega/o) Damping: C (c Omega/o) cxx cyy cyx cxy Modified Sommerfeld number igure 17: Dimenionle Siffne and Damping Coefficien v. Sommerfeld Number (σ) for Shor Journal Bearing RTO-EN-AVT-14

21 Noe ha he dimenionle force coefficien do no repreen he acual phyical rend. or example, a e o 0, K K 0, bu he dimenionle k k have non zero value. Thi peculiariy follow from he definiion of dimenionle force coefficien uing he applied load ( o ). Recall ha, a e o 0, he aic load o i alo naugh. 4. Dynamic orce Coefficien for Journal Cenered Operaion, i.e. No Applied Load A he journal cener approache he bearing cener, e o 0, φ o 90 o, and from he formula preened, K µ Ω R L π Ω µ R L π K k c; C C c (40) c 4 c Thu, a he cenered journal poiion a hydrodynamic bearing offer no direc (uppor) iffne bu only cro-coupled force. A mall load applied on he bearing will caue a journal diplacemen in a direcion orhogonal (ranvere) o he load, a hown in he chemaic view below. Thi behaviour i common o all fluid film journal bearing of rigid geomery. Non-roaing rucure Roaing rucure Ω The ignificance of he cro-coupled effec in fluid film bearing 5.0 ROTORDNAMIC STABILIT O RIGID ROTOR SUPPORTED ON SHORT LENGTH BEARINGS The linearized equaion of moion are wrien in dimenionle form a [5, 9] p x c c x k k x + + co( τ ) p δ (41) y c c y k k y in( τ ) where u x, y, τ Ω, δ, d C M ( ' ) ; p c c c dτ c ij are he dimenionle iffne and damping force coefficien. Ω o i a dimenionle ma, and k ij and I i of inere o deermine if he roor-bearing yem i able for mall ampliude journal cener moion (perurbaion) abou he equilibrium poiion. To hi end, e he imbalance parameer δ 0 in he equaion above o obain, RTO-EN-AVT

22 p x c c x k k x (4) y c c y k k y 0 If he roor-bearing yem i o become unable, hi will occur a a hrehold peed of roaion (Ω ) and he roor will perform (undamped) orbial moion a a whirl frequency (ω ). Thee moion, aifying equaion (4), are of he form: jω jωτ jω jωτ x Ae Ae ; y B e B e ; j 1 (4) where ω ω Ω i known a he whirl frequency raio, i.e. he raio of he roor whirl or preceional frequency o he roor one peed of inabiliy. Subiuion of (4) ino equaion (4) lead o [5]: p ω k + k + + jω c jω c k p ω + + k jω c + jω c A 0 B 0 (44) In Eqn. (44), he deerminan, mu be zero for a non-rivial oluion of he homogenou yem. Afer algebraic manipulaion, he real and imaginary par of render [5,9] p C M (45) k c + k c c k c k ω S ω keq c + c o and ( k k )( k k ) k k eq eq ω ω (46) c c c c Ω 5.1 Threhold Speed, Criical Ma, Equivalen Siffne and Whirl requency Raio or a given value of journal eccenriciy (ε o ), i.e. a given Sommerfeld number (σ ), one evaluae Eqn. (45) o obain he equivalen iffne k eq, and hen Eqn. (46) o ge he whirl frequency raioω S. Thi ubiuion hen yield p ω (~criical ma), which in urn render he one peed of inabiliy Ω. k eq S igure 18 and 19 depic he whirl frequency raio (ω/ω) and he dimenionle hrehold peed of inabiliy (p ) veru equilibrium journal eccenriciy and modified Sommerfeld number, repecively. Noe ha for near cenered journal operaion, i.e. large Sommerfeld number, he whirl frequency i 0.50, i.e. half-ynchronou whirl. Recall ha in a mechanical yem, an equivalen damping raio > 0 caue he aenuaion of moion induced by mall perurbaion from an equilibrium poiion. A null damping raio bring he yem ino uained periodic moion wihou decay or growh, hu denoing he hrehold beween abiliy and inabiliy (ampliude growing moion) RTO-EN-AVT-14

23 1 Whirl frequency raio Whirl frequency raio Modified Sommerfeld number 1 Whirl frequency raio Whirl frequency raio journal eccenriciy (raio) high peed low peed mall load < > large load large vicoiy mall vicoiy igure 18: Whirl requency Raio v. Sommerfeld Number (σ) and Journal Eccenriciy (ε). RTO-EN-AVT

24 Dim [-] hrehold peed (p) 6 (p) hrehold peed inabiliy 4 unable Modified Sommerfeld number Dim [-] hrehold peed (p) (p) hrehold peed inabiliy 6 4 unable journal eccenriciy (e/c) high peed low peed mall load < > large load large vicoiy mall vicoiy igure 19: Threhold Speed of Inabiliy (p ) v. Sommerfeld Number (σ) and Journal Eccenriciy (ε). On he oher hand, if one aume ha he curren roaional peed (Ω) i he one peed of inabiliy, hen from he relaion above i follow he large magniude of ½ yem ma (M) o make he roorbearing yem unable. Thi ma i known a he criical ma, M c, and correpond o he limi ma which he yem can carry dynamically. If he roor ma i equal o or larger han wice M c, hen he yem will become unable a he raed peed Ω 4. The whirl frequency raio (WR), ω Ω, i he raio of he roor whirl frequency o he one peed of inabiliy. Noe ha hi raio, a given in equaion (46), depend only on he fluid film bearing 4 Recall ha each bearing carrie half he aic load, and alo half he dynamic or ineria load (.M c C Ω ) RTO-EN-AVT-14

25 characeriic and he equilibrium eccenriciy. The WR i independen of he roor characeriic (roor ma and flexibiliy) [5]. Reference [10] preen an analyi including fluid ineria effec, more applicable o annular preure eal and bearing handling proce fluid of large deniy. The parameer k eq i a journal bearing (dimenionle) equivalen iffne, alo depiced in igure 16 and 17. rom he definiion of hrehold peed and whirl raio, p M Ω ( C ) andω ω Ω, hen o Thu, he whirl or preceional frequency i ω K eq (47) C o M keq K eq ω ωn (48) M i.e., he whirl frequency equal he naural frequency of he rigid roor uppored on journal bearing. or operaion cloe o he concenric poiion, ε o 0, i.e. large Sommerfeld number (no load condiion), he force coefficien are, ee equaion (40), k k 0 ; c c ; k k ; c c 0 (49) eq ( k c c k ) c k + 0 ω k and 0.50 a ε 0 Ω c (50) Thi value of whirl frequency raio (WR) i a characeriic of hydrodynamic plain journal bearing. The WR how ha a he one ped of inabiliy he roor whirl a i naural frequency equal o 50% of he hrehold roaional peed. urhermore, under no exernally applied load, o 0, a in verically urbomachinery, he bearing poee no uppor iffne, i.e. K eq 0 and he yem naural frequency (ω n ) i zero, i.e. he roor-bearing yem mu whirl a all operaing peed. Noe ha if K 0, i.e. he bearing doe no have cro-coupled effec, hen he WR 0, i.e. no whirl occur and he yem i alway dynamically able. Cro-coupled effec are hen reponible for he inabiliie o commonly oberved in roor mouned on journal bearing. If he whirl frequency raio i 0.50, hen he maximum roaional peed ha he roor-bearing yem can aain i ju, i.e., wice he naural frequency (or oberved rigid roor criical peed). ω Ω max ω ω n (51) 0.50 igure 18, 19 and 0 depic he whirl frequency raio, he dimenionle hrehold peed (p ) and he criical ma (p ) veru he Sommerfeld number and equilibrium journal eccenriciy. The reul demonrae ha a rigid-roor uppored on plain journal bearing i STABLE for journal eccenriciy raio ε > 0.75 (mall Sommerfeld number) for all L/D raio. Noe ha he criical ma and whirl frequency raio are nearly invarian for operaion wih journal eccenriciie (ε o ) below RTO-EN-AVT

26 0 Criical ma Dimenionle Criical Ma unable Modified Sommerfeld number 0 Dimenionle Criical ma Dimenionle Criical Ma unable journal eccenriciy (e/c) high peed low peed mall load < > large load large vicoiy mall vicoiy igure 0: Criical Ma (m c p ) v. Sommerfeld Number (σ) and Journal Eccenriciy (ε). Keep in mind ha increaing he roaional peed of he roor-bearing yem deermine larger Sommerfeld number, and conequenly, operaion a maller journal eccenriciie for he ame applied aic load. Thu, operaion a ever increaing peed will evenually lead o a roor dynamically unable yem a he analyi reul how. 5. Effec of Roor lexibiliy on Sabiliy of Syem A imilar analyi can be performed conidering roor flexibiliy [5, 11]. Thi analyi i more laboriou hough raighforward. The analyi how ha roor flexibiliy doe no affec he whirl frequency raio. However, he one peed of inabiliy i dramaically reduced ince he naural frequency of he roorbearing yem i much lower. The relaionhip for he hrehold peed of inabiliy of a flexible roor i: 10-6 RTO-EN-AVT-14

27 K ro M bearing p f p T 1+ keq C (5) where he ub index f denoe a flexible roor. K ro i he roor iffne on each ide of he mid dik hown in he graph, and T o K i he roor aic ag or elaic deformaion a midpan. ro The elaic haf and bearing are mouned in erie, i.e. he bearing and haf flexibiliie add (reciprocal of iffnee), and hu he equivalen yem iffne i lower han ha of he bearing, and herefore he yem naural frequency decreae ignificanly. igure 1 depic he hrehold peed of inabiliy (p f ) for a flexible roor mouned on plain hor lengh journal bearing. Noe ha he more flexible he roor i, he lower he hrehold peed of inabiliy. If he fluid film bearing are deigned oo iff (low Sommerfeld number), hen he naural frequency of he roor-bearing yem i ju (K ro /M) 1/, irrepecive of he bearing configuraion. Threhold peed (p) for flexible roor 6 Threhold peed (p) 4 unable rigid T/c0.1 T/c1 T/c10 Modified Sommerfeld number low peed high peed large load < > mall load low vicoiy large vicoiy igure 1: Threhold Speed of Inabiliy (p ) for lexible Roor veru Sommerfeld Number (σ). Saic Sag (Γ/c) Varie. 5. Phyical Inerpreaion of Dynamic orce for Circular Cenered Whirl The bearing dynamic force in he radial and angenial are RTO-EN-AVT

28 r d K K rr r K K r e C e φ C rr C C 0 r e 0 r e φ (5) Conider circular journal moion of ampliude e a a forward frequency (ω), a hown in igure. A he cenered poiion, he bearing ha no direc iffnee, only cro-coupled iffne and direc damping, i.e. K rr K C r C r 0 (54) K K r K r Ω C ; C C C rr µ R L C π whirl orbi (ω) -(C ω + K r ) e Roor pin e r -(C r ω + K rr ) e igure : orce Diagram for Circular Cenered Whirl Moion. And hu, he radial and angenial force become r d 0; ( C ω K ) e (55) d r A deabilizing force will drive he journal in he direcion of he forward whirl moion, i.e. >0 if he equivalen damping (C eq ) i negaive (ee igure ), i.e. 1 ( C K r ) Ceq < 0 ω (56) 10-8 RTO-EN-AVT-14

29 RTO-EN-AVT whirl orbi Cro-coupled force K r e Damping force - C ω e Roor pin igure : orce Driving and Rearding Roor Whirl Moion. A he hrehold peed of inabiliy, K r C Ω. Thu, unable forward whirl moion occur for roor peed Ω ω. In he (,) coordinae yem, e co(ω) and e in(ω). Thu, he bearing dynamic force become Ω + Ω Ω + C e C e C d d ω ω ω ω ω ω ω ω 1 ) co( ) in( 1 ) co( ) in( (57) Noe ha (, ) oppoe he forward whirl moion for journal peed Ω < ω. or larger roor peed, he bearing force become poiive and aid o he growh of he forward whirl ampliude of moion, a hown graphically in igure 4.

30 whirl orbi -K -K Energy of orbi: Area_orbi (K -K ) K >0, K <0 <0, >0 igure 4: Repreenaion of ollower orce from Cro-Coupled Siffnee. The work performed by he bearing force during a full orbial period (Tπ/ω) equal [9] E ( ) d ( ω e) d ( ) + d ( ) E ( π e )( C ω K ) Area C ω r orbi d ( C ω K eq r ) ω e T (58) Noe ha E<0 i equivalen o negaive work, i.e. energy removed or diipaed from he roor-bearing yem. However when E>0, i.e. for Ω ω, he fluid film bearing add "energy" ino he roor-bearing yem hu driving he whirl moion forward. rom hi dicuion one can eaily deduce ha roor-bearing evidencing whirl orbi wih kewed area (harp ellipoid) will be le prone o roordynamic inabiliy, ee igure 5. Thi ype of dynamic repone i obained by deign and conrucion of a bearing generaing (direc) iffne aymmery, a given in muliple pad bearing configuraion (ellipic, muliple-lobe wih preload, preure-dam bearing). However, hee bearing are limied o fixed orienaion aic load and roor pin in only one direcion RTO-EN-AVT-14

31 Energy of orbi: Area x (K -K ) igure 5: Influence of Bearing Aymmery on Whirl Orbi. 5.4 Experimenal Meauremen of Roor-Bearing Syem Inabiliy The archival lieraure i abundan in experimenal and field decripion of evere inabiliie induced by fluid film bearing on roaing machinery. A an example of e conduced a he auhor laboraory on a high peed e rig, igure 6 depic recorded ampliude of moion veru haf peed in a rigid roor uppored on plain journal bearing. The diplacemen meauremen correpond o roor moion along he verical and horizonal plane (LV, LH). The curve wih larger ampliude denoe he oal ampliude of moion while he oher in ligh color how he filered ynchronou (1) moion wih low roll compenaion. The paage hrough a well-damped criical peed i eviden a ~ 8.5 krpm. A he haf peed increae, he ampliude of moion decreae. However, a a haf peed ~ wice he criical peed, he roor become violenly unable wih large ampliude moion nearly equalling he bearing clearance. igure 6: Ampliude of Roor Moion veru Shaf Speed. Experimenal Evidence of Roordynamic Inabiliy. RTO-EN-AVT

32 igure 7 depic he waerfall of he verical haf moion. The graph how he frequency conen of he vibraion ignal a he roor accelerae. The ynchronou moion are denoed by he 1 line. The whirl frequency raio i 0.50 a he one of he evere ubynchronou moion. A he peed increae, he whirl frequency lock a he yem naural frequency. Thi phenomenon i known a oil-whip. The roor wa everely damaged upon compleion of he experimen. igure 7: Waerfall of Recorded Roor Moion Demonraing Subynchronou Whirl. 6.0 CLOSURE Compreor, urbine, pump, elecric moor, elecric generaor and oher roaing machine are commonly uppored on fluid film bearing. In he pa, mo applicaion implemened common cylindrical plain journal bearing. A machine have achieved higher peed and larger power, roor dynamic inabiliy problem uch a oil whirl have brough he need o implemen oher bearing configuraion. Cuing axial groove in he bearing o upply oil flow ino he lubricaed urface generae ome of hee geomerie. Oher bearing ype have variou paern of variable clearance (preload and offe) o creae a pad film hickne ha ha rongly converging and diverging region, hu generaing a direc iffne for operaion even a he journal cenered poiion. Variou oher geomerie have evolved a well, uch a he iling pad bearing, which allow each pad o pivo, and hu o ake i own equilibrium poiion. Thi feaure uually reul in a rongly converging film region for each loaded pad and he near abence of cro-coupled iffne coefficien. Table and ummarize ome of he advanage and diadvanage of variou bearing in condened form. igure 8 how graphical keche for ome of he bearing configuraion below. Reference [1, 1, 14] offer imporan echnical informaion on he deign, operaion and abiliy conideraion for he mo common fluid film bearing ued in indurial applicaion, wih emphai in pump and compreor RTO-EN-AVT-14

33 Table : ixed Pad Non-Pre Loaded Journal Bearing Bearing Type Advanage Diadvanage Commen Plain Journal 1. Eay o make. Low Co 1. Mo prone o oil whirl Round bearing are nearly alway cruhed o make ellipical bearing Parial Arc 1. Eay o make. Low Co. Low horepower lo 1. Poor vibraion reiance. Oil upply no eaily conained Bearing ued only on raher old machine Axial Groove 1. Eay o make. Low Co 1. Subjec o oil whirl Round bearing are nearly alway cruhed o make ellipical or muli-lobe loaing Ring 1. Relaively eay o make. Low Co 1. Subjec o oil whirl (wo whirl frequencie from inner and ouer film (50% haf peed, 50% [haf + ring] peed) Ued primarily on high peed urbocharger for PV and CV engine Ellipical 1. Eay o make. Low Co. Good damping a criical peed 1. Subjec o oil whirl a high peed. Load direcion mu be known Probably mo widely ued bearing a low or moderae roor peed Offe Half (Wih Horizonal Spli) 1. Excellen uppreion of whirl a high peed. Low Co. Eay o make 1. air uppreion of whirl a moderae peed. Load direcion mu be known High horizonal iffne and low verical iffne - may become popular - ued ouide U.S. Three and our Lobe 1. Good uppreion of whirl. Overall good performance. Moderae co 1. Expenive o make properly. Subjec o whirl a high peed Currenly ued by ome manufacurer a a andard bearing deign RTO-EN-AVT

34 Table : Pad Journal Bearing wih Sep, Dam or Pocke, Tiling Pad Bearing Bearing Type Advanage Diadvanage Commen Preure Dam (Single Dam) 1. Good uppreion of whirl. Low co. Good damping a criical peed 4. Eay o make 1. Goe unable wih lile warning. Dam may be ubjec o wear or build up over ime. Load direcion mu be known Very popular in he perochemical indury. Eay o conver ellipical over o preure dam Muli-Dam Axial Groove or Muliple-Lobe 1. Dam are relaively eay o place in exiing bearing. Good uppreion of whirl. Relaively low co 4. Good overall performance 1. Complex bearing requiring deailed analyi. May no uppre whirl due o non bearing caue Ued a andard deign by ome manufacurer Hydroaic 1. Good uppreion of oil whirl. Wide range of deign parameer. Moderae co 1. Poor damping a criical peed. Require careful deign. Require high preure lubrican upply Generally high iffne properie ued for high preciion roor NON-IED PAD JOURNAL BEARINGS Bearing Type Advanage Diadvanage Commen Tiling Pad Journal Bearing lexure Pivo, Tiling Pad Bearing 1. Will no caue whirl (no cro coupling) 1. High Co. Require careful deign. Poor damping a criical peed 4. Hard o deermine acual clearance 5. Load direcion mu be known Widely ued bearing o abilize machine wih ubynchronou nonbearing relaed exciaion oil Bearing 1. Tolerance o mialignmen.. Oil-free 1. High co. Dynamic performance no well known for heavily loaded machinery. Prone o ubynchronou whirl Ued mainly for low load uppor on high peed machinery (APU uni) 10-4 RTO-EN-AVT-14

35 θ θ 0 W W (a) Ellipical bearing (b) -pad bearing θ θ W (c) 4-pad bearing W (d) 5-pad bearing Tiling pad bearing lexure pivo hydroaic pad bearing igure 8: Schemaic View of Variou Radial luid ilm Bearing Configuraion. REERENCES [1] Turbulence in luid ilm Bearing, L. San André, Lecure Noe (#8) in Modern Lubricaion, hp://phn.amu.edu/tribgroup, 00. [] Tribology ricion, Lubricaion & Wear, A. Szeri, Hemiphere Pub, RTO-EN-AVT

36 [] Caviaion in Liquid ilm Bearing, L. San André, Lecure Noe (#6) in Modern Lubricaion, hp://phn.amu.edu/tribgroup, 00. [4] Effec of luid Ineria on inie Lengh Sealed Squeeze ilm Damper, L. San André & J.M. Vance, ASLE Tranacion, 0,, pp. 84-9, [5] Turbomachinery Roordynamic, (chaper ), D. Child, John Wiley & Son, Inc., 199. [6] A Table of he Journal Bearing Inegral, J.. Booker, ASME Journal of Baic Engineering, pp. 5-55, [7] Roordynamic of Turbomachinery, J.M. Vance, J., Wiley Iner-Science Pub., [8] Self-Excied, Saionary Whirl Orbi of a Journal in a Sleeve Bearing, J. Lund, Ph.D. Thei, Renelaer Polyechnic Iniue, Troy, N.., [9] Dynamic of Simple Roor-luid ilm Bearing Syem, L. San André, Lecure Noe (#5) in Modern Lubricaion, hp://phn.amu.edu/tribgroup, 00. [10] Effec of Eccenriciy on he orce Repone of a Hybrid Bearing, L. San André, STLE Tribology Tranacion, 4, 4, pp , [11] The Sabiliy of an Elaic Roor in Journal Bearing wih lexible Suppor, J, Lund, ASME Journal of Applied Mechanic, pp , [1] Deign of Journal Bearing for Roaing Machinery, P. Allaire & R.D. lack, Proc. of he 10 h Turbomachinery Sympoium, TAMU, pp. 5-45, [1] luid ilm Bearing undamenal and ailure,. Zeidan & B. Herbage, Proc. of he 0 h Turbomachinery Sympoium, TAMU, pp [14] undamenal of luid ilm Journal Bearing Operaion and Modeling, M. He & J. Byrne, Proc. of he 4 h Turbomachinery Sympoium, TAMU, pp , RTO-EN-AVT-14

Buckling of a structure means failure due to excessive displacements (loss of structural stiffness), and/or

Buckling of a structure means failure due to excessive displacements (loss of structural stiffness), and/or Buckling Buckling of a rucure mean failure due o exceive diplacemen (lo of rucural iffne), and/or lo of abiliy of an equilibrium configuraion of he rucure The rule of humb i ha buckling i conidered a mode

More information

CHAPTER 7: SECOND-ORDER CIRCUITS

CHAPTER 7: SECOND-ORDER CIRCUITS EEE5: CI RCUI T THEORY CHAPTER 7: SECOND-ORDER CIRCUITS 7. Inroducion Thi chaper conider circui wih wo orage elemen. Known a econd-order circui becaue heir repone are decribed by differenial equaion ha

More information

TP B.2 Rolling resistance, spin resistance, and "ball turn"

TP B.2 Rolling resistance, spin resistance, and ball turn echnical proof TP B. olling reiance, pin reiance, and "ball urn" upporing: The Illuraed Principle of Pool and Billiard hp://billiard.coloae.edu by Daid G. Alciaore, PhD, PE ("Dr. Dae") echnical proof originally

More information

u(t) Figure 1. Open loop control system

u(t) Figure 1. Open loop control system Open loop conrol v cloed loop feedbac conrol The nex wo figure preen he rucure of open loop and feedbac conrol yem Figure how an open loop conrol yem whoe funcion i o caue he oupu y o follow he reference

More information

Lecture 5 Buckling Buckling of a structure means failure due to excessive displacements (loss of structural stiffness), and/or

Lecture 5 Buckling Buckling of a structure means failure due to excessive displacements (loss of structural stiffness), and/or AOE 204 Inroducion o Aeropace Engineering Lecure 5 Buckling Buckling of a rucure mean failure due o exceive diplacemen (lo of rucural iffne), and/or lo of abiliy of an equilibrium configuraion of he rucure

More information

To become more mathematically correct, Circuit equations are Algebraic Differential equations. from KVL, KCL from the constitutive relationship

To become more mathematically correct, Circuit equations are Algebraic Differential equations. from KVL, KCL from the constitutive relationship Laplace Tranform (Lin & DeCarlo: Ch 3) ENSC30 Elecric Circui II The Laplace ranform i an inegral ranformaion. I ranform: f ( ) F( ) ime variable complex variable From Euler > Lagrange > Laplace. Hence,

More information

Physics 240: Worksheet 16 Name

Physics 240: Worksheet 16 Name Phyic 4: Workhee 16 Nae Non-unifor circular oion Each of hee proble involve non-unifor circular oion wih a conan α. (1) Obain each of he equaion of oion for non-unifor circular oion under a conan acceleraion,

More information

Angular Motion, Speed and Velocity

Angular Motion, Speed and Velocity Add Imporan Angular Moion, Speed and Velociy Page: 163 Noe/Cue Here Angular Moion, Speed and Velociy NGSS Sandard: N/A MA Curriculum Framework (006): 1.1, 1. AP Phyic 1 Learning Objecive: 3.A.1.1, 3.A.1.3

More information

Chapter 7: Inverse-Response Systems

Chapter 7: Inverse-Response Systems Chaper 7: Invere-Repone Syem Normal Syem Invere-Repone Syem Baic Sar ou in he wrong direcion End up in he original eady-ae gain value Two or more yem wih differen magniude and cale in parallel Main yem

More information

2. VECTORS. R Vectors are denoted by bold-face characters such as R, V, etc. The magnitude of a vector, such as R, is denoted as R, R, V

2. VECTORS. R Vectors are denoted by bold-face characters such as R, V, etc. The magnitude of a vector, such as R, is denoted as R, R, V ME 352 VETS 2. VETS Vecor algebra form he mahemaical foundaion for kinemaic and dnamic. Geomer of moion i a he hear of boh he kinemaic and dnamic of mechanical em. Vecor anali i he imehonored ool for decribing

More information

5.2 GRAPHICAL VELOCITY ANALYSIS Polygon Method

5.2 GRAPHICAL VELOCITY ANALYSIS Polygon Method ME 352 GRHICL VELCITY NLYSIS 52 GRHICL VELCITY NLYSIS olygon Mehod Velociy analyi form he hear of kinemaic and dynamic of mechanical yem Velociy analyi i uually performed following a poiion analyi; ie,

More information

EECE 301 Signals & Systems Prof. Mark Fowler

EECE 301 Signals & Systems Prof. Mark Fowler EECE 31 Signal & Syem Prof. Mark Fowler Noe Se #27 C-T Syem: Laplace Tranform Power Tool for yem analyi Reading Aignmen: Secion 6.1 6.3 of Kamen and Heck 1/18 Coure Flow Diagram The arrow here how concepual

More information

Discussion Session 2 Constant Acceleration/Relative Motion Week 03

Discussion Session 2 Constant Acceleration/Relative Motion Week 03 PHYS 100 Dicuion Seion Conan Acceleraion/Relaive Moion Week 03 The Plan Today you will work wih your group explore he idea of reference frame (i.e. relaive moion) and moion wih conan acceleraion. You ll

More information

13.1 Accelerating Objects

13.1 Accelerating Objects 13.1 Acceleraing Objec A you learned in Chaper 12, when you are ravelling a a conan peed in a raigh line, you have uniform moion. However, mo objec do no ravel a conan peed in a raigh line o hey do no

More information

EE Control Systems LECTURE 2

EE Control Systems LECTURE 2 Copyrigh F.L. Lewi 999 All righ reerved EE 434 - Conrol Syem LECTURE REVIEW OF LAPLACE TRANSFORM LAPLACE TRANSFORM The Laplace ranform i very ueful in analyi and deign for yem ha are linear and ime-invarian

More information

6.8 Laplace Transform: General Formulas

6.8 Laplace Transform: General Formulas 48 HAP. 6 Laplace Tranform 6.8 Laplace Tranform: General Formula Formula Name, ommen Sec. F() l{ f ()} e f () d f () l {F()} Definiion of Tranform Invere Tranform 6. l{af () bg()} al{f ()} bl{g()} Lineariy

More information

EE202 Circuit Theory II

EE202 Circuit Theory II EE202 Circui Theory II 2017-2018, Spring Dr. Yılmaz KALKAN I. Inroducion & eview of Fir Order Circui (Chaper 7 of Nilon - 3 Hr. Inroducion, C and L Circui, Naural and Sep epone of Serie and Parallel L/C

More information

Chapter 6. Laplace Transforms

Chapter 6. Laplace Transforms 6- Chaper 6. Laplace Tranform 6.4 Shor Impule. Dirac Dela Funcion. Parial Fracion 6.5 Convoluion. Inegral Equaion 6.6 Differeniaion and Inegraion of Tranform 6.7 Syem of ODE 6.4 Shor Impule. Dirac Dela

More information

1 Motivation and Basic Definitions

1 Motivation and Basic Definitions CSCE : Deign and Analyi of Algorihm Noe on Max Flow Fall 20 (Baed on he preenaion in Chaper 26 of Inroducion o Algorihm, 3rd Ed. by Cormen, Leieron, Rive and Sein.) Moivaion and Baic Definiion Conider

More information

Introduction to SLE Lecture Notes

Introduction to SLE Lecture Notes Inroducion o SLE Lecure Noe May 13, 16 - The goal of hi ecion i o find a ufficien condiion of λ for he hull K o be generaed by a imple cure. I urn ou if λ 1 < 4 hen K i generaed by a imple curve. We will

More information

Problem Set If all directed edges in a network have distinct capacities, then there is a unique maximum flow.

Problem Set If all directed edges in a network have distinct capacities, then there is a unique maximum flow. CSE 202: Deign and Analyi of Algorihm Winer 2013 Problem Se 3 Inrucor: Kamalika Chaudhuri Due on: Tue. Feb 26, 2013 Inrucion For your proof, you may ue any lower bound, algorihm or daa rucure from he ex

More information

Chapter 6. Laplace Transforms

Chapter 6. Laplace Transforms Chaper 6. Laplace Tranform Kreyzig by YHLee;45; 6- An ODE i reduced o an algebraic problem by operaional calculu. The equaion i olved by algebraic manipulaion. The reul i ranformed back for he oluion of

More information

, the. L and the L. x x. max. i n. It is easy to show that these two norms satisfy the following relation: x x n x = (17.3) max

, the. L and the L. x x. max. i n. It is easy to show that these two norms satisfy the following relation: x x n x = (17.3) max ecure 8 7. Sabiliy Analyi For an n dimenional vecor R n, he and he vecor norm are defined a: = T = i n i (7.) I i eay o how ha hee wo norm aify he following relaion: n (7.) If a vecor i ime-dependen, hen

More information

CHAPTER HIGHER-ORDER SYSTEMS: SECOND-ORDER AND TRANSPORTATION LAG. 7.1 SECOND-ORDER SYSTEM Transfer Function

CHAPTER HIGHER-ORDER SYSTEMS: SECOND-ORDER AND TRANSPORTATION LAG. 7.1 SECOND-ORDER SYSTEM Transfer Function CHAPTER 7 HIGHER-ORDER SYSTEMS: SECOND-ORDER AND TRANSPORTATION LAG 7. SECOND-ORDER SYSTEM Tranfer Funcion Thi ecion inroduce a baic yem called a econd-order yem or a quadraic lag. Second-order yem are

More information

Let. x y. denote a bivariate time series with zero mean.

Let. x y. denote a bivariate time series with zero mean. Linear Filer Le x y : T denoe a bivariae ime erie wih zero mean. Suppoe ha he ime erie {y : T} i conruced a follow: y a x The ime erie {y : T} i aid o be conruced from {x : T} by mean of a Linear Filer.

More information

Curvature. Institute of Lifelong Learning, University of Delhi pg. 1

Curvature. Institute of Lifelong Learning, University of Delhi pg. 1 Dicipline Coure-I Semeer-I Paper: Calculu-I Leon: Leon Developer: Chaianya Kumar College/Deparmen: Deparmen of Mahemaic, Delhi College of r and Commerce, Univeriy of Delhi Iniue of Lifelong Learning, Univeriy

More information

Structural Dynamics and Earthquake Engineering

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

More information

Sample Final Exam (finals03) Covering Chapters 1-9 of Fundamentals of Signals & Systems

Sample Final Exam (finals03) Covering Chapters 1-9 of Fundamentals of Signals & Systems Sample Final Exam Covering Chaper 9 (final04) Sample Final Exam (final03) Covering Chaper 9 of Fundamenal of Signal & Syem Problem (0 mar) Conider he caual opamp circui iniially a re depiced below. I LI

More information

Introduction to Congestion Games

Introduction to Congestion Games Algorihmic Game Theory, Summer 2017 Inroducion o Congeion Game Lecure 1 (5 page) Inrucor: Thoma Keelheim In hi lecure, we ge o know congeion game, which will be our running example for many concep in game

More information

Non-uniform circular motion *

Non-uniform circular motion * OpenSax-CNX module: m14020 1 Non-uniform circular moion * Sunil Kumar Singh This work is produced by OpenSax-CNX and licensed under he Creaive Commons Aribuion License 2.0 Wha do we mean by non-uniform

More information

Laplace Transform. Inverse Laplace Transform. e st f(t)dt. (2)

Laplace Transform. Inverse Laplace Transform. e st f(t)dt. (2) Laplace Tranform Maoud Malek The Laplace ranform i an inegral ranform named in honor of mahemaician and aronomer Pierre-Simon Laplace, who ued he ranform in hi work on probabiliy heory. I i a powerful

More information

INTRODUCTION TO INERTIAL CONFINEMENT FUSION

INTRODUCTION TO INERTIAL CONFINEMENT FUSION INTODUCTION TO INETIAL CONFINEMENT FUSION. Bei Lecure 7 Soluion of he imple dynamic igniion model ecap from previou lecure: imple dynamic model ecap: 1D model dynamic model ρ P() T enhalpy flux ino ho

More information

18 Extensions of Maximum Flow

18 Extensions of Maximum Flow Who are you?" aid Lunkwill, riing angrily from hi ea. Wha do you wan?" I am Majikhie!" announced he older one. And I demand ha I am Vroomfondel!" houed he younger one. Majikhie urned on Vroomfondel. I

More information

Single Phase Line Frequency Uncontrolled Rectifiers

Single Phase Line Frequency Uncontrolled Rectifiers Single Phae Line Frequency Unconrolle Recifier Kevin Gaughan 24-Nov-03 Single Phae Unconrolle Recifier 1 Topic Baic operaion an Waveform (nucive Loa) Power Facor Calculaion Supply curren Harmonic an Th

More information

Motion In One Dimension. Graphing Constant Speed

Motion In One Dimension. Graphing Constant Speed Moion In One Dimenion PLATO AND ARISTOTLE GALILEO GALILEI LEANING TOWER OF PISA Graphing Conan Speed Diance v. Time for Toy Car (0-5 ec.) be-fi line (from TI calculaor) d = 207.7 12.6 Diance (cm) 1000

More information

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

!!#$%&#'()!#&'(*%)+,&',-)./0)1-*23) "#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

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

More information

NEUTRON DIFFUSION THEORY

NEUTRON DIFFUSION THEORY NEUTRON DIFFUSION THEORY M. Ragheb 4//7. INTRODUCTION The diffuion heory model of neuron ranpor play a crucial role in reacor heory ince i i imple enough o allow cienific inigh, and i i ufficienly realiic

More information

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

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

More information

IB Physics Kinematics Worksheet

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

More information

13.1 Circuit Elements in the s Domain Circuit Analysis in the s Domain The Transfer Function and Natural Response 13.

13.1 Circuit Elements in the s Domain Circuit Analysis in the s Domain The Transfer Function and Natural Response 13. Chaper 3 The Laplace Tranform in Circui Analyi 3. Circui Elemen in he Domain 3.-3 Circui Analyi in he Domain 3.4-5 The Tranfer Funcion and Naural Repone 3.6 The Tranfer Funcion and he Convoluion Inegral

More information

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

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

More information

Randomized Perfect Bipartite Matching

Randomized Perfect Bipartite Matching Inenive Algorihm Lecure 24 Randomized Perfec Biparie Maching Lecurer: Daniel A. Spielman April 9, 208 24. Inroducion We explain a randomized algorihm by Ahih Goel, Michael Kapralov and Sanjeev Khanna for

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

From Particles to Rigid Bodies

From Particles to Rigid Bodies Rigid Body Dynamics From Paricles o Rigid Bodies Paricles No roaions Linear velociy v only Rigid bodies Body roaions Linear velociy v Angular velociy ω Rigid Bodies Rigid bodies have boh a posiion and

More information

KINEMATICS IN ONE DIMENSION

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

More information

1. VELOCITY AND ACCELERATION

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

More information

NECESSARY AND SUFFICIENT CONDITIONS FOR LATENT SEPARABILITY

NECESSARY AND SUFFICIENT CONDITIONS FOR LATENT SEPARABILITY NECESSARY AND SUFFICIENT CONDITIONS FOR LATENT SEPARABILITY Ian Crawford THE INSTITUTE FOR FISCAL STUDIES DEPARTMENT OF ECONOMICS, UCL cemmap working paper CWP02/04 Neceary and Sufficien Condiion for Laen

More information

S. S. Patil, S.B.Javheri, C.G.Konapure. Main structure

S. S. Patil, S.B.Javheri, C.G.Konapure. Main structure ISSN: 77-3754 ISO 900:008 Ceriied Inernaional Journal o Engineering and Innovaive Technology (IJEIT) Volume, Iue 6, June0 Eecivene o Muliple Tuned Ma Damper S. S. Pail, S.B.Javheri, C.G.Konapure Abrac

More information

Let us start with a two dimensional case. We consider a vector ( x,

Let us start with a two dimensional case. We consider a vector ( x, Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our

More information

Linear Algebra Primer

Linear Algebra Primer Linear Algebra rimer And a video dicuion of linear algebra from EE263 i here (lecure 3 and 4): hp://ee.anford.edu/coure/ee263 lide from Sanford CS3 Ouline Vecor and marice Baic Mari Operaion Deerminan,

More information

Linear Motion, Speed & Velocity

Linear Motion, Speed & Velocity Add Iporan Linear Moion, Speed & Velociy Page: 136 Linear Moion, Speed & Velociy NGSS Sandard: N/A MA Curriculu Fraework (2006): 1.1, 1.2 AP Phyic 1 Learning Objecive: 3.A.1.1, 3.A.1.3 Knowledge/Underanding

More information

Flow networks. Flow Networks. A flow on a network. Flow networks. The maximum-flow problem. Introduction to Algorithms, Lecture 22 December 5, 2001

Flow networks. Flow Networks. A flow on a network. Flow networks. The maximum-flow problem. Introduction to Algorithms, Lecture 22 December 5, 2001 CS 545 Flow Nework lon Efra Slide courey of Charle Leieron wih mall change by Carola Wenk Flow nework Definiion. flow nework i a direced graph G = (V, E) wih wo diinguihed verice: a ource and a ink. Each

More information

Rectilinear Kinematics

Rectilinear Kinematics Recilinear Kinemaic Coninuou Moion Sir Iaac Newon Leonard Euler Oeriew Kinemaic Coninuou Moion Erraic Moion Michael Schumacher. 7-ime Formula 1 World Champion Kinemaic The objecie of kinemaic i o characerize

More information

Algorithmic Discrete Mathematics 6. Exercise Sheet

Algorithmic Discrete Mathematics 6. Exercise Sheet Algorihmic Dicree Mahemaic. Exercie Shee Deparmen of Mahemaic SS 0 PD Dr. Ulf Lorenz 7. and 8. Juni 0 Dipl.-Mah. David Meffer Verion of June, 0 Groupwork Exercie G (Heap-Sor) Ue Heap-Sor wih a min-heap

More information

Admin MAX FLOW APPLICATIONS. Flow graph/networks. Flow constraints 4/30/13. CS lunch today Grading. in-flow = out-flow for every vertex (except s, t)

Admin MAX FLOW APPLICATIONS. Flow graph/networks. Flow constraints 4/30/13. CS lunch today Grading. in-flow = out-flow for every vertex (except s, t) /0/ dmin lunch oday rading MX LOW PPLIION 0, pring avid Kauchak low graph/nework low nework direced, weighed graph (V, ) poiive edge weigh indicaing he capaciy (generally, aume ineger) conain a ingle ource

More information

18.03SC Unit 3 Practice Exam and Solutions

18.03SC Unit 3 Practice Exam and Solutions Sudy Guide on Sep, Dela, Convoluion, Laplace You can hink of he ep funcion u() a any nice mooh funcion which i for < a and for > a, where a i a poiive number which i much maller han any ime cale we care

More information

Mon Apr 2: Laplace transform and initial value problems like we studied in Chapter 5

Mon Apr 2: Laplace transform and initial value problems like we studied in Chapter 5 Mah 225-4 Week 2 April 2-6 coninue.-.3; alo cover par of.4-.5, EP 7.6 Mon Apr 2:.-.3 Laplace ranform and iniial value problem like we udied in Chaper 5 Announcemen: Warm-up Exercie: Recall, The Laplace

More information

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

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

More information

6 December 2013 H. T. Hoang - www4.hcmut.edu.vn/~hthoang/ 1

6 December 2013 H. T. Hoang - www4.hcmut.edu.vn/~hthoang/ 1 Lecure Noe Fundamenal of Conrol Syem Inrucor: Aoc. Prof. Dr. Huynh Thai Hoang Deparmen of Auomaic Conrol Faculy of Elecrical & Elecronic Engineering Ho Chi Minh Ciy Univeriy of Technology Email: hhoang@hcmu.edu.vn

More information

6.302 Feedback Systems Recitation : Phase-locked Loops Prof. Joel L. Dawson

6.302 Feedback Systems Recitation : Phase-locked Loops Prof. Joel L. Dawson 6.32 Feedback Syem Phae-locked loop are a foundaional building block for analog circui deign, paricularly for communicaion circui. They provide a good example yem for hi cla becaue hey are an excellen

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

Module 2: Analysis of Stress

Module 2: Analysis of Stress Module/Leon Module : Anali of Sre.. INTRODUCTION A bod under he acion of eernal force, undergoe diorion and he effec due o hi em of force i ranmied hroughou he bod developing inernal force in i. To eamine

More information

A STUDY ON COMPLICATED ROLL MOTION OF A SHIP EQUIPPED WITH AN ANTI-ROLLING TANK

A STUDY ON COMPLICATED ROLL MOTION OF A SHIP EQUIPPED WITH AN ANTI-ROLLING TANK 607 A STUDY ON COMPLICATED ROLL MOTION OF A SHIP EQUIPPED WITH AN ANTI-ROLLING TANK Harukuni Taguchi, Hirohi Sawada and Kauji Tanizawa Naional Mariime Reearch Iniue (JAPAN) Abrac Nonlinear roll moion of

More information

Wall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing.

Wall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing. MECHANICS APPLICATIONS OF SECOND-ORDER ODES 7 Mechanics applicaions of second-order ODEs Second-order linear ODEs wih consan coefficiens arise in many physical applicaions. One physical sysems whose behaviour

More information

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan Ground Rules PC11 Fundamenals of Physics I Lecures 3 and 4 Moion in One Dimension A/Prof Tay Seng Chuan 1 Swich off your handphone and pager Swich off your lapop compuer and keep i No alking while lecure

More information

Traveling Waves. Chapter Introduction

Traveling Waves. Chapter Introduction Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from

More information

ECE Spring Prof. David R. Jackson ECE Dept. Notes 39

ECE Spring Prof. David R. Jackson ECE Dept. Notes 39 ECE 6341 Spring 2016 Prof. David R. Jackon ECE Dep. Noe 39 1 Finie Source J ( y, ) For a phaed curren hee: p p j( k kyy) J y = J e + (, ) 0 The angenial elecric field ha i produced i: ( p ˆ ˆ) j( k + k

More information

2001 November 15 Exam III Physics 191

2001 November 15 Exam III Physics 191 1 November 15 Eam III Physics 191 Physical Consans: Earh s free-fall acceleraion = g = 9.8 m/s 2 Circle he leer of he single bes answer. quesion is worh 1 poin Each 3. Four differen objecs wih masses:

More information

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

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

More information

Additional Methods for Solving DSGE Models

Additional Methods for Solving DSGE Models Addiional Mehod for Solving DSGE Model Karel Meren, Cornell Univeriy Reference King, R. G., Ploer, C. I. & Rebelo, S. T. (1988), Producion, growh and buine cycle: I. he baic neoclaical model, Journal of

More information

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

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

More information

LAPLACE TRANSFORM AND TRANSFER FUNCTION

LAPLACE TRANSFORM AND TRANSFER FUNCTION CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions

More information

Math 221: Mathematical Notation

Math 221: Mathematical Notation Mah 221: Mahemaical Noaion Purpose: One goal in any course is o properly use he language o ha subjec. These noaions summarize some o he major conceps and more diicul opics o he uni. Typing hem helps you

More information

Today in Physics 218: radiation reaction

Today in Physics 218: radiation reaction Today in Physics 18: radiaion reacion Radiaion reacion The Abraham-Lorenz formula; radiaion reacion force The pah of he elecron in oday s firs example (radial decay grealy exaggeraed) 6 March 004 Physics

More information

Design of Controller for Robot Position Control

Design of Controller for Robot Position Control eign of Conroller for Robo oiion Conrol Two imporan goal of conrol: 1. Reference inpu racking: The oupu mu follow he reference inpu rajecory a quickly a poible. Se-poin racking: Tracking when he reference

More information

Notes on cointegration of real interest rates and real exchange rates. ρ (2)

Notes on cointegration of real interest rates and real exchange rates. ρ (2) Noe on coinegraion of real inere rae and real exchange rae Charle ngel, Univeriy of Wiconin Le me ar wih he obervaion ha while he lieraure (mo prominenly Meee and Rogoff (988) and dion and Paul (993))

More information

Chapter 12: Velocity, acceleration, and forces

Chapter 12: Velocity, acceleration, and forces To Feel a Force Chaper Spring, Chaper : A. Saes of moion For moion on or near he surface of he earh, i is naural o measure moion wih respec o objecs fixed o he earh. The 4 hr. roaion of he earh has a measurable

More information

Echocardiography Project and Finite Fourier Series

Echocardiography Project and Finite Fourier Series Echocardiography Projec and Finie Fourier Series 1 U M An echocardiagram is a plo of how a porion of he hear moves as he funcion of ime over he one or more hearbea cycles If he hearbea repeas iself every

More information

4. Electric field lines with respect to equipotential surfaces are

4. Electric field lines with respect to equipotential surfaces are Pre-es Quasi-saic elecromagneism. The field produced by primary charge Q and by an uncharged conducing plane disanced from Q by disance d is equal o he field produced wihou conducing plane by wo following

More information

y z P 3 P T P1 P 2. Werner Purgathofer. b a

y z P 3 P T P1 P 2. Werner Purgathofer. b a Einführung in Viual Compuing Einführung in Viual Compuing 86.822 in co T P 3 P co in T P P 2 co in Geomeric Tranformaion Geomeric Tranformaion W P h f Werner Purgahofer b a Tranformaion in he Rendering

More information

What is maximum Likelihood? History Features of ML method Tools used Advantages Disadvantages Evolutionary models

What is maximum Likelihood? History Features of ML method Tools used Advantages Disadvantages Evolutionary models Wha i maximum Likelihood? Hiory Feaure of ML mehod Tool ued Advanage Diadvanage Evoluionary model Maximum likelihood mehod creae all he poible ree conaining he e of organim conidered, and hen ue he aiic

More information

The Residual Graph. 11 Augmenting Path Algorithms. Augmenting Path Algorithm. Augmenting Path Algorithm

The Residual Graph. 11 Augmenting Path Algorithms. Augmenting Path Algorithm. Augmenting Path Algorithm Augmening Pah Algorihm Greedy-algorihm: ar wih f (e) = everywhere find an - pah wih f (e) < c(e) on every edge augmen flow along he pah repea a long a poible The Reidual Graph From he graph G = (V, E,

More information

s-domain Circuit Analysis

s-domain Circuit Analysis Domain ircui Analyi Operae direcly in he domain wih capacior, inducor and reior Key feaure lineariy i preerved c decribed by ODE and heir I Order equal number of plu number of Elemenbyelemen and ource

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

Applications of the Basic Equations Chapter 3. Paul A. Ullrich

Applications of the Basic Equations Chapter 3. Paul A. Ullrich Applicaions of he Basic Equaions Chaper 3 Paul A. Ullrich paullrich@ucdavis.edu Par 1: Naural Coordinaes Naural Coordinaes Quesion: Why do we need anoher coordinae sysem? Our goal is o simplify he equaions

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

Mathematische Annalen

Mathematische Annalen Mah. Ann. 39, 33 339 (997) Mahemaiche Annalen c Springer-Verlag 997 Inegraion by par in loop pace Elon P. Hu Deparmen of Mahemaic, Norhweern Univeriy, Evanon, IL 628, USA (e-mail: elon@@mah.nwu.edu) Received:

More information

Effects of Coordinate Curvature on Integration

Effects of Coordinate Curvature on Integration Effecs of Coordinae Curvaure on Inegraion Chrisopher A. Lafore clafore@gmail.com Absrac In his paper, he inegraion of a funcion over a curved manifold is examined in he case where he curvaure of he manifold

More information

1 CHAPTER 14 LAPLACE TRANSFORMS

1 CHAPTER 14 LAPLACE TRANSFORMS CHAPTER 4 LAPLACE TRANSFORMS 4 nroducion f x) i a funcion of x, where x lie in he range o, hen he funcion p), defined by p) px e x) dx, 4 i called he Laplace ranform of x) However, in hi chaper, where

More information

Numerical Dispersion

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

More information

ū(e )(1 γ 5 )γ α v( ν e ) v( ν e )γ β (1 + γ 5 )u(e ) tr (1 γ 5 )γ α ( p ν m ν )γ β (1 + γ 5 )( p e + m e ).

ū(e )(1 γ 5 )γ α v( ν e ) v( ν e )γ β (1 + γ 5 )u(e ) tr (1 γ 5 )γ α ( p ν m ν )γ β (1 + γ 5 )( p e + m e ). PHY 396 K. Soluion for problem e #. Problem (a: A poin of noaion: In he oluion o problem, he indice µ, e, ν ν µ, and ν ν e denoe he paricle. For he Lorenz indice, I hall ue α, β, γ, δ, σ, and ρ, bu never

More information

Geology 229 Engineering Geology. Lecture 12. Elementary Soil Mechanics (cont. 2) (West, Ch. 7)

Geology 229 Engineering Geology. Lecture 12. Elementary Soil Mechanics (cont. 2) (West, Ch. 7) Geology 229 Engineering Geology Lecure 12 Elemenary Soil echanic (con. 2) (We, Ch. 7) Ouline of hi Lecure 1. Compacion/Conolidaion Soil denificaion include compacion and conolidaion. Denificaion compacion

More information

Roller-Coaster Coordinate System

Roller-Coaster Coordinate System Winer 200 MECH 220: Mechanics 2 Roller-Coaser Coordinae Sysem Imagine you are riding on a roller-coaer in which he rack goes up and down, wiss and urns. Your velociy and acceleraion will change (quie abruply),

More information

Chapter 9 - The Laplace Transform

Chapter 9 - The Laplace Transform Chaper 9 - The Laplace Tranform Selece Soluion. Skech he pole-zero plo an region of convergence (if i exi) for hee ignal. ω [] () 8 (a) x e u = 8 ROC σ ( ) 3 (b) x e co π u ω [] ( ) () (c) x e u e u ROC

More information

ANALYSIS OF SOME SAFETY ASSESSMENT STANDARD ON GROUNDING SYSTEMS

ANALYSIS OF SOME SAFETY ASSESSMENT STANDARD ON GROUNDING SYSTEMS ANAYSIS OF SOME SAFETY ASSESSMENT STANDARD ON GROUNDING SYSTEMS Shang iqun, Zhang Yan, Cheng Gang School of Elecrical and Conrol Engineering, Xi an Univeriy of Science & Technology, 710054, Xi an, China,

More information

The Residual Graph. 12 Augmenting Path Algorithms. Augmenting Path Algorithm. Augmenting Path Algorithm

The Residual Graph. 12 Augmenting Path Algorithms. Augmenting Path Algorithm. Augmenting Path Algorithm Augmening Pah Algorihm Greedy-algorihm: ar wih f (e) = everywhere find an - pah wih f (e) < c(e) on every edge augmen flow along he pah repea a long a poible The Reidual Graph From he graph G = (V, E,

More information

Graphs III - Network Flow

Graphs III - Network Flow Graph III - Nework Flow Flow nework eup graph G=(V,E) edge capaciy w(u,v) 0 - if edge doe no exi, hen w(u,v)=0 pecial verice: ource verex ; ink verex - no edge ino and no edge ou of Aume every verex v

More information