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

Size: px
Start display at page:

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

Transcription

1 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) Dpartmnt of Mchanical Enginring, Indian Institut of Tchnology Guwahati. Guwahati , India. * Corrsponding author: sunpbokdi@gmail.com ABSTRACT Gas foil barings (GFBs) satisfy many of th ruirmnts notd for novl oil-fr turbomachinry. Howvr, GFBs hav a limitd load carrying capacity. This papr prsnts a numrical modl in ordr to assss th prformanc charactristics of gas foil barings. Th finit diffrnc schm has bn usd to discrtiz th govrning Rynolds uation and th prssur is calculatd by solving non-linar matrix uation using Nwton-Raphson tchniu. Th static prformanc analysis has bn carrid out. Th computational analysis hav bn compard with th xprimntal and thortical rsults availabl in th litratur and th ffcts of bump foil thicknss, numbr of bumps and bump complianc cofficint on th load carrying capacity at diffrnt rotor spd hav bn invstigatd. Th rsults of th study show that too thin bump foil thicknss may lad to a significant dcras in th load capacity. Howvr for accurat prdictions of th foil baring prformancs, mor dtails foil structur of 1D and 2D finit lmnt modl should b considrd. KEYWORDS Gas Foil Barings (GFBs), Nwton-Raphson tchniu, bump foil thicknss, bump foil complianc 1. INTRODUCTION Gas foil barings hav bn succssfully usd in various typs of rotating machinry, such as turbochargrs, auxiliary powr units, and gas turbin ngins. It offrs svral advantags ovr th convntional rigid surfac barings including highr load capacity, lowr powr loss and bttr stability. Ths advantagous charactristics of foil barings hav bn dmonstratd through many xprimntal invstigations [1 4]. Th foil structur of a gas foil baring consists of two parts: A smooth top foil that provids th baring surfac and a corrugatd bump foil that provids rsilint support for th top foil. Undr th action of th hydrodynamic prssur, th foil structur dforms. This yilds an apprciabl chang in th lubricant-film profil. Thrfor, th fluid film prssur must b coupld to th dformation of th foil structur in ordr to ascrtain th charactristics of th foil baring DOI : /ijmch

2 prformanc. From this point of viw, many analytical studis hav bn conductd basd on a rang of structural modls. Walowit and Anno [5] first prsntd an lastic modl for a singl bump, whr in th bnding and tnsil stiffnss of th bump wr considrd. This structural modl was usd by Hshmat t al. [6] in an invstigation of th stady stat bhavior of a foil baring. Ku and Hshmat [7, 8] also dvlopd a thortical modl of th corrugatd foil strip dformation introducing friction forc btwn th bump foils and th baring housing or top foil, and th ffct of bump gomtry on th foil strip complianc. An xprimntal procdur has bn prsntd to invstigat th foil strip dflction undr static loads. Through an optical track systm, bump dflction imags ar capturd indicating that th horizontal dflction of th sgmnt btwn bumps is ngligibl compard to th transvrsal dflction of th bumps. Rubio and San Andrés [9] also prsntd an xprimntal and analytical procdur aimd to idntify th structural stiffnss for an ntir bump-typ foil baring. Ruscitto t al. [10] prformd a sris of load capacity tsts of bump typ GFBs. Th static prformanc analysis of GFBs considring thrdimnsional shap of th foil structur was prformd by L t al. [11]. Using this modl, th dflctions of intrconnctd bumps ar compard to thos of sparatd bumps, and th minimum film thicknss is compard to thos of prvious modls. In addition, th ffcts of th top foil and bump foil thicknss on th foil baring static prformanc ar also valuatd. L t al. [12] studid oprating charactristics of th bump GFBs considring top foil bnding phnomnon and corrlation among bump foils. Kim and San Andr s [13, 14] in comparisons with tst data [10] validat a GFB modl that implmnts th simpl lastic foundation modl with formulas for bump stiffnss takn from [15]. Thy did analysis of GFBs intgrating 1D and 2D FE top foil modls. 2D FE modl prdictions ovrstimat th minimum film thicknss at th baring cntrlin, but slightly undrstimat it at th baring dgs. Prdictions from th 1D FE modl compar bst to th limitd tsts data, rproducing closly th xprimntal circumfrntial profil of minimum film thicknss. Svral approachs hav bn put forward to calculat th stady stat charactristics of gas foil barings by diffrnt rsarchrs. Kospl [16] dvlopd a forward itration mthod. Finit diffrnc mthod has bn usd for solving th fluid modl and finit lmnt mthod has bn usd for th structural modl. Hshmat t al. [6] has dvlopd an invrs mthod whr finit diffrnc mthod has bn usd for solving fluid modl and th structural modl. A Nwton- Raphson mthod has bn usd to solv th finit diffrnc mthod. A modifid forward itration mthod has bn put forwardd by Carpino and Png [17] whr finit lmnt mthod has bn usd for both fluid modl and structural modl. Subsuntly, this modifid itration mthod has bn usd by Carpino and Png [17] for calculating stady stat charactristics of gas foil barings whil calculating stiffnss and damping cofficints for lastically supportd gas foil barings. L t al. [11] dvlopd an itrativ mthod for calculating th static prformanc analysis of foil journal barings. A thr dimnsional modl has bn considrd for th foil structur and finit lmnt mthod has bn usd for solving fluid modl and th thr dimnsional foil modl. In this papr finit diffrnc mthod has bn usd to discrtiz th govrning Rynolds uation. Th discrtizd uation has bn writtn in th form of a uadratic function for vry nod and Nwton-Raphson tchniu hav bn usd to solv th non-linar matrix uation. Th stady stat charactristic of th GFB ar thn assssd and th stady stat charactristics rsults hav bn compard with th xprimntal and thortical rsults availabl in th litratur. Th 16

3 ffcts of various baring paramtrs on th load carrying capacity of th baring hav bn nvisagd 2. GAS FOIL BEARING DESCRIPTION 2.1 Govrning Euations Figur 1 shows a schmatic viw of bump typ GFB. Th Rynolds uation dscribs th gnration of th gas prssur (p) within th film thicknss (h) and for an isothrmal, isoviscous idal gas this uation is, 3 p 3 p ( ph) ( ph) ph ph + = 6µω R + 12µ x x z z x t (1) whr, ( x, z ) ar th circumfrntial and axial coordinats on th plan of baring,. Th prssur taks ambint valu p on th sid boundaris of th baring. a Fig.1 Schmatic viw of bump typ GFB Th boundary conditions for th solution of En. (1) ar p = pa at θ = 0 and 2π p = p at z = 0 and L (2) a Th film thicknss ( h ) for a prfctly alignd journal configuration is givn by h = c + cos( θ φ) + wt (3) Whr,c, and wt ar th assmbld claranc, journal ccntricity and th lastic dformation of th foil structur rspctivly. 17

