DOUBLE PENDULUM VIBRATION MOTION IN FLUID FLOW

Size: px
Start display at page:

Download "DOUBLE PENDULUM VIBRATION MOTION IN FLUID FLOW"

Transcription

1 ENGINEERING FOR RURA DEEOPMENT Jelgava, DOUBE PENDUUM IBRATION MOTION IN FUID FOW Janis iba, Maris Eiuks, Martins Irbe Riga Technical University, atvia Abstract. The paper analyses the ouble penulum motion in the vertical plane interacting with the lui (air or water) low. The penulum consists o two plates, astene together with each other an with the ounation by the axes. The case was investigate where the lui lows in a horizontal irection but the plates move aroun the pivots in the vertical plane. It is assume that there is a laminar lui low. Double penulum interaction with the low is escribe as localize particles interaction with the plates, using the mass centre an angular momentum exchange theorems o mechanical systems. A system with two egrees o reeom is investigate. Interaction with lui is approximate by a square step epenence on the low relative velocity taking into account the irection o motion. The resulting system o equations is use to simulate the motion. The three cases o motion are analyse: ree oscillations against a stationary lui to etermine experimentally all physical parameters o the interaction; auto vibration motion at a constant low velocity; time variable low inuce oscillations. Some interesting results o the plate motion graphics are given. Keywors: ouble penulum, plate interaction with lui, low-inuce vibrations. Introuction Object vibrations in surrouning environment inclue motion insie lui at low spee region with a range when moving parts o elements change the velocity irection. Thereore, it is necessary to analyze the mechanical system motion as a ouble penulum () vibration in lui low () (Fig. ). Nature observation shows that some o the leaves o the trees start penulum motion (rattle) in winy weather. We can also observe water plant vibrations in the river by water stream lows. It is important to check up such vibrations or untraitional energy prouction (without propeller system) []. For that it is necessary investigating some simpler systems like a ouble penulum motion using the classical mechanics theory [-]. This report eals with analysis o two plates astene together with each other an with the ounation connecte to the pivots (that is, with the axis perpenicular to the irection o movement). The ouble penulum an lui interaction orces are analyze an their reuction to the penulum pivot points is oun out. Ater that, consiering win interaction orces ierential equations or the ouble penulum motion analysis are given. Some results o moelling are shown. Fig.. Double penulum in lui low: ouble penulum; horizontally low Moel o object The object consists o two rectangular plates an (Fig..). The plates are rotate aroun the horizontal axes Oz an Az. As generalize coorinates angles φ an φ are use. In the plates or the mass centers an vertically gravity orces mg an mg are applie respectively. The rotary axes Oz an Az can be place asymmetrically with istances an rom the plates ens. Flui low is horizontal with a spee. The analysis o the motion o the ouble penulum system by

2 ENGINEERING FOR RURA DEEOPMENT Jelgava, ierential equations requires ining the lui low interaction orces. They are obtaine in the next section. Fig.. Plane moel o penulum: penulum rotating part; penulum plane motion part; horizontal low; φ, φ generalize coorinates Reuction o lui low interaction orces Flui low interaction orces epen on the relative spee o the contacts points. The main interaction is the square unction o transverse spee components, perpenicular to the plate. Thereore, it is necessary to in out the point velocities projections on the perpenicular irections AO an AB (Fig..). It means that or the plate AO expression rom two graphs an (Fig. ) can be oun. Accoringly, or the plate AB all three graphs, an can be use. In this ouble penulum task or orces reuction it is convenient to use the points O an A as reuction centers. Then, the lui low principal vectors R, R an the principal moments M, M can be oun rom the ollowing integrals (or Fig. ): R M, Fin M R ; () R M, Fin M R. ()

3 ENGINEERING FOR RURA DEEOPMENT Jelgava, Here, & ; & ; where, length o plates;, angular velocities;, constants, incluing rag coeicients, ensity o lui; with o plates. I the rotation axis o the plates are not at the ens, the low principal vectors R, R an the principal moments M, M can be oun by separate length (- till, an - till ) rom ormulas () an (), taking into account eviations an (see Fig..). Fig.. Graphics o ouble penulum points velocities components: graphic o irst penulum rotation velocity component; secon penulum relative rotation velocity component; horizontal low component; R, R principal vectors; M, M principal moments Dierential equations o penulum motion Dierential equations o motion can be oun, or example, by the D Alembert s principle o the given mechanical system (see Fig..). Sum o orces moments against the point A or only the secon plate gives: M φaon rsin ( J m r ) && m && r mg rsin. Sum o orces moment against the point O or both plates gives: () 6

