Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
|
|
- Brooke Riley
- 5 years ago
- Views:
Transcription
1 Prof. Dr. Eleni Chatzi Lecture March, 2016
2 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural Response
3 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural Response
4 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural Response
5 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural Response
6 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural Response
7 Multiple DOF Systems Let s revisit the simple 2DOF system k 1 F 1 ( t ) k 2 F ( t 2 ) FBD k1x1 F 1 ( t ) k ( x x ) k2( x2 x1) F ( t 2 ) c 1 m 1 x ( t 1 ) c 1 m 2 x 2 ( t ) cx 1 1 m 1 c2( x 2 x 1) c2( x 2 x 1) (Lumped Mass System) m 2 The equations of motion can be written as m 1 ẍ 1 + (c 1 + c 2 )ẋ 1 c 2 ẋ 2 + (k 1 + k 2 )x 1 k 2 x 2 = F 1 (t) m 2 ẍ 2 + c 2 ẋ 2 c 2 ẋ 1 + k 2 x 2 k 2 x 1 = F 2 (t) The system can be written in matrix form as follows: [ ] [ ] [ ] [ ] [ m1 0 ẍ1 c1 + c 2 c 2 ẋ1 k1 + k 2 k m 2 ẍ 2 c 2 c 2 ẋ 2 k 2 k 2 ] [ x1 x 2 ] = [ F1(t) F 2(t) Eq. (1) ]
8 State Space Equation Formulation for MDOF systems 2dof Mass Spring System or otherwise more compactly: Mx d + Cx d + Kx d = F where x d = [ x 1 x 2 ] T We now introduce the augmented state vector: x = [ x 1 x 2 ẋ 1 ẋ 2 ] T. Then, ẋ = [ M 1 K ] [ M 1 C ] x 1 x 2 ẋ 1 ẋ [ M 1 ] [ F1 F 2 ]
9 State Space Equation Formulation 2dof Mass Spring System We obtain the following equivalent 1st order ODE system where u = [ F 1 F 2 ] T ẋ = Ax + Bu Assume you would like to monitor (measure) both displacements x 1, x 2. Then the observation vector is: y = x 1 x 2 x 3 x 4 + O 4 2 u(t) y = Cx + Du
10 State Space Equation Formulation Note Using the state space representation we have converted a 2nd order ODE into an equivalent 1st order ODE system. We can now use any of the aforementioned 1st order ODE integration methods in order to convert the continuous system into a discrete one and obtain an approximate solution For instance MATLAB s ode45, which is a Runge Kutta integration scheme may be used. What are other integration schemes that may be utilized?
11 Numerical Integration for 1st order ODEs Using these methods a continuous system is brought into an equivalent discrete formulation and an approximative solution is sought. 1st order ODE Integration Methods Assume dy dt = f (t, y(t)), y(t 0) = 0 Forward Euler Method y n+1 = y n + hf (t n, y n ) where h is the integration time step. This explicit expression is obtained from the truncated Taylor Expansion of y(t n + h) Backward Euler Method y n+1 = y n + hf (t n+1, y n+1 ) This implicit expression (since y n+1 is on the right hand side) is obtained from the truncated Taylor Expansion of y(t n+1 h)
12 Numerical Integration for 1st order ODEs 2nd Order Runge Kutta (RK2) k 1 = hf (t n, y n ), k 2 = hf (t n h, y n k 1) y n+1 = y n + k 2 + O(h 3 ) 4th Order Runge Kutta (RK4) - MATLAB ode45 funcction k 1 = hf (t n, y n ), k 2 = hf (t n h, y n k 1) k 3 = hf (t n h, y n k 2), k 4 = hf (t n + h, y n + k 3 ) y n+1 = y n k k k k 4 + O(h 5 )
13 Multiple DOF Systems The Eigenvalue Problem The state-space formulation provides an easy way for solving the equations of motion and is a method that is directly implementable in the time domain. An alternative however exists for solving the problem and for additionally extracting significant information on the properties of the system in the frequency domain, meaning its natural frequencies, modal shapes and damping characteristics.
14 Multiple DOF Systems The Eigenvalue Problem The state-space formulation provides an easy way for solving the equations of motion and is a method that is directly implementable in the time domain. An alternative however exists for solving the problem and for additionally extracting significant information on the properties of the system in the frequency domain, meaning its natural frequencies, modal shapes and damping characteristics.
15 Multiple DOF Systems The Eigenvalue Problem The system of equations (1) is coupled in the sense that each equation involves more than one coordinate (DOF) furthermore, it can be rewritten in matrix - vector notation as: Mẍ + Cẋ + Kx = F(t) k For an SDOF system the natural frequency is obtained as ω n = m. In a similar manner for an n-dof system we consider the following: The undamped - free vibration case can be written as: Mẍ + Kx = 0 n n
16 Multiple DOF Systems The Eigenvalue Problem As for the SDOF system, the general solution for this 2nd order homogeneous (free response) problem is of the following type (harmonic): [ x1 x 2 ] = [ X1 X 2 ] e iωt We can plug this in the homogeneous differential equation: { [ ] [ ]} [ ] [ ω 2 m1 0 k1 + k + 2 k 2 X1 e iωt 0 = 0 m 2 k 2 k 2 X 2 0 ]
17 Multiple DOF Systems The Eigenvalue Problem For x 1 = x 2 = 0 we obtain the trivial solution x(t) = 0. So, instead we demand: { [ ] [ ]} [ ] [ ] ω 2 m1 0 k1 + k + 2 k 2 X1 0 = 0 m 2 k 2 k 2 X 2 0 From det(k Mω 2 ) = 0 we obtain the eigenfrequencies (or natural frequencies) ω 1,2. Then for ω = ω 1 we have that: [ ω 2 1m 1 + (k 1 + k 2 )]X 1 k 2 X 2 = 0 k 2 X 1 + [ ω 2 1m 2 + (k 2 + k 3 )]X 2 = 0 The two equations are linearly dependent, hence we solve one to finally obtain: [ ] [ ] X1 1 k 2 = ω1 2m X 1 X 1 2 φ 1 = 1 τ1 1 + (k 1 + k 2 ) X 1 = X st eigenvector }{{} τ 1
18 Multiple DOF Systems The Eigenvalue Problem Similarly we solve for ω = ω 2 to obtain the 2nd eigenvector: [ ] τ2 φ 2 = X1 2 1 The total solution occurs through the superposition of the two modes: [ ] x1 (t) = C x 2 (t) 1 φ 1 e iω1t + C 1 φ 2 e iω 2t Example: for m 1 = 2m, m 2 + m, k 1 = k 2 = k 1 st mode 2 nd mode 2 2 k 1 F 1 ( t ) k 2 m 1 1 (out of (in phase) phase) Institute of Structural FBD Engineering F 1 ( t ) Identification Methods for Structural F 2Systems ( t ) 1 c 1 x ( t 1 ) c 1
19 Properties of Eigenvectors The previously described process is known as: Solution of the Eigenproblem ω 2 Mφ + Kφ = 0 (eigenproblem) The above homogeneous linear system of equations can only have a solution if the determinant is equal to 0: ω 2 M + K = 0 The solution of the above equation will yield: n eigenvalues ωi 2, i = 1,..., n with 0 ω 1...ω n (the eigenfrequencies) and The solution of ω 2 Mφ = Kφ (eigenproblem) will yield: n eigenvectors (or modal vectors) φ i
20 Properties of Eigenvectors M-Orthonormality From Matrix Properties we know that (AB) T = B T A T, Thus, for two eigenvectors φ n, φ r we obtain: (φ T nkφ r ) T = φ T r K T φ n Since the Stiffness matrix is symmetric, K = K T, hence (φ T nkφ r ) T = φ T r Kφ n The same applies for the Mass matrix, M = M T, yielding: (φ T nmφ r ) T = φ T r Mφ n
21 Properties of Eigenvectors M-Orthonormality Then eigenproblem for vector r can then be written as: (multiply with φ T n on the left) ωr 2 φ T nmφ r = φ T transpose nkφ r ω 2 r φ T r Mφ n = φ T r Kφ n and the eigenproblem for vector n is: (multiply with φ T r on the left) ω 2 nφ T r Mφ n = φ T r Kφ n By subtracting the two previous formulas we obtain: (ω 2 r ω 2 n)φ T r Mφ n = 0 We conclude that for n r φ T r Mφ n = 0, φ T r Mφ n = 0,
22 Properties of Eigenvectors M-Orthonormality We choose φ n, such that: Φ T MΦ = I therefore from MΦΩ 2 = KΦ we have that: Φ T KΦ = Ω 2 The principle of M-orthonormality can then be written as: { 1, n = r φ T r Mφ n = 0, n r { ω 2 φ T r Kφ n = n, n = r 0, n r
23 Mode Superposition Method Matrix Form Defining a matrix Φ whose columns are the eigenvectors φ i and a diagonal matrix Ω 2 which stores the eigenvalues ωi 2 on its diagonal, i.e: ω 1 Φ = [ ] φ 1, φ 2,... φ n ; Ω 2 = ω 2... ; we can write the n solutions to the eigenproblem as: MΦΩ 2 = KΦ ω n
24 Mode Superposition Method Decoupling the Equation of Motion Assume a regular coupled system where we note the displacement vector as U(t) and therefore we can write the following equation of motion: MÜ + C U + KU = F(t) We can now use Φ as a transformation matrix by defining: U(t) = ΦX(t) Which leads to the transformed equilibrium equation: Ẍ(t) + Φ T CΦẊ(t) + Ω 2 X(t) = Φ T R(t) where using the property of M-orthonormality, the initial conditions will be: 0 X = Φ T M 0 U; 0 Ẋ = Φ T M 0 U
25 Dynamic Response with Damping Neglected If we neglect the velocity dependent damping effects the equilibrium equation reduces to: Ẍ (t) + Ω 2 X(t) = Φ T R(t) i.e n individual equations of the form (since Ω 2 is diagonal): ẍ i (t) + ωi 2x } i(t) = r i (t) r i (t) = φ T i R(t) i = 1, 2,..., n with 0 x i = φ T i M 0 U; 0ẋ i = φ T i M 0 U
26 Dynamic Response with Damping Neglected Individual SDOF system response Each equation of the previous system describes a single degree of freedom system with unit mass and stiffness ωi 2. The solution to this equation for a random input excitation can either be obtained by using the SDOF integration methods (like direct integration) or by using the Duhamel Integral: x i (t) = 1 t r i (τ)sinω i (t τ) dτ + α i sinω i t + β i cosω i t ω i 0 where α i, β i are determined from the initial conditions. Therefore, the SDOF response id owed to two contributions A dynamic (steady - state) response obtained by multiplying the static response by a dynamic load factor (this is the particular solution of the governing differential equation), and An additional dynamic response called the transient response
27 Dynamic Response with Damping Neglected Complete Response The solution of all n SDOF equations are calculated and the finite element nodal point displacements are obtained by superposition of the response in each mode: U(t) = ΦX(t) U(t) = n φ i x i (t) i=1 Therefore the solution scheme is: Solve for the eigenvalues and eigenvectors of problem Solve for the response of the decoupled SDOF equations Use superposition to find the total response.
28 Superposition Principle An alternative View Superposition A dynamic load can be designed as Fourier series of harmonic sine and cosine contributions The total solution of such a problem is equal to the superposition of solution of the Fourier terms. A dynamic load can be designed as Fourier series The total solution of such a problem is equal to the superposition of solution of the Fourier terms. U(t) = ΦX(t) U(t) = n φ i x i (t) i= Mode Superposition 11
29 Swiss Federal Institute of Technology Example - 2 DOF system Direct Integration ti Methods Calculate the displacement response of the system Pag k 1 = 4 U, U, U m = 2 R = k 2 = m 2 = 1 R 2 = U U1 0 + = 0 1 U U 2 10 U, U, U k 3 = 2 Method of Finite Elements II
30 Example - 2 DOF system Eigenproblem Setup ω 2 nmφ n = Kφ n ω 2 n Eigenvalue Calculation [ ] = [ ] φ n K ω 2 n M 6 2ω 2 = 0 n ωn 2 = 2ω4 n 14ωn = 0 ω1 2 = Eigenvector Calculation = 2, ω 2 2 = = 5 ( K ω 2 1 M ) φ 1 = 0 [ ] [ φ11 φ 12 ] = 0 [ φ11 φ 12 ] [ 1 = λ 1 ]
31 Example - 2 DOF system Use of M - orthonormality φ T 1Mφ 1 = 1 λ 2 (m 1 + m 2 ) = 1 λ = 1 3 Hence, φ 1 = Similarly, φ 2 = Transformed Equilibrium Equation [ 2, 0 Ẍ(t) Ẍ(t) + Ω 2 X(t) = Φ T R(t) ] X(t) = [ 0 10 ]
32 Example - 2 DOF system Decoupled SDOF equations Initial Conditions ẍ 1 + 2x 1 = 10 3 U(0) = 0 U(0) = 0 and 2 ẍ 2 + 5x 2 = 10 3 x i (0) = φ T i MU(0) ẋ i (0) = φ T i M U(0) x 1 (0) = 0, ẋ 1 (0) = 0, x 2 (0) = 0, ẋ 2 (0) = 0 Then the exact solutions to the ODEs are: x 1 = 5 (1 cos 2 2t) x 2 = ( 1 + cos 5t) And the total solution is: U(t) = n φ i x i (t) i=1
33 Example - 2 DOF system Note In this case we dealt with a 2DOF example and used an analytical solution to solve the two decoupled SDOF system equations In practice for multi degree of freedom systems, higher modes are neglected and instead of an analytical solution we use a numerical scheme in order to solve for each SDOF ODE Neglecting higher modes offers the tremendous advantage of obtaining a Reduced Order System (ROM).
34 The effect of Damping Problems with neglected dampi The presence of damping reduces the dynamic load factor (which then cannot be infinite) and damps out the transient response The response in the modes with ˆω large is ω negligible For ˆω close to zero the ω system follows the loads statically Effectively only the first p modes need to be used p n, in order to obtain a good approximate solution Mode Superposition
35 Analysis Including Damping Transformed Equilibrium Equation Ẍ(t) + Φ T CΦẊ(t) + Ω 2 X(t) = Φ T R(t) Proportional Damping Assumption φ T i Cφ j = 2ω i ξ i δ ij (1) where ξ i is a modal damping parameter and δ ij is the Kronecker delta (δ ij = 1 for i = j, δ ij = 0 for i j)
36 Analysis Including Damping Therefore, we still end up with Decoupled SDOF equations for each x i : ẍ i (t) + 2ω i ξ i ẋ i (t) + ω 2 i x i (t) = r i (t) with the Duhamel Integral now being: x i (t) = 1 ω t { r i (τ) e ξ i ω i (t τ) sin ω i (t τ) dτ i 0 } + e ξ i ω i t (α i sin ω i t + β i cos ω i t) where ω i = ω i 1 ξ 2 i
37 Analysis Including Damping Rayleigh Damping If there are only two different damping ratios ξ i, Damping can be used: i = 1, 2 Rayleigh C = αm + βk (2) Eqns (1), (2) now yield: φ T i (αm + βk)φ i = 2ω i ξ i α + βω 2 i = 2ω i ξ i The 2 2 system can be solved to obtain α, β. In actual analysis it may well be that the damping ratios are known for many more than two frequencies. In that case two average values say ξ 1, ξ 2 are used to evaluate α, β.
38 Example - Damping for an MDOF system Assume that the approximate damping to be specified for a multiple degree of freedom system is as follows:
39 Example - Damping for an MDOF system Damping as function of frequency A problem with Rayleigh damping is that higher modes are much more damped than lower modes. However, in general Rayleigh damping provides a good approximation
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 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 informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture 6-5 November, 2015
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Lecture 6-5 November, 015 Institute of Structural Engineering Method of Finite Elements II 1 Introduction
More informationAA242B: MECHANICAL VIBRATIONS
AA242B: MECHANICAL VIBRATIONS 1 / 50 AA242B: MECHANICAL VIBRATIONS Undamped Vibrations of n-dof Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations:
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 informationFREE VIBRATION RESPONSE OF UNDAMPED SYSTEMS
Lecture Notes: STRUCTURAL DYNAMICS / FALL 2011 / Page: 1 FREE VIBRATION RESPONSE OF UNDAMPED SYSTEMS : : 0, 0 As demonstrated previously, the above Equation of Motion (free-vibration equation) has a solution
More information3. Mathematical Properties of MDOF Systems
3. Mathematical Properties of MDOF Systems 3.1 The Generalized Eigenvalue Problem Recall that the natural frequencies ω and modes a are found from [ - ω 2 M + K ] a = 0 or K a = ω 2 M a Where M and K are
More information1. Multiple Degree-of-Freedom (MDOF) Systems: Introduction
1. Multiple Degree-of-Freedom (MDOF) Systems: Introduction Lesson Objectives: 1) List examples of MDOF structural systems and state assumptions of the idealizations. 2) Formulate the equation of motion
More informationIdentification Methods for Structural Systems
Prof. Dr. Eleni Chatzi System Stability Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can be defined from
More informationAA 242B / ME 242B: Mechanical Vibrations (Spring 2016)
AA 242B / ME 242B: Mechanical Vibrations (Spring 206) Solution of Homework #3 Control Tab Figure : Schematic for the control tab. Inadequacy of a static-test A static-test for measuring θ would ideally
More informationDr.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 informationReduction in number of dofs
Reduction in number of dofs Reduction in the number of dof to represent a structure reduces the size of matrices and, hence, computational cost. Because a subset of the original dof represent the whole
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 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 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 informationMultiple Degree of Freedom Systems. The Millennium bridge required many degrees of freedom to model and design with.
Multiple Degree of Freedom Systems The Millennium bridge required many degrees of freedom to model and design with. The first step in analyzing multiple degrees of freedom (DOF) is to look at DOF DOF:
More informationTransient Response Analysis of Structural Systems
. 21 Transient Response Analysis of Structural Systems 21 1 Chapter 21: TRANSIENT RESPONSE ANALYSIS OF STRUCTURAL SYSTEMS 21 2 21.1 MODAL APPROACH TO TRANSIENT ANALYSIS Consider the following large-order
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 informationThis work is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 4.0 License.
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 015 14. Oscillations Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License This wor
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 informationSolution of Vibration and Transient Problems
. 9 Solution of Vibration and Transient Problems 9 Chapter 9: SOLUTION OF VIBRATION AND TRANSIENT PROBLEMS 9 9. MODAL APPROACH TO TRANSIENT ANALYSIS Consider the following large-order finite element model
More informationSTRUCTURAL DYNAMICS BASICS:
BASICS: STRUCTURAL DYNAMICS Real-life structures are subjected to loads which vary with time Except self weight of the structure, all other loads vary with time In many cases, this variation of the load
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 information22.2. Applications of Eigenvalues and Eigenvectors. Introduction. Prerequisites. Learning Outcomes
Applications of Eigenvalues and Eigenvectors 22.2 Introduction Many applications of matrices in both engineering and science utilize eigenvalues and, sometimes, eigenvectors. Control theory, vibration
More informationTHE subject of the analysis is system composed by
MECHANICAL VIBRATION ASSIGNEMENT 1 On 3 DOF system identification Diego Zenari, 182160, M.Sc Mechatronics engineering Abstract The present investigation carries out several analyses on a 3-DOF system.
More informationDynamic Response of Structures With Frequency Dependent Damping
Dynamic Response of Structures With Frequency Dependent Damping Blanca Pascual & S Adhikari School of Engineering, Swansea University, Swansea, UK Email: S.Adhikari@swansea.ac.uk URL: http://engweb.swan.ac.uk/
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014
Prof. Dr. Eleni Chatzi System Stability - 26 March, 24 Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can
More informationL = 1 2 a(q) q2 V (q).
Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More informationUsing Simulink to analyze 2 degrees of freedom system
Using Simulink to analyze 2 degrees of freedom system Nasser M. Abbasi Spring 29 page compiled on June 29, 25 at 4:2pm Abstract A two degrees of freedom system consisting of two masses connected by springs
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 informationAppendix A Equations of Motion in the Configuration and State Spaces
Appendix A Equations of Motion in the Configuration and State Spaces A.1 Discrete Linear Systems A.1.1 Configuration Space Consider a system with a single degree of freedom and assume that the equation
More informationAdvanced Vibrations. Elements of Analytical Dynamics. By: H. Ahmadian Lecture One
Advanced Vibrations Lecture One Elements of Analytical Dynamics By: H. Ahmadian ahmadian@iust.ac.ir Elements of Analytical Dynamics Newton's laws were formulated for a single particle Can be extended to
More informationStructural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian
Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:
More informationLecture 27: Structural Dynamics - Beams.
Chapter #16: Structural Dynamics and Time Dependent Heat Transfer. Lectures #1-6 have discussed only steady systems. There has been no time dependence in any problems. We will investigate beam dynamics
More informationMathematics Qualifying Exam Study Material
Mathematics Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering mathematics topics. These topics are listed below for clarification. Not all instructors
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 informationCollocated versus non-collocated control [H04Q7]
Collocated versus non-collocated control [H04Q7] Jan Swevers September 2008 0-0 Contents Some concepts of structural dynamics Collocated versus non-collocated control Summary This lecture is based on parts
More information! 4 4! o! +! h 4 o=0! ±= ± p i And back-substituting into the linear equations gave us the ratios of the amplitudes of oscillation:.»» = A p e i! +t»»
Topic 6: Coupled Oscillators and Normal Modes Reading assignment: Hand and Finch Chapter 9 We are going to be considering the general case of a system with N degrees of freedome close to one of its stable
More informationIntroduction to Vibration. Mike Brennan UNESP, Ilha Solteira São Paulo Brazil
Introduction to Vibration Mike Brennan UNESP, Ilha Solteira São Paulo Brazil Vibration Most vibrations are undesirable, but there are many instances where vibrations are useful Ultrasonic (very high
More informationCE 6701 Structural Dynamics and Earthquake Engineering Dr. P. Venkateswara Rao
CE 6701 Structural Dynamics and Earthquake Engineering Dr. P. Venkateswara Rao Associate Professor Dept. of Civil Engineering SVCE, Sriperumbudur Difference between static loading and dynamic loading Degree
More informationChapter 5 Design. D. J. Inman 1/51 Mechanical Engineering at Virginia Tech
Chapter 5 Design Acceptable vibration levels (ISO) Vibration isolation Vibration absorbers Effects of damping in absorbers Optimization Viscoelastic damping treatments Critical Speeds Design for vibration
More informationSHOCK RESPONSE OF MULTI-DEGREE-OF-FREEDOM SYSTEMS Revision F By Tom Irvine May 24, 2010
SHOCK RESPONSE OF MULTI-DEGREE-OF-FREEDOM SYSTEMS Revision F By Tom Irvine Email: tomirvine@aol.com May 4, 010 Introduction The primary purpose of this tutorial is to present the Modal Transient method
More informationProblem 1: Find the Equation of Motion from the static equilibrium position for the following systems: 1) Assumptions
Problem 1: Find the Equation of Motion from the static equilibrium position for the following systems: 1) Assumptions k 2 Wheels roll without friction k 1 Motion will not cause block to hit the supports
More informationME 375 EXAM #1 Friday, March 13, 2015 SOLUTION
ME 375 EXAM #1 Friday, March 13, 2015 SOLUTION PROBLEM 1 A system is made up of a homogeneous disk (of mass m and outer radius R), particle A (of mass m) and particle B (of mass m). The disk is pinned
More informationThe Finite Element Method
Page he Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. o. Dr. Michael Havbro bo Faber Dr. Nebojsa Mojsilovic Swiss Federal Institute of EH Zurich, Switzerland Contents of
More informationDamping Matrix. Donkey2Ft
1 Damping Matrix DonkeyFt Damping in a single-degree-of-freedom (SDOF) system is well studied. Whether the system is under-damped, over-damped, or critically damped is well known. For an under-damped system,
More informationStochastic Dynamics of SDOF Systems (cont.).
Outline of Stochastic Dynamics of SDOF Systems (cont.). Weakly Stationary Response Processes. Equivalent White Noise Approximations. Gaussian Response Processes as Conditional Normal Distributions. Stochastic
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 information1 Review of simple harmonic oscillator
MATHEMATICS 7302 (Analytical Dynamics YEAR 2017 2018, TERM 2 HANDOUT #8: COUPLED OSCILLATIONS AND NORMAL MODES 1 Review of simple harmonic oscillator In MATH 1301/1302 you studied the simple harmonic oscillator:
More informationMATH 2250 Final Exam Solutions
MATH 225 Final Exam Solutions Tuesday, April 29, 28, 6: 8:PM Write your name and ID number at the top of this page. Show all your work. You may refer to one double-sided sheet of notes during the exam
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 informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY SPRING SEMESTER 207 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationOutline. Structural Matrices. Giacomo Boffi. Introductory Remarks. Structural Matrices. Evaluation of Structural Matrices
Outline in MDOF Systems Dipartimento di Ingegneria Civile e Ambientale, Politecnico di Milano May 8, 014 Additional Today we will study the properties of structural matrices, that is the operators that
More informationAppendix C. Modal Analysis of a Uniform Cantilever with a Tip Mass. C.1 Transverse Vibrations. Boundary-Value Problem
Appendix C Modal Analysis of a Uniform Cantilever with a Tip Mass C.1 Transverse Vibrations The following analytical modal analysis is given for the linear transverse vibrations of an undamped Euler Bernoulli
More informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear
More informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture 9-23 April, 2013
Prof. Dr. Eleni Chatzi Lecture 9-23 April, 2013 Identification Methods The work is done either in the frequency domain - Modal ID methods using the Frequency Response Function (FRF) information or in the
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 informationEQUIVALENT SINGLE-DEGREE-OF-FREEDOM SYSTEM AND FREE VIBRATION
1 EQUIVALENT SINGLE-DEGREE-OF-FREEDOM SYSTEM AND FREE VIBRATION The course on Mechanical Vibration is an important part of the Mechanical Engineering undergraduate curriculum. It is necessary for the development
More informationPeriodic functions: simple harmonic oscillator
Periodic functions: simple harmonic oscillator Recall the simple harmonic oscillator (e.g. mass-spring system) d 2 y dt 2 + ω2 0y = 0 Solution can be written in various ways: y(t) = Ae iω 0t y(t) = A cos
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 informationMore chapter 3...linear dependence and independence... vectors
More chapter 3...linear dependence and independence... vectors It is important to determine if a set of vectors is linearly dependent or independent Consider a set of vectors A, B, and C. If we can find
More informationIn which we count degrees of freedom and find the normal modes of a mess o masses and springs, which is a lovely model of a solid.
Coupled Oscillators Wednesday, 26 October 2011 In which we count degrees of freedom and find the normal modes of a mess o masses and springs, which is a lovely model of a solid. Physics 111 θ 1 l k θ 2
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 informationKeble College - Hilary 2014 CP3&4: Mathematical methods I&II Tutorial 5 - Waves and normal modes II
Tomi Johnson 1 Keble College - Hilary 2014 CP3&4: Mathematical methods I&II Tutorial 5 - Waves and normal modes II Prepare full solutions to the problems with a self assessment of your progress on a cover
More informationConsistency and Convergence
Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained
More informationCoupled Oscillators Monday, 29 October 2012 normal mode Figure 1:
Coupled Oscillators Monday, 29 October 2012 In which we count degrees of freedom and find the normal modes of a mess o masses and springs, which is a lovely model of a solid. Physics 111 θ 1 l m k θ 2
More informationTWO-STAGE ISOLATION FOR HARMONIC BASE EXCITATION Revision A. By Tom Irvine February 25, 2008
TWO-STAGE ISOLATION FOR HARMONIC BASE EXCITATION Revision A By Tom Irvine Email: tomirvine@aol.com February 5, 008 Introduction Consider a base plate mass m and an avionics mass m modeled as two-degree-of-freedom.
More informationME8230 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 informationANNEX A: ANALYSIS METHODOLOGIES
ANNEX A: ANALYSIS METHODOLOGIES A.1 Introduction Before discussing supplemental damping devices, this annex provides a brief review of the seismic analysis methods used in the optimization algorithms considered
More informationChapter 4 Analysis of a cantilever
Chapter 4 Analysis of a cantilever Before a complex structure is studied performing a seismic analysis, the behaviour of simpler ones should be fully understood. To achieve this knowledge we will start
More informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationFinal Exam Solution Dynamics :45 12:15. Problem 1 Bateau
Final Exam Solution Dynamics 2 191157140 31-01-2013 8:45 12:15 Problem 1 Bateau Bateau is a trapeze act by Cirque du Soleil in which artists perform aerial maneuvers on a boat shaped structure. The boat
More informationThe Phasor Analysis Method For Harmonically Forced Linear Systems
The Phasor Analysis Method For Harmonically Forced Linear Systems Daniel S. Stutts, Ph.D. April 4, 1999 Revised: 10-15-010, 9-1-011 1 Introduction One of the most common tasks in vibration analysis is
More informationReview: 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 informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More informationdx n a 1(x) dy
HIGHER ORDER DIFFERENTIAL EQUATIONS Theory of linear equations Initial-value and boundary-value problem nth-order initial value problem is Solve: a n (x) dn y dx n + a n 1(x) dn 1 y dx n 1 +... + a 1(x)
More informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationResponse to Periodic and Non-periodic Loadings. Giacomo Boffi. March 25, 2014
Periodic and Non-periodic Dipartimento di Ingegneria Civile e Ambientale, Politecnico di Milano March 25, 2014 Outline Introduction Fourier Series Representation Fourier Series of the Response Introduction
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 informationComputational Physics (6810): Session 8
Computational Physics (6810): Session 8 Dick Furnstahl Nuclear Theory Group OSU Physics Department February 24, 2014 Differential equation solving Session 7 Preview Session 8 Stuff Solving differential
More informationMatrix Iteration. Giacomo Boffi.
http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April 12, 2016 Outline Second -Ritz Method Dynamic analysis of MDOF
More informationStructural Dynamics A Graduate Course in Aerospace Engineering
Structural Dynamics A Graduate Course in Aerospace Engineering By: H. Ahmadian ahmadian@iust.ac.ir The Science and Art of Structural Dynamics What do all the followings have in common? > A sport-utility
More informationTruncation Errors Numerical Integration Multiple Support Excitation
Errors Numerical Integration Multiple Support Excitation http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April 10,
More information2C9 Design for seismic and climate changes. Jiří Máca
2C9 Design for seismic and climate changes Jiří Máca List of lectures 1. Elements of seismology and seismicity I 2. Elements of seismology and seismicity II 3. Dynamic analysis of single-degree-of-freedom
More informationCHAPTER 6 STATE SPACE: FREQUENCY RESPONSE, TIME DOMAIN
CHAPTER 6 STATE SPACE: FREQUENCY RESPONSE, TIME DOMAIN 6. Introduction Frequency Response This chapter will begin with the state space form of the equations of motion. We will use Laplace transforms to
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 informationL = 1 2 a(q) q2 V (q).
Physics 3550 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6, 11.1 11.3 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More informationProblem 1: Lagrangians and Conserved Quantities. Consider the following action for a particle of mass m moving in one dimension
105A Practice Final Solutions March 13, 01 William Kelly Problem 1: Lagrangians and Conserved Quantities Consider the following action for a particle of mass m moving in one dimension S = dtl = mc dt 1
More informationDesign of Structures for Earthquake Resistance
NATIONAL TECHNICAL UNIVERSITY OF ATHENS Design of Structures for Earthquake Resistance Basic principles Ioannis N. Psycharis Lecture 3 MDOF systems Equation of motion M u + C u + K u = M r x g(t) where:
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 informationLinear Algebra Review (Course Notes for Math 308H - Spring 2016)
Linear Algebra Review (Course Notes for Math 308H - Spring 2016) Dr. Michael S. Pilant February 12, 2016 1 Background: We begin with one of the most fundamental notions in R 2, distance. Letting (x 1,
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 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 informationApplied Math for Engineers
Applied Math for Engineers Ming Zhong Lecture 15 March 28, 2018 Ming Zhong (JHU) AMS Spring 2018 1 / 28 Recap Table of Contents 1 Recap 2 Numerical ODEs: Single Step Methods 3 Multistep Methods 4 Method
More informationAutonomous system = system without inputs
Autonomous system = system without inputs State space representation B(A,C) = {y there is x, such that σx = Ax, y = Cx } x is the state, n := dim(x) is the state dimension, y is the output Polynomial representation
More informationTime integration and the Trefftz Method, Part II Second-order and hyperbolic problems
ime integration and the refftz Method, Part II Second-order and hyperbolic problems JA eixeira de Freitas Departamento de Engenharia Civil e Arquitectura, Instituto Superior écnico echnical University
More informationA Guide to linear dynamic analysis with Damping
A Guide to linear dynamic analysis with Damping This guide starts from the applications of linear dynamic response and its role in FEA simulation. Fundamental concepts and principles will be introduced
More informationCIVL 8/7117 Chapter 12 - Structural Dynamics 1/75. To discuss the dynamics of a single-degree-of freedom springmass
CIV 8/77 Chapter - /75 Introduction To discuss the dynamics of a single-degree-of freedom springmass system. To derive the finite element equations for the time-dependent stress analysis of the one-dimensional
More informationChapter 2 Basics of Vibration Dynamics
Chapter 2 Basics of Vibration Dynamics Vibrations are mechanical oscillations about an equilibrium position. There are cases when vibrations are desirable, such as in certain types of machine tools or
More information