4 Thrfor En. (1) and En. (3) ar coupld through th hydrodynamic prssur and lastic dformation of th foil structur as a function of th hydrodynamic prssur. Th abov govrning uations can b normalizd by using th following substitutions as follows: α p x S = a H = h z = Z = θ C C R R Y p w y = P = ε = W = t L p C C a (4) Hnc stady stat non-dimnsional Rynolds uation as givn in En. (1) is givn by whr 3 P 3 P PH PH PH + Λ = 0 Z Z θ θ θ Λ 6µω = p a ( R ) 2 C (5) And similarly th non-dimnsional film thicknss is givn as, H = 1 + ε cos( θ φ ) + W (6) whr ε is th ccntricity ratio. 2.2 Stady Stat Formulation of Rynolds Euation Th non-dimnsionalizd Rynolds En. (5) has bn discrtizd using a finit diffrnc approach. For all th drivativs prsnt in En. (5) cntral diffrnc formula hav bn usd. Figur 2 shows an xfoliatd viw of a baring showing th msh siz ( θ Z ). m and n ar th divisions along θ and Z rspctivly. Th dimnsion of th prssur in th domain in matrix m + 1 n + 1 rprsnt by P. Th vctor rprsntation of P is p and its lngth form will b ( ) ( ) is( m 1)( n 1) + +. Similarly foil dflction and film thicknss in matrix forms ar rprsntd by W and H rspctivly. If (i, j) is any arbitrary point in th matrix rprsntation, its position in th = j 1 m i. Th four adjacnt nods associatd vctor rprsntation can b writtn as ( )( ) with nod (i, j) ar ( i 1, j), ( i 1, j), ( i, j 1) and ( i, j 1) + + as shown in Fig 2. Th discrtizd form of uation (5) for (i, j) th nod or th nod can b writtn in a uadratic form as ( ) T p A p + b p = 0 (7) whr = { P( i, j 1) P( i 1, j) P( i, j) P( i+ 1, j) P( i, j+ 1) } p and its siz is 5 1, T 18

5 A H 3 3 ( i, j 1) ( i, j 1) ( Z ) 4( Z) 3 3 H( i 1, j) H( i 1, j) ( θ ) 4( θ ) H H 2H 2H H H ( i, j) ( i, j) ( i, j) ( i, j) ( i, j) ( i, j) = ( Z) ( θ ) ( θ ) ( Z) ( θ ) ( Z) 3 3 H( i+ 1, j) H( i+ 1, j) ( θ ) 4( θ ) 3 3 H( i, j+ 1) H(, 1) i j + 2 4( 2 Z) 4( Z ) H A is a 5 5 non-symmtric matrix, ΛH ΛH ( i 1, j) ( i+ 1, j) b = 0 0 0, 2( θ ) 2( θ ) b is a row vctor of siz 1 5. Fig. 2: An xfoliatd viw of a baring showing th msh siz ( θ Z ) Th forc of th gas film acting on th journal can b computd by intgrating th prssur ovr th baring surfac. According to th coordinat systm illustratd in Fig. 3, this intgration can b writtn as: L D 2π Wx Pcosθdθ dz = L D 0 19

6 L D 2π Wy PsinθdθdZ = (8) L D 0 For numrical intgration, Simpson s on third rul has bn usd. Finally th total non-dimnsional load is givn by 2 2 W = Wx + Wy (9) dθ θ 180-dθ 1. Y P sinθ rdθdz P rdθdz -P cosθ rdθdz X Fig. 3: Th coordinat systm and th sign convntion of th journal forcs 3. SOLUTION METHOD In th prsnt algorithm, a fixd coordinat systm has bn considrd (X and Y). Thrfor, for stady stat uilibrium, th horizontal load componnt should bcom zro. Th flowchart of th solution procss is shown in Fig. 4. For a givn ccntricity ratio, th attitud angl is varid till th horizontal load componnt approximatly bcoms zro. Bisction mthod has bn usd for brackting and also to calculat th corrct attitud angl. Th load capacity is simply ual to th vrtical load componnt. Initially, th non-dimnsional top foil dflction W is assumd as zro. To start with th itration, th non-dimntionalisd prssurs at all th msh points ar assumd as 1 (i.. ambint prssur) and thos at th boundaris ar also st to th ambint prssur. Nwton-Raphson mthod has bn usd to solv th discrtizd Rynolds En. (7) for 20

7 air prssur p in this papr. Onc th prssur is found out th non-dimnsional top foil dflction W is calculatd. Aftr that, substituting th nw valu of non-dimnsional top foil dflction W in th uation of film thicknss, th uadratic uation is solvd for all th msh points to stimat th prssur at all ths points. This uadratic uation would not b satisfid for th prssurs which ar assumd to b constant in th bginning. Hnc, th itrativ procss is carrid out until th following convrgnc critrion is satisfid. 3.1 Nwton-Raphson Mthod ( k 1) ( k ) ( p ) p ( ) ( k ) ( p ) 10 In ordr to apply Nwton-Raphson mthod w dfin th rsidu vctor r in En. (7) as ( T ) whr 1,2,..., ( 1)( 1) r = p A p + b p = m + n + (11) whr r is th rsidu vctor for th nod and its lngth is 5 1. Th lngth of th rsidu vctor r is ( m 1)( n 1) ( ) ( ) ( ) + +. Th rsidu vctor of th th nod for th k-th itration is ( k ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) whr 1,2,..., ( 1)( 1) r k T k k k k = + = m + n + p A p b p (12) Th th row of th snsitivity matrix, G is calculatd as ( ) ( k ) ( ) ( k) 6 { } ( k ) ( ) ( ) whr 1,2,..., ( 1)( 1) G T ( ) ( k ) T k = + ( ) + = m + n + p A A b (13) Th siz of th matrix G is ( m + 1)( n + 1) ( m + 1)( n + 1) ( k + 1) ( k ) ( k ) ( 1) ( ) ( k ) k. Th prssur is updatd as p = p λ G r (14) whr ( k +1) p and ( k ) p ar th prssurs at ( k + 1) th and k th itrations rspctivly. ( k) G (10) is th snsitivity matrix at k th ( k) itration. λ is th k th stp lngth. Th computation is stoppd whn th 6 rror btwn two succssiv itration stps is smallr than

