Formulation of the Stress Distribution Due to a Concentrated Force Acting on the Boundary of Viscoelastic Half-Space

Size: px
Start display at page:

Download "Formulation of the Stress Distribution Due to a Concentrated Force Acting on the Boundary of Viscoelastic Half-Space"

Transcription

1 Formulaion of he Sress Disribuion Due o a Concenraed Force Acing on he Boundar of Viscoelasic Half-Space Yun eng and Debao Zhou Deparmen of Mechanical and Indusrial Engineering Universi of Minnesoa, Duluh MN 8 USA peng04@d.umn.edu and dhou@d.umn.edu. Absrac - This paper presens he derivaive of he sress disribuion in he condiion ha a poin (or concenraed) force acs on he surface of a viscoelasic half-space. This soluion is obained hrough he combinaion of he soluion o he sress disribuion from he forces in boh normal and angenial direcions. The angenial-force viscoelasic problem is defined as he applicaion of a angenial poin force acing on he surface of a viscoelasic half-space and he normal-force viscoelasic problem is relaed o a normal force in he same siuaion. The elasic-viscoelasic correspondence principle is used in he derivaion of he sress disribuion. Two ime-dependen funcions, also called deermining funcions, are uilied o epress he muliplicaion of he eernal force and he ime-dependen viscoelasic maerial properies, i.e. he volumeric and deviaoric relaaion funcions and oisson s raio. Using he boundar condiions and he equilibrium condiions of a half-space, he mahemaical epressions of he deermining funcions can be eplicil obained. The sress disribuion in a viscoelasic maerial under general poin force (wih boh normal and angenial componens) is hen obained. These soluions can be used o formulae he sress disribuion in he crushing or cuing of linear viscoelasic maerials. Ke words: Sress disribuion linear viscoelasic angenial poin force half-space. I INTODUCTION Cuing operaions are involved in he process of man viscoelasic maerials. Modeling he relaionship beween he cuing force and sress disribuion can help predic he fracure due o cuing and idenif differen maerials along he cuing pah. Cuing problem can be modeled as a bel-shaped area force acing on he surface of a half-space. The sud of he sress disribuion due o a poin force will be he firs sep in modeling he cuing sress disribuion. A cuing force can be considered as he resulan of a normal and a angenial componen acing on he conac surface beween he ool and he maerial. The sress disribuion in an elasic half-space due o a normal poin force and ha due o a angenial poin force has been modeled b Boussinesq and Cerrui, respecivel. Correspondingl, he are called Boussinesq s problem and Cerrui s problem. For viscoelasic maerials, Talbl (00) raised he mehod b subsiuing he muliplicaion of eernal force and viscoelasic maerial funcions wih ime-dependen funcions. This paper will concenrae on he formulaion of sress disribuion in a viscoelasic under a angenial poin force. The simplified model is shown in Fig.. Due o he asmmer resuling from he angenial force, onl Caresian coordinaes could be used. This makes our formulaion much more complicaed. Fig. : A angenial force ( ) acs on he surface of a half-space Furhermore, using he soluion formulaed in his paper, a soluion o he problem in viscoelasic maerials under boh normal and angenial poin forces can be obained. This generalied soluion can be used in man furher problems, for eample, he food slicing cu problems formulaed b Zhou and McMurra (0). The remainder of his paper is as follows: angenial-force viscoelasic problem is eplained in Secion II and he soluion o his problem is derived in Secion III. The general soluion is shown in Secion IV. The conclusions are drawn in Secion V. II STATEMENT OF A VISCOELASTIC HALF-SACE UNDE TANENTIAL-FOCE A. Tangenial oin Force In his problem, a reference frame O- is defined on he half-space as shown in Fig., where he boundar of he half-space is a = 0, and a angenial poin force, which could be a funcion of ime, is applied a he origin along he -ais, and he posiive -ais poins owards he inerior of he half-space. The sress soluion o his problem for elasic maerials was given b Cerrui in 88 hrough he use of singulariies from poenial heor. The resuls were also presened b Love (97). The displacemen disribuions a poin (,, ) inside he half-space bod are:

