A NEW FREQUENCY METHOD APPLIED TO CALCULATE CABLE FORCE FOR CABLES WITH THE INTERMEDIATE MULTI-SUPPORT

Size: px
Start display at page:

Download "A NEW FREQUENCY METHOD APPLIED TO CALCULATE CABLE FORCE FOR CABLES WITH THE INTERMEDIATE MULTI-SUPPORT"

Transcription

1 Blucher Mechancal Engneerng Proceedngs Ma 014, vol. 1, num. 1 A NEW FREQUENCY METHOD APPLIED TO CALCULATE CABLE FORCE FOR CABLES WITH THE INTERMEDIATE MULTI-SUPPORT Pngje L Ronghu Wang Nujng Ma School of Cvl Engneerng and Transportaton, South Chna Unverst of Technolog(pngjel@163.com) Abstract. Ths paper proposes a new frequenc method appled to calculate cable force for cables wth the ntermedate mult-support from measured natural frequences. The proposed method has establshed new relatonshp between cable frequenc and tenson force, takng cable force, the ntermedate mult-support, flexural rgdt, boundar condtons of two ends nto account. Followng the classcal gudelnes of cable dnamcs, and b makng some smplfcatons, frequenc equaton has been obtaned n non-dmensonal form. When natural frequences have measured, cable force can be calculated b solvng frequenc equaton. Comparson wth Matharmethod and the proposed method n the feld test, the dfference of cable force result between two methods s just onl 1%, so the proposed method has been made to valdate. Kewords:frequenc method, cable force, mult-support, cable 1. INTRODUCTION Frequenc method s the wa to calculate the cable force b measured natural frequences of cables. In recent ears, Frequenc method has receved ncreasng attenton because of ts smplct and speedness. In practce, frequenc method s wdel used for cables wth boundar condtons of constrants on both ends, and rarel used for cables wth boundar condtons of constrants on both ends and the ntermedate. As cables are appled more and more wdel, the boundar condton of cables s more and more complex, such as the damper at suspenders n the suspenson brdge and ntermedate support at the te bar n the arch brdge. So a new frequenc method s useful for applcaton to calculate the tenson force of cables wth the ntermedate mult-support. Nowadas, there are a lot of researches on the vbraton of cables. Irvne H M [1] has proposed an analtcal soluton for the frequenc equaton for sagged cable wthout bendng stffness. For such long cables, wth such large dameters, flexural stffness cannot be neglected n determnng dnamcal propertes of cables. So Zu H [],Geer Rand Km B H [3-4] has proposed methods whch allow us to obtan cable force n a smple manner, whch are based on the dentfcaton of modal propertes of cables. Flamand O and Matsumoto M [5-6] have observed the nfluence of ran wnd acton to cables. These studes are mostl for cables wth the boundar condton of constrants on both ends, and there are rare for cables wth the boundar condton of 1

2 constrants on both ends and the ntermedate. R H wang [7] has put forward to a fnte element soluton for tenson force of cables wth elastc mult-support, whch has been expermentall verfed. However, there s not an analtc method for tenson force of cables wth elastc mult-support. Ths paper proposes a new frequenc method for applcaton to calculate the tenson force of cables wth the ntermedate mult-support from measured natural frequences. The proposed method has establshed new relatonshp between cable frequenc and tenson force, takng cable force, the ntermedate mult-support, flexural rgdt, boundar condtons of two ends nto account. Followng the classcal gudelnes of cable dnamcs, and b makng some smplfcatons, frequenc equaton has been obtaned n non-dmensonal form. When natural frequences have measured, cable force can be calculated b solvng frequenc equaton..dynamic MODEL OF CABLES.1Vbraton model for cables The cable wth the boundar condton of restrants on both ends and the ntermedate s equvalent for the mult-span beam wth axal force, when flexural stffness of cables has not been neglected., the calculated parameters can be descrbed as show n fgure 1: L s the undeformed cable length; H s the horzontal component of ntal cable tenson; n s the number of constrants on the ntermedat of cables; K (=1,,...,n) s stffness of constrants; X (=1,,...,n) s the dstance between the adjacent constrants. The boundar condton of constrants on both ends s consder as clamped. Fgure 1. the cable wth the ntermedate mult-support In accordance wth [1], the lnearzed equatons of moton are: 4 v d v H h() t EI (,) 4 v f x t mv (1) x dx x v s the transverse vbraton dsplacement of cables changng over tme, ht () s the addtonal horzontal component of tenson changng over tme, ( xt, ) s the transverse coordnate of cables, x s the coordnate along the axs of beam, t s

