Dynamics of Non-Smooth Systems

Size: px
Start display at page:

Download "Dynamics of Non-Smooth Systems"

Transcription

1 Dynamics of Non-Smooth Systems P Chandramouli Indian Institute of Technology Madras October 5, 2013 Mouli, IIT Madras Non-Smooth System Dynamics 1/39

2 Introduction Filippov Systems Behaviour around hyper-surface Numerical Implementation Case Studies Mouli, IIT Madras Non-Smooth System Dynamics 2/39

3 What is a non-smooth system Physical systems often operate in different modes Transition from one mode to another idealized as instantaneous Time scales of transition usually small Examples Dry friction Switches Switch from slip to stick or vice-versa Systems with impact/backlash Switch from zero to finite stiffness Diodes in electrical circuits Switch from short-circuit when forward biased to open-circuit when reverse biased Mouli, IIT Madras Non-Smooth System Dynamics 3/39

4 Continuous and Smooth Functions A function f (x) is continuous at x X if for all ɛ > 0 there exists δ > 0 such that for y B(x; δ) implies f (x) f (y) ɛ Roughly speaking this means for a single-valued function one can draw the function on a piece of paper without taking our hands off the paper A function f (x) is smooth if it is continuously differentible upto any order in x Mouli, IIT Madras Non-Smooth System Dynamics 4/39

