Quintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation

Size: px
Start display at page:

Download "Quintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation"

Transcription

1 General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation spring has a known displacement,, which is almost always zero. The differential equation of equilibrium of a beam, of flexural stiffness EI, resting on a foundation, of stiffness k per unit length, with a transverse load per unit length of w, subjected to a tensile axial load N is: 0. (1) where v(x) is the transverse displacement of the beam. The beam equation is a fourth order ordinary differential equation. Therefore it will generally need four boundary conditions. Figure 1 A beam-column on an elastic foundation The distributed load per unit length can include point transverse shear loads, V, by using the Dirac Delta distribution. Likewise, employing a doublet distribution in defining w(x) allows for the inclusion of point couples, or moments, M. Related physical quantities are the slope,, the bending moment,, and the transverse shear force,. The sign conventions are that the position, x, is positive to the right, the deflection, v, point forces, P, and the line load, w, are positive in the y-direction (upward), and the slopes and moments are positive in the counterclockwise direction. Since the foundation springs are restrained, there can be no rigid body motions of the beam. In other words, the assembled equations should never be singular and it is acceptable to just have non-essential boundary conditions applied to the beam. A simple Winkler foundation model like this one can push or pull on the beam as needed and no gaps can occur. Engineers designing beams are usually interested in the last two quantities since they define the stress levels and the material failure criteria. The exact solution of the homogeneous (w==0) general form given above is given in terms of hyperbolic sines and cosines. Based on Tong s Theorem, exact solutions at the nodes are obtained if such functions are used in a finite element model. Advanced elements of that type have been applied with excellent results. To apply the Galerkin weak form the governing ODE is multiplied by v(x) and the integral over the length of the beam is set to zero. The highest derivative term is always present, and needs to be integrated twice by parts to reduce the inter-element continuity requirement. That integral becomes Copyright J.E. Akin. All rights reserved. Page 1 of 11