3 vbraton tme, E s the Young modulus; I s the nerta moment; s vscous dampng coeffcents per untar cable length, v s frst tme dervatve of the vbraton dsplacement, v s second tme dervatve of the vbraton dsplacement, m s the cable mass per untar length, f ( xt, ) s mpressed force changng over tme. The vbraton model of the cable has consdered elastc support onl, but hasn t consdered the effect of dampng, so 0.accordng to Irvne H M [1], It can be neglgble to the addtonal horzontal component of tenson [1], so ht () =0. Assumng v v x, t x qt ( ), f ( xt, )= f ( xqt ) ( ),the partal dfferental equaton (1) s transformed nto two ordnar dfferental equatons, one governng the tme evoluton and one governng the spatal dstrbuton of the moton: qt () qt () 0 (a) IV H m 1 ( x) ( x) ( x) f ( x) (b) EI EI EI Expresson (b) s the fourth order nonhomogeneous lnear dfferental equaton, whch can be solved as: 1 ( x) A1cos( x) Asn( x) A3cosh( x) A4snh( x) f ( ) x (3) m In whch Aj ( j = 1,..., 4) are constants and the followng quanttes have been ntroduced: ( ) ( H ) m H EI EI EI (4a) ( ) ( H ) m H EI EI EI (4b).The cable frequenc equaton Boundar condtons for clamped restrant at cable anchorages are: v( L/, t) 0 (5a) vl ( /, t) 0 (5b) v ( L/, t) 0 x v ( L /, t ) 0 (5d) x The boundar condtons of cables wth clamped ends, gven b Eqs. (5a)- (5d), 3 (5c)

4 substtute nto Eq (3); And then, an algebrac set of four equatons, wth the unknown quanttes Bj ( j = 1,..., 4), can be solved; And Bj ( j = 1,..., 4) can be substtuted nto Eq. (3) and the fnal expresson of the modal shapes s obtaned followng: fˆ ( ) ˆcosh( ˆ ˆ)sn( ˆ / ) ˆcos( ˆ ˆ)snh( ˆ x x x / ) ˆ ( xˆ ) 1 ˆ ˆ cosh( ˆ / )sn( ˆ / ) ˆ cos( ˆ / )snh( ˆ / ) (6) The parameters of expresson (6) can be obtaned b Eqs. (7a)- (7f). ( x x ) s the Drac delta functon. x xˆ L (7a) L ˆ (7b) EI / m L H / EI (7c) ˆ ( ˆ) ( ) ˆ (7d) ˆ ( ˆ) ( ) ˆ (7e) ˆ KL f ( x) ( x ) ( xx ) ( 1,,..., n) (7f) 4 n ˆ ˆ ˆ ˆ ˆ EI 1 So frequenc equaton of cables can be gotten b Eq (6) and Eq (7f), the smplfed expressons s wrtten as: C I 0 (8) The parameters of expresson (8) can be gotten b Eqs. (9a)- (9b) EI ˆ ( 1,,..., n) (9a) 4 KL ˆ cosh( ˆ xˆ )sn( ˆ / ) ˆ cos( ˆ ˆ )snh( ˆ x / ) Cj 1 ( j 1,,..., n 1,,..., n) ˆ cosh( ˆ / )sn( ˆ / ) ˆ cos( ˆ / )snh( ˆ / ) (9b) 3 CABLE FORCE CALCULATION 3.1 Cable force equaton Cable frequenc equaton (8) contans two non-dmensonal parameters, one s 4

