Dynamic Response of Structures With Frequency Dependent Damping
|
|
- Charlene Dora Conley
- 6 years ago
- Views:
Transcription
1 Dynamic Response of Structures With Frequency Dependent Damping Blanca Pascual & S Adhikari School of Engineering, Swansea University, Swansea, UK S.Adhikari@swansea.ac.uk URL: adhikaris Frequency Dependent Damping p.1/35
2 Outline of the presentation Damping models in structural dynamics Review of current approaches Dynamic response of frequency dependent damped systems Non-linear eigenvalue problem SDOF systems: real & complex solutions MDOF systems: real & complex solutions Conclusions & discussions Frequency Dependent Damping p.2/35
3 Damping models In general a physically realistic model of damping may not be a viscous damping model. Damping models in which the dissipative forces depend on any quantity other than the instantaneous generalized velocities are non-viscous damping models. Possibly the most general way to model damping within the linear range is to use non-viscous damping models which depend on the past history of motion via convolution integrals over kernel functions. Frequency Dependent Damping p.3/35
4 Equation of motion The equations of motion of a N-DOF linear system: Mü(t) + t 0 G(t τ) u(τ) dτ + Ku(t) = f(t) (1) together with the initial conditions u(t = 0) = u 0 R N and u(t = 0) = u 0 R N. (2) u(t): displacement vector, f(t): forcing vector, M, K: mass and stiffness matrices. In the limit when G(t τ) = Cδ(t τ), where δ(t) is the Diracdelta function, this reduces to viscous damping. Frequency Dependent Damping p.4/35
5 Damping functions - 1 Model Damping function Author and year of publication Number 1 G(s) = n k=1 a k s s + b k Biot[1] G(s) = E 1s α E 0 bs β 1 + bs β (0 < α, β < 1) Bagley and Torvik[2] sg(s) = G [1 + ] k α s 2 + 2ξ k ω k s k s 2 + 2ξ k ω k s + ωk 2 Golla and Hughes[3] and McTavish and Hughes[4] G(s) = 1 + n k s k=1 Lesieutre and Mingori[5] s + β k 5 G(s) = c 1 e st 0 Adhikari[6] st 0 6 G(s) = c 1 + 2(st 0 /π) 2 e st0 st (st 0 /π) 2 Adhikari[6] [ ( )] 7 G(s) = c e s2 /4µ s 1 erf 2 Adhikari and Woodhouse[7] µ Some damping functions in the Laplace domain. Frequency Dependent Damping p.5/35
6 Damping functions - 2 We use a damping model for which the kernel function matrix: n G(t) = µ k e µ kt C k (3) k=1 The constants µ k R + are known as the relaxation parameters and n denotes the number relaxation parameters. When µ k, k this reduces to the viscous damping model: n C = C k. (4) k=1 Frequency Dependent Damping p.6/35
7 State-space Approach Eq (1) can be transformed to B ż(t) = A z(t) + r(t) (5) where the m (m = 2N + nn) dimensional matrices and vectors are: B = n C k M C 1 /µ 1 C n /µ n k=1 M O O O O C 1 /µ 1 O C 1 /µ 2 1 O O, r(t) =.... O O O C n /µ n O O O C n /µ 2 n f(t) (6) A = K O O O O u(t) O M O O O v(t) O O C 1 /µ 1 O O, z(t) = y 1 (t)... O O O O. O O O O C n /µ n y n (t) (7) Frequency Dependent Damping p.7/35
8 Main issues The main reasons for not using a frequency dependent non-viscous damping model include, but not limited to: although exact in nature, the state-space approach usually needed for this type of damped systems is computationally very intensive for real-life systems; the physical insights offered by methods in the original space (eg, the modal analysis) is lost in a state-space based approach the experimental identification of the parameters of a frequency dependent damping model is difficult. Frequency Dependent Damping p.8/35
9 Dynamic Response - 1 Taking the Laplace transform of equation (1) and considering the initial conditions in (2) we have s 2 M q smq 0 M q 0 + sg(s) q G(s)q 0 + K q = f(s) or D(s) q = f(s) + M q 0 + [sm + G(s)]q 0. The dynamic stiffness matrix is defined as D(s) = s 2 M + sg(s) + K C N N. (8) The inverse of the dynamics stiffness matrix, known as the transfer function matrix, is given by H(s) = D 1 (s) C N N. (9) Frequency Dependent Damping p.9/35
10 Dynamic Response - 2 Using the residue-calculus the transfer function matrix can be expressed like a viscously damped system as H(s) = m j=1 R j s s j ; R j = res s=s j [H(s)] = z T j z j z T j (10) D(s j ) s j z j where m is the number of non-zero eigenvalues (order) of the system, s j and z j are respectively the eigenvalues and eigenvectors of the system, which are solutions of the non-linear eigenvalue problem [s 2 jm + s j G(s j ) + K]z j = 0, for j = 1,,m (11) Frequency Dependent Damping p.10/35
11 Dynamic Response - 3 The expression of H(s) allows the response to be expressed as modal summation as q(s) = where m j=1 z T f(s) j + z T j M q 0 + sz T j Mq 0 + z T j G(s)q 0 (s) γ j z j s s j γ j = z T j 1 (12) D(s j ) s j z j. (13) We aim to derive the eigensolutions in N -space by solving the non-linear eigenvalue problem. Frequency Dependent Damping p.11/35
12 Non-linear Eigenvalue Problem The eigenvalue problem associated with a linear system with exponential damping model: [ ] n s 2 µ k jm + s j C k + K z j = 0, s j + µ k k=1 Two types of eigensolutions: 2N complex conjugate solutions - underdamped/vibrating modes for j = 1,,m. p real solutions [p = n k=1 rank (C k)] - overdamped modes (14) Frequency Dependent Damping p.12/35
13 Non-linear Eigenvalue Problem The following four cases are considered: single-degree-of-freedom system with single exponential kernel (N = 1,n = 1) single-degree-of-freedom system with multiple exponential kernels (N = 1,n > 1) multiple-degree-of-freedom system with single exponential kernel (N > 1,n = 1) multiple-degree-of-freedom system with multiple exponential kernels (N > 1,n > 1) Frequency Dependent Damping p.13/35
14 SDOF systems Computational cost and other relevant issues identified before do not strictly affect the eigenvalue problem of a single-degree-of-freedom system (SDOF) with exponential damping. The main reason for considering a SDOF system is that in many cases the underlying approximation method can be extended to MDOF systems in a relatively straight-forward manner. Frequency Dependent Damping p.14/35
15 Complex-conjugate solutions The main motivation of the approximations is that the approximate solution can be constructed from the solution of equivalent viscously damped system. The solution of equivalent viscously damped system can in turn be expressed in terms of the undamped eigensolutions. Combining these together, one can therefore obtain the eigensolutions of frequency-dependent systems by simple post-processing of undamped solutions only. Frequency Dependent Damping p.15/35
16 Complex-conjugate solutions The eigenvalues of the equivalent viscously damped system: s 0 = ζ n ω n ± iω n 1 ζ 2 n ζ n ω n ± iω n (15) ω n = k u /m u and ζ n = c/2 k u m u. Viscous damped system is a special case when the function g(s) is replaced by g(s ). For that case this solution would have been the exact solution of the characteristic equation. The difference between the viscous solution and the true solution is essentially arising due to the varying nature of the function g(s). Frequency Dependent Damping p.16/35
17 Complex-conjugate solutions The central idea here is that the actual solution can be obtained by expanding the solution in a Taylor series around s 0. We assume s = s 0 + δ, (δ is small). Substituting this into the characteristic equation we have (s 0 + δ) 2 m u + (s 0 + δ)g(s 0 + δ) + k u = 0. (16) First-order approximation δ (1) = s 0(s 0 m u + g(s 0 )) + k u s 0 (2m u + g (s 0 )) + g(s 0 ) (17) Frequency Dependent Damping p.17/35
18 Complex-conjugate solutions Second-order approximation: δ (2) = B B 2 4AC 2A (18) where A = (m u + g (s 0 ) s 0 + g (s 0 )) 2! (19) B = (2m u s 0 + s 0 g (s 0 ) + g(s 0 )) (20) and C = (s 2 0m u + s 0 g(s 0 ) + k). (21) g (s 0 ) and g (s 0 ) are respectively the first and second order derivative of g(s) evaluated at s = s 0. Frequency Dependent Damping p.18/35
19 Real solutions When N = 1,n = 1 the eigenvalue equation: s 2 m u + sg(s) + k u = 0 where g(s) = µ c. (22) s + µ While the complex-conjugate solution can be expected to be close to the solution of the equivalent viscously damped system, no such analogy can be made for the real solution as the equivalent viscously damped system doesn t have one. Rewrite the characteristic equation as (s 2 m u + k u )(µ + s) + scµ = 0. (23) Frequency Dependent Damping p.19/35
20 Real solutions Consider that the damping is small so that scµ 0. Since (s 2 m u + k u ) 0 as we are considering the real solution only, the first guess is obtained as µ + s + 0 = 0 or s 0 = µ. We take the first approximation of the real root as s = s 0 + = µ + Substituting into the characteristic equation and neglecting all the higher-order terms: µ 2 c µ 2 m u + k u + µc (24) Frequency Dependent Damping p.20/35
21 Real solutions: General case The characteristic equation: s 2 m u + s n k=1 µ k s + µ k c k + k u = 0. (25) This a polynomial in s of order (n + 2) - it has (n + 2) roots [2 complex conjugate and n real]. Multiplying the characteristic equation by the product n j (s + µ j): (s 2 m u + k u ) n (s + µ j ) + s j=1 n n (µ k c k (µ j + s)) = 0. (26) k=1 j=1 j k Frequency Dependent Damping p.21/35
22 Real solutions: General case Use the approximation s k = µ k + k, k = 1, 2, n. Substituting into the characteristic equation and retaining only the first-order terms in k, after some simplifications: k where θ k = c k µ 2 k n (µ j µ k ) j=1 j k (µ 2 k m + k) n µ k c k nj=1 j k r=1 n k nm=1 m j m k (µ r µ k ) + µ k θ k for k = 1, 2,, n. (µ m µ k ) n r=1 r k µ r c r nj=1 j n j k (µ j µ k ) + n j=1 j k c k (µ j µ k ) Frequency Dependent Damping p.22/35
23 Numerical Results SDOF system with eight exponential kernels: m = 1 kg, k = 2 N/m, µ k for k = 1, 2.., 8 are , , , , , , , Frequency Dependent Damping p.23/35
24 Numerical Results µ Exact solution (statespace) Proposed solution approximate Percentage error Complex solution Conjugate ± i ±1.4718i ±0i Exact and approximate eigenvalues of the SDOF system. Frequency Dependent Damping p.24/35
25 Frequency response function 10 1 exact solution approximation 10 2 log(abs(h)) ω (rad/s) Frequency response function of the SDOF system obtained using exact and approximate eigenvalues. Frequency Dependent Damping p.25/35
26 Complex solutions: MDOF Complex modes can be expanded as a complex linear combination of undamped modes z j = N l=1 α(j) l x l Substituting in to the eigenvalue equation (11): N l=1 s 2 jα (j) l Mx l + s j α (j) l G(s j )x l + α (j) l Kx l = 0. (27) Premultiplying by x T k and using the mass-orthogonality property of the undamped eigenvectors: s 2 jα (j) k +s j N l=1 α (j) l G kl(s j )+ωkα 2 (j) k = 0, G kl(s j ) = x T kg(s j )x l (28) Frequency Dependent Damping p.26/35
27 Complex solutions: MDOF Considering the j-th set of equation and neglecting the second-order terms involving α (j) k and G kl (s j), k l: s 2 j + s j G jj(s j ) + ω 2 j 0 (29) Similar to the SDOF case (replace m u by 1, k u by ω 2 j and g(s) by G jj(s)). To obtain the eigenvectors rewrite Eq. (28) for j k as s 2 j α(j) k + s j G kj (s j) + α (j) k G kk (s j) + N l k j l G kl (s j) +ωk 2 α(j) k = 0, α (j) k = 1,, N; j. (30) Frequency Dependent Damping p.27/35
28 Complex solutions: MDOF Retaining only the product terms α (j) l G kl : z j x j N k=1 k j s j G kj (s j)x k ω 2 k + s2 j + s jg kk (s j). (31) A second-order expressions is given in the paper. The results derived here are based on small non-proportional damping. Frequency Dependent Damping p.28/35
29 Real solutions: MDOF For systems with single exponential kernel assume s l = µ + l (32) Substituting in the approximate characteristic Eq (29) l µ 2 C ll ; l = 1, 2,,N (33) µ 2 + ωl 2 + µc ll Frequency Dependent Damping p.29/35
30 Real solutions: MDOF Assuming all coefficient matrices are of full rank, for systems with n kernels there are in general nn number of purely real eigenvalues. The approximate eigenvalues can be obtained as s lk = µ k + lk (34) where lk C kll µ 2 k (µ 2 k + ω2 l ) n θ lk = µ k C kll nj=1 j k r=1 n k nm=1 n (µ j µ k ) j=1 j k m j m k and C kll = x T l C kx l. (µ r µ k ) + µ k θ lk k = 1, 2,, n; l = 1, 2,, N (µ m µ k ) n r=1 r k µ r c r nj=1 j n j k (µ j µ k )+ n j=1 C kll (µ j µ k ) j k Frequency Dependent Damping p.30/35
31 Numerical example We consider a three degree-of-freedom system: M = m u m u 0, K = 2k u k u 0 k u 2k u k u (35) 0 0 m u 0 k u 2k u G(s) = C 6 k=1 µ k s + µ k, where C = (36) m u = 1 kg, k u = 1 N/m and µ k are , , , , , Schaumburg, Illinois, 9 April 2008 Frequency Dependent Damping p.31/35
32 Numerical results Exact solution (state-space) Proposed approximate solution Real solutions Percentage error Frequency Dependent Damping p.32/35
33 Numerical results Exact solution (state-space) Proposed approximate solution Real solutions Complex Conjugate solutions Percentage error ± i ± i ± i ± i ± i ± i ± i ± i ± i Frequency Dependent Damping p.33/35
34 Conclusions - 1 Multiple degree-of-freedom linear systems with frequency depended damping kernels is considered. The transfer function matrix of the system was derived in terms of the eigenvalues and eigenvectors of the second-order system. The response can be expressed as a sum of two parts, one that arises in usual viscously damped systems and the other that occurs due to non-viscous damping. The calculation of the eigensolutions of frequency-depended damped systems requires the solution of a non-linear eigenvalue problem. Frequency Dependent Damping p.34/35
35 Conclusions - 2 Approximate expressions are derived for the complex and real eigenvalues of the SDOF system with single and multiple exponential kernels. These results are extended to MDOF systems. These approximations allow one to obtain the dynamic response of general frequency-depended damped systems by simple post-processing of undamped eigensolutions. The accuracy of the proposed approximations were verified using numerical examples. The complex conjugate eigensolutions turn out to be more accurate compared to the real eigensolutions. Frequency Dependent Damping p.35/35
36 References [1] Biot, M. A., Variational principles in irreversible thermodynamics with application to viscoelasticity, Physical Review, Vol. 97, No. 6, 1955, pp [2] Bagley, R. L. and Torvik, P. J., Fractional calculus a different approach to the analysis of viscoelastically damped structures, AIAA Journal, Vol. 21, No. 5, May 1983, pp [3] Golla, D. F. and Hughes, P. C., Dynamics of viscoelastic structures - a time domain finite element formulation, Transactions of ASME, Journal of Applied Mechanics, Vol. 52, December 1985, pp [4] McTavish, D. J. and Hughes, P. C., Modeling of linear viscoelastic space structures, Transactions of ASME, Journal of Vibration and Acoustics, Vol. 115, January 1993, pp [5] Lesieutre, G. A. and Mingori, D. L., Finite element modeling of frequency-dependent material properties using augmented thermodynamic fields, AIAA Journal of Guidance, Control and Dynamics, Vol. 13, 1990, pp [6] Adhikari, S., Energy Dissipation in Vibrating Structures, Master s thesis, Cambridge University Engineering Department, Cambridge, UK, May 1998, First Year Report. 35-1
37 [7] Adhikari, S. and Woodhouse, J., Identification of damping: part 1, viscous damping, Journal of Sound and Vibration, Vol. 243, No. 1, May 2001, pp
Calculation of Eigensolution Derivatives for Nonviscously Damped Systems Using Nelson s Method
AIAA JOURNAL Vol. 44, No. 8, August 006 Calculation of Eigensolution Derivatives for Nonviscously Damped Systems Using Nelson s Method S. Adhikari and M. I. Friswell University of Bristol, Bristol, England
More informationDAMPING MODELLING AND IDENTIFICATION USING GENERALIZED PROPORTIONAL DAMPING
DAMPING MODELLING AND IDENTIFICATION USING GENERALIZED PROPORTIONAL DAMPING S. Adhikari Department of Aerospace Engineering, University of Bristol, Queens Building, University Walk, Bristol BS8 1TR (U.K.)
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 informationIdentification of Damping Using Proper Orthogonal Decomposition
Identification of Damping Using Proper Orthogonal Decomposition M Khalil, S Adhikari and A Sarkar Department of Aerospace Engineering, University of Bristol, Bristol, U.K. Email: S.Adhikari@bristol.ac.uk
More informationCOMPLEX MODES IN LINEAR STOCHASTIC SYSTEMS
Proceedings of VETOMAC-I October 25-27, 2, Bangalore, INDIA COMPLEX MODES IN LINEAR STOCHASTIC SYSTEMS S. Adhikari Cambridge University Engineering Department Trumpington Street Cambridge CB2 1PZ (U.K.)
More informationDamping Modelling and Identification Using Generalized Proportional Damping
Damping Modelling and Identification Using Generalized Proportional Damping S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. Email: S.Adhikari@bristol.ac.uk Generalized
More informationTOWARDS IDENTIFICATION OF A GENERAL MODEL OF DAMPING
TOWARDS IDENTIFICATION OF A GENERAL MODEL OF DAMPING S. Adhikari and J. Woodhouse Department of Engineering, University of Cambridge, Trumpington Street, Cambridge CB PZ (U.K.) ABSTRACT Characterization
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 informationRandom Eigenvalue Problems Revisited
Random Eigenvalue Problems Revisited S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. Email: S.Adhikari@bristol.ac.uk URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html
More informationThe finite element method, coupled with model updating techniques, allow us to accurately model the
46th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics & Materials Conference, 18-21 Apr 25, Austin, Texas Dynamics of Non-viscously Damped Distributed Parameter Systems Sondipon Adhikari University
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 informationRandom Matrix Eigenvalue Problems in Probabilistic Structural Mechanics
Random Matrix Eigenvalue Problems in Probabilistic Structural Mechanics S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html
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 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 informationME scope Application Note 28
App Note 8 www.vibetech.com 3/7/17 ME scope Application Note 8 Mathematics of a Mass-Spring-Damper System INTRODUCTION In this note, the capabilities of ME scope will be used to build a model of the mass-spring-damper
More informationTime domain analysis of viscoelastic models
Time domain analysis of viscoelastic models L. Rouleau 1, J.-F. Deü 1 1 LMSSC, Conservatoire national des arts et métiers, 2D6R10, 292 rue Saint-Martin, 75141 Paris cedex 03, France e-mail: lucie.rouleau@lecnam.net
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 informationDynamic response of structures with uncertain properties
Dynamic response of structures with uncertain properties S. Adhikari 1 1 Chair of Aerospace Engineering, College of Engineering, Swansea University, Bay Campus, Fabian Way, Swansea, SA1 8EN, UK International
More informationRayleigh s Classical Damping Revisited
Rayleigh s Classical Damping Revisited S. Adhikari and A. Srikantha Phani Department of Aerospace Engineering, University of Bristol, Bristol, U.K. Email: S.Adhikari@bristol.ac.uk URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html
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 informationExact Solution of Time History Response for Dynamic Systems with Coupled Damping using Complex Mode Superposition
Clemson University TigerPrints All Theses Theses 8-2007 Exact Solution of Time History Response for Dynamic Systems with Coupled Damping using Complex Mode Superposition Manoj kumar Chinnakonda Clemson
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 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 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 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 informationCopyright 2009 Tech Science Press CMC, vol.298, no.1, pp.1-29, 2009
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 Copyright 2009 Tech Science Press CMC, vol.298, no.1, pp.1-29, 2009 Identification of dynamical systems with fractional derivative
More informationTime Response Analysis (Part II)
Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary
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 informationDamping modelling using generalized proportional damping
Journal of Sound and Vibration 9 (006) 156 10 JOURNAL OF SOUND AND VIBRATION www.elsevier.com/locate/jsvi Damping modelling using generalized proportional damping S. Adhikari Department of Aerospace Engineering,
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 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 informationSpectral methods for fuzzy structural dynamics: modal vs direct approach
Spectral methods for fuzzy structural dynamics: modal vs direct approach S Adhikari Zienkiewicz Centre for Computational Engineering, College of Engineering, Swansea University, Wales, UK IUTAM Symposium
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 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 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 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 informationEigensolution Derivatives for Arbitrarily Normalized Modes
B. P. Wang Department of Mechanical Aerospace Engineering The University of Texas at Arlington Box 19023, Arlington, TX 76019 Eigensolution Derivatives for Arbitrarily Normalized Modes Methods for computing
More informationIdentification of modal parameters from ambient vibration data using eigensystem realization algorithm with correlation technique
Journal of Mechanical Science and Technology 4 (1) (010) 377~38 www.springerlink.com/content/1738-494x DOI 107/s106-010-1005-0 Identification of modal parameters from ambient vibration data using eigensystem
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 2.6 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS
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 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 informationVISCO-ELASTIC ANALYSIS OF CANTILEVER BEAM USING AUGMENTED THERMODYNAMIC FIELD METHOD IN MATLAB
VISCO-ELASTIC ANALYSIS OF CANTILEVER BEAM USING AUGMENTED THERMODYNAMIC FIELD METHOD IN MATLAB 1 G.N.V.V. Manikanta, 2 P. Jaikishan, 3 E. Raghavendra Yadav 1 Research Scholar, 2,3 Asst. Professor, Aditya
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 informationOn Dominant Poles and Model Reduction of Second Order Time-Delay Systems
On Dominant Poles and Model Reduction of Second Order Time-Delay Systems Maryam Saadvandi Joint work with: Prof. Karl Meerbergen and Dr. Elias Jarlebring Department of Computer Science, KULeuven ModRed
More informationRandom Eigenvalue Problems in Structural Dynamics: An Experimental Investigation
Random Eigenvalue Problems in Structural Dynamics: An Experimental Investigation S. Adhikari, A. Srikantha Phani and D. A. Pape School of Engineering, Swansea University, Swansea, UK Email: S.Adhikari@swansea.ac.uk
More informationAn Improved Computational Strategy for Vibration- Proof Structures Equipped with Nano-Enhanced Viscoelastic Devices
An Improved Computational Strategy for Vibration- Proof Structures Equipped with Nano-Enhanced Viscoelastic Devices E. Ntotsios & A. Palmeri Loughborough University, United Kingdom SUMMARY: Viscoelastic
More informationNonlinear system identification with the use of describing functions a case study
Nonlinear system identification with the use of describing functions a case study Zhongge Zhao 1, Chuanri Li 2, Kjell Ahlin 3, Huan Du 4 1, 2, 4 School of Reliability and System Engineering, Beihang University,
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 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 informationChapter 8: Frequency Domain Analysis
Chapter 8: Frequency Domain Analysis Samantha Ramirez Preview Questions 1. What is the steady-state response of a linear system excited by a cyclic or oscillatory input? 2. How does one characterize the
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 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 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 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 informatione st f (t) dt = e st tf(t) dt = L {t f(t)} s
Additional operational properties How to find the Laplace transform of a function f (t) that is multiplied by a monomial t n, the transform of a special type of integral, and the transform of a periodic
More informationProbabilistic Structural Dynamics: Parametric vs. Nonparametric Approach
Probabilistic Structural Dynamics: Parametric vs. Nonparametric Approach S Adhikari School of Engineering, Swansea University, Swansea, UK Email: S.Adhikari@swansea.ac.uk URL: http://engweb.swan.ac.uk/
More informationFinite Element Modules for Demonstrating Critical Concepts in Engineering Vibration Course
Finite Element Modules for Demonstrating Critical Concepts in Engineering Vibration Course Shengyong Zhang Assistant Professor of Mechanical Engineering College of Engineering and Technology Purdue University
More informationTuning TMDs to Fix Floors in MDOF Shear Buildings
Tuning TMDs to Fix Floors in MDOF Shear Buildings This is a paper I wrote in my first year of graduate school at Duke University. It applied the TMD tuning methodology I developed in my undergraduate research
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 informationIn this lecture you will learn the following
Module 9 : Forced Vibration with Harmonic Excitation; Undamped Systems and resonance; Viscously Damped Systems; Frequency Response Characteristics and Phase Lag; Systems with Base Excitation; Transmissibility
More informationAn Efficient Computational Solution Scheme of the Random Eigenvalue Problems
50th AIAA SDM Conference, 4-7 May 2009 An Efficient Computational Solution Scheme of the Random Eigenvalue Problems Rajib Chowdhury & Sondipon Adhikari School of Engineering Swansea University Swansea,
More informationStructural Health Monitoring using Shaped Sensors
Structural Health Monitoring using Shaped Sensors Michael I. Friswell and Sondipon Adhikari Swansea University, UK This paper is concerned with distributed sensors to measure the response of beam and plate
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 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 informationDifferential Equations, Math 315 Midterm 2 Solutions
Name: Section: Differential Equations, Math 35 Midterm 2 Solutions. A mass of 5 kg stretches a spring 0. m (meters). The mass is acted on by an external force of 0 sin(t/2)n (newtons) and moves in a medium
More informationMath 308 Exam II Practice Problems
Math 38 Exam II Practice Problems This review should not be used as your sole source for preparation for the exam. You should also re-work all examples given in lecture and all suggested homework problems..
More informationSENSOR DESIGN FOR PIEZOELECTRIC CANTILEVER BEAM ENERGY HARVESTERS
SENSOR DESIGN FOR PIEZOELECTRIC CANTILEVER BEAM ENERGY HARVESTERS Michael I. Friswell and Sondipon Adhikari School of Engineering Swansea University Singleton Park, Swansea SA2 8PP, UK E-mail: m.i.friswell@swansea.ac.uk;
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 informationLecture 7: Laplace Transform and Its Applications Dr.-Ing. Sudchai Boonto
Dr-Ing Sudchai Boonto Department of Control System and Instrumentation Engineering King Mongkut s Unniversity of Technology Thonburi Thailand Outline Motivation The Laplace Transform The Laplace Transform
More information93 Analytical solution of differential equations
1 93 Analytical solution of differential equations 1. Nonlinear differential equation The only kind of nonlinear differential equations that we solve analytically is the so-called separable differential
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 information742. Time-varying systems identification using continuous wavelet analysis of free decay response signals
74. Time-varying systems identification using continuous wavelet analysis of free decay response signals X. Xu, Z. Y. Shi, S. L. Long State Key Laboratory of Mechanics and Control of Mechanical Structures
More informationSystem Parameter Identification for Uncertain Two Degree of Freedom Vibration System
System Parameter Identification for Uncertain Two Degree of Freedom Vibration System Hojong Lee and Yong Suk Kang Department of Mechanical Engineering, Virginia Tech 318 Randolph Hall, Blacksburg, VA,
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 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 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 informationThe quantum mechanics approach to uncertainty modeling in structural dynamics
p. 1/3 The quantum mechanics approach to uncertainty modeling in structural dynamics Andreas Kyprianou Department of Mechanical and Manufacturing Engineering, University of Cyprus Outline Introduction
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 informationTHE APPLICATION OF THE IRS AND BALANCED REALIZATION METHODS TO OBTAIN REDUCED MODELS OF STRUCTURES WITH LOCAL NON-LINEARITIES
Journal of Sound and Vibration (1996) 196(4), 453 468 THE APPLICATION OF THE IRS AND BALANCED REALIZATION METHODS TO OBTAIN REDUCED MODELS OF STRUCTURES WITH LOCAL NON-LINEARITIES M. I. FRISWELL Department
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 informationStructural Dynamics Modification and Modal Modeling
SEM Handboo of Experimental Structural Dynamics - Structural Dynamics Modification and Modal Modeling Structural Dynamics Modification and Modal Modeling Structural Dynamics Modification (SDM) also nown
More informationOn the Nature of Random System Matrices in Structural Dynamics
On the Nature of Random System Matrices in Structural Dynamics S. ADHIKARI AND R. S. LANGLEY Cambridge University Engineering Department Cambridge, U.K. Nature of Random System Matrices p.1/20 Outline
More informationDamped Harmonic Oscillator
Damped Harmonic Oscillator Wednesday, 23 October 213 A simple harmonic oscillator subject to linear damping may oscillate with exponential decay, or it may decay biexponentially without oscillating, or
More informationProblem 4.4 Derive the ramp response by nding the particular and homogeneous solutions of the ordinary di. equation.
Homework #4 Name: HyoSoo Kim UFID:8989-3565 Date:/6/6 Problem 4.4 Derive the ramp response by nding the particular and homogeneous solutions of the ordinary di. equation. The unit ramp fuction is r(t)
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 informationSensitivity analysis of viscoelastic structures
Shock and Vibration 13 (26) 545 558 545 IOS Press Sensitivity analysis of viscoelastic structures A.M.G. de Lima a, M.H. Stoppa b, D.A. Rade a, and V. Steffen. Jr. a a Federal University of Uberlândia,
More informationStructural Matrices in MDOF Systems
in MDOF Systems http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April 9, 2016 Outline Additional Static Condensation
More informationOrder Reduction of Parametrically Excited Linear and Nonlinear Structural Systems
Venkatesh Deshmukh 1 Eric A. Butcher e-mail: ffeab@uaf.ediu Department of Mechanical Engineering, University of Alaska Fairbanks, Fairbanks, AK 99775 S. C. Sinha Nonlinear Systems Research Laboratory,
More informationIn-Structure Response Spectra Development Using Complex Frequency Analysis Method
Transactions, SMiRT-22 In-Structure Response Spectra Development Using Complex Frequency Analysis Method Hadi Razavi 1,2, Ram Srinivasan 1 1 AREVA, Inc., Civil and Layout Department, Mountain View, CA
More informationDifferential equations
Differential equations Math 7 Spring Practice problems for April Exam Problem Use the method of elimination to find the x-component of the general solution of x y = 6x 9x + y = x 6y 9y Soln: The system
More informationABSTRACT Modal parameters obtained from modal testing (such as modal vectors, natural frequencies, and damping ratios) have been used extensively in s
ABSTRACT Modal parameters obtained from modal testing (such as modal vectors, natural frequencies, and damping ratios) have been used extensively in system identification, finite element model updating,
More informationTheory of Vibrations in Stewart Platforms
Theory of Vibrations in Stewart Platforms J.M. Selig and X. Ding School of Computing, Info. Sys. & Maths. South Bank University London SE1 0AA, U.K. (seligjm@sbu.ac.uk) Abstract This article develops a
More informationEE102 Homework 2, 3, and 4 Solutions
EE12 Prof. S. Boyd EE12 Homework 2, 3, and 4 Solutions 7. Some convolution systems. Consider a convolution system, y(t) = + u(t τ)h(τ) dτ, where h is a function called the kernel or impulse response of
More informationA Vector-Space Approach for Stochastic Finite Element Analysis
A Vector-Space Approach for Stochastic Finite Element Analysis S Adhikari 1 1 Swansea University, UK CST2010: Valencia, Spain Adhikari (Swansea) Vector-Space Approach for SFEM 14-17 September, 2010 1 /
More informationModal Decomposition and the Time-Domain Response of Linear Systems 1
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING.151 Advanced System Dynamics and Control Modal Decomposition and the Time-Domain Response of Linear Systems 1 In a previous handout
More informationREVIEW NOTES FOR MATH 266
REVIEW NOTES FOR MATH 266 MELVIN LEOK 1.1: Some Basic Mathematical Models; Direction Fields 1. You should be able to match direction fields to differential equations. (see, for example, Problems 15-20).
More informationA Multi-Step Hybrid Method for Multi-Input Partial Quadratic Eigenvalue Assignment with Time Delay
A Multi-Step Hybrid Method for Multi-Input Partial Quadratic Eigenvalue Assignment with Time Delay Zheng-Jian Bai Mei-Xiang Chen Jin-Ku Yang April 14, 2012 Abstract A hybrid method was given by Ram, Mottershead,
More information+ + LAPLACE TRANSFORM. Differentiation & Integration of Transforms; Convolution; Partial Fraction Formulas; Systems of DEs; Periodic Functions.
COLOR LAYER red LAPLACE TRANSFORM Differentiation & Integration of Transforms; Convolution; Partial Fraction Formulas; Systems of DEs; Periodic Functions. + Differentiation of Transforms. F (s) e st f(t)
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 information