8 Fig. 4: Flowchart of solution algorithm 4. MODELLING OF FOIL SUPPORT STRUCTURE Th lastic support structur considrd in th prsnt analysis is a simpl foundation modl, th original work of Hshmat t al [6] which most of th publishd modls ar basd on. This analysis rlis on svral assumptions which othr rsarchrs hav also rproduc: 22

9 (1) Th stiffnss of a bump strip is uniformly distributd throughout th baring surfac, i.. th bump strip is rgardd as a uniform lastic foundation. (2) A bump stiffnss is constant, indpndnt of th actual bump dflction, not rlatd to or constraind by adjacnt bumps. (3) Th top foil dos not sag btwn adjacnt bumps. Th top foil dos not hav ithr bnding or mmbran stiffnss, and its dflction follows that of th bump With this considrations, th dformation ( w t ) dpnds on th bump complianc (α ) and th avrag prssur across th baring width as shown Fig. 5, w = α( p p ) (15) t whr, p and α ar th arithmtic man prssur in th axial dirction and complianc of th bump foil. a Fig.5 Singl sgmnt of Bump foil [18] Th dimnsionlss dflction of th bump is givn by W = S( P 1) (16) Coupling of th simpl modl uation (16) with th solution of Rynolds uation (5) is straightforward for th prdiction of th static prformanc of GFBs. 5. VALIDATION, RESULTS AND PREDICTIONS OF THE MODEL 5.1. Validation Th validity of th prsnt analysis and computational program ar assssd by comparison of prdictions with publishd data availabl in opn litratur. As shown in figur 6 and 7, th 23

10 prsnt computation analysis ar assssd with Hshmat t al [6] and Png and Carpino[5]. From th comparison, it has bn obsrvd that th prsnt rsults ar fairly in good agrmnt with thos from th rfrncs. Not that, a baring rducs to a ordinary gas baring for complianc cofficint S=0. Figur 8 shows th prssur profil of th GFB for L/D =1.0, Λ=1.0, S=1 and ε = 0.6. Fig. 6: Validation of load carrying capacity with Hshmat t al [6] and Png and Carpino[5]. 24

11 Fig. 7: Validation of attitud angl with Hshmat t al [6] and Png and Carpino[5]. Fig 8: Prssur Profil of GFB for L/D =1.0, Λ=1.0, S=1 and ε = Configuration of Tst GFB Tabl 2 provids th paramtrs for th tst foil baring givn in [10] and Figur 5 dpicts th configuration of th bump foil strip. Th foil baring is a first gnration typ with on 360º top foil and on bump strip layr, both mad of Inconl X-750. Th top foil and bump layr ar spot wldd at on nd to th baring slv. Th othr nd of th top foil is fr as wll as th nd of th bump strip layr. Th journal rotational dirction is from th fr nd of th top foil towards its fixd nd and all tsts in wr prformd with air at ambint condition. Tabl 2: Dsign dtails of foil baring, rfrnc [10] Baring radius, R=D/2 Baring lngth, L Foil arc circumfrntial lngth, l x mm 38.1 mm 120 mm Radial journal travl, c ( ~ Claranc, C ) 31.8 µm J Top foil thicknss t t µm Bump foil thicknss, t b µm Bump pitch, s mm Half bump lngth, l mm Bump hight, h b mm Numbr of Bumps 26 Bump foil Young s modulus, E 214 GPa b 25

12 Bump foil Poisson s ratio, υ Rsults Effct of Bump Foil Thicknss on th Static Prformancs Th GFB computational tools for simpl lastic foundation modl hav bn dvlopd which prdict th static prformanc of th GFB. Th prdiction for th simpl lastic modl uss a msh of 90 and 10 lmnts in th circumfrntial and axial dirctions rspctivly. Th sam msh siz wr usd for th finit diffrnc numrical schm solving Rynolds En. (5) and calculating th hydrodynamic gas film prssur. Th bump foils which is a masur of th foil structur has bn varis and its charactristics has bn obsrvd on th prformanc of load capacity of th baring. (a) 26

13 (b) Fig 9: Effct of bump foil thicknss on th load carrying capacity of GFBs at diffrnt ccntricity whn th rotating vlocity is (a) 30,000 rpm and (b) 45,000 rpm Figur 9 (a) and (b) shows th prdiction of th load carrying capacity for various bump foil thicknss at diffrnt ccntricity ratios and at th rotating vlocity of 30,000 rpm and 45,000 rpm rspctivly. As th ccntricity incrass, th load capacity incras linarly for all valus of bump foil thicknss, which bcoms significant at highr ccntricity ratios and at highr rpm. Th incras though is not so prominnt for smallr bump foil thicknss at tb = 0.05µ m compard to highr foil thicknss at tb = 0.20 µ m oprating at diffrnt spd which is minnt from figur 9 (a) and (b). Ths rsults hnc indicat that too thin a bump foil thicknss can lad to a significant dcras in th load capacity du to th xcssiv foil structur dflction. In addition, a larg incrmnt of th load capacity from variation of th foil stiffnss may not b xpctd undr a simpl foil structur Effct of numbr of bumps and bump complianc cofficint on th load carrying capacity of GFBs Fig 10: Numbr of bumps vrsus load carrying capacity of GFBs 27

14 Fig 11: Complianc cofficint vrsus load carrying capacity of GFBs Figur 10 dpicts th variation of th load capacity of th baring as th numbr of bump foils varis. Th numbr of bumps ar calculatd by dividing th foil arc circumfrntial lngth ( l x ) by th bump pitch (s), and as th bump numbr incrass th bumps bcom mor crowdd and stiffr which rsultd in an incras in th load capacity of th baring as dpictd from fig. 10. Th ffct of complianc cofficint (S) which is a function of th bump complianc, on th load capacity has bn portrait in figur 11. Th load capacity gnrally dcrass as th complianc cofficint incrass, which prdicts that as th bumps complianc incrass th stiffnss tnds to dcras, as bump complianc is th invrs of stiffnss hnc a dcras in th load prformancs of th baring is minnt. 6. CONCLUSION In this study, a simpl lastic foundation modl of th foil has bn considrd whr in th stiffnss of th bump strip ar uniformly distributd throughout th baring surfac, i.. th bump strip is rgardd as a uniform lastic foundation and constant, indpndnt of th actual bump dflction and ar not rlatd to or constraind by adjacnt bumps. Also th top foil dos not sag btwn adjacnt bumps and dos not hav ithr bnding or mmbran stiffnss and its dflction follows that of th bump. A numrical modl is thn dvlopd in ordr to find out th prformanc charactristics of th gas foil barings. Th finit diffrnc schm has bn usd to discrtiz th govrning Rynolds uation and th prssurs ar calculatd by solving a nonlinar matrix uation using Nwton-Raphson tchniu and th ffcts of bump foil thicknss, numbr of bumps and bump complianc cofficint on th foil baring static prformancs wr valuatd. Th rsults of th study show that too thin bump foil thicknss may lad to a significant dcras in th load capacity. Howvr for accurat prdictions of th foil baring prformancs, th prsnt computational analysis can b furthr xtndd by incorporating 1D and 2D finit lmnt modl of th foil considring th dflction of th top foil as wll as th bump foils. 28

