VALIDATION OF FINITE ELEMENT PROGRAM FOR JOURNAL BEARINGS -- STATIC AND DYNAMIC PROPERTIES

Size: px
Start display at page:

Download "VALIDATION OF FINITE ELEMENT PROGRAM FOR JOURNAL BEARINGS -- STATIC AND DYNAMIC PROPERTIES"

Transcription

1 Univrsity of Kntucky UKnowldg Univrsity of Kntucky Mastr's Thss Graduat School 2004 VALIDATION OF FINITE ELEMENT PROGRAM FOR JOURNAL BEARINGS -- STATIC AND DYNAMIC PROPERTIES Raja Shkar Balupari Univrsity of Kntucky, Click hr to lt us know how accss to this documnt bnfits you. Rcommndd Citation Balupari, Raja Shkar, "VALIDATION OF FINITE ELEMENT PROGRAM FOR JOURNAL BEARINGS -- STATIC AND DYNAMIC PROPERTIES" Univrsity of Kntucky Mastr's Thss This Thsis is brought to you for fr and opn accss by th Graduat School at UKnowldg. It has bn accptd for inclusion in Univrsity of Kntucky Mastr's Thss by an authorizd administrator of UKnowldg. For mor information, plas contact

2 ABSTRACT OF THESIS VALIDATION OF FINITE ELEMENT PROGRAM FOR JOURNAL BEARINGS STATIC AND DYNAMIC PROPERTIES Th analysis of baring systms involvs th prdiction of thir static and dynamic charactristics. Th capability to comput th dynamic charactristics for hydrodynamic barings has bn addd to Baring Dsign Systm BRGDS, a finit lmnt program dvlopd by Dr. R.W. Stphnson, and th rsults obtaind wr validatd. In this softwar, a standard finit lmnt implmntation of th Rynolds quation is usd to modl th land rgion of th baring with prssur dgrs of frdom. Th assumptions of incomprssibl flow, constant viscosity, and no fluid inrtia trms ar mad. Th prssur solution is intgratd to giv th baring load, and th stiffnss and damping charactristics wr calculatd by a prturbation mthod. Th static and dynamic charactristics of 60, 120 and 180 partial barings wr vrifid and compard for a lngth to diamtr L/D ratio of 0.5. A comparison has also bn obtaind for th 120 baring with L/D ratios of 0.5, 0.75 and 1.0. A 360 -journal baring was vrifid for an L/D ratio of 0.5 and also compard to an L/D ratio of 1.0. Th rsults ar in good agrmnt with othr vrifid rsults. Th ffct of providing lubricant to th rcsss has bn shown for a 120 hybrid hydrostatic baring with a singl and doubl rcss. Kywords: Rynolds quation, Journal barings, Stiffnss, Damping, Sommrfld Numbr. Raja Shkar Balupari 07/28/2004 Copyright Raja Shkar Balupari 2004

3 VALIDATION OF FINITE ELEMENT PROGRAM FOR JOURNAL BEARINGS STATIC AND DYNAMIC PROPERTIES By Raja Shkar Balupari Dr. Kith Rouch Dirctor of Thsis Dr. Gorg Huang Dirctor of Graduat Studis 07/28/2004

4 RULES FOR THE USE OF THESIS Unpublishd thsis submittd for th Mastr s dgr and dpositd in th Univrsity of Kntucky Library ar as a rul opn for inspction, but ar to b usd only with du rgard to th rights of th authors. Bibliographical rfrncs may b notd, but quotations or summaris of parts may b publishd only with th prmission of th author, and with th usual scholarly acknowldgmnts. Extnsiv copying or publication of th thsis in whol or in part also rquirs th consnt of th Dan of th Graduat School of th Univrsity of Kntucky.

5 THESIS Raja Shkar Balupari Th Graduat School Univrsity of Kntucky 2004

6 VALIDATION OF FINITE ELEMENT PROGRAM FOR JOURNAL BEARINGS STATIC AND DYNAMIC PROPERTIES THESIS A thsis submittd in partial fulfillmnt of th rquirmnts for th dgr of Mastr of Scinc in Mchanical Enginring in th Collg of Enginring at th Univrsity of Kntucky By Raja Shkar Balupari Lxington, Kntucky Dirctor: Dr. Kith Rouch, Profssor of Mchanical Enginring Lxington, Kntucky 2004 Copyright Raja Shkar Balupari 2004

7 Ddication To my family, for thir vrlasting lov, support and ncouragmnt

8 ACKNOWLEDGEMENTS I would lik to xprss dp gratitud to Dr. Kith Rouch for having givn m this opportunity to work with him. Th prsnt work would not hav bn possibl without th consistnt guidanc and patinc shown by him. I would lik to thank Dr. R. W. Stphnson from Mchcon, Inc., for his guidanc and most importantly for providing m with important data. I would also lik to thank Dr. Raymond LBau and Dr. David Hrrin for srving on my thsis committ. I would lik to thank my family for thir trmndous support and ncouragmnt to pursu highr studis. I would particularly lik to thank my brothr Ravindra Balupari for his support and guidanc through out my graduat study. I would also lik to thank all my frinds for thir assistanc during my stay hr at th Univrsity of Kntucky. iii

9 TABLE OF CONTENTS Acknowldgmnts...iii List of Tabls. vii List of Figurs...viii List of Fils.xi Chaptr 1. Introduction Introduction Typs of Journal Barings Partial Arc Barings Plain Journal Barings 360º Axial Groov Barings Prssur Dam Barings Multi-lob Barings Tilting Pad Barings Hybrid Barings Litratur Rviw Scop of Thsis Solution for Fluid Film Barings Introduction Viscosity Solution of Fluid Film Barings Assumptions Rynolds Equation Finit Elmnt Mthod Finit Elmnt Solution of Rynolds Equation Baring Dsign Program Introduction Ovrviw Baring Coordinat Systm.26 iv

10 3.4 Solution Procdur Gnral Systm of Equations Eccntricity Film Thicknss Calculation Baring Forcs and Momnts Hydrostatic Solution Convrgnc Conditions Load Convrgnc Prssur Convrgnc Rcss Flow Convrgnc Damping and Stiffnss Calculation Rsults and Discussion Introduction Sommrfld Numbr Partial Arc Baring Partial Baring Partial Baring Partial Baring Effct of Baring Arc Effct of Lngth to Diamtr Ratio Full Journal Baring Effct of Lngth to Diamtr Ratio Partial Arc Hybrid Baring Conclusions and Futur Work Conclusions Futur Work. 70 Appndix 71 A1. Nomnclatur.. 71 A2.STIFFCALC Subroutin.. 72 A3. DAMPCALC Subroutin 74 A4. Sampl input fil.76 v

11 A5. Sampl output fil...78 Rfrncs..81 Vita...85 vi

12 LIST OF TABLES Tabl 4.1: Baring dsign data for 60º, 120º and 180º partial arc barings...45 Tabl 4.2: Stiffnss and damping cofficints for a 120 baring with no rcss.45 Tabl 4.3: Stiffnss and damping cofficints for a 120 baring with singl rcss...45 Tabl 4.4: Stiffnss and damping cofficints for a 120 baring with two rcsss...45 vii

13 LIST OF FIGURES Figur 1.1: Schmatic diagram of full plain journal barings Allair t al., Figur 1.2: Schmatic diagram of Partial arc baring Allair t al., Figur 2.1: Vlocity profil btwn two paralll plats to quantify viscosity Szri, Figur 2.2: Fluid film gomtry Allair t. al., Figur 2.3: Boundary conditions on unrolld baring surfac Allair t. al., Figur 2.4: Division of Fluid film into lmnts Nicholas, Figur 2.5: Th procss of finit lmnt analysis Bath, Figur 2.6: Th th two dimnsional simplx triangular lmnt Nicholas, Figur 3.1: Baring Coordinat systm Stphnson, Figur 3.2: Sid viw and prssur profil of journal baring Allair, Figur 3.3: Flow chart of BRGDS Baring analysis program Stphnson, Figur 3.4: Th physical rprsntation of Dynamic forc cofficints of a fluid film Baring San Andrs, Figur 4.1: Dimnsionlss dirct stiffnss cofficints for a 60º baring with L/D= Figur 4.2: Dimnsionlss dirct damping cofficints for a 60º baring with L/D= Figur 4.3: Dimnsionlss cross-coupld stiffnss for a 60º baring with L/D= Figur 4.4: Dimnsionlss cross-coupld damping for a 60º baring with L/D= Figur 4.5: Dimnsionlss stiffnss cofficints vs. Sommrfld numbr for a 120º baring with L/D= Figur 4.6: Dimnsionlss damping cofficints vs. Sommrfld numbr for a 120º baring with L/D= Figur 4.7: Dimnsionlss stiffnss cofficints vs. Sommrfld numbr for a 180º baring with L/D= Figur 4.8: Dimnsionlss damping cofficints vs. Sommrfld numbr for a 180º baring with L/D= Figur 4.9: Sommrfld numbr vs. ccntricity ratio for 60º, 120º and 180º barings 52 viii

14 Figur 4.10: Sommrfld numbr vs. ccntricity angl for 60º, 120º and 180º barings 52 Figur 4.11: Prssur at circumfrnc vs. baring angl for 60º, 120º and 180º barings...53 Figur 4.12: Load capacity vs. ccntricity ratio for 60º, 120º and 180º barings 53 Figur 4.13: Baring flow vs. Sommrfld numbr for 60º, 120º and 180º barings 54 Figur 4.14: Dimnsionlss dirct stiffnss proprtis for 60º, 120º and 180º barings...54 Figur 4.15: Dimnsionlss dirct damping proprtis for 60º, 120º and 180º barings...55 Figur 4.16: Cross-coupld stiffnss proprtis for 60º, 120º and 180º barings..55 Figur 4.17: Cross-coupld damping proprtis for 60º, 120º and 180º barings.56 Figur 4.18: Sommrfld numbr vs. ccntricity ratio for L/D=0.5, 0.75 and Figur 4.19: Sommrfld numbr vs. ccntricity angl for L/D=0.5, 0.75 and Figur 4.20: Load capacity vs. ccntricity ratio for L/D=0.5, 0.75 and Figur 4.21: Total flow vs. ccntricity ratio for L/D=0.5, 0.75 and Figur 4.22: Dimnsionlss dirct stiffnss charactristics for L/D=0.5, 0.75 and Figur 4.23: Dimnsionlss dirct damping charactristics for L/D=0.5, 0.75 and Figur 4.24: Dimnsionlss cross-coupld stiffnss for L/D=0.5, 0.75 and Figur 4.25: Dimnsionlss cross-coupld damping for L/D=0.5, 0.75 and Figur 4.26: Dimnsionlss dirct stiffnss for a 360º baring with L/D= Figur 4.27: Dimnsionlss dirct damping for a 360º baring with L/D= Figur 4.28: Dimnsionlss cross-coupld stiffnss for a 360º baring with L/D= Figur 4.29: Dimnsionlss cross-coupld damping for a 360º baring with L/D= Figur 4.30: Sommrfld numbr vs. ccntricity ratio for 360º baring L/D=0.5 and Figur 4.31: Eccntricity angl vs. Sommrfld numbr for a 360º baring L/D=0.5 and Figur 4.32: Nodal prssurs vs. angl of baring for a 360º baring L/D=0.5 and Figur 4.33: Load capacity vs. ccntricity ratio for a 360º baring L/D=0.5 and Figur 4.34: Dimnsionlss flow vs. Sommrfld numbr for a 360º baring ix

