THE STATIC SUBSTRUCTURE METHOD FOR DYNAMIC ANALYSIS OF STRUCTURES. Lou Menglin* SUMMARY
|
|
- Eugenia Farmer
- 6 years ago
- Views:
Transcription
1 264 THE STATIC SUBSTRUCTURE METHOD FOR DYNAMIC ANALYSIS OF STRUCTURES Lou Mengl* SUMMARY In this paper, the static substructure method based on the Ritz vector direct superposition method is suggested for analyzg the dynamic response of structures. The advantage of this algorithm is that the computer cost can be reduced and the static analysis and the dynamic analysis of large structures can be simplified by usg the identical static substructure method. INTRODUCTION Recently, the dynamic substructure mode synthesis method (1-6) was developed for dynamic analysis of large structures. Although there are many different algorithms for the dynamic substructure method, the prciple of them all is similar. In general, the total system is divided to several substructures and a smaller eigenvalue problem is solved for each substructure. Then, by truncatg the higher modes, the lowest mode shapes or pseudo-static displacements of the substructure are used to form the Ritz vectors of the total system. The size of the total system is reduced significantly by usg these Ritz vectors. As several smaller eigenvalue problems are evaluated stead of evaluation of a large eigenvalue problem, the dynamic substructure method can reduce the computer cost and computer storage required. The dynamic substructure method is based on the modal superposition method. If the dynamic substructure method is used for dynamic analysis of large structures, the followg shortcomgs exist: 1. Solution for a set of eigenvalue problems of the substructures is still expensive, especially for large structures. 2. The error produced by truncatg higher modes of substructure cannot be avoided. The purpose of this paper is to propose a new substructure method. When this method is applied to the dynamic analysis of structures, the solution of a set of algebraic equations is required stead of a solution of an eigenvalue problem the substructure and the shortcomgs mentioned above do not exist. TWO METHODS FOR LINEAR DYNAMIC ANALYSIS OF OF STRUCTURES Dalian Institute of Technology, Dalian, P.R. Cha. The dynamic equilibrium equation for a structure, modelled by a fite number of discrete members and lumped masses, is written terms of jot displacements as [M]{u(t)}+[C]{u(t)}+[]{u(t) } = {P(t)} (1) where [M], [C] and [] are the given (nxn) mass, dampg and stiffness matrices. The time-dependent vectors {u(t)}, (u(t)} and {u(t)} are the jot accelerations, velocities and displacements respectively. If the time-varyg load can be further factored to a sum of space vectors multiplied by time functions, Eq. (1) can be written as [M]{ii(t) }+[C]{u(t) }+[]{u(t) } = {f }g(t) (2) For lear structures, the displacements vector u(t) can be approximately expressed by lear superposition of several orthogonal base vectors m (u(t)}= I {t.}q.(t) j=l 3 3 = [T]{q(t)} (3) where (q(t)} is an array of m unknown time functions and [T] is a (nxm) matrix. The substitution of Eq. 3 to Eq. 2 and the premultiplication by [T] 1^ yields [M*]{q(t) }+[C*]{<4(t) }+[*]{q(t) } [M*]=[T] T [M][T] [C*]=[T] T [C][T] = {f*}g(t) (4) [*]=[T] T [][T] (5) {f*}=[t] T {f} When m<<n, solution of Eq. (2) can be simplified significantly. BULLETIN OF THE NEW ZEALAND NATIONAL SOCIETY FOR EARTHQUAE ENGINEERING, Vol. 20 No. 4, December 1987
2 265 For lear analysis, the transformation matrix [T] can be obtaed by two methods. 1. The Model Superposition Method This is a well known method for dynamic analysis of structures. The matrix [T] is generated by solvg the followg eigenvalue problem: ( []-x [M]) lt>{0}. (3 = 1,2,...,m) (6) The matrix [T] is the assembly of the lowest eigenvectors of the system. The matrix [M*] and [*] are both diagonal. If proportional dampg is assumed, the matrix [C*] is also diagonal. Therefore, Eq. (4) is reduced to a set of uncoupled, lear, second-order, ordary differential equations. However, solvg for the eigenvalue Eq. (6) is expensive for large structures with hundreds or thousands degrees of freedom. 2. Ritz Vectors Direct Superposition Method (R-method) As the R-method turns evaluation of the eigenvalue equation to a set of evaluation of the algebraic equations, the computer cost can be reduced. The new substructure method presented this paper is based on the R-method. THE STATIC SUBSTRUCTURE METHOD In the R-method, the ma effort lies evaluatg the set of algebraic equations (7) and (9), can be written a unified form []{y}={f} (18) when {y} dicates {y. }, if) dicates {f } such as Eq. (7) and Eq. (9). A static substructure method is suggested for solvg Eq. (18) to reduce the computer cost and computer store for large structures. The prciple of the method is described as follows. Recently, E.L. Wilson et al (7) suggested a new approach to generate the matrix [T]. The ma steps of this method are described as follows: (1) Solve for the first vector {x^} []{y 1}={f 1} {x 1}=a 1{y 1> (7) (8) {f x}={f} a 1=({y 1} T [M]{y 1})"" - 5 (2) Solve for additional vectors {x^}(k=l, 2,...,m) []{y k}={f k> (9) k-1 {x*} = {y, }- I c.{x.} M orthognalized (10) {x k}=a k{x.} { f k } = [ M ] { x k-l } c k j ={y k^[mhx.} a k = ({x*} T [M]{x }r - 5 M normalized (11) (3) Reduce the stiffness matrix of the total system (12) (13) (14) []=[X] T [][X] (15) (4) Solve the generalized eigenvalue problem ([]-A.[I] ) {z.}={0} J 3 (5) Obta the matrix [T] [T]=[X][Z] (j=l,2,...,m) (16) (17) FIG. 1 - TRIANGULAR CROSS SECTION OF GRAVITY DAM AND ITS DIVISION OF SUB STRUCTURES The cross section of gravity dam shown Fig. 1 is divided to three substructures by the common boundaries AB and BC. The jots on the common boundaries AB or BC are called boundary jots and the other jots are free are called ternal jots. The variables referrg to boundary jots or ternal jots are denoted by subscripts "B" or "I" respectively. The procedure of the static method is summarized as follows: 1. Generate the mass and stiffness matrix for each substructure [ ( i ) ] = [M ( i ) ]= II BI II IB BB (19a) (19b). th, superscript "i" dicates the i substructure. 2. Condense [^] to the boundary jots to form the condensed stiffness matrix
3 266 [^] of the i*"* 1 substructure. P 2 =[]{x 2 l [ ( i ) ] = [^ } 3-[^ } ] [^ } ]" 1 [^ ) ] (20) =[](a 2 (x 2 )) 3. Assemble [^J_ to the total condensed stiffness matrix [] of the system, accordg to the order of DOF of boundary jots. []^[ ( i ) ] 4. Calculate the condensed load vector for each substructure and assemble them to form the total condensed load vector (f R ) of the system { I B i } = } { i } fb }~l ii 1 ^ ij r 1^ 1' > (21) {f } ={f ( i ) } U B ; B * 5. Calculate the displacements of the boundary jots under the load 6. Calculate the displacements of the ternal jots of each substructure <yfb-[<ih- 1 (f< i ' ) -, l Bl'l' II, l~ 1 > T After fishg step 1 to step 4 for each substructure, [][<i>]" x and [ ii } 3" 1 must rema for further use step 6 and [f 1 )] is no longer needed so that the requirements for computer storage are reduced. (22) (23) It must be noted that when Ritz vector {xj^} has been obtaed, {fii) } and {fj 1 '} Eq. (21) and Eq. (23) must be calculated as follows (24) (25) Ix.^ } and {x^ } are displacements of ternal and boundary jots of the i ^ substructure {x^. As the stiffness matrix of each substructure is not available to assemble the stiffness matrix of the total system, the reduced stiffness matrix [ ] of the total system cannot be formed from Eq. (15), but it can be generated another way. Let {p k > =[]{x k > (26) The vector (p k > is an elastic force vector of the structure makes the structure deform as {x }. Due to the recurrence relationship between the Ritz vectors {x.} (k=l,2,...,m), a recurrence relationship between (p,) exists too. Accordg to Eq. 24 and Eqs, (7,8) or Eqs. (9-11), {p^ and {p 2 > can be written as { P l )= []{ X l } = [](a 1 {y 1 }) = a l { f l } (27) = a 2 ( [ ] { y 2 } " C 21 [ ] { x l } ) =a 2 ({f 2 }-c 2 1 { P l }) (28) Similarly, the recurrence relationship can be expressed as {p v )=[]{x } k-1 - a k ( { f k } c - j l 1 kj { x j 1 > Eq. (29) dicates that {p^} is concerned with previous force vectors {pi} to {p^-^} and also with the load (fv)* After the Ritz vector {x k } is generated, the (p^) should be calculated immediately by Eq. (27) Thus, Eq. 15 becomes []=[x 1,x 2,...,x m ][p 1,p 2,...,p m ] =[X] T [P] (30) In the static substructure method,^the order of the reduced stiffness matrix [] is only dependent on the number of Ritz vectors and is dependent of the number of DOF of the boundary between substructures. APPLICATIONS Two examples are solved to illustrate the application of the static substructure method for dynamic analysis of structures. 1. Simple Shear Buildg The buildg shown Fig. 2a is simplified as a shear cantilever model with lumped mass. For simplification, it is assumed that all sprg factors k and lumped mass m are equal to unity. Four methods were adopted and compared with each other: Method A: the static substructure method presented this paper. The model was divided to four substructures as shown Fig. 2b. Three Ritz vectors were generated (m=3). The calculation was not performed by computer. Method B: the dynamic substructure method (non-terface method for substructural mode synthesis (5)). The model was divided to two substructures as shown Fig. 2c. Two modes were taken for each substructure. Four Ritz vectors were generated (m=4). Method C: the dynamic substructure method (Craig's constraed mode method (2)). The model was divided to two substructures as shown Fig. 2d. One mode was taken for each substructure. Includg one pseudostatic displacement vector, three Ritz vectors were generated (m=3). Method D: Craig's constraed mode method. The model was also divided to four substructures as shown Fig. 2b. One mode was taken for each substructure. Includg three pseudo-static displacement vectors.
4 substructure 1 substructure 2 substructure A triangular cross section of concrete gravity dam was divided to three substructures as shown Fig. 1. Its height was 100 metres and the width of its base was 80 metres. Table 2 gives a comparison of the frequencies obtaed by the static substructure method and the subspace iteration method. When usg ten Ritz vectors, the ratio of computer time needed for the two methods was about 1:1.5. TABLE 2 FIRST FIVE FREQUENCIES OF THE DAM (Hz) (a) substructure 4 (b) Mode The static substructure method the subspace iteration method CONCLUSIONS substructure 1 substructure 1 substructure 2 substructure 2 FIG. (c) (d) > SIMPLE SHEAR BUILDING AND ITS DIVISION OF SUBSTRUCTURES seven Ritz vectors were generated (m=7). The exact frequencies and the approximate frequencies of the system by usg four substructure algorithms are given Table 1. TABLE 1 - FIRST THREE FREQUENCIES OF THE BUILDING Mode exact method method method method A B C D This example shows that the static substructure method possesses higher precision and needs less computation than the dynamic substructure method as the number of Ritz vectors is the same. 2. Two-Dimensional Concrete Gravity Dam The static substructure method is suggested for analyzg the dynamic response of the structures; it possesses three advantages: 1. No algorithm error exists the analysis stage of the substructures as the static substructure method is itself an accurate algorithm. 2. An eigenvalue problem does not occur the analysis of the substructures, so the computer cost can be reduced. 3. The number of reduced eigenvalue equations for system is equal to the number of Ritz vectors and is dependent of the number of DOF of the terface between the substructures. This is very meangful analyzg the dynamic response of massive structures such as arch dams. When the Lanczos method (8) is used for the dynamic analysis of structures, its ma work is also solvg a set of algebraic equations. Consequently, the static substructure method could also be used. REFERENCES 1. W.C. Hurty, "Dynamic analysis of structural systems usg component modes", AIAA J. 3(4): 678: R.R. Craig, Jr., and M.C.C. Bampton, "Couplg of substructures for dynamic analyses", AIAA J. 6(7): S. Rub, "Improved componentmode representation for structural analysis", AIAA J. 13(8): R.R. Craig, Jr. and C-J. Chang, "Substructure couplg for dynamic analysis and testg", NASA CR R. Zhang and M. Dong, "Nonterface method for substructural mode synthesis", Computational Structural Mechanics & Applications (Cha) 1(3):
5 Z. Ma and Y. Lu, "Takg 'residual stiffness 1 and 'residual mass' to consideration for free-terface method of component mode synthesis", Jnl. of Vibration and Shock (Cha) (3): E.L. Wilson et al, "Dynamic analysis by direct superposion of Ritz vectors", Earthquake Engeerg and Structural Dynamics, 10(6): B. Nour-Omid and R.W. Clough, "Dynamic analysis of structures usg Lanczos co-ordates", Earthquake and Structural Dynamics 12(4):
Reduction 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 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 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 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 informationEffect of Mass Matrix Formulation Schemes on Dynamics of Structures
Effect of Mass Matrix Formulation Schemes on Dynamics of Structures Swapan Kumar Nandi Tata Consultancy Services GEDC, 185 LR, Chennai 600086, India Sudeep Bosu Tata Consultancy Services GEDC, 185 LR,
More informationANALYSIS OF HIGHRISE BUILDING STRUCTURE WITH SETBACK SUBJECT TO EARTHQUAKE GROUND MOTIONS
ANALYSIS OF HIGHRISE BUILDING SRUCURE WIH SEBACK SUBJEC O EARHQUAKE GROUND MOIONS 157 Xiaojun ZHANG 1 And John L MEEK SUMMARY he earthquake response behaviour of unframed highrise buildings with setbacks
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 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 informationAutomated Multi-Level Substructuring CHAPTER 4 : AMLS METHOD. Condensation. Exact condensation
Automated Multi-Level Substructuring CHAPTER 4 : AMLS METHOD Heinrich Voss voss@tu-harburg.de Hamburg University of Technology AMLS was introduced by Bennighof (1998) and was applied to huge problems of
More informationThis appendix gives you a working knowledge of the theory used to implement flexible bodies in ADAMS. The topics covered include
Appendix D Theoretical Background This appendix gives you a working knowledge of the theory used to implement flexible bodies in ADAMS. The topics covered include modal superposition component mode synthesis,
More informationMechanical Vibrations Chapter 6 Solution Methods for the Eigenvalue Problem
Mechanical Vibrations Chapter 6 Solution Methods for the Eigenvalue Problem Introduction Equations of dynamic equilibrium eigenvalue problem K x = ω M x The eigensolutions of this problem are written in
More informationCodal Provisions IS 1893 (Part 1) 2002
Abstract Codal Provisions IS 1893 (Part 1) 00 Paresh V. Patel Assistant Professor, Civil Engineering Department, Nirma Institute of Technology, Ahmedabad 38481 In this article codal provisions of IS 1893
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 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 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 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 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 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 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 informationEfficient Reduced Order Modeling of Low- to Mid-Frequency Vibration and Power Flow in Complex Structures
Efficient Reduced Order Modeling of Low- to Mid-Frequency Vibration and Power Flow in Complex Structures Yung-Chang Tan Graduate Student Research Assistant Matthew P. Castanier Assistant Research Scientist
More informationElementary Row Operations on Matrices
King Saud University September 17, 018 Table of contents 1 Definition A real matrix is a rectangular array whose entries are real numbers. These numbers are organized on rows and columns. An m n matrix
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More information18. FAST NONLINEAR ANALYSIS. The Dynamic Analysis of a Structure with a Small Number of Nonlinear Elements is Almost as Fast as a Linear Analysis
18. FAS NONLINEAR ANALYSIS he Dynamic Analysis of a Structure with a Small Number of Nonlinear Elements is Almost as Fast as a Linear Analysis 18.1 INRODUCION he response of real structures when subjected
More informationDynamic Stress Analysis of a Bus Systems
Dynamic Stress Analysis of a Bus Systems *H. S. Kim, # Y. S. Hwang, # H. S. Yoon Commercial Vehicle Engineering & Research Center Hyundai Motor Company 772-1, Changduk, Namyang, Whasung, Kyunggi-Do, Korea
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 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 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 informationSubstructure model updating through iterative minimization of modal dynamic residual
Substructure model updating through iterative minimization of modal dynamic residual Dapeng Zhu, Xinun Dong, Yang Wang* School of Civil and Environmental Eng., Georgia Inst. of Technology, Atlanta, GA
More informationComputational Stiffness Method
Computational Stiffness Method Hand calculations are central in the classical stiffness method. In that approach, the stiffness matrix is established column-by-column by setting the degrees of freedom
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 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 informationNormal modes. where. and. On the other hand, all such systems, if started in just the right way, will move in a simple way.
Chapter 9. Dynamics in 1D 9.4. Coupled motions in 1D 491 only the forces from the outside; the interaction forces cancel because they come in equal and opposite (action and reaction) pairs. So we get:
More informationDynamic Loads CE 543. Examples. Harmonic Loads
CE 543 Structural Dynamics Introduction Dynamic Loads Dynamic loads are time-varying loads. (But time-varying loads may not require dynamic analysis.) Dynamics loads can be grouped in one of the following
More 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 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 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 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 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 informationA pragmatic approach to including complex natural modes of vibration in aeroelastic analysis
A pragmatic approach to including complex natural modes of vibration in aeroelastic analysis International Aerospace Symposium of South Africa 14 to 16 September, 215 Stellenbosch, South Africa Louw van
More informationNONLINEAR DYNAMIC ANALYSIS OF COMPLEX STRUCTURES
EARTHQUAKE ENGINEERING AND STRUCTURAL DYNAMICS, VOL. 1, 241-252 (1973) NONLINEAR DYNAMIC ANALYSIS OF COMPLEX STRUCTURES E. L. WILSON~ University of California, Berkeley, California I. FARHOOMAND$ K. J.
More informationMultilevel Methods for Eigenspace Computations in Structural Dynamics
Multilevel Methods for Eigenspace Computations in Structural Dynamics Ulrich Hetmaniuk & Rich Lehoucq Sandia National Laboratories, Computational Math and Algorithms, Albuquerque, NM Joint work with Peter
More informationSoftware Verification
EXAMPLE 6-003 LINK GAP ELEMENT PROBLEM DESCRIPTION This example uses a single-bay, single-story rigid frame to test the gap link element. This link element carries compression loads only; it has zero stiffness
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 informationEFFECTIVE MODAL MASS & MODAL PARTICIPATION FACTORS Revision F
EFFECTIVE MODA MASS & MODA PARTICIPATION FACTORS Revision F By Tom Irvine Email: tomirvine@aol.com March 9, 1 Introduction The effective modal mass provides a method for judging the significance of a vibration
More informationAPPLICATION OF RESPONSE SPECTRUM METHOD TO PASSIVELY DAMPED DOME STRUCTURE WITH HIGH DAMPING AND HIGH FREQUENCY MODES
3 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 4 Paper No. 5 APPLICATION OF RESPONSE SPECTRUM METHOD TO PASSIVELY DAMPED DOME STRUCTURE WITH HIGH DAMPING AND HIGH FREQUENCY
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 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 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 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 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 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 informationhave invested in supercomputer systems, which have cost up to tens of millions of dollars each. Over the past year or so, however, the future of vecto
MEETING THE NVH COMPUTATIONAL CHALLENGE: AUTOMATED MULTI-LEVEL SUBSTRUCTURING J. K. Bennighof, M. F. Kaplan, y M. B. Muller, y and M. Kim y Department of Aerospace Engineering & Engineering Mechanics The
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 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 informationNonlinear Analysis of Reinforced Concrete Bridges under Earthquakes
6 th International Conference on Advances in Experimental Structural Engineering 11 th International Workshop on Advanced Smart Materials and Smart Structures Technology August 1-2, 2015, University of
More informationon the figure. Someone has suggested that, in terms of the degrees of freedom x1 and M. Note that if you think the given 1.2
1) A two-story building frame is shown below. The mass of the frame is assumed to be lumped at the floor levels and the floor slabs are considered rigid. The floor masses and the story stiffnesses are
More informationPreliminary Examination - Dynamics
Name: University of California, Berkeley Fall Semester, 2018 Problem 1 (30% weight) Preliminary Examination - Dynamics An undamped SDOF system with mass m and stiffness k is initially at rest and is then
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 informationStress analysis of a stepped bar
Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has cross-sectional areas of A ) and A ) over the lengths l ) and l ), respectively.
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 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 informationNONLINEAR STRUCTURAL DYNAMICS USING FE METHODS
NONLINEAR STRUCTURAL DYNAMICS USING FE METHODS Nonlinear Structural Dynamics Using FE Methods emphasizes fundamental mechanics principles and outlines a modern approach to understanding structural dynamics.
More informationResponse Analysis for Multi Support Earthquake Excitation
Chapter 5 Response Analysis for Multi Support Earthquake Excitation 5.1 Introduction It is very important to perform the dynamic analysis for the structure subjected to random/dynamic loadings. The dynamic
More informationME 475 Modal Analysis of a Tapered Beam
ME 475 Modal Analysis of a Tapered Beam Objectives: 1. To find the natural frequencies and mode shapes of a tapered beam using FEA.. To compare the FE solution to analytical solutions of the vibratory
More informationDynamics of Ocean Structures Prof. Dr. Srinivasan Chandrasekaran Department of Ocean Engineering Indian Institute of Technology, Madras
Dynamics of Ocean Structures Prof. Dr. Srinivasan Chandrasekaran Department of Ocean Engineering Indian Institute of Technology, Madras Lecture 25 Continuous System In the last class, in this, we will
More informationNumerical Simulation of Vibroacoustic Problems. by the Modal Synthesis
Adv. Theor. Appl. Mech., Vol. 6, 2013, no. 1, 1-12 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/atam.2013.21237 Numerical Simulation of Vibroacoustic Problems by the Modal Synthesis M. Mansouri
More informationLast Time. Social Network Graphs Betweenness. Graph Laplacian. Girvan-Newman Algorithm. Spectral Bisection
Eigenvalue Problems Last Time Social Network Graphs Betweenness Girvan-Newman Algorithm Graph Laplacian Spectral Bisection λ 2, w 2 Today Small deviation into eigenvalue problems Formulation Standard eigenvalue
More informationVIBRATION PROBLEMS IN ENGINEERING
VIBRATION PROBLEMS IN ENGINEERING FIFTH EDITION W. WEAVER, JR. Professor Emeritus of Structural Engineering The Late S. P. TIMOSHENKO Professor Emeritus of Engineering Mechanics The Late D. H. YOUNG Professor
More informationINELASTIC RESPONSES OF LONG BRIDGES TO ASYNCHRONOUS SEISMIC INPUTS
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 24 Paper No. 638 INELASTIC RESPONSES OF LONG BRIDGES TO ASYNCHRONOUS SEISMIC INPUTS Jiachen WANG 1, Athol CARR 1, Nigel
More informationAdaptation of the Lanczos Algorithm for the Solution of Buckling Eigenvalue Problems
Original Article Adaptation of the Lanczos Algorithm for the Solution of Buckling Eigenvalue Problems Abstract An adaptation of the conventional Lanczos algorithm is proposed to solve the general symmetric
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 informationstiffness to the system stiffness matrix. The nondimensional parameter i is introduced to allow the modeling of damage in the ith substructure. A subs
A BAYESIAN PROBABILISTIC DAMAGE DETECTION USING LOAD-DEPENDENT RIT VECTORS HOON SOHN Λ and KINCHO H. LAW y Department of Civil and Environmental Engineering, Stanford University, Stanford, CA 9435-42,U.S.A.
More informationEXPERIMENTAL CONFIRMATION FOR THE VALIDITY OF RITZ METHOD IN STRUCTURAL DYNAMIC ANALYSIS
JOURNAL OF THEORETICAL AND APPLIED MECHANICS 52, 4, pp. 981-993, Warsaw 2014 EXPERIMENTAL CONFIRMATION FOR THE VALIDITY OF RITZ METHOD IN STRUCTURAL DYNAMIC ANALYSIS Gultekin Aktas, Abdulhalim Karasin
More informationNumerical Solution of Equation of Motion
Class Notes: Earthquake Engineering, Ahmed Elgamal, September 25, 2001 (DRAFT) Numerical Solution of Equation of Motion Average Acceleration Method (Trapezoidal method) m a + c v + k d = f (t) In the above
More 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 informationSoftware Verification
SAP000 EXAMPLE 6-00 LINK LINEAR LINK WITH RAMP LOADING PROBLEM DESCRIPTION In this example a ramp load is applied to an undamped single degree of freedom structure. The ramp loading has a finite rise time
More informationStochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore
Stochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore Lecture No. # 32 Probabilistic Methods in Earthquake Engineering-1 (Refer Slide
More informationSubstructuring using Impulse Response Functions for Impact Analysis
Proceedings of the IMAC-XXVIII February 1 4, 21, Jacksonville, Florida USA 21 Society for Experimental Mechanics Inc. Substructuring using Impulse Response Functions for Impact Analysis Daniel J. Rixen
More informationLarge-scale eigenvalue problems
ELE 538B: Mathematics of High-Dimensional Data Large-scale eigenvalue problems Yuxin Chen Princeton University, Fall 208 Outline Power method Lanczos algorithm Eigenvalue problems 4-2 Eigendecomposition
More informationPart D: Frames and Plates
Part D: Frames and Plates Plane Frames and Thin Plates A Beam with General Boundary Conditions The Stiffness Method Thin Plates Initial Imperfections The Ritz and Finite Element Approaches A Beam with
More informationCRAIG-BAMPTON METHOD FOR A TWO COMPONENT SYSTEM Revision C
CRAIG-BAMPON MEHOD FOR A WO COMPONEN SYSEM Revision C By om Irvine Email: tom@vibrationdata.com May, 03 Introduction he Craig-Bampton method is method for reducing the size of a finite element model, particularly
More informationPart 6: Dynamic design analysis
Part 6: Dynamic design analysis BuildSoft nv All rights reserved. No part of this document may be reproduced or transmitted in any form or by any means, electronic or manual, for any purpose, without written
More informationPreliminaries: Beam Deflections Virtual Work
Preliminaries: Beam eflections Virtual Work There are several methods available to calculate deformations (displacements and rotations) in beams. They include: Formulating moment equations and then integrating
More information3.4 Analysis for lateral loads
3.4 Analysis for lateral loads 3.4.1 Braced frames In this section, simple hand methods for the analysis of statically determinate or certain low-redundant braced structures is reviewed. Member Force Analysis
More informationx 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7
Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)
More informationMODAL PARAMETER TRACKING FOR SHAPE-CHANGING GEOMETRIC OBJECTS
MODAL PARAMETER TRACKING FOR SHAPE-CHANGING GEOMETRIC OBJECTS Cynthia Bruyns Maxwell University of California at Berkeley Computer Science Department, CNMAT Berkeley, California USA cbruyns@cs.berkeley.edu
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 informationParametric Identification of a Cable-stayed Bridge using Substructure Approach
Parametric Identification of a Cable-stayed Bridge using Substructure Approach *Hongwei Huang 1), Yaohua Yang 2) and Limin Sun 3) 1),3) State Key Laboratory for Disaster Reduction in Civil Engineering,
More informationResponse Spectrum Analysis Methods for Spatial Structures
8 September October 009, Universidad Politecnica de Valencia, Spain Alberto DOMINGO and Carlos LAZARO (eds. Response Spectrum Analysis Methods for Spatial Structures Yongfeng LUO *, Lei WANG Department
More informationMODELLING DAMAGE AND PROGRESSIVE COLLAPSE OF FRAMES USING A GAUSSIAN SPRINGS BASED APPLIED ELEMENT METHOD
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 018, Glasgow, UK MODELLING DAMAGE AND PROGRESSIVE COLLAPSE OF FRAMES
More informationMath Camp II. Basic Linear Algebra. Yiqing Xu. Aug 26, 2014 MIT
Math Camp II Basic Linear Algebra Yiqing Xu MIT Aug 26, 2014 1 Solving Systems of Linear Equations 2 Vectors and Vector Spaces 3 Matrices 4 Least Squares Systems of Linear Equations Definition A linear
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 informationIntroduction to Finite Element Method. Dr. Aamer Haque
Introduction to Finite Element Method 4 th Order Beam Equation Dr. Aamer Haque http://math.iit.edu/~ahaque6 ahaque7@iit.edu Illinois Institute of Technology July 1, 009 Outline Euler-Bernoulli Beams Assumptions
More informationBack Matter Index The McGraw Hill Companies, 2004
INDEX A Absolute viscosity, 294 Active zone, 468 Adjoint, 452 Admissible functions, 132 Air, 294 ALGOR, 12 Amplitude, 389, 391 Amplitude ratio, 396 ANSYS, 12 Applications fluid mechanics, 293 326. See
More informationResidual Force Equations
3 Residual Force Equations NFEM Ch 3 Slide 1 Total Force Residual Equation Vector form r(u,λ) = 0 r = total force residual vector u = state vector with displacement DOF Λ = array of control parameters
More informationInvestigation of Stress Analysis in Gear Simulations Using Different Model Order Reduction Techniques
Master Thesis DIPL-MSC-215 Investigation of Stress Analysis in Gear Simulations Using Different Model Order Reduction Techniques from Iosu Oloriz del Rosario Advisor: Prof. Dr. Ing. Prof. E.h. P. Eberhard
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 informationLinear Algebra. Linear Equations and Matrices. Copyright 2005, W.R. Winfrey
Copyright 2005, W.R. Winfrey Topics Preliminaries Systems of Linear Equations Matrices Algebraic Properties of Matrix Operations Special Types of Matrices and Partitioned Matrices Matrix Transformations
More informationContents i. Contents
Contents i Contents 7 SEISMIC LOADS Commentary 7. Estimation of Seismic Loads.............................. 7.. Seismic load and design earthquake motion.................. 4 7..2 Idealization of building
More information