5 the non-dmensonal frequenc ˆ, and the other s the non-dmensonal cable force. So cable force equaton can be obtaned b eqs (8), (7d), (7e): f ( ˆ, ) 0 As cable frequenc s tested out, the non-dmensonal frequenc ˆ can be gotten. So expresson (10) s a sngle parameter equaton as show n expresson (11). f ( ) 0 (10) (11) 3. Iteratve algorthm for the cable force equaton As we known, the Eq(11) s a transcendental equaton. It can t be solve drectl. So t must adopt numercal methods, lke Nowton-Raphson, for the purpose of obtanng the roots of Eq. (11).It s convenent to manpulate ths expresson to elmnate some dfferences, whch, when the are made between close approxmate numbers, can lead to round off errors as show n fgure. f mleik,,,, xˆ, ˆ, ˆ( ˆ), ˆ ( ˆ), fˆ ( xˆ),, C j H 1 ˆn f n H n 1 H n H H n 1 H n H f f f,suppose error as = 10 n 4 fn f 0 f f 0 n H x Fgure. flowsheet for cable force calculaton 4 APPLICATION EXAMPLE 4.1 background 5

6 Sanshanx Brdge s concrete flled steel tubular ted arch brdge. Te bar, consstng of hundreds of steel strands, s located on the both end of brdge under the brdge deck wth the length of 60m,whch s used for resstng horzontal thrust of man arch. A bunch of te ncludes 10 to 1 steel strands, and these steel strands are contact each other. There are 50 bunches of te separated. Steel strands have been separated to three parts b ponts slk plate. The length of parts of steel strands s 49.7m, 160.6m, 49.7m respectvel. Ponts slk plate s weldng nsde steel box, whch s usng for protectng the te bar from corroson. Steel strand and ponts slk plate contact tghtl, so each bunch of te bar can be consdered to be anchored at the place of ponts slk plate. It can t be completed to determne the axal force of te bar b applng frequenc method drectl f we have not solved the followng two problems. One s how to test frequenc of te bar. Steel strands are contact each other, so we can appl some fllng blocks to separate the steel strand, what s mportant s determne the compresson stffness of fllng blocks. The other s how to calculate the axal force of te bar b measured frequenc. There s changed boundar condton of steel strand b ncreased fllng blocks. So t s essental to deduce a new equaton of cable force and frequenc for cables wth the ntermedate mult-support. 4. Test parameter Te bar and steel strand both can be consder to be the cable, but t just can be test one steel strand frequenc, so the cable s stand for steel strand n the followng. The cable parameter s descrbng at table 1, E s the modulus of elastct, I s the nerta moment, s the mass denst, A s area, d s dameter of secton, L s the length of cables. E(pa) I(m 4 ) Table 1. test parameter (kg/m 3 A(m ) d(mm) L(m) ) 1.95E E E The fllng block has three laers. There s steel plate at the top and bottom laer wth the thckness of 5mm and block rubber at the mddle laer wth the thckness of 40mm. The sze of fllng blocks s 5cm 5cm.The compresson stffness of fllng blocks s determned b traxal compresson test. Two fllng blocks have been made 6 and the compresson stffness of them are respectl : K N/m 1 K N/m. 4.3 Cable frequenc and cable force Cable frequenc cannot be tested drectl because of the contact between cables, so fllng blocks are used for solatng cables as show n fgure 3. And cable frequenc 6

7 can be tested b dnamc testng nstrument when optmum vbraton sstem set at the mddle of two fllng blocks as show n fgure 4, whch are nstalled under the cable as show n fgure 5. Fgure 3. The model for testng cable frequenc Fgure 4. The nstallng place of optmum Fgure 5. The nstallng place of the fllng vbraton sstem block It s fnd that the 4th order frequenc of cables can be eas tested as show n table. Table. Tested cable frequenc number The 4th order cable frequenc(hz) B flowsheet for cable force calculaton (fg.), cable force has calculated as show n table 3. The average cable force s 97.4kN, whch can evaluate the whole level of te bar force. Usng non-destructve testng (NDT) methods of Matharmethod [8], te bar force can be tested to be 98.3kN. So the dfference of te bar force tested b two methods s small, just wth the error of 1%. Table 3. Calculated cable force number average cable force(kn)

