Transactions on Modelling and Simulation vol 20, 1998 WIT Press, ISSN X
|
|
- Beryl Harris
- 6 years ago
- Views:
Transcription
1 A BEM-formulation applied to the solution of nonlinear dynamic problems H.J. Holl, A.K. Belyaev & H. Irschik Division of Technical Mechanics, University Linz, Austria uni-linz. ac. at Abstract A semi-analytical time-integration procedure is presented for the integration of discretized dynamic mechanical systems. This method utilizes the advantages of the boundary element method, well known from quasi-static field problems Motivated by these spatial formulations, the present dynamic method is based on influence functions in time, and gives exact solutions in the linear time-invariant case. Similar to domain-type BEM's for nonlinear field problems, the method is extended for different nonlinear dynamic systems having nonclassical damping The numerical stability and accuracy of the semi-analytical method is discussed in two steps for the nonclassical damping and for the nonlinear restoring forces The Duffing oscillator and a system with nonclassical damping are used as representative model problems. A comparison is given to other conventionally used time integration procedures, which shows the efficiency of the present method. 1 Introduction During the last decades, various powerful computational methods have been developed for nonlinear structural dynamics. A main stage of the computational analysis is the space-wise discretization of the governing field equations leading to a system of ordinary differential equations in time. Especially, a semidiscretization by the Finite Element Method (FEM) leads to a dynamic system of second-order equations of the form '* (1) compare Zienkiewicz & Taylor [1], where X denotes the vector of generalised coordinates, M stands for the mass matrix, K and C are matrices of stiffness and
2 300 Boundary Elements damping, respectively, and F gathers the imposed forcing terms. In many practical situations nonlinearities do occur, the matrices in eqn (1) then depending on the generalised coordinates: M = m + M(X,X), C = c + C(X,X), KX = kx + Q(X,X), (2) where m, c and k denote the symmetric and constant portions of the matrices in eqn (1). The damping matrix c is of the Rayleigh type. It has been shown in Zienkiewicz, Taylor [1] that the numerical integration of eqn (1) by a finite element discretization in time provides a strategy unifying many existing singleand multi-step algorithms and providing a variety of new ones. When discussing stability problems of these schemes, it is useful to revert to modally uncoupled equations, see [1] and Hughes [2]: Xk-f2CkO^Xk + o^xk = pk = fk + rk(x,x,x), (3) where the index (k) refers to the number of the mode, which is supressed in the following. The modal analysis transforming eqn (1) into eqn (3) is performed with respect to the constant symmetric matrices m, c and k of eqn (2), while the remaining unsymmetric and nonlinear portions are gathered in the k-th modal force r. The modal force due to F is denoted by f, co is for the natural frequency, and C denotes slight modal damping. With respect to eqn (3), it is noted that our group has been successful to solve various structural elasto-viscoplastic problems by utilizing an analogy between the inelastic strains and a distribution of eigenstrains in a linear-elastic background-structure and by taking the continuous modes of that linear structure for a modal reduction to equations of the type of eqn (3). The nonlinear restoring forces q have been calculated from the constitutive equations in a time-stepping procedure, and a semi-analytic timeintegration of eqn (3) has been performed, see Irschik & Ziegler [3], [4] and Fotiu, Irschik & Ziegler [5] for a recent application relating modal analysis to integral equations for elastic -plastic plate vibrations. For the discretized case of nonlinear rotordynamics, an analogous semi-analytic time-integration scheme has been developed by Holl [6], [7] and extended to systems with time-dependent mass in [8]. Seen from the point of the FEM, however, the time-integration of eqns (1) or (3) by means offiniteelements in time of course is a concise strategy. One should remind, nevertheless, that a formulation of the type of eqn (1) can also be obtained by utilising the Boundary Element Method (BEM) for the space-wise semi-discretization, Nardini & Brebbia [9], see Manolis & Beskos [10] for a review. In that connection, it is the scope of the present contribution to discuss the integration of eqn (3) by means of a BEM in time. It is shown that a timewise BEM corresponds to an exact integration of the modal equations in the linear case, while a domain-boundary element formulation (DBEM) is obtained for the case of nonlinear modal restoring forces, compare [11] for DBEMs of quasi-static inelastic problems. The time-domain integrals of the present
3 Boundary Elements 301 formulation coincide with Duhamel's convolution integral, well-known from the theory of structural vibrations, see e.g. Ziegler [12]. In contrast to an approximation with FEM in time, the present BEM does not approximate the generalised coordinates, but only the integrands of the domain integrals occurring in the nonlinear case are approximated. In the present paper the numerical stability is discussed for the nonlinear case, where a representative approximation of the domain-integrals is used. Results are shown for the Duffing oscillator and are compared to the outcomes of the Newmark method. As a further example a system with time-varying mass and nonclassical damping is analysed. 2 Boundary Integral Formulation in Time Following the fundamental concepts of BEM presented by Brebbia & Dominguez [13], we start with a weighted residual formulation for eqn (3) within the time interval: At [X(t) + 2 C co X(t) + 0)2 X(t) - p(t)] W(t) dt = 0; 0<t<At. (4) W(t) denotes a sufficiently continuous weight function. Integrating by parts twice gives - 2 c co W(t) + of w(t)) x(t) - p(t) w(t)] dt + x(t) w(t) * + 2(toX(t)W(t) '-X(t)\V(t) At = 0- (5) In order to relate eqn (5) to BEM, we follow the singular approach indicated by Brebbia & Dominguez [13]. In the sense of this latter approach, the weight function W is identified by the fundamental solution W, of the modal oscillator due to an unit impulse acting at the inner site of the boundary t=0. W, is governed by the differential equation Vv-i-2Cco\V,+co*Wi=S(O, where (6) Dirac's delta function is denoted by 8, and t+ = lim^ <# - e). W, satisfies homogeneous initial conditions at t=0. As a second fundamental solution, the solution Wj of the modal oscillator due to a unit dipole 6 is introduced:
4 302 Boundary Elements, where [cos (# + ^ sin co^t ] H(Q. (8) % Applying both eqns (7) and (8) to eqn (5), considering the properties of 8 and S, r X(t) 5(t) dt = X(t=0) = Xo, (^ X(t) S(t) dt = - X(t=0) = -*o, (9) Jo Jo using some trigonometric transformation rules and solving for the state-vector at the end of the time-interval, the following transition matrix formulation is obtained: It is thus seen that the evaluation of the weighted residual statement of eqn (5) in the sense of BEM leads to the exact integration of the modal equations (3) by means of Duhamel's integral, since A(At) is the exact transition matrix of the linear oscillator, compare Hughes [2], and see Refs. [3]- [6], where a semianalytical time-stepping procedure analogous to eqn (10) has been used. It is certainly not surprising that the use of BEM in time leads to exact formulations, since the boundaries of the time-interval reduce to a two-point problem, with algebraic relations instead of boundary integrals. In the case of nonlinear problems, the portion r of pin eqn (3) is not known in advance. An iterative procedure has to be established in the sense of the DBEM by approximating the time-evolution of q in eqn (10), see the following Section. The solution for the load transfer matrix of the k-th modal equations of motion due to the excitation forces can be computed using eqn (10), when inserting P(T) = tmco- The excitation force in the integral is approximated in the form of a parabolic interpolation function of the force within the m-th time step according i ( i \* to fkmw = 4mO + ^km^ + 4AM^I 1 We can define a corresponding forcing vector with the three parameters of the interpolation function in the form ^km = { 4m A4m ^4m) ^ the following the index k, which refers to the number of the generalized coordinate, is omitted as it is obvious in these equations. The evaluation of the convolution integral of eqn (10) for the above stated interpolation function (10) gives the corresponding solution of the elements of the load transfer matrix B(At) involving dynamic Green's functions. If we write the result for the generalized coordinate at the end of the m-th time step, we get for a linear system
5 Boundary Elements 393 A'»-PBfp. p = l This is the exact solution if the parabolic interpolation of the excitation force is exact. The suitable application of modal reduction gives additional advantages for the computational effort and can be included with some small changes in the present algorithm. 3 Analysis of Systems with Nonlinear Restoring Forces * The transition-matrix formulation of eqn (10) provides exact solutions for linear vibration^ In order to obtain a measure for the quality of various approximations of the modal force r in the nonlinear case, however, free damped linear vibrations are studied at first, where the damping force formally is treated as a nonlinear restoring force: X + co*x = r(x) = -2CX Using the undamped matrix A(%., the solution of eqn (12) is studied for the approximation of the integrand in eqn (10): VY-\ _ Y j. V. XCD-XO + XOT 2 At Inserting eqn (13) into (10), integrating and rearranging terms gives the form /X(At)\ \*w>r*"\*.,. 04) where A is the approximate transfer matrix. This equation is appropriate for a stability discussion. The stability analysis is performed for a linear initial value problem which is extended to the consideration of the above described effects. The analysis steps for the linear case are taken from [2]. So it can be stated, that the absolute value of the maximum eigenvalues has to be less than one for every numerical algorithm if it is said to be unconditionally stable. The eigenvalues of the exact transfer matrix A are ^=e( C»±i«*)t ^absolute value of the eigenvalues is computedtop max Ap j -/Aj - e * «", which is also called spectral radius. The presented semi-analytic procedure has an approximate transfer mairix A and
6 304 Boundary Elements can be analysed with respect to its numerical behaviour. In [6] the equations for the numerical dissipation, the numerical dispersion and the accuracy are given for the approximate transfer matrix. When considering a linear modal oscillator it is confirmed, that the present BEM-type method is exact. The numerical dissipation is the sum of the modal damping of the oscillator and an additional damping due to the assumptions of the algorithm. The numerical dispersion is characterized by the resulting relative period error. The accuracy is computed from the local error for one time step. The advantage of the present method for the linear system results due to the application of the exact transfer matrix, so the numerical dissipation and dispersion and the local error are zero for the linear symmetric system. 1,2 "S <o 0,8 --* 0,6 '- 0,4-0,2 -, 0 0,5..., 1., _ 1,5 2 dunensionsless lime step At/I Fig. 1: Absolute values of the eigenvalue of the approximate transition matrix The absolute values X for the case of eqn (13) are shown in Fig. 1 as a function of At / T, where T denotes the period of the undamped vibration. Also shown are the exact solutions in the form of dashed lines. The errors accumulated within the time-interval while compared to the exact solution become: (15) It is seen that the accuracy is of second order. (Note that these results are purely artificial, since in practice one would treat linear oscillations by means of the exact transition matrix.) Similar estimates can be obtained for the errors in the vibration period and for the phenomenon of numerical damping (see Holl [6] for details): Here for the above algorithm some additional results are presented for the case of the equation of motion
7 Boundary Elements 305 * (17) The exact transfer matrix formulation for this case is knwon. The eqn (10) is calculated using the convolution integral and inserting the approximation _ Y - A + - T. The approximate transfer matrix A results after some algebraic manipulations. Fig. 2 shows the absolute values of the eigenvalues of the resulting algorithm with the parameters a = 0.3 and = 0.1. An analogous analysis was performed for the Newmark method. Based on the characteristic equations for 7= 1 and 0 = 1, see [2], [6], an application to eqn (17) was computed. In Fig. 2 the corresponding result for the Newmark method can be seen. 1 S? 3 a Newrr ark me od J\ / E act valu *s^ Pres ^s^ ent meth od ^ 0.2 k C=o.i 0=0.3 0 _L J L _L 0 0, Dimensionsless time step At/T 1 Fig. 2 Absolute values of the eigenvalues for the system of eqn (17) The result for the accuracy of the above presented algorithm is where XQ and XQ are the initial values at the beginning of the time step For the Newmark method the local error for = 0 results to (18) (19)
8 306 Boundary Elements We see, that y= ^ gives an algorithm, which has a convergence of second order. This result for the convergence does not change for the linear case a = 0, whereas from eqn (18) it is obvious, that for a = 0 the local error disappears, i.e. K = 0. 4 Examples 4.1 Free Vibrations of an Undamped Duffing Oscillator The above BEM-type time-integration scheme is applied to free vibrations of a Duffing oscillator, r = ~px^ with initial conditions 3 X = 4 of and %o = 0. The approximation (13) is used in a time-stepping procedure with CD At = 2 7t / 100. The increments of the state vector were calculated solving a cubic equation. Fig. 3 shows the error in the displacement e(t) = c exact versus the error in velocity a(t) = for the present method and for the Newmark method. It is seen that the present method is more accurate, while the CPU-time of the present method was smaller by a factor of 1/2. Comparisons with other direct time-integration schemes turned out to give an even better benefit for our method. 100 Present Method Newmark Method e(t) Fig. 3 Phase-plane of the non-dimensional error for the Duffing oscillator
9 Boundary Elements Examples of Systems with Time-dependent Mass As a further example a system with two DOF is analysed, where the matrices are taken from [14]: M = 110] -«<o,t c_[l j [3.8 Ol loir ' M ' *= 0 4 ' = lsin(vt) withv/g>j=l Itfe assumed, that this system has a time varying mass according to the function e with a = 1 /10, which reduces all elements of the mass matrix proporrt^j ^ S»ven system the computed results are shown in Fig. 4 for the first DOF and for homogeneous initial conditions in a dimensionless form. X,, is the static displacements. The four computed cases are: A) with full damping matrix; B) decoupling approximation when neglecting off-diagonal elements- C) as A but with time-variable mass; D) as B but withtime-variablemass I -1,5 0,5 1 1,5 2 Dimensionless time t/t Fig.4 Converged solution for the first DOF 2,5 5 Conclusion The presented semi-analytic algorithm makes use of the methods of linear structural dynamics and applies them in a consistent way to the analysis of transient vibrations of nonlinear rotordynamic systems. It is demonstrated that the resulting algorithm is efficient and accurate. The computation time for solving example problems shows an advantage for the present method when
10 308 Boundary Elements compared to the Newmark method. When using the same time step, the error to a converged solution in the computed results showed an advantage as well. References [1] Zienkiewicz, O.C.; Taylor, R.L.: The Finite Element Method, 4th edition, Vol. 2, McGraw Hill [2] Hughes, T.J.R.: Analysis of Transient Algorithms with Particular Reference to Stability Behavior. In: Computational Methods for Transient Analysis, eds.: Belytschko, T.; Hughes, T.J.R.; North-Holland, Elsevier Science Publishers B.V., [3] Irschik, H.; Ziegler, R: Dynamics of Linear Elastic Structures with Selfstress: A Unified Treatment for Linear and Nonlinear Problems. ZAMM 68 (6), [4] Irschik, H.; Ziegler, F.: Dynamic Processes in Structural Thermo- Viscoplasticity. Applied Mechanics Review 48 (6), , [5] Fotiu, P.; Irschik, H.; Ziegler, F.: Modal Analysis of Elastic-plastic Plate Vibrations by Integral Equations. Engineering Analysis with Boundary Elements 14, Elsevier, [6] Holl, H.J.: Ein effizienter Algorithmus fur nichtlineare Probleme der Strukturdynamik mit Anwendung in der Rotordynamik. Dissertation, University Linz, [7] Holl, H.J.: An Time Integration Algorithm For Time-Varying Systems With Nonclassical Damping Based on Modal Methods. Proceedings of the 15th IMAC, Orlando, ,1997. [8] Holl, H.J.: An Effizient Time Domain Formulation For Nonlinear Rotordynamic Systems Using Modal Reduction. Proceedings of the 13th IMAC, Nashville, [9] Nardini, D.; Brebbia, C.A.: Transient Dynamic Analysis by the Boundary Element Method. In: Boundary Elements, eds.: Brebbia, C.A.; Futagami, T.; Tanaka, M.; Springer Verlag , [10] Manolis, G.D.; Beskos, D.E.: Boundary Element Methods in Elastodynamics. Unwin Hyman [11] Mukherjee, S.: Boundary Element Methods in Creep and Fracture. Applied Science Publisher [12] Ziegler, F.: Mechanics of Solids and Fluids. Springer Verlag [13] Brebbia, C.A.; Dominguez, J.: Boundary Elements. 2nd edition, McGraw Hill [14] Hwang, J.H.; Ma, F.: On the Approximate Solution of Nonclassically Damped Linear Systems. Journal of Applied Mechanics 60, , 1993.
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 informationThe Finite Element Method for Solid and Structural Mechanics
The Finite Element Method for Solid and Structural Mechanics Sixth edition O.C. Zienkiewicz, CBE, FRS UNESCO Professor of Numerical Methods in Engineering International Centre for Numerical Methods in
More informationTransactions on Modelling and Simulation vol 8, 1994 WIT Press, ISSN X
Model analysis of plates using the dual reciprocity boundary element method T.W. Davies & F.A. Moslehy Department of Mechanical and Aerospace Engineering, University of Central Florida, Orlando, Florida,
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 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 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 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 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 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 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 informationCOPYRIGHTED MATERIAL. Index
Index A Admissible function, 163 Amplification factor, 36 Amplitude, 1, 22 Amplitude-modulated carrier, 630 Amplitude ratio, 36 Antinodes, 612 Approximate analytical methods, 647 Assumed modes method,
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
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 information4.4 1) 단순지지된깊은보 선형동적해석검증예제 ANALYSIS REFERENCE. REFERENCE NAFEMS 1 Beam elements, solid elements
그림 5.4.3 가진방향에따른응답변화예시 Reaction Sum. Total Strain Energy 0 30 60 90 120 150 180 Excitation ngle 4.4 선형동적해석검증예제 1) 단순지지된깊은보 REFERENCE NFEMS 1 ELEMENTS Beam elements, solid elements MODEL FILENME LinearDynamic01.mpb
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 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 informationEXTENDED ABSTRACT. Dynamic analysis of elastic solids by the finite element method. Vítor Hugo Amaral Carreiro
EXTENDED ABSTRACT Dynamic analysis of elastic solids by the finite element method Vítor Hugo Amaral Carreiro Supervisor: Professor Fernando Manuel Fernandes Simões June 2009 Summary The finite element
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 informationDYNAMIC RESPONSE OF EARTHQUAKE EXCITED INELASTIC PRIMARY- SECONDARY SYSTEMS
DYNAMIC RESPONSE OF EARTHQUAKE EXCITED INELASTIC PRIMARY- SECONDARY SYSTEMS Christoph ADAM 1 And Peter A FOTIU 2 SUMMARY The objective of the paper is to investigate numerically the effect of ductile material
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 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 informationTransactions on Modelling and Simulation vol 18, 1997 WIT Press, ISSN X
A simple coupling of 2D BEM and FEM bar model applied to mass matrix elastodynamic analysis Coda,H.B*; Venturini W.S* and Aliabadi, M.HT * Sao Carlos School of Eng., University of Sao Paulo, Br, venturin@usp.
More informationVibration Dynamics and Control
Giancarlo Genta Vibration Dynamics and Control Spri ringer Contents Series Preface Preface Symbols vii ix xxi Introduction 1 I Dynamics of Linear, Time Invariant, Systems 23 1 Conservative Discrete Vibrating
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 informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Simplifying integration for logarithmic singularities R.N.L. Smith Department ofapplied Mathematics & OR, Cranfield University, RMCS, Shrivenham, Swindon, Wiltshire SN6 SLA, UK Introduction Any implementation
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 informationA new simple method of implicit time integration for dynamic problems of engineering structures
Acta Mech Sin 2007 23:91 99 DOI 10.1007/s1009-006-0050-5 RESEARCH PAPER A new simple method of implicit time integration for dynamic problems of engineering structures Jun Zhou Youhe Zhou Received: 17
More informationmidas Civil Dynamic Analysis
Edgar De Los Santos Midas IT August 23 rd 2017 Contents: Introduction Eigen Value Analysis Response Spectrum Analysis Pushover Analysis Time History Analysis Seismic Analysis Seismic Analysis The seismic
More informationDISPENSA FEM in MSC. Nastran
DISPENSA FEM in MSC. Nastran preprocessing: mesh generation material definitions definition of loads and boundary conditions solving: solving the (linear) set of equations components postprocessing: visualisation
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 informationENGI 9420 Lecture Notes 1 - ODEs Page 1.01
ENGI 940 Lecture Notes - ODEs Page.0. Ordinary Differential Equations An equation involving a function of one independent variable and the derivative(s) of that function is an ordinary differential equation
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 informationShape Optimization of Revolute Single Link Flexible Robotic Manipulator for Vibration Suppression
15 th National Conference on Machines and Mechanisms NaCoMM011-157 Shape Optimization of Revolute Single Link Flexible Robotic Manipulator for Vibration Suppression Sachindra Mahto Abstract In this work,
More informationA STRATEGY FOR IDENTIFICATION OF BUILDING STRUCTURES UNDER BASE EXCITATIONS
A STRATEGY FOR IDENTIFICATION OF BUILDING STRUCTURES UNDER BASE EXCITATIONS G. Amato and L. Cavaleri PhD Student, Dipartimento di Ingegneria Strutturale e Geotecnica,University of Palermo, Italy. Professor,
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 informationPreprocessor Geometry Properties )Nodes, Elements(, Material Properties Boundary Conditions(displacements, Forces )
در برنامه يك تدوين براي بعدي دو يك سازه محيط MATLAB Preprocessor Geometry Properties )Nodes, Elements(, Material Properties Boundary Conditions(displacements, Forces ) Definition of Stiffness Matrices
More informationA MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC PROBLEMS IN SEMI-INFINITE MEDIA
8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 A MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC
More informationApplication of pseudo-symmetric technique in dynamic analysis of concrete gravity dams
Application of pseudo-symmetric technique in dynamic analysis of concrete gravity dams V. Lotfi Department of Civil and Environmental Engineering, Amirkabir University, Iran Abstract A new approach is
More informationON THE ACCURACY AND COMPUTATIONAL COST OF TIME HISTORY ANALYSIS OF RESIDENTIAL BUILDINGS BY THE QUASI- WILSON-THETA METHOD
11 th International Conference on Vibration Problems Z. Dimitrovová et al. (eds.) Lisbon, Portugal, 9-12 September 2013 ON THE ACCURACY AND COMPUTATIONAL COST OF TIME HISTORY ANALYSIS OF RESIDENTIAL BUILDINGS
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 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 informationON THE INTEGRATION OF EQUATIONS OF MOTION: FEM AND MOLECULAR DYNAMICS PROBLEMS
8th International Congress on Computational Mechanics, Volos, 1-15 July 015 ON THE INTEGRATION OF EQUATIONS OF MOTION: FEM AND MOLECULAR DYNAMICS PROBLEMS E.G. Kakouris, V.K. Koumousis Institute of Structural
More informationSPECIAL DYNAMIC SOIL- STRUCTURE ANALYSIS PROCEDURES DEMONSTATED FOR TWO TOWER-LIKE STRUCTURES
2010/2 PAGES 1 8 RECEIVED 21. 9. 2009 ACCEPTED 20. 1. 2010 Y. KOLEKOVÁ, M. PETRONIJEVIĆ, G. SCHMID SPECIAL DYNAMIC SOIL- STRUCTURE ANALYSIS PROCEDURES DEMONSTATED FOR TWO TOWER-LIKE STRUCTURES ABSTRACT
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 informationDynamic analysis of railway bridges by means of the spectral method
Dynamic analysis of railway bridges by means of the spectral method Giuseppe Catania, Silvio Sorrentino DIEM, Department of Mechanical Engineering, University of Bologna, Viale del Risorgimento, 436 Bologna,
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 informationelastoplastic contact problems D. Martin and M.H. Aliabadi Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK
Non-conforming BEM elastoplastic contact problems D. Martin and M.H. Aliabadi discretisation in Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK Abstract In this paper,
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 informationHydro-elastic Wagner impact using variational inequalities
Hydro-elastic Wagner impact using variational inequalities Thomas GAZZOLA, Alexander KOROBKIN, Šime MALENICA Introduction A simple model of water impact has been introduced by Wagner [6]. This model is
More 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 informationDYNAMIC RESPONSE OF THIN-WALLED GIRDERS SUBJECTED TO COMBINED LOAD
DYNAMIC RESPONSE OF THIN-WALLED GIRDERS SUBJECTED TO COMBINED LOAD P. WŁUKA, M. URBANIAK, T. KUBIAK Department of Strength of Materials, Lodz University of Technology, Stefanowskiego 1/15, 90-924 Łódź,
More informationFinite Element Analysis Lecture 1. Dr./ Ahmed Nagib
Finite Element Analysis Lecture 1 Dr./ Ahmed Nagib April 30, 2016 Research and Development Mathematical Model Mathematical Model Mathematical Model Finite Element Analysis The linear equation of motion
More informationAdvanced Vibrations. Distributed-Parameter Systems: Approximate Methods Lecture 20. By: H. Ahmadian
Advanced Vibrations Distributed-Parameter Systems: Approximate Methods Lecture 20 By: H. Ahmadian ahmadian@iust.ac.ir Distributed-Parameter Systems: Approximate Methods Rayleigh's Principle The Rayleigh-Ritz
More informationCOMPUTATION OF MASS OSCILLATIONS IN A SURGE TANK BY FINITE ELEMENT TECHNIQUE
Larhyss Journal, ISSN -368, n 5, Septembre 3, pp. 39-49 3 Tous droits réservés COMPUTATION OF MASS OSCILLATIONS IN A SURGE TANK BY FINITE ELEMENT TECHNIQUE AMARA L., BERREKSI A., ABDOUNE K. 3 Research
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 informationEffects of damping matrix in the response of structures with added linear viscous dampers
Effects of damping matrix in the response of structures with added linear viscous dampers J.R. Arroyo', J. Marte2 l Department of General Engineering, University of Puerto Rico, Mayagiiez Campus 'Department
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 informationTransactions on Engineering Sciences vol 6, 1994 WIT Press, ISSN
A computational method for the analysis of viscoelastic structures containing defects G. Ghazlan," C. Petit," S. Caperaa* " Civil Engineering Laboratory, University of Limoges, 19300 Egletons, France &
More informationA Fast Newton-Raphson Method in Stochastic Linearization
A Fast Newton-Raphson Method in Stochastic Linearization Thomas Canor 1, Nicolas Blaise 1, Vincent Denoël 1 1 Structural Engineering Division, University of Liège, Chemin des Chevreuils, B52/3, 4000 Liège,
More informationAPPLICATION OF VARIATIONAL METHODS AND GALERKIN METHOD IN SOLVING ENGINEERING PROBLEMS REPRESENTED BY ORDINARY DIFFERENTIAL EQUATIONS
APPLICATION OF VARIATIONAL METHODS AND GALERKIN METHOD IN SOLVING ENGINEERING PROBLEMS REPRESENTED BY ORDINARY DIFFERENTIAL EQUATIONS 1 B.V. SIVA PRASAD REDDY, 2 K. RAJESH BABU 1,2 Department of Mechanical
More informationNotes for Expansions/Series and Differential Equations
Notes for Expansions/Series and Differential Equations In the last discussion, we considered perturbation methods for constructing solutions/roots of algebraic equations. Three types of problems were illustrated
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 informationThursday, August 4, 2011
Chapter 16 Thursday, August 4, 2011 16.1 Springs in Motion: Hooke s Law and the Second-Order ODE We have seen alrealdy that differential equations are powerful tools for understanding mechanics and electro-magnetism.
More informationChapter 23: Principles of Passive Vibration Control: Design of absorber
Chapter 23: Principles of Passive Vibration Control: Design of absorber INTRODUCTION The term 'vibration absorber' is used for passive devices attached to the vibrating structure. Such devices are made
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 informationa) Find the equation of motion of the system and write it in matrix form.
.003 Engineering Dynamics Problem Set Problem : Torsional Oscillator Two disks of radius r and r and mass m and m are mounted in series with steel shafts. The shaft between the base and m has length L
More informationResearch Article A Note on the Solutions of the Van der Pol and Duffing Equations Using a Linearisation Method
Mathematical Problems in Engineering Volume 1, Article ID 693453, 1 pages doi:11155/1/693453 Research Article A Note on the Solutions of the Van der Pol and Duffing Equations Using a Linearisation Method
More informationSeismic analysis of structural systems with uncertain damping
Seismic analysis of structural systems with uncertain damping N. Impollonia, G. Muscolino, G. Ricciardi Dipartirnento di Costruzioni e Tecnologie Avanzate, Universita di Messina Contrada Sperone, 31 -
More informationDue Monday, November 16 th, 12:00 midnight
Due Monday, November 16 th, 12:00 midnight This homework is considering the finite element analysis of transient and dynamic FEM analysis. You are asked to include transient and/dynamic effects to MatLab
More information22.3. Repeated Eigenvalues and Symmetric Matrices. Introduction. Prerequisites. Learning Outcomes
Repeated Eigenvalues and Symmetric Matrices. Introduction In this Section we further develop the theory of eigenvalues and eigenvectors in two distinct directions. Firstly we look at matrices where one
More informationStatic & Dynamic. Analysis of Structures. Edward L.Wilson. University of California, Berkeley. Fourth Edition. Professor Emeritus of Civil Engineering
Static & Dynamic Analysis of Structures A Physical Approach With Emphasis on Earthquake Engineering Edward LWilson Professor Emeritus of Civil Engineering University of California, Berkeley Fourth Edition
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 information1 Introduction IPICSE-2016
(06) DOI: 0.05/ matecconf/06860006 IPICSE-06 Numerical algorithm for solving of nonlinear problems of structural mechanics based on the continuation method in combination with the dynamic relaxation method
More informationFOURIER TRANSFORMS. At, is sometimes taken as 0.5 or it may not have any specific value. Shifting at
Chapter 2 FOURIER TRANSFORMS 2.1 Introduction The Fourier series expresses any periodic function into a sum of sinusoids. The Fourier transform is the extension of this idea to non-periodic functions by
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 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 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 informationSTATIC AND DYNAMIC ANALYSIS OF A BISTABLE PLATE FOR APPLICATION IN MORPHING STRUCTURES
STATIC AND DYNAMIC ANALYSIS OF A BISTABLE PLATE FOR APPLICATION IN MORPHING STRUCTURES A. Carrella 1, F. Mattioni 1, A.A. Diaz 1, M.I. Friswell 1, D.J. Wagg 1 and P.M. Weaver 1 1 Department of Aerospace
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 informationRandom Vibrations & Failure Analysis Sayan Gupta Indian Institute of Technology Madras
Random Vibrations & Failure Analysis Sayan Gupta Indian Institute of Technology Madras Lecture 1: Introduction Course Objectives: The focus of this course is on gaining understanding on how to make an
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 informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Plate-soil elastodynamic coupling using analysis S.F.A. Baretto, H.B. Coda, W.S. Venturini Sao Carlos School of Engineering, University ofsao Paulo, Sao Carlos - SP, Brazil BEM Abstract The aim of this
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 informationStructural Dynamics Lecture 7. Outline of Lecture 7. Multi-Degree-of-Freedom Systems (cont.) System Reduction. Vibration due to Movable Supports.
Outline of Multi-Degree-of-Freedom Systems (cont.) System Reduction. Truncated Modal Expansion with Quasi-Static Correction. Guyan Reduction. Vibration due to Movable Supports. Earthquake Excitations.
More informationComputation Time Assessment of a Galerkin Finite Volume Method (GFVM) for Solving Time Solid Mechanics Problems under Dynamic Loads
Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 215 Paper o. 31 Computation Time Assessment of a Galerkin Finite Volume
More informationFinal Exam December 11, 2017
Final Exam Instructions: You have 120 minutes to complete this exam. This is a closed-book, closed-notes exam. You are NOT allowed to use a calculator with communication capabilities during the exam. Usage
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 informationTHE REFLECTION PHENOMENA OF SV-WAVES IN A GENERALIZED THERMOELASTIC MEDIUM
Internat. J. Math. & Math. Sci. Vol., No. 8 () 59 56 S67 Hindawi Publishing Corp. THE REFLECTION PHENOMENA OF SV-WAVES IN A GENERALIZED THERMOELASTIC MEDIUM ABO-EL-NOUR N. ABD-ALLA and AMIRA A. S. AL-DAWY
More informationAdditive resonances of a controlled van der Pol-Duffing oscillator
Additive resonances of a controlled van der Pol-Duffing oscillator This paper has been published in Journal of Sound and Vibration vol. 5 issue - 8 pp.-. J.C. Ji N. Zhang Faculty of Engineering University
More informationResponse Analysis of thin-walled Structure under Non-uniform temperature field and Acoustic Loads Yundong Sha, Xiyang Zheng
International Conference on Advances in Mechanical Engineering and Industrial Informatics (AMEII 2015) Response Analysis of thin-walled Structure under Non-uniform temperature field and Acoustic Loads
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 informationDYNAMIC FAILURE ANALYSIS OF LAMINATED COMPOSITE PLATES
Association of Metallurgical Engineers of Serbia AMES Scientific paper UDC:669.1-419:628.183=20 DYNAMIC FAILURE ANALYSIS OF LAMINATED COMPOSITE PLATES J. ESKANDARI JAM 1 and N. GARSHASBI NIA 2 1- Aerospace
More informationINTRODUCCION AL ANALISIS DE ELEMENTO FINITO (CAE / FEA)
INTRODUCCION AL ANALISIS DE ELEMENTO FINITO (CAE / FEA) Title 3 Column (full page) 2 Column What is Finite Element Analysis? 1 Column Half page The Finite Element Method The Finite Element Method (FEM)
More informationTheoretical Manual Theoretical background to the Strand7 finite element analysis system
Theoretical Manual Theoretical background to the Strand7 finite element analysis system Edition 1 January 2005 Strand7 Release 2.3 2004-2005 Strand7 Pty Limited All rights reserved Contents Preface Chapter
More informationREAL-TIME REVERB SIMULATION FOR ARBITRARY OBJECT SHAPES
REAL-TIME REVERB SIMULATION FOR ARBITRARY OBJECT SHAPES Cynthia Bruyns Maxwell University of California at Berkeley Computer Science Department Center for New Music and Audio Technologies Berkeley, California
More informationMODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD
MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD Andrzej Skrzat, Rzeszow University of Technology, Powst. Warszawy 8, Rzeszow, Poland Abstract: User-defined material models which can be used
More informationVirtual distortions applied to structural modelling and sensitivity analysis. Damage identification testing example
AMAS Workshop on Smart Materials and Structures SMART 03 (pp.313 324) Jadwisin, September 2-5, 2003 Virtual distortions applied to structural modelling and sensitivity analysis. Damage identification testing
More informationTransactions on Engineering Sciences vol 14, 1997 WIT Press, ISSN
On the Computation of Elastic Elastic Rolling Contact using Adaptive Finite Element Techniques B. Zastrau^, U. Nackenhorst*,J. Jarewski^ ^Institute of Mechanics and Informatics, Technical University Dresden,
More informationGlobal Format for Energy Momentum based Time Integration in Nonlinear Dynamics
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 204; 00: 2 Published online in Wiley InterScience (www.interscience.wiley.com). Global Format for Energy Momentum based
More informationSuppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber
Suppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber J.C. Ji, N. Zhang Faculty of Engineering, University of Technology, Sydney PO Box, Broadway,
More information