15 REFERENCES [1] Hshmat, H., Shapiro, W., and Gray, S., 1982, Dvlopmnt of Foil Journal Barings for High Load Capacity and High Spd Whirl Stability, ASME J. Lubr. Tchnol., 104, pp [2] Hshmat, H., 1994, Advancmnts in th Prformanc of Arodynamic Foil Journal Barings: High Spd and Load Capability, ASME J. Tribol., 116, pp [3] Hshmat, H., 2000, Opration of Foil Barings Byond th Bnding Critical Mod, ASME J. Tribol., 122, pp [4] L, Y. B., Kim, T. H., Kim, C. H., L, N. S., and Choi, D. H., 2004, Dynamic Charactristics of a Flxibl Rotor Systm Supportd by a Vis-colastic Foil Baring, Tribol. Int., 37, pp [5] Walowit, J. A., and Anno, J. N., 1975, Modrn Dvlopmnt of Lubrication Mchanics Applid Scinc, London, Chap. 7. [6] Hshmat, H., Walowit, J. A., and Pinkus, O., 1983, Analysis of Gas Lubricatd Foil Journal Barings ASME J. Lubr. Tchnol., 105, pp [7] Ku, C. P., and Hshmat, H. 1992, Complaint Foil Baring Structural Stiffnss Analysis Part I: Thortical Modl Including Strip and Variabl Bump Foil Gomtry ASME J. Tribol pp [8] Ku, C. P., and Hshmat, H., 1993, Complaint Foil Baring Structural Stiffnss Analysis Part II: Exprimntal Invstigation ASME J. Tribol. 113, pp [9] Rubio, D., and San Andrés, L., 2004, Bump-Typ Foil Baring Structural Stiffnss: Exprimnts and Prdictions ASME Papr GT [10] Ruscitto, D., McCormick, J., and Gray, S., 1978, Hydrodynamic Air Lubricatd Compliant Surfac Baring for an Automotiv Gas Turbin Engin I-Journal Baring Prformanc. NASA CR [11] L, D., Kim, Y., and Kim. K., 2008, Th Static Prformanc Analysis of foil journal barings considring thr-dimnsional shap of th foil structur ASME J. of Tribol. 130/ [12] L, Y., Park, D., Kim, C., Kim, S., 2008, Oprating charactristics of th bump foil journal barings with top foil bnding phnomnon and corrlation among bump foils Tribol. Int., 41, pp [13] Kim, T. H., and San Andrés, L., 2008, Havily Loadd Gas Foil Barings: A Modl Anchord to Tst Data ASME J.Eng Gas Turb Powr pp [14] San Andrés, L., and Kim, T., 2009, Analysis of gas foil barings intgrating FE top foil modls. Tribol. Int., 42, pp [15] Iordanoff, I., 1999, Analysis of an Arodynamic Compliant Foil Thrust Baring: Mthod for a Rapid Dsign J. Tribol., 121, pp [16] Kopsl, W., 1977, Gas Lubricatd Foil Baring Dvlopmnt for Advancd Turbomachins, U. S. Air Forc Aro Propulsion Laboratory, Rp. AFAPL-TR , Vol. 1, 2. [17] Carpino, M., and Png, J. P, 1993, Thortical Prformanc Foil Journal Barings, prsntd at th AIAA/SAE/ASME/ASME. 27th Joint Propulsion Confrnc, papr no. AIAA [18] Andr San Luis, Kim Ho Ta, 2008, Forcd nonlinar rspons of gas foil baring supportd rotors Tribology Intrnational, 41, pp NOMENCLATURE C Baring radial claranc (m) D Diamtr of journal (m) Baring ccntricity (m) 2 E Young s modulus for bump foil ( N/m ) b 2 E Young s modulus for top foil ( N/m ) t h H i, j Film thicknss (m) Non-dimnsional minimum film thicknss Grid location in circumfrntial and axial dirctions of FDM msh 29

16 3 K Bump foil structural stiffnss pr unit ara ( N/m ) l 0 f Half bump lngth (m) L Baring lngth (m) m Numbr of divisions along j dirction of FDM msh n Numbr of divisions along j dirction of FDM msh O Cntr of baring O ' Cntr of journal p 2 Hydrodynamic prssur in gas film ( N/m ) 2 p a Atmosphric prssur ( N/m ) p 2 Arithmtic man prssur along baring lngth ( N/m ) P Non-dimnsional hydrodynamic prssur R Radius of journal (m) s Bump foil pitch (m) S pa Complianc cofficint of bump foil : CK t b t t w t W W x, y, z Z α Bump foil thicknss (m) Top foil thicknss (m) Top foil transvrs dflction (m) Non-dimnsional top foil transvrs dflction Non-dimnsional stady stat load carrying capacity Coordinat systm on th plan of baring Non-dimnsional axial coordinat of baring : z R 3 1 Complianc of th bump foil ( m /N ) : ε Eccntricity ratio Λ 2 6µω R Baring numbr : pa C µ 2 Gas viscosity ( N-s/m ) φ Attitud angl (rad) θ Angular coordinat of baring (rad) : x / R υ Poisson s ratio ω Rotor angular vlocity ( rad/s ) θ, Z Non-dimnsional msh siz of FDM msh f K f 30