4 ENGINEERING FOR RURA DEEOPMENT Jelgava, Here (see Fig..): M M R ( J m r ) && ( J m r m r ) && φan sin m( r ) φaon rsin φaon mg ( r sin Mφ J && ; Mφ J && ; φaon m & ; mg r sin rsin( ). φot m r && ; φat m r&& ; φaot m &&, where J, J central moments o inertia mass m, m; r, r istance rom center masses till pins; length o irst plate; &, & angular accelerations o plates. φon m r & ; && φan m r & ; () Fig.. D Alembert s principle orces: (mg, mg, R, R, M, M) active orces an their moments;, M components o orces o inertia an their moments From two equations () an () the angular accelerations&, & can be oun. Then, by integrations &, & the motion parameters can be calculate, using numerical computer calculation methos (or example, Euler metho). There exists an aitional possibility o analysis o this motion by special lui low analysis computer programs, or example, Working Moel. 7

5 ENGINEERING FOR RURA DEEOPMENT Jelgava, Motion moelling In this research the investigation was mae by soli boies plane motion calculation computer program Working Moel D. Some interesting graphics o motion are shown in Fig omments are given below or each igure. 5 Fig. 5. Transitional ampe vibration motion o penulum: angular velocity o irst plate as time unction; initial position o ouble penulum; mile position o penulum; close en rest position o penulum; 5 angular velocity o secon plate (see Fig. ) 5 Fig. 6. Stationary vibration motion o penulum with large amplitue: angular velocity o irst plate as time unction; angular velocity o secon plate as time unction; initial position o ouble penulum; positions o rotation axis are with eccentricities an (see Fig. ); 5 experiment in win tunnel with low velocity about 6 m sec - Fig. 7. Biurcation o vibration motion (with small amplitue): angular velocities; initial position o ouble penulum 8

6 ENGINEERING FOR RURA DEEOPMENT Jelgava, Results an iscussion Analysis o the moeling results shows that it is possible to obtain a stationary oscillations motion in the ouble penulum system. For that the axis o rotation can be isplace by eccentricity, against the en points o the plates. This means that lui lows at a constant spee in the observe system generating perioically vibration motion. In this system the existence o vibration biurcations shoul be checke out. This question requires aitional investigations an is not covere by the given article. onclusions. Double penulum vibration system by using lui low interaction can be use or untraitional energy prouction (without propeller system).. Stationary vibration motion epens on special penulum rotation axis isplacements.. Theory valiation was checke by the win tunnel. Reerences. iba J., itols D., Grischenko M., Eiuks M., Kulikovskis G. Investigation o unerwater robotic system or object motion propulsion an energy generation. Proceeings o International conerence Engineering or rural evelopment, 9-.5., Jelgava, atvia, pp Martynyuk A. A. an Nikitina N.. The Theory o Motion o a Double Mathematical Penulum, International Applie Mechanics, ol. 6, No. 9,, pp Double Penulum rom Eric Weisstein s. Worl o Physics scienceworl.wolram.com Mechanics Penula [online][..] Available at: physics/doublepenulum.html. Double Penulum. [online][..] Available at: penulum.html 5. Arnol,. I. Problem in Mathematical Methos o lassical Mechanics, n e. New York: Springer-erlag, p. 9, Wells, D. A. Theory an Problems o agrangian Dynamics. New York: McGraw-Hill, pp. -,, an -, arry M. Bates; James M. Nester. On D Alembert s Principle. ommunications in Mathematics (). olume: 9, Issue:, page ISSN:

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

International Conference KNOWLEDGE-BASED ORGANIZATION Vol. XXIII No

International Conference KNOWLEDGE-BASED ORGANIZATION Vol. XXIII No International Conference KNOWLEDGE-BAED ORGANIZATION Vol. XXIII No 3 2017 METHOD FOR DETERMINATION OF THE PARAMETER OF A MACHINE GUN UPENION MOUNTED ON AN ARMOURED VEHICLE vilen PIRDONOV, vilen TEFANOV

More information

DOUBLE PENDULUM MOTION ANALYSIS IN VARIABLE FLUID FLOW

DOUBLE PENDULUM MOTION ANALYSIS IN VARIABLE FLUID FLOW ENGINEERING FOR RURAL DEVELOPMENT Jelgava, 5.-7.5.16. DOUBLE PENDULUM MOTION ANALYSIS IN VARIABLE FLUID FLOW Oskars Ozolins, Valters Cipruss, Janis Viba, Olegs Jakovlevs Riga Technical University, Latvia

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

AN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A

AN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A AN INTRODUCTION TO AIRCRAFT WIN FLUTTER Revision A By Tom Irvine Email: tomirvine@aol.com January 8, 000 Introuction Certain aircraft wings have experience violent oscillations uring high spee flight.

More information

Product and Quotient Rules and Higher-Order Derivatives. The Product Rule

Product and Quotient Rules and Higher-Order Derivatives. The Product Rule 330_003.q 11/3/0 :3 PM Page 119 SECTION.3 Prouct an Quotient Rules an Higher-Orer Derivatives 119 Section.3 Prouct an Quotient Rules an Higher-Orer Derivatives Fin the erivative o a unction using the Prouct

More information

Numerical Integrator. Graphics

Numerical Integrator. Graphics 1 Introuction CS229 Dynamics Hanout The question of the week is how owe write a ynamic simulator for particles, rigi boies, or an articulate character such as a human figure?" In their SIGGRPH course notes,

More information

MATHEMATICAL MODELLING OF INVERTED PENDULUM WITH DISTURBANCE INPUT

MATHEMATICAL MODELLING OF INVERTED PENDULUM WITH DISTURBANCE INPUT Mathematical Moelling O Inverte Penulum With Disturbance Input MATHEMATICAL MODELLING OF INVERTED PENDULUM WITH DISTURBANCE INPUT PANKAJ WARAK Stuent, ME Control Systems, KKW COE &Research, Nashik Abstract

More information

u t v t v t c a u t b a v t u t v t b a

u t v t v t c a u t b a v t u t v t b a Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying

More information

Numerical Investigation of Non-Stationary Parameters on Effective Phenomena of a Pitching Airfoil at Low Reynolds Number

Numerical Investigation of Non-Stationary Parameters on Effective Phenomena of a Pitching Airfoil at Low Reynolds Number Journal of Applie Flui Mechanics, Vol. 9, No. 2, pp. 643-651, 2016. Available online at www.jafmonline.net, ISSN 1735-3572, EISSN 1735-3645. Numerical Investigation of Non-Stationary Parameters on Effective

More information

Chapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies

Chapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies hapter- Steay Stoes flow aroun eforme sphere. class of oblate axi-symmetric boies. General In physical an biological sciences, an in engineering, there is a wie range of problems of interest lie seimentation

More information

Gyroscopic matrices of the right beams and the discs

Gyroscopic matrices of the right beams and the discs Titre : Matrice gyroscopique es poutres roites et es i[...] Date : 15/07/2014 Page : 1/16 Gyroscopic matrices of the right beams an the iscs Summary: This ocument presents the formulation of the matrices

More information

f sends 1 to 2, so f 1 sends 2 back to 1. f sends 2 to 4, so f 1 sends 4 back to 2. f 1 = { (2,1), (4,2), (6,3), (1,4), (3,5), (5,6) }.

f sends 1 to 2, so f 1 sends 2 back to 1. f sends 2 to 4, so f 1 sends 4 back to 2. f 1 = { (2,1), (4,2), (6,3), (1,4), (3,5), (5,6) }. .3 Inverse Functions an their Derivatives In this unit you ll review inverse unctions, how to in them, an how the graphs o unctions an their inverses are relate geometrically. Not all unctions can be unone,

More information

2.3 Product and Quotient Rules and Higher-Order Derivatives

2.3 Product and Quotient Rules and Higher-Order Derivatives Chapter Dierentiation. Prouct an Quotient Rules an Higher-Orer Derivatives Fin the erivative o a unction using the Prouct Rule. Fin the erivative o a unction using the Quotient Rule. Fin the erivative

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity

1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity AP Physics Multiple Choice Practice Electrostatics 1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity. A soli conucting sphere is given a positive charge Q.

More information

Longitudinal Waves in a Rotating Solid Cylinder Immersed in an Inviscid Fluid Using Chebyshev Polynomial

Longitudinal Waves in a Rotating Solid Cylinder Immersed in an Inviscid Fluid Using Chebyshev Polynomial Longituinal Waves in a Rotating Soli Cyliner Immerse in an Invisci Flui Using Chebyshev Polynomial R. Selvamani *1, P.Ponnusamy 1 Department o Mathematics, Karunya University, Coimbatore, TamilNau, Inia

More information

Research on Thin Film Thickness Uniformity for Deposition of Rectangular Planar Sputtering Target

Research on Thin Film Thickness Uniformity for Deposition of Rectangular Planar Sputtering Target Available online at www.scienceirect.com Physics Proceia 3 (01 ) 903 913 18 th International Vacuum Congress (IVC-18) Research on Thin Film Thickness Uniormity or Deposition o Rectangular Planar puttering

More information

A-level PHYSICS A PHYA4/1. Unit 4 Fields and Further Mechanics. Section A. Monday 20 June 2016 Morning

A-level PHYSICS A PHYA4/1. Unit 4 Fields and Further Mechanics. Section A. Monday 20 June 2016 Morning Please write clearly in block capitals. entre number aniate number Surname Forename(s) aniate signature -level PHYSIS Unit 4 Fiels an Further Mechanics Section Monay 20 June 2016 Morning Materials In aition

More information

12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes

12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes Differentiation of vectors 12.5 Introuction The area known as vector calculus is use to moel mathematically a vast range of engineering phenomena incluing electrostatics, electromagnetic fiels, air flow

More information

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann Gravitation as the result of the reintegration of migrate electrons an positrons to their atomic nuclei. Osvalo Domann oomann@yahoo.com (This paper is an extract of [6] liste in section Bibliography.)

More information

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

More information

6. Friction and viscosity in gasses

6. Friction and viscosity in gasses IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner

More information

Class Notes for Advanced Dynamics (MEAM535)

Class Notes for Advanced Dynamics (MEAM535) Class Notes f Avance Dynamics MEAM535 Michael A. Carchii December 9, 9 Equilibrium Points, Small Perturbations Linear Stability Analysis The following notes were initially base aroun the text entitle:

More information

Average value of position for the anharmonic oscillator: Classical versus quantum results

Average value of position for the anharmonic oscillator: Classical versus quantum results verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive

More information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 4. Electrostatics of Macroscopic Media Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1

More information

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann Gravitation as the result of the reintegration of migrate electrons an positrons to their atomic nuclei. Osvalo Domann oomann@yahoo.com (This paper is an extract of [6] liste in section Bibliography.)

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

CIRCULAR MOTION AND SHM: Solutions to Higher Level Questions

CIRCULAR MOTION AND SHM: Solutions to Higher Level Questions CIRCULAR MOTION AND SHM: Solutions to Higher Level Questions ****ALL QUESTIONS-ANSWERS ARE HIGHER LEVEL*** Solutions 015 Question 6 (i) Explain what is meant by centripetal force. The force - acting in

More information

Week 8 Lecture: Concepts of Quantum Field Theory (QFT) Klein-Gordon Green s Functions and Raising/Lowering Operators

Week 8 Lecture: Concepts of Quantum Field Theory (QFT) Klein-Gordon Green s Functions and Raising/Lowering Operators Week 8 Lecture: Concepts of Quantum Fiel Theory (QFT) Anrew Forrester February 29, 2008 Klein-Goron Green s Functions an aising/lowering Operators This Week s Questions How o the Green s functions of the

More information

CENTURION UNIVERSITY OF TECHNOLOGY & MANAGEMENT,ODISHA CUEE-2015

CENTURION UNIVERSITY OF TECHNOLOGY & MANAGEMENT,ODISHA CUEE-2015 CENTURION UNIVERSITY OF TECHNOLOGY & MANAGEMENT,ODISHA CUEE-015 PHYSICS 1. The imensional formula of angular momentum is a) ML T - b) MLT - c) MLT -1 ) ML T -1. If A B = B A, then the angle between A an