8 5 CONCLUSION The computatonal formula of frequenc method cannot calculate cable force for cables wth boundar condtons of constrants on the ntermedate. so a new computatonal formula of frequenc method s establshed for cables wth the ntermedate mult-support. Some conclusons are made as followng: 1)The proposal computatonal formula can be useful for applcaton to calculate the tenson force of cables wth the ntermedate mult-support, such as the damper at suspenders n the suspenson brdge and ntermedate support beam at the te bar n the arch brdge. )The proposal computatonal formula s a transcendental equaton, so numercal method should be adopted, but t s eas to make a programme, followng fg,to calculatng cable force. 3)The feld test shows that the dfference of cable force, testng b the proposal method and the matharmethod,s small, so the proposal method s relable. Acknowledgement Ths work s supported b Program of Natonal Natural Scence Foundaton of Chna ( ) REFERENCE [1]Irvne H M. Cable Structures [M]. Cambrge:The MIT Press, [] Zu H, Shnke T, Namta Y. Practcal formulas for estmaton of cable tenson b vbraton method [J]. J Struct Eng 1996, 1(6): [3] Geer R, De Roeck G, Flesch R. Accurate cable force determnaton usng ambent vbraton measurements. Struct Infrastruct Eng 006;(1):43 5. [4] Km B H, Park T. Estmaton of cable tenson force usng the frequenc-based sstem dentfcaton method [J]. Journal of Sound and Vbraton. 007, 304(3-5): [5] Flamand O. Ran-wnd nduced vbraton of cables. J Wnd Eng Ind Aerodn 1995;57: [6] Matsumoto M, Dato Y, Kanamura T, Shgemura Y, Sakuma S,Ishzak H. Wnd-nduced vbraton of cables of cable-staed brdges.j of Wnd Eng and Ind Aerodn 1998;74 76: [7]Wang Rong-hu,Gan Quan,Gong Lng-lng,et al.a fnte element soluton for tenson of cable wth elastc mult-support[j]. Journal of Vbraton Engneerng, 010, 3(5): [8]BA50/93. Post-tensoned Concrete Brdges. Plannng, Organsaton and Methods for Carrng out Specal Inspectons. Department of Transport. Jul

Indeterminate pin-jointed frames (trusses)

Indeterminate pin-jointed frames (trusses) Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all

More information

829. An adaptive method for inertia force identification in cantilever under moving mass