15 L/D=0.5 and Figur 4.35: Dimnsionlss dirct stiffnss charactristics for a 360º baring L/D=0.5 and Figur 4.36: Dimnsionlss dirct damping charactristics for a 360º baring L/D=0.5 and Figur 4.37: Dimnsionlss cross-coupld stiffnss for a 360º baring L/D=0.5 and Figur 4.38: Prssur profil at various rcss flows for a 120 baring with singl rcss Figur 4.39: Prssur profil at various journal spds for a 120 baring with singl rcss Figur 4.40: Prssur profil at various rcss flows for a 120 baring with two rcsss Figur 4.41: Prssur profil at various journal spd for a 120 baring with two rcsss...68 x

16 LIST OF FILES Filnam Fil siz 1. RajaThs.pdf MB xi

17 CHAPTER ONE INTRODUCTION 1.1 Introduction A fluid film baring may b dfind as a baring in which th opposing or mating surfacs ar compltly sparatd by a layr of fluid lubricant Wilcock, A widly usd baring typ is th plain journal baring with application in comprssors, turbins, pumps, lctric motors and lctric gnrators. Journal Barings consists of two cylindrs rotating rlativ to ach othr. Th outr cylindr baring is stationary and th innr ring shaft rotating with an angular vlocity is calld th Journal as shown in figur 1.1. Th main purpos of th journal baring is to support th rotating machinry by providing sufficint lubrication to sparat th moving parts and to minimiz th friction du to rotation Allair, Th high-prssur fluid film in th claranc btwn th journal and th baring du to rotation of th journal provids th hydrodynamic film lubrication, and th load capacity to th baring. Th displacmnt of th shaft cntr rlativ to th baring cntr is known as ccntricity. Narly all havy industrial turbomachins us fluid film barings of som typ to support th shaft wight and control motions causd by unbalanc forcs. Th two primary advantags of fluid film barings ovr rolling lmnt barings ar thir suprior ability to absorb nrgy to dampn vibrations, and thir longvity du to th absnc of rolling contact strsss. Th damping is vry important in many typs of rotating machins whr th fluid film barings ar oftn th primary sourc of th nrgy absorption ndd to control vibrations. Fluid film journal barings also play a major rol in dtrmining rotordynamic stability, making thir carful slction and application is a crucial stp in th dvlopmnt of suprior rotor-baring systms. Journal barings ar incrdibly long-livd providd th lubricant is contaminant fr and sufficintly supplid. Thus dynamic analysis of hydrodynamic barings is important du to th forcs imposd on th shaft from machin unbalanc forcs, arodynamic forcs, and xtrnal xcitations from sals and couplings Allair,

18 1.2 Typs of Journal Barings Journal barings may b groupd into two main classs - full baring in which th baring arc compltly surrounds th journal and partial baring in which th arc is lss than 360º Raimondi and Boyd, Barings ar furthr classifid as plain journal barings, axial groov barings, prssur dam barings, multi-lob journal barings, tilting pad journal barings, and hybrid hydrostatic barings Partial Arc Barings Partial barings as shown in figur 1.2 may b furthr classifid as cntrally or ccntrically loadd dpnding upon whthr th load lin biscts th baring pad arc or divids it ccntrically. Typically th radius of th baring is gratr than th radius of th journal by an amount qual to th radial claranc. Whn th load applid is primarily unidirctional thr is no nd of full journal barings and instad a singl partial arc baring may b usd Allair, Partial arc barings ar usd in rlativly low spd applications. Th partial arc pad may also b dsignd with a hydrostatic rcss thus nabling fluid film lubrication du to both hydrostatic and hydrodynamic ffcts. Partial arc barings ar important bcaus thy form th building blocks for axial groov, multi-lob and tilting pad barings. Th prsnt work dals with th computation of dynamic charactristics of partial arc barings Plain Journal Barings 360º Plain barings ar composd of a cylindrical shaft of radius R, calld a journal, rotating with an angular vlocity, ω, about its axis in a cylindrical bushing of radius R+c and lngth L. Th cntr of th baring is labld as O b and th cntr of journal as O j as shown in figur 1.1. Undr stady stat conditions th journal cntr rmains at a constant ccntricity and attitud angl Φ for a givn load W acting on th shaft Allair t al., 1975, Allair t al., Th plain journal baring is th simplst and most common radial baring dsign, whr a plain cylindrical shll ncircls th shaft. Plain slv barings hav th highst cross coupling of all barings and ar suitabl for highly loadd or low-spd shafts. Advantags ar low cost and as of manufactur. Exampls includ automotiv crankshaft barings and low-spd or highly loadd turbo machinry applications. 2

19 1.2.3 Axial Groov Barings Th axial groov baring is similar to th plain journal baring, but with two or mor axial groovs addd for oil supply. As in th plain journal baring, thr is no prload and has a high tndncy for instability Allair t. al., Advantags ar low cost and as of manufactur. Axial groov barings ar vry common in many typs of commrcial machinry including turbins, gnrators, motors, pumps, and comprssors, and hav slightly bttr dynamic proprtis than do plain journal barings Prssur Dam Barings Th prssur dam baring is a fixd-gomtry baring improvs rotor dynamic stability than plain journal barings. A pockt is milld in th uppr unloadd half of th baring that nds in a "dam". A sharp prssur pak is cratd at th dam du to fluid inrtia ffcts. Th prssur pak imposs a downward load on th journal, forcing th shaft down to a gratr ccntricity that inhrntly improvs stability bcaus of th stiffnss and damping asymmtry inducd. Prssur dam barings hav rlativly high powr losss bcaus of th built-in load. Thy ar not suitabl for applications whr th load dirction changs bcaus th top half is pocktd for th dam. Th dam also rstricts thm to unidirctional opration and th oil must b kpt clan to prvnt sludg accumulation in th pockt. Manufacturing of prssur dam barings is mor difficult and xpnsiv than plain barings, as th dam gomtry is milld sparatly from th bor and must b prcis. Prssur dam barings ar usd primarily in highr spd applications to incras stability ovr othr typs of fixd-gomtry barings Multi-lob Barings Multi-lob barings ar ssntially barings with mor than on baring pad that nabl a combination of numbr of pads, rotation of baring, claranc, prload, and offst. This producs a stabilizing ffct on th shaft and can incras load capacity. Th cntr of curvatur of th pad radius R+c may b offst from th bushing/baring cntr prload. Th prload factor in this cas is givn by th ratio of distanc of cntr of curvatur from bushing cntr to th claranc in th barings Allair t al., 1975, Allair t al., Howvr, it can also consum mor powr du to th prload. Multi-lob barings can b ithr bidirctional or unidirctional, dpnding on whthr th lobs hav symmtric or asymmtric lobs. Multi-lob barings ar 3

20 difficult and xpnsiv to manufactur bcaus of th prcis machining oprations rquird. Thy ar commonly usd in smallr, high-spd machins rquiring high load capacity or high stability Tilting Pad barings Th tilting pad baring diffrs from th multi-lobd baring in that ach pad rotats about a pivot, nabling ach pad to hav highr dgrs of frdom corrsponding to movmnt about th pivot point. Th pad tilts such that its cntr of curvatur movs to crat a convrging pad film thicknss. Ths pads ar abl to follow th shaft motion rsulting in littl cross-coupld stiffnss and damping Allair t. al., Thy ar widly usd to stabiliz machins that hav sub synchronous vibrations. Ths barings hav highr powr losss and highr cost to manufactur Hybrid Barings In hydrostatic barings, high-prssur lubricant is fd to a rcss in th pad. Thy hav improvd high-spd load carrying capacity compard to th hydrodynamic opration. Thir dsign dpnds on th hydrodynamic ffct in addition to th hydrostatic ffct to achiv ncssary load support. Hybrid barings ar advantagous ovr pur hydrodynamic barings in that war can b avoidd at starting or stopping, thy tolrat substantial loads abov normal dsign load, withstand havy dynamic loadings which vary widly in dirction of rotation, and may allow for dsign with a smallr journal diamtr, rducing initial cost and oprating powr consumption. All th barings apart from th tilting pad barings discussd abov ar composd of fixd pad plain baring gomtry. Thus study of th plain baring is significant, as th majority of th barings possss similar charactristics as th plain baring in gomtry as wll as opration. Full barings and cntrally loadd partial arc barings of th claranc typ ar common, and it is to ths barings that th analysis prsntd in this dissrtation is applicabl. Also th dynamic analyss hav bn prformd on th fixd partial arc barings. 4

21 1.3 Litratur Rviw As alrady statd th dynamic analysis of hydrodynamic barings is ssntial as thy provid stability and control mchanical vibrations that occur in rotating machinry. Svral numrical tchniqus hav bn proposd to provid th solution of th fluid film lubrication problm. Raimondi and Boyd 1958 providd th numrical solution of finit journal barings for application in analysis and dsign using finit diffrnc mthod. Li 1977 and Hays 1959 providd th analytical dynamics of partial journal barings and finit journal barings rspctivly using variational mthods to solv th Rynolds quation. Th finit lmnt mthod has bn usd prominntly for som yars to continuum and fild problms Zinkiwicz, Rddi 1969 prsntd th finit lmnt solution for incomprssibl lubrication problms of complx gomtris without th loss of accuracy as th finit diffrnc mthod. Wada, t al has succssfully usd th finit lmnt mthod for finit width and infinit width lubrication problms. Thy also showd that finit lmnts might produc mor accurat solutions for journal barings than finit diffrnc mthods. Hubnr and Bookr 1972 applid th finit lmnt mthod to th gnral lubrication problm with a systmatic dscription of procdurs. Hubnr 1975 providd th most rcnt and xtnsiv tratmnt of fluid film lubrication. Allair, t al dvlopd a systmatic matrix approach for finit lmnts which automatically produc minimum bandwidth of algbraic quations. By organizing th labling in matrix form throughout th analysis, th solution could b asily obtaind using Gaussian limination mthods or othr mthods. Latr, Allair, t al providd a prssur paramtr mthod for th finit lmnt solution of Rynolds quation to improv th accuracy significantly. Rao 1982 providd th mathmatical formulation of th finit lmnt mthod to th hydrodynamic lubrication problm govrnd by th Rynolds quation. Allair, t al dvlopd a finit lmnt program FINBRG1 to comput th baring charactristics of plain, partial arc, axial groov, multi-lob and tilting pad barings. Allair, t al also prformd th dynamic analysis of incomprssibl fluid film barings. Stphnson 1997 dvlopd Baring Dsign Systm, a program to prform baring analysis using th finit lmnt mthod. 5