2 which brings the non-essential boundary conditions into the integral form. The first term is the product of the displacement and the point shear force at both ends, while the second term is the product of the slope and the point moment at both ends. Thus, the contributions from the highest derivative term in the ODE are Therefore, the integral form contains second derivatives ( as its highest derivative term. The presence of second derivatives in the integral form means that calculus requires the elements to have inter-element continuity of the deflection and the slope (v(x) and x ). Such elements are said to have an inter-element continuity of C. Therefore, each shared beam node must have two degrees of freedom (dof), at least. Now that there is more than one dof per node the scatter/gather operations previously used for axial models have become more generalized. The usual convention for numbering the unknowns is to count all of the unknowns at a node before moving on to another node. Therefore, the displacements are associated with the odd number of rows and columns in a matrix, while the slopes contributions are occurring in the even ones. For previous one-dimensional scalar unknowns scatter simply meant adding a scalar term at any node into the system equilibrium matrix and source vector. Now, that changes to adding a two by two sub-matrix from each element node into the system square matrix and a two by one column matrix into the system source vector. The prior algorithms for doing the scatter and gather still work because they had the number of degrees of freedom per node (n_g) as an input argument. It is just for hand solutions that you need to clearly understand this change. Figure 2 graphically illustrates the sub-matrices scatter operations for a line element with two unknowns per node. The general numbering scheme that relates local element equation numbers to the system equation numbers refer to Figure 3. As a detailed example of the change in scatter operations, consider a long pipeline with hundreds of three node beam element that are numbered in a random fashion. Figure 4 shows typical element connectivity data for the hypothetical example. Since there are two dof per node there are six local equations associated with each quintic beam element. They have to be scattered to a much much larger set of system equilibrium equations. Figure 5 illustrates how all of the local generalized force vector terms are scattered, as well as showing the same information for six typical terms selected for the element local stiffness matrix. Copyright J.E. Akin. All rights reserved. Page 2 of 11

3 Figure 2 Scattering sub-matrices for line elements with 2 dof per node Copyright J.E. Akin. All rights reserved. Page 3 of 11

4 Figure 3 Calculating equation numbers for general scatter operations Figure 4 Example data for gather/scatter of element number 21 Copyright J.E. Akin. All rights reserved. Page 4 of 11

5 Figure 5 Scattering typical matrix terms from example element number 21 The presence of the second derivatives in the integral form also means that the essential boundary conditions (EBC) will specify or or both at a point. The non-essential boundary conditions are that the shear,, and/or the moment,, is given at a point. When the deflection is given, then the shear force is a reaction quantity. When the slope is specified, then the moment is a reaction. If an axial force is present, then this is usually called a Beam-Column model. If the axial load is an unknown constant, then it is a buckling load that has to be found by an eigen-solution rather than by solving a linear algebraic system. Beams on elastic foundations are fairly common (like pipelines or rail lines). However, most introductions to beams assume no axial load, nor any foundation support. Then, the above general beam equation reduces to that covered in an introduction to solid mechanics: EI wx. (2) The homogeneous solution of this equation is simply a cubic polynomial. A cubic polynomial has four constants. Calculus continuity requires the elements share their slope and deflection values. Thus, classic beam elements approximated by at least a cubic polynomial having C continuity can give nodally exact deflections and slopes. One-dimensional polynomial interpolation with C or C continuity usually is based on Hermite interpolation polynomials. The classic beam element is often said to be a Hermite cubic. Some typical Hermite polynomials are given in the Appendix. A line element with only two nodes will provide the four constants needed to define a cubic polynomial. Such elements are referred to as classic beam elements. But since designers are Copyright J.E. Akin. All rights reserved. Page 5 of 11

6 interested in the second and third derivatives of the solution, classic beam elements give very poor design results unless a large number of them are used. For a beam with a cross-sectional moment of inertia, I, length, L, a depth of h, and material with an elastic modulus of E and a coefficient of thermal expansion of α, the corresponding matrices for the equilibrium of a single classic beam element, are , (3) where and are the load per unit length at the first (left) and second end, respectively, and is the temperature increase, over h, from the top to the bottom of the beam. Also, the and terms are point shear forces and point moments applied to the nodes as external loadings and/or unknown reactions. In matrix symbol notation this equation is written as. (4) The classic beam element has a deflection given by (for rx L): (5) Its moment and shear are linear and constant along the beam length, respectively. Specifically, they are: / (6) /. (7) The Hermite family of interpolation polynomials can be increased in order by increasing the number of nodes and/or by increasing the continuity level at each node. Many exact solutions for beams not on an elastic foundation are tabulated in typical structural design manuals. They show that for the most common load and support conditions the exact solution is a third, fourth or fifth degree polynomial. Thus, a fifth degree polynomial will generally give a very accurate or exact beam solution with only a few elements, perhaps with only one element. However, if a foundation is present polynomial approximations at the nodes will no longer be exact. Here, a more accurate three-node C beam element will be presented. This will be created by adding a third (middle) node to the element. That adds two more degrees of freedom and raises the polynomial to a fifth degree (a quintic polynomial). For simplicity, all of the coefficients in Eq. (1) are taken as constants evaluated at the center of the element length. That allows the corresponding finite element matrices of the fifth-degree element to be expressed in closed form in terms of the element deflection and slope at each of its three nodes [v θ v θ v θ. The element interpolation relations and the corresponding slope, moment and shear distributions are: Copyright J.E. Akin. All rights reserved. Page 6 of 11

7 (8) / / / (9) / / / / (10) 396 1,632 1,440 / / ,248 1,350 / / (11) The 6 by 6 element matrix equations of equilibrium of a single element are: (12) where { } T = [v, θ,v,θ,v,θ ] are the generalized nodal displacements and rotations, and the are the generalized resultant nodal loading forces and moments. The terms are stiffness matrices arising from the first three terms in the differential equation. Usually, only the first is nonzero. This system is singular until enough boundary conditions are applied to prevent rigid body translations and rotations. Since this element has three nodes, the loading per unit length can be input at each of those nodes. The resultant element external loading vectors due to point forces,, and couples,, and the quadratic distributed load values,, at each node (or constant, or linear 2 values) are,, 1,995 1, , (13) ,540 1, where the are the six Hermite beam interpolations, and the are the three quadratic Lagrangian interpolations for the distributed loads. The integral of their product forms a rectangular six by three array to convert distributed loads to concentrated shears and moments at the beam nodes. Likewise, it can be shown that a temperature difference through the depth of a beam causes only end kmoment loadings given by A constant transverse load reduces the resultant load and moment vector to,3, , , , ( ) Copyright J.E. Akin. All rights reserved. Page 7 of 11

8 The symmetric flexural stiffness matrix for the three noded quintic element is 5,092 1,1383,584 1, , ,168 1, , , ,9201, , ,280 1, ,920 5,092 1, , (14) The symmetric beam-column (or geometric stiffness) matrix due to any axial load is, 15, ,824 2,160 1, , , , , ,304 2, , ,824 2,160 15, (15) The element foundation stiffness is, 2, , , (16) Since it is common for beams to have piecewise constant properties, the classic beam element and the more accurate quintic beam element can be programmed in closed form. Using either the classic or quintic beam element some beam analysis problems can be solved in closed form. As an example of the three-node beam, consider a fixed-fixed beam with a constant line load. The essential boundary conditions are 0, 0, 0, 0. The external point force and moment at center node 2 are zero ( 0, 0). The middle two rows define the remaining unknown center point generalized displacements: 0 7, , ,7, Multiplying by the inverse of the square matrix gives the middle node solutions: 351 7, , ,280 7, which are exact. The slope was expected to be zero due to symmetry. It could have been used as a boundary condition to compute from one equation. Here, the reactions on the left are found from the first two rows of the equilibrium equations (since all the displacements are now known): Copyright J.E. Akin. All rights reserved. Page 8 of 11

9 35 3,584 1, ,700 3, Likewise, at the right end, utilizing the last two rows of the equilibrium equations 1 6. Note that the net resultant force is wl, which is equal and opposite to the applied transverse load. The net external moment is zero. Usually static equilibrium is taught in undergraduate classes with Newton s Laws. You can use them to verify that the computed reactions do indeed satisfy that the sum of the forces is zero, and that the sum of the moments, taken at any reference point, is zero. Since a 5-th degree element was used this problem is exact everywhere. That includes the moment and shear diagrams. The exact shear diagrams of can have a discontinuity (due to point loads P) and the exact moment diagram of can have a discontinuity (due to point moments M). The mesh should be constricted such that any point source occurs at an interface node between elements. Those items can be discontinuous at element interfaces because the calculus requirement for splitting integrals on says that v and v must be continuous there. That explains why a two node Hermite interpolation was not chosen to create a quintic beam. That choice could be made only if the beam was not allowed to have point sources. Inside any element the solution is continuous and thus is at any non-interface node. In other words, inside any standard element you can calculate an infinite number of continuous derivatives (but the vast majority is zero). If you were to apply a point source at an interior node then the second and third derivatives would still be continuous and the shear and moment diagrams in that element would be wrong. Next consider the involvement of the foundation in the post-processing needed to recover the reactions on each beam element. Before the deflections of a beam on an elastic foundation (BOEF) are computed, a typical element free body diagram would be similar to Figure 1. As before, any external point load, P, and/or any point couple, M, should be placed at an element interface node and they go directly into the first or second equilibrium row of the node at that point. A distributed line load is integrated to become a consistent load vector. Those two effects are detailed in Eq. 13. Both the beam stiffness and the foundation stiffness resist the external applied loads. The element equilibrium equation (for no axial force) is Likewise, the assembled equations have a similar form and are solved for the system displacements, v. Returning to the element to get the beam member reactions, the now known element degrees of freedom,, are gathered. The free body diagram of the element now has a significant change Copyright J.E. Akin. All rights reserved. Page 9 of 11

10 because the foundation is now applying a known pressure that opposes the displacement of the beam. That pressure,, is a polynomial of the same degree as the assumed displacements. That state is shown in Figure 6, which defines the sought element reaction vector,. Figure 6 Beam element reaction vector, including foundation pressure Now, there is a new consistent load vector to evaluate coming from the foundation pressure:. The foundation stiffness caused the new load, but that load and others must be resisted by only the beam stiffness. Now the beam member equilibrium equations define the reactions on the nodes of the beam, : or Copyright J.E. Akin. All rights reserved. Page 10 of 11

11 Appendix Figure 7 The first three Hermite interpolation options for elements with two nodes Copyright J.E. Akin. All rights reserved. Page 11 of 11

General elastic beam with an elastic foundation

General elastic beam with an elastic foundation General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation

More information

10. Applications of 1-D Hermite elements

10. Applications of 1-D Hermite elements 10. Applications of 1-D Hermite elements... 1 10.1 Introduction... 1 10.2 General case fourth-order beam equation... 3 10.3 Integral form... 5 10.4 Element Arrays... 7 10.5 C1 Element models... 8 10.6

More information

Contents as of 12/8/2017. Preface. 1. Overview...1

Contents as of 12/8/2017. Preface. 1. Overview...1 Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite

More information

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 11

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 11 Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Module - 01 Lecture - 11 Last class, what we did is, we looked at a method called superposition

More information

k 21 k 22 k 23 k 24 k 31 k 32 k 33 k 34 k 41 k 42 k 43 k 44

k 21 k 22 k 23 k 24 k 31 k 32 k 33 k 34 k 41 k 42 k 43 k 44 CE 6 ab Beam Analysis by the Direct Stiffness Method Beam Element Stiffness Matrix in ocal Coordinates Consider an inclined bending member of moment of inertia I and modulus of elasticity E subjected shear

More information

COPYRIGHTED MATERIAL. Index

COPYRIGHTED 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 information

Discontinuous Distributions in Mechanics of Materials

Discontinuous Distributions in Mechanics of Materials Discontinuous Distributions in Mechanics of Materials J.E. Akin, Rice University 1. Introduction The study of the mechanics of materials continues to change slowly. The student needs to learn about software

More information

3. Numerical integration

3. Numerical integration 3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...

More information

CIVL 7/8117 Chapter 4 - Development of Beam Equations - Part 2 1/34. Chapter 4b Development of Beam Equations. Learning Objectives

CIVL 7/8117 Chapter 4 - Development of Beam Equations - Part 2 1/34. Chapter 4b Development of Beam Equations. Learning Objectives CIV 7/87 Chapter 4 - Development of Beam Equations - Part /4 Chapter 4b Development of Beam Equations earning Objectives To introduce the work-equivalence method for replacing distributed loading by a

More information

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13 Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras (Refer Slide Time: 00:25) Module - 01 Lecture - 13 In the last class, we have seen how

More information

BEAM DEFLECTION THE ELASTIC CURVE

BEAM DEFLECTION THE ELASTIC CURVE BEAM DEFLECTION Samantha Ramirez THE ELASTIC CURVE The deflection diagram of the longitudinal axis that passes through the centroid of each cross-sectional area of a beam. Supports that apply a moment

More information

Methods of Analysis. Force or Flexibility Method

Methods of Analysis. Force or Flexibility Method INTRODUCTION: The structural analysis is a mathematical process by which the response of a structure to specified loads is determined. This response is measured by determining the internal forces or stresses

More information

Chapter 5 Structural Elements: The truss & beam elements

Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations

More information

FLEXIBILITY METHOD FOR INDETERMINATE FRAMES

FLEXIBILITY METHOD FOR INDETERMINATE FRAMES UNIT - I FLEXIBILITY METHOD FOR INDETERMINATE FRAMES 1. What is meant by indeterminate structures? Structures that do not satisfy the conditions of equilibrium are called indeterminate structure. These

More information

Transactions of the VŠB Technical University of Ostrava, Mechanical Series. article No Roland JANČO *

Transactions of the VŠB Technical University of Ostrava, Mechanical Series. article No Roland JANČO * Transactions of the VŠB Technical University of Ostrava, Mechanical Series No. 1, 013, vol. LIX article No. 1930 Roland JANČO * NUMERICAL AND EXACT SOLUTION OF BUCKLING LOAD FOR BEAM ON ELASTIC FOUNDATION

More information

Introduction to Finite Element Method

Introduction to Finite Element Method Introduction to Finite Element Method Dr. Rakesh K Kapania Aerospace and Ocean Engineering Department Virginia Polytechnic Institute and State University, Blacksburg, VA AOE 524, Vehicle Structures Summer,

More information

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method Module 2 Analysis of Statically Indeterminate Structures by the Matrix Force Method Lesson 8 The Force Method of Analysis: Beams Instructional Objectives After reading this chapter the student will be

More information

CIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen

CIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall

More information

Computational Stiffness Method

Computational 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 information

Converting Plane Stress Statics to 2D Natural Frequencies, changes in red

Converting Plane Stress Statics to 2D Natural Frequencies, changes in red The following illustrates typical modifications for converting any plane stress static formulation into a plane stress natural frequency and mode shape calculation. The changes and additions to a prior

More information

Structural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.

Structural 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 information

Chapter 3 Higher Order Linear ODEs

Chapter 3 Higher Order Linear ODEs Chapter 3 Higher Order Linear ODEs Advanced Engineering Mathematics Wei-Ta Chu National Chung Cheng University wtchu@cs.ccu.edu.tw 1 2 3.1 Homogeneous Linear ODEs 3 Homogeneous Linear ODEs An ODE is of

More information

Section 6: PRISMATIC BEAMS. Beam Theory

Section 6: PRISMATIC BEAMS. Beam Theory Beam Theory There are two types of beam theory aailable to craft beam element formulations from. They are Bernoulli-Euler beam theory Timoshenko beam theory One learns the details of Bernoulli-Euler beam

More information

Structural 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 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 information

Lecture 7: The Beam Element Equations.

Lecture 7: The Beam Element Equations. 4.1 Beam Stiffness. A Beam: A long slender structural component generally subjected to transverse loading that produces significant bending effects as opposed to twisting or axial effects. MECH 40: Finite

More information

9.1 Introduction to bifurcation of equilibrium and structural

9.1 Introduction to bifurcation of equilibrium and structural Module 9 Stability and Buckling Readings: BC Ch 14 earning Objectives Understand the basic concept of structural instability and bifurcation of equilibrium. Derive the basic buckling load of beams subject

More information

2. Polynomial interpolation

2. Polynomial interpolation 2. Polynomial interpolation Contents 2. POLYNOMIAL INTERPOLATION... 1 2.1 TYPES OF INTERPOLATION... 1 2.2 LAGRANGE ONE-DIMENSIONAL INTERPOLATION... 2 2.3 NATURAL COORDINATES*... 15 2.4 HERMITE ONE-DIMENSIONAL

More information

Stress analysis of a stepped bar

Stress 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 information

UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES

UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the

More information

FINITE GRID SOLUTION FOR NON-RECTANGULAR PLATES

FINITE GRID SOLUTION FOR NON-RECTANGULAR PLATES th International Conference on Earthquake Geotechnical Engineering June 5-8, 7 Paper No. 11 FINITE GRID SOLUTION FOR NON-RECTANGULAR PLATES A.Halim KARAŞĐN 1, Polat GÜLKAN ABSTRACT Plates on elastic foundations

More information

M5 Simple Beam Theory (continued)

M5 Simple Beam Theory (continued) M5 Simple Beam Theory (continued) Reading: Crandall, Dahl and Lardner 7.-7.6 In the previous lecture we had reached the point of obtaining 5 equations, 5 unknowns by application of equations of elasticity

More information

Chapter 2: Rigid Bar Supported by Two Buckled Struts under Axial, Harmonic, Displacement Excitation..14

Chapter 2: Rigid Bar Supported by Two Buckled Struts under Axial, Harmonic, Displacement Excitation..14 Table of Contents Chapter 1: Research Objectives and Literature Review..1 1.1 Introduction...1 1.2 Literature Review......3 1.2.1 Describing Vibration......3 1.2.2 Vibration Isolation.....6 1.2.2.1 Overview.

More information

Chapter 11. Displacement Method of Analysis Slope Deflection Method

Chapter 11. Displacement Method of Analysis Slope Deflection Method Chapter 11 Displacement ethod of Analysis Slope Deflection ethod Displacement ethod of Analysis Two main methods of analyzing indeterminate structure Force method The method of consistent deformations

More information

Lecture 8: Assembly of beam elements.

Lecture 8: Assembly of beam elements. ecture 8: Assembly of beam elements. 4. Example of Assemblage of Beam Stiffness Matrices. Place nodes at the load application points. Assembling the two sets of element equations (note the common elemental

More information

Institute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I

Institute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I Institute of Structural Engineering Page 1 Chapter 2 The Direct Stiffness Method Institute of Structural Engineering Page 2 Direct Stiffness Method (DSM) Computational method for structural analysis Matrix

More information

Lecture 11: The Stiffness Method. Introduction

Lecture 11: The Stiffness Method. Introduction Introduction Although the mathematical formulation of the flexibility and stiffness methods are similar, the physical concepts involved are different. We found that in the flexibility method, the unknowns

More information

Vibration of Thin Beams by PIM and RPIM methods. *B. Kanber¹, and O. M. Tufik 1

Vibration of Thin Beams by PIM and RPIM methods. *B. Kanber¹, and O. M. Tufik 1 APCOM & ISCM -4 th December, 23, Singapore Vibration of Thin Beams by PIM and RPIM methods *B. Kanber¹, and O. M. Tufik Mechanical Engineering Department, University of Gaziantep, Turkey. *Corresponding

More information

Advanced Structural Analysis Prof. Devdas Menon Department of Civil Engineering Indian Institute of Technology, Madras

Advanced Structural Analysis Prof. Devdas Menon Department of Civil Engineering Indian Institute of Technology, Madras Advanced Structural Analysis Prof. Devdas Menon Department of Civil Engineering Indian Institute of Technology, Madras Module No. # 5.1 Lecture No. # 27 Matrix Analysis of Beams and Grids Good morning,

More information

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 17

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 17 Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil engineering Indian Institute of Technology, Madras Module - 01 Lecture - 17 In the last class, we were looking at this general one-dimensional

More information

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).

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). 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 information

Interpolation Functions for General Element Formulation

Interpolation Functions for General Element Formulation CHPTER 6 Interpolation Functions 6.1 INTRODUCTION The structural elements introduced in the previous chapters were formulated on the basis of known principles from elementary strength of materials theory.

More information

CHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES

CHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES CHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES 14.1 GENERAL REMARKS In structures where dominant loading is usually static, the most common cause of the collapse is a buckling failure. Buckling may

More information

Bilinear Quadrilateral (Q4): CQUAD4 in GENESIS

Bilinear Quadrilateral (Q4): CQUAD4 in GENESIS Bilinear Quadrilateral (Q4): CQUAD4 in GENESIS The Q4 element has four nodes and eight nodal dof. The shape can be any quadrilateral; we ll concentrate on a rectangle now. The displacement field in terms

More information

Chapter 4 Deflection and Stiffness

Chapter 4 Deflection and Stiffness Chapter 4 Deflection and Stiffness Asst. Prof. Dr. Supakit Rooppakhun Chapter Outline Deflection and Stiffness 4-1 Spring Rates 4-2 Tension, Compression, and Torsion 4-3 Deflection Due to Bending 4-4 Beam

More information

Eigenvalues of Trusses and Beams Using the Accurate Element Method

Eigenvalues of Trusses and Beams Using the Accurate Element Method Eigenvalues of russes and Beams Using the Accurate Element Method Maty Blumenfeld Department of Strength of Materials Universitatea Politehnica Bucharest, Romania Paul Cizmas Department of Aerospace Engineering

More information

FEA CODE WITH MATLAB. Finite Element Analysis of an Arch ME 5657 FINITE ELEMENT METHOD. Submitted by: ALPAY BURAK DEMIRYUREK

FEA CODE WITH MATLAB. Finite Element Analysis of an Arch ME 5657 FINITE ELEMENT METHOD. Submitted by: ALPAY BURAK DEMIRYUREK FEA CODE WITH MATAB Finite Element Analysis of an Arch ME 5657 FINITE EEMENT METHOD Submitted by: APAY BURAK DEMIRYUREK This report summarizes the finite element analysis of an arch-beam with using matlab.

More information

LATERAL STABILITY OF BEAMS WITH ELASTIC END RESTRAINTS

LATERAL STABILITY OF BEAMS WITH ELASTIC END RESTRAINTS LATERAL STABILITY OF BEAMS WITH ELASTIC END RESTRAINTS By John J. Zahn, 1 M. ASCE ABSTRACT: In the analysis of the lateral buckling of simply supported beams, the ends are assumed to be rigidly restrained

More information

Chapter Objectives. Copyright 2011 Pearson Education South Asia Pte Ltd

Chapter Objectives. Copyright 2011 Pearson Education South Asia Pte Ltd Chapter Objectives To generalize the procedure by formulating equations that can be plotted so that they describe the internal shear and moment throughout a member. To use the relations between distributed

More information

ENGI Second Order Linear ODEs Page Second Order Linear Ordinary Differential Equations

ENGI Second Order Linear ODEs Page Second Order Linear Ordinary Differential Equations ENGI 344 - Second Order Linear ODEs age -01. Second Order Linear Ordinary Differential Equations The general second order linear ordinary differential equation is of the form d y dy x Q x y Rx dx dx Of

More information

Module 3. Analysis of Statically Indeterminate Structures by the Displacement Method

Module 3. Analysis of Statically Indeterminate Structures by the Displacement Method odule 3 Analysis of Statically Indeterminate Structures by the Displacement ethod Lesson 21 The oment- Distribution ethod: rames with Sidesway Instructional Objectives After reading this chapter the student

More information

Tensor Transformations and the Maximum Shear Stress. (Draft 1, 1/28/07)

Tensor Transformations and the Maximum Shear Stress. (Draft 1, 1/28/07) Tensor Transformations and the Maximum Shear Stress (Draft 1, 1/28/07) Introduction The order of a tensor is the number of subscripts it has. For each subscript it is multiplied by a direction cosine array

More information

Chapter 12 Plate Bending Elements. Chapter 12 Plate Bending Elements

Chapter 12 Plate Bending Elements. Chapter 12 Plate Bending Elements CIVL 7/8117 Chapter 12 - Plate Bending Elements 1/34 Chapter 12 Plate Bending Elements Learning Objectives To introduce basic concepts of plate bending. To derive a common plate bending element stiffness

More information

MECE 3321: Mechanics of Solids Chapter 6

MECE 3321: Mechanics of Solids Chapter 6 MECE 3321: Mechanics of Solids Chapter 6 Samantha Ramirez Beams Beams are long straight members that carry loads perpendicular to their longitudinal axis Beams are classified by the way they are supported

More information

Consider an elastic spring as shown in the Fig.2.4. When the spring is slowly

Consider an elastic spring as shown in the Fig.2.4. When the spring is slowly .3 Strain Energy Consider an elastic spring as shown in the Fig..4. When the spring is slowly pulled, it deflects by a small amount u 1. When the load is removed from the spring, it goes back to the original

More information

Part D: Frames and Plates

Part 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 information

Deflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering

Deflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering Deflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering 008/9 Dr. Colin Caprani 1 Contents 1. Introduction... 3 1.1 General... 3 1. Background... 4 1.3 Discontinuity Functions...

More information

CH. 4 BEAMS & COLUMNS

CH. 4 BEAMS & COLUMNS CH. 4 BEAMS & COLUMNS BEAMS Beams Basic theory of bending: internal resisting moment at any point in a beam must equal the bending moments produced by the external loads on the beam Rx = Cc + Tt - If the

More information

CRITERIA FOR SELECTION OF FEM MODELS.

CRITERIA FOR SELECTION OF FEM MODELS. CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.

More information

Deflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering

Deflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering Deflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering 009/10 Dr. Colin Caprani 1 Contents 1. Introduction... 4 1.1 General... 4 1. Background... 5 1.3 Discontinuity Functions...

More information

A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS

A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,

More information

EML4507 Finite Element Analysis and Design EXAM 1

EML4507 Finite Element Analysis and Design EXAM 1 2-17-15 Name (underline last name): EML4507 Finite Element Analysis and Design EXAM 1 In this exam you may not use any materials except a pencil or a pen, an 8.5x11 formula sheet, and a calculator. Whenever

More information

Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.

Introduction 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 information

1859. Forced transverse vibration analysis of a Rayleigh double-beam system with a Pasternak middle layer subjected to compressive axial load

1859. Forced transverse vibration analysis of a Rayleigh double-beam system with a Pasternak middle layer subjected to compressive axial load 1859. Forced transverse vibration analysis of a Rayleigh double-beam system with a Pasternak middle layer subjected to compressive axial load Nader Mohammadi 1, Mehrdad Nasirshoaibi 2 Department of Mechanical

More information

Mechanics of Inflatable Fabric Beams

Mechanics of Inflatable Fabric Beams Copyright c 2008 ICCES ICCES, vol.5, no.2, pp.93-98 Mechanics of Inflatable Fabric Beams C. Wielgosz 1,J.C.Thomas 1,A.LeVan 1 Summary In this paper we present a summary of the behaviour of inflatable fabric

More information

CHAPTER -6- BENDING Part -1-

CHAPTER -6- BENDING Part -1- Ishik University / Sulaimani Civil Engineering Department Mechanics of Materials CE 211 CHAPTER -6- BENDING Part -1-1 CHAPTER -6- Bending Outlines of this chapter: 6.1. Chapter Objectives 6.2. Shear and

More information

Analysis III for D-BAUG, Fall 2017 Lecture 11

Analysis III for D-BAUG, Fall 2017 Lecture 11 Analysis III for D-BAUG, Fall 2017 Lecture 11 Lecturer: Alex Sisto (sisto@math.ethz.ch) 1 Beams (Balken) Beams are basic mechanical systems that appear over and over in civil engineering. We consider a

More information

Chapter 8 Supplement: Deflection in Beams Double Integration Method

Chapter 8 Supplement: Deflection in Beams Double Integration Method Chapter 8 Supplement: Deflection in Beams Double Integration Method 8.5 Beam Deflection Double Integration Method In this supplement, we describe the methods for determining the equation of the deflection

More information

202 Index. failure, 26 field equation, 122 force, 1

202 Index. failure, 26 field equation, 122 force, 1 Index acceleration, 12, 161 admissible function, 155 admissible stress, 32 Airy's stress function, 122, 124 d'alembert's principle, 165, 167, 177 amplitude, 171 analogy, 76 anisotropic material, 20 aperiodic

More information

Stiffness Matrices, Spring and Bar Elements

Stiffness Matrices, Spring and Bar Elements CHAPTER Stiffness Matrices, Spring and Bar Elements. INTRODUCTION The primary characteristics of a finite element are embodied in the element stiffness matrix. For a structural finite element, the stiffness

More information

UNIT IV FLEXIBILTY AND STIFFNESS METHOD

UNIT IV FLEXIBILTY AND STIFFNESS METHOD SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road 517583 QUESTION BANK (DESCRIPTIVE) Subject with Code : SA-II (13A01505) Year & Sem: III-B.Tech & I-Sem Course & Branch: B.Tech

More information

Forced Vibration Analysis of Timoshenko Beam with Discontinuities by Means of Distributions Jiri Sobotka

Forced Vibration Analysis of Timoshenko Beam with Discontinuities by Means of Distributions Jiri Sobotka 21 st International Conference ENGINEERING MECHANICS 2015 Svratka, Czech Republic, May 11 14, 2015 Full Text Paper #018, pp. 45 51 Forced Vibration Analysis of Timoshenko Beam with Discontinuities by Means

More information

Mechanics PhD Preliminary Spring 2017

Mechanics PhD Preliminary Spring 2017 Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n

More information

Final Exam Solution Dynamics :45 12:15. Problem 1 Bateau

Final 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 information

Workshop 8. Lateral Buckling

Workshop 8. Lateral Buckling Workshop 8 Lateral Buckling cross section A transversely loaded member that is bent about its major axis may buckle sideways if its compression flange is not laterally supported. The reason buckling occurs

More information

Mechanics of Materials II. Chapter III. A review of the fundamental formulation of stress, strain, and deflection

Mechanics of Materials II. Chapter III. A review of the fundamental formulation of stress, strain, and deflection Mechanics of Materials II Chapter III A review of the fundamental formulation of stress, strain, and deflection Outline Introduction Assumtions and limitations Axial loading Torsion of circular shafts

More information

2 marks Questions and Answers

2 marks Questions and Answers 1. Define the term strain energy. A: Strain Energy of the elastic body is defined as the internal work done by the external load in deforming or straining the body. 2. Define the terms: Resilience and

More information

Advanced Structural Analysis Prof. Devdas Menon Department of Civil Engineering Indian Institute of Technology, Madras

Advanced Structural Analysis Prof. Devdas Menon Department of Civil Engineering Indian Institute of Technology, Madras Advanced Structural Analysis Prof. Devdas Menon Department of Civil Engineering Indian Institute of Technology, Madras Module - 5.2 Lecture - 28 Matrix Analysis of Beams and Grids (Refer Slide Time: 00:23)

More information

Back Matter Index The McGraw Hill Companies, 2004

Back 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 information

Free vibration analysis of elastically connected multiple-beams with general boundary conditions using improved Fourier series method

Free vibration analysis of elastically connected multiple-beams with general boundary conditions using improved Fourier series method Free vibration analysis of elastically connected multiple-beams with general boundary conditions using improved Fourier series method Jingtao DU*; Deshui XU; Yufei ZHANG; Tiejun YANG; Zhigang LIU College

More information

Advanced Structural Analysis Prof. Devdas Menon Department of Civil Engineering Indian Institute of Technology, Madras

Advanced Structural Analysis Prof. Devdas Menon Department of Civil Engineering Indian Institute of Technology, Madras Advanced Structural Analysis Prof. Devdas Menon Department of Civil Engineering Indian Institute of Technology, Madras Module No. # 5.4 Lecture No. # 30 Matrix Analysis of Beams and Grids (Refer Slide

More information

Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method

Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method 9210-203 Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method You should have the following for this examination one answer book No additional data is attached

More information

[8] Bending and Shear Loading of Beams

[8] Bending and Shear Loading of Beams [8] Bending and Shear Loading of Beams Page 1 of 28 [8] Bending and Shear Loading of Beams [8.1] Bending of Beams (will not be covered in class) [8.2] Bending Strain and Stress [8.3] Shear in Straight

More information

STRESS STRAIN AND DEFORMATION OF SOLIDS, STATES OF STRESS

STRESS STRAIN AND DEFORMATION OF SOLIDS, STATES OF STRESS 1 UNIT I STRESS STRAIN AND DEFORMATION OF SOLIDS, STATES OF STRESS 1. Define: Stress When an external force acts on a body, it undergoes deformation. At the same time the body resists deformation. The

More information

Lecture 15 Strain and stress in beams

Lecture 15 Strain and stress in beams Spring, 2019 ME 323 Mechanics of Materials Lecture 15 Strain and stress in beams Reading assignment: 6.1 6.2 News: Instructor: Prof. Marcial Gonzalez Last modified: 1/6/19 9:42:38 PM Beam theory (@ ME

More information

Discrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method

Discrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method 131 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 21 (2008) Published online (http://hdl.handle.net/10114/1532) Discrete Analysis for Plate Bending Problems by Using

More information

ACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY THE SUPERPOSITION METHOD

ACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY THE SUPERPOSITION METHOD Journal of Sound and Vibration (1999) 219(2), 265 277 Article No. jsvi.1998.1874, available online at http://www.idealibrary.com.on ACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY

More information

NON LINEAR BUCKLING OF COLUMNS Dr. Mereen Hassan Fahmi Technical College of Erbil

NON LINEAR BUCKLING OF COLUMNS Dr. Mereen Hassan Fahmi Technical College of Erbil Abstract: NON LINEAR BUCKLING OF COLUMNS Dr. Mereen Hassan Fahmi Technical College of Erbil The geometric non-linear total potential energy equation is developed and extended to study the behavior of buckling

More information

Development of Beam Equations

Development of Beam Equations CHAPTER 4 Development of Beam Equations Introduction We begin this chapter by developing the stiffness matrix for the bending of a beam element, the most common of all structural elements as evidenced

More information

MODULE C: COMPRESSION MEMBERS

MODULE C: COMPRESSION MEMBERS MODULE C: COMPRESSION MEMBERS This module of CIE 428 covers the following subjects Column theory Column design per AISC Effective length Torsional and flexural-torsional buckling Built-up members READING:

More information

Geometric Stiffness Effects in 2D and 3D Frames

Geometric Stiffness Effects in 2D and 3D Frames Geometric Stiffness Effects in D and 3D Frames CEE 41. Matrix Structural Analsis Department of Civil and Environmental Engineering Duke Universit Henri Gavin Fall, 1 In situations in which deformations

More information

ENGN2340 Final Project: Implementation of a Euler-Bernuolli Beam Element Michael Monn

ENGN2340 Final Project: Implementation of a Euler-Bernuolli Beam Element Michael Monn ENGN234 Final Project: Implementation of a Euler-Bernuolli Beam Element Michael Monn 12/11/13 Problem Definition and Shape Functions Although there exist many analytical solutions to the Euler-Bernuolli

More information

Vibrationdata FEA Matlab GUI Package User Guide Revision A

Vibrationdata FEA Matlab GUI Package User Guide Revision A Vibrationdata FEA Matlab GUI Package User Guide Revision A By Tom Irvine Email: tom@vibrationdata.com March 25, 2014 Introduction Matlab Script: vibrationdata_fea_preprocessor.zip vibrationdata_fea_preprocessor.m

More information

6. Bending CHAPTER OBJECTIVES

6. Bending CHAPTER OBJECTIVES CHAPTER OBJECTIVES Determine stress in members caused by bending Discuss how to establish shear and moment diagrams for a beam or shaft Determine largest shear and moment in a member, and specify where

More information

Investigation of Slopes and Deflections of a Stepped Beam Using a Global Formula for

Investigation of Slopes and Deflections of a Stepped Beam Using a Global Formula for Investigation of Slopes and Deflections of a Stepped Beam Using a Global Formula for An Undergraduate Honors College Thesis in the Departmanet of Mechnical Engineering College of Engineering University

More information

M.S Comprehensive Examination Analysis

M.S Comprehensive Examination Analysis UNIVERSITY OF CALIFORNIA, BERKELEY Spring Semester 2014 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials Name:......................................... M.S Comprehensive

More information

Chapter 12 Elastic Stability of Columns

Chapter 12 Elastic Stability of Columns Chapter 12 Elastic Stability of Columns Axial compressive loads can cause a sudden lateral deflection (Buckling) For columns made of elastic-perfectly plastic materials, P cr Depends primarily on E and

More information

MAE 323: Chapter 6. Structural Models

MAE 323: Chapter 6. Structural Models Common element types for structural analyis: oplane stress/strain, Axisymmetric obeam, truss,spring oplate/shell elements o3d solid ospecial: Usually used for contact or other constraints What you need

More information

3. BEAMS: STRAIN, STRESS, DEFLECTIONS

3. BEAMS: STRAIN, STRESS, DEFLECTIONS 3. BEAMS: STRAIN, STRESS, DEFLECTIONS The beam, or flexural member, is frequently encountered in structures and machines, and its elementary stress analysis constitutes one of the more interesting facets

More information

Prepared by M. GUNASHANKAR AP/MECH DEPARTMENT OF MECHANICAL ENGINEERING

Prepared by M. GUNASHANKAR AP/MECH DEPARTMENT OF MECHANICAL ENGINEERING CHETTINAD COLLEGE OF ENGINEERING AND TECHNOLOGY-KARUR FINITE ELEMENT ANALYSIS 2 MARKS QUESTIONS WITH ANSWER Prepared by M. GUNASHANKAR AP/MECH DEPARTMENT OF MECHANICAL ENGINEERING FINITE ELEMENT ANALYSIS

More information

Advanced Vibrations. Distributed-Parameter Systems: Approximate Methods Lecture 20. By: H. Ahmadian

Advanced 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 information