C.-H. Lamarque. University of Lyon/ENTPE/LGCB & LTDS UMR CNRS 5513
|
|
- Marilynn Dorsey
- 5 years ago
- Views:
Transcription
1 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 Structures Sperlonga (Italy), September 213 C.-H. Lamarque University of Lyon/ENTPE/LGCB & LTDS UMR CNRS 5513
2 1. Numerical simulations and question of identification 2. Identification process: Theoretical examples 3. Identification: Numerical examples
3 Introduction Numerical simulations and question of identification : Responses of some studied examples: Numerical examples, Transient, cycles, responses under loading solicitation Numerics... = "experimental" results Identification: Theoretical examples via modelling Assuming type for non smooth elements: E.g. dry friction, Saint-Venant Theoretical location of "non smooth" points of cycles Identification: Numerical examples For discrete data (numerical, experimental), how to process? Adapted numerica tools
4 Schedule 1. Numerical simulations and question of identification 2. Identification process: Theoretical examples 3. Identification: Numerical examples
5 Description of the mechanical systems We studied problems with the general mathematical formulation: Proposition: Under assumptions, for all ξ D( φ), there exists a unique function X in W 1,1 (, T, R p ) such that { Ẋ (t) + M φ(x (t)) G(t, X (t)) a.e. on ], T [, (1) X () = ξ. M symmetric definite positive matrix, φ maximal monotone... In every case, models are based on masses, springs (stiffness), dashpots (damping coefficient), Saint-Venant elements (threshold). These parameters have to be identified.
6 Thus, all the systems we have seen can be written under the form (1) and have a unique solution. We can observe that there are three classes of mechanical systems : in the first class, the function φ is a linear combinations of x i ; in the second class, the function φ is a sum of functions ψ [ ηi,η i ] (indicatrix function); in the third class, φ involves the both functions x and ψ [ ηi,η i ]. But a model must be given to be identified. Let us look at a few numerical results...
7 Numerical simulations We present some numerical simulations for : - the Prandtl model, - the Prandtl model with linear hardening, - the generalized Prandtl model with linear hardening - several viscoelastoplastic models, For all these numerical simulations, we solve the differential inclusion (1) by using the implicit Euler scheme We denote by X h the linear interpolation of the X n s. The function X h converges to the solution X of system (1) in C ([, T ], R p ).
8 Study of the rheological Prandtl models Let us consider the rheological Prandtl models: k α 1 1 k k α 2 2 k 1 α 1 m F k α 2 2 kn αn m F (a) kn αn The generalized Prandtl model (a) and the generalized Prandtl model + linear hardening (b). Displacement x and mass m. We choose m = 1. We study a harmonic forcing and other periodic forcings. (b)
9 Harmonic forcing For this section, we choose F (t) = f cos(ωt). (2) We first present essential difference between the Prandtl model and the Prandtl model with linear hardening : we compute x, y and z for the Prandtl model on the interval [, T ] with T = 2, F (t) = 2 cos(.1t), η = k = 1, x = y = u =, (3) and for the Prandtl model with linear hardening onl [, T ] with T = 3, F (t) = 2 cos(.1t), k = 1, η = k = 1, Sx = y = u =. (4) Results :
10 x(t) x (t) u(t) t t t x(t) x (t) u(t) t t t The functions x, y and u for the Prandtl model The functions x, y and u for the Prandtl model with linear hardening,
11 We can see that the amplitude of the functions x and ẋ is larger for the Prandtl model, which corresponds to k = than for Prandtl model with hardening, which is the same system with k. We sketch a physical explanation: If the St-Venant element slips, then and u εη, ẍ(t) + k m x(t) = f εkη cos(ωt) m m, (5) where ε { 1, 1}. If k =, and if the slip phase starts at t, then x(t) = 2f sin(ω(t + t mω 2 )) sin(ω(t t )) εkη 2m (t t ) 2 + ( ẋ(t ) f mω sin(ωt ) ) (6) (t t ) + x(t ). For ω 1 and t t, we have x(t) 2f mω 2 sin(ω(t t )); (7) thus, the approximate amplitude of x is roughly 2f = 4, (8) mω2
12 If k = 1 and if the slip phase starts at t, (5) implies x(t) = D cos(ω t + φ) + f εkη m(ω 2 cos(ωt) ω2 ) mω 2, (9) where ω = k /m ω and (D, φ) depends of (x(t ), ẋ(t )). Here, with ω =.1, the amplitude of the sinusoidal component of period ω of the function x is roughly as we can see on the right. - Coherent results. f m(ω 2 2, (1) ω2 ) - We start observing also "smoothness" of the solutions.
13 To examine "smoothness", see in next Fig. a) curve {x(t), F (t)} t [35,8] for the Prandtl model with F (t) = 2 cos(.1t), η = k = 1, and x = y = u =. We observe a limit cycle, but it gives no information on the physically relevant parameters of the system. If we plot as in Fig. b), {x(t), F (t) mẍ(t)} t [35,8], we observe that the slope of the oblique parts of the cycle is equal to k (a) (b) F(t) 5 1 F(t) mx (t) x(t) x(t) {x(t), F (t)} t [35,8] (a) and {x(t), F (t) mẍ(t)} t [35,8] (b) for the Prandtl model defined by F (t) = 2 cos(.1t), η = k = 1, and x = y = u =. These two figures differ only by the choice of coordinates.
14 t (a) F(t) mx (t) F(t) mx"(t) x(t) x(t) {x(t), F (t) mẍ(t), t} t [,18] for the Prandtl model, defined by F (t) = 2 cos(.1t), η = k = 1, x = y = u =. {x(t), F (t) mẍ(t)} t [5,2] for the generalized Prandtl model with linear hardening F (t) = cos(.5t), n = 1, k =, η 1 = k 1 = 1, (x = y = u,1 = ) We can observe (left Figure) transients followed by a periodic regime (t 15)of hysteresis cycles. Without transient (right Figure a), the cycle does not contain (, ) in the plane (x, F mẍ). Slopes and location of slopes changes...
15 F(t) mx (t) 2 F(t) mx (t) x(t) x(t) {x(t), F (t) mẍ(t)} t [75,1] for the generalized Prandtl rheological model with linear hardening, with mechanical parameters and F (t) = f cos(.5t), n = 5, m = 1, x = y =, i {1,..., 5}, u,i =, f = 6.6. {x(t), F (t) mẍ(t)} t [96,1] for the generalized Prandtl rheological model with linear hardening, defined by F (t) = 8H 1 (.5t), n = 3, k = 1/n, x = y =, i {1,..., 3}, k i = 1/3, η i = i, u,i =.
16 Viscoelastoplastic model Adding classical dissipation cẋ in our one dof model... 1 (a) (b) F(t) mx (t) F(t) mx (t) x(t) x(t) {x(t), F (t) mẍ(t)} t [4,2] for the viscoelastoplastic model defined by F (t) = cos(.5t), η = k = m = 1, x = y = u =, c = 1 (a) and c = 1 (b). These two figures differ only by (mathematical) value of damping c.
17 Schedule 1. Numerical simulations and question of identification 2. Identification process: Theoretical examples 3. Identification: Numerical examples
18 Process for identification To get such cycle or at least half cycle: Continuous data (theory), Discrete data (experiments), F known, To deal with damping... Find slope s changes (non smooth), Identify from a given model. Let us start first from theoretical point of view. Remark: - sometimes hysteresis limit cycles periodic but not convex. - elastic part may be nonlinear... - in both cases: Similar process works.
19 Most of the responses of the generalized Prandtl model to a sinusoidal forcing or to a periodic (but non harmonic) forcing tend to hysteresis limit cycles. We observe that these hysteresis cycles have a center of symmetry. Therefore, we study a loading phase corresponding to a half-cycle in the (x, F mẍ) plan : F(t)-mx"(t) A n+1 A n A n+2 A i+1 A i A 2 A 1 d 1 d i d n-1 d n x(t) Loading curve for the generalized Prandtl rheological model with linear hardening.
20 Let us denote by A 1, A 2,...,A n+1 and A n+2 the ends of segments which constitute the hysteresis half-cycle. For i {1,..., n}, d i denotes the difference between abscissa of A i+1 and A i (positive real). For i {1,..., n + 1}, p i denotes the slope of the segment [A i, A i+1 ]. F(t)-mx"(t) A n A n+1 A n+2 A i+1 A i A 2 A 1 d 1 d i d n-1 d n x(t)
21 If the forcing F is not constant on any open non empty subinterval of [, T ], then x has the same property; assume indeed that x is constant on (t 1, t 2 ); then ẋ and ẍ vanish on (t 1, t 2 ) and by uniqueness of the solution of { u i + β(u i /η i ), (11) u i (t 1 ) given, u i (t) is equal to u i (t 1 ) on (t 1, t 2 ). This implies that the expression F = mẍ + k x + n k i u i i=1 is constant over (t 1, t 2 ), which contradicts our assumption.
22 We are now able to describe the shape of the representation of a trajectory of the system in the x and F mẍ coordinates, in some special case. We assume that the η i = α i /k i are all distinct; otherwise two elements with identical ratio η i would enter the plasticity phase simultaneously. We reorder the indices so that η 1 < η 2 <... < η n 1 < η n. (12) Assume that F vanishes on no open non empty subinterval of [, T ], and the following properties: 1) u j () = η j, j = 1,..., n, 2) x is increasing on [, T ].
23 Then there exists a increasing sequence t 1 t 2... t n such that for some j {1,..., n} and we have moreover t 1 < t 2 <... < t j = t j+1 =... = T (13) x(t j ) x() = 2η j, (14) n F mẍ k + x is constant on ]t j, t j+1 ]. (15) i=j+1 k i
24 Observe first that x is strictly increasing on [, T ], thanks to our assumption on F. Let û j (t) = η j + x(t) x(), (16) and denote by ξ the inverse function of x and define { ξ (2η j + x()) if 2η j + x() < x(t ), t j = T if 2η j + x() x(t ). (17) We check immediately that the function u j (t) = is the (unique) solution of { {ûj (t) if t < t j, η j if t t j, u j + β(u j /η j ) ẋ, u j () = η j. (18) (19)
25 Therefore, on ]t j, t j+1 ], relation (15) holds. Moreover, relation (17) is equivalent to (14). We write equations (15) and (14) as j {1,..., n + 1}, p j = k + n k l, (2) i {1,..., n}, d i = 2η i. (21) On the other hand, we remark that the assumptions k, k i > and equation (2) imply equations (12) and (21) imply l=j p n+1 < p n <... < p 2 < p 1 ; (22) d 1 < d 2 <... < d n 1 < d n. (23)
26 From (2) and (21) we have a one to one correspondence between the parameters of generalized Prandtl model k i and η i and geometrical parameters p j and d j of the hysteresis cycle in the (x, F mẍ) plane. Thus, a partial identification of the model is possible. But we cannot infer from the values of (p i ) 1 i n+1 and (d i ) 1 i n the abscissa and the ordinate of the lower left point of the cycle. They depend on the initial conditions and on the forcing:
27 (a) (b) F(t) mx (t) F(t) mx (t) x(t) x(t) x 1 4 {x(t), F (t) mẍ(t)} t [5,2] for the generalized Prandtl rheological model with linear hardening, with F(t) = cos(.5t), n = 1, k =, η 1 = k 1 = 1, (x = y = u,1 = ) (a), (x = 1, y = 2, u,1 =.99) (b). These two figures differ only by the initial conditions.
28 There are more or less parameters to identify...: F(t) mx (t) 2 2 F(t) mx (t) x(t) x(t) {x(t), F (t) mẍ(t)} t [75,1] for the generalized Prandtl rheological model with linear hardening, with mechanical parameters and F (t) = f cos(.5t), n = 5, m = 1, x = y =, i {1,..., 5}, u,i =, f = phases... (n) {x(t), F (t) mẍ(t)} t [96,1] for the generalized Prandtl rheological model with linear hardening, defined by F (t) = 8H 1 (.5t), n = 3, k = 1/n, x = y =, i {1,..., 3}, k i = 1/3, η i = i, u,i =. 3 phases... (n)
29 Schedule 1. Numerical simulations and question of identification 2. Identification process: Theoretical examples 3. Identification: Numerical examples
30 Numerical identification of a model from its limit cycles If limit cycle or at least half cycle is known from continuous data, with no damping or known damping, F known: OK (calculation of parameters...) If limit cycle or at least half cycle is known from discrete data with no noise, high quality sampling, with no damping or known damping, F known (discrete): OK (calculation of parameters, mean squares,...). If limit cycle or at least half cycle is known from discrete data with poor sampling, noise etc. (real data), with no damping or known damping, F known (discrete): OK with wavelet analysis. To deal with damping...also wavelet analysis... Identification from a given model and high enough solicitation: n = 1 rather than n = 2 leads to an error. Unlocking only one Saint-Venant element rather than two because of too weak forcing F also leads to an error... obviously.
31 Damping with wavelets Continuous wavelet transform: f signal, g analyzing function (b, a) R R +, W g f (b, a) = 1 + a f (t)g( t b )dt, (24) a For f (t) = A exp( ct) cos(ωt + φ), focusing on b j = jt, j N, T = 2π/ω, one can prove c = 1 (m n)t log( W g f (nt, a) (25) (mt, a)), W g f for m, n integers, at the analyzing scale a > small enough. C.-H. LAMARQUE, S. PERNOT, A. CUER, Damping identification in multi-degree-of-freedom systems via a wavelet-logarithmic decrement- Part I: Theory, Journal of Sound and Vibration, 235 (3), , 2.
32 Multiresolution wavelet analysis: f signal, φ scaling function, φ jk (x) = r Z φ jk (x + r), (26) f f J = 2 J 1 k= s Jk φjk (27) For f (t) = A exp( ct) cos(ωt + φ), focusing at scale a = 2 ja on b k = k/2 ja and b l = l/2 ja, one can prove for k, l, j a integers. c = 2ja (k l) log( s j al s jak ), (28) C.-H. LAMARQUE, S. PERNOT, A. CUER, Damping identification in multi-degree-of-freedom systems via a wavelet-logarithmic decrement- Part I: Theory, Journal of Sound and Vibration, 235 (3), , 2.
33 Theoretical and signal processing validation: C.-H. LAMARQUE, S. PERNOT, A. CUER, Damping identification in multi-degree-of-freedom systems via a wavelet-logarithmic decrement- Part I: Theory, Journal of Sound and Vibration, 235 (3), , 2. Experimental validation: S. HANS, E. IBRAIM, S.PERNOT, C. BOUTIN, C.-H. LAMARQUE, Damping identification in multi-degree-of-freedom systems via a wavelet-logarithmic decrement- Part II: Study of a civil engineering building, Journal of Sound and Vibration, 235 (3), , 2. It permits to process identification of damping from more general signal, discrete, with noise, from experiments. Better than several classical algorithms (logarithmic decrement, mean squares methods, etc.). It includes filters...
34 Identification of non locally Lipschitz with damping mẍ + G(x) + G(ẋ) + i {1,..., n}, n k i u i = F, i=1 u i + β(u i /η i ) ẋ (29) where F is deterministic or stochastic. We choose n = 5, k =, m = 1, σ = 3, T = 1 and G(y) =.5y, G(x) =.5x 3. We focus on a monotonic phase of x assumed to be hold between and T. On the graph (x, F mẍ) t T = (X, G(X )) some regularity s phases can be observed. In fact, in the deterministic case, ẋ is a derivable function with respect to x and we observe only n discontinuities on the graph (X, G (X )). We can locate this discontinuities by using wavelet s tools.
35 Identification of non locally Lipschitz with damping Slope s detections in next figure for n = 5, k =, m = 1, σ = 3, T = 1 and the functions G 1 (x, y) = G(y) =.5y and G 2 (x, y) =. Horizontal one: we represent on the horizontal axis not the value of the discretization (x i ) 1 i N of x but the number i. Vertical axis = scales Color = highest amplitude of wavelets coefficients
36 Identification of non locally Lipschitz with damping x x 1 4 Amplitude of wavelet s coefficients versus displacement x. High amplitudes separation of phases lest square identification... Remark: Identification may be still difficult (external forcing, model of damping, choice of model)...
Dr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum
STRUCTURAL DYNAMICS Dr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum Overview of Structural Dynamics Structure Members, joints, strength, stiffness, ductility Structure
More informationDynamics 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 informationDESIGN OF A NONLINEAR VIBRATION ABSORBER
DESIGN OF A NONLINEAR VIBRATION ABSORBER Maxime Geeroms, Laurens Marijns, Mia Loccufier and Dirk Aeyels Ghent University, Department EESA, Belgium Abstract Linear vibration absorbers can only capture certain
More informationSome 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 informationOSE801 Engineering System Identification. Lecture 09: Computing Impulse and Frequency Response Functions
OSE801 Engineering System Identification Lecture 09: Computing Impulse and Frequency Response Functions 1 Extracting Impulse and Frequency Response Functions In the preceding sections, signal processing
More informationIdentification 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 informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationIntroduction to Vibration. Professor Mike Brennan
Introduction to Vibration Professor Mie Brennan Introduction to Vibration Nature of vibration of mechanical systems Free and forced vibrations Frequency response functions Fundamentals For free vibration
More information2. Determine whether the following pair of functions are linearly dependent, or linearly independent:
Topics to be covered on the exam include: Recognizing, and verifying solutions to homogeneous second-order linear differential equations, and their corresponding Initial Value Problems Recognizing and
More informationSolutions for homework 5
1 Section 4.3 Solutions for homework 5 17. The following equation has repeated, real, characteristic roots. Find the general solution. y 4y + 4y = 0. The characteristic equation is λ 4λ + 4 = 0 which has
More informationProblem Set Number 01, MIT (Winter-Spring 2018)
Problem Set Number 01, 18.377 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Thursday, March 8, 2018. Turn it in (by 3PM) at the Math.
More informationVibrations: Second Order Systems with One Degree of Freedom, Free Response
Single Degree of Freedom System 1.003J/1.053J Dynamics and Control I, Spring 007 Professor Thomas Peacock 5//007 Lecture 0 Vibrations: Second Order Systems with One Degree of Freedom, Free Response Single
More informationME 563 HOMEWORK # 7 SOLUTIONS Fall 2010
ME 563 HOMEWORK # 7 SOLUTIONS Fall 2010 PROBLEM 1: Given the mass matrix and two undamped natural frequencies for a general two degree-of-freedom system with a symmetric stiffness matrix, find the stiffness
More informationDynamics of Structures: Theory and Analysis
1. Free vibrations 2. Forced vibrations 3. Transient response 4. Damping mechanisms Dynamics of Structures: Theory and Analysis Steen Krenk Technical University of Denmark 5. Modal analysis I: Basic idea
More informationDamping Identification with the Morlet-Wave
Damping Identification with the Morlet-Wave Janko Slavič, Miha Boltežar Laboratory for Dynamics of Machines and Structures, Faculty of Mechanical Engineering, University of Ljubljana, Aškerčeva 6, 1000
More information1-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 informationThe student will experimentally determine the parameters to represent the behavior of a damped oscillatory system of one degree of freedom.
Practice 3 NAME STUDENT ID LAB GROUP PROFESSOR INSTRUCTOR Vibrations of systems of one degree of freedom with damping QUIZ 10% PARTICIPATION & PRESENTATION 5% INVESTIGATION 10% DESIGN PROBLEM 15% CALCULATIONS
More informationTHEORY OF VIBRATION ISOLATION
CHAPTER 30 THEORY OF VIBRATION ISOLATION Charles E. Crede Jerome E. Ruzicka INTRODUCTION Vibration isolation concerns means to bring about a reduction in a vibratory effect. A vibration isolator in its
More informationDamped harmonic motion
Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,
More information4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
More informationMulti Degrees of Freedom Systems
Multi Degrees of Freedom Systems MDOF s http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March 9, 07 Outline, a System
More informationDynamics 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 informationOutline of parts 1 and 2
to Harmonic Loading http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March, 6 Outline of parts and of an Oscillator
More informationBSc/MSci EXAMINATION. Vibrations and Waves. Date: 4 th May, Time: 14:30-17:00
BSc/MSci EXAMINATION PHY-217 Vibrations and Waves Time Allowed: 2 hours 30 minutes Date: 4 th May, 2011 Time: 14:30-17:00 Instructions: Answer ALL questions in section A. Answer ONLY TWO questions from
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic angular frequency
More informationSection 3.7: Mechanical and Electrical Vibrations
Section 3.7: Mechanical and Electrical Vibrations Second order linear equations with constant coefficients serve as mathematical models for mechanical and electrical oscillations. For example, the motion
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ small angle approximation θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic
More informationMath 1302, Week 8: Oscillations
Math 302, Week 8: Oscillations T y eq Y y = y eq + Y mg Figure : Simple harmonic motion. At equilibrium the string is of total length y eq. During the motion we let Y be the extension beyond equilibrium,
More informationDifferential 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 informationMechanics IV: Oscillations
Mechanics IV: Oscillations Chapter 4 of Morin covers oscillations, including damped and driven oscillators in detail. Also see chapter 10 of Kleppner and Kolenkow. For more on normal modes, see any book
More informationDynamics of Non-Smooth Systems
Dynamics of Non-Smooth Systems P Chandramouli Indian Institute of Technology Madras October 5, 2013 Mouli, IIT Madras Non-Smooth System Dynamics 1/39 Introduction Filippov Systems Behaviour around hyper-surface
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t
More informationMath Assignment 5
Math 2280 - Assignment 5 Dylan Zwick Fall 2013 Section 3.4-1, 5, 18, 21 Section 3.5-1, 11, 23, 28, 35, 47, 56 Section 3.6-1, 2, 9, 17, 24 1 Section 3.4 - Mechanical Vibrations 3.4.1 - Determine the period
More informationSolutions of Spring 2008 Final Exam
Solutions of Spring 008 Final Exam 1. (a) The isocline for slope 0 is the pair of straight lines y = ±x. The direction field along these lines is flat. The isocline for slope is the hyperbola on the left
More informationGiacomo 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 informationEN40: Dynamics and Vibrations. Final Examination Wed May : 2pm-5pm
EN40: Dynamics and Vibrations Final Examination Wed May 10 017: pm-5pm School of Engineering Brown University NAME: General Instructions No collaboration of any kind is permitted on this examination. You
More informationPHYSICS 110A : CLASSICAL MECHANICS HW 2 SOLUTIONS. Here is a sketch of the potential with A = 1, R = 1, and S = 1. From the plot we can see
PHYSICS 11A : CLASSICAL MECHANICS HW SOLUTIONS (1) Taylor 5. Here is a sketch of the potential with A = 1, R = 1, and S = 1. From the plot we can see 1.5 1 U(r).5.5 1 4 6 8 1 r Figure 1: Plot for problem
More informationOscillations. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring of constant k is
Dr. Alain Brizard College Physics I (PY 10) Oscillations Textbook Reference: Chapter 14 sections 1-8. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring
More informationVector hysteresis models
Euro. Jnl. Appl. Math. 2 (99), 28 292 Vector hysteresis models Pavel Krejčí Matematický ústav ČSAV, Žitná 25 5 67 Praha, Czechoslovakia Key words: vector hysteresis operator, hysteresis potential, differential
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationWEEKS 8-9 Dynamics of Machinery
WEEKS 8-9 Dynamics of Machinery References Theory of Machines and Mechanisms, J.J.Uicker, G.R.Pennock ve J.E. Shigley, 2011 Mechanical Vibrations, Singiresu S. Rao, 2010 Mechanical Vibrations: Theory and
More informationDynamic Loads CE 543. Examples. Harmonic Loads
CE 543 Structural Dynamics Introduction Dynamic Loads Dynamic loads are time-varying loads. (But time-varying loads may not require dynamic analysis.) Dynamics loads can be grouped in one of the following
More informationVibrations Qualifying Exam Study Material
Vibrations Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering vibrations topics. These topics are listed below for clarification. Not all instructors
More informationPhysics 8 Monday, December 4, 2017
Physics 8 Monday, December 4, 2017 HW12 due Friday. Grace will do a review session Dec 12 or 13. When? I will do a review session: afternoon Dec 17? Evening Dec 18? Wednesday, I will hand out the practice
More informationNumerical Solution of Equation of Motion
Class Notes: Earthquake Engineering, Ahmed Elgamal, September 25, 2001 (DRAFT) Numerical Solution of Equation of Motion Average Acceleration Method (Trapezoidal method) m a + c v + k d = f (t) In the above
More informationExplosion Protection of Buildings
Explosion Protection of Buildings Author: Miroslav Mynarz Explosion Protection of Buildings Introduction to the Problems of Determination of Building Structure's Response 3 Classification of actions According
More informationIntroduction to structural dynamics
Introduction to structural dynamics p n m n u n p n-1 p 3... m n-1 m 3... u n-1 u 3 k 1 c 1 u 1 u 2 k 2 m p 1 1 c 2 m2 p 2 k n c n m n u n p n m 2 p 2 u 2 m 1 p 1 u 1 Static vs dynamic analysis Static
More informationTheoretical Basis of Modal Analysis
American Journal of Mechanical Engineering, 03, Vol., No. 7, 73-79 Available online at http://pubs.sciepub.com/ajme//7/4 Science and Education Publishing DOI:0.69/ajme--7-4 heoretical Basis of Modal Analysis
More informationThe 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 informationFree Vibration of Single-Degree-of-Freedom (SDOF) Systems
Free Vibration of Single-Degree-of-Freedom (SDOF) Systems Procedure in solving structural dynamics problems 1. Abstraction/modeling Idealize the actual structure to a simplified version, depending on the
More information5th-order differentiation
ARBITRARY-ORDER REAL-TIME EXACT ROBUST DIFFERENTIATION A. Levant Applied Mathematics Dept., Tel-Aviv University, Israel E-mail: levant@post.tau.ac.il Homepage: http://www.tau.ac.il/~levant/ 5th-order differentiation
More informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More informationSoftware Verification
POGAM NAME: EVISION NO.: 0 EXAMPLE 6-005 LINK DAMPE ELEMENT UNDE HAMONIC LOADING POBLEM DESCIPTION In this single degree of freedom example a spring-mass-damper system is subjected to a harmonic load.
More informationExercises 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 informationComparison between the visco-elastic dampers And Magnetorheological dampers and study the Effect of temperature on the damping properties
Comparison between the visco-elastic dampers And Magnetorheological dampers and study the Effect of temperature on the damping properties A.Q. Bhatti National University of Sciences and Technology (NUST),
More informationChapter 3 Mathematical Methods
Chapter 3 Mathematical Methods Slides to accompany lectures in Vibro-Acoustic Design in Mechanical Systems 0 by D. W. Herrin Department of Mechanical Engineering Lexington, KY 40506-0503 Tel: 859-8-0609
More informationFinal Exam December 20, 2011
Final Exam December 20, 2011 Math 420 - Ordinary Differential Equations No credit will be given for answers without mathematical or logical justification. Simplify answers as much as possible. Leave solutions
More informationThe 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 informationLaboratory notes. Torsional Vibration Absorber
Titurus, Marsico & Wagg Torsional Vibration Absorber UoB/1-11, v1. Laboratory notes Torsional Vibration Absorber Contents 1 Objectives... Apparatus... 3 Theory... 3 3.1 Background information... 3 3. Undamped
More informationANALYTICAL MODEL FOR HIGH DAMPING ELASTOMERS APPLIED TO ENERGY DISSIPATING DEVICES. NUMERICAL STUDY AND EXPERIMENTAL VALIDATION
Proceedings of the Third European Conference on Structural Control, 3ECSC, 12-15 July 2004, Vienna University of Technology, Vienna, Austria Austria ANALYTICAL MODEL FOR HIGH DAMPING ELASTOMERS APPLIED
More informationODE 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 informationContinuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms
Continuum mechanics office Math 0.107 ales.janka@unifr.ch http://perso.unifr.ch/ales.janka/mechanics Mars 16, 2011, Université de Fribourg 1. Constitutive equation: definition and basic axioms Constitutive
More informationSlide 1. Slide 2. Slide 3 Remark is a new function derived from called derivative. 2.2 The derivative as a Function
Slide 1 2.2 The derivative as a Function Slide 2 Recall: The derivative of a function number : at a fixed Definition (Derivative of ) For any number, the derivative of is Slide 3 Remark is a new function
More information21.55 Worksheet 7 - preparation problems - question 1:
Dynamics 76. Worksheet 7 - preparation problems - question : A coupled oscillator with two masses m and positions x (t) and x (t) is described by the following equations of motion: ẍ x + 8x ẍ x +x A. Write
More informationDynamics of structures
Dynamics of structures 2.Vibrations: single degree of freedom system Arnaud Deraemaeker (aderaema@ulb.ac.be) 1 Outline of the chapter *One degree of freedom systems in real life Hypothesis Examples *Response
More informationPhysics 2101 S c e t c i cti n o 3 n 3 March 31st Announcements: Quiz today about Ch. 14 Class Website:
Physics 2101 Section 3 March 31 st Announcements: Quiz today about Ch. 14 Class Website: http://www.phys.lsu.edu/classes/spring2010/phys2101 3/ http://www.phys.lsu.edu/~jzhang/teaching.html Simple Harmonic
More informationSubharmonic Oscillations and Chaos in Dynamic Atomic Force Microscopy
Subharmonic Oscillations and Chaos in Dynamic Atomic Force Microscopy John H. CANTRELL 1, Sean A. CANTRELL 2 1 NASA Langley Research Center, Hampton, Virginia 23681, USA 2 NLS Analytics, LLC, Glencoe,
More informationModeling and Experimentation: Mass-Spring-Damper System Dynamics
Modeling and Experimentation: Mass-Spring-Damper System Dynamics Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin July 20, 2014 Overview 1 This lab is meant to
More informationCMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation
CMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University September 26, 2005 1 Sinusoids Sinusoids
More informationTorsion Spring Oscillator with Dry Friction
Torsion Spring Oscillator with Dry Friction Manual Eugene Butikov Annotation. The manual includes a description of the simulated physical system and a summary of the relevant theoretical material for students
More informationDynamics of structures
Dynamics of structures 2.Vibrations: single degree of freedom system Arnaud Deraemaeker (aderaema@ulb.ac.be) 1 One degree of freedom systems in real life 2 1 Reduction of a system to a one dof system Example
More informationSeminar 6: COUPLED HARMONIC OSCILLATORS
Seminar 6: COUPLED HARMONIC OSCILLATORS 1. Lagrangian Equations of Motion Let consider a system consisting of two harmonic oscillators that are coupled together. As a model, we will use two particles attached
More informationSession 1 : Fundamental concepts
BRUFACE Vibrations and Acoustics MA1 Academic year 17-18 Cédric Dumoulin (cedumoul@ulb.ac.be) Arnaud Deraemaeker (aderaema@ulb.ac.be) Exercise 1 Session 1 : Fundamental concepts Consider the following
More informationChaotic motion. Phys 750 Lecture 9
Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to
More informationIntroduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.
Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.
More informationSHAKING TABLE DEMONSTRATION OF DYNAMIC RESPONSE OF BASE-ISOLATED BUILDINGS ***** Instructor Manual *****
SHAKING TABLE DEMONSTRATION OF DYNAMIC RESPONSE OF BASE-ISOLATED BUILDINGS ***** Instructor Manual ***** A PROJECT DEVELOPED FOR THE UNIVERSITY CONSORTIUM ON INSTRUCTIONAL SHAKE TABLES http://wusceel.cive.wustl.edu/ucist/
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare http://ocw.mit.edu 5.74 Introductory Quantum Mechanics II Spring 9 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Andrei Tokmakoff,
More informationCopyright 2009, August E. Evrard.
Unless otherwise noted, the content of this course material is licensed under a Creative Commons BY 3.0 License. http://creativecommons.org/licenses/by/3.0/ Copyright 2009, August E. Evrard. You assume
More informationSEISMIC RESPONSE OF SINGLE DEGREE OF FREEDOM STRUCTURAL FUSE SYSTEMS
3 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 4 Paper No. 377 SEISMIC RESPONSE OF SINGLE DEGREE OF FREEDOM STRUCTURAL FUSE SYSTEMS Ramiro VARGAS and Michel BRUNEAU
More informationStructural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma).
Structural Dynamics Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). We will now look at free vibrations. Considering the free
More informationDamping of materials and members in structures
Journal of Physics: Conference Series Damping of materials and members in structures To cite this article: F Orban 0 J. Phys.: Conf. Ser. 68 00 View the article online for updates and enhancements. Related
More informationIn-class exercises Day 1
Physics 4488/6562: Statistical Mechanics http://www.physics.cornell.edu/sethna/teaching/562/ Material for Week 11 Exercises due Mon Apr 16 Last correction at April 16, 2018, 11:19 am c 2018, James Sethna,
More informationMethod of Finite Elements II Modeling ohysteresis
Prof. Dr. Eleni Chatzi, Dr. K. Agathos, Dr. G. Abbiati Institute of Structural Engineering (IBK) Department of Civil, Environmental, and Geomatic Engineering (DBAUG) ETH Zürich Method of Finite Elements
More informationThe 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 informationApplication of bouc-wen model to frequency-dependent nonlinear hysteretic friction damper
Journal of Mechanical Science and Technology 24 (6) () 3~37 www.springerlink.com/content/738-494x DOI.7/s26--4-6 Application of bouc-wen model to frequency-dependent nonlinear hysteretic friction damper
More informationSIMPLE HARMONIC MOTION
SIMPLE HARMONIC MOTION PURPOSE The purpose of this experiment is to investigate simple harmonic motion. We will determine the elastic spring constant of a spring first and then study small vertical oscillations
More information1.1 OBJECTIVE AND CONTENTS OF THE BOOK
1 Introduction 1.1 OBJECTIVE AND CONTENTS OF THE BOOK Hysteresis is a nonlinear phenomenon exhibited by systems stemming from various science and engineering areas: under a low-frequency periodic excitation,
More informationNew implicit method for analysis of problems in nonlinear structural dynamics
Applied and Computational Mechanics 5 (2011) 15 20 New implicit method for analysis of problems in nonlinear structural dynamics A. A. Gholampour a,, M. Ghassemieh a a School of Civil Engineering, University
More informationNumerical Methods for Initial Value Problems; Harmonic Oscillators
Lab 1 Numerical Methods for Initial Value Problems; Harmonic Oscillators Lab Objective: Implement several basic numerical methods for initial value problems (IVPs), and use them to study harmonic oscillators.
More informationMath 221 Topics since the second exam
Laplace Transforms. Math 1 Topics since the second exam There is a whole different set of techniques for solving n-th order linear equations, which are based on the Laplace transform of a function. For
More informationLyapunov exponent calculation of a two-degreeof-freedom vibro-impact system with symmetrical rigid stops
Chin. Phys. B Vol. 20 No. 4 (2011) 040505 Lyapunov exponent calculation of a two-degreeof-freedom vibro-impact system with symmetrical rigid stops Li Qun-Hong( ) and Tan Jie-Yan( ) College of Mathematics
More informationSecond order linear equations
Second order linear equations Samy Tindel Purdue University Differential equations - MA 266 Taken from Elementary differential equations by Boyce and DiPrima Samy T. Second order equations Differential
More informationIntroduction 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 informationMathematics for Engineers II. lectures. Differential Equations
Differential Equations Examples for differential equations Newton s second law for a point mass Consider a particle of mass m subject to net force a F. Newton s second law states that the vector acceleration
More informationSecond Order Linear ODEs, Part II
Craig J. Sutton craig.j.sutton@dartmouth.edu Department of Mathematics Dartmouth College Math 23 Differential Equations Winter 2013 Outline Non-homogeneous Linear Equations 1 Non-homogeneous Linear Equations
More informationMath 216 Second Midterm 20 March, 2017
Math 216 Second Midterm 20 March, 2017 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material
More information18.03SC Practice Problems 14
1.03SC Practice Problems 1 Frequency response Solution suggestions In this problem session we will work with a second order mass-spring-dashpot system driven by a force F ext acting directly on the mass:
More informationa. Follow the Start-Up Procedure in the laboratory manual. Note the safety rules.
Lab #1 - Free Vibration Name: Date: Section / Group: Procedure Steps (from lab manual): a. Follow the Start-Up Procedure in the laboratory manual. Note the safety rules. b. Locate the various springs and
More informationVibrations and waves: revision. Martin Dove Queen Mary University of London
Vibrations and waves: revision Martin Dove Queen Mary University of London Form of the examination Part A = 50%, 10 short questions, no options Part B = 50%, Answer questions from a choice of 4 Total exam
More informationEffects of Damping Ratio of Restoring force Device on Response of a Structure Resting on Sliding Supports with Restoring Force Device
Effects of Damping Ratio of Restoring force Device on Response of a Structure Resting on Sliding Supports with Restoring Force Device A. Krishnamoorthy Professor, Department of Civil Engineering Manipal
More information