22 1.4 Scop of Thsis 'BRGDS is an acronym for Baring Dsign Systm. Th initial program dvlopd by Dr. Stphnson has th capability to comput th baring prssur distribution, film thicknss, baring forcs and momnts for a singl partial arc baring with multipl rcsss in th pad. Th program includs turbulnc as an option. Eithr ccntricity or load can b input to th barings. If load is input, th ccntricity and attitud angl ar solvd itrativly. Th motivation for th prsnt work has bn to includ th capability to comput th damping charactristics for plain, partial arc and hybrid hydrodynamic barings and to furthr validat th rsults. Th finit lmnt solution for th partial arc and full plain journal baring including th baring damping calculations ar th highlight of this thsis. Th baring program is basd on Fortran 77. Th function STIFFCALC to comput stiffnss and DAMPCALC that has bn introducd to comput th damping charactristics of th barings ar shown in Appndix A2 and A3. Th rsults obtaind wr compard and vrifid with othr rsults. Most of th rsults usd for comparison wr obtaind from Dynamic Rotor Baring Systm program, DyRoBs. DyRoBS is a powrful and sophisticatd softwar tool authord by Wn Jng Chn for rotor dynamics including comprhnsiv baring analysis, with rsults courtsy of Stphnson R. W. from Mchcon, Inc. Th nxt chaptr of this work dals with th solution of fluid film barings, viscosity, assumptions that ar usd in th analysis and thir validity. Th Rynolds quation for fluid barings has bn shown and rvisd for finit width barings. Also th finit lmnt solution of th Rynolds quation that is th backbon of th prsnt work has bn dtaild. Chaptr thr dals with th baring dsign program BRGDS. Th gnral systms of quation, baring coordinat systm, solution procdur and convrgnc conditions usd in th program ar documntd. Howvr, th main focus is on th computation of stiffnss and damping charactristics for fluid film barings. Th rsults that ar obtaind by th analysis ar shown in chaptr four. Comparisons hav bn mad for th stiffnss and damping charactristics to succd th prsnt program for furthr analysis. Plots showing th analysis of various othr baring typs ar also shown and commnts on th rsults hav bn providd. Finally, th scop for futur work has bn outlind basd on th prsnt work. 6

23 Y x Baring ω o b W o j hө Ө y X Journal Ф Lin of Cntrs Load Lin Ψ Figur 1.1: Schmatic diagram of a full plain journal baring Allair t al.,

24 Y Lin of cntrs W ω hθ R+c O b R O j X Ф Figur 1.2: Schmatic diagram of partial arc baring Allair t. al., Copyright Raja Shkar Balupari

25 CHAPTER TWO SOLUTION FOR FLUID FILM BEARINGS 2.1 Introduction Fluid film barings hav bn dfind as barings in which th opposing or mating surfacs ar compltly sparatd from ach othr by a layr of fluid lubricant. Thus it is ssntial to undrstand th fundamntal procss by which th fluid is maintaind whil supporting th load. This chaptr dals with th basic principls of fluid film barings, viscosity of fluids, th Rynolds quation and th finit lmnt intrprtation of th Rynolds quation. Th basic principl can b undrstood by considring high-spd fluid flow btwn two plan surfacs with a convrging wdg and lubricant btwn th surfacs. Th convrgnc coupld with high-spd fluid flow and fluid viscosity gnrats a high-prssur fluid film that supports th load. In th cas of a plain journal baring, th convrging wdg is formd at th bottom of th baring du to th wight of th shaft and any othr xtrnal applid load. Th fluid is pulld into th rgion undr th shaft du to th shar forcs gnratd by th shaft rotation. Th fluid is forcd into th convrging film thicknss at th bottom of th shaft producing a high-prssur film Allair, This film supports th wight of th rotating machin componnts and prvnts th shaft from touching th bushing surfac. For a givn ccntricity th fluid film has convrging and divrging gomtry, such that cavitations may occur in th divrging portion. It is thus vry important to b abl to prdict th prssur distribution and load capacity of th baring. 2.2 Viscosity Th viscosity of th fluid lubricant is its most significant physical charactristic as far as fluid film barings ar concrnd. In this analysis of hydrodynamic barings th fluid is incomprssibl, and th viscosity rmains constant throughout th flow. Th amount of prssur gnratd dpnds on th viscosity and dnsity of th fluid. Viscosity is th proprty that dfins th rsistanc of th fluid to motion. Th viscosity is du to th molcular attraction btwn th adjacnt layrs of fluid. Any applid shar forc will caus th fluid to mov, rsulting in a 9

26 rsistanc to flow that dpnds on th viscosity of th fluid. Whn a fluid is at rst, thr is no rsistanc to an applid shar forc. In a positiv sns, viscosity controls th fluid flow out of th baring but it is th caus of powr consumption of th baring du to lubricating shar Wilcock, A mathmatical quation to quantify viscosity can b obtaind from figur 2.1. A fluid film of thicknss H sparats a lowr plat that is stationary from an uppr plat. Th uppr plat movs with vlocity U du to application of forc F. Th fluid particl O in contact with th lowr plat has zro vlocity, and th fluid particl Q attachd to th uppr plat movs with th vlocity U. For a Nwtonian fluid th fluid vlocity incrass linarly from th stationary plat to th moving plat as shown in th figur 2.1. Th viscosity µ of any intrmdiat particl P locatd at a distanc y from th lowr plat moving with vlocity v can b obtaind from quation 2.1. Rarranging th trms w obtain quation 2.2, whr τ is th shar strss dvlopd in th fluid, A is th ara of moving plat in contact with th fluid, and U/H is th vlocity gradint btwn th plats. Th flow causd by viscosity is known as Coutt flow. UA F = µ 2.1 H τ τ µ = = 2.2 v y U H Thus dynamic or absolut viscosity has bn dfind as th forc rquird for moving a flat unit plat locatd at a unit distanc from a stationary plat by a unit vlocity whr th spac btwn th plats is filld by th fluid. Th units of absolut viscosity ar lb.s/ft 2. Fluids whos viscosity rmains constant with chang in vlocity gradint ar calld Nwtonian fluids, whil fluids that do not hav a linar rlationship ar calld non-nwtonian fluids. Th ratio of absolut viscosity to th dnsity of th fluid is calld kinmatic viscosity Szri, 1999 as in quation 2.3. Th units of kinmatic viscosity ar ft 2 /sc. µ ν = 2.3 ρ 10

27 Solution of fluid film barings Th war fr transfr of forc in hydrodynamic lubricatd barings is basd on th application of prssur in th lubricating film that balancs th xtrnal forcs on th baring. Osborn Rynolds and othrs Pkn, t al, 1983 drivd from th common Navir-Stoks quations 2.4, 2.5 and 2.6 and th continuity quation 2.9, a diffrntial quation for th calculation of th prssur distribution Szri, Th Laplacian and th divrgnc oprators ar shown in quations 2.7 and 2.8. u x u x p f z u w y u v x u u t u x + + = µ µ ρ ρ X-dirction 2.4 u y v y p f z v w y v v x v u t v y + + = µ µ ρ ρ Y-dirction 2.5 u z w z p f z w w y w v x w u t w z + + = µ µ ρ ρ Z-dirction z y x + + = 2.7 z u y u x u u + + = = w z v y u x t ρ ρ ρ ρ 2.9 Th Navir-Stoks quations in which th fluid inrtia, body, prssur and viscous forcs ar includd ar too complicatd to yild an analytical solution for most problms. Th Rynolds quation is a simplifid vrsion of Navir-Stoks quation. Th analytical solution of th journal baring is obtaind by using th most popular form of th Rynolds quation. Th following assumptions hav to b considrd in ordr to obtain th solution.

28 2.3.1 Assumptions Th xact rprsntation of th fluid flow in mathmatical trms using Navir-Stoks quation is complicatd for fluid film barings. Crtain assumptions may b mad to obtain a standard quation known as Rynolds quation Th following assumptions ar typical in th analysis of fluid film barings Camron, Body forcs ar nglctd. This mans that thr ar no inrtial or othr distributd forcs acting on th fluid. Th ffct of gravity of th fluid is nglctd 2. Th prssur of th fluid is assumd constant across th thicknss of th film. 3. Th curvatur of th surfacs is larg compard to th film thicknss. Surfac vlocitis ar considrd unidirctional. 4. Thr is no slip at th boundaris. Th vlocity of th oil layr adjacnt to th boundary is th sam as that of th boundary. Th fluid attachd to th baring surfac is stationary whil th fluid nar th rotor or shaft has th sam angular vlocity as th shaft itslf. Apart from ths assumptions thr ar crtain assumptions that ar considrd to nabl simplification of th mathmatical quations. Ths assumptions ar not ssntial but mak th quations simplr Allair, Th fluid lubricant usd is Nwtonian. This mans th fluid obys Nwton s law of viscosity and th fluid strss is proportional to th rat of shar. 6. Fluid film flow in th barings is laminar. 7. Th viscosity of th fluid dos not chang. 8. Th dnsity of th fluid rmains constant Rynolds Equation Rynolds quation with appropriat boundary conditions dscribs th prssur btwn two surfacs sparatd by a thin fluid film. Considr two surfacs moving rlativ to ach othr sparatd by a constant fluid film of thicknss as shown in figur 2.2. Th most gnralizd rprsntation of Rynolds quation for fluid film lubrication, dvlopd by Dowson 1962, can b rprsntd as shown in quation

29 3 ρh p + x k x x µ 3 ρh p = y k x y µ hu 2 ρ + hv x 2 ρ ρ + h y t x U U ρ h y V V ρ h W + ρ W 1 Whr, U 2 =U 2 i+v 2 j+w 2 k, is th vlocity of th uppr plat U 1 =U 1 i+v 1 j+w 1 k, is th vlocity of th lowr plat i, j and k ar unit vlocity vctors in th x, y, z dirctions rspctivly. U, V and W ar vlocity componnts in th x, y, z dirctions rspctivly. ρ x,y = fluid film dnsity h x,y = fluid film thicknss µ x,y = fluid film viscosity p x,y = fluid film prssur v x,y = diffusion vlocity t = tim k x = turbulnt constant k x = 12, for laminar flow k x = *R 0.9, for turbulnt flow Th lft hand sid of quation 2.10 can b rwrittn as shown in quation LHS Hr, 3 ρh = p 12 µ = i + j and for incomprssibl fluid flow x y ρ = x ρ = y ρ =0 t 2.11 Thus th right hand sid can b rwrittn as shown in quation 2.12; for slidr barings w considr quation

30 RHS ρ h x h y t U U + ρ V V + ρh = h = W 2 W 1 t Nglcting strtch ffcts th right hand sid can b rprsntd as quation RHS x y x x ρhu + ρhv ρhu ρhv + ρh = Thus from quations 2.12 and 2.14, th full Rynolds quation can b writtn as quation t ρh p = [ ρh U U ] + [ ρh V V ] + ρh k x µ x y t With U= 1 U U, th avrag vlocity of uppr and lowr surfacs and nglcting th body 2 [ ] 1 2 forcs F with diffusion ffcts, th Rynolds quation can b rwrittn as in quation 2.16 Rddi, ρh p = k xµ t ρhu + ρh + ρv 2.16 Thus, quation 2.16 dscribs th Rynolds quation for incomprssibl and isoviscous fluid films. Body forcs ar includd in th quation. Onc th Rynolds quation is solvd and th prssur is known, th avrag fluid vlocity, u is givn by quation 2.17 and th mass flow rat q is givn by th quation h u = U P k µ x 2.17 q = ρhu 2.18 Spcifying ithr th prssur or th avrag normal fluid vlocity at ach point along th thin fluid film boundary complts th formulation of th fluid film boundary valu problm. Lt C P dnot th portion of boundary C whr th prssur is considrd known and C q rprsnt th rmaining portion of boundary whr th normal fluid vlocity or normal mass flow rat is known as shown in figur 2.3. This implis P is known on C P and q.n is known on C q, whr n is 14