2 e ( ) ( ) u e v e 4 ( ) ( ) v w ( ) ( ) where u e, v e and w e denoe he displacemen in he posiive -, - and -ais direcions in elasic case,, is he eernal angenial force, is he elasic shear modulus and is he oisson s raio. The sress disribuions can hen be obained using he kinemaic equaions and consiuive equaions for elasic maerials. The sresses a poin (,, ) saisf he following boundar condiions for his problem: dd 0 ( ) dd 0 dd 0 ( ) dd 0 dd 0 ( ) dd 0, where he hree equaions in lef column represen he force balance, and he hree equaions in righ column represen he orque balance. B. Viscoelasic Half-Space under Tangenial oin Force In our derivaion, he viscoelasic effec is aken ino consideraion based on he original angenial-force elasic problem. The difference beween elasic and viscoelasic problems lies in he consiuive equaions. In an elasic maerial, he relaionship beween sress and srain can be described b Hooke s Law. However, in a viscoelasic problem, he maerial will show elasic behaviors like solids and also show viscous behaviors like fluids. This changes he relaionship beween sress and srain. During he calculaion of he sress, insead of mulipling srain wih maerial properies, such as Young s Modulus and oisson s raio, he convoluions of srain wih relaaion funcions are used. These relaionships are shown as: d d d d d d d d d, where and are he deviaoric and volumeric relaaion funcions, and / is he mean srain. Therefore, we sae he angenial-force viscoelasic problem as follows: finding ou he soluions o he sress disribuion of a half-space bod under a poin angenial force applied o he surface, wih boundar condiions and sress-srain relaionships saisfied. III SOLUTIONS TO TANENTIAL-FOCE VISCOELASTIC OBLEM A. Displacemens For he linear viscoelasic case, b appling he elasic-viscoelasic correspondence principle and replacing and wih and, The following epressions can be obained: u u u v v v 4 w w w where u, v and w denoe he displacemens in he posiive direcions of he -, - and -ais wih: u u ( ) ( ) v v ( ) ( ) w w. ( ) ( ) In following secions, all derivaions will be conduced wih he erms of and carried ou individuall. B. Normal Srains The normal srains are obained via he derivaive of displacemens wih coordinaes. The obained normal srains are: u u u v v v = 4 w w w where: u ( ) ( ) ( ) u ( ) ( ) ( ) v ( ) ( ) ( ) v ( ) ( ) ( ) w ( ) ( ) w ( ) ( ) ( ) ( ) Thus he mean srain is: 4 and he mean sress is d d d 4

3 C. Shear Srains The shear srains, and can be obained as: v w v w v w = u w u w u w 4 u v u v u v where v w v w 0. u w u w 0. u v ( ) ( ) ( ) u v. ( ) ( ) ( ) Noed is ha he coefficiens of in and are ero, hen he srains can be wrien as 0 = 0 4 u v u v D. Sresses Wih srains known, he sresses can be obained from viscoelasic consiuive equaions (Zhang 994, p. 6) as: d d 4 d 4 ( ) ( ) ( ) d ( ) ( ) ( ) d d 4 d 4 ( ) ( ) ( ) d ( ) ( ) ( ) d d 4 d 4 d d d 4 d d 4 d d 4 ( ) ( ) ( ) d. ( ) ( ) ( ) E. Boundar Condiions Now consider he boundar condiions. Inegraing he sress componen and noicing ha convoluion will no ake par in spaial inegrals, he firs boundar condiion becomes: d dd d 0 4 The above equaion means ha he resulan force in direcion on he O-- plane cancels he eernal force. I is worh o menion ha he inegrals of he oher sress componens are all eros because all of hem are odd funcions of wih he inegraions over,. The value of he above inegral should be ero as required b he balance of force. Then he following epression can be obained: d F. Equilibrium Equaion Bringing all he sresses ino he firs equilibrium equaion E 0 will ield: E d 4 u u v u w d 4 u u v 0, where he coefficiens of d and are: u u v u w u u v u w u u v u u v. ()

4 Thus here is: E d d B eliminaing he same coefficiens conaining he coordinaes, here is: d d 0 () The second equilibrium equaion will ield he same equaion relaionship as given in (). The hird equilibrium equaion 0 has been saisfied auomaicall.. Soluion o Tangenial-Force Viscoelasic roblem and are linearl independen in Equaions () and (). Thus, and can be uniquel solved b solving he Volerra inegral equaions of he second kind (Zhang 994, p.0-6). iven he iniial condiion (0) = 0, all sress and srain componens are ero a = 0 and 0 0, 0 0, he soluion o can be obained. Noed is ha he deailed derivaions are no included in his par for he purpose of concision since he deailed derivaions and he proof of he uniqueness can be found in he paper b Talbl (00). Manipulaing equaions () and (), here is: d 6 Viewing in whole and solving his Volerra inegral equaion, can be obained wih previousl solved. Taking advanage of () and (), he erm of can be eliminaed and he sress disribuions can be wrien as follows: ( ) ( ) ( ) ( ) ( ) ( ) d. ( ) ( ) ( ) d d () Comparing wih classic Cerrui s soluions (Johnson, 98 p.69-70), we found ha he erm ()- *dψ in viscoelasic soluions plas he same role as he erm (-μ) in elasic soluions. IV STESS DISTIBUTION IN A VISCOELASTIC HALF- SACE UNDE ENEAL OINT FOCES We now consider a general force wih boh normal and angenial componens, in which wo angenial forces () and (), along - and - ais respecivel and one normal force () along -ais are applied a he origin. The displacemen and sress disribuions are obained using he superposiion of he displacemen and sress disribuion b (), () and (). The sress and displacemen disribuions for he problem when onl one angenial force ()= () along -ais has been discussed in previous discussions. For disambiguaion, we rewrie he soluion wih proper superscrip (i) (i =, or ) o denoe he applied forces. Therefore we have:,,,,,. where he epressions of he righ hand side are he same as in () onl wih () replaced b (). Fig. : Model for coordinae s ransformaion When here is a single angenial force () applied on he O-- plane a poin O along -ais, he sress disribuion is obained hrough he calculaion of he frame roaion. In his mehod, as shown in Fig., we firs roae he frames around -ais b 90 and denoe he new coordinae ssem as O-. The corresponding Jacob mari is: cos(90 ) sin(90 ) sin(90 ) cos(90 ) The coordinaes are ransformed as ' ' ' ' ' and ' ' ' ' The corresponding sress ensor is ransformed as ' ' ' ' ' ' ' ' ' ' T ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' Correspondingl, he deermining funcions can be rewrien as (Talbl, 00):