More information

Nonlinear Backstepping Control of Permanent Magnet Synchronous Motor (PMSM)

Nonlinear Backstepping Control of Permanent Magnet Synchronous Motor (PMSM) International ournal o Systems Control (Vol.1-010/Iss.1) Merzoug an Benalla / Nonlinear Backstepping Control o ermanent... / pp. 30-34 Nonlinear Backstepping Control o ermanent Magnet Synchronous Motor

More information

Tutorial Test 5 2D welding robot

Tutorial Test 5 2D welding robot Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

More information

Vectors in two dimensions

Vectors in two dimensions Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication

More information

Computed Tomography Notes, Part 1. The equation that governs the image intensity in projection imaging is:

Computed Tomography Notes, Part 1. The equation that governs the image intensity in projection imaging is: Noll 6 CT Notes : Page Compute Tomograph Notes Part Challenges with Projection X-ra Sstems The equation that governs the image intensit in projection imaging is: z I I ep μ z Projection -ra sstems are

More information

INFLUENCE OF SURFACE ROUGHNESS THROUGH A SERIES OF FLOW FACTORS ON THE PERFORMANCE OF A LONGITUDINALLY ROUGH FINITE SLIDER BEARING

INFLUENCE OF SURFACE ROUGHNESS THROUGH A SERIES OF FLOW FACTORS ON THE PERFORMANCE OF A LONGITUDINALLY ROUGH FINITE SLIDER BEARING ANNALS of Faculty Engineering Huneoara International Journal of Engineering Tome XIV [2016] Fascicule 2 [May] ISSN: 1584-2665 [print; online] ISSN: 1584-2673 [CD-Rom; online] a free-accessmultiisciplinarypublication