31 th outward normal to C q. It can b shown that th prssur must b spcifid at a minimum of on point along th boundary for th solution of P to th boundary valu problm to b uniqu. It has bn shown that th prssur P that minimizs this functional whil satisfying th boundary conditions on C P is also a solution of th Rynolds quation. Th variational principl for incomprssibl, isoviscous films of any thin film gomtry is to find a prssur P x,y which minimizs th function shown in quation 2.19 and satisfis th boundary condition on C P Allair t. al., Th prssur thus obtaind is th sam as that which is th solution to th Rynolds quation With ths quations, th thin fluid film problm is compltly formulatd xcpt for inrtia and nrgy ffcts. J 3 ρh 2kxµ t P P ρhu P + ρh + ρv P da+ q n = A C q PdC 2.19 For a fluid film journal baring th lowr surfac is fixd, fluid film viscosity and dnsity ar constant, and body forcs and diffusion forcs ar nglctd. Thus, th Rynolds quation in cylindrical coordinat systm for laminar flow of a circular baring can b rprsntd as in quation 2.20 and th corrsponding functional as quation r 2 J 3 h p [ ] + θ µ Θ 3 h p [ ] = z µ z h 6 ϖ + j θ h t ρh 1 P P ρ P ρh P = + ϖ b + ϖ j h + A 2k xµ R θ y 2 θ t P Rdθdy Finit Elmnt Mthod Finit lmnt procdurs ar vry widly usd in nginring analysis. Th procdurs ar mployd xtnsivly in th analysis of solids and structurs and of hat transfr and fluids. Finit lmnt mthods ar usful in virtually vry fild of nginring analysis. Th finit lmnt mthod is usd to solv physical problms in nginring analysis and dsign. Th physical problm typically involvs th structur or structural componnt subjctd to 15

32 crtain loads. Th idalization of th physical problm to th mathmatical modl rquirs crtain assumptions that togthr lad to diffrntial quations govrning th mathmatical modl. Th finit lmnt analysis provids an approximat solution to th mathmatical modl. Sinc th finit lmnt procdur is a numrical procdur it is ncssary to assss th solution accuracy. If th accuracy critria ar not mt, th numrical solution has to b rpatd with rfind solution paramtrs such as a finr msh until a sufficint accuracy is rachd. Also th choic of an approximat mathmatical modl is crucial and compltly dtrmins th insight into th actual physical problm that w can obtain by th analysis. A flow chart rprsntation of th gnral procdur for finit lmnt analysis is shown in figur 2.5 Bath, In th fluid film baring application th finit lmnt mthod rducs th fild problm in which th prssur is a continuous function of spac to a problm in which th prssur is valuatd only at crtain spcific nodal points. Btwn th nods th prssur is assumd to vary according to a simpl function calld an intrpolation function. Th problm thn bcoms on of valuating a systm of algbraic quations for unknown nodal prssurs. Th cofficints of algbraic quations dpnd upon th nodal coordinats, th intrpolation polynomials chosn and othr factors. Th finit lmnt mthod is convnintly drivd from a variational principl. In th cas of fluid film barings th variational principl has bn drivd from th Rynolds quation. Th variational principl consists of a functional in th form of an intgral ovr th baring ara. Th baring ara to b analyzd is dividd up into a finit numbr of small sub rgions, calld lmnts, of simpl gomtric shap. Th lmnts must b chosn such that thr ar no gaps or ovrlapping occurs btwn th lmnts. On th boundary of ach lmnt, a numbr of nods ar chosn whr prssur is to b valuatd. Intrior nods can b chosn but thy ar not usually mployd. An intrpolation function is chosn to approximat th prssur within ach lmnt. Th function for a givn lmnt is not valid outsid th lmnt. It is found in trms of nodal prssur and nodal coordinats. Lt us considr a polynomial of th form shown in quation 2.22 Hubnr, 1975 and Zinkiwicz, Hr it is ncssary that th lmnt has six nods placd, such that ach nod corrsponds to a constant in th polynomial. In ordr to 16

33 insur accuracy of th solution, th intrpolation function is usually chosn such that th prssur is continuous across th intr lmnt boundaris compatibl formulation. a a1x + a2 y + a3x + a4 xy a5 y Each lmnt will mak a contribution to th intgral in th functional. Ths contributions must b valuatd for ach lmnt and addd togthr. Th solution to th problm is obtaind whn th nodal prssurs ar chosn such that th functional is a minimum. Th algbraic systm of quations rsulting from a minimization procss yilds th nodal prssur. In summary th stps involvd ar 1. Exprss th problm as a functional. 2. Divid th ara into lmnts. 3. Choos intrpolation functions and nods. 4. Evaluat lmnt contributions to functional. 5. Minimiz functional Finit lmnt solution of Rynolds Equation This sction dvlops th basic thory of th finit lmnt mthod for two-dimnsional thin fluid films. Th thory is gnralizd to a finit slidr baring with rlativ surfac vlocitis as wll as a squz vlocity. Th analysis applis to any baring that can b unrolld onto a horizontal surfac such as a journal baring, partial arc baring and tilt pad barings. Considr a gnral finit lmnt slidr shown in figur 2.2. Th prssur must minimiz th function in quation 2.23 nglcting th body forcs Allair t. al., J 3 ρh 2k x µ t P = P ρhu P + ρh P da + q n A C q PdC 2.23 Hr th subscripts x and y dnot th vctor componnts in th X and Y dirctions rspctivly, th functional may b writtn as show in quation Th boundary conditions along part of th boundary, C p, th prssur is spcifid as in quation 2.25 and along th rmaindr of th boundary, C q, th mass rat of flow is spcifid as shown in quations 2.26 and Th valus 17

34 18 of p a and q a may vary around th boundary. If th prssur is spcifid around th ntir boundary, C p = C, th last intgral in quation 2.24 is simply ignord Nicholas, = A y x x da P h t y P U x P U h y P x P k h P J ρ ρ µ ρ dc P n q n q q C y y x x [P] on Cp = [P a ] on Cp 2.25 [q x ] on Cq = [q ax ] on Cq 2.26 [q y ] on Cq = [q ay ] on Cq 2.27 Th fluid is dividd into a global systm of lmnts in a mannr such as that shown in figurs 2.4 and 2.6. Lt th rgion of intrst in th X-Y plan b dividd into E finit lmnts. First th baring surfac is dividd into M quadrilatrals in th y dirction and N quadrilatrals in th x dirction. This prmits th nods at th intrsctions to b labld in a rctangular matrix of ordr MxN. By xprssing th prssur and th distribution of th various forcing functions in trms of thir rspctiv nodal valus through intrpolation function, y x L i, w obtain quations 2.28 and Hr, P r and U r ar th vctors of nodal prssurs and x-componnt vlocitis rspctivly, of lmnt. [ ],,, i i i P y x L P y x L y x p r = = 2.28 [ ],, i i i x U y x L L U y x U r = = 2.29 Each of th nodal prssurs P i, P j, P k ar labld corrspondingly and th coordinats of th nods ar dnotd by x i,y i, x j,y j, x k,y k. Within ach lmnt a linar intrpolation function is rprsntd as in quation k k j j i i P y x L P y x L P y x L y x P,,,, + + = 2.30 Whr, [ ] 2, n n n n A y c x b a y x L + + = a n, b n, c n ar constants A = ara of th lmnt Sinc th cofficints a n, b n, c n and th ara of th lmnts ar known, only th nodal prssurs ar unknown in th lmnt. Th function is brokn up into intgrals ovr th individual lmnts as shown in quations 2.31 and 2.24.

35 19 = = E P J P J Th scond intgral is valuatd only for th portion of th lmnt boundary that is locatd along C q. If no portion of boundary falls along C q, thn th last intgral is simply ignord. Th prssur, shar in x-dirction, shar in y-dirction, squz and boundary flow constants ar show in quations 2.32, 2.33, 2.34, 2.35 and 2.36 Rao, 1999 and Nicholas, Thus, th lmnt functional in compact form and th total functional can thus b rprsntd as in quations 2.37 and 2.38 rspctivly Bookr and Hubnr, A m m n m x pmn da y L y L x L x L k h K + = µ ρ 2.32 A n X Uxn da x L hu K = ρ 2.33 A n y Uyn da y L hu K = ρ 2.34 = A n n h da L t h K ρ & 2.35 [ ] + = C q n y y x x n dc L n q n q q 2.36 { } [ ] { } { } { } { } { } 2 1 Uy T Ux T P T K P K P P K P P J = { } { } { } { } T h T q P K P + & 2.37 { } [ ] { } { } [ ] { } [ ] = = = = E Uy T E Ux T E P T K P K P P K P P J { } [ ] { } [] = = + E T E h T q P K P 1 1 & 2.38 Th functional must b minimizd with rspct to all of th unknown nodal prssurs. W assum that all th nodal prssurs ar unknown vn though thr ar a fw known nodal prssurs. Thus minimization of quation 2.38, as shown in quation 2.39 rsults in quation This quation is known as th finit lmnt formulation of Rynolds quation for fluid film journal barings Nicholas, This quation can b convnintly solvd in a rcognizd

36 20 matrix form as in quation Th lmnts thus obtaind can b assmbld into an ovrall systm as in = = = N J M I P P J J I,..., 1,2 1,2,..., 0,, 2.39 [ ] { } [ ] [ ] [ ] [ ] = = = = = + = E E h E Uy E Ux E P q K K K P K & 2.40 [ ] [ ] [ ] h u p F H K U K Q P K r r & r r r & = = 2.41 [ ] F P K P r r = 2.42 Whr, th ordr of matrix P K is qual to M. Applying th prssur boundary conditions can solv th quations thus obtaind. Onc th nodal prssur P r and flow U r ar found, th baring load capacity W can b computd from = S ds y x p W, 2.43

37 F Moving Plat Q U U H y P v Vlocity Profil O Stationary Plat Figur 2.1: Vlocity profil btwn two paralll plats to quantify viscosity Szri, 1999 X Y h t Squz vlocity U 2 Surfac vlocity Film thicknss, h U 1 Surfac vlocity U 1 Atmosphric prssur, P a ρhu*n=q a Mass flow rat Figur 2.2: Fluid film gomtry Allair t al.,

38 Lading dg Boundary C q whr P =0 x Tailing dg X b RӨ dirction P=0 L/2 P=0 Y b Z dirction Boundary C p whr P=0 Figur 2.3: Boundary conditions on unrolld baring surfac Allair t al., 1975 X Y N Nods M Nods 32 M x N quadrilatral lmnts Figur 2.4: Division of fluid film into lmnts Nicholas,

39 Physical Problm Chang of physical problm Mathmatical Modl Govrnd by diffrntial quations considring assumptions on Gomtry Kinmatics Matrial Law Loading Boundary conditions tc. Improv Mathmatical Modl Finit Elmnt Solution Choic of Finit Elmnts Msh Dnsity Solution Paramtrs Rprsntation of Loading Boundary conditions Rfin msh, solution paramtrs tc. Assssmnt of accuracy of finit lmnt solution of mathmatical modl Intrprtation of Rsults Rfin Analysis Dsign Improvmnts Structural optimization Figur 2.5: Th procss of finit lmnt analysis Bath,

40 Z Elmnt Ara = A nod k P k P i nod i nod j nodal prssurs P j Y X i, Y i X k, Y k X X j, Y j Linar intrpolation function Lx,y = a 0 +a 1 x+a 2 y Figur 2.6: Th th two-dimnsional simplx triangular lmnt Nicholas, Copyright Raja Shkar Balupari