5 i i i d, i,,, 0 0 i i i U i d, i,, In coordinaes O-, () is applied a he origin along -ais. Based on he soluion we obained in (), he sress in O- ssem can be obained as: ' ' ' ' ' ' d '( ' ') ' ( ' ') ' ( ' ') ' ' ' ' ' ' ' ' ' ' ' d '( ' ') ' ( ' ') ' ( ' ') ' ' ' ' ' ' (4) ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' d. '( ' ') ' ( ' ') ' ( ' ') ' where ' ' ' ' Afer ransformaion o O- coordinaes, sresses in (4) can be wrien as: d ( ) ( ) ( ) () d. ( ) ( ) ( ) When a single normal force () is applied a poin O and along he posiive -ais, he sress disribuion is shown in (6). Noed is ha his epression has been given b Talbl (00). d ( ) ( ) ( ) d 4 d 0 0. Based on (4), () and (6), we provide he soluion of sress disribuions o he general problem as follows: i i,, i i i i,, i i i i,. i i (6) Similarl, he displacemen disribuions o viscoelasic problem u, v and w as follows: i u u, i i v v, i i w w, i where he componens are: u ( ) ( ) v 4 ( ) ( ) w ( ) ( ) ( ) ( ) u v 4 w ( ) ( ) ( ) ( ) u v 4 ( ) ( ) w The resuls will be used in our fuure research abou he formulaion of he cuings for linear viscoelasic maerials, in which a disribuive force is used insead of a ij r, poin force. We could obain he sress response o a cerain disribuive force b considering he poin force as a funcion of and, replacing b, b, and hen inegraing he corresponding sress or displacemen componens for and in - plane. and b V CONCLUSION The soluions o he sress and displacemen disribuions in a viscoelasic half-space under a poin force are presened in his paper. The mehod is based on he elasic-viscoelasic corresponding principle. The soluion o he displacemen of an elasic half-space under a angenial poin force was used as he displacemen soluion o he viscoelasic problem. Based on he equilibrium equaions and boundar condiions of a viscoelasic half-space, we obained he soluion o he wo ime-dependen deermining funcions via he Volerra inegral equaions of he second kind. Then, he sress disribuion due o a angenial force in -ais direcion is solved based on a frame roaion mehod.

6 Finall, b combining he soluions when here is onl a normal force componen, we ge he resuls for he generalied case, where an arbirar force wih hree nonero componens in, and direcions, is applied. These resuls could be furher used o solve he sress disribuion in a viscoelasic half-space under a disribuive force, b aking he inegral over he area where he disribuive force is applied. VI EFEENCES Boussinesq, J., Applicaion des oeniels a I Eude de l Euqilibre e du Mouvemen des Solides Elasiques, auhier-villars, aris, (88). Love, A. E. H., A Treaise on he Mahemaical Theor of Elasici, 4h Ediion. Cambridge Universi ress (97). Talbl, L. K., Boussinesq's viscoelasic problem on normal concenraed force on a half-space surface, Mechanics of Time- Dependen Maerials. :-9, (00). Xu, Z., Elasici Mechanics (Tan Xing Li Xue) h Ediion, Higher-level Educaion ublicaion, in Chinese, (996). Zhang, C. Y., Viscoelasic Fracure Mechanics (Nian Tan Xing Duan Lie Li Xue), Huahong Universi of Science and Technolog ress, in Chinese, (994). Zhou, D. and McMurra,., Slicing Cus on Food Maerials Using oboic Conrolled aor Blade, Modelling and Simulaion in Engineering, Acceped and o be published in Nov 0, Hindawi ublishing Corporaion, (0). Johnson, K. L., Conac Mechanics, Cambridge Universi ress, (98).