17 AUTHORS 1. T. Moasunp Jamir Rsarch Scholar (PhD) Dpartmnt of Mchanical Enginring Indian Institut of Tchnology Guwahati [ sunpbokdi@gmail.com ; t.jamir@iitg.rnt.in 2. Dr. S. K. Kakoty Profssor Dpartmnt of Mchanical Enginring Indian Institut of Tchnology Guwahati [sashin@iitg.rnt.in] 3. Dr. Karuna Kalita Assistant Profssor Dpartmnt of Mchanical Enginring Indian Institut of Tchnology Guwahati [karuna.kalita@iitg.rnt.in] 31

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

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

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

Dynamic response of a finite length euler-bernoulli beam on linear and nonlinear viscoelastic foundations to a concentrated moving force

Dynamic response of a finite length euler-bernoulli beam on linear and nonlinear viscoelastic foundations to a concentrated moving force Journal of Mchanical Scinc and Tchnology 2 (1) (21) 1957~1961 www.springrlink.com/contnt/1738-9x DOI 1.17/s1226-1-7-x Dynamic rspons of a finit lngth ulr-brnoulli bam on linar and nonlinar viscolastic

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

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

Effect of Aspect Ratio on the Performance of Hydrodynamic Journal Bearing Operating Under Wear

Effect of Aspect Ratio on the Performance of Hydrodynamic Journal Bearing Operating Under Wear Intrnational ournal of Thortical and Applid Mchanics. ISSN 0973-6085 Volum 12, Numbr 3 (2017) pp. 497-522 Rsarch India Publications http://.ripublication.com Effct of Aspct Ratio on th Prformanc of Hydrodynamic

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

Dynamic Characteristics Analysis of Blade of Fan Based on Ansys

Dynamic Characteristics Analysis of Blade of Fan Based on Ansys Powr and Enrgy Enginring Confrnc 1 Dynamic Charactristics Analysis of Blad of Fan Basd on Ansys Junji Zhou, Bo Liu, Dingbiao Wang, Xiaoqian li School of Chmical Enginring Zhngzhou Univrsity Scinc Road

More information

A Comparative study of Load Capacity and Pressure Distribution of Infinitely wide Parabolic and Inclined Slider Bearings

A Comparative study of Load Capacity and Pressure Distribution of Infinitely wide Parabolic and Inclined Slider Bearings Procdings of th World Congrss on Enginring 2 Vol II WCE 2, Jun 3 - July 2, 2, London, U.K. A Comparativ study of Load Capacity and Prssur Distribution of Infinitly wid Parabolic and Inclind Slidr Barings

More information

Koch Fractal Boundary Single feed Circularly Polarized Microstrip Antenna

Koch Fractal Boundary Single feed Circularly Polarized Microstrip Antenna 1 Journal of Microwavs, Optolctronics and Elctromagntic Applications, Vol. 6, No. 2, Dcmbr 2007 406 Koch Fractal Boundary Singl fd Circularly Polarizd Microstrip Antnna P. Nagswara Rao and N. V. S.N Sarma

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

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

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

Ultimate lateral load resistance of laterally loaded pile

Ultimate lateral load resistance of laterally loaded pile Ultimat latral load rsistanc of latrally loadd pil Md. M. Rahman Assistant Profssor, Dpartmnt of Civil Enginring, RUET, Rajshahi, Bangladsh Md. R. arim, A. L. Baki & D.. Paul Lctr, Dpartmnt of Civil Enginring,

More information

Elastic Analysis of Functionally Graded Variable Thickness Rotating Disk by Element Based Material Grading

Elastic Analysis of Functionally Graded Variable Thickness Rotating Disk by Element Based Material Grading Journal of Solid Mchanics ol. 9, No. 3 (017) pp. 650-66 Elastic Analysis of Functionally Gradd ariabl hicknss Rotating Disk by Elmnt Basd Matrial Grading A.K. hawait 1,*, L. Sondhi 1, Sh. Sanyal, Sh. Bhowmick

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

Rotor Stationary Control Analysis Based on Coupling KdV Equation Finite Steady Analysis Liu Dalong1,a, Xu Lijuan2,a

Rotor Stationary Control Analysis Based on Coupling KdV Equation Finite Steady Analysis Liu Dalong1,a, Xu Lijuan2,a 204 Intrnational Confrnc on Computr Scinc and Elctronic Tchnology (ICCSET 204) Rotor Stationary Control Analysis Basd on Coupling KdV Equation Finit Stady Analysis Liu Dalong,a, Xu Lijuan2,a Dpartmnt of

More information

843. Efficient modeling and simulations of Lamb wave propagation in thin plates by using a new spectral plate element

843. Efficient modeling and simulations of Lamb wave propagation in thin plates by using a new spectral plate element 843. Efficint modling and simulations of Lamb wav propagation in thin plats by using a nw spctral plat lmnt Chunling Xu, Xinwi Wang Stat Ky Laboratory of Mchanics and Control of Mchanical Structurs aning

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

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

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

INFLUENCE OF GROUND SUBSIDENCE IN THE DAMAGE TO MEXICO CITY S PRIMARY WATER SYSTEM DUE TO THE 1985 EARTHQUAKE

INFLUENCE OF GROUND SUBSIDENCE IN THE DAMAGE TO MEXICO CITY S PRIMARY WATER SYSTEM DUE TO THE 1985 EARTHQUAKE 13 th World Confrnc on Earthquak Enginring Vancouvr, B.C., Canada August 1-6, 2004 Papr No. 2165 INFLUENCE OF GROUND SUBSIDENCE IN THE DAMAGE TO MEXICO CITY S PRIMARY WATER SYSTEM DUE TO THE 1985 EARTHQUAKE

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

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

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

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

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

Ultimate strength analysis & design of residential slabs on reactive soil

Ultimate strength analysis & design of residential slabs on reactive soil Ultimat strngth analysis & dsign of rsidntial slabs on ractiv soil This documnt prsnts an ovrviw of thory undrlying ultimat strngth analysis and dsign of stiffnd raft and waffl raft slabs, as commonly

More information

Effects of Electron Model on Three-Grid Ion Engine Analyses

Effects of Electron Model on Three-Grid Ion Engine Analyses Effcts of Elctron Modl on Thr-Grid Ion Engin Analyss IEPC-2011-205 Prsntd at th 32nd Intrnational Elctric Propulsion Confrnc, Wisbadn Grmany Takshi Miyasaka 1 and Katsuo Asato 2 Gifu Univrsity, Gifu, 501-1193,

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

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

The Autonomous Underwater Vehicle (AUV) MAYA: General Description

The Autonomous Underwater Vehicle (AUV) MAYA: General Description Introduction h ocans and rivrs always hav bn and still ar an important sourc of rvnu and prosprity for mankind. Du to th grat importanc of ocans and rivrs, th scintific community maks us of Autonomous

More information

Procdings of IC-IDC0 ( and (, ( ( and (, and (f ( and (, rspctivly. If two input signals ar compltly qual, phas spctra of two signals ar qual. That is

Procdings of IC-IDC0 ( and (, ( ( and (, and (f ( and (, rspctivly. If two input signals ar compltly qual, phas spctra of two signals ar qual. That is Procdings of IC-IDC0 EFFECTS OF STOCHASTIC PHASE SPECTRUM DIFFERECES O PHASE-OLY CORRELATIO FUCTIOS PART I: STATISTICALLY COSTAT PHASE SPECTRUM DIFFERECES FOR FREQUECY IDICES Shunsu Yamai, Jun Odagiri,

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

Sara Godoy del Olmo Calculation of contaminated soil volumes : Geostatistics applied to a hydrocarbons spill Lac Megantic Case

Sara Godoy del Olmo Calculation of contaminated soil volumes : Geostatistics applied to a hydrocarbons spill Lac Megantic Case wwwnvisol-canadaca Sara Godoy dl Olmo Calculation of contaminatd soil volums : Gostatistics applid to a hydrocarbons spill Lac Mgantic Cas Gostatistics: study of a PH contamination CONTEXT OF THE STUDY

More information

A New Approach to the Fatigue Life Prediction for Notched Components Under Multiaxial Cyclic Loading. Zhi-qiang TAO and De-guang SHANG *

A New Approach to the Fatigue Life Prediction for Notched Components Under Multiaxial Cyclic Loading. Zhi-qiang TAO and De-guang SHANG * 2017 2nd Intrnational Conrnc on Applid Mchanics, Elctronics and Mchatronics Enginring (AMEME 2017) ISBN: 978-1-60595-497-4 A Nw Approach to th Fatigu Li Prdiction or Notchd Componnts Undr Multiaxial Cyclic

More information

Rational Approximation for the one-dimensional Bratu Equation

Rational Approximation for the one-dimensional Bratu Equation Intrnational Journal of Enginring & Tchnology IJET-IJES Vol:3 o:05 5 Rational Approximation for th on-dimnsional Bratu Equation Moustafa Aly Soliman Chmical Enginring Dpartmnt, Th British Univrsity in

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

Modeling of the Plasma Flow in High-Power TAL

Modeling of the Plasma Flow in High-Power TAL 39th AIAA/ASME/SAE/ASEE Joint Propulsion Confrnc and Exhibit 20-23 July 2003, Huntsvill, Alabama AIAA 2003-4701 AIAA-2003-4701 Modling of th Plasma Flow in High-Powr TAL Michal Kidar, Iain D. Boyd and

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

Two Products Manufacturer s Production Decisions with Carbon Constraint

Two Products Manufacturer s Production Decisions with Carbon Constraint Managmnt Scinc and Enginring Vol 7 No 3 pp 3-34 DOI:3968/jms9335X374 ISSN 93-34 [Print] ISSN 93-35X [Onlin] wwwcscanadant wwwcscanadaorg Two Products Manufacturr s Production Dcisions with Carbon Constraint

More information

Evaluating Reliability Systems by Using Weibull & New Weibull Extension Distributions Mushtak A.K. Shiker

Evaluating Reliability Systems by Using Weibull & New Weibull Extension Distributions Mushtak A.K. Shiker Evaluating Rliability Systms by Using Wibull & Nw Wibull Extnsion Distributions Mushtak A.K. Shikr مشتاق عبذ الغني شخير Univrsity of Babylon, Collg of Education (Ibn Hayan), Dpt. of Mathmatics Abstract

More information

2017 Water Reactor Fuel Performance Meeting September 10 (Sun) ~ 14 (Thu), 2017 Ramada Plaza Jeju Jeju Island, Korea

2017 Water Reactor Fuel Performance Meeting September 10 (Sun) ~ 14 (Thu), 2017 Ramada Plaza Jeju Jeju Island, Korea 2017 Watr Ractor Ful Prformanc Mting Sptmbr 10 (Sun) ~ 14 (Thu), 2017 Ramada Plaza Jju Jju Island, Kora Study of Ful Rod Bhavior with Missing Pllt Surfac Dfct Zhnhai Liu 1, Yi Zhou 1, Ping Chn 1, Yuanming

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

MEASURING HEAT FLUX FROM A COMPONENT ON A PCB

MEASURING HEAT FLUX FROM A COMPONENT ON A PCB MEASURING HEAT FLUX FROM A COMPONENT ON A PCB INTRODUCTION Elctronic circuit boards consist of componnts which gnrats substantial amounts of hat during thir opration. A clar knowldg of th lvl of hat dissipation

More information

Learning Spherical Convolution for Fast Features from 360 Imagery

Learning Spherical Convolution for Fast Features from 360 Imagery Larning Sphrical Convolution for Fast Faturs from 36 Imagry Anonymous Author(s) 3 4 5 6 7 8 9 3 4 5 6 7 8 9 3 4 5 6 7 8 9 3 3 3 33 34 35 In this fil w provid additional dtails to supplmnt th main papr

More information

A Novel FE Bulk-Flow Model Improving Predictions of Force Coefficients in Off-Centered Grooved Oil Seals

A Novel FE Bulk-Flow Model Improving Predictions of Force Coefficients in Off-Centered Grooved Oil Seals Txas A&M Univrsity Mchanical Enginring Dpartmnt Turbomachinry laboratory A Novl FE Bulk-Flow Modl Improving Prdictions of Forc Cofficints in Off-Cntrd Groovd Oil Sals Rsarch Progrss Rport to th Turbomachinry

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

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

ANALYTICAL MODEL FOR CFRP SHEETS BONDED TO CONCRETE

ANALYTICAL MODEL FOR CFRP SHEETS BONDED TO CONCRETE ANALYTICAL MODEL FOR CFRP SHEETS BONDED TO CONCRETE Brian Millr and Dr. Antonio Nanni Univrsity of Missouri Rolla Dpartmnt of Civil Enginring 5 ERL 1870 Minr Circl Rolla, MO 65401, USA Dr. Charls E. Bakis

More information

NUMERICAL SIMULATION OF THERMAL WARPING AND BUCKLING IN ENAMELLED STEEL PARTS

NUMERICAL SIMULATION OF THERMAL WARPING AND BUCKLING IN ENAMELLED STEEL PARTS NUMERICAL SIMULATION OF THERMAL WARPING AND BUCKLING IN ENAMELLED STEEL PARTS 337 XXI Intrnational Enamllrs Congrss Numrical Simulation of Thrmal Warping and Buckling in Enamlld Stl Parts Filip Van dn

More information

Twist analysis of piezoelectric laminated composite plates

Twist analysis of piezoelectric laminated composite plates wist analysis of pizolctric laminatd composit plats Mchatronics Enginring Dpartmnt, Faculty of Enginring, Intrnational Islamic Univrsity Malaysia, Malaysia raisuddin@iiu.du.my ABSAC cntly scintists ar

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

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

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

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

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

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

arxiv: v1 [physics.comp-ph] 30 Jun 2016

arxiv: v1 [physics.comp-ph] 30 Jun 2016 On anisotropy function in crystal growth simulations using Lattic Boltzmann quation AMINA YOUNSI a,1, ALAIN CARTALADE a, a Dn DM2S, STMF, LMSF, CEA, Univrsité d Paris-Saclay, F-91191, Gif-sur-Yvtt, Franc.

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

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

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

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

Optimal Design Schemes of Centrifugal Pump Compound Impeller at Super-Low Specific Speed

Optimal Design Schemes of Centrifugal Pump Compound Impeller at Super-Low Specific Speed 811 A publication of CHEMICAL ENGINEERING TRANSACTIONS VOL., 15 Gust Editors: Piyu Rn, Yancang Li, Huiping Song Copyright 15, AIDIC Srvizi S.r.l., ISBN 78-88-58-7-; ISSN 8-1 Th Italian Association of Chmical

More information

Review Statistics review 14: Logistic regression Viv Bewick 1, Liz Cheek 1 and Jonathan Ball 2

Review Statistics review 14: Logistic regression Viv Bewick 1, Liz Cheek 1 and Jonathan Ball 2 Critical Car Fbruary 2005 Vol 9 No 1 Bwick t al. Rviw Statistics rviw 14: Logistic rgrssion Viv Bwick 1, Liz Chk 1 and Jonathan Ball 2 1 Snior Lcturr, School of Computing, Mathmatical and Information Scincs,

More information

Full Waveform Inversion Using an Energy-Based Objective Function with Efficient Calculation of the Gradient

Full Waveform Inversion Using an Energy-Based Objective Function with Efficient Calculation of the Gradient Full Wavform Invrsion Using an Enrgy-Basd Objctiv Function with Efficint Calculation of th Gradint Itm yp Confrnc Papr Authors Choi, Yun Sok; Alkhalifah, ariq Ali Citation Choi Y, Alkhalifah (217) Full

More information

Dynamic behaviour of a rotating cracked beam

Dynamic behaviour of a rotating cracked beam Journal of Physics: Confrnc Sris PAPER OPEN ACCESS Dynamic bhaviour of a rotating crackd bam To cit this articl: Ahmd Yashar t al 6 J. Phys.: Conf. Sr. 744 57 Viw th articl onlin for updats and nhancmnts.

More information

FINITE BEAM ELEMENT WITH PIEZOELECTRIC LAYERS AND FUNCTIONALLY GRADED MATERIAL OF CORE

FINITE BEAM ELEMENT WITH PIEZOELECTRIC LAYERS AND FUNCTIONALLY GRADED MATERIAL OF CORE ECCOMAS Congrss 20 II Europan Congrss on Computational Mthods in Applid Scincs and Enginring M. Papadrakakis,. Papadopoulos, G. Stfanou,. Plvris (ds.) Crt Island, Grc, 5 0 Jun 20 FINITE BEAM ELEMENT WITH

More information

Numerical Analysis of Transient Responses for Elastic Structures Connected to a Viscoelastic Shock Absorber Using FEM with a Nonlinear Complex Spring

Numerical Analysis of Transient Responses for Elastic Structures Connected to a Viscoelastic Shock Absorber Using FEM with a Nonlinear Complex Spring Numrical Analysis of Transint Rsponss for Elastic Structurs Connctd to a Viscolastic Shock Absorbr Using FEM with a Nonlinar Complx Spring Takao Yamaguchi, Yusaku Fujii, Toru Fukushima, Akihiro Takita,

More information

Observer Bias and Reliability By Xunchi Pu

Observer Bias and Reliability By Xunchi Pu Obsrvr Bias and Rliability By Xunchi Pu Introduction Clarly all masurmnts or obsrvations nd to b mad as accuratly as possibl and invstigators nd to pay carful attntion to chcking th rliability of thir

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

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

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

Effects of Wave Non-Linearity on Residual Pore Pressures in Marine Sediments

Effects of Wave Non-Linearity on Residual Pore Pressures in Marine Sediments Th Opn Civil Enginring Journal, 8,, 63-74 63 Effcts of Wav Non-Linarity on Rsidual Por Prssurs in Marin Sdimnts Dong-Shng Jng* Opn Accss Division of Civil Enginring, School of Enginring, Physics and Mathmatics,

More information

High Energy Physics. Lecture 5 The Passage of Particles through Matter

High Energy Physics. Lecture 5 The Passage of Particles through Matter High Enrgy Physics Lctur 5 Th Passag of Particls through Mattr 1 Introduction In prvious lcturs w hav sn xampls of tracks lft by chargd particls in passing through mattr. Such tracks provid som of th most

More information

FE modeling of inelastic behavior of reinforced high-strength concrete continuous beams

FE modeling of inelastic behavior of reinforced high-strength concrete continuous beams Structural Enginring and Mchanics, Vol. 49, No. 3 (214) 373-393 DOI: http://dx.doi.org/1.12989/sm.214.49.3.373 373 FE modling of inlastic bhavior of rinforcd high-strngth concrt continuous bams Tijiong

More information

MCE503: Modeling and Simulation of Mechatronic Systems Discussion on Bond Graph Sign Conventions for Electrical Systems

MCE503: Modeling and Simulation of Mechatronic Systems Discussion on Bond Graph Sign Conventions for Electrical Systems MCE503: Modling and Simulation o Mchatronic Systms Discussion on Bond Graph Sign Convntions or Elctrical Systms Hanz ichtr, PhD Clvland Stat Univrsity, Dpt o Mchanical Enginring 1 Basic Assumption In a

More information

Performance Prediction of the Single-Sided Linear. Induction Motors for Transportation Considers. Longitudinal End Effect by Using Analytic Method

Performance Prediction of the Single-Sided Linear. Induction Motors for Transportation Considers. Longitudinal End Effect by Using Analytic Method Contmporary Enginring Scincs, Vol., 9, no., 95-14 Prformanc Prdiction of th Singl-Sidd Linar Induction Motors for Transportation Considrs Longitudinal End Effct by Using Analytic Mthod Ali Suat Grçk Univrsity

More information

MCB137: Physical Biology of the Cell Spring 2017 Homework 6: Ligand binding and the MWC model of allostery (Due 3/23/17)

MCB137: Physical Biology of the Cell Spring 2017 Homework 6: Ligand binding and the MWC model of allostery (Due 3/23/17) MCB37: Physical Biology of th Cll Spring 207 Homwork 6: Ligand binding and th MWC modl of allostry (Du 3/23/7) Hrnan G. Garcia March 2, 207 Simpl rprssion In class, w drivd a mathmatical modl of how simpl

More information

SME 3033 FINITE ELEMENT METHOD. Bending of Prismatic Beams (Initial notes designed by Dr. Nazri Kamsah)

SME 3033 FINITE ELEMENT METHOD. Bending of Prismatic Beams (Initial notes designed by Dr. Nazri Kamsah) Bnding of Prismatic Bams (Initia nots dsignd by Dr. Nazri Kamsah) St I-bams usd in a roof construction. 5- Gnra Loading Conditions For our anaysis, w wi considr thr typs of oading, as iustratd bow. Not:

More information

Forces. Quantum ElectroDynamics. α = = We have now:

Forces. Quantum ElectroDynamics. α = = We have now: W hav now: Forcs Considrd th gnral proprtis of forcs mdiatd by xchang (Yukawa potntial); Examind consrvation laws which ar obyd by (som) forcs. W will nxt look at thr forcs in mor dtail: Elctromagntic

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

VII. Quantum Entanglement

VII. Quantum Entanglement VII. Quantum Entanglmnt Quantum ntanglmnt is a uniqu stat of quantum suprposition. It has bn studid mainly from a scintific intrst as an vidnc of quantum mchanics. Rcntly, it is also bing studid as a basic

More information

A nonequilibrium molecular dynamics simulation of evaporation

A nonequilibrium molecular dynamics simulation of evaporation Intrnational Confrnc Passiv and Low Enrgy Cooling 543 A nonquilibrium molcular dynamics simulation of vaporation Z.-J. Wang, M. Chn and Z.-Y. Guo Dpartmnt of Enginring Mchanics, Tsinghua Univrsity, Bijing

More information

Sliding Mode Flow Rate Observer Design

Sliding Mode Flow Rate Observer Design Sliding Mod Flow Rat Obsrvr Dsign Song Liu and Bin Yao School of Mchanical Enginring, Purdu Univrsity, Wst Lafaytt, IN797, USA liu(byao)@purdudu Abstract Dynamic flow rat information is ndd in a lot of

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

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

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

Determination of Vibrational and Electronic Parameters From an Electronic Spectrum of I 2 and a Birge-Sponer Plot

Determination of Vibrational and Electronic Parameters From an Electronic Spectrum of I 2 and a Birge-Sponer Plot 5 J. Phys. Chm G Dtrmination of Vibrational and Elctronic Paramtrs From an Elctronic Spctrum of I 2 and a Birg-Sponr Plot 1 15 2 25 3 35 4 45 Dpartmnt of Chmistry, Gustavus Adolphus Collg. 8 Wst Collg

More information

4037 ADDITIONAL MATHEMATICS

4037 ADDITIONAL MATHEMATICS CAMBRIDGE INTERNATIONAL EXAMINATIONS GCE Ordinary Lvl MARK SCHEME for th Octobr/Novmbr 0 sris 40 ADDITIONAL MATHEMATICS 40/ Papr, maimum raw mark 80 This mark schm is publishd as an aid to tachrs and candidats,

More information

A System Identification Algorithm for Vehicle Lumped Parameter Model in Crash Analysis

A System Identification Algorithm for Vehicle Lumped Parameter Model in Crash Analysis Intrnational Journal of Modling and Optimization, Vol., No., Jun A Systm Idntification Algorithm for Vhicl Lumpd Paramtr Modl in Crash Analysis Javad Marzbanrad and Mostafa Pahlavani Abstract This study

More information

EFFECT OF BALL PROPERTIES ON THE BALL-BAT COEFFICIENT OF RESTITUTION

EFFECT OF BALL PROPERTIES ON THE BALL-BAT COEFFICIENT OF RESTITUTION EFFECT OF BALL PROPERTIES ON THE BALL-BAT COEFFICIENT OF RESTITUTION A. M. NATHAN 1 AND L. V. SMITH 2 1 Univrsity of Illinois, 1110 W. Grn Strt, Urbana, IL 61801, USA, E-mail: a-nathan@illinois.du 2 Washington

More information

Introduction to Multicopter Design and Control

Introduction to Multicopter Design and Control Introduction to Multicoptr Dsign and Control Lsson 05 Coordinat Systm and Attitud Rprsntation Quan Quan, Associat Profssor _uaa@uaa.du.cn BUAA Rlial Flight Control Group, http://rfly.uaa.du.cn/ Bihang

More information

Recursive Estimation of Dynamic Time-Varying Demand Models

Recursive Estimation of Dynamic Time-Varying Demand Models Intrnational Confrnc on Computr Systms and chnologis - CompSysch 06 Rcursiv Estimation of Dynamic im-varying Dmand Modls Alxandr Efrmov Abstract: h papr prsnts an implmntation of a st of rcursiv algorithms

More information

Dynamic analysis of a Timoshenko beam subjected to moving concentrated forces using the finite element method

Dynamic analysis of a Timoshenko beam subjected to moving concentrated forces using the finite element method Shock and Vibration 4 27) 459 468 459 IOS Prss Dynamic analysis of a Timoshnko bam subjctd to moving concntratd forcs using th finit lmnt mthod Ping Lou, Gong-lian Dai and Qing-yuan Zng School of Civil

More information

KINEMATIC SOIL-STRUCTURE INTERACTION EFFECTS ON MAXIMUM INELASTIC DISPLACEMENT DEMANDS OF SDOF SYSTEMS

KINEMATIC SOIL-STRUCTURE INTERACTION EFFECTS ON MAXIMUM INELASTIC DISPLACEMENT DEMANDS OF SDOF SYSTEMS Th 14 th World Confrnc on Earthquak Enginring Octobr 12-17, 2008, Bijing, China KINEMATIC SOIL-STRUCTURE INTERACTION EFFECTS ON MAXIMUM INELASTIC DISPLACEMENT DEMANDS OF SDOF SYSTEMS Y.Y. Lin 1 1 Associat

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