41 CHAPTER THREE BEARING DESIGN PROGRAM 3.1 Introduction Th program nam BRGDS is an acronym for Baring Dsign. It is a finit lmnt basd computr program for th analysis of hydrodynamic partial arc barings and barings with hydrostatic rcsss. Dr. Robrt W. Stphnson initially dvlopd th program. Additional work has bn don to includ th capability to comput th dynamic charactristics for partial arc, full and hybrid barings. In th prsnt chaptr a dtaild ovrviw of th BRGDS program, coordinat systm, flow charts for finit lmnt mthod and convrgnc conditions ar prsntd. Th main focus is on th computation of th dynamic charactristics. 3.2 Ovrviw A standard finit lmnt implmntation is usd to modl th land rgion of th baring. Th finit lmnt modl of th baring is gnratd using four-nod quadrilatral and/or thr-nod triangular lmnts. Th finit lmnt solution of th Rynolds quation is usd to obtain th prssur ovr th finit baring surfac. Th assumptions usd in th solution of Rynolds quation in this program ar incomprssibl flow, constant lubricant viscosity, and no fluid inrtia trms. Th baring pad is assumd to b unrolld or flat, so that a two-dimnsional msh rprsntation in th axial-circumfrntial plan can b usd. Th baring dg is subjctd to th boundary conditions. Th nods hav ithr prssur or flow as th only dgr of frdom. Th systm matrics ar thus dvlopd from th baring gomtry and oprating paramtrs such as journal rotational spd, lubricant viscosity and claranc. Eithr a spcifid load or a spcifid ccntricity is usd to apply loading to th baring. For a spcifid load th program itrats on th ccntricity,, and th ccntricity angl, Ф, until load quilibrium is rachd. Howvr, for spcifid ccntricity th itrations ar prformd on ccntricity angl and th baring load is calculatd. A convntional Nwton-Raphson mthod is 25

42 usd for adjusting th journal displacmnt to find quilibrium position. Prssurs lowr than th cavitation prssur ar st qual to th spcifid cavitation prssur. Turbulnc is includd as an option and is implmntd as a corrction factor to th laminar solution. For turbulnt solution, itrations ar prformd until prssur convrgnc is achivd sinc th corrction factors in th systm quations ar prssur dpndnt. Th hydrostatic capability of th baring is includd in th analysis by spcifying ithr flow or prssur boundary conditions on th baring land around th dg of th hydrostatic rcss. Hr, all nods around th rcss priphry ar assumd to b at th rcss prssur itslf. For spcifid rcss prssur, th rcss nods at th priphry ar st to hav th boundary condition prssur and solution is achivd. Th flow into and out of th rcss is calculatd by back substitution of th prssur solution into th systm quations. In th othr cas of rcss flow bing dfind, an itrativ solution is prformd to adjust th rcss prssur, until th flow calculatd by back substitution of prssur solution agrs with th spcifid rcss flow within a spcifid tolranc. Th prssur is thn intgratd ovr th baring to giv th load capacity and momnt on th baring about a spcifid axial location. A 2X2 Gauss quadratur routin is usd to numrically intgrat th lmnts. A thr-point intgration schm is usd to numrically intgrat th triangular lmnts Zinkiwicz and Bath. Linar stiffnss and damping cofficints ar obtaind from small prturbations of displacmnt and vlocity about th quilibrium position rspctivly. In addition, th baring film thicknss, torqu loss and fluid flow into th baring from hydrodynamic action ar calculatd. Th ffcts of th axial misalignmnt of th shaft of th baring ar not considrd in th currnt implmntation. 3.3 Baring Coordinat Systm For th purpos of th program a right handd coordinatd systm fixd in spac is st up as shown in Figurs 3.1 and 3.2. X and Y ar th horizontal and vrtical coordinats of th gomtrical cntr of th shaft. Th angular coordinat, Θ, is masurd positiv in th countrclockwis dirction, starting always from th positiv x-axis in th dirction of rotation of th shaft. Th two dimnsional nodal coordinats rquird for th baring modl ar in trms of 26

43 Z axial and Θ circumfrntial dgrs. Th baring radius is spcifid as input and it is usd along with Θ to corrctly dfin th circumfrntial distanc. 3.4 Solution Procdur A block diagram flow chart of th program solution procdur is show in Figur 3.3. Th squnc of stps prformd during baring analysis by BRGDS can b outlind as follows Stphnson, Rad th baring information from th input fil and gnrat th baring finit lmnt modl. 2. Assign th prssur boundary conditions to th dg nod and rcss nods if any. 3. Assum initial valus of hydrostatic rcss prssurs, turbulnt corrction factors and journal position. 4. Solv th finit lmnt quation for unknown prssur. Comput th unknown flows by back substitution of prssur solution into systm quations. 5. Intgrat th prssur solution ovr th lmnts to comput th baring forcs and baring momnts. 6. Chck for forc convrgnc. If convrgnc conditions ar not satisfactory thn modify journal position and prform solution for unknown prssur stp 4 7. Rst any prssur blow th cavitation prssur to th spcifid cavitation prssur and prform stp Vrify prssur convrgnc for turbulnt solution if includd. If not convrgd thn updat turbulnt corrction factors and prform stp In th vnt of a hydrostatic rcss with spcifid rcss flow, chck rcss flow convrgnc. If convrgnc condition is not achivd updat rcss prssur and prform stp Display th convrgd baring solution. Prform stp 3 for mor load cass. 3.5 Gnral Systm Equations Th basic systm quations for th hydrodynamic problm in matrix form can b rprsntd as, [K p ]{p} = {q} - [K u ]{u}

44 Whr, [K p ] = systm prssur fluidity matrix {p} = nodal prssur vctor {q} = nodal flow vctor [K u ] = systm vlocity fluidity matrix {u} = nodal vlocity vctor Th systm matrics [K p ] and [K u ] ar assmbld from th individual lmnt matrics, which ar dpndnt on th film thicknss, lubricant viscosity, and lmnt gomtry Bookr t. al., Each ntry of th vctor {u} is a known valu dpndnt of th rotational spd of th journal and is qual to th rotational spd of th journal -Rω/2. Each nod contains only on unknown dgr of frdom in 3.1, ithr p or q. Not that th prssurs and flows ar scalar valus. For nods with unknown prssurs, th nt flow valus ar always zro. Dgrs of frdom ar rmovd from 3.1 at nods with spcifid prssurs. Thrfor, th right hand sid bcoms a known vctor and all unknowns in th systm ar in {p}. Th unknown flow valus at nods with spcifid prssurs can b back calculatd from th systm quations aftr th complt systm prssur solution is obtaind Eccntricity Th position of th journal in th baring is dirctly rlatd to amount of loading on th baring. Whn th journal baring is sufficintly supplid with oil and with no load acting, th journal rotats concntrically with th baring. Howvr, whn th load is applid, th journal movs to th quilibrium position or an incrasingly ccntric position, thus forming a wdg shapd oil film to support th load. As show in Figur 1.1, th distanc btwn th baring cntr and th shaft cntr is known as ccntricity. Th ratio of ccntricity to th baring claranc is dfind as ccntricity ratio, ε. If ε=0, thn thr is no load on th pad and for ε=1 th journal would touch th baring undr largr xtrnal loads Film Thicknss Calculation Th lmnts in th matrics [K p ] and [K u ] ar functions of th film thicknss, h, which in turn varis as a function of th angular location, θ, at ach nod and is calculatd from quation

AS 5850 Finite Element Analysis

AS 5850 Finite Element Analysis AS 5850 Finit Elmnt Analysis Two-Dimnsional Linar Elasticity Instructor Prof. IIT Madras Equations of Plan Elasticity - 1 displacmnt fild strain- displacmnt rlations (infinitsimal strain) in matrix form

More information

3 Finite Element Parametric Geometry

3 Finite Element Parametric Geometry 3 Finit Elmnt Paramtric Gomtry 3. Introduction Th intgral of a matrix is th matrix containing th intgral of ach and vry on of its original componnts. Practical finit lmnt analysis rquirs intgrating matrics,

More information

Finite element discretization of Laplace and Poisson equations

Finite element discretization of Laplace and Poisson equations Finit lmnt discrtization of Laplac and Poisson quations Yashwanth Tummala Tutor: Prof S.Mittal 1 Outlin Finit Elmnt Mthod for 1D Introduction to Poisson s and Laplac s Equations Finit Elmnt Mthod for 2D-Discrtization

More information

VSMN30 FINITA ELEMENTMETODEN - DUGGA

VSMN30 FINITA ELEMENTMETODEN - DUGGA VSMN3 FINITA ELEMENTMETODEN - DUGGA 1-11-6 kl. 8.-1. Maximum points: 4, Rquird points to pass: Assistanc: CALFEM manual and calculator Problm 1 ( 8p ) 8 7 6 5 y 4 1. m x 1 3 1. m Th isotropic two-dimnsional

More information

CHAPTER 1. Introductory Concepts Elements of Vector Analysis Newton s Laws Units The basis of Newtonian Mechanics D Alembert s Principle

CHAPTER 1. Introductory Concepts Elements of Vector Analysis Newton s Laws Units The basis of Newtonian Mechanics D Alembert s Principle CHPTER 1 Introductory Concpts Elmnts of Vctor nalysis Nwton s Laws Units Th basis of Nwtonian Mchanics D lmbrt s Principl 1 Scinc of Mchanics: It is concrnd with th motion of matrial bodis. odis hav diffrnt

More information

ME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002

ME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002 3.4 Forc Analysis of Linkas An undrstandin of forc analysis of linkas is rquird to: Dtrmin th raction forcs on pins, tc. as a consqunc of a spcifid motion (don t undrstimat th sinificanc of dynamic or

More information

A Propagating Wave Packet Group Velocity Dispersion

A Propagating Wave Packet Group Velocity Dispersion Lctur 8 Phys 375 A Propagating Wav Packt Group Vlocity Disprsion Ovrviw and Motivation: In th last lctur w lookd at a localizd solution t) to th 1D fr-particl Schrödingr quation (SE) that corrsponds to

More information

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches. Subjct Chmistry Papr No and Titl Modul No and Titl Modul Tag 8/ Physical Spctroscopy / Brakdown of th Born-Oppnhimr approximation. Slction ruls for rotational-vibrational transitions. P, R branchs. CHE_P8_M

More information

Division of Mechanics Lund University MULTIBODY DYNAMICS. Examination Name (write in block letters):.

Division of Mechanics Lund University MULTIBODY DYNAMICS. Examination Name (write in block letters):. Division of Mchanics Lund Univrsity MULTIBODY DYNMICS Examination 7033 Nam (writ in block lttrs):. Id.-numbr: Writtn xamination with fiv tasks. Plas chck that all tasks ar includd. clan copy of th solutions

More information

Homotopy perturbation technique

Homotopy perturbation technique Comput. Mthods Appl. Mch. Engrg. 178 (1999) 257±262 www.lsvir.com/locat/cma Homotopy prturbation tchniqu Ji-Huan H 1 Shanghai Univrsity, Shanghai Institut of Applid Mathmatics and Mchanics, Shanghai 272,