New effective moduli of isotropic viscoelastic composites. Part I. Theoretical justification

New effective moduli of isotropic viscoelastic composites. Part I. Theoretical justification IOP Conference Series: Maerials Science and Engineering PAPE OPEN ACCESS New effecive moduli of isoropic viscoelasic composies. Par I. Theoreical jusificaion To cie his aricle: A A Sveashkov and A A akurov

More information

CH.7. PLANE LINEAR ELASTICITY. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

CH.7. PLANE LINEAR ELASTICITY. Continuum Mechanics Course (MMC) - ETSECCPB - UPC CH.7. PLANE LINEAR ELASTICITY Coninuum Mechanics Course (MMC) - ETSECCPB - UPC Overview Plane Linear Elasici Theor Plane Sress Simplifing Hpohesis Srain Field Consiuive Equaion Displacemen Field The Linear

More information

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow 1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering

More information

Lecture 10: Wave equation, solution by spherical means

Lecture 10: Wave equation, solution by spherical means Lecure : Wave equaion, soluion by spherical means Physical modeling eample: Elasodynamics u (; ) displacemen vecor in elasic body occupying a domain U R n, U, The posiion of he maerial poin siing a U in

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

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

More information

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec

More information

Flow-Induced Vibration Analysis of Supported Pipes with a Crack

Flow-Induced Vibration Analysis of Supported Pipes with a Crack Flow-Induced Vibraion Analsis of Suppored Pipes wih a Crack Jin-Huk Lee, Samer Masoud Al-Said Deparmen of Mechanical Engineering American Universi of Sharjah, UAE Ouline Inroducion and Moivaion Aeroacousicall

More information

Second Order Linear Differential Equations

Second Order Linear Differential Equations Second Order Linear Differenial Equaions Second order linear equaions wih consan coefficiens; Fundamenal soluions; Wronskian; Exisence and Uniqueness of soluions; he characerisic equaion; soluions of homogeneous

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

6.2 Transforms of Derivatives and Integrals.

6.2 Transforms of Derivatives and Integrals. SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

Structural Dynamics and Earthquake Engineering

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

More information

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

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

More information

v A Since the axial rigidity k ij is defined by P/v A, we obtain Pa 3

v A Since the axial rigidity k ij is defined by P/v A, we obtain Pa 3 The The rd rd Inernaional Conference on on Design Engineering and Science, ICDES 14 Pilsen, Czech Pilsen, Republic, Czech Augus Republic, 1 Sepember 1-, 14 In-plane and Ou-of-plane Deflecion of J-shaped

More information

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion

More information

ME 391 Mechanical Engineering Analysis

ME 391 Mechanical Engineering Analysis Fall 04 ME 39 Mechanical Engineering Analsis Eam # Soluions Direcions: Open noes (including course web posings). No books, compuers, or phones. An calculaor is fair game. Problem Deermine he posiion of

More information

Differential Equations

Differential Equations Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding

More information

4.6 One Dimensional Kinematics and Integration