More information

DYNAMIC PERFORMANCE OF RELUCTANCE SYNCHRONOUS MACHINES

DYNAMIC PERFORMANCE OF RELUCTANCE SYNCHRONOUS MACHINES Annals of the University of Craiova, Electrical Engineering series, No 33, 9; ISSN 184-485 7 TH INTERNATIONAL CONFERENCE ON ELECTROMECHANICAL AN POWER SYSTEMS October 8-9, 9 - Iaşi, Romania YNAMIC PERFORMANCE

More information

Day 4: Motion Along a Curve Vectors

Day 4: Motion Along a Curve Vectors Day 4: Motion Along a Curve Vectors I give my stuents the following list of terms an formulas to know. Parametric Equations, Vectors, an Calculus Terms an Formulas to Know: If a smooth curve C is given

More information

Debond crack growth in fatigue along fiber in UD composite with broken fibers

Debond crack growth in fatigue along fiber in UD composite with broken fibers Debon crack growth in atigue along iber in UD composite with broken ibers Johannes Eitzenberger an Janis arna Dept o Applie Physics an Mechanical Engineering Lulea University o Technology, Sween SE 97

More information

Estimation amount of snow deposits on the road

Estimation amount of snow deposits on the road Estimation amount of snow eposits on the roa Olga. Glaysheva oronezh State University of Architecture an Civil Engineering, Russia Email: glaov@ox.vsi.ru ABSTRACT The article uner consieration gives the