More information

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis Middl East Tchnical Univrsity Dpartmnt of Mchanical Enginring ME 43 Introduction to Finit Elmnt Analysis Chaptr 3 Computr Implmntation of D FEM Ths nots ar prpard by Dr. Cünyt Srt http://www.m.mtu.du.tr/popl/cunyt

More information

Elements of Statistical Thermodynamics

Elements of Statistical Thermodynamics 24 Elmnts of Statistical Thrmodynamics Statistical thrmodynamics is a branch of knowldg that has its own postulats and tchniqus. W do not attmpt to giv hr vn an introduction to th fild. In this chaptr,

More information

1 Isoparametric Concept

1 Isoparametric Concept UNIVERSITY OF CALIFORNIA BERKELEY Dpartmnt of Civil Enginring Spring 06 Structural Enginring, Mchanics and Matrials Profssor: S. Govindj Nots on D isoparamtric lmnts Isoparamtric Concpt Th isoparamtric

More information

Dynamic Modelling of Hoisting Steel Wire Rope. Da-zhi CAO, Wen-zheng DU, Bao-zhu MA *

Dynamic Modelling of Hoisting Steel Wire Rope. Da-zhi CAO, Wen-zheng DU, Bao-zhu MA * 17 nd Intrnational Confrnc on Mchanical Control and Automation (ICMCA 17) ISBN: 978-1-6595-46-8 Dynamic Modlling of Hoisting Stl Wir Rop Da-zhi CAO, Wn-zhng DU, Bao-zhu MA * and Su-bing LIU Xi an High

More information

5.80 Small-Molecule Spectroscopy and Dynamics

5.80 Small-Molecule Spectroscopy and Dynamics MIT OpnCoursWar http://ocw.mit.du 5.80 Small-Molcul Spctroscopy and Dynamics Fall 008 For information about citing ths matrials or our Trms of Us, visit: http://ocw.mit.du/trms. Lctur # 3 Supplmnt Contnts

More information

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis Middl East Tchnical Univrsity Dpartmnt of Mchanical Enginring ME Introduction to Finit Elmnt Analysis Chaptr 5 Two-Dimnsional Formulation Ths nots ar prpard by Dr. Cünyt Srt http://www.m.mtu.du.tr/popl/cunyt

More information

The van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012

The van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012 Th van dr Waals intraction D. E. Sopr 2 Univrsity of Orgon 20 pril 202 Th van dr Waals intraction is discussd in Chaptr 5 of J. J. Sakurai, Modrn Quantum Mchanics. Hr I tak a look at it in a littl mor

More information

Background: We have discussed the PIB, HO, and the energy of the RR model. In this chapter, the H-atom, and atomic orbitals.

Background: We have discussed the PIB, HO, and the energy of the RR model. In this chapter, the H-atom, and atomic orbitals. Chaptr 7 Th Hydrogn Atom Background: W hav discussd th PIB HO and th nrgy of th RR modl. In this chaptr th H-atom and atomic orbitals. * A singl particl moving undr a cntral forc adoptd from Scott Kirby

More information

Higher order derivatives

Higher order derivatives Robrto s Nots on Diffrntial Calculus Chaptr 4: Basic diffrntiation ruls Sction 7 Highr ordr drivativs What you nd to know alrady: Basic diffrntiation ruls. What you can larn hr: How to rpat th procss of

More information

Direct Approach for Discrete Systems One-Dimensional Elements

Direct Approach for Discrete Systems One-Dimensional Elements CONTINUUM & FINITE ELEMENT METHOD Dirct Approach or Discrt Systms On-Dimnsional Elmnts Pro. Song Jin Par Mchanical Enginring, POSTECH Dirct Approach or Discrt Systms Dirct approach has th ollowing aturs:

More information

1.2 Faraday s law A changing magnetic field induces an electric field. Their relation is given by:

1.2 Faraday s law A changing magnetic field induces an electric field. Their relation is given by: Elctromagntic Induction. Lorntz forc on moving charg Point charg moving at vlocity v, F qv B () For a sction of lctric currnt I in a thin wir dl is Idl, th forc is df Idl B () Elctromotiv forc f s any

More information

Difference -Analytical Method of The One-Dimensional Convection-Diffusion Equation

Difference -Analytical Method of The One-Dimensional Convection-Diffusion Equation Diffrnc -Analytical Mthod of Th On-Dimnsional Convction-Diffusion Equation Dalabav Umurdin Dpartmnt mathmatic modlling, Univrsity of orld Economy and Diplomacy, Uzbistan Abstract. An analytical diffrncing

More information

Addition of angular momentum

Addition of angular momentum Addition of angular momntum April, 07 Oftn w nd to combin diffrnt sourcs of angular momntum to charactriz th total angular momntum of a systm, or to divid th total angular momntum into parts to valuat

More information

u 3 = u 3 (x 1, x 2, x 3 )

