Triangular Plate Displacement Elements
|
|
- Lily Linette Hopkins
- 5 years ago
- Views:
Transcription
1 Triangular Plate Displacement Elements
2 Chapter : TRIANGULAR PLATE DISPLACEMENT ELEMENTS TABLE OF CONTENTS Page. Introduction Triangular Element Properties Triangle Geometry Triangular Coordinates Triangular CoordinatesProperties Linear Interpolation Coordinate Transformations Partial Derivatives Six Node Quadratic Interpolation Cubic Interpolation Cubic Interpolation Choice : 0 Nodes Cubic Interpolation Choice : 4 Nodes plus 6 Side Slopes Cubic Interpolation Choice : 4 Nodes plus 6 Cartesian Slopes Cubic Interpolation Choice 4: 4 Nodes plus 6 Cartesian Rotations 9.4 The BCIZ Plate Bending Element The Kinematic Constraint The Curvature Displacement Matrix The Element Stiffness Matrix Notes and Bibliography Exercises
3 . TRIANGULAR ELEMENT PROPERTIES.. Introduction This Chapter presents an overview of finite element models for thin plates using the Kirchhoff Plate Bending (KPB) model and the Total Potential Energy (TPE) vaariational principle.the derivation of shape functions for the triangle geometry is covered in the next Chapter. The purpose of this Chapter is to explain the construction of displacement-based triangular Kirchhoff plate bending elements through the development of the necessary interpolation formulas. The resulting element is complete but does not satisfy full normal-slope conformity... Triangular Element Properties We recall the following properties of a straight-sided triangular element, taken from Chapter of IFEM.... Triangle Geometry The geometry of the -node triangle shown in Figure.(a) is specified by the location of its three corner nodes on the {x, y} plane. The nodes are labeled,, while traversing the sides in counterclockwise fashion. The location of the corners is defined by their coordinates: x i, y i, i =,, (.) The area of the triangle is denoted by A and is given by [ A = det x x x = (x y x y ) + (x y x y ) + (x y x y ). (.) y y y It is important to realize that the area given by formula (.) is a signed quantity. It is positive if the corners are numbered in counterclockwise order as shown in Figure.(b). This convention is followed in the sequel. (a) (x,y ) (b) Area A > 0 y (x,y ) x z up, towards you (x,y ) Figure.. Geometry of stright-sided triangular element.
4 Chapter : TRIANGULAR PLATE DISPLACEMENT ELEMENTS... Triangular Coordinates Points of the triangle may also be located in terms of a parametric coordinate system: ζ,ζ,ζ. (.) In the literature these three parameters receive many names. In the sequel the name triangular coordinates will be used to stress its close association with this particular geometry.... Triangular CoordinatesProperties Equations ζ i = constant (.4) represent a set of straight lines parallel to the side opposite to the i th corner. See Figure.. The equation of sides, and are ζ = 0, ζ = 0 and ζ = 0, respectively. The three corners have coordinates (,0,0), (0,,0) and (0,0,). The three midpoints of the sides have coordinates (,, 0), (0,, ) and (, 0, ), the centroid (,, ), and so on. The coordinates are not independent because their sum is unity: ζ + ζ + ζ =. (.) /4 ζ = 0 /4 ζ = / / /4 / /4 /4 /4 ζ = ζ = 0 ζ = ζ = 0 Figure.. Triangular coordinates...4. Linear Interpolation Consider a function w(x, y) that varies linearly over the triangle domain. In terms of Cartesian coordinates it may be expressed as w(x, y) = a 0 + a x + a y, (.6) where a 0, a and a are coefficients to be determined from three conditions. In finite element work such conditions are often the nodal values taken by w at the corners: w,w,w (.7) 4
5 . TRIANGULAR ELEMENT PROPERTIES (a) (b) 6 4 Figure.. Nodal configurations for: (a) linear interpolation of w by three values w i, at corners i =,, ; (b) quadratic interpolation of w by six values w i at corners i =,, and midpoints i = 4,, 6. The expression in triangular coordinates makes direct use of these three values: w(ζ,ζ,ζ ) = w ζ + w ζ + w ζ = [ w w w [ ζ ζ ζ = [ ζ ζ ζ Expression (.8) is called a linear interpolant for w. See Figure.(a).... Coordinate Transformations [ w w w. (.8) Quantities which are closely linked with the element geometry are naturally expressed in triangular coordinates. On the other hand, quantities such as displacements, strains and stresses are often expressed in the Cartesian system x, y. We therefore need transformation equations through which we can pass from one coordinate system to the other. Cartesian and triangular coordinates are linked by the relation [ x y [ = x x x y y y [ ζ ζ. (.9) ζ The first equation says that the sum of the three coordinates is one. The second and third express x and y linearly as homogeneous forms in the triangular coordinates. These simply apply the linear interpolant formula (.8) to the Cartesian coordinates: x = x ζ + x ζ + x ζ and y = y ζ + y ζ + y ζ. Inversion of (.9) yields [ ζ ζ ζ [ [ = x y x y y y x x x y x y y y x x x = A x y x y y y x x y A [ [ A y x A y x x. A y x y (.0) Here x jk = x j x k, y jk = y j y k, A is the triangle area given by (.) and A jk denotes the area subtended by corners j, k and the origin of the x y system. If this origin is taken at the centroid of the triangle, A = A = A = A/.
6 Chapter : TRIANGULAR PLATE DISPLACEMENT ELEMENTS..6. Partial Derivatives From equations (.9) and (.0) we immediately obtain the following relations between partial derivatives: x y = x ζ i, = y i ζ i, (.) i A ζ i x = y jk, A ζ i y = x kj. (.) In (.) j and k denote the cyclic permutations of i. For example, if i =, then j = and k =. The derivatives of a function w(ζ,ζ,ζ ) with respect to x or y follow immediately from (.) and application of the chain rule: w x = ( w y + w y + w ) y A ζ ζ ζ w y = ( w x + w x + w ) (.) x A ζ ζ ζ which in matrix form is..7. Six Node Quadratic Interpolation w w ζ x w = [ y y y w A x x x ζ. (.4) y w ζ Consider next the six-node triangle shown in Figure.(b). This element still has straight sides, and nodes 4, and 6 are located at the midpoint of the sides. In Chapter 6 of the IFEM course it was shown a function w(x, y) that varies quadratically over the element and takes on the node values w i, i =,,, 4,, 6, can be interpolated in terms of the triangular coordinates by the formula N q N q N q w = [ w w w w 4 w w 6 N q. (.) 4 N q N q 6 where the N q i are the quadratic shape functions N q = ζ (ζ ) N q = ζ (ζ ) N q = ζ (ζ ) N q 4 = 4ζ ζ N q = 4ζ ζ N q 6 = 4ζ (.6) ζ The geometry is still defined by (.9) because we are not permitting nodes 4,,6 to be away from the midpoint positions. Similarly, the partial derivative expressions of..6 remain valid. 6
7 . CUBIC INTERPOLATION s (a) 9 8 s 7 (b) s s s s Figure.4. Nodal configurations for two forms of cubic interpolation of w(x, y) over triangle: (a) interpolation by ten values w i at corners i =,,, thirdpoints i = 4,, 6, 7, 8, 9 and centroid i = 0; (b) interpolation by four w i values at corners i =,, and centroid i = 0, plus six side slopes ( w/ s ), ( w/ s ), ( w/ s ), ( w/ s ), ( w/ s ) and ( w/ s )... Cubic Interpolation Cubic interpolants for w are fundamental in the construction of the simplest Kirchhoff plate bending elements using the TPE functional. The complete two-dimensional cubic polynomial has 0 terms. There are several choices for the selection of nodal values that determine that interpolation. Several important ones are examined next.... Cubic Interpolation Choice : 0 Nodes Next consider the ten-node triangle shown in Figure.4(a). This element still has straight sides. Nodes 4 through 9 are placed at the thirdpoints of the sides as indicated. The tenth node is placed at the centroid and is labeled 0. In a Homework of the IFEM course it was shown that a function w(x, y) that varies cubically over the element and takes on the node values w i, i =,,, 4,, 6, 7, 8, 9, 0, can be interpolated in terms of the triangular coordinates by the formula w = [ w w w w 4 w w 6 w 7 w 8 w 9 w 0 7 N c N c N c N c 4 N c N c 6. N0 c. (.7)
8 Chapter : TRIANGULAR PLATE DISPLACEMENT ELEMENTS where the N c i are the 0-node cubic shape functions N c = ζ (ζ )(ζ ) N c = ζ (ζ )(ζ ) N c = ζ (ζ )(ζ ) N4 c = 9 ζ ζ (ζ ) N c = 9 ζ ζ (ζ ) N6 c = 9 ζ ζ (ζ ) N7 c = 9 ζ ζ (ζ ) N8 c = 9 ζ ζ (ζ ) N9 c = 9 ζ ζ (ζ ) N0 c = 7ζ ζ ζ (.8) The geometry is still defined by (.9) because we are not permitting nodes 4 through 9 to be away from the thirdpoint positions as well as forcing 0 to be exactly at the centroid. Similarly, the partial derivative expressions of..6 remain valid.... Cubic Interpolation Choice : 4 Nodes plus 6 Side Slopes Consider now a variant of the cubic interpolation over the triangle as shown in Figure.4(b). This element still has straight sides. In addition to values at the corners,, and the centroid 0, we we specify the side-slope corner derivatives w jk = ( w/ s j ) k. There are six combinations: w = ( w/ s ), w = ( w/ s ), w = ( w/ s ), w = ( w/ s ), w = ( w/ s ) and w = ( w/ s ). The resulting interpolation is w = [ w w w w w w w w w w 0 N c N c N c N c 4 N c N c 6. N0 c. (.9) in which the N c i are the shape functions N c = ζ (ζ + ζ + ζ ) 7ζ ζ ζ N c = L (ζ ζ ζ ζ ζ ) N c = L (ζ ζ ζ ζ ζ ) N4 c = ζ (ζ + ζ + ζ ) 7ζ ζ ζ N c = L (ζ ζ ζ ζ ζ ) N6 c = L (ζ ζ ζ ζ ζ ) N c 7 = ζ (ζ + ζ + ζ ) 7ζ ζ ζ N c 8 = L (ζ ζ ζ ζ ζ ) N c 9 = L (ζ ζ ζ ζ ζ ) N c 0 = 7ζ ζ ζ (.0) and L jk denotes the length of the side joining corners j and k. The function of the corrective terms ζ ζ ζ is to make the first nine functions (.0) vanish at 0. 8
9 The c and c formulas are related by the transformation matrix. CUBIC INTERPOLATION w w w w 4 w w 6 w 7 w 8 w 9 w 0 = L 7 L L 0 4L L 7 L L 0 4L L L 0 0 4L L w w w w w w w w w w 0 (.)... Cubic Interpolation Choice : 4 Nodes plus 6 Cartesian Slopes This choice is similar to the previous one but the corner derivatives are taken with respect to the directions x, y of the Cartesian reference system. The notation used for the slopes is w x = ( w/ x), w y = ( w/ y), w x = ( w/ x), w y = ( w/ y), w x = ( w/ x) and w y = ( w/ y). The resulting interpolation is in which the N c i w = [ w w x w y w w x w y w w x w y w 0 are the shape functions N c N c N c N c 4 N c N c 6. N0 c. (.) N c = ζ (ζ + ζ + ζ ) 7ζ ζ ζ N c = ζ (x ζ x ζ ) + (x x )ζ ζ ζ N c = ζ (y ζ y ζ ) + (y y )ζ ζ ζ N c 4 = ζ (ζ + ζ + ζ ) 7ζ ζ ζ N c = ζ (x ζ x ζ ) + (x x )ζ ζ ζ N c 6 = ζ (y ζ y ζ ) + (y y )ζ ζ ζ N c 7 = ζ (ζ + ζ + ζ ) 7ζ ζ ζ N c 8 = ζ (x ζ x ζ ) + (x x )ζ ζ ζ N9 c = ζ (y ζ y ζ ) + (y y )ζ ζ ζ N0 c = 7ζ ζ ζ (.) in which x ij = x i x j and y ij = y i y j...4. Cubic Interpolation Choice 4: 4 Nodes plus 6 Cartesian Rotations This interpolation is the same as the previous one, with some rearrangements. The rotational degrees of freedom at the corners are introduced using the relations θ x = w/ y and θ y = w/ x.. The 9
10 Chapter : TRIANGULAR PLATE DISPLACEMENT ELEMENTS resulting interpolation is in which the i w = [ w θ x θ y w θ x θ y w θ x θ y w (.4) are the same shape functions of Choice with some sign and index changes: = N c, = N c, = N c,, 8, 4 = N c 4, = N c 6, 6 = N c 7 = N c 7, 8 = N c 9, 9 = N c 0 = N c The BCIZ Plate Bending Element (.) One of the simplest Kirchhoff plate bending elements was presented in 966 by Bazeley, Cheung, Irons and Zienkiewicz [66. This is called the BCIZ element after the authors initials. This element can be derived from the cubic interpolation choice 4 (.4). The technique appears a bit mysterious at first. Basically one can construct a 0 0 plate stiffness matrix. Freedom w 0 at the centroid can be statically condensed out and the resulting 9 9 stiffness used in finite element analysis. Unfortunately the static condensation destroys curvature completeness and the solutions will not generally converge, as emperically shown in Tocher s thesis [770. The idea behing the BCIZ element is that elimination of the centroidal DOF is done in such a way that completeness is maintained, using a kinemtic constraint..4.. The Kinematic Constraint As explained above, we seek a kinematic constraint to eliminate w 0 in terms of the nine connector DOFs: w 0 = [ a w a θ x a θy a w a θ x a θy a w a θ x a θy w θ x θ y w θ x θ y w θ x θ y (.6) Historically the first triangular plate bending element that satisfied completeness and invariance. This element is also important as the motivation for development of the patch test. 0
11 .4 THE BCIZ PLATE BENDING ELEMENT Substituting this into any of the 4 choices of the previous section, w 0 is eliminated. For choice 4: w = [ w θ x θ y w θ x θ y w θ x θ y + a wn 0 + a θ xn 0 + a θyn 0. N9 c4 + a θyn 0. (.7) in which N 0 = 7ζ ζ ζ is the same for all choices. To determine the a s in (.7) we impose the condition that all constant curvature states must be exactly represented. This is a completeness condition. Constant curvature states are associated with quadratic variations of w, which are defined by the 6-node interpolation (.6). The procedure is as follows: (i) Use (.6) to express w 0 in terms of w through w 6 by setting ζ = ζ = ζ = /. (ii) Take the cubic interpolation choice 4 (.4) and express the value of w 4, w and w 6 in terms of the nine connection DOF w through θ y by setting ζ, ζ and ζ to appropriate midpoint coordinates. (iii) Eliminate w 4, w and w 6 from (i) and (ii) to get (.6) and hence (.7). The resulting formula is called the BCIZ interpolation. The associated N s are the BCIZ shape functions. It can be shown that the constraint is q = [/,( x + x + x )/8,( y + y + y )/8, /,(x x + x )/8,(y y + y )/8, /,(x + x x )/8,(y + y y )/8 (.8) The derivation of the entries of (.8) is the subject of a homework Exercise.
12 Chapter : TRIANGULAR PLATE DISPLACEMENT ELEMENTS (a) (b) Figure.. The two -point Gauss integration rules for triangles, which are useful in the computation of the plate element stiffness matrix K. Numbers annotated near sample points are the weights..4.. The Curvature Displacement Matrix Once we have an interpolation formula for w, as in (.7), the curvatures over the plate element can be evaluated by double differentiation of the shape functions with respect to x and y. The resulting relation can be expressed in matrix form as κ xx w/ x B B B B 4 B B 6 B 7 B 8 B 9 κ yy = w/ y = B B B B 4 B B 6 B 7 B 8 B 9 κ xy w/ y y B B B B 4 B B 6 B 7 B 8 B 9 w θ x θ y w θ x θ y w θ x θ y (.9) or κ = Bu (.0) The entries B ij are linear in the triangular coordinates ζ, ζ and ζ. The derivation of the entries of the B matrix is the subject of a homework Exercise..4.. The Element Stiffness Matrix For this displacement element, the stiffness matrix is given by the usual formula K = B T DB d (.) (e) where (e) is the element area, D is the plate rigidity matrix that relates moment to curvatures (see 9..) and B is the curvature-displacement matrix defined in (.0). If D is constant over the element, as frequently assumed, the integrand of K is quadratic in the triangular coordinates. If so the integral can be done exactly using one of the two -point Gauss rule presented in Chapter of IFEM. The two rules are reproduced below for convenience: F(ζ,ζ,ζ ) d = A F(, 6, 6 ) + F( 6,, 6 ) + F( 6, 6, ). (.) (e)
13 .4 THE BCIZ PLATE BENDING ELEMENT F(ζ,ζ,ζ ) d = A F(,, 0) + F(0,, ) + F(, 0, ). (.) (e) The second one is called the midpoint rule because the three sample points are at the triangle midpoints. These rules are depicted in Figures.(a) and (b), respectively. Both are exact up to quadratic polynomials in the triangular coordinates, which is what is needed for the plate stiffness matrix.
14 Chapter : TRIANGULAR PLATE DISPLACEMENT ELEMENTS Notes and Bibliography The foregoing shape functions were first developed and presented in terms of triangular coordinates in the 966 thesis [. 4
15 Exercises Homework Exercises for Chapter Triangular Plate Displacement Elements EXERCISE. [A: Check that the shape functions given for cubic interpolation choice 4 satisfy the conditions of being for the associated DOF, and zero for all others. For example N, which is associated with θ x, must satisfy N (, 0, 0) = 0, N y (, 0, 0) = = θ x, N x (, 0, 0) = 0 = θ y, etc. (E.) Note: It is sufficient to verify N, N, and N because the next six shape functions are obtained by cyclic permutation of subscripts. N 0 is easy to check since its corner value and corner slopes are zero. EXERCISE. [A: Using the chain rule and the results of..6, show that the following formula can be used to compute the plate curvatures from an interpolation formula in w: w x w y w x y y x x y = y x x y y x x y 4A y y x x (x y + x y ) y y x x (x y + x y ) y y x x (x y + x y ) T w ζ w ζ w ζ w ζ ζ w ζ ζ w ζ ζ, (E.) EXERCISE. [A: Derive the expression for the a s in (.7) following the procedure outline in that section, and derive the modified shape functions appearing in (.8). EXERCISE.4 [A: Program the BCIZ element starting from the interpolation formula provided by the previous exercise, the expression (E.) to form the curvature-displacement matrix B and one of the -point Gauss rules given in. to evaluate K (e) = (e) B T DBd. Assume that D = (h /)E in which both the thickness h and the stress-strain matrix E are constant over the element, and supplied through arguments. EXERCISE. [A: Evaluate K (e) of the BCIZ element for a triangle with corners at (0, 0), (, ) and (, ), fabricated with isotropic material of E = 0 and ν = 0, and constant thickness h =. As a check, compute the 9 eigenvalues; three of them should be zero and six positive.
Thin Plate Bending Elements: Overview
Thin Plate Bending Elements: Overview Chapter : THIN PLATE BENDING ELEMENTS: OVERVIEW TABLE OF CONTENTS Page. Introduction..................... 3. Overview of KTP FE Models.............. 3.. Triangles...................
More informationNonconservative Loading: Overview
35 Nonconservative Loading: Overview 35 Chapter 35: NONCONSERVATIVE LOADING: OVERVIEW TABLE OF CONTENTS Page 35. Introduction..................... 35 3 35.2 Sources...................... 35 3 35.3 Three
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 information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationUNIVERSITY OF HAWAII COLLEGE OF ENGINEERING DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING
UNIVERSITY OF HAWAII COLLEGE OF ENGINEERING DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING ACKNOWLEDGMENTS This report consists of the dissertation by Ms. Yan Jane Liu, submitted in partial fulfillment
More informationImplementation of Iso-P Triangular Elements
24 Implementation of Iso-P Triangular Elements 24 1 Chapter 24: IMPLEMENTATION OF ISO-P TRIANGULAR ELEMENTS TABLE OF CONTENTS Page 24.1 Introduction..................... 24 3 24.2 Gauss Quadrature for
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 information3. 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 informationMMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
More informationVariational Formulation of Plane Beam Element
. Variational Formulation of Plane Beam Element Chapter : VARIATIONAL FORMULATION OF PLANE BEAM ELEMENT TABLE OF CONTENTS Page.. Introduction 3.. What is a Beam? 3... Terminology................. 3...
More informationLecture 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 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 informationTRIANGULAR, NINE-DEGREES-OF-FREEDOM, C PLATE BENDING ELEMENT OF QUADRATIC ACCURACY*
QUARTERLY OF APPLIED MATHEMATICS 303 OCTOBER, 1973 TRIANGULAR, NINE-DEGREES-OF-FREEDOM, C PLATE BENDING ELEMENT OF QUADRATIC ACCURACY* By ISAAC FRIED and SHOK KENG YANG (Boston University) Introduction.
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 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 informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
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 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 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 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 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 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 informationZero Energy Modes in One Dimension: An Introduction to Hourglass Modes
Zero Energy Modes in One Dimension: An Introduction to Hourglass Modes David J. Benson March 9, 2003 Reduced integration does a lot of good things for an element: it reduces the computational cost, it
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 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 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 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 informationChapter 10 INTEGRATION METHODS Introduction Unit Coordinate Integration
Chapter 0 INTEGRATION METHODS 0. Introduction The finite element analysis techniques are always based on an integral formulation. At the very minimum it will always be necessary to integrate at least an
More informationCurve Fitting. Objectives
Curve Fitting Objectives Understanding the difference between regression and interpolation. Knowing how to fit curve of discrete with least-squares regression. Knowing how to compute and understand the
More informationMIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING
144 MIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING J. N. Reddy* and Chen-Shyh-Tsay School of Aerospace, Mechanical and Nuclear Engineering, University of Oklahoma, Norman, Oklahoma The paper describes
More informationChapter 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 informationFitting Element Fields
F Fitting Element Fields F Appendix F: FITTING ELEMENT FIELDS TABLE OF CONTENTS Page F. Introduction..................... F F. Least-Square Fitting.................. F F. Line Elements.....................
More informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
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 informationLINEAR AND NONLINEAR SHELL THEORY. Contents
LINEAR AND NONLINEAR SHELL THEORY Contents Strain-displacement relations for nonlinear shell theory Approximate strain-displacement relations: Linear theory Small strain theory Small strains & moderate
More informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
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 informationInterpolation (Shape Functions)
Mètodes Numèrics: A First Course on Finite Elements Interpolation (Shape Functions) Following: Curs d Elements Finits amb Aplicacions (J. Masdemont) http://hdl.handle.net/2099.3/36166 Dept. Matemàtiques
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 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 informationContent. Department of Mathematics University of Oslo
Chapter: 1 MEK4560 The Finite Element Method in Solid Mechanics II (January 25, 2008) (E-post:torgeiru@math.uio.no) Page 1 of 14 Content 1 Introduction to MEK4560 3 1.1 Minimum Potential energy..............................
More informationProfessor Terje Haukaas University of British Columbia, Vancouver The M4 Element. Figure 1: Bilinear Mindlin element.
Professor Terje Hakaas University of British Colmbia, ancover www.inrisk.bc.ca The M Element variety of plate elements exist, some being characterized as Kirchhoff elements, i.e., for thin plates, and
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 informationINTERNATIONAL JOURNAL OF PURE AND APPLIED RESEARCH IN ENGINEERING AND TECHNOLOGY
INTERNATIONAL JOURNAL OF PURE AND APPLIED RESEARCH IN ENGINEERING AND TECHNOLOGY A PATH FOR HORIZING YOUR INNOVATIVE WORK SPECIAL ISSUE FOR INTERNATIONAL LEVEL CONFERENCE "ADVANCES IN SCIENCE, TECHNOLOGY
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 informationExtra Problems for Math 2050 Linear Algebra I
Extra Problems for Math 5 Linear Algebra I Find the vector AB and illustrate with a picture if A = (,) and B = (,4) Find B, given A = (,4) and [ AB = A = (,4) and [ AB = 8 If possible, express x = 7 as
More informationOBJECTIVES & ASSUMPTIONS
OBJECTIVES & ASSUMPTIONS Z Project Objectives: Construct a template for the 5-node pyramid element Conduct higher-order patch tests Verify element formulation Determine the optimal element in bending X
More informationMEC-E8001 FINITE ELEMENT ANALYSIS
MEC-E800 FINIE EEMEN ANAYSIS 07 - WHY FINIE EEMENS AND IS HEORY? Design of machines and structures: Solution to stress or displacement by analytical method is often impossible due to complex geometry,
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 informationPlate bending approximation: thin (Kirchhoff) plates and C 1 continuity requirements
4 Plate bending approximation: thin (Kirchhoff) plates and C 1 continuity requirements 4.1 Introduction The subject of bending of plates and indeed its extension to shells was one of the rst to which the
More informationMeasurement of deformation. Measurement of elastic force. Constitutive law. Finite element method
Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation
More informationComputational Stiffness Method
Computational Stiffness Method Hand calculations are central in the classical stiffness method. In that approach, the stiffness matrix is established column-by-column by setting the degrees of freedom
More informationLevel 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 informationMath 462: Homework 2 Solutions
Math 46: Homework Solutions Paul Hacking /6/. Consider the matrix A = (a) Check that A is orthogonal. (b) Determine whether A is a rotation or a reflection / rotary reflection. [Hint: What is the determinant
More informationContents. Prologue Introduction. Classical Approximation... 19
Contents Prologue........................................................................ 15 1 Introduction. Classical Approximation.................................. 19 1.1 Introduction................................................................
More informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
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 informationSymmetric Bending of Beams
Symmetric Bending of Beams beam is any long structural member on which loads act perpendicular to the longitudinal axis. Learning objectives Understand the theory, its limitations and its applications
More informationENGN 2340 Final Project Report. Optimization of Mechanical Isotropy of Soft Network Material
ENGN 2340 Final Project Report Optimization of Mechanical Isotropy of Soft Network Material Enrui Zhang 12/15/2017 1. Introduction of the Problem This project deals with the stress-strain response of a
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 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 informationNumerical Mathematics & Computing, 7 Ed. 4.1 Interpolation
Numerical Mathematics & Computing, 7 Ed. 4.1 Interpolation Ward Cheney/David Kincaid c UT Austin Engage Learning: Thomson-Brooks/Cole www.engage.com www.ma.utexas.edu/cna/nmc6 November 7, 2011 2011 1 /
More informationIntroduction Linear system Nonlinear equation Interpolation
Interpolation Interpolation is the process of estimating an intermediate value from a set of discrete or tabulated values. Suppose we have the following tabulated values: y y 0 y 1 y 2?? y 3 y 4 y 5 x
More informationTruss Structures: The Direct Stiffness Method
. Truss Structures: The Companies, CHAPTER Truss Structures: The Direct Stiffness Method. INTRODUCTION The simple line elements discussed in Chapter introduced the concepts of nodes, nodal displacements,
More informationTHE UPPER BOUND PROPERTY FOR SOLID MECHANICS OF THE LINEARLY CONFORMING RADIAL POINT INTERPOLATION METHOD (LC-RPIM)
International Journal of Computational Methods Vol. 4, No. 3 (2007) 521 541 c World Scientific Publishing Company THE UPPER BOUND PROPERTY FOR SOLID MECHANICS OF THE LINEARLY CONFORMING RADIAL POINT INTERPOLATION
More informationAssumptions: beam is initially straight, is elastically deformed by the loads, such that the slope and deflection of the elastic curve are
*12.4 SLOPE & DISPLACEMENT BY THE MOMENT-AREA METHOD Assumptions: beam is initially straight, is elastically deformed by the loads, such that the slope and deflection of the elastic curve are very small,
More informationChapter 11 Three-Dimensional Stress Analysis. Chapter 11 Three-Dimensional Stress Analysis
CIVL 7/87 Chapter - /39 Chapter Learning Objectives To introduce concepts of three-dimensional stress and strain. To develop the tetrahedral solid-element stiffness matri. To describe how bod and surface
More informationAdvanced 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 informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
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 informationDirectional Field. Xiao-Ming Fu
Directional Field Xiao-Ming Fu Outlines Introduction Discretization Representation Objectives and Constraints Outlines Introduction Discretization Representation Objectives and Constraints Definition Spatially-varying
More informationStructural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian
Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:
More 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 informationKirchhoff Plates: Field Equations
20 Kirchhoff Plates: Field Equations AFEM Ch 20 Slide 1 Plate Structures A plate is a three dimensional bod characterized b Thinness: one of the plate dimensions, the thickness, is much smaller than the
More informationLectures 9-10: Polynomial and piecewise polynomial interpolation
Lectures 9-1: Polynomial and piecewise polynomial interpolation Let f be a function, which is only known at the nodes x 1, x,, x n, ie, all we know about the function f are its values y j = f(x j ), j
More informationDispersion relation for transverse waves in a linear chain of particles
Dispersion relation for transverse waves in a linear chain of particles V. I. Repchenkov* It is difficult to overestimate the importance that have for the development of science the simplest physical and
More informationLecture 21: Isoparametric Formulation of Plane Elements.
6.6. Rectangular Plane Stress/Strain Element. he CS element uses a triangular shape. he 3 nodes of the CS element allow us to employ linear multivariate approximations to the u and v displacements. he
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 informationEFFECTS OF THERMAL STRESSES AND BOUNDARY CONDITIONS ON THE RESPONSE OF A RECTANGULAR ELASTIC BODY MADE OF FGM
Proceedings of the International Conference on Mechanical Engineering 2007 (ICME2007) 29-31 December 2007, Dhaka, Bangladesh ICME2007-AM-76 EFFECTS OF THERMAL STRESSES AND BOUNDARY CONDITIONS ON THE RESPONSE
More informationMechanical Properties of Materials
Mechanical Properties of Materials Strains Material Model Stresses Learning objectives Understand the qualitative and quantitative description of mechanical properties of materials. Learn the logic of
More informationOutline. Structural Matrices. Giacomo Boffi. Introductory Remarks. Structural Matrices. Evaluation of Structural Matrices
Outline in MDOF Systems Dipartimento di Ingegneria Civile e Ambientale, Politecnico di Milano May 8, 014 Additional Today we will study the properties of structural matrices, that is the operators that
More informationMA 323 Geometric Modelling Course Notes: Day 07 Parabolic Arcs
MA 323 Geometric Modelling Course Notes: Day 07 Parabolic Arcs David L. Finn December 9th, 2004 We now start considering the basic curve elements to be used throughout this course; polynomial curves and
More informationCIVL4332 L1 Introduction to Finite Element Method
CIVL L Introduction to Finite Element Method CIVL L Introduction to Finite Element Method by Joe Gattas, Faris Albermani Introduction The FEM is a numerical technique for solving physical problems such
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 informationStiffness 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 informationMEC-E8001 Finite Element Analysis, Exam (example) 2018
MEC-E8 inite Element Analysis Exam (example) 8. ind the transverse displacement wx ( ) of the structure consisting of one beam element and point forces and. he rotations of the endpoints are assumed to
More informationA new closed-form model for isotropic elastic sphere including new solutions for the free vibrations problem
A new closed-form model for isotropic elastic sphere including new solutions for the free vibrations problem E Hanukah Faculty of Mechanical Engineering, Technion Israel Institute of Technology, Haifa
More informationModule-6: Laminated Composites-II. Learning Unit-1: M6.1. M 6.1 Structural Mechanics of Laminates
Module-6: Laminated Composites-II Learning Unit-1: M6.1 M 6.1 Structural Mechanics of Laminates Classical Lamination Theory: Laminate Stiffness Matrix To this point in the development of classical lamination
More informationChapter 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 informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MATERIALS SCIENCE AND ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139
MASSACHUSETTS NSTTUTE OF TECHNOLOGY DEPARTMENT OF MATERALS SCENCE AND ENGNEERNG CAMBRDGE, MASSACHUSETTS 39 3. MECHANCAL PROPERTES OF MATERALS PROBLEM SET SOLUTONS Reading Ashby, M.F., 98, Tensors: Notes
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 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 informationInterpolation Models CHAPTER 3
CHAPTER 3 CHAPTER OUTLINE 3. Introduction 75 3.2 Polynomial Form of Interpolation Functions 77 3.3 Simple, Comple, and Multiple Elements 78 3.4 Interpolation Polynomial in Terms of Nodal Degrees of Freedom
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 informationFIXED BEAMS IN BENDING
FIXED BEAMS IN BENDING INTRODUCTION Fixed or built-in beams are commonly used in building construction because they possess high rigidity in comparison to simply supported beams. When a simply supported
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 informationAnalytical Mechanics: Elastic Deformation
Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 1 / 60 Agenda Agenda
More informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
More informationGlobal polynomial interpolants suffer from the Runge Phenomenon if the data sites (nodes) are not chosen correctly.
Piecewise polynomial interpolation Global polynomial interpolants suffer from the Runge Phenomenon if the data sites (nodes) are not chosen correctly. In many applications, one does not have the freedom
More informationMathematics FINITE ELEMENT ANALYSIS AS COMPUTATION. What the textbooks don't teach you about finite element analysis. Chapter 3
Mathematics FINITE ELEMENT ANALYSIS AS COMPUTATION What the textbooks don't teach you about finite element analysis Chapter 3 Completeness and continuity: How to choose shape functions? Gangan Prathap
More information