829. An adaptive method for inertia force identification in cantilever under moving mass 89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(

More information

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

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

More information

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS

More information

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

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

More information

DUE: WEDS FEB 21ST 2018

DUE: WEDS FEB 21ST 2018 HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant

More information

UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 2017/2018 FINITE ELEMENT AND DIFFERENCE SOLUTIONS

UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 2017/2018 FINITE ELEMENT AND DIFFERENCE SOLUTIONS OCD0 UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 07/08 FINITE ELEMENT AND DIFFERENCE SOLUTIONS MODULE NO. AME6006 Date: Wednesda 0 Ma 08 Tme: 0:00

More information

PART 8. Partial Differential Equations PDEs

PART 8. Partial Differential Equations PDEs he Islamc Unverst of Gaza Facult of Engneerng Cvl Engneerng Department Numercal Analss ECIV 3306 PAR 8 Partal Dfferental Equatons PDEs Chapter 9; Fnte Dfference: Ellptc Equatons Assocate Prof. Mazen Abualtaef

More information

Grid Generation around a Cylinder by Complex Potential Functions

Grid Generation around a Cylinder by Complex Potential Functions Research Journal of Appled Scences, Engneerng and Technolog 4(): 53-535, 0 ISSN: 040-7467 Mawell Scentfc Organzaton, 0 Submtted: December 0, 0 Accepted: Januar, 0 Publshed: June 0, 0 Grd Generaton around

More information

New Method for Solving Poisson Equation. on Irregular Domains

New Method for Solving Poisson Equation. on Irregular Domains Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad

More information

THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES

THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES W. C. Lao Department of Cvl Engneerng, Feng Cha Unverst 00 Wen Hwa Rd, Tachung, Tawan SUMMARY: The ndentaton etween clndrcal ndentor

More information

Numerical Solution of Ordinary Differential Equations

Numerical Solution of Ordinary Differential Equations Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples

More information

Module 3: Element Properties Lecture 1: Natural Coordinates

Module 3: Element Properties Lecture 1: Natural Coordinates Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers

More information

GEO-SLOPE International Ltd, Calgary, Alberta, Canada Vibrating Beam

GEO-SLOPE International Ltd, Calgary, Alberta, Canada   Vibrating Beam GEO-SLOPE Internatonal Ltd, Calgary, Alberta, Canada www.geo-slope.com Introducton Vbratng Beam Ths example looks at the dynamc response of a cantlever beam n response to a cyclc force at the free end.

More information

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,

More information

Numerical Heat and Mass Transfer

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

More information

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion Assessment of Ste Amplfcaton Effect from Input Energy Spectra of Strong Ground Moton M.S. Gong & L.L Xe Key Laboratory of Earthquake Engneerng and Engneerng Vbraton,Insttute of Engneerng Mechancs, CEA,

More information

An identification algorithm of model kinetic parameters of the interfacial layer growth in fiber composites

An identification algorithm of model kinetic parameters of the interfacial layer growth in fiber composites IOP Conference Seres: Materals Scence and Engneerng PAPER OPE ACCESS An dentfcaton algorthm of model knetc parameters of the nterfacal layer growth n fber compostes o cte ths artcle: V Zubov et al 216

More information

Finite Element Modelling of truss/cable structures

Finite Element Modelling of truss/cable structures Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

A new Approach for Solving Linear Ordinary Differential Equations

A new Approach for Solving Linear Ordinary Differential Equations , ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of

More information

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc

More information

Chapter Newton s Method

Chapter Newton s Method Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve

More information

Numerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method

Numerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method Appled Mathematcs, 6, 7, 5-4 Publshed Onlne Jul 6 n ScRes. http://www.scrp.org/journal/am http://.do.org/.436/am.6.77 umercal Solutons of a Generalzed th Order Boundar Value Problems Usng Power Seres Approxmaton

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

More information

MODELLING OF ELASTO-STATICS OF POWER LINES BY NEW COMPOSITE BEAM FINITE ELEMENT Bratislava

MODELLING OF ELASTO-STATICS OF POWER LINES BY NEW COMPOSITE BEAM FINITE ELEMENT Bratislava ODING OF ASTO-STATICS OF POW INS BY NW COPOSIT BA FINIT NT urín Justín 1 rabovský Jura 1 Gogola oman 1 utš Vladmír 1 Paulech Jura 1 1 Insttute of Automotve echatroncs FI STU n Bratslava Ilkovčova 3 812

More information

Δ x. u(x,t) Fig. Schematic view of elastic bar undergoing axial motions

Δ x. u(x,t) Fig. Schematic view of elastic bar undergoing axial motions ME67 - Handout 4 Vbratons of Contnuous Systems Axal vbratons of elastc bars The fgure shows a unform elastc bar of length and cross secton A. The bar materal propertes are ts densty ρ and elastc modulus

More information

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1

More information

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project SE 8 Fnal Project Story Shear Frame u m Gven: u m L L m L L EI ω ω Solve for m Story Bendng Beam u u m L m L Gven: m L L EI ω ω Solve for m 3 3 Story Shear Frame u 3 m 3 Gven: L 3 m m L L L 3 EI ω ω ω

More information

A Mechanics-Based Approach for Determining Deflections of Stacked Multi-Storey Wood-Based Shear Walls

A Mechanics-Based Approach for Determining Deflections of Stacked Multi-Storey Wood-Based Shear Walls A Mechancs-Based Approach for Determnng Deflectons of Stacked Mult-Storey Wood-Based Shear Walls FPINNOVATIONS Acknowledgements Ths publcaton was developed by FPInnovatons and the Canadan Wood Councl based

More information

Application to Plane (rigid) frame structure

Application to Plane (rigid) frame structure Advanced Computatonal echancs 18 Chapter 4 Applcaton to Plane rgd frame structure 1. Dscusson on degrees of freedom In case of truss structures, t was enough that the element force equaton provdes onl

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

Computers and Mathematics with Applications

Computers and Mathematics with Applications Computers and Mathematcs wth Applcatons 56 (2008 3204 3220 Contents lsts avalable at ScenceDrect Computers and Mathematcs wth Applcatons journal homepage: www.elsever.com/locate/camwa An nnovatve egenvalue

More information

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

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

More information

Δ x. u(x,t) Fig. Schematic view of elastic bar undergoing axial motions

Δ x. u(x,t) Fig. Schematic view of elastic bar undergoing axial motions ME67 - Handout 4 Vbratons of Contnuous Systems Axal vbratons of elastc bars The fgure shows a unform elastc bar of length and cross secton A. The bar materal propertes are ts densty ρ and elastc modulus

More information

Week 9 Chapter 10 Section 1-5

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

More information

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI

More information

Free torsional vibration analysis of H-section hangers attached with horizontal wind-resistant cables

Free torsional vibration analysis of H-section hangers attached with horizontal wind-resistant cables Free torsonal vbraton analyss of H-secton hangers attached wth horzontal wnd-resstant cables * Ka Xu ), Yang Zhao ), Zhhao Wang 3) ), ), 3) Department of Cvl Engneerng, North Chna Unversty of Water Resources