u 3 = u 3 (x 1, x 2, x 3 ) Lctur 23: Curvilinar Coordinats (RHB 8.0 It is oftn convnint to work with variabls othr than th Cartsian coordinats x i ( = x, y, z. For xampl in Lctur 5 w mt sphrical polar and cylindrical polar coordinats.

More information

4.2 Design of Sections for Flexure

4.2 Design of Sections for Flexure 4. Dsign of Sctions for Flxur This sction covrs th following topics Prliminary Dsign Final Dsign for Typ 1 Mmbrs Spcial Cas Calculation of Momnt Dmand For simply supportd prstrssd bams, th maximum momnt

More information

Instantaneous Cutting Force Model in High-Speed Milling Process with Gyroscopic Effect

Instantaneous Cutting Force Model in High-Speed Milling Process with Gyroscopic Effect Advancd Matrials sarch Onlin: -8-6 ISS: 66-8985, Vols. 34-36, pp 389-39 doi:.48/www.scintific.nt/am.34-36.389 rans ch Publications, Switzrland Instantanous Cutting Forc Modl in High-Spd Milling Procss

More information

Addition of angular momentum

Addition of angular momentum Addition of angular momntum April, 0 Oftn w nd to combin diffrnt sourcs of angular momntum to charactriz th total angular momntum of a systm, or to divid th total angular momntum into parts to valuat th

More information

Supplementary Materials

Supplementary Materials 6 Supplmntary Matrials APPENDIX A PHYSICAL INTERPRETATION OF FUEL-RATE-SPEED FUNCTION A truck running on a road with grad/slop θ positiv if moving up and ngativ if moving down facs thr rsistancs: arodynamic

More information

Exam 1. It is important that you clearly show your work and mark the final answer clearly, closed book, closed notes, no calculator.

Exam 1. It is important that you clearly show your work and mark the final answer clearly, closed book, closed notes, no calculator. Exam N a m : _ S O L U T I O N P U I D : I n s t r u c t i o n s : It is important that you clarly show your work and mark th final answr clarly, closd book, closd nots, no calculator. T i m : h o u r

More information

MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design

MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design MAE4700/5700 Finit Elmnt Analysis for Mchanical and Arospac Dsign Cornll Univrsity, Fall 2009 Nicholas Zabaras Matrials Procss Dsign and Control Laboratory Sibly School of Mchanical and Arospac Enginring

More information

Einstein Equations for Tetrad Fields

Einstein Equations for Tetrad Fields Apiron, Vol 13, No, Octobr 006 6 Einstin Equations for Ttrad Filds Ali Rıza ŞAHİN, R T L Istanbul (Turky) Evry mtric tnsor can b xprssd by th innr product of ttrad filds W prov that Einstin quations for

More information

On the Hamiltonian of a Multi-Electron Atom

On the Hamiltonian of a Multi-Electron Atom On th Hamiltonian of a Multi-Elctron Atom Austn Gronr Drxl Univrsity Philadlphia, PA Octobr 29, 2010 1 Introduction In this papr, w will xhibit th procss of achiving th Hamiltonian for an lctron gas. Making

More information

Topology Optimization of Suction Muffler for Noise Attenuation

Topology Optimization of Suction Muffler for Noise Attenuation Purdu Univrsity Purdu -Pubs Intrnational Comprssor Enginring Confrnc School of Mchanical Enginring 2012 Topology Optimization of Suction Mufflr for Nois Attnuation Jin Woo L jinwool@ajou.ac.kr Dong Wook

More information

That is, we start with a general matrix: And end with a simpler matrix:

That is, we start with a general matrix: And end with a simpler matrix: DIAGON ALIZATION OF THE STR ESS TEN SOR INTRO DUCTIO N By th us of Cauchy s thorm w ar abl to rduc th numbr of strss componnts in th strss tnsor to only nin valus. An additional simplification of th strss

More information

Quasi-Classical States of the Simple Harmonic Oscillator

Quasi-Classical States of the Simple Harmonic Oscillator Quasi-Classical Stats of th Simpl Harmonic Oscillator (Draft Vrsion) Introduction: Why Look for Eignstats of th Annihilation Oprator? Excpt for th ground stat, th corrspondnc btwn th quantum nrgy ignstats

More information

Search sequence databases 3 10/25/2016

Search sequence databases 3 10/25/2016 Sarch squnc databass 3 10/25/2016 Etrm valu distribution Ø Suppos X is a random variabl with probability dnsity function p(, w sampl a larg numbr S of indpndnt valus of X from this distribution for an

More information

COHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim.

COHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim. MTH rviw part b Lucian Mitroiu Th LOG and EXP functions Th ponntial function p : R, dfind as Proprtis: lim > lim p Eponntial function Y 8 6 - -8-6 - - X Th natural logarithm function ln in US- log: function

More information

Hydrogen Atom and One Electron Ions

Hydrogen Atom and One Electron Ions Hydrogn Atom and On Elctron Ions Th Schrödingr quation for this two-body problm starts out th sam as th gnral two-body Schrödingr quation. First w sparat out th motion of th cntr of mass. Th intrnal potntial

More information

Brief Introduction to Statistical Mechanics

Brief Introduction to Statistical Mechanics Brif Introduction to Statistical Mchanics. Purpos: Ths nots ar intndd to provid a vry quick introduction to Statistical Mchanics. Th fild is of cours far mor vast than could b containd in ths fw pags.

More information

The pn junction: 2 Current vs Voltage (IV) characteristics

The pn junction: 2 Current vs Voltage (IV) characteristics Th pn junction: Currnt vs Voltag (V) charactristics Considr a pn junction in quilibrium with no applid xtrnal voltag: o th V E F E F V p-typ Dpltion rgion n-typ Elctron movmnt across th junction: 1. n

More information

EXST Regression Techniques Page 1

EXST Regression Techniques Page 1 EXST704 - Rgrssion Tchniqus Pag 1 Masurmnt rrors in X W hav assumd that all variation is in Y. Masurmnt rror in this variabl will not ffct th rsults, as long as thy ar uncorrlatd and unbiasd, sinc thy

More information

16. Electromagnetics and vector elements (draft, under construction)

16. Electromagnetics and vector elements (draft, under construction) 16. Elctromagntics (draft)... 1 16.1 Introduction... 1 16.2 Paramtric coordinats... 2 16.3 Edg Basd (Vctor) Finit Elmnts... 4 16.4 Whitny vctor lmnts... 5 16.5 Wak Form... 8 16.6 Vctor lmnt matrics...

More information

SCALING OF SYNCHROTRON RADIATION WITH MULTIPOLE ORDER. J. C. Sprott

SCALING OF SYNCHROTRON RADIATION WITH MULTIPOLE ORDER. J. C. Sprott SCALING OF SYNCHROTRON RADIATION WITH MULTIPOLE ORDER J. C. Sprott PLP 821 Novmbr 1979 Plasma Studis Univrsity of Wisconsin Ths PLP Rports ar informal and prliminary and as such may contain rrors not yt

More information

Title: Vibrational structure of electronic transition

Title: Vibrational structure of electronic transition Titl: Vibrational structur of lctronic transition Pag- Th band spctrum sn in th Ultra-Violt (UV) and visibl (VIS) rgions of th lctromagntic spctrum can not intrprtd as vibrational and rotational spctrum

More information

Finite Element Model of a Ferroelectric

Finite Element Model of a Ferroelectric Excrpt from th Procdings of th COMSOL Confrnc 200 Paris Finit Elmnt Modl of a Frrolctric A. Lópz, A. D Andrés and P. Ramos * GRIFO. Dpartamnto d Elctrónica, Univrsidad d Alcalá. Alcalá d Hnars. Madrid,

More information

Finite Strain Elastic-Viscoplastic Model

Finite Strain Elastic-Viscoplastic Model Finit Strain Elastic-Viscoplastic Modl Pinksh Malhotra Mchanics of Solids,Brown Univrsity Introduction Th main goal of th projct is to modl finit strain rat-dpndnt plasticity using a modl compatibl for

More information

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation. Lur 7 Fourir Transforms and th Wav Euation Ovrviw and Motivation: W first discuss a fw faturs of th Fourir transform (FT), and thn w solv th initial-valu problm for th wav uation using th Fourir transform

More information

2.3 Matrix Formulation

2.3 Matrix Formulation 23 Matrix Formulation 43 A mor complicatd xampl ariss for a nonlinar systm of diffrntial quations Considr th following xampl Exampl 23 x y + x( x 2 y 2 y x + y( x 2 y 2 (233 Transforming to polar coordinats,

More information

San José State University Aerospace Engineering AE 138 Vector-Based Dynamics for Aerospace Applications, Fall 2016

San José State University Aerospace Engineering AE 138 Vector-Based Dynamics for Aerospace Applications, Fall 2016 San José Stat Univrsity Arospac Enginring AE 138 Vctor-Basd Dynamics for Arospac Applications, Fall 2016 Instructor: Offic Location: Email: Offic Hours: Class Days/Tim: Classroom: Prof. J.M. Huntr E272F

More information

Outline. Thanks to Ian Blockland and Randy Sobie for these slides Lifetimes of Decaying Particles Scattering Cross Sections Fermi s Golden Rule

Outline. Thanks to Ian Blockland and Randy Sobie for these slides Lifetimes of Decaying Particles Scattering Cross Sections Fermi s Golden Rule Outlin Thanks to Ian Blockland and andy obi for ths slids Liftims of Dcaying Particls cattring Cross ctions Frmi s Goldn ul Physics 424 Lctur 12 Pag 1 Obsrvabls want to rlat xprimntal masurmnts to thortical

More information

DIFFERENTIAL EQUATION

DIFFERENTIAL EQUATION MD DIFFERENTIAL EQUATION Sllabus : Ordinar diffrntial quations, thir ordr and dgr. Formation of diffrntial quations. Solution of diffrntial quations b th mthod of sparation of variabls, solution of homognous

More information

Massachusetts Institute of Technology Department of Mechanical Engineering

Massachusetts Institute of Technology Department of Mechanical Engineering Massachustts Institut of Tchnolog Dpartmnt of Mchanical Enginring. Introduction to Robotics Mid-Trm Eamination Novmbr, 005 :0 pm 4:0 pm Clos-Book. Two shts of nots ar allowd. Show how ou arrivd at our

More information

Data Assimilation 1. Alan O Neill National Centre for Earth Observation UK

Data Assimilation 1. Alan O Neill National Centre for Earth Observation UK Data Assimilation 1 Alan O Nill National Cntr for Earth Obsrvation UK Plan Motivation & basic idas Univariat (scalar) data assimilation Multivariat (vctor) data assimilation 3d-Variational Mthod (& optimal

More information

(Upside-Down o Direct Rotation) β - Numbers

(Upside-Down o Direct Rotation) β - Numbers Amrican Journal of Mathmatics and Statistics 014, 4(): 58-64 DOI: 10593/jajms0140400 (Upsid-Down o Dirct Rotation) β - Numbrs Ammar Sddiq Mahmood 1, Shukriyah Sabir Ali,* 1 Dpartmnt of Mathmatics, Collg

More information

2008 AP Calculus BC Multiple Choice Exam

2008 AP Calculus BC Multiple Choice Exam 008 AP Multipl Choic Eam Nam 008 AP Calculus BC Multipl Choic Eam Sction No Calculator Activ AP Calculus 008 BC Multipl Choic. At tim t 0, a particl moving in th -plan is th acclration vctor of th particl

More information

Principles of Humidity Dalton s law

Principles of Humidity Dalton s law Principls of Humidity Dalton s law Air is a mixtur of diffrnt gass. Th main gas componnts ar: Gas componnt volum [%] wight [%] Nitrogn N 2 78,03 75,47 Oxygn O 2 20,99 23,20 Argon Ar 0,93 1,28 Carbon dioxid

More information

Unfired pressure vessels- Part 3: Design

Unfired pressure vessels- Part 3: Design Unfird prssur vssls- Part 3: Dsign Analysis prformd by: Analysis prformd by: Analysis vrsion: According to procdur: Calculation cas: Unfird prssur vssls EDMS Rfrnc: EF EN 13445-3 V1 Introduction: This

More information

Analysis of potential flow around two-dimensional body by finite element method

Analysis of potential flow around two-dimensional body by finite element method Vol. 7(2), pp. 9-22, May, 2015 DOI: 10.5897/JMER2014.0342 rticl Numbr: 20E80053033 ISSN 2141 2383 Copyright 2015 uthor(s) rtain th copyright of this articl http://www.acadmicjournals.org/jmer Journal of

More information

COMPUTATIONAL NUCLEAR THERMAL HYDRAULICS

COMPUTATIONAL NUCLEAR THERMAL HYDRAULICS COMPUTTIONL NUCLER THERML HYDRULICS Cho, Hyoung Kyu Dpartmnt of Nuclar Enginring Soul National Univrsity CHPTER4. THE FINITE VOLUME METHOD FOR DIFFUSION PROBLEMS 2 Tabl of Contnts Chaptr 1 Chaptr 2 Chaptr

More information

GAS FOIL BEARING ANALYSIS AND THE EFFECT OF BUMP FOIL THICKNESS ON ITS PERFORMANCE CHARACTERISTICS USING A NON-LINEAR MATRIX EQUATION SOLVER

GAS FOIL BEARING ANALYSIS AND THE EFFECT OF BUMP FOIL THICKNESS ON ITS PERFORMANCE CHARACTERISTICS USING A NON-LINEAR MATRIX EQUATION SOLVER GAS FOIL BEARING ANALYSIS AND THE EFFECT OF BUMP FOIL THICKNESS ON ITS PERFORMANCE CHARACTERISTICS USING A NON-LINEAR MATRIX EQUATION SOLVER T. Moasunp. Jamir 1)*, S. K. Kakoty 1), Karuna Kalita 1) 1)

More information

CE 530 Molecular Simulation

CE 530 Molecular Simulation CE 53 Molcular Simulation Lctur 8 Fr-nrgy calculations David A. Kofk Dpartmnt of Chmical Enginring SUNY Buffalo kofk@ng.buffalo.du 2 Fr-Enrgy Calculations Uss of fr nrgy Phas quilibria Raction quilibria

More information

SECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero.

SECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero. SETION 6. 57 6. Evaluation of Dfinit Intgrals Exampl 6.6 W hav usd dfinit intgrals to valuat contour intgrals. It may com as a surpris to larn that contour intgrals and rsidus can b usd to valuat crtain

More information

DIELECTRIC AND MAGNETIC PROPERTIES OF MATERIALS

DIELECTRIC AND MAGNETIC PROPERTIES OF MATERIALS DILCTRIC AD MAGTIC PROPRTIS OF MATRIALS Dilctric Proprtis: Dilctric matrial Dilctric constant Polarization of dilctric matrials, Typs of Polarization (Polarizability). quation of intrnal filds in liquid

More information

Complex Powers and Logs (5A) Young Won Lim 10/17/13

Complex Powers and Logs (5A) Young Won Lim 10/17/13 Complx Powrs and Logs (5A) Copyright (c) 202, 203 Young W. Lim. Prmission is grantd to copy, distribut and/or modify this documnt undr th trms of th GNU Fr Documntation Licns, Vrsion.2 or any latr vrsion

More information

Nusselt number correlations for simultaneously developing laminar duct flows of liquids with temperature dependent properties

Nusselt number correlations for simultaneously developing laminar duct flows of liquids with temperature dependent properties Journal of Physics: Confrnc Sris OPEN ACCESS Nusslt numbr corrlations for simultanously dvloping laminar duct flows of liquids with tmpratur dpndnt proprtis To cit this articl: Stfano Dl Giudic t al 2014

More information

Construction of asymmetric orthogonal arrays of strength three via a replacement method

Construction of asymmetric orthogonal arrays of strength three via a replacement method isid/ms/26/2 Fbruary, 26 http://www.isid.ac.in/ statmath/indx.php?modul=prprint Construction of asymmtric orthogonal arrays of strngth thr via a rplacmnt mthod Tian-fang Zhang, Qiaoling Dng and Alok Dy

More information

Problem Set 6 Solutions

Problem Set 6 Solutions 6.04/18.06J Mathmatics for Computr Scinc March 15, 005 Srini Dvadas and Eric Lhman Problm St 6 Solutions Du: Monday, March 8 at 9 PM in Room 3-044 Problm 1. Sammy th Shark is a financial srvic providr

More information

Preliminary Fundamentals

Preliminary Fundamentals 1.0 Introduction Prliminary Fundamntals In all of our prvious work, w assumd a vry simpl modl of th lctromagntic torqu T (or powr) that is rquird in th swing quation to obtain th acclrating torqu. This

More information

Lecture 13: Conformational Sampling: MC and MD

Lecture 13: Conformational Sampling: MC and MD Statistical Thrmodynamics Lctur 13: Conformational Sampling: MC and MD Dr. Ronald M. Lvy ronlvy@tmpl.du Contributions from Mik Andrc and Danil Winstock Importanc Sampling and Mont Carlo Mthods Enrgy functions

More information

A Sub-Optimal Log-Domain Decoding Algorithm for Non-Binary LDPC Codes

A Sub-Optimal Log-Domain Decoding Algorithm for Non-Binary LDPC Codes Procdings of th 9th WSEAS Intrnational Confrnc on APPLICATIONS of COMPUTER ENGINEERING A Sub-Optimal Log-Domain Dcoding Algorithm for Non-Binary LDPC Cods CHIRAG DADLANI and RANJAN BOSE Dpartmnt of Elctrical

More information

cycle that does not cross any edges (including its own), then it has at least

cycle that does not cross any edges (including its own), then it has at least W prov th following thorm: Thorm If a K n is drawn in th plan in such a way that it has a hamiltonian cycl that dos not cross any dgs (including its own, thn it has at last n ( 4 48 π + O(n crossings Th

More information

ph People Grade Level: basic Duration: minutes Setting: classroom or field site

ph People Grade Level: basic Duration: minutes Setting: classroom or field site ph Popl Adaptd from: Whr Ar th Frogs? in Projct WET: Curriculum & Activity Guid. Bozman: Th Watrcours and th Council for Environmntal Education, 1995. ph Grad Lvl: basic Duration: 10 15 minuts Stting:

More information

The Matrix Exponential

The Matrix Exponential Th Matrix Exponntial (with xrciss) by D. Klain Vrsion 207.0.05 Corrctions and commnts ar wlcom. Th Matrix Exponntial For ach n n complx matrix A, dfin th xponntial of A to b th matrix A A k I + A + k!

More information

Answer Homework 5 PHA5127 Fall 1999 Jeff Stark

Answer Homework 5 PHA5127 Fall 1999 Jeff Stark Answr omwork 5 PA527 Fall 999 Jff Stark A patint is bing tratd with Drug X in a clinical stting. Upon admiion, an IV bolus dos of 000mg was givn which yildd an initial concntration of 5.56 µg/ml. A fw

More information

Mor Tutorial at www.dumblittldoctor.com Work th problms without a calculator, but us a calculator to chck rsults. And try diffrntiating your answrs in part III as a usful chck. I. Applications of Intgration

More information

Machine Detector Interface Workshop: ILC-SLAC, January 6-8, 2005.

Machine Detector Interface Workshop: ILC-SLAC, January 6-8, 2005. Intrnational Linar Collidr Machin Dtctor Intrfac Workshop: ILCSLAC, January 68, 2005. Prsntd by Brtt Parkr, BNLSMD Mssag: Tools ar now availabl to optimiz IR layout with compact suprconducting quadrupols

More information

The Matrix Exponential

The Matrix Exponential Th Matrix Exponntial (with xrciss) by Dan Klain Vrsion 28928 Corrctions and commnts ar wlcom Th Matrix Exponntial For ach n n complx matrix A, dfin th xponntial of A to b th matrix () A A k I + A + k!

More information

A General Thermal Equilibrium Discharge Flow Model

A General Thermal Equilibrium Discharge Flow Model Journal of Enrgy and Powr Enginring 1 (216) 392-399 doi: 1.17265/1934-8975/216.7.2 D DAVID PUBLISHING A Gnral Thrmal Equilibrium Discharg Flow Modl Minfu Zhao, Dongxu Zhang and Yufng Lv Dpartmnt of Ractor

More information

CO-ORDINATION OF FAST NUMERICAL RELAYS AND CURRENT TRANSFORMERS OVERDIMENSIONING FACTORS AND INFLUENCING PARAMETERS

CO-ORDINATION OF FAST NUMERICAL RELAYS AND CURRENT TRANSFORMERS OVERDIMENSIONING FACTORS AND INFLUENCING PARAMETERS CO-ORDINATION OF FAST NUMERICAL RELAYS AND CURRENT TRANSFORMERS OVERDIMENSIONING FACTORS AND INFLUENCING PARAMETERS Stig Holst ABB Automation Products Swdn Bapuji S Palki ABB Utilitis India This papr rports

More information

Solution: APPM 1360 Final (150 pts) Spring (60 pts total) The following parts are not related, justify your answers:

Solution: APPM 1360 Final (150 pts) Spring (60 pts total) The following parts are not related, justify your answers: APPM 6 Final 5 pts) Spring 4. 6 pts total) Th following parts ar not rlatd, justify your answrs: a) Considr th curv rprsntd by th paramtric quations, t and y t + for t. i) 6 pts) Writ down th corrsponding

More information

COMPUTER GENERATED HOLOGRAMS Optical Sciences 627 W.J. Dallas (Monday, April 04, 2005, 8:35 AM) PART I: CHAPTER TWO COMB MATH.

COMPUTER GENERATED HOLOGRAMS Optical Sciences 627 W.J. Dallas (Monday, April 04, 2005, 8:35 AM) PART I: CHAPTER TWO COMB MATH. C:\Dallas\0_Courss\03A_OpSci_67\0 Cgh_Book\0_athmaticalPrliminaris\0_0 Combath.doc of 8 COPUTER GENERATED HOLOGRAS Optical Scincs 67 W.J. Dallas (onday, April 04, 005, 8:35 A) PART I: CHAPTER TWO COB ATH

More information

Basic Polyhedral theory

Basic Polyhedral theory Basic Polyhdral thory Th st P = { A b} is calld a polyhdron. Lmma 1. Eithr th systm A = b, b 0, 0 has a solution or thr is a vctorπ such that π A 0, πb < 0 Thr cass, if solution in top row dos not ist

More information

NONLINEAR ANALYSIS OF PLATE BENDING

NONLINEAR ANALYSIS OF PLATE BENDING NONLINEAR ANALYSIS OF PLATE BENDING CONTENTS Govrning Equations of th First-Ordr Shar Dformation thor (FSDT) Finit lmnt modls of FSDT Shar and mmbran locking Computr implmntation Strss calculation Numrical

More information

4 x 4, and. where x is Town Square

4 x 4, and. where x is Town Square Accumulation and Population Dnsity E. A city locatd along a straight highway has a population whos dnsity can b approimatd by th function p 5 4 th distanc from th town squar, masurd in mils, whr 4 4, and

More information

Introduction to Arithmetic Geometry Fall 2013 Lecture #20 11/14/2013

Introduction to Arithmetic Geometry Fall 2013 Lecture #20 11/14/2013 18.782 Introduction to Arithmtic Gomtry Fall 2013 Lctur #20 11/14/2013 20.1 Dgr thorm for morphisms of curvs Lt us rstat th thorm givn at th nd of th last lctur, which w will now prov. Thorm 20.1. Lt φ:

More information

The graph of y = x (or y = ) consists of two branches, As x 0, y + ; as x 0, y +. x = 0 is the

The graph of y = x (or y = ) consists of two branches, As x 0, y + ; as x 0, y +. x = 0 is the Copyright itutcom 005 Fr download & print from wwwitutcom Do not rproduc by othr mans Functions and graphs Powr functions Th graph of n y, for n Q (st of rational numbrs) y is a straight lin through th

More information

Pipe flow friction, small vs. big pipes

Pipe flow friction, small vs. big pipes Friction actor (t/0 t o pip) Friction small vs larg pips J. Chaurtt May 016 It is an intrsting act that riction is highr in small pips than largr pips or th sam vlocity o low and th sam lngth. Friction

More information

What are those βs anyway? Understanding Design Matrix & Odds ratios

What are those βs anyway? Understanding Design Matrix & Odds ratios Ral paramtr stimat WILD 750 - Wildlif Population Analysis of 6 What ar thos βs anyway? Undrsting Dsign Matrix & Odds ratios Rfrncs Hosmr D.W.. Lmshow. 000. Applid logistic rgrssion. John Wily & ons Inc.

More information

Thermodynamical insight on the role of additives in shifting the equilibrium between white and grey tin

Thermodynamical insight on the role of additives in shifting the equilibrium between white and grey tin hrmodynamical insight on th rol of additivs in shifting th quilibrium btwn whit and gry tin Nikolay Dmntv Dpartmnt of Chmistry, mpl Univrsity, Philadlphia, PA 19122 Abstract In this study mthods of statistical

More information

TOPOLOGY DESIGN OF STRUCTURE LOADED BY EARTHQUAKE. Vienna University of Technology

TOPOLOGY DESIGN OF STRUCTURE LOADED BY EARTHQUAKE. Vienna University of Technology Bluchr Mchanical Enginring Procdings May 2014, vol. 1, num. 1 www.procdings.bluchr.com.br/vnto/10wccm TOPOLOGY DESIG OF STRUCTURE LOADED BY EARTHQUAKE P. Rosko 1 1 Cntr of Mchanics and Structural Dynamics,

More information

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis Middl East Tchnical Univrsity Dpartmnt of Mchanical Enginring ME 4 Introduction to Finit Elmnt Analysis Chaptr 4 Trusss, Bams and Frams Ths nots ar prpard by Dr. Cünyt Srt http://www.m.mtu.du.tr/popl/cunyt

More information

Simulated Analysis of Tooth Profile Error of Cycloid Steel Ball Planetary Transmission

Simulated Analysis of Tooth Profile Error of Cycloid Steel Ball Planetary Transmission 07 4th Intrnational Matrials, Machinry and Civil Enginring Confrnc(MATMCE 07) Simulatd Analysis of Tooth Profil Error of Cycloid Stl Ball Plantary Transmission Ruixu Hu,a, Yuquan Zhang,b,*, Zhanliang Zhao,c,

More information

Introduction to Condensed Matter Physics

Introduction to Condensed Matter Physics Introduction to Condnsd Mattr Physics pcific hat M.P. Vaughan Ovrviw Ovrviw of spcific hat Hat capacity Dulong-Ptit Law Einstin modl Dby modl h Hat Capacity Hat capacity h hat capacity of a systm hld at

More information

Math 34A. Final Review

Math 34A. Final Review Math A Final Rviw 1) Us th graph of y10 to find approimat valus: a) 50 0. b) y (0.65) solution for part a) first writ an quation: 50 0. now tak th logarithm of both sids: log() log(50 0. ) pand th right

