Abstract. keywords: finite elements, beams, Hermitian interpolation, shape functions
|
|
- Arlene Gallagher
- 5 years ago
- Views:
Transcription
1 Generation of shape functions for straight beam elements Charles E. Augarde Originally published in Computers and Structures, 68 (1998) Corrections highlighted. Abstract Straight beam finite elements with greater than two nodes are used for edge stiffening in plane stress analyses and elsewhere. It is often necessary to match the number of nodes on the edge stiffener to the number on a whole plane stress element side. Beam elements employ shape functions which are recognised to be level one Hermitian polynomials. An alternative to the commonly adopted method for determining these shape functions is given in this note, using a formula widely reported in mathematical texts which has hitherto not been applied to this task in the finite element literature. The procedure derives shape functions for beams entirely from the set of Lagrangian interpolating polynomials. Examples are given for the derivation of functions for a three and four-noded beam elements. keywords: finite elements, beams, Hermitian interpolation, shape functions 1
2 Introduction Analysis of structures using the finite element method is well-established. Many formulations exist for complex elements but simple elements remain popular since they are usually well-tested and easy to implement into an analysis program. Two-dimensional plane stress analysis, for thin structures subject to in-plane loading, may employ continuum elements, such as the fifteen-node triangle, having a large number of nodes along a side. Where edge stiffening is required, beam elements can be connected to continua edge nodes. There is then a requirement for formulations of beam elements having more than two nodes. Conventional two-dimensional beam elements have two degrees of freedom at each node: one lateral displacement and one rotation. Unless the structure is loaded entirely laterally, axial stiffness must also be incorporated, by an additional degree of freedom at each node. With this amendment, beam elements are usually referred to as frame elements. The axial effects are uncoupled for straight beams and the determination of suitable shape functions is straightforward [1]. The author has recently contributed to the development of a complex threedimensional finite element model at the Department of Engineering Science, Oxford University to study the damage accruing to surface structures from adjacent tunnelling [2]. Modelling was also undertaken in two-dimensions (to validate the complex model) where tunnel linings were represented with simple beam elements. The procedure outlined in this note was used in the implementation of these beam elements into an existing in-house finite element code (OXFEM). In this note we are concerned with the generation of the shape functions which interpolate the lateral displacements along beam elements having more than two nodes. Bernoulli-Euler beam theory is assumed where transverse shear deformation is zero. While most finite element texts describe the simple twonoded beam [1,4] few explain how more complex elements may be formulated [3]. 2
3 Standard procedure The beam element of Figure 1a is defined by two nodes, a distance l apart along the element local x axis. A common convention, and that adopted here, is to make the nodal coordinates non-dimensional that is x = x l (1) so that the two-noded beam has nodes at x = 0, +1. Using shape functions, lateral displacement w( x) of the two-noded beam element of Figure 1a is w( x) = Nd (2) where N = { N 1 N 2 N 3 N 4 } & d T = { w 1 θ 1 w 2 θ 2 } (3) N j, (j = 1, 4) are the bending shape functions and w i, θ i, (i = 1, 2) are the displacements and rotations at the nodes. Determination of the former, using the method to be found in many finite element texts [4,5], proceeds by first writing w( x) as an n-termed polynomial with unknown coefficients, n being the number of degrees of freedom in the element where, w( x) = Xα (4) X = { 1 x x 2 x 3 } & α T = { α 1 α 2 l α 3 l 2 α 4 l 3 } (5) Taking equation (4) and its first derivative with respect to x, four further equations can be formed, one for each degree of freedom: d = Aα (6) where A = (7)
4 Solving for α in Equation (6) and substituting into Equations (4) and (2) gives N = XA 1 (8) The three-node beam of Figure 1b, of overall length 2l, has nodes at x = 1, 0, +1. The same procedure as above yields N = { N 1 N 2 N 3 N 4 N 5 N 6 }, (9) d T = { w 1 θ 1 w 2 θ 2 w 3 θ 3 }, (10) X = { 1 x x 2 x 3 x 4 x 5 } (11) α T = { α 1 α 2 l α 3 l 2 α 4 l 3 α 5 l 4 α 6 l 5 } (12) and to solve for N requires the inversion of A = (13) Hermitian interpolation The shape functions in Equation (2) are Hermitian polynomials since the displacement w(x) is interpolated from nodal rotations as well as nodal displacements. This contrasts with Lagrangian interpolation, used for continuum elements shape functions and for the axial effects in frame elements. Considering small displacements, the nodal rotations are the first derivatives of the unknown real displacement function at the nodes thus fulfilling the definition of Hermitian interpolation. This property allows an alternative procedure to be used to determine the shape functions to that outlined above The generation of Hermitian (or Hermite) polynomials from Lagrangian interpolation polynomials is described in many mathematical texts [6,7,8]. Despite the clear understanding that bending shape functions are equivalent to Hermitian polynomials the technique described below has not, to the author s knowledge, been linked to the generation of those shape functions. 4
5 One-dimensional interpolation is required for straight beam elements. The single dimension is along the element centreline, defined as the x-axis. Onedimensional Hermitian interpolation for an unknown, w(x) proceeds as [6,7] w(x) = nnod i=1 [ H r 0iw i + H r 1i ( ) dw +... Hri r dx i ( d r ) ] w dx r i (14) where H r ji is a Hermite polynomial of level r, relating to node i and to derivative order j of w. The sum is over the number of nodes, nnod where values of w and its derivatives are available. The level of the polynomial indicates the highest order derivative used in the interpolation. Comparison of Equations (2) and (14) reveals that the bending shape functions are level one Hermitian polynomials as follows: N 1 = H 1 01 N 2 = H 1 11 N 3 = H 1 02 N 4 = H 1 12 (15) Level one Hermitian polynomials are derived from Lagrangian polynomials by the following formulae: H 1 0i = [1 2(x x i )L i(x i )][L i (x)] 2 (16) H 1 1i = (x x i )[L i (x)] 2 (17) where L i (x) is the one-dimensional Langrangian polynomial of degree (nnod 1) calculated at node i, given by L i (x) = nnod j=1,j i x x j x i x j (18) and L i(x) is its first derivative with respect to x. A polynomial of order (nnod 2 1) is required to interpolate over nnod points, each contributing two values. Inspection of Equations (16) and (17) shows that the Hermite polynomials are of the correct order for interpolation. Re-writing Equations (16) and (17) in terms of the non-dimensional coordinate, x H 1 0i = [1 2l( x x i )L i( x i )][L i ( x)] 2 (19) H 1 1i = l( x x i )[L i ( x)] 2 (20) 5
6 The advantage of a derivation based on the Lagrangian polynomials and their first derivatives is that these are already likely to be present in a program code. The former are required for continuum elements and the derivatives are required for isoparametric elements. The use of this procedure also provides a systematic approach to allow simpler coding. Examples Three-node beam A three-node element with six bending degrees of freedom and total length 2l is shown in Figure 1b. This element has nodes at x = ( 1, 0, 1). The axial degrees of freedom are omitted from this element as in the derivations above. Shape functions for these degrees of freedom are the Lagrangian polynomials of order 2. From the preceding section, it is clear that these are also required for the derivation of the bending shape functions. The Lagrangian polynomials are L 1 (x) = (x x 2)(x x 3 ) (x 1 x 2 )(x 1 x 3 ) L 2 (x) = (x x 1)(x x 3 ) (x 2 x 1 )(x 2 x 3 ) L 3 (x) = (x x 1)(x x 2 ) (x 3 x 1 )(x 3 x 2 ) (21) (22) (23) Rewriting in terms of the non-dimensional coordinate x and substituting for values of x at nodes (i.e. ( x 1, x 2, x 3 ) = ( 1, 0, 1) ), gives The derivatives are L 1 ( x) = x ( x 1) 2 (24) L 2 ( x) = 1 x 2 (25) L 3 ( x) = x ( x + 1) (26) 2 L 1( x) = 1 (2 x 1) (27) 2l 6
7 L 2( x) = 2 ( x) (28) l L 3( x) = 1 (2 x + 1) (29) 2l The bending shape functions are equivalent to the following Hermite polynomials: N 1 = H 1 01 N 2 = H 1 11 N 3 = H 1 02 (30) N 4 = H 1 12 N 5 = H 1 03 N 6 = H 1 13 (31) Note that the first four are not the same as those in Equation (15) since each set is based on different order Lagrangian polynomials. From (24) and (27) and with substitution for x 1 we obtain [ ( H01 1 = 1 2l( x + 1) 3 )] [ ] x 2 2l 2 ( x 1) = x x3 1 2 x x5 (32) Similarly, [ ] x 2 H11 1 = l( x + 1) 2 ( x 1) = l [ x 2 x 3 x 4 + x 5] (33) 4 H 1 02 = [1] [ (1 x 2 ) ] 2 = 1 2 x 2 + x 4 (34) H12 1 = l x [ x( x + 1)] 2 = l [ x 2 x 3 + x 5] (35) [ ( )] [ ] 3 x 2 H03 1 = 1 2l( x 1) 2l 2 ( x + 1) = x x3 1 2 x4 3 (36) 4 x5 [ ] x 2 H13 1 = l( x 1) 2 ( x + 1) = l [ x 2 x 3 + x 4 + x 5] (37) 4 The suitability of the above as shape functions can be demonstrated by plotting each polynomial (Figure 2). The polynomials related to displacements have value 1 at the associated node and zero elsewhere. Those related to derivatives have slopes of 1 at the associated node. Four-node beam Figure 1c shows a beam with four nodes having the same overall length as the three-node beam. The non-dimensional coordinates ( x) of the nodes are 7
8 1, 1/3, +1/3 and +1. Eight shape functions are required, derived from the following Lagrangian polynomials and derivatives: L 1 ( x) = 9 [ x ] 9 + x2 x 3 (38) L 2 ( x) = 27 [ ] 1 x2 x x3 (39) L 3 ( x) = 27 [ ] 1 x2 + x x3 (40) L 4 ( x) = 9 [ x ] 9 + x2 + x 3 (41) L 1( x) = 9 [ ] 1 16l x 3 x2 (42) L 2( x) = 27 [ 1 2 ] 16l 3 x + 3 x2 (43) L 3( x) = 27 [1 2 ] 16l 3 x 3 x2 (44) L 4( x) = 9 [ 1 ] 16l x + 3 x2 (45) The first two shape functions, relating to lateral displacement and rotation at node 1, can then be derived as N 1 = 1 [ x 243 x x x x x x (46) 7] 512 N 2 = l [ 1 x 19 x x x 4 99 x 5 81 x x 7] (47) 256 Plotting these two functions (Figures 3a & b) shows that they also satisfy the basic requirement of shape functions as outlined above. Conclusion This work presents an alternative derivation of bending shape functions for simple beam elements, for implementation of many-noded straight beam elements within a finite element analysis code. While the elements described are simple the theory will be of interest to developers of other C 1 continuous elements such as rectangular plates. 8
9 Acknowledgements This work was carried out as part of an EPSRC funded project examining numerical modelling of tunnelling at the Department of Engineering Science, University of Oxford, UK under the guidance of Dr. H.J. Burd. The author would like to acknowledge the help and contribution of Dr. Burd to this work. References 1. Astley, R.J., Finite Elements in Solids and Structures, Chapman & Hall, London, Augarde, C.E., Numerical modelling of tunnelling processes for assessment of damage to buildings, D.Phil thesis, University of Oxford, Dawe, D.J., Matrix and Finite Element Displacement Analysis of Structures, Clarendon Press, Oxford, Cook, R.D., Concepts and Applications of Finite Element Analysis, 2nd edn., John Wiley & Sons, Chichester, Mohr, G.A., Finite Elements for Solids, Fluids and Optimisation, Oxford University Press, Oxford, Jacques, I. and Judd, C., Numerical Analysis, Chapman & Hall, London, Morris, J.L., Computational Methods in Elementary Numerical Analysis, John Wiley & Sons, Chichester, Spanier, J. and Oldham, K.B., An Atlas of Functions, Hemisphere, London,
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 informationMITOCW MITRES2_002S10linear_lec07_300k-mp4
MITOCW MITRES2_002S10linear_lec07_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
More informationBilinear 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 informationCRITERIA 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 informationThe Finite Element Method for Solid and Structural Mechanics
The Finite Element Method for Solid and Structural Mechanics Sixth edition O.C. Zienkiewicz, CBE, FRS UNESCO Professor of Numerical Methods in Engineering International Centre for Numerical Methods in
More informationLecture 2: Finite Elements
Materials Science & Metallurgy Master of Philosophy, Materials Modelling, Course MP7, Finite Element Analysis, H. K. D. H. Bhadeshia Lecture 2: Finite Elements In finite element analysis, functions of
More informationLecture 12: Finite Elements
Materials Science & Metallurgy Part III Course M6 Computation of Phase Diagrams H. K. D. H. Bhadeshia Lecture 2: Finite Elements In finite element analysis, functions of continuous quantities such as temperature
More informationUNCONVENTIONAL 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 informationEsben Byskov. Elementary Continuum. Mechanics for Everyone. With Applications to Structural Mechanics. Springer
Esben Byskov Elementary Continuum Mechanics for Everyone With Applications to Structural Mechanics Springer Contents Preface v Contents ix Introduction What Is Continuum Mechanics? "I Need Continuum Mechanics
More informationQuintic beam closed form matrices (revised 2/21, 2/23/12) 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 informationMAE 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 informationContents 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 informationFinite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems
Finite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems Introduction Till now we dealt only with finite elements having straight edges.
More informationFinite 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 informationLecture 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 informationGeneral 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 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 informationThe CR Formulation: BE Plane Beam
6 The CR Formulation: BE Plane Beam 6 Chapter 6: THE CR FORMUATION: BE PANE BEAM TABE OF CONTENTS Page 6. Introduction..................... 6 4 6.2 CR Beam Kinematics................. 6 4 6.2. Coordinate
More informationFLEXIBILITY 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 informationUniversity of Groningen
University of Groningen Nature-inspired microfluidic propulsion using magnetic actuation Khaderi, S. N.; Baltussen, M. G. H. M.; Anderson, P. D.; Ioan, D.; den Toonder, J.M.J.; Onck, Patrick Published
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationDue Tuesday, September 21 st, 12:00 midnight
Due Tuesday, September 21 st, 12:00 midnight The first problem discusses a plane truss with inclined supports. You will need to modify the MatLab software from homework 1. The next 4 problems consider
More informationGyroscopic matrixes of the straight beams and the discs
Titre : Matrice gyroscopique des poutres droites et des di[...] Date : 29/05/2013 Page : 1/12 Gyroscopic matrixes of the straight beams and the discs Summarized: This document presents the formulation
More informationPLAT DAN CANGKANG (TKS 4219)
PLAT DAN CANGKANG (TKS 4219) SESI I: PLATES Dr.Eng. Achfas Zacoeb Dept. of Civil Engineering Brawijaya University INTRODUCTION Plates are straight, plane, two-dimensional structural components of which
More informationCOMPARATIVE STUDY OF VARIOUS BEAMS UNDER DIFFERENT LOADING CONDITION USING FINITE ELEMENT METHOD
COMPARATIVE STUDY OF VARIOUS BEAMS UNDER DIFFERENT LOADING CONDITION USING FINITE ELEMENT METHOD A THESIS SUBMITTED IN PARTIAL FULLFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF Bachelor of Technology
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 informationFinite 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 informationCIVL 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 informationENGN2340 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 informationSteps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen
Steps in the Finite Element Method Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen General Idea Engineers are interested in evaluating effects such as deformations, stresses,
More informationMethods 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 informationInterpolation 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 information6. 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 informationModelling and Finite Element Analysis of Double Wishbone Suspension
Modelling and Finite Element Analysis of Double Wishbone Suspension Amol Patil, Varsha Patil, Prashant Uhle P.G. Student, Dept. of Mechanical Engineering, S S B T S College of Engineering, Jalgaon, Maharastra,
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 03 Quiz # April 8, 03 Name: SOUTION ID#: PS.: A the
More informationCIV-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 informationLecture 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 informationAEROELASTIC ANALYSIS OF SPHERICAL SHELLS
11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
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 informationSection 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 informationGeometry-dependent MITC method for a 2-node iso-beam element
Structural Engineering and Mechanics, Vol. 9, No. (8) 3-3 Geometry-dependent MITC method for a -node iso-beam element Phill-Seung Lee Samsung Heavy Industries, Seocho, Seoul 37-857, Korea Hyu-Chun Noh
More informationChapter 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 informationEML4507 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 informationInstitute 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 informationFinite element modelling of structural mechanics problems
1 Finite element modelling of structural mechanics problems Kjell Magne Mathisen Department of Structural Engineering Norwegian University of Science and Technology Lecture 10: Geilo Winter School - January,
More informationLecture 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 informationAE3160 Experimental Fluid and Solid Mechanics
AE3160 Experimental Fluid and Solid Mechanics Cantilever Beam Bending Claudio Di Leo 1 Learning Objectives 1. On Structural Mechanics: a) Mechanics of Slender Beams b) Strain Transformation Theory c) Principal
More informationCE 715: Advanced Strength of Materials
CE 715: Advanced Strength of Materials Lecture 1 CE 715 Course Information Instructor: Tasnim Hassan Office: Mann Hall 419 Office Hours: TTh 2:00-4:00 pm Phone: 515-8123 Email: thassan@eos.ncsu.edu 1 Course
More informationVirtual Work and Variational Principles
Virtual Work and Principles Mathematically, the structural analysis problem is a boundary value problem (BVP). Forces, displacements, stresses, and strains are connected and computed within the framework
More informationCHAPTER -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 informationPrinciples of Finite Element for Design Engineers and Analysts. Ayman Shama, Ph.D., P.E., F.ASCE
Principles of Finite Element for Design Engineers and Analysts Ayman Shama, Ph.D., P.E., F.ASCE Outline Principles of Engineering Analysis The development of the finite element method Types of elements
More informationACCURATE MODELLING OF STRAIN DISCONTINUITIES IN BEAMS USING AN XFEM APPROACH
VI International Conference on Adaptive Modeling and Simulation ADMOS 213 J. P. Moitinho de Almeida, P. Díez, C. Tiago and N. Parés (Eds) ACCURATE MODELLING OF STRAIN DISCONTINUITIES IN BEAMS USING AN
More informationBHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I
BHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I 635 8 54. Third Year M E C H A NICAL VI S E M ES TER QUE S T I ON B ANK Subject: ME 6 603 FIN I T E E LE ME N T A N A L YSIS UNI T - I INTRODUCTION
More informationReduction of Random Variables in Structural Reliability Analysis
Reduction of Random Variables in Structural Reliability Analysis S. Adhikari and R. S. Langley Department of Engineering University of Cambridge Trumpington Street Cambridge CB2 1PZ (U.K.) February 21,
More informationA 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 informationBending of Simply Supported Isotropic and Composite Laminate Plates
Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,
More informationChapter 6 2D Elements Plate Elements
Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda
More informationToward a novel approach for damage identification and health monitoring of bridge structures
Toward a novel approach for damage identification and health monitoring of bridge structures Paolo Martino Calvi 1, Paolo Venini 1 1 Department of Structural Mechanics, University of Pavia, Italy E-mail:
More informationPLAXIS. Scientific Manual
PLAXIS Scientific Manual 2016 Build 8122 TABLE OF CONTENTS TABLE OF CONTENTS 1 Introduction 5 2 Deformation theory 7 2.1 Basic equations of continuum deformation 7 2.2 Finite element discretisation 8 2.3
More informationCode No: RT41033 R13 Set No. 1 IV B.Tech I Semester Regular Examinations, November - 2016 FINITE ELEMENT METHODS (Common to Mechanical Engineering, Aeronautical Engineering and Automobile Engineering)
More informationChapter 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 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 information1 Static Plastic Behaviour of Beams
1 Static Plastic Behaviour of Beams 1.1 Introduction Many ductile materials which are used in engineering practice have a considerable reserve capacity beyond the initial yield condition. The uniaxial
More informationDynamic and buckling analysis of FRP portal frames using a locking-free finite element
Fourth International Conference on FRP Composites in Civil Engineering (CICE8) 22-24July 8, Zurich, Switzerland Dynamic and buckling analysis of FRP portal frames using a locking-free finite element F.
More informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More informationNonlinear bending analysis of laminated composite stiffened plates
Nonlinear bending analysis of laminated composite stiffened plates * S.N.Patel 1) 1) Dept. of Civi Engineering, BITS Pilani, Pilani Campus, Pilani-333031, (Raj), India 1) shuvendu@pilani.bits-pilani.ac.in
More informationFREE VIBRATION OF AXIALLY LOADED FUNCTIONALLY GRADED SANDWICH BEAMS USING REFINED SHEAR DEFORMATION THEORY
FREE VIBRATION OF AXIALLY LOADED FUNCTIONALLY GRADED SANDWICH BEAMS USING REFINED SHEAR DEFORMATION THEORY Thuc P. Vo 1, Adelaja Israel Osofero 1, Marco Corradi 1, Fawad Inam 1 1 Faculty of Engineering
More informationCourse No: (1 st version: for graduate students) Course Name: Continuum Mechanics Offered by: Chyanbin Hwu
Course No: (1 st version: for graduate students) Course Name: Continuum Mechanics Offered by: Chyanbin Hwu 2011. 11. 25 Contents: 1. Introduction 1.1 Basic Concepts of Continuum Mechanics 1.2 The Need
More informationCorrection of local-linear elasticity for nonlocal residuals: Application to Euler-Bernoulli beams
Correction of local-linear elasticity for nonlocal residuals: Application to Euler-Bernoulli beams Mohamed Shaat* Engineering and Manufacturing Technologies Department, DACC, New Mexico State University,
More informationInstitute 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 informationCHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES
CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric
More informationAN INTEGRATED KIRCHHOFF PLATE ELEMENT BY GALERKIN METHOD FOR THE ANALYSIS OF PLATES ON ELASTIC FOUNDATION
AN INTEGRATED IRCHHOFF PLATE ELEMENT BY GALERIN METHOD FOR THE ANALYSIS OF PLATES ON ELASTIC FOUNDATION Ragesh.P.P, V.Mustafa 2, T.P.Somasundaran 3 (Research Scholar, National Institute of Technology Calicut,
More informationAn Introduction to the Finite Element Method
An Introduction to the Finite Element Method Third Edition J. N. REDDY Department 01 Mechanical Engineering Texas A&M University College Station, Texas, USA 77843 11 Boston Burr Ridge, IL Dubuque, IA Madison,
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 informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationDynamic Model of a Badminton Stroke
ISEA 28 CONFERENCE Dynamic Model of a Badminton Stroke M. Kwan* and J. Rasmussen Department of Mechanical Engineering, Aalborg University, 922 Aalborg East, Denmark Phone: +45 994 9317 / Fax: +45 9815
More informationFREE VIBRATIONS OF UNIFORM TIMOSHENKO BEAMS ON PASTERNAK FOUNDATION USING COUPLED DISPLACEMENT FIELD METHOD
A R C H I V E O F M E C H A N I C A L E N G I N E E R I N G VOL. LXIV 17 Number 3 DOI: 1.1515/meceng-17- Key words: free vibrations, Coupled Displacement Field method, uniform Timoshenko beam, Pasternak
More information3 Relation between complete and natural degrees of freedom
Stiffness matrix for D tapered beams by ouie. Yaw, PhD, PE, SE Walla Walla University March 9, 9 Introduction This article presents information necessary for the construction of the stiffness matrix of
More informationMechanics 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 informationUNIT 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 informationPost Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method
9210-220 Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method You should have the following for this examination one answer book scientific calculator No
More informationk 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 informationUsing MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,
Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business
More informationThe Effect of Distribution for a Moving Force
Paper Number 66, Proceedings of ACOUSTICS 2011 2-4 November 2011, Gold Coast, Australia The Effect of Distribution for a Moving Force Ahmed M. Reda (1,2), Gareth L. Forbes (2) (1) Atkins, Perth, Australia
More informationFORMULA FOR FORCED VIBRATION ANALYSIS OF STRUCTURES USING STATIC FACTORED RESPONSE AS EQUIVALENT DYNAMIC RESPONSE
FORMULA FOR FORCED VIBRATION ANALYSIS OF STRUCTURES USING STATIC FACTORED RESPONSE AS EQUIVALENT DYNAMIC RESPONSE ABSTRACT By G. C. Ezeokpube, M. Eng. Department of Civil Engineering Anambra State University,
More informationFinite Element Method
Finite Element Method Finite Element Method (ENGC 6321) Syllabus Objectives Understand the basic theory of the FEM Know the behaviour and usage of each type of elements covered in this course one dimensional
More informationTable of Contents. Preface... 13
Table of Contents Preface... 13 Chapter 1. Vibrations of Continuous Elastic Solid Media... 17 1.1. Objective of the chapter... 17 1.2. Equations of motion and boundary conditions of continuous media...
More informationUsing the finite element method of structural analysis, determine displacements at nodes 1 and 2.
Question 1 A pin-jointed plane frame, shown in Figure Q1, is fixed to rigid supports at nodes and 4 to prevent their nodal displacements. The frame is loaded at nodes 1 and by a horizontal and a vertical
More informationCOPYRIGHTED MATERIAL. Index
Index A Admissible function, 163 Amplification factor, 36 Amplitude, 1, 22 Amplitude-modulated carrier, 630 Amplitude ratio, 36 Antinodes, 612 Approximate analytical methods, 647 Assumed modes method,
More information2D Kirchhoff Thin Beam Elements
2D Kirchhoff Thin Beam Elements General Element Name Y,v X,u BM3 2 3 1 Element Group Element Subgroup Element Description Number Of Nodes 3 Freedoms Node Coordinates Geometric Properties Kirchhoff Parabolically
More informationKINEMATICS OF CONTINUA
KINEMATICS OF CONTINUA Introduction Deformation of a continuum Configurations of a continuum Deformation mapping Descriptions of motion Material time derivative Velocity and acceleration Transformation
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationMITOCW MITRES2_002S10nonlinear_lec20_300k-mp4
MITOCW MITRES2_002S10nonlinear_lec20_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high-quality educational resources
More informationInternational Journal of Advanced Engineering Technology E-ISSN
Research Article INTEGRATED FORCE METHOD FOR FIBER REINFORCED COMPOSITE PLATE BENDING PROBLEMS Doiphode G. S., Patodi S. C.* Address for Correspondence Assistant Professor, Applied Mechanics Department,
More informationME 1401 FINITE ELEMENT ANALYSIS UNIT I PART -A. 2. Why polynomial type of interpolation functions is mostly used in FEM?
SHRI ANGALAMMAN COLLEGE OF ENGINEERING AND TECHNOLOGY (An ISO 9001:2008 Certified Institution) SIRUGANOOR, TIRUCHIRAPPALLI 621 105 Department of Mechanical Engineering ME 1401 FINITE ELEMENT ANALYSIS 1.
More informationMATERIAL PROPERTIES. Material Properties Must Be Evaluated By Laboratory or Field Tests 1.1 INTRODUCTION 1.2 ANISOTROPIC MATERIALS
. MARIAL PROPRIS Material Properties Must Be valuated By Laboratory or Field ests. INRODUCION he fundamental equations of structural mechanics can be placed in three categories[]. First, the stress-strain
More information2. 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 informationDynamic Analysis of Coupling Vehicle-Bridge System Using Finite Prism Method
Dynamic Analysis of Coupling Vehicle-Bridge System Using Finite Prism Method A. T. Saeed and Zhongfu Xiang Abstract To investigate the transient responses of bridges under moving vehicles, Finite Prism
More informationA New Method of Analysis of Continuous Skew Girder Bridges
A New Method of Analysis of Continuous Skew Girder Bridges Dr. Salem Awad Ramoda Hadhramout University of Science & Technology Mukalla P.O.Box : 50511 Republic of Yemen ABSTRACT One three girder two span
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 information