More information

APPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam

APPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam APPENDIX F A DISPACEMENT-BASED BEAM EEMENT WITH SHEAR DEFORMATIONS Never use a Cubc Functon Approxmaton for a Non-Prsmatc Beam F. INTRODUCTION { XE "Shearng Deformatons" }In ths appendx a unque development

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

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

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

More information

Power law and dimension of the maximum value for belief distribution with the max Deng entropy

Power law and dimension of the maximum value for belief distribution with the max Deng entropy Power law and dmenson of the maxmum value for belef dstrbuton wth the max Deng entropy Bngy Kang a, a College of Informaton Engneerng, Northwest A&F Unversty, Yanglng, Shaanx, 712100, Chna. Abstract Deng

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur Module Random Processes Lesson 6 Functons of Random Varables After readng ths lesson, ou wll learn about cdf of functon of a random varable. Formula for determnng the pdf of a random varable. Let, X be

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

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

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

More information

Torsion Stiffness of Thin-walled Steel Beams with Web Holes

Torsion Stiffness of Thin-walled Steel Beams with Web Holes Torson Stffness of Thn-walled Steel Beams wth Web Holes MARTN HORÁČEK, JNDŘCH MELCHER Department of Metal and Tmber Structures Brno Unversty of Technology, Faculty of Cvl Engneerng Veveří 331/95, 62 Brno

More information

Newton s Method for One - Dimensional Optimization - Theory

Newton s Method for One - Dimensional Optimization - Theory Numercal Methods Newton s Method for One - Dmensonal Optmzaton - Theory For more detals on ths topc Go to Clck on Keyword Clck on Newton s Method for One- Dmensonal Optmzaton You are free to Share to copy,

More information

Lecture 21: Numerical methods for pricing American type derivatives

Lecture 21: Numerical methods for pricing American type derivatives Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)

More information

Kinematics in 2-Dimensions. Projectile Motion

Kinematics in 2-Dimensions. Projectile Motion Knematcs n -Dmensons Projectle Moton A medeval trebuchet b Kolderer, c1507 http://members.net.net.au/~rmne/ht/ht0.html#5 Readng Assgnment: Chapter 4, Sectons -6 Introducton: In medeval das, people had

More information

ANALYSIS OF TIMOSHENKO BEAM RESTING ON NONLINEAR COMPRESSIONAL AND FRICTIONAL WINKLER FOUNDATION

ANALYSIS OF TIMOSHENKO BEAM RESTING ON NONLINEAR COMPRESSIONAL AND FRICTIONAL WINKLER FOUNDATION VOL. 6, NO., NOVEMBER ISSN 89-668 6- Asan Research Publshng Network (ARPN). All rghts reserved. ANALYSIS OF TIMOSHENKO BEAM RESTING ON NONLINEAR COMPRESSIONAL AND FRICTIONAL WINKLER FOUNDATION Adel A.

More information

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

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

More information

Image Processing for Bubble Detection in Microfluidics

Image Processing for Bubble Detection in Microfluidics Image Processng for Bubble Detecton n Mcrofludcs Introducton Chen Fang Mechancal Engneerng Department Stanford Unverst Startng from recentl ears, mcrofludcs devces have been wdel used to buld the bomedcal

More information

2 Finite difference basics

2 Finite difference basics Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T

More information

Lecture 8 Modal Analysis

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

More information

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

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

More information

6 Supplementary Materials

6 Supplementary Materials 6 Supplementar Materals 61 Proof of Theorem 31 Proof Let m Xt z 1:T : l m Xt X,z 1:t Wethenhave mxt z1:t ˆm HX Xt z 1:T mxt z1:t m HX Xt z 1:T + mxt z 1:T HX We consder each of the two terms n equaton

More information

Formulas for the Determinant

Formulas for the Determinant page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use