More information

Lecture 37 (Schrödinger Equation) Physics Spring 2018 Douglas Fields

Lecture 37 (Schrödinger Equation) Physics Spring 2018 Douglas Fields Lctur 37 (Schrödingr Equation) Physics 6-01 Spring 018 Douglas Filds Rducd Mass OK, so th Bohr modl of th atom givs nrgy lvls: E n 1 k m n 4 But, this has on problm it was dvlopd assuming th acclration

More information

Ch. 24 Molecular Reaction Dynamics 1. Collision Theory

Ch. 24 Molecular Reaction Dynamics 1. Collision Theory Ch. 4 Molcular Raction Dynamics 1. Collision Thory Lctur 16. Diffusion-Controlld Raction 3. Th Matrial Balanc Equation 4. Transition Stat Thory: Th Eyring Equation 5. Transition Stat Thory: Thrmodynamic

More information

Extraction of Doping Density Distributions from C-V Curves

Extraction of Doping Density Distributions from C-V Curves Extraction of Doping Dnsity Distributions from C-V Curvs Hartmut F.-W. Sadrozinski SCIPP, Univ. California Santa Cruz, Santa Cruz, CA 9564 USA 1. Connction btwn C, N, V Start with Poisson quation d V =

More information

Differential Equations

Differential Equations UNIT I Diffrntial Equations.0 INTRODUCTION W li in a world of intrrlatd changing ntitis. Th locit of a falling bod changs with distanc, th position of th arth changs with tim, th ara of a circl changs

More information

Estimation of apparent fraction defective: A mathematical approach

Estimation of apparent fraction defective: A mathematical approach Availabl onlin at www.plagiarsarchlibrary.com Plagia Rsarch Library Advancs in Applid Scinc Rsarch, 011, (): 84-89 ISSN: 0976-8610 CODEN (USA): AASRFC Estimation of apparnt fraction dfctiv: A mathmatical

More information

Numerical considerations regarding the simulation of an aircraft in the approaching phase for landing

Numerical considerations regarding the simulation of an aircraft in the approaching phase for landing INCAS BULLETIN, Volum, Numbr 1/ 1 Numrical considrations rgarding th simulation of an aircraft in th approaching phas for landing Ionl Cristinl IORGA ionliorga@yahoo.com Univrsity of Craiova, Alxandru

More information