More information

and from it produce the action integral whose variation we set to zero:

and from it produce the action integral whose variation we set to zero: Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine

More information

Wheel Dynamic Response to Rolling Excitation on Road Weiji Wang1, a*

Wheel Dynamic Response to Rolling Excitation on Road Weiji Wang1, a* Secon International Conference on Mechanics, Materials an Structural Engineering (ICMMSE 7) Wheel Dynamic Response to Rolling Excitation on Roa Weiji Wang, a* Department of Engineering an Design, University

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek

More information

ARCH 614 Note Set 5 S2012abn. Moments & Supports

ARCH 614 Note Set 5 S2012abn. Moments & Supports RCH 614 Note Set 5 S2012abn Moments & Supports Notation: = perpenicular istance to a force from a point = name for force vectors or magnitue of a force, as is P, Q, R x = force component in the x irection

More information

Analysis of Instantaneous Center of Zero Acceleration of Rigid Body in Planar Motion

Analysis of Instantaneous Center of Zero Acceleration of Rigid Body in Planar Motion Moern pplie Science pril, 009 nalysis of Instantaneous Center of Zero cceleration of Rigi Boy in Planar Motion Haibin Sun (Corresponing author) & Tingting Liu Department of Physics an Electronics Taishan

More information

PHYS1169: Tutorial 8 Solutions

PHYS1169: Tutorial 8 Solutions PHY69: Tutorial 8 olutions Wae Motion ) Let us consier a point P on the wae with a phase φ, so y cosϕ cos( x ± ωt) At t0, this point has position x0, so ϕ x0 ± ωt0 Now, at some later time t, the position

More information

Lecture 27: Generalized Coordinates and Lagrange s Equations of Motion

Lecture 27: Generalized Coordinates and Lagrange s Equations of Motion Lecture 27: Generalize Coorinates an Lagrange s Equations of Motion Calculating T an V in terms of generalize coorinates. Example: Penulum attache to a movable support 6 Cartesian Coorinates: (X, Y, Z)

More information

Vibration Analysis of Railway Tracks Forced by Distributed Moving Loads

Vibration Analysis of Railway Tracks Forced by Distributed Moving Loads IJR International Journal of Railway Vol. 6, No. 4 / December 13, pp. 155-159 The Korean Society for Railway Vibration Analysis of Railway Tracks Force by Distribute Moving Loas Sinyeob Lee*, Dongkyu Kim*,

More information

Second Major Solution Q1. The three capacitors in the figure have an equivalent capacitance of 2.77 µf. What is C 2?

Second Major Solution Q1. The three capacitors in the figure have an equivalent capacitance of 2.77 µf. What is C 2? Secon Major Solution Q1. The three capacitors in the figure have an equivalent capacitance of.77 µf. What is C? C 4.0 µf.0 µf A) 7 µf B) µf C) 4 µf D) 3 µf E) 6 µf Q. When the potential ifference across

More information

Chapter 14: Vibration Generations Mechanisms: Self Excited Vibration

Chapter 14: Vibration Generations Mechanisms: Self Excited Vibration hapter 14: ibration Generations echanisms: Self Ecite ibration Introuction: Self-ecite systems begin to vibrate of their own accor spontaneously, the amplitue increasing until some nonlinear effect limits

More information

MATH section 2.3 Basic Differentiation Formulas Page 1 of 5

MATH section 2.3 Basic Differentiation Formulas Page 1 of 5 MATH 0100 section. Basic Dierentiation Formulas Page 1 o The tetbook is using Leibniz notation or this section. I ll continue to use this notation when we are using ormulas an rules. Notation: This means,

More information

A Short Note on Self-Similar Solution to Unconfined Flow in an Aquifer with Accretion

A Short Note on Self-Similar Solution to Unconfined Flow in an Aquifer with Accretion Open Journal o Flui Dynamics, 5, 5, 5-57 Publishe Online March 5 in SciRes. http://www.scirp.org/journal/oj http://x.oi.org/.46/oj.5.57 A Short Note on Sel-Similar Solution to Unconine Flow in an Aquier