4.6 One Dimensional Kinematics and Integration 4.6 One Dimensional Kinemaics and Inegraion When he acceleraion a( of an objec is a non-consan funcion of ime, we would like o deermine he ime dependence of he posiion funcion x( and he x -componen of

More information

From Particles to Rigid Bodies

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

More information

and v y . The changes occur, respectively, because of the acceleration components a x and a y

and v y . The changes occur, respectively, because of the acceleration components a x and a y Week 3 Reciaion: Chaper3 : Problems: 1, 16, 9, 37, 41, 71. 1. A spacecraf is raveling wih a veloci of v0 = 5480 m/s along he + direcion. Two engines are urned on for a ime of 84 s. One engine gives he

More information

Integration of the equation of motion with respect to time rather than displacement leads to the equations of impulse and momentum.

Integration of the equation of motion with respect to time rather than displacement leads to the equations of impulse and momentum. Inegraion of he equaion of moion wih respec o ime raher han displacemen leads o he equaions of impulse and momenum. These equaions greal faciliae he soluion of man problems in which he applied forces ac

More information

Torsion CHAPTER 7.1 INTRODUCTION

Torsion CHAPTER 7.1 INTRODUCTION CHAPTE 7 Torsion 7. INTODUCTION The orsion of circular shafs has been discussed in elemenar srengh of maerials. There, we were able o obain a soluion o his problem under he assumpion ha he cross-secions

More information

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15. SMT Calculus Tes Soluions February 5,. Le f() = and le g() =. Compue f ()g (). Answer: 5 Soluion: We noe ha f () = and g () = 6. Then f ()g () =. Plugging in = we ge f ()g () = 6 = 3 5 = 5.. There is a

More information

And the solution to the PDE problem must be of the form Π 1

And the solution to the PDE problem must be of the form Π 1 5. Self-Similar Soluions b Dimensional Analsis Consider he diffusion problem from las secion, wih poinwise release (Ref: Bluman & Cole, 2.3): c = D 2 c x + Q 0δ(x)δ() 2 c(x,0) = 0, c(±,) = 0 Iniial release

More information

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

Second Order Linear Differential Equations

Second Order Linear Differential Equations Second Order Linear Differenial Equaions Second order linear equaions wih consan coefficiens; Fundamenal soluions; Wronskian; Exisence and Uniqueness of soluions; he characerisic equaion; soluions of homogeneous

More information

Elementary Differential Equations and Boundary Value Problems

Elementary Differential Equations and Boundary Value Problems Elemenar Differenial Equaions and Boundar Value Problems Boce. & DiPrima 9 h Ediion Chaper 1: Inroducion 1006003 คณ ตศาสตร ว ศวกรรม 3 สาขาว ชาว ศวกรรมคอมพ วเตอร ป การศ กษา 1/2555 ผศ.ดร.อร ญญา ผศ.ดร.สมศ

More information

CHAPTER 12 DIRECT CURRENT CIRCUITS

CHAPTER 12 DIRECT CURRENT CIRCUITS CHAPTER 12 DIRECT CURRENT CIUITS DIRECT CURRENT CIUITS 257 12.1 RESISTORS IN SERIES AND IN PARALLEL When wo resisors are conneced ogeher as shown in Figure 12.1 we said ha hey are conneced in series. As

More information

The Contradiction within Equations of Motion with Constant Acceleration

The Contradiction within Equations of Motion with Constant Acceleration The Conradicion wihin Equaions of Moion wih Consan Acceleraion Louai Hassan Elzein Basheir (Daed: July 7, 0 This paper is prepared o demonsrae he violaion of rules of mahemaics in he algebraic derivaion

More information

( ) 2. Review Exercise 2. cos θ 2 3 = = 2 tan. cos. 2 x = = x a. Since π π, = 2. sin = = 2+ = = cotx. 2 sin θ 2+

( ) 2. Review Exercise 2. cos θ 2 3 = = 2 tan. cos. 2 x = = x a. Since π π, = 2. sin = = 2+ = = cotx. 2 sin θ 2+ Review Eercise sin 5 cos sin an cos 5 5 an 5 9 co 0 a sinθ 6 + 4 6 + sin θ 4 6+ + 6 + 4 cos θ sin θ + 4 4 sin θ + an θ cos θ ( ) + + + + Since π π, < θ < anθ should be negaive. anθ ( + ) Pearson Educaion

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

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

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

More information

Equivalence Problem of the Painlevé Equations

Equivalence Problem of the Painlevé Equations Advances in Pure Mahemaics 0 97-0 hp://ddoiorg/06/apm00 Pulished Online March 0 (hp://wwwscirporg/journal/apm) Equivalence Prolem of he Painlevé Equaions Sopia Khamrod Deparmen of Mahemaics Facul of Science

More information

ASTR415: Problem Set #5

ASTR415: Problem Set #5 ASTR45: Problem Se #5 Curran D. Muhlberger Universi of Marland (Daed: April 25, 27) Three ssems of coupled differenial equaions were sudied using inegraors based on Euler s mehod, a fourh-order Runge-Kua

More information

ln 2 1 ln y x c y C x

ln 2 1 ln y x c y C x Lecure 14 Appendi B: Some sample problems from Boas Here are some soluions o he sample problems assigned for Chaper 8 8: 6 Soluion: We wan o find he soluion o he following firs order equaion using separaion

More information

Optimal Path Planning for Flexible Redundant Robot Manipulators

Optimal Path Planning for Flexible Redundant Robot Manipulators 25 WSEAS In. Conf. on DYNAMICAL SYSEMS and CONROL, Venice, Ialy, November 2-4, 25 (pp363-368) Opimal Pah Planning for Flexible Redundan Robo Manipulaors H. HOMAEI, M. KESHMIRI Deparmen of Mechanical Engineering

More information

4.1.1 Mindlin plates: Bending theory and variational formulation

4.1.1 Mindlin plates: Bending theory and variational formulation Chaper 4 soropic fla shell elemens n his chaper, fia shell elemens are formulaed hrough he assembly of membrane and plae elemens. The exac soluion of a shell approximaed by fia faces compared o he exac

More information

Finite Element Analysis of Structures

Finite Element Analysis of Structures KAIT OE5 Finie Elemen Analysis of rucures Mid-erm Exam, Fall 9 (p) m. As shown in Fig., we model a russ srucure of uniform area (lengh, Area Am ) subjeced o a uniform body force ( f B e x N / m ) using

More information

EE243 Advanced Electromagnetic Theory Lec # 13: Waveguides and sources

EE243 Advanced Electromagnetic Theory Lec # 13: Waveguides and sources Applied M Fall 6, Neureuher Lecure #3 er /8/6 43 Advanced lecromagneic Theor Lec # 3: Waveguides and sources Source Free Region: ecor Poenials A and F Single direcion componen of A and F Give TM and T

More information

The expectation value of the field operator.

The expectation value of the field operator. The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining

More information

Unsteady Flow Problems

Unsteady Flow Problems School of Mechanical Aerospace and Civil Engineering Unseady Flow Problems T. J. Craf George Begg Building, C41 TPFE MSc CFD-1 Reading: J. Ferziger, M. Peric, Compuaional Mehods for Fluid Dynamics H.K.

More information

Outline of Topics. Analysis of ODE models with MATLAB. What will we learn from this lecture. Aim of analysis: Why such analysis matters?

Outline of Topics. Analysis of ODE models with MATLAB. What will we learn from this lecture. Aim of analysis: Why such analysis matters? of Topics wih MATLAB Shan He School for Compuaional Science Universi of Birmingham Module 6-3836: Compuaional Modelling wih MATLAB Wha will we learn from his lecure Aim of analsis: Aim of analsis. Some

More information

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx. . Use Simpson s rule wih n 4 o esimae an () +. Soluion: Since we are using 4 seps, 4 Thus we have [ ( ) f() + 4f + f() + 4f 3 [ + 4 4 6 5 + + 4 4 3 + ] 5 [ + 6 6 5 + + 6 3 + ]. 5. Our funcion is f() +.

More information

Chapter 5 Kinematics

Chapter 5 Kinematics Chaper 5 Kinemaics In he firs place, wha do we mean b ime and space? I urns ou ha hese deep philosophical quesions have o be analzed ver carefull in phsics, and his is no eas o do. The heor of relaivi

More information

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem) Week 1 Lecure Problems, 5 Wha if somehing oscillaes wih no obvious spring? Wha is ω? (problem se problem) Sar wih Try and ge o SHM form E. Full beer can in lake, oscillaing F = m & = ge rearrange: F =

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

LAPLACE TRANSFORM AND TRANSFER FUNCTION

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

More information

Second quantization and gauge invariance.

Second quantization and gauge invariance. 1 Second quanizaion and gauge invariance. Dan Solomon Rauland-Borg Corporaion Moun Prospec, IL Email: dsolom@uic.edu June, 1. Absrac. I is well known ha he single paricle Dirac equaion is gauge invarian.

More information

STRESS ANALYSIS OF POLARIZATION MAINTAINING OPTICAL FIBERS BY THE FINITE ELEMENT METHOD

STRESS ANALYSIS OF POLARIZATION MAINTAINING OPTICAL FIBERS BY THE FINITE ELEMENT METHOD IIUM Engineering Journal, Vol. 1, No. 1, Januar STRESS ANALYSIS OF POLARIZATION MAINTAINING OPTICAL FIBERS BY THE FINITE ELEMENT METHOD M. H. Al Facul of Engineering, Universi of Aleandria, Aleandria 1,

More information

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1. Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of

More information

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

More information

KINEMATICS IN ONE DIMENSION

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

More information

This is an example to show you how SMath can calculate the movement of kinematic mechanisms.

This is an example to show you how SMath can calculate the movement of kinematic mechanisms. Dec :5:6 - Kinemaics model of Simple Arm.sm This file is provided for educaional purposes as guidance for he use of he sofware ool. I is no guaraeed o be free from errors or ommissions. The mehods and

More information

Review of EM and Introduction to FDTD

Review of EM and Introduction to FDTD 1/13/016 5303 lecromagneic Analsis Using Finie Difference Time Domain Lecure #4 Review of M and Inroducion o FDTD Lecure 4These noes ma conain coprighed maerial obained under fair use rules. Disribuion

More information

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation Hea (iffusion) Equaion erivaion of iffusion Equaion The fundamenal mass balance equaion is I P O L A ( 1 ) where: I inpus P producion O oupus L losses A accumulaion Assume ha no chemical is produced or

More information

Fractional Method of Characteristics for Fractional Partial Differential Equations

Fractional Method of Characteristics for Fractional Partial Differential Equations Fracional Mehod of Characerisics for Fracional Parial Differenial Equaions Guo-cheng Wu* Modern Teile Insiue, Donghua Universiy, 188 Yan-an ilu Road, Shanghai 51, PR China Absrac The mehod of characerisics

More information

Numerical Dispersion

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

More information

MATH 31B: MIDTERM 2 REVIEW. x 2 e x2 2x dx = 1. ue u du 2. x 2 e x2 e x2] + C 2. dx = x ln(x) 2 2. ln x dx = x ln x x + C. 2, or dx = 2u du.

MATH 31B: MIDTERM 2 REVIEW. x 2 e x2 2x dx = 1. ue u du 2. x 2 e x2 e x2] + C 2. dx = x ln(x) 2 2. ln x dx = x ln x x + C. 2, or dx = 2u du. MATH 3B: MIDTERM REVIEW JOE HUGHES. Inegraion by Pars. Evaluae 3 e. Soluion: Firs make he subsiuion u =. Then =, hence 3 e = e = ue u Now inegrae by pars o ge ue u = ue u e u + C and subsiue he definiion

More information

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

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

More information

Theory of! Partial Differential Equations!

Theory of! Partial Differential Equations! hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

Two Coupled Oscillators / Normal Modes

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

More information

MECHANICS OF MATERIALS Poisson s Ratio

MECHANICS OF MATERIALS Poisson s Ratio Poisson s Raio For a slender bar subjeced o axial loading: ε x x y 0 The elongaion in he x-direcion i is accompanied by a conracion in he oher direcions. Assuming ha he maerial is isoropic (no direcional

More information

Chapter 8 The Complete Response of RL and RC Circuits

Chapter 8 The Complete Response of RL and RC Circuits Chaper 8 The Complee Response of RL and RC Circuis Seoul Naional Universiy Deparmen of Elecrical and Compuer Engineering Wha is Firs Order Circuis? Circuis ha conain only one inducor or only one capacior

More information

The Arcsine Distribution

The Arcsine Distribution The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

Check in: 1 If m = 2(x + 1) and n = find y when. b y = 2m n 2

Check in: 1 If m = 2(x + 1) and n = find y when. b y = 2m n 2 7 Parameric equaions This chaer will show ou how o skech curves using heir arameric equaions conver arameric equaions o Caresian equaions find oins of inersecion of curves and lines using arameric equaions

More information

Polymer Engineering (MM3POE)

Polymer Engineering (MM3POE) Polymer Engineering (MM3POE) VISCOELASTICITY hp://www.noingham.ac.uk/~eazacl/mm3poe Viscoelasiciy 1 Conens Wha is viscoelasiciy? Fundamenals Creep & creep recovery Sress relaxaion Modelling viscoelasic

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details! MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his

More information

Theory of Elasticity Ct 5141 Direct Methods

Theory of Elasticity Ct 5141 Direct Methods Delf Universi of Technolog Facul of Civil Engineering and Geosciences Theor of Elasici C 54 Direc Mehods Prof.dr.ir. J. Blaauwendraad June 003 Las updae: April 004 C 54 Acknowledgemen In wriing and updaing

More information

Motion along a Straight Line

Motion along a Straight Line chaper 2 Moion along a Sraigh Line verage speed and average velociy (Secion 2.2) 1. Velociy versus speed Cone in he ebook: fer Eample 2. Insananeous velociy and insananeous acceleraion (Secions 2.3, 2.4)

More information

Combined Bending with Induced or Applied Torsion of FRP I-Section Beams

Combined Bending with Induced or Applied Torsion of FRP I-Section Beams Combined Bending wih Induced or Applied Torsion of FRP I-Secion Beams MOJTABA B. SIRJANI School of Science and Technology Norfolk Sae Universiy Norfolk, Virginia 34504 USA sirjani@nsu.edu STEA B. BONDI

More information

Second-Order Differential Equations

Second-Order Differential Equations WWW Problems and Soluions 3.1 Chaper 3 Second-Order Differenial Equaions Secion 3.1 Springs: Linear and Nonlinear Models www m Problem 3. (NonlinearSprings). A bod of mass m is aached o a wall b means

More information

Variational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations

Variational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 1, Issue 6 Ver. II (Nov - Dec. 214), PP 48-54 Variaional Ieraion Mehod for Solving Sysem of Fracional Order Ordinary Differenial

More information

6. 6 v ; degree = 7; leading coefficient = 6; 7. The expression has 3 terms; t p no; subtracting x from 3x ( 3x x 2x)

6. 6 v ; degree = 7; leading coefficient = 6; 7. The expression has 3 terms; t p no; subtracting x from 3x ( 3x x 2x) 70. a =, r = 0%, = 0. 7. a = 000, r = 0.%, = 00 7. a =, r = 00%, = 7. ( ) = 0,000 0., where = ears 7. ( ) = + 0.0, where = weeks 7 ( ) =,000,000 0., where = das 7 = 77. = 9 7 = 7 geomeric 0. geomeric arihmeic,

More information

Non-uniform circular motion *

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

More information

A Note on the Equivalence of Fractional Relaxation Equations to Differential Equations with Varying Coefficients

A Note on the Equivalence of Fractional Relaxation Equations to Differential Equations with Varying Coefficients mahemaics Aricle A Noe on he Equivalence of Fracional Relaxaion Equaions o Differenial Equaions wih Varying Coefficiens Francesco Mainardi Deparmen of Physics and Asronomy, Universiy of Bologna, and he

More information

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4. PHY1 Elecriciy Topic 7 (Lecures 1 & 11) Elecric Circuis n his opic, we will cover: 1) Elecromoive Force (EMF) ) Series and parallel resisor combinaions 3) Kirchhoff s rules for circuis 4) Time dependence