5 Non-smooth Type Continuous vector field but non-smooth mẍ + cẋ + kx = f 0 cos(ωt) f (x) { f (x) = 0, if x < x c k c (x x c ), else Support Derivatives of f (x) discontinuous and hence f (x) non-smooth Mouli, IIT Madras Non-Smooth System Dynamics 5/39

6 Discontinuous Continuous vector field but non-smooth mẍ + kx = f 0 cos(ωt) f (ẋ) µ, ẋ > 0 f (ẋ) [ µ, µ], ẋ = 0 µ, ẋ < 0 Friction f (ẋ) discontinuous Called Filippov type or differential inclusion Mouli, IIT Madras Non-Smooth System Dynamics 6/39

7 Aircraft Engine Vibration Imbalance forces due to fan blade loss Vibration of engine caused Blades start to rub against stator Contact with friction occurs Will the engine survive for rest of flight? Simulation needs to provide the answer Mouli, IIT Madras Non-Smooth System Dynamics 7/39

8 Example Let us look at an example ẋ = f (x) = 1 2 sgn(x) For initial conditions x(0) 0 one can obtain a solution { 3t + C 1, x < 0 x(t) = t + C 2, x > 0 C 1 and C 2 are determined from IC Both these solutions reach x = 0 in finite time If it arrives at x = 0 it cannot leave x = 0 Why? PhasePlane Mouli, IIT Madras Non-Smooth System Dynamics 8/39

9 Example System Dynamics When x < 0 the velocity ẋ > 0 However if x > 0 the velocity ẋ < 0 Hence when the system trajectory reaches x = 0 it implies that ẋ = 0 There is no discontinuity of the variables as a function of time Therefore the system continues to remain at x = 0 This does not agree with our usual notion that sgn(0) = 0 Then ẋ(0) = 1 0 Mouli, IIT Madras Non-Smooth System Dynamics 9/39

10 Set Valued Functions We can deal with the situation if we replace the function concept with a set value function { 1}, x < 0 Sgn(x) = [ 1, 1], x = 0 {1}, x > 0 In the above set [a, b] implies the interval {x R a x b} The set {a, b} implies the elements are a and b The differential equation is then replaced by the differential inclusion ẋ F(x) or ẋ 1 2 Sgn(x) Mouli, IIT Madras Non-Smooth System Dynamics 10/39

11 Solution With the use of set-valued functions x(t) = 0 is a unique global solution to the differential inclusion ẋ 1 2 Sgn(x) With initial condition x(0) = 0 The example used is a one-dimensional one Let us now try to generalize to a multi-dimensional case Mouli, IIT Madras Non-Smooth System Dynamics 11/39

12 Multi-dimensional Case We assume that the discontinuous right hand side is limited to one hyper-surface The space R n is split into two sub-spaces V and V + by a hyper-surface S R n = V S V + The hyper-surface S is defined by a scalar indicator function h(x(t)) When x(t) is on S then h(x) = 0 The normal to the hypersurface is given by n(x) = h(x) Mouli, IIT Madras Non-Smooth System Dynamics 12/39

13 Different Spaces The different spaces can now be written in terms of the indicator function V = {x R n h(x(t)) < 0} S = {x R n h(x(t)) = 0} V + = {x R n h(x(t)) > 0} The vector field f(x, t) is assumed to be smooth and locally continuous provided x / S { f (x, t), x V ẋ(t) = f(x, t) = f + (x, t), x V + The above system of differential equations does not define the dynamics if x S Mouli, IIT Madras Non-Smooth System Dynamics 13/39

14 Multi-dimensional Inclusion We do a set-value extension like the single dimension case f (x, t), x V ẋ(t) F(x, t) = co(f, f + ), x S f + (x, t), x V + The co represents a convex set which means co(f, f +) = {qf + (1 q)f +, q [0, 1]} The solution of Filippov systems are not unique Two different types of behaviour around S will be demonstrated Mouli, IIT Madras Non-Smooth System Dynamics 14/39

15 Example 2 Let us look at the following system x 1 = sgn(x 2 c) x 2 = sgn(x 2 c) The phase plane for this case is shown here In this case the following holds (n T f ) (n T f + ) > 0 for x S Transversal Mouli, IIT Madras Non-Smooth System Dynamics 15/39

16 Example 2 Let us look at the following system x 1 = sgn(x 2 c) x 2 = 2 4 sgn(x 2 c) The phase plane for this case is shown here Aslide In this case the solution is attracted to the hyper-surface and move along the surface It is called a sliding mode This implies n T f = 0 One can express f = αf + + (1 α)f Calculate α = nt f n T (f f + ) Mouli, IIT Madras Non-Smooth System Dynamics 16/39

17 Final remarks Solution to the differential inclusion with x 0 S is locally unique in forward time if Projections of the vector fields are on the same side of S (n T f ) (n T f +) > 0 Projections point to S (n T f ) > 0 and (n T f +) < 0 Mouli, IIT Madras Non-Smooth System Dynamics 17/39

18 Numerical Implementation: Switch Model One could use smoothening of the sgn(x) function Arctan approximation 2 π arctan γx with γ quite large Often leads to stiff differential equations Use of switch model helps solve these systems without encountering stiff equations A tolerance band near the indicator function hyper-surface is made Band Mouli, IIT Madras Non-Smooth System Dynamics 18/39

19 Procedure Check if h(x) < tol If it is so then we are in the band Now compute the projections of the vector fields n T f and n T f + If n T f > 0 and n T f + < 0 we have an attractive sliding mode The opposite of this is the repulsive sliding mode which I have not discussed To ensure convergence from the edge of the band to the middle for the sliding mode we let ḣ = 1 τ h Mouli, IIT Madras Non-Smooth System Dynamics 19/39

20 Calculations Now since ḣ = h x x t one can simplify as ḣ = n T (αf + + (1 α)f ) From this one can solve for α as α = nt f + h τ n T (f f +) Choose the parameter tol sufficiently small If the projections are of the same sign then simply continue integrating Mouli, IIT Madras Non-Smooth System Dynamics 20/39

21 Other Methods Switch model is not suitable if many friction interfaces or contact/friction interfaces are present Event driven integration methods are used Contact/friction is solved using Linear Complementarity Programming (LCP) approach Solution gives the next mode of the trajectory of the system Mouli, IIT Madras Non-Smooth System Dynamics 21/39

22 3 DOF Brake Model We now look at an example problem to apply the Filippov Method Mouli, IIT Madras Non-Smooth System Dynamics 22/39

23 Result velocity displacement Mouli, IIT Madras Non-Smooth System Dynamics 23/39

24 SD Oscillator Realization of the Oscillator Non-dimensional equation is x + 2ζx + x(1 1 x 2 +α 2 ) = f 0 cos Ωτ RestForce Mouli, IIT Madras Non-Smooth System Dynamics 24/39

25 Typical Result Results obtained for α = 0 using switch model and event based time integration Mouli, IIT Madras Non-Smooth System Dynamics 25/39

26 Friction Interface of Turbine Blade Variable normal load model Four states: Stick, two slip states and separation possible mü + c u + ku + f u = F u (t); m v + f v = F v (t) Mouli, IIT Madras Non-Smooth System Dynamics 26/39

27 Nonlinear Forces k u (u u 0 ) + f 0, stick µ(n 0 + k v v), positive slip f u = µ(n 0 + k v v), negative slip 0, separation { n0 + k v v, contact f v = 0, separation Mouli, IIT Madras Non-Smooth System Dynamics 27/39

28 Four States Mouli, IIT Madras Non-Smooth System Dynamics 28/39

29 Response with different Pre-loads Mouli, IIT Madras Non-Smooth System Dynamics 29/39

30 Concluding Remarks I have not talked about stability of periodic solutions here This is done from the fundamental solution matrix Φ(x, t) The eigenvalues of this matrix are the Floquet multipliers More recent work available using shooting methods for Filippov systems We are now looking at using a time-variational method for such systems Faster computations as we seek periodic solutions in time domain Mouli, IIT Madras Non-Smooth System Dynamics 30/39

31 Major References Leine, R. I. and Nijmeijer, H., 2004, Dynamics and Bifurcations of Non-Smooth Mechanical Systems, Springer-Verlag: Berlin Moreau, J. J. and Panagiotopoulos, P. D., 1988, Non-smooth Mechanics and Applications, CISM Courses and Lectures, vol. 302, Springer: Wien. Glocker, Ch., 2001, Set-Valued Force Laws, Dynamics of Non-smooth Systems, Lecture Notes in Applied Mechanics, vol. 1, Springer-Verlag:Berlin. Mouli, IIT Madras Non-Smooth System Dynamics 31/39

32 Switch Types Return Mouli, IIT Madras Non-Smooth System Dynamics 32/39

33 Support with Clearance Back Mouli, IIT Madras Non-Smooth System Dynamics 33/39

34 Support with Clearance Return Mouli, IIT Madras Non-Smooth System Dynamics 34/39

35 Phase Plane Trajectory Back Mouli, IIT Madras Non-Smooth System Dynamics 35/39

36 Transversal Intersection Return Mouli, IIT Madras Non-Smooth System Dynamics 36/39

37 Attractive Sliding Mode Back Mouli, IIT Madras Non-Smooth System Dynamics 37/39

38 Switch Model Band Return Mouli, IIT Madras Non-Smooth System Dynamics 38/39

39 Restoring Force Back Mouli, IIT Madras Non-Smooth System Dynamics 39/39

Sliding Modes in Filippov Systems.

Sliding Modes in Filippov Systems. in Filippov Systems. School of Mathematics, Georgia Tech Joint work with Luciano Lopez (Univ. Bari) Geneva Conference In honor of Ernst Hairer 60th birthday References Filippov [DEs with discontinuous

More information

Nonsmooth systems: synchronization, sliding and other open problems

Nonsmooth systems: synchronization, sliding and other open problems John Hogan Bristol Centre for Applied Nonlinear Mathematics, University of Bristol, England Nonsmooth systems: synchronization, sliding and other open problems 2 Nonsmooth Systems 3 What is a nonsmooth

More information

AXIAL STICK-SLIP LIMIT CYCLING IN DRILL-STRING DYNAMICS WITH DELAY

AXIAL STICK-SLIP LIMIT CYCLING IN DRILL-STRING DYNAMICS WITH DELAY ENOC-2005, Eindhoven, Netherlands, 7-12 August 2005 AXIAL STICK-SLIP LIMIT CYCLING IN DRILL-STRING DYNAMICS WITH DELAY Christophe Germay Department of Electrical Engineering and Computer Science University

More information

Controlling Attractivity of Friction-Induced Equilibrium Sets

Controlling Attractivity of Friction-Induced Equilibrium Sets Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 2005 Seville, Spain, December 12-15, 2005 TuA16.5 Controlling Attractivity of Friction-Induced Equilibrium

More information

1-DOF Vibration Characteristics. MCE371: Vibrations. Prof. Richter. Department of Mechanical Engineering. Handout 7 Fall 2017

1-DOF Vibration Characteristics. MCE371: Vibrations. Prof. Richter. Department of Mechanical Engineering. Handout 7 Fall 2017 MCE371: Vibrations Prof. Richter Department of Mechanical Engineering Handout 7 Fall 2017 Free Undamped Vibration Follow Palm, Sect. 3.2, 3.3 (pp 120-138), 3.5 (pp 144-151), 3.8 (pp. 167-169) The equation

More information

8 Example 1: The van der Pol oscillator (Strogatz Chapter 7)

8 Example 1: The van der Pol oscillator (Strogatz Chapter 7) 8 Example 1: The van der Pol oscillator (Strogatz Chapter 7) So far we have seen some different possibilities of what can happen in two-dimensional systems (local and global attractors and bifurcations)

More information

Derivation of border-collision maps from limit cycle bifurcations

Derivation of border-collision maps from limit cycle bifurcations Derivation of border-collision maps from limit cycle bifurcations Alan Champneys Department of Engineering Mathematics, University of Bristol Mario di Bernardo, Chris Budd, Piotr Kowalczyk Gabor Licsko,...

More information

C.-H. Lamarque. University of Lyon/ENTPE/LGCB & LTDS UMR CNRS 5513

C.-H. Lamarque. University of Lyon/ENTPE/LGCB & LTDS UMR CNRS 5513 Nonlinear Dynamics of Smooth and Non-Smooth Systems with Application to Passive Controls 3rd Sperlonga Summer School on Mechanics and Engineering Sciences on Dynamics, Stability and Control of Flexible

More information

Lehrstuhl B für Mechanik Technische Universität München D Garching Germany

Lehrstuhl B für Mechanik Technische Universität München D Garching Germany DISPLACEMENT POTENTIALS IN NON-SMOOTH DYNAMICS CH. GLOCKER Lehrstuhl B für Mechanik Technische Universität München D-85747 Garching Germany Abstract. The paper treats the evaluation of the accelerations

More information

LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)

LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,

More information

LMI Methods in Optimal and Robust Control

LMI Methods in Optimal and Robust Control LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear

More information

The reduction to invariant cones for nonsmooth systems. Tassilo Küpper Mathematical Institute University of Cologne

The reduction to invariant cones for nonsmooth systems. Tassilo Küpper Mathematical Institute University of Cologne The reduction to invariant cones for nonsmooth systems Tassilo Küpper Mathematical Institute University of Cologne kuepper@math.uni-koeln.de International Workshop on Resonance Oscillations and Stability

More information

Lecture 9 Nonlinear Control Design

Lecture 9 Nonlinear Control Design Lecture 9 Nonlinear Control Design Exact-linearization Lyapunov-based design Lab 2 Adaptive control Sliding modes control Literature: [Khalil, ch.s 13, 14.1,14.2] and [Glad-Ljung,ch.17] Course Outline

More information

ME8230 Nonlinear Dynamics

ME8230 Nonlinear Dynamics ME8230 Nonlinear Dynamics Lecture 1, part 1 Introduction, some basic math background, and some random examples Prof. Manoj Srinivasan Mechanical and Aerospace Engineering srinivasan.88@osu.edu Spring mass

More information

Simulating Two-Dimensional Stick-Slip Motion of a Rigid Body using a New Friction Model

Simulating Two-Dimensional Stick-Slip Motion of a Rigid Body using a New Friction Model Proceedings of the 2 nd World Congress on Mechanical, Chemical, and Material Engineering (MCM'16) Budapest, Hungary August 22 23, 2016 Paper No. ICMIE 116 DOI: 10.11159/icmie16.116 Simulating Two-Dimensional

More information

Complex Dynamic Systems: Qualitative vs Quantitative analysis

Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic

More information

Advanced Dynamics. - Lecture 4 Lagrange Equations. Paolo Tiso Spring Semester 2017 ETH Zürich

Advanced Dynamics. - Lecture 4 Lagrange Equations. Paolo Tiso Spring Semester 2017 ETH Zürich Advanced Dynamics - Lecture 4 Lagrange Equations Paolo Tiso Spring Semester 2017 ETH Zürich LECTURE OBJECTIVES 1. Derive the Lagrange equations of a system of particles; 2. Show that the equation of motion

More information

System Control Engineering 0

System Control Engineering 0 System Control Engineering 0 Koichi Hashimoto Graduate School of Information Sciences Text: Nonlinear Control Systems Analysis and Design, Wiley Author: Horacio J. Marquez Web: http://www.ic.is.tohoku.ac.jp/~koichi/system_control/

More information

A Remark on IVP and TVP Non-Smooth Viscosity Solutions to Hamilton-Jacobi Equations

A Remark on IVP and TVP Non-Smooth Viscosity Solutions to Hamilton-Jacobi Equations 2005 American Control Conference June 8-10, 2005. Portland, OR, USA WeB10.3 A Remark on IVP and TVP Non-Smooth Viscosity Solutions to Hamilton-Jacobi Equations Arik Melikyan, Andrei Akhmetzhanov and Naira

More information

LIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS

LIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS 5 th GRACM International Congress on Computational Mechanics Limassol, 29 June 1 July, 2005 LIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS G.A. Drosopoulos I, G.E.

More information

Implicit Euler Numerical Simulation of Sliding Mode Systems

Implicit Euler Numerical Simulation of Sliding Mode Systems Implicit Euler Numerical Simulation of Sliding Mode Systems Vincent Acary and Bernard Brogliato, INRIA Grenoble Rhône-Alpes Conference in honour of E. Hairer s 6th birthday Genève, 17-2 June 29 Objectives

More information

IMES - Center of Mechanics ETH Zentrum, Tannenstrasse 3 CH-8092 Zürich, Switzerland

IMES - Center of Mechanics ETH Zentrum, Tannenstrasse 3 CH-8092 Zürich, Switzerland Chapter 16 THE GEOMETRY OF NEWTON S CRADLE Christoph Glocker and Ueli Aeberhard IMES - Center of Mechanics ETH Zentrum, Tannenstrasse 3 CH-8092 Zürich, Switzerland glocker@imes.mavt.ethz.ch Abstract A

More information

BIBO STABILITY AND ASYMPTOTIC STABILITY

BIBO STABILITY AND ASYMPTOTIC STABILITY BIBO STABILITY AND ASYMPTOTIC STABILITY FRANCESCO NORI Abstract. In this report with discuss the concepts of bounded-input boundedoutput stability (BIBO) and of Lyapunov stability. Examples are given to

More information

Vibration Control Prof. Dr. S. P. Harsha Department of Mechanical & Industrial Engineering Indian Institute of Technology, Roorkee

Vibration Control Prof. Dr. S. P. Harsha Department of Mechanical & Industrial Engineering Indian Institute of Technology, Roorkee Vibration Control Prof. Dr. S. P. Harsha Department of Mechanical & Industrial Engineering Indian Institute of Technology, Roorkee Module - 1 Review of Basics of Mechanical Vibrations Lecture - 2 Introduction

More information

Thursday, August 4, 2011

Thursday, August 4, 2011 Chapter 16 Thursday, August 4, 2011 16.1 Springs in Motion: Hooke s Law and the Second-Order ODE We have seen alrealdy that differential equations are powerful tools for understanding mechanics and electro-magnetism.

More information

3D revolute joint with clearance in multibody systems

3D revolute joint with clearance in multibody systems 3D revolute joint with clearance in multibody systems Narendra Akhadkar 1, Vincent Acary 2, and Bernard Brogliato 2 1 Schneider Electric, 37 Quai Paul-Louis Merlin, 38000 Grenoble, France, e-mail: narendra.akhadkar@schneider-electric.com

More information

Introductory Physics. Week 2015/05/29

Introductory Physics. Week 2015/05/29 2015/05/29 Part I Summary of week 6 Summary of week 6 We studied the motion of a projectile under uniform gravity, and constrained rectilinear motion, introducing the concept of constraint force. Then

More information

CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD

CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD Vietnam Journal of Mechanics, VAST, Vol. 34, No. 3 (2012), pp. 157 167 CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD Nguyen Van Khang, Hoang Manh Cuong, Nguyen

More information

Figure 1: Schematic of ship in still water showing the action of bouyancy and weight to right the ship.

Figure 1: Schematic of ship in still water showing the action of bouyancy and weight to right the ship. MULTI-DIMENSIONAL SYSTEM: In this computer simulation we will explore a nonlinear multi-dimensional system. As before these systems are governed by equations of the form x 1 = f 1 x 2 = f 2.. x n = f n

More information

Explore Kapitza s Pendulum Behavior via Trajectory Optimization. Yifan Hou

Explore Kapitza s Pendulum Behavior via Trajectory Optimization. Yifan Hou 1 Introduction 12 course Explore Kapitza s Pendulum Behavior via Trajectory Optimization Project for: Mechanics of Manipulation, 16 741 Yifan Hou Andrew id: yifanh houyf11@gmail.com or yifanh@andrew.cmu.edu

More information

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Non-Linear Dynamics Part I

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Non-Linear Dynamics Part I The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Non-Linear Dynamics Part I Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 5/Part A - 23 November,

More information

The problem of singularity in impacting systems

The problem of singularity in impacting systems The problem of singularity in impacting systems Soumitro Banerjee, Department of Physics Indian Institute of Science Education & Research, Kolkata, India The problem of singularity in impacting systems

More information

Dynamics of Ocean Structures Prof. Dr. Srinivasan Chandrasekaran Department of Ocean Engineering Indian Institute of Technology, Madras

Dynamics of Ocean Structures Prof. Dr. Srinivasan Chandrasekaran Department of Ocean Engineering Indian Institute of Technology, Madras Dynamics of Ocean Structures Prof. Dr. Srinivasan Chandrasekaran Department of Ocean Engineering Indian Institute of Technology, Madras Module - 01 Lecture - 09 Characteristics of Single Degree - of -

More information

Exercises Lecture 15

Exercises Lecture 15 AM1 Mathematical Analysis 1 Oct. 011 Feb. 01 Date: January 7 Exercises Lecture 15 Harmonic Oscillators In classical mechanics, a harmonic oscillator is a system that, when displaced from its equilibrium

More information

1. Consider the 1-DOF system described by the equation of motion, 4ẍ+20ẋ+25x = f.

1. Consider the 1-DOF system described by the equation of motion, 4ẍ+20ẋ+25x = f. Introduction to Robotics (CS3A) Homework #6 Solution (Winter 7/8). Consider the -DOF system described by the equation of motion, ẍ+ẋ+5x = f. (a) Find the natural frequency ω n and the natural damping ratio

More information

Hydro-elastic Wagner impact using variational inequalities

Hydro-elastic Wagner impact using variational inequalities Hydro-elastic Wagner impact using variational inequalities Thomas GAZZOLA, Alexander KOROBKIN, Šime MALENICA Introduction A simple model of water impact has been introduced by Wagner [6]. This model is

More information

Review: control, feedback, etc. Today s topic: state-space models of systems; linearization

Review: control, feedback, etc. Today s topic: state-space models of systems; linearization Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered

More information

1 2 predators competing for 1 prey

1 2 predators competing for 1 prey 1 2 predators competing for 1 prey I consider here the equations for two predator species competing for 1 prey species The equations of the system are H (t) = rh(1 H K ) a1hp1 1+a a 2HP 2 1T 1H 1 + a 2

More information

Nonlinear effects on the rotor driven by a motor with limited power

Nonlinear effects on the rotor driven by a motor with limited power Applied and Computational Mechanics 1 (007) 603-61 Nonlinear effects on the rotor driven by a motor with limited power L. Pst Institute of Thermomechanics, Academy of Sciences of CR, Dolejškova 5,18 00

More information

5 Linear Algebra and Inverse Problem

5 Linear Algebra and Inverse Problem 5 Linear Algebra and Inverse Problem 5.1 Introduction Direct problem ( Forward problem) is to find field quantities satisfying Governing equations, Boundary conditions, Initial conditions. The direct problem

More information

PDE and Boundary-Value Problems Winter Term 2014/2015

PDE and Boundary-Value Problems Winter Term 2014/2015 PDE and Boundary-Value Problems Winter Term 2014/2015 Lecture 13 Saarland University 5. Januar 2015 c Daria Apushkinskaya (UdS) PDE and BVP lecture 13 5. Januar 2015 1 / 35 Purpose of Lesson To interpretate

More information

The dynamics of a nonautonomous oscillator with friction memory

The dynamics of a nonautonomous oscillator with friction memory The dynamics of a nonautonomous oscillator with friction memory L. A. Igumnov 1, V. S. Metrikin 2, M. V. Zaytzev 3 1, 2 Research Institute for mechanics, National Research Lobachevsky State University

More information

Damped Oscillation Solution

Damped Oscillation Solution Lecture 19 (Chapter 7): Energy Damping, s 1 OverDamped Oscillation Solution Damped Oscillation Solution The last case has β 2 ω 2 0 > 0. In this case we define another real frequency ω 2 = β 2 ω 2 0. In

More information

Differential Equations

Differential Equations Differential Equations A differential equation (DE) is an equation which involves an unknown function f (x) as well as some of its derivatives. To solve a differential equation means to find the unknown

More information

Problem Set Number 02, j/2.036j MIT (Fall 2018)

Problem Set Number 02, j/2.036j MIT (Fall 2018) Problem Set Number 0, 18.385j/.036j MIT (Fall 018) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139) September 6, 018 Due October 4, 018. Turn it in (by 3PM) at the Math. Problem Set

More information

Poincaré Map, Floquet Theory, and Stability of Periodic Orbits

Poincaré Map, Floquet Theory, and Stability of Periodic Orbits Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct

More information

(x) = lim. dx = f f(x + h) f(x)

(x) = lim. dx = f f(x + h) f(x) Chapter 4 The Derivative In our investigation so far, we have defined the notion of an instantaneous rate of change, and called this the derivative. We have also identified this mathematical concept with

More information

Nonlinear Control Lecture 1: Introduction

Nonlinear Control Lecture 1: Introduction Nonlinear Control Lecture 1: Introduction Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2011 Farzaneh Abdollahi Nonlinear Control Lecture 1 1/15 Motivation

More information

Lecture - 11 Bendixson and Poincare Bendixson Criteria Van der Pol Oscillator

Lecture - 11 Bendixson and Poincare Bendixson Criteria Van der Pol Oscillator Nonlinear Dynamical Systems Prof. Madhu. N. Belur and Prof. Harish. K. Pillai Department of Electrical Engineering Indian Institute of Technology, Bombay Lecture - 11 Bendixson and Poincare Bendixson Criteria

More information

A Geometric Interpretation of Newtonian Impacts with Global Dissipation Index

A Geometric Interpretation of Newtonian Impacts with Global Dissipation Index A Geometric Interpretation of Newtonian Impacts with Global Dissipation Index Christoph Glocker Institute B of Mechanics, Technical University of Munich, D-85747 Garching, Germany (glocker@lbm.mw.tu-muenchen.de)

More information

Szalai, R., & Osinga, H. M. (2007). Unstable manifolds of a limit cycle near grazing.

Szalai, R., & Osinga, H. M. (2007). Unstable manifolds of a limit cycle near grazing. Szalai, R., & Osinga, H. M. (2007). Unstable manifolds of a limit cycle near grazing. Early version, also known as pre-print Link to publication record in Explore Bristol Research PDF-document University

More information

Dynamics of Ocean Structures Prof. Dr. Srinivasan Chandrasekaran Department of Ocean Engineering Indian Institute of Technology, Madras

Dynamics of Ocean Structures Prof. Dr. Srinivasan Chandrasekaran Department of Ocean Engineering Indian Institute of Technology, Madras Dynamics of Ocean Structures Prof. Dr. Srinivasan Chandrasekaran Department of Ocean Engineering Indian Institute of Technology, Madras Module - 1 Lecture - 13 Undamped and Damped Systems II (Refer Slide

More information

Some Aspects of Structural Dynamics

Some Aspects of Structural Dynamics Appendix B Some Aspects of Structural Dynamics This Appendix deals with some aspects of the dynamic behavior of SDOF and MDOF. It starts with the formulation of the equation of motion of SDOF systems.

More information

(Refer Slide Time: 00:32)

(Refer Slide Time: 00:32) Nonlinear Dynamical Systems Prof. Madhu. N. Belur and Prof. Harish. K. Pillai Department of Electrical Engineering Indian Institute of Technology, Bombay Lecture - 12 Scilab simulation of Lotka Volterra

More information

Simulation of the Stick-Slip Friction between Steering Shafts Using ADAMS/PRE

Simulation of the Stick-Slip Friction between Steering Shafts Using ADAMS/PRE Simulation of the Stick-Slip Friction between Steering Shafts Using ADAMS/PRE Dexin Wang and Yuting Rui Research & Vehicle Technology Ford Motor Company ABSTRACT Cyclic stick-slip friction is a well-known

More information

AME 301: Differential Equations, Control and Vibrations

AME 301: Differential Equations, Control and Vibrations AME 301: Differential Equations, Control and Vibrations Introduction, Overview and Motivation Some boring catalog stuff: First of a two-course sequence, the course introduces methods of differential-equation

More information

Lecture 1 Notes: 06 / 27. The first part of this class will primarily cover oscillating systems (harmonic oscillators and waves).

Lecture 1 Notes: 06 / 27. The first part of this class will primarily cover oscillating systems (harmonic oscillators and waves). Lecture 1 Notes: 06 / 27 The first part of this class will primarily cover oscillating systems (harmonic oscillators and waves). These systems are very common in nature - a system displaced from equilibrium

More information

CDS 101/110a: Lecture 2.1 Dynamic Behavior

CDS 101/110a: Lecture 2.1 Dynamic Behavior CDS 11/11a: Lecture.1 Dynamic Behavior Richard M. Murray 6 October 8 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium

More information

Lecture 3 Dynamics 29

Lecture 3 Dynamics 29 Lecture 3 Dynamics 29 30 LECTURE 3. DYNAMICS 3.1 Introduction Having described the states and the observables of a quantum system, we shall now introduce the rules that determine their time evolution.

More information

The Stick-Slip Vibration and Bifurcation of a Vibro-Impact System with Dry Friction

The Stick-Slip Vibration and Bifurcation of a Vibro-Impact System with Dry Friction Send Orders for Reprints to reprints@benthamscience.ae 308 The Open Mechanical Engineering Journal, 2014, 8, 308-313 Open Access The Stick-Slip Vibration and Bifurcation of a Vibro-Impact System with Dry

More information

Dynamics of Structures

Dynamics of Structures Dynamics of Structures Elements of structural dynamics Roberto Tomasi 11.05.2017 Roberto Tomasi Dynamics of Structures 11.05.2017 1 / 22 Overview 1 SDOF system SDOF system Equation of motion Response spectrum

More information

Handout 2: Invariant Sets and Stability

Handout 2: Invariant Sets and Stability Engineering Tripos Part IIB Nonlinear Systems and Control Module 4F2 1 Invariant Sets Handout 2: Invariant Sets and Stability Consider again the autonomous dynamical system ẋ = f(x), x() = x (1) with state

More information

On Stability of Linear Complementarity Systems

On Stability of Linear Complementarity Systems On Stability of Linear Complementarity Systems Alen Agheksanterian University of Maryland Baltimore County Math 710A: Special Topics in Optimization On Stability of Linear Complementarity Systems Alen

More information

Effect of various periodic forces on Duffing oscillator

Effect of various periodic forces on Duffing oscillator PRAMANA c Indian Academy of Sciences Vol. 67, No. 2 journal of August 2006 physics pp. 351 356 Effect of various periodic forces on Duffing oscillator V RAVICHANDRAN 1, V CHINNATHAMBI 1, and S RAJASEKAR

More information

Computers and Mathematics with Applications. Chaos suppression via periodic pulses in a class of piece-wise continuous systems

Computers and Mathematics with Applications. Chaos suppression via periodic pulses in a class of piece-wise continuous systems Computers and Mathematics with Applications 64 (2012) 849 855 Contents lists available at SciVerse ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa

More information

CDS 101/110a: Lecture 2.1 Dynamic Behavior

CDS 101/110a: Lecture 2.1 Dynamic Behavior CDS 11/11a: Lecture 2.1 Dynamic Behavior Richard M. Murray 6 October 28 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium

More information

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016 Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural

More information

Converse Lyapunov theorem and Input-to-State Stability

Converse Lyapunov theorem and Input-to-State Stability Converse Lyapunov theorem and Input-to-State Stability April 6, 2014 1 Converse Lyapunov theorem In the previous lecture, we have discussed few examples of nonlinear control systems and stability concepts

More information

Stability analysis of sliding-grazing phenomenon in dry-friction oscillator using TS fuzzy approach

Stability analysis of sliding-grazing phenomenon in dry-friction oscillator using TS fuzzy approach Stability analysis of sliding-grazing phenomenon in dry-friction oscillator using TS fuzzy approach Kamyar Mehran School of Engineering University of Warwick Coventry, CV4 7LW, UK Email: k.mehran@warwick.ac.uk

More information

SCALAR FORCE POTENTIALS IN RIGID MULTIBODY SYSTEMS. Ch. Glocker Technical University of Munich, Germany

SCALAR FORCE POTENTIALS IN RIGID MULTIBODY SYSTEMS. Ch. Glocker Technical University of Munich, Germany SCALAR FORCE POTENTIALS IN RIGID MULTIBODY SYSTEMS Ch. Glocker Technical University of Munich, Germany ABSTRACT These lectures treat the motion of finite-dimensional mechanical systems under the influence

More information

Kinetics of Particles: Work and Energy

Kinetics of Particles: Work and Energy Kinetics of Particles: Work and Energy Total work done is given by: Modifying this eqn to account for the potential energy terms: U 1-2 + (-ΔV g ) + (-ΔV e ) = ΔT T U 1-2 is work of all external forces

More information

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution

More information

Continuous Threshold Policy Harvesting in Predator-Prey Models

Continuous Threshold Policy Harvesting in Predator-Prey Models Continuous Threshold Policy Harvesting in Predator-Prey Models Jon Bohn and Kaitlin Speer Department of Mathematics, University of Wisconsin - Madison Department of Mathematics, Baylor University July

More information

Nonlinear Control. Nonlinear Control Lecture # 2 Stability of Equilibrium Points

Nonlinear Control. Nonlinear Control Lecture # 2 Stability of Equilibrium Points Nonlinear Control Lecture # 2 Stability of Equilibrium Points Basic Concepts ẋ = f(x) f is locally Lipschitz over a domain D R n Suppose x D is an equilibrium point; that is, f( x) = 0 Characterize and

More information

Physics 326 Lab 6 10/18/04 DAMPED SIMPLE HARMONIC MOTION

Physics 326 Lab 6 10/18/04 DAMPED SIMPLE HARMONIC MOTION DAMPED SIMPLE HARMONIC MOTION PURPOSE To understand the relationships between force, acceleration, velocity, position, and period of a mass undergoing simple harmonic motion and to determine the effect

More information

DETC2001/VIB NONLINEAR DYNAMICS OF THE WOODPECKER TOY

DETC2001/VIB NONLINEAR DYNAMICS OF THE WOODPECKER TOY Proceedings of DETC' ASME Design Engineering Technical Conference and Computers and Information in Engineering Conference Pittsburgh, PA, September 9-, DETC/VIB-68 DETC/VIB-68 NONLINEAR DYNAMICS OF THE

More information

AP Physics C Mechanics Objectives

AP Physics C Mechanics Objectives AP Physics C Mechanics Objectives I. KINEMATICS A. Motion in One Dimension 1. The relationships among position, velocity and acceleration a. Given a graph of position vs. time, identify or sketch a graph

More information

Theory and Practice of Rotor Dynamics Prof. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati

Theory and Practice of Rotor Dynamics Prof. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati Theory and Practice of Rotor Dynamics Prof. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati Module - 7 Instability in Rotor Systems Lecture - 2 Fluid-Film Bearings

More information

Dynamics of structures

Dynamics of structures Dynamics of structures 1.2 Viscous damping Luc St-Pierre October 30, 2017 1 / 22 Summary so far We analysed the spring-mass system and found that its motion is governed by: mẍ(t) + kx(t) = 0 k y m x x

More information

Nonlinear dynamics of dry friction

Nonlinear dynamics of dry friction J. Phys. A: Math. Gen. 30 (1997) 6057 6063. Printed in the UK PII: S0305-4470(97)79385-5 Nonlinear dynamics of dry friction Franz-Josef Elmer Institut für Physik, Universität Basel, CH-4056 Basel, Switzerland

More information

VTU-NPTEL-NMEICT Project

VTU-NPTEL-NMEICT Project MODULE-II --- SINGLE DOF FREE S VTU-NPTEL-NMEICT Project Progress Report The Project on Development of Remaining Three Quadrants to NPTEL Phase-I under grant in aid NMEICT, MHRD, New Delhi SME Name : Course

More information

Giacomo Boffi. Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano

Giacomo Boffi.  Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April 21, 2017 Outline of Structural Members Elastic-plastic Idealization

More information

Modelling biological oscillations

Modelling biological oscillations Modelling biological oscillations Shan He School for Computational Science University of Birmingham Module 06-23836: Computational Modelling with MATLAB Outline Outline of Topics Van der Pol equation Van

More information

Passivity-based Stabilization of Non-Compact Sets

Passivity-based Stabilization of Non-Compact Sets Passivity-based Stabilization of Non-Compact Sets Mohamed I. El-Hawwary and Manfredi Maggiore Abstract We investigate the stabilization of closed sets for passive nonlinear systems which are contained

More information

ODE Homework 1. Due Wed. 19 August 2009; At the beginning of the class

ODE Homework 1. Due Wed. 19 August 2009; At the beginning of the class ODE Homework Due Wed. 9 August 2009; At the beginning of the class. (a) Solve Lẏ + Ry = E sin(ωt) with y(0) = k () L, R, E, ω are positive constants. (b) What is the limit of the solution as ω 0? (c) Is

More information

Small-scale instabilities in dynamical systems with sliding. Sieber, J and Kowalczyk, P. MIMS EPrint:

Small-scale instabilities in dynamical systems with sliding. Sieber, J and Kowalczyk, P. MIMS EPrint: Small-scale instabilities in dynamical systems with sliding Sieber, J and Kowalczyk, P 2010 MIMS EPrint: 2010.2 Manchester Institute for Mathematical Sciences School of Mathematics The University of Manchester

More information

1 Trajectory Generation

1 Trajectory Generation CS 685 notes, J. Košecká 1 Trajectory Generation The material for these notes has been adopted from: John J. Craig: Robotics: Mechanics and Control. This example assumes that we have a starting position

More information

DYNAMICAL SYSTEMS

DYNAMICAL SYSTEMS 0.42 DYNAMICAL SYSTEMS Week Lecture Notes. What is a dynamical system? Probably the best way to begin this discussion is with arguably a most general and yet least helpful statement: Definition. A dynamical

More information

Lecture 9 Nonlinear Control Design. Course Outline. Exact linearization: example [one-link robot] Exact Feedback Linearization

Lecture 9 Nonlinear Control Design. Course Outline. Exact linearization: example [one-link robot] Exact Feedback Linearization Lecture 9 Nonlinear Control Design Course Outline Eact-linearization Lyapunov-based design Lab Adaptive control Sliding modes control Literature: [Khalil, ch.s 13, 14.1,14.] and [Glad-Ljung,ch.17] Lecture

More information

Tracking Control for a Class of Measure Differential Inclusions

Tracking Control for a Class of Measure Differential Inclusions Proceedings of the 47th IEEE Conference on Decision and Control Cancun, Mexico, Dec. 9-, 2008 Tracking Control for a Class of Measure Differential Inclusions Nathan van de Wouw, Remco I. Leine Abstract

More information

2 Discrete growth models, logistic map (Murray, Chapter 2)

2 Discrete growth models, logistic map (Murray, Chapter 2) 2 Discrete growth models, logistic map (Murray, Chapter 2) As argued in Lecture 1 the population of non-overlapping generations can be modelled as a discrete dynamical system. This is an example of an

More information

DETC ON THE COUPLING BETWEEN TORSIONAL AND LATERAL VIBRATIONS IN A ROTOR DYNAMIC SYSTEM WITH SET-VALUED FRICTION

DETC ON THE COUPLING BETWEEN TORSIONAL AND LATERAL VIBRATIONS IN A ROTOR DYNAMIC SYSTEM WITH SET-VALUED FRICTION Proceedings of the ASME 7 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference IDETC/CIE 7 September -7, 7, Las Vegas, Nevada, USA DETC7-91 ON THE

More information

Lecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. Poincare-Bendixson theorem Hopf bifurcations Poincare maps

Lecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. Poincare-Bendixson theorem Hopf bifurcations Poincare maps Lecture 5 Outline: Limit Cycles Definition and examples How to rule out limit cycles Gradient systems Liapunov functions Dulacs criterion Poincare-Bendixson theorem Hopf bifurcations Poincare maps Limit

More information

CHAPTER 3 THE MAXIMUM PRINCIPLE: MIXED INEQUALITY CONSTRAINTS. p. 1/73

CHAPTER 3 THE MAXIMUM PRINCIPLE: MIXED INEQUALITY CONSTRAINTS. p. 1/73 CHAPTER 3 THE MAXIMUM PRINCIPLE: MIXED INEQUALITY CONSTRAINTS p. 1/73 THE MAXIMUM PRINCIPLE: MIXED INEQUALITY CONSTRAINTS Mixed Inequality Constraints: Inequality constraints involving control and possibly

More information

Rotor-stator contact in multiple degree of freedom systems

Rotor-stator contact in multiple degree of freedom systems Rotor-stator contact in multiple degree of freedom systems A. D. Shaw 1, A. R. Champneys 2, M. I. Friswell 1 1 Swansea University, College of Engineering, UK 2 University of Bristol, Department of Engineering

More information

AEROELASTICITY IN AXIAL FLOW TURBOMACHINES

AEROELASTICITY IN AXIAL FLOW TURBOMACHINES von Karman Institute for Fluid Dynamics Lecture Series Programme 1998-99 AEROELASTICITY IN AXIAL FLOW TURBOMACHINES May 3-7, 1999 Rhode-Saint- Genèse Belgium STRUCTURAL DYNAMICS: BASICS OF DISK AND BLADE

More information

Engineering Mechanics Prof. Siva Kumar Department of Civil Engineering Indian Institute of Technology, Madras Statics - 5.2

Engineering Mechanics Prof. Siva Kumar Department of Civil Engineering Indian Institute of Technology, Madras Statics - 5.2 Engineering Mechanics Prof. Siva Kumar Department of Civil Engineering Indian Institute of Technology, Madras Statics - 5.2 Now what we want to do is given a surface which is let s assume the surface is

More information

KEELE UNIVERSITY PHYSICS/ASTROPHYSICS MODULE PHY OSCILLATIONS AND WAVES PRACTICE EXAM

KEELE UNIVERSITY PHYSICS/ASTROPHYSICS MODULE PHY OSCILLATIONS AND WAVES PRACTICE EXAM KEELE UNIVERSITY PHYSICS/ASTROPHYSICS MODULE PHY-10012 OSCILLATIONS AND WAVES PRACTICE EXAM Candidates should attempt ALL of PARTS A and B, and TWO questions from PART C. PARTS A and B should be answered

More information

Introduction to Geotechnical Earthquake Engineering

Introduction to Geotechnical Earthquake Engineering Module 1 Introduction to Geotechnical Earthquake Engineering by Dr. Deepankar Choudhury Professor Department of Civil Engineering IIT Bombay, Powai, Mumbai 400 076, India. Email: dc@civil.iitb.ac.in URL:

More information

APPPHYS 217 Tuesday 6 April 2010

APPPHYS 217 Tuesday 6 April 2010 APPPHYS 7 Tuesday 6 April Stability and input-output performance: second-order systems Here we present a detailed example to draw connections between today s topics and our prior review of linear algebra

More information