More information

Solutions to the Exercises of Chapter 9

Solutions to the Exercises of Chapter 9 9A. Vectors an Forces Solutions to the Exercises of Chapter 9. F = 5 sin 5.9 an F = 5 cos 5 4.8.. a. By the Pythagorean theorem, the length of the vector from to (, ) is + = 5. So the magnitue of the force

More information

LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR

LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR 5 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR Tatsunori Kaneko* an Yasuto Asano* * Department of Mechanical Engineering,

More information

This section outlines the methodology used to calculate the wave load and wave wind load values.

This section outlines the methodology used to calculate the wave load and wave wind load values. COMPUTERS AND STRUCTURES, INC., JUNE 2014 AUTOMATIC WAVE LOADS TECHNICAL NOTE CALCULATION O WAVE LOAD VALUES This section outlines the methoology use to calculate the wave loa an wave win loa values. Overview

More information

PERMANENT MAGNETS CHAPTER MAGNETIC POLES AND BAR MAGNETS

PERMANENT MAGNETS CHAPTER MAGNETIC POLES AND BAR MAGNETS CHAPTER 6 PERAET AGET 6. AGETIC POLE AD BAR AGET We have seen that a small current-loop carrying a current i, prouces a magnetic fiel B o 4 ji ' at an axial point. Here p ia is the magnetic ipole moment

More information

Improved Geoid Model for Syria Using the Local Gravimetric and GPS Measurements 1

Improved Geoid Model for Syria Using the Local Gravimetric and GPS Measurements 1 Improve Geoi Moel for Syria Using the Local Gravimetric an GPS Measurements 1 Ria Al-Masri 2 Abstract "The objective of this paper is to iscuss recent research one by the author to evelop a new improve

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects IOP Conference Series: Materials Science an Engineering 3-D FEM Moeling of fiber/matrix interface eboning in UD composites incluing surface effects To cite this article: A Pupurs an J Varna 2012 IOP Conf.

More information

DAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1

DAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1 11 th National Congress on Theoretical an Applie Mechanics, 2-5 Sept. 2009, Borovets, Bulgaria DAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1 E. MANOACH Institute of Mechanics,

More information

Evaporating droplets tracking by holographic high speed video in turbulent flow

Evaporating droplets tracking by holographic high speed video in turbulent flow Evaporating roplets tracking by holographic high spee vieo in turbulent flow Loïc Méès 1*, Thibaut Tronchin 1, Nathalie Grosjean 1, Jean-Louis Marié 1 an Corinne Fournier 1: Laboratoire e Mécanique es

More information

VUMAT for Fabric Reinforced Composites

VUMAT for Fabric Reinforced Composites VUMAT or Fabric Reinorce Composites. Introuction This ocument escribes a constitutive mo or abric reinorce composites that was introuce in Abaqus/Exicit 6.8. The mo has been imemente as a built-in VUMAT

More information

A Model of Electron-Positron Pair Formation

A Model of Electron-Positron Pair Formation Volume PROGRESS IN PHYSICS January, 8 A Moel of Electron-Positron Pair Formation Bo Lehnert Alfvén Laboratory, Royal Institute of Technology, S-44 Stockholm, Sween E-mail: Bo.Lehnert@ee.kth.se The elementary

More information

Sensors & Transducers 2015 by IFSA Publishing, S. L.

Sensors & Transducers 2015 by IFSA Publishing, S. L. Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity

More information

FINAL EXAM 1 SOLUTIONS Below is the graph of a function f(x). From the graph, read off the value (if any) of the following limits: x 1 +

FINAL EXAM 1 SOLUTIONS Below is the graph of a function f(x). From the graph, read off the value (if any) of the following limits: x 1 + FINAL EXAM 1 SOLUTIONS 2011 1. Below is the graph of a function f(x). From the graph, rea off the value (if any) of the following its: x 1 = 0 f(x) x 1 + = 1 f(x) x 0 = x 0 + = 0 x 1 = 1 1 2 FINAL EXAM

More information

The Application of Euler Lagrange Method of Optimization for Electromechanical Motion Control

The Application of Euler Lagrange Method of Optimization for Electromechanical Motion Control HE ANNAS OF "DUNAREA DE OS" UNIVERSIY OF GAAI FASCICE III, ISSN -454X EECROECHNICS, EECRONICS, AUOMAIC CONRO, INFORMAICS he Application o Euler agrange Metho o Optiization or Electroechanical Motion Control

More information

G j dq i + G j. q i. = a jt. and

G j dq i + G j. q i. = a jt. and Lagrange Multipliers Wenesay, 8 September 011 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine

More information

Problem 1 (20 points)