More information

Summary of shear rate kinematics (part 1)

Summary of shear rate kinematics (part 1) InroToMaFuncions.pdf 4 CM465 To proceed o beer-designed consiuive equaions, we need o know more abou maerial behavior, i.e. we need more maerial funcions o predic, and we need measuremens of hese maerial

More information

Theory of! Partial Differential Equations-I!

Theory of! Partial Differential Equations-I! hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-I! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his

More information

NMR Spectroscopy: Principles and Applications. Nagarajan Murali 1D - Methods Lecture 5

NMR Spectroscopy: Principles and Applications. Nagarajan Murali 1D - Methods Lecture 5 NMR pecroscop: Principles and Applicaions Nagarajan Murali D - Mehods Lecure 5 D-NMR To full appreciae he workings of D NMR eperimens we need o a leas consider wo coupled spins. omeimes we need o go up

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

More information

Chapter 6. Systems of First Order Linear Differential Equations

Chapter 6. Systems of First Order Linear Differential Equations Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh

More information

Math 334 Fall 2011 Homework 11 Solutions

Math 334 Fall 2011 Homework 11 Solutions Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then

More information

Math 333 Problem Set #2 Solution 14 February 2003

Math 333 Problem Set #2 Solution 14 February 2003 Mah 333 Problem Se #2 Soluion 14 February 2003 A1. Solve he iniial value problem dy dx = x2 + e 3x ; 2y 4 y(0) = 1. Soluion: This is separable; we wrie 2y 4 dy = x 2 + e x dx and inegrae o ge The iniial