More information

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam. Please revew the followng statement: I certfy that I have not gven unauthorzed ad nor have I receved ad n the completon of ths exam. Sgnature: Instructor s Name and Secton: (Crcle Your Secton) Sectons:

More information

SCALARS AND VECTORS All physical quantities in engineering mechanics are measured using either scalars or vectors.

SCALARS AND VECTORS All physical quantities in engineering mechanics are measured using either scalars or vectors. SCALARS AND ECTORS All phscal uanttes n engneerng mechancs are measured usng ether scalars or vectors. Scalar. A scalar s an postve or negatve phscal uantt that can be completel specfed b ts magntude.

More information

2. PROBLEM STATEMENT AND SOLUTION STRATEGIES. L q. Suppose that we have a structure with known geometry (b, h, and L) and material properties (EA).

2. PROBLEM STATEMENT AND SOLUTION STRATEGIES. L q. Suppose that we have a structure with known geometry (b, h, and L) and material properties (EA). . PROBEM STATEMENT AND SOUTION STRATEGIES Problem statement P, Q h ρ ρ o EA, N b b Suppose that we have a structure wth known geometry (b, h, and ) and materal propertes (EA). Gven load (P), determne the

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

SIMPLE LINEAR REGRESSION

SIMPLE LINEAR REGRESSION Smple Lnear Regresson and Correlaton Introducton Prevousl, our attenton has been focused on one varable whch we desgnated b x. Frequentl, t s desrable to learn somethng about the relatonshp between two

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

11. Dynamics in Rotating Frames of Reference

11. Dynamics in Rotating Frames of Reference Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons

More information

CHAPTER 9 CONCLUSIONS

CHAPTER 9 CONCLUSIONS 78 CHAPTER 9 CONCLUSIONS uctlty and structural ntegrty are essentally requred for structures subjected to suddenly appled dynamc loads such as shock loads. Renforced Concrete (RC), the most wdely used

More information

NONLINEAR NATURAL FREQUENCIES OF A TAPERED CANTILEVER BEAM

NONLINEAR NATURAL FREQUENCIES OF A TAPERED CANTILEVER BEAM Advanced Steel Constructon Vol. 5, No., pp. 59-7 (9) 59 NONLINEAR NATURAL FREQUENCIES OF A TAPERED CANTILEVER BEAM M. Abdel-Jaber, A.A. Al-Qasa,* and M.S. Abdel-Jaber Department of Cvl Engneerng, Faculty

More information

ANALYSIS OF ELECTROMAGNETIC FIELD USING THE CONSTRAINED INTERPOLATION PROFILE METHOD PHÂN TÍCH TRƯỜNG ĐIỆN TỪ SỬ DỤNG PHƯƠNG PHÁP CIP

ANALYSIS OF ELECTROMAGNETIC FIELD USING THE CONSTRAINED INTERPOLATION PROFILE METHOD PHÂN TÍCH TRƯỜNG ĐIỆN TỪ SỬ DỤNG PHƯƠNG PHÁP CIP ANALYSIS OF ELECTROMAGNETIC FIELD USING THE CONSTRAINED INTERPOLATION PROFILE METHOD PHÂN TÍCH TRƯỜNG ĐIỆN TỪ SỬ DỤNG PHƯƠNG PHÁP CIP LÊ VŨ HƯNG Cao đẳng kỹ thuật quốc ga Kushro, Nhật Bản Kushro Natonal

More information

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

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

More information

6.3.4 Modified Euler s method of integration

6.3.4 Modified Euler s method of integration 6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from

More information

THE EFFECT OF BEAM TO COLUMN CONNECTION IN ARC PORTAL FRAME

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

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

In this section is given an overview of the common elasticity models.

In this section is given an overview of the common elasticity models. Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp

More information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information

Solution Thermodynamics

Solution Thermodynamics Soluton hermodynamcs usng Wagner Notaton by Stanley. Howard Department of aterals and etallurgcal Engneerng South Dakota School of nes and echnology Rapd Cty, SD 57701 January 7, 001 Soluton hermodynamcs

More information

CHAPTER 14 GENERAL PERTURBATION THEORY

CHAPTER 14 GENERAL PERTURBATION THEORY CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves

More information

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites 7 Asa-Pacfc Engneerng Technology Conference (APETC 7) ISBN: 978--6595-443- The Two-scale Fnte Element Errors Analyss for One Class of Thermoelastc Problem n Perodc Compostes Xaoun Deng Mngxang Deng ABSTRACT

