DOUBLE PENDULUM VIBRATION MOTION IN FLUID FLOW

Similar documents
Chapter 2 Lagrangian Modeling

International Conference KNOWLEDGE-BASED ORGANIZATION Vol. XXIII No

DOUBLE PENDULUM MOTION ANALYSIS IN VARIABLE FLUID FLOW

The Principle of Least Action

AN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A

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

Numerical Integrator. Graphics

MATHEMATICAL MODELLING OF INVERTED PENDULUM WITH DISTURBANCE INPUT

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

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

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

Gyroscopic matrices of the right beams and the discs

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) }.

2.3 Product and Quotient Rules and Higher-Order Derivatives

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

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

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

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

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

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

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

Lagrangian and Hamiltonian Mechanics

6. Friction and viscosity in gasses

Class Notes for Advanced Dynamics (MEAM535)

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

Chapter 4. Electrostatics of Macroscopic Media

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

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

CIRCULAR MOTION AND SHM: Solutions to Higher Level Questions

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

CENTURION UNIVERSITY OF TECHNOLOGY & MANAGEMENT,ODISHA CUEE-2015

Nonlinear Backstepping Control of Permanent Magnet Synchronous Motor (PMSM)

Tutorial Test 5 2D welding robot

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

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

Vectors in two dimensions

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

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

DYNAMIC PERFORMANCE OF RELUCTANCE SYNCHRONOUS MACHINES

Day 4: Motion Along a Curve Vectors

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

Estimation amount of snow deposits on the road

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

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

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

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

ARCH 614 Note Set 5 S2012abn. Moments & Supports

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

PHYS1169: Tutorial 8 Solutions

Lecture 27: Generalized Coordinates and Lagrange s Equations of Motion

Vibration Analysis of Railway Tracks Forced by Distributed Moving Loads

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

Chapter 14: Vibration Generations Mechanisms: Self Excited Vibration

MATH section 2.3 Basic Differentiation Formulas Page 1 of 5

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

Solutions to the Exercises of Chapter 9

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

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

PERMANENT MAGNETS CHAPTER MAGNETIC POLES AND BAR MAGNETS

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

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

Table of Common Derivatives By David Abraham

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

DAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1

Evaporating droplets tracking by holographic high speed video in turbulent flow

VUMAT for Fabric Reinforced Composites

A Model of Electron-Positron Pair Formation

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

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 +

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

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

Problem 1 (20 points)

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

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

EXPERIMENTAL AND NUMERICAL STUDY ON FREE SPANNING PIPELINES SUBJECTED TO EARTHQUAKES

Dusty Plasma Void Dynamics in Unmoving and Moving Flows

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

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

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.

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

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

Separation of Variables

Stability-Guaranteed Impedance Control of Hydraulic Robotic Manipulators

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

ME 141. Engineering Mechanics

Goldstein Chapter 1 Exercises

Patterns in bistable resonant-tunneling structures

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

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

State-Space Model for a Multi-Machine System

Ordinary Differential Equations: Homework 1

Angles-Only Orbit Determination Copyright 2006 Michel Santos Page 1

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

Solutions for Homework #8. Landing gear

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

Modeling of vegetated rivers for inbank and overbank ows

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

The Energy Flow Approach for Oscillation Source Location and Damping Evaluation

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

Applications of First Order Equations

Transcription:

ENGINEERING FOR RURA DEEOPMENT Jelgava,.-.5.5. DOUBE PENDUUM IBRATION MOTION IN FUID FOW Janis iba, Maris Eiuks, Martins Irbe Riga Technical University, atvia janis.viba@rtu.lv, maris.eiuks@rtu.lv, martins.irbe@rtu.lv 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

ENGINEERING FOR RURA DEEOPMENT Jelgava,.-.5.5. 5 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. ()

ENGINEERING FOR RURA DEEOPMENT Jelgava,.-.5.5. 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

ENGINEERING FOR RURA DEEOPMENT Jelgava,.-.5.5. 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

ENGINEERING FOR RURA DEEOPMENT Jelgava,.-.5.5. 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. 5-7. 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

ENGINEERING FOR RURA DEEOPMENT Jelgava,.-.5.5. 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. 6-65.. Martynyuk A. A. an Nikitina N.. The Theory o Motion o a Double Mathematical Penulum, International Applie Mechanics, ol. 6, No. 9,, pp. 5-58.. Double Penulum rom Eric Weisstein s. Worl o Physics scienceworl.wolram.com Mechanics Penula [online][..] Available at: http://scienceworl.wolram.com/ physics/doublepenulum.html. Double Penulum. [online][..] Available at: http://www.myphysicslab.com/bl_ penulum.html 5. Arnol,. I. Problem in Mathematical Methos o lassical Mechanics, n e. New York: Springer-erlag, p. 9, 989. 6. Wells, D. A. Theory an Problems o agrangian Dynamics. New York: McGraw-Hill, pp. -,, an -, 967. 7. arry M. Bates; James M. Nester. On D Alembert s Principle. ommunications in Mathematics (). olume: 9, Issue:, page 57-7. ISSN: 8-88. 9