More information

Boundedness and Stability of Solutions of Some Nonlinear Differential Equations of the Third-Order.

Boundedness and Stability of Solutions of Some Nonlinear Differential Equations of the Third-Order. Boundedness Sabili of Soluions of Some Nonlinear Differenial Equaions of he Third-Order. A.T. Ademola, M.Sc. * P.O. Arawomo, Ph.D. Deparmen of Mahemaics Saisics, Bowen Universi, Iwo, Nigeria. Deparmen

More information

CHAPTER 2: Mathematics for Microeconomics

CHAPTER 2: Mathematics for Microeconomics CHAPTER : Mahemaics for Microeconomics The problems in his chaper are primarily mahemaical. They are inended o give sudens some pracice wih he conceps inroduced in Chaper, bu he problems in hemselves offer

More information

3.3 Internal Stress. Cauchy s Concept of Stress

3.3 Internal Stress. Cauchy s Concept of Stress INTERNL TRE 3.3 Inernal ress The idea of sress considered in 3.1 is no difficul o concepualise since objecs ineracing wih oher objecs are encounered all around us. more difficul concep is he idea of forces

More information

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

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

More information

Modified one dimensional elastic wave equations that retain time synchronization under spatial coordinate transformations

Modified one dimensional elastic wave equations that retain time synchronization under spatial coordinate transformations Modified one dimensional elasic wave equaions ha reain ime synchronizaion under spaial coordinae ransformaions Ruiwen Yao, Zhihai Xiang * AML, Deparmen of Engineering Mechanics, Tsinghua Universiy, Beijing

More information

Solutions to Assignment 1

Solutions to Assignment 1 MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we

More information

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits EEE25 ircui Analysis I Se 4: apaciors, Inducors, and Firs-Order inear ircuis Shahriar Mirabbasi Deparmen of Elecrical and ompuer Engineering Universiy of Briish olumbia shahriar@ece.ubc.ca Overview Passive

More information

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

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

More information

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,

More information

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

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

More information

Experiment 123 Determination of the sound wave velocity with the method of Lissajous figures

Experiment 123 Determination of the sound wave velocity with the method of Lissajous figures perimen 3 Deerminaion of he sound wave veloci wih he mehod of Lissajous figures The aim of he eercise To sud acousic wave propagaion in he air To deermine of he sound wave veloci in he air Mehodolog of

More information