Problem 1 (20 points) ME 309 Fall 01 Exam 1 Name: C Problem 1 0 points Short answer questions. Each question is worth 5 points. Don t spen too long writing lengthy answers to these questions. Don t use more space than is given.

More information

Study on aero-acoustic structural interactions in fan-ducted system

Study on aero-acoustic structural interactions in fan-ducted system Stuy on aero-acoustic structural interactions in fan-ucte system Yan-kei CHIANG 1 ; Yat-sze CHOY ; Li CHENG 3 ; Shiu-keung TANG 4 1,, 3 Department of Mechanical Engineering, The Hong Kong Polytechnic University,

More information

Objective: To introduce the equations of motion and describe the forces that act upon the Atmosphere

Objective: To introduce the equations of motion and describe the forces that act upon the Atmosphere Objective: To introuce the equations of motion an escribe the forces that act upon the Atmosphere Reaing: Rea pp 18 6 in Chapter 1 of Houghton & Hakim Problems: Work 1.1, 1.8, an 1.9 on pp. 6 & 7 at the

More information

EXPERIMENTAL AND NUMERICAL STUDY ON FREE SPANNING PIPELINES SUBJECTED TO EARTHQUAKES

EXPERIMENTAL AND NUMERICAL STUDY ON FREE SPANNING PIPELINES SUBJECTED TO EARTHQUAKES The th Worl Conerence on Earthquake Enineerin October -7, 8, Beijin, China EXPERIMENTAL AND NUMERICAL STUDY ON FREE SPANNING PIPELINES SUBJECTED TO EARTHQUAKES X. Li, Q. Jin, R.B. Don 3 an J. Zhou Associate

More information

Dusty Plasma Void Dynamics in Unmoving and Moving Flows

Dusty Plasma Void Dynamics in Unmoving and Moving Flows 7 TH EUROPEAN CONFERENCE FOR AERONAUTICS AND SPACE SCIENCES (EUCASS) Dusty Plasma Voi Dynamics in Unmoving an Moving Flows O.V. Kravchenko*, O.A. Azarova**, an T.A. Lapushkina*** *Scientific an Technological

More information

Simple Electromagnetic Motor Model for Torsional Analysis of Variable Speed Drives with an Induction Motor

Simple Electromagnetic Motor Model for Torsional Analysis of Variable Speed Drives with an Induction Motor DOI: 10.24352/UB.OVGU-2017-110 TECHNISCHE MECHANIK, 37, 2-5, (2017), 347-357 submitte: June 15, 2017 Simple Electromagnetic Motor Moel for Torsional Analysis of Variable Spee Drives with an Inuction Motor

More information

Outline. Calculus for the Life Sciences II. Introduction. Tides Introduction. Lecture Notes Differentiation of Trigonometric Functions

Outline. Calculus for the Life Sciences II. Introduction. Tides Introduction. Lecture Notes Differentiation of Trigonometric Functions Calculus for the Life Sciences II c Functions Joseph M. Mahaffy, mahaffy@math.ssu.eu Department of Mathematics an Statistics Dynamical Systems Group Computational Sciences Research Center San Diego State

More information

for which the garbage will fall within the can. If the range is large, the shot is easy; if the range is small, the target is hard to hit.

for which the garbage will fall within the can. If the range is large, the shot is easy; if the range is small, the target is hard to hit. EN40: Dynamics an Vibrations Homework : Soling equations of motion for particles Due Friay Feb 18 011 School of Engineering Brown Uniersity 1. Suppose that you wish to throw a piece of garbage into a garbage

More information

fv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n

fv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n Chapter 11 Rossby waves Supplemental reaing: Pelosky 1 (1979), sections 3.1 3 11.1 Shallow water equations When consiering the general problem of linearize oscillations in a static, arbitrarily stratifie

More information

Modeling and analysis of parallel connected permanent magnet synchronous generators in a small hydropower plant

Modeling and analysis of parallel connected permanent magnet synchronous generators in a small hydropower plant Proceeings of the 2006 IASME/WSEAS International Conference on Energy & Environmental Systems, Chalkia, Greece, May 8-10, 2006 (pp83-88) Moeling an analysis of parallel connecte permanent magnet synchronous

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

Stability-Guaranteed Impedance Control of Hydraulic Robotic Manipulators

Stability-Guaranteed Impedance Control of Hydraulic Robotic Manipulators Tampere University o Technology Stability-Guarantee Impeance Control o Hyraulic Robotic Manipulators Citation Koivumäki, J., & Mattila, J. (217). Stability-Guarantee Impeance Control o Hyraulic Robotic

More information

The proper definition of the added mass for the water entry problem

The proper definition of the added mass for the water entry problem The proper efinition of the ae mass for the water entry problem Leonaro Casetta lecasetta@ig.com.br Celso P. Pesce ceppesce@usp.br LIE&MO lui-structure Interaction an Offshore Mechanics Laboratory Mechanical