More information

Tensor Smooth Length for SPH Modelling of High Speed Impact

Tensor Smooth Length for SPH Modelling of High Speed Impact Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru

More information

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Blucher Mechancal Engneerng Proceedngs May 0, vol., num. www.proceedngs.blucher.com.br/evento/0wccm STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Takahko Kurahash,

More information

FUZZY FINITE ELEMENT METHOD

FUZZY FINITE ELEMENT METHOD FUZZY FINITE ELEMENT METHOD RELIABILITY TRUCTURE ANALYI UING PROBABILITY 3.. Maxmum Normal tress Internal force s the shear force, V has a magntude equal to the load P and bendng moment, M. Bendng moments

More information

Chapter 11: Simple Linear Regression and Correlation

Chapter 11: Simple Linear Regression and Correlation Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

One Dimensional Axial Deformations

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

More information

Chapter Eight. Review and Summary. Two methods in solid mechanics ---- vectorial methods and energy methods or variational methods

Chapter Eight. Review and Summary. Two methods in solid mechanics ---- vectorial methods and energy methods or variational methods Chapter Eght Energy Method 8. Introducton 8. Stran energy expressons 8.3 Prncpal of statonary potental energy; several degrees of freedom ------ Castglano s frst theorem ---- Examples 8.4 Prncpal of statonary

More information

Cubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two

Cubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two World Academy of Scence Engneerng and echnology Internatonal Journal of Mathematcal and omputatonal Scences Vol: No:0 00 ubc rgonometrc B-Splne Appled to Lnear wo-pont Boundary Value Problems of Order

More information

Constitutive Modelling of Superplastic AA-5083

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

More information

Handout # 6 (MEEN 617) Numerical Integration to Find Time Response of SDOF mechanical system. and write EOM (1) as two first-order Eqs.

Handout # 6 (MEEN 617) Numerical Integration to Find Time Response of SDOF mechanical system. and write EOM (1) as two first-order Eqs. Handout # 6 (MEEN 67) Numercal Integraton to Fnd Tme Response of SDOF mechancal system State Space Method The EOM for a lnear system s M X + DX + K X = F() t () t = X = X X = X = V wth ntal condtons, at

More information

Chapter 3. Estimation of Earthquake Load Effects

Chapter 3. Estimation of Earthquake Load Effects Chapter 3. Estmaton of Earthquake Load Effects 3.1 Introducton Sesmc acton on chmneys forms an addtonal source of natural loads on the chmney. Sesmc acton or the earthquake s a short and strong upheaval

More information

EN40: Dynamics and Vibrations. Homework 7: Rigid Body Kinematics

EN40: Dynamics and Vibrations. Homework 7: Rigid Body Kinematics N40: ynamcs and Vbratons Homewor 7: Rgd Body Knematcs School of ngneerng Brown Unversty 1. In the fgure below, bar AB rotates counterclocwse at 4 rad/s. What are the angular veloctes of bars BC and C?.

More information

FINITE DIFFERENCE ANALYSIS OF CURVED DEEP BEAMS ON WINKLER FOUNDATION

FINITE DIFFERENCE ANALYSIS OF CURVED DEEP BEAMS ON WINKLER FOUNDATION VOL. 6, NO. 3, MARCH 0 ISSN 89-6608 006-0 Asan Research Publshng Network (ARPN). All rghts reserved. FINITE DIFFERENCE ANALYSIS OF CURVED DEEP BEAMS ON WINKLER FOUNDATION Adel A. Al-Azzaw and Al S. Shaker

More information

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS Unversty of Oulu Student Laboratory n Physcs Laboratory Exercses n Physcs 1 1 APPEDIX FITTIG A STRAIGHT LIE TO OBSERVATIOS In the physcal measurements we often make a seres of measurements of the dependent

More information

I have not received unauthorized aid in the completion of this exam.

I have not received unauthorized aid in the completion of this exam. ME 270 Sprng 2013 Fnal Examnaton Please read and respond to the followng statement, I have not receved unauthorzed ad n the completon of ths exam. Agree Dsagree Sgnature INSTRUCTIONS Begn each problem

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

Second Order Analysis

Second Order Analysis Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information