More information

ME 141. Engineering Mechanics

ME 141. Engineering Mechanics ME 141 Engineering Mechanics Lecture : Statics of particles Ahma Shahei Shakil Lecturer, Dept. of Mechanical Engg, BUET E-mail: sshakil@me.buet.ac.b, shakil6791@gmail.com Website: teacher.buet.ac.b/sshakil

More information

Goldstein Chapter 1 Exercises

Goldstein Chapter 1 Exercises Golstein Chapter 1 Exercises Michael Goo July 17, 2004 1 Exercises 11. Consier a uniform thin isk that rolls without slipping on a horizontal plane. A horizontal force is applie to the center of the isk

More information

Patterns in bistable resonant-tunneling structures

Patterns in bistable resonant-tunneling structures PHYSICAL REVIEW B VOLUME 56, NUMBER 20 Patterns in bistable resonant-tunneling structures 15 NOVEMBER 1997-II B. A. Glavin an V. A. Kochelap Institute of Semiconuctor Physics, Ukrainian Acaemy of Sciences,

More information

Computed Tomography Notes, Part 1. The equation that governs the image intensity in projection imaging is:

Computed Tomography Notes, Part 1. The equation that governs the image intensity in projection imaging is: Noll 3 CT Notes : Page Compute Tomograph Notes Part Challenges with Projection X-ra Sstems The equation that governs the image intensit in projection imaging is: z I I ep µ z Projection -ra sstems are

More information

Derivation of angular momentum balance law using the cauchy stress tensor measure. (HW#4, MAE 295. UCI)

Derivation of angular momentum balance law using the cauchy stress tensor measure. (HW#4, MAE 295. UCI) Derivation of angular momentum balance law using the cauchy stress tensor measure. (HW#4, MAE 295. UCI) by Nasser Abbasi February 28, 2006 Problem Derive the angular momentum balance (AMB) equation for

More information

State-Space Model for a Multi-Machine System

State-Space Model for a Multi-Machine System State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal

More information

Ordinary Differential Equations: Homework 1

Ordinary Differential Equations: Homework 1 Orinary Differential Equations: Homework 1 M. Gameiro, J.-P. Lessar, J.D. Mireles James, K. Mischaikow January 12, 2017 2 Chapter 1 Motivation 1.1 Exercises Exercise 1.1.1. (Frictionless spring) Consier

More information

Angles-Only Orbit Determination Copyright 2006 Michel Santos Page 1

Angles-Only Orbit Determination Copyright 2006 Michel Santos Page 1 Angles-Only Orbit Determination Copyright 6 Michel Santos Page 1 Abstract This ocument presents a re-erivation of the Gauss an Laplace Angles-Only Methos for Initial Orbit Determination. It keeps close

More information

With the Chain Rule. y x 2 1. and. with respect to second axle. dy du du dx. Rate of change of first axle. with respect to third axle

With the Chain Rule. y x 2 1. and. with respect to second axle. dy du du dx. Rate of change of first axle. with respect to third axle 0 CHAPTER Differentiation Section The Chain Rule Fin the erivative of a composite function using the Chain Rule Fin the erivative of a function using the General Power Rule Simplif the erivative of a function

More information

Solutions for Homework #8. Landing gear

Solutions for Homework #8. Landing gear Solutions or Homewor #8 PROBEM. (P. 9 on page 78 in the note) An airplane is modeled as a beam with masses as shown below: m m m m π [rad/sec] anding gear m m.5 Find the stiness and mass matrices. Find

More information

SOLUTION & ANSWER FOR KCET-2009 VERSION A1 [PHYSICS]

SOLUTION & ANSWER FOR KCET-2009 VERSION A1 [PHYSICS] SOLUTION & ANSWER FOR KCET-009 VERSION A [PHYSICS]. The number of significant figures in the numbers.8000 ---- 5 an 7.8000 5 significant igits 8000.50 7 significant igits. β-ecay means emission of electron

More information

Modeling of vegetated rivers for inbank and overbank ows

Modeling of vegetated rivers for inbank and overbank ows Loughborough University Institutional Repository Moeling o vegetate rivers or inbank an overbank ows This item was submitte to Loughborough University's Institutional Repository by the/an author. Citation:

More information

The effect of nonvertical shear on turbulence in a stably stratified medium

The effect of nonvertical shear on turbulence in a stably stratified medium The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:

More information

The Energy Flow Approach for Oscillation Source Location and Damping Evaluation

The Energy Flow Approach for Oscillation Source Location and Damping Evaluation 1 The Energy Flow Approach for Oscillation Source Location an Damping Evaluation Lei CHEN PhD, Associate Professor Department of Electrical Engineering Tsinghua University chenlei08@tsinghua.eu.cn Backgroun

More information

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate

More information

Applications of First Order Equations

Applications of First Order Equations Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more

More information