MAE 323: Chapter 6. Structural Models
|
|
- Samuel Black
- 5 years ago
- Views:
Transcription
1 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 to know before selecting an element: owhat is the dominant mechanism by which work is dissipated? owhat DoF s and derived quantities does the element support? owhat special inputs (if any are required) owhat parent shape(s) does the element support? These considerations make some elements more suitable for certain problems than others 2011 Alex Grishin MAE 323 Chapter 6 1
2 Beams and Shells θ z1 u z1 1 z 3 θ y1 θ x1 u y1 y X-X x ux1 θ z2 z u z2 2 Cross section X-X The user must input cross sectional properties. θ x2 θ y2 y u y2 u x2 v 1,y v 2 1 θ 1 2 θ 2 x u u 1 2 The traditional Euler-Bernoulli beam element has two nodes, with a third optional orientation node in three dimensions. This is a reduced continuum element. Supported nodal DoF s are: displacement, and slope. Element DoF s are bending stress, strain, The element dissipates work by bending and linear extension exclusively. This means that only bending and compressive/tensile stresses and strains are calculated for derived quantities 2011 Alex Grishin MAE 323 Chapter 6 2
3 Beam and Shells z x y u z1 θ z1 1 u z1 θ z1 4 t θ y1 u y1 u x1 θ x1 θ y1 u y1 u x1 θ x1 2 s r u z θ z 1 1 θ y1 u y1 u x1 θ x1 u z1 θ z1 3 θ y1 u y1 u x1 θ x1 1 t 3 s r 2 Traditional (Mindlin) shell elements come in rectangular (4-node) shapes, as well as triangular (3-node) shapes. These elements are also dominated by bending deformation. Nodal and element DoF s are the same as beams but with an additional spatial dimension Again, the user must supply cross-sectional information 2011 Alex Grishin MAE 323 Chapter 6 3
4 Beam and Shells z,u z x C P x,u x z,u z C P θ y x,u x Kirchoff Assumption: Straight lines normal to the mid-plane before bending remain straight and normal to midplane after bending. Suitable for thin shells z,u z x C P x,u x z,u z γ xz C P x,u x Mindlin Assumption: straight lines normal to midplane before bending remain straight BUT not necessarily normal after bending. This is due to transverse shear. Suitable for thick shells 2011 Alex Grishin MAE 323 Chapter 6 4
5 Beam and Shells Today, some commercial FE software products offer fully parametric beam and shell elements like those shown below. The advantage of these types of elements over traditional small-deflection beams and shells is that they carry transverse shear components (and so are better at modeling thick beams and shells), and they more accurately represent large rotations and small curvatures They also offer a fully integrated solution through the cross-section. This is helpful when modeling composites or assessing the affects of cross section warping. The disadvantage of these element types is that they may experience shear locking and/or localized spurious deflections. Also, care must be taken to reproduce Euler- Bernoulli (predominantly bending) behavior 4 4 t t s r 1 2 r s Alex Grishin MAE 323 Chapter
6 2D Continuum Elements s r s r All the above isoparametric elements (whose shape functions we have already seen) can be used to model plane stress, plane strain, and axisymmetric problems in structural mechanics. The work dissipation mechanism is governed by the elastic internal strain energy relationship (Chapter 4, equation 3). Nodal DoF s are displacement x and y displacement. Element DoF s are stress and strain At this point, the reader should notice a relationship between the number of nodes an element has and the number of shape functions. When we solve the algebraic system equations, the solution vector represents shape function coefficients, which in turn represents solution values at nodal locations. Because of this, an element is said to have nodal degrees of freedom Alex Grishin MAE 323 Chapter 6 6
7 3D Continuum Elements s s t r t r All the above isoparametric elements (whose shape functions are obtained the same way as the 2D elements) may be used to model structural problems over any 3D continuum. Nodal and element DoF s are the same as for plane stress and strain, but with an additional spatial dimension The energy dissipation mechanism is again supplied by the elastic strain energy density relationship (Hooke s Law) 2011 Alex Grishin MAE 323 Chapter 6 7
8 Spring Elements y z 1 u z1 u y1 U y2 ux x1 2 u x2 u z2 Type 1: Longitudinal Spring 1 2 u z1 u z2 u y2 u y1 u x1 u x2 Type 2: Coincident-Node Spring We have already seen spring elements (Chapter 2). They may be defined in two spatial dimensions or three. The longitudinal spring is usually fully defined by two nodal locations and a longitudinal stiffness, k (as in Chapter 2). The second type of spring (the Coincident-Node spring) is fully defined by two stiffness components in 2D (three in 3D) and two coincident nodal locations. This type of spring may also have either translational OR rotational DoF s (the latter are identical to the former except that the stiffness values are defined in terms of force x distance divided by degree or radian) The energy dissipation mechanism is supplied by Hooke s Law Alex Grishin MAE 323 Chapter 6 8
9 Canonical DoF s and Numerical DoF s In Chapter 2, the student (hopefully) got a feel for what nodal degrees of freedom are when we spoke of spring, truss and beam elements. In that chapter, we saw that they are simply the primary variable (the one we want to solve for) at nodal locations. If we were to count the total number of nodal DoF s an element has, we would see a pattern that looks something m*n, where m is the spatial dimension and n is the number of nodes. For example, if we have a two-dimensional spring element, it appears that we should have 2 x 2 = 4 nodal DoF s. However, in the element coordinate system, this reduces to 1 x 2 =2 nodal DoF s (spring elements don t use the transverse, or y-element direction). Furthermore, a two-dimensional beam element appears to have (with extension) 6 DoF s instead of 4! So, it seems that the DoF count per node for an element depends on : - coordinate reference -degree and form of the governing equation(s) 2011 Alex Grishin MAE 323 Chapter 6 9
10 Canonical DoF s and Numerical DoF s In an attempt to make this all a little less bewildering, we ll introduce some general principles. It helps to make a distinction between what we ll call canonical DoF s and numerical DoF s. Numerical DoF s will refer to the total number of DoF s per element that we use in constructing the algebraic system. This will usually determine the matrix system size during solution. Canonical DoF s will refer to the theoretical minimum number of DoF s needed to solve the problem (this will be an important idea if you re ever asked to solve a matrix structural problem manually ). The number of canonical DoF s of an element solving Ndifferential equations may N be found by: c ndof = n pi (1) Where p i is the order of the ith differential equation divided by 2 The number of numerical DoF s will be equal to or larger than this depending on the number of independent spatial directions, m: i= 1 m N n ndof = n pij j= 1 i= 1 (2) 2011 Alex Grishin MAE 323 Chapter 6 10
11 Canonical DoF s and Numerical DoF s The number of canonical degrees of freedom, ndof c is important in physics and relates to the idea of generalized coordinates. This notion helps us to understand what an element can do when it is used to approximate a structure. Example 1: A bar/truss element utilizes a single differential equation: 2 u EA + b 2 x = 0 So, p 1 = 2/2=1 x And it has two nodes, n=2 So, by equation (1), ndof c =1*2=2 It is instructive to compare this number to m*n (the number of spatial dimensions times the number of nodes), because this is the number of rigid body DoF s. It is less because the element can only dissipate energy in the direction of its length 2011 Alex Grishin MAE 323 Chapter 6 11
12 Canonical DoF s and Numerical DoF s Thus, this simple assessment of ndof c tells us something very important about truss/bar elements. Without any further calculations or evidence, we can immediately deduce what will happen in the following situation: y x u x =0 1 2 u y =0 F y A truss element fixed at node 1 with a transverse load at node 2 We already know from the previous example that this element cannot dissipate energy in the y-direction (only along it s length). So, what does this mean in terms of the finite element solution? 2011 Alex Grishin MAE 323 Chapter 6 12
13 Canonical DoF s and Numerical DoF s To reinforce this line of reasoning, we can always check it by actually solving the finite element system! We utilize equation 12b of Chapter 2: F 2 2 c cs c cs 1 1x 2 2 v F cs s cs s 1 1y = 2 2 L c cs c cs u F 2 2x 2 2 cs s cs s v F 2 2 y EA u 2011 Alex Grishin MAE 323 Chapter 6 13
14 Canonical DoF s and Numerical DoF s Using the technique we learned there (slide 12 of Chapter 2) of striking out DoF s that are set to zero: F 2 2 c cs c cs 1 1x 2 2 v F cs s cs s 1 1y = 2 2 L c cs c cs u cs s cs s v F 2 y EA u EA c 2 cs 0 = 2 L cs s Fy But recall that c and s are the direction cosines the element makes with the coordinate axes: 2011 Alex Grishin MAE 323 Chapter 6 14
15 Canonical DoF s and Numerical DoF s So, we have: E A L = Fy This equation has no solution (k is singular)! On the other hand, if we had used a beam element (slide 28 of Chapter 2) and set rotation at node 1 equal to zero, we wouldget a solution. So, if we need a line element that carries a transverse load, use beams (this is the difference between trusses and frames)! This is a major theme of ROM (Reduced Order Modeling - using reduced continuum elements) We ll leave the student with a final thought (and exercise). Count the DoF s of the beam of slide 28 in Chapter 2 and compare to equation (1). Then compare this number to the rigid-body DoF s m*n and give an explanation 2011 Alex Grishin MAE 323 Chapter 6 15
16 So far, we have discussed nodal DoF s, which usually represent the primary variable(s) at nodal locations. However, just as the slope of the transverse deflection of beam elements constitutes yet another DoF, quantities involving other deflection derivatives (often called derivedquantities) also constitute DoF s of their own. The most notable examples of such quantities are stress and strain. However, in continuum elements, numerical integration results in an odd phenomenon. Recall that in the FEM, we are actually solving integral equations and the integrals are performed numerically with quadrature rules. This produces some error. This error gets amplified when we take derivatives*, which we must do to calculate stresses and strains (see slides 15 and 16 of Chapter 4). This error is not great (and diminishes with increased order of shape functions used), but is enough to give stress and strain contours unphysical and spurious patterns *This is usually the case when taking derivatives of approximate quantities 2011 Alex Grishin MAE 323 Chapter 6 16
17 As an example, consider a two-element (2 nd order Serendipity) plane stress model with an applied uniform shear stress at the free end: y s r s r v=-1 x Now, let s look at the internal shear stress. Before we continue, students should be able to predict what the analytical result would be 2011 Alex Grishin MAE 323 Chapter 6 17
18 But look at the shear stress contours that result if one simply plots the stresses as calculated naively from equation (8) of Chapter 5 (with second degree isoparametric shape functions): 0 r 0 N 0 u σ = C s 0 N v r s σ xy 2011 Alex Grishin MAE 323 Chapter 6 18
19 Specifically, let s focus on the shear stress distribution at s =1/ 3 σ ( s = 1/ 3) xy The location s = 1/ 3 happens to coincide with one of the Gauss points of a two-point rule The figure marks the x-locations of the same Gauss rule for both elements. Note that the averagestress at Gauss point locations is closer to the analytical result. As we refine the elements over this domain, these points would converge faster to the analytical result than other points Alex Grishin MAE 323 Chapter 6 19
20 The fact that derived quantity averages are closer to the analytical solution at Gauss points reflects the superconvergentproperties of such sites. Even though taking derivatives of approximate quantities amplifies their error, it turns out that Gauss points of a rule 1 order less than that used to perform the integration have the same error as the displacements!* Most commercial FE code calculates and stores stress and strains at these Gauss point locations OR at element centroids. Such locations may be thought of as providing the element DoF s (and are often referred to as such) And because these sites provide such increased accuracy for derived quantities, the value of these quantities at other locations is usually extrapolated from these sites *R. D. Cook, D. S. Malkus, and M. E. Plesha, Concepts and Applications of Finite Element Analysis, 3rd ed. New York, NY, USA: John Wiley & Sons, Alex Grishin MAE 323 Chapter 6 20
21 Thus, to obtain smoother, more accurate stress and strain contour plots, most modern FE codes extrapolate these quantities according to a scheme like the one shown below for 2 nd order isoparametric elements D 4 4 s q A B p C Here, r and s are the usual isoparametric coordinates, such that Gauss Points 1, thru 4 lie at r, s = ± 1/ 3 To make things a little easier, let a new parametrization, p, q range from ±1 over the Gauss points Thus, the two coordinate systems are related by a scale factor of sqrt(3). In other words: p = q = 3r 3s Alex Grishin MAE 323 Chapter 6 21 r
22 To find element stresses anywhere within the domain, use: σ 4 p = i= 1 N So, if we wanted to calculate σ x (say)at point A (node 1: p=q=-sqrt(3)), we would have: p i Where σ i is any stress component at the four Gauss point locations, and the N i are now bilinear shape functions in p and q given by: σ ( 1 p)( 1 q) ( 1+ p)( 1 q) ( 1+ p)( 1+ q) ( 1 p)( 1+ q) 1 N( p, q) = 4 σ = 1.866σ 0.500σ σ 0.500σ i xa x1 x2 x3 x Alex Grishin MAE 323 Chapter 6 22
23 Similar extrapolation formulas can be used for triangular, as well three dimensional tetrahedral and hexahedral elements Once the extrapolation is performed for each element, coincident stress/strain values at corner nodes are averaged. In commercial systems, such plots are often referred to as averaged, or nodal stresses (unmodified stress values at Gauss points are referred to as element stresses). With schemes such as this, only corner node stresses are calculated and stored. Below is an averaged shear stress plot of the problem on slide 17 Even with just two elements, the stress plot is now much closer to the analytical value. This plot at least now starts to make physical sense 2011 Alex Grishin MAE 323 Chapter 6 23
24 Unaveraged element stresses and strains also provide a good measure of mesh quality. This usually manifests itself as elevated stress or strain levels in poorly shaped elements, as shown below Ω 3 Ω 2 Poorly shaped element Leads to spurious strain result Ω 4 Ω Alex Grishin MAE 323 Chapter 6 24
25 Structural Symmetry We have already encountered two types of symmetry found in structural elastostatic problems: planar and axisymmetry. Both of these types of symmetry exploit the fact that a domain s loads, boundary conditions, and geometry admit a planar section in which all other parallel planar sections yield an identical structural response. Other types of symmetry one can encounter in structural problems involve boundary conditions. These are: Symmetry Boundary Condition Anti-Symmetry Boundary Condition Cyclic Symmetry Boundary Condition 2011 Alex Grishin MAE 323 Chapter 6 25
26 Structural Symmetry osymmetry Boundary Condition A symmetry boundary condition is a surface on which primary solution variables (displacement in structural models) are set to zero normal to the surface. Such a boundary condition usually reflects a theoretical dividing line through a domain both sides of which behave identically. The figure below shows a planar model (one type of symmetry), which is identically loaded at opposite ends. The geometry and loading admit a plane of symmetry which can be used to cut the model in half u x =0 y Half-symmetry model x 2011 Alex Grishin MAE 323 Chapter 6 26
27 Structural Symmetry oanti-symmetry Boundary Condition An anti-symmetry boundary condition is a surface on which primary solution variables (displacement in structural models) are set to zero within the plane of the surface. Such a boundary condition also reflects a hypothetical dividing line through a domain, but this time, it usually represents a half of structure the other half of which deforms in the opposite direction along the symmety plane. The figure below shows a planar model whose loading suggests a plane of symmetry parallel to it y u x =0 Half-anti- symmetry model x 2011 Alex Grishin MAE 323 Chapter 6 27
28 Structural Symmetry ocyclic Symmetry Boundary Condition A cyclic symmetry boundary condition can be imposed on a structure if it has repeated or identical features (or parts) around a central axis of revolution. S 2 θ S 1 y θ S = θ S 1 2 r x 2011 Alex Grishin MAE 323 Chapter 6 28
MITOCW 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 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 informationTheoretical Manual Theoretical background to the Strand7 finite element analysis system
Theoretical Manual Theoretical background to the Strand7 finite element analysis system Edition 1 January 2005 Strand7 Release 2.3 2004-2005 Strand7 Pty Limited All rights reserved Contents Preface Chapter
More 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 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 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 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 informationAircraft Structures Kirchhoff-Love Plates
University of Liège erospace & Mechanical Engineering ircraft Structures Kirchhoff-Love Plates Ludovic Noels Computational & Multiscale Mechanics of Materials CM3 http://www.ltas-cm3.ulg.ac.be/ Chemin
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 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 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 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 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 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 informationMatrix Assembly in FEA
Matrix Assembly in FEA 1 In Chapter 2, we spoke about how the global matrix equations are assembled in the finite element method. We now want to revisit that discussion and add some details. For example,
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 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 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 informationMAE 323: Lecture 1. Review
This review is divided into two parts. The first part is a mini-review of statics and solid mechanics. The second part is a review of matrix/vector fundamentals. The first part is given as an refresher
More informationAxisymmetric Modeling. This tutorial gives an overview of axisymmetric modeling. Learn how to:
Axisymmetric Modeling I-DEAS Tutorials: Simulation Projects This tutorial gives an overview of axisymmetric modeling. Learn how to: sketch on the XZ plane apply boundary conditions mesh axisymmetric elements
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 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 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 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 informationFinite Element Modeling and Analysis. CE 595: Course Part 2 Amit H. Varma
Finite Element Modeling and Analysis CE 595: Course Part 2 Amit H. Varma Discussion of planar elements Constant Strain Triangle (CST) - easiest and simplest finite element Displacement field in terms of
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 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 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 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 informationIV B.Tech. I Semester Supplementary Examinations, February/March FINITE ELEMENT METHODS (Mechanical Engineering) Time: 3 Hours Max Marks: 80
www..com www..com Code No: M0322/R07 Set No. 1 IV B.Tech. I Semester Supplementary Examinations, February/March - 2011 FINITE ELEMENT METHODS (Mechanical Engineering) Time: 3 Hours Max Marks: 80 Answer
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 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 informationDue Monday, September 14 th, 12:00 midnight
Due Monday, September 14 th, 1: midnight This homework is considering the analysis of plane and space (3D) trusses as discussed in class. A list of MatLab programs that were discussed in class is provided
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 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 informationDHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS PART A (2 MARKS)
DHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS UNIT I : FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PART A (2 MARKS) 1. Write the types
More informationME 475 Modal Analysis of a Tapered Beam
ME 475 Modal Analysis of a Tapered Beam Objectives: 1. To find the natural frequencies and mode shapes of a tapered beam using FEA.. To compare the FE solution to analytical solutions of the vibratory
More 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 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 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 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 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 informationSemiloof Curved Thin Shell Elements
Semiloof Curved Thin Shell Elements General Element Name Y,v,θy X,u,θx Z,w,θz Element Group Element Subgroup Element Description Number Of Nodes Freedoms Node Coordinates TSL 1 2 Semiloof 3 QSL8 7 8 1
More informationChapter 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 informationME751 Advanced Computational Multibody Dynamics
ME751 Advanced Computational Multibody Dynamics November 2, 2016 Antonio Recuero University of Wisconsin-Madison Quotes of the Day The methods which I set forth do not require either constructions or geometrical
More informationLecture 4: PRELIMINARY CONCEPTS OF STRUCTURAL ANALYSIS. Introduction
Introduction In this class we will focus on the structural analysis of framed structures. We will learn about the flexibility method first, and then learn how to use the primary analytical tools associated
More informationAdvanced Structural Analysis EGF Section Properties and Bending
Advanced Structural Analysis EGF316 3. Section Properties and Bending 3.1 Loads in beams When we analyse beams, we need to consider various types of loads acting on them, for example, axial forces, shear
More informationPresented By: EAS 6939 Aerospace Structural Composites
A Beam Theory for Laminated Composites and Application to Torsion Problems Dr. BhavaniV. Sankar Presented By: Sameer Luthra EAS 6939 Aerospace Structural Composites 1 Introduction Composite beams have
More informationMAE 323: Lecture 1. Review
This review is divided into two parts. The first part is a mini-review of statics and solid mechanics. The second part is a review of matrix/vector fundamentals. The first part is given as an refresher
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 informationModule 5: Laminate Theory Lecture 17: Laminate Constitutive Relations. The Lecture Contains: Laminate Constitutive Relations
Lecture 17: Laminate Constitutive Relations The Lecture Contains: Laminate Constitutive Relations Classification of Laminates Cross-Ply Laminates Specially Orthotropic Laminates Examples Homework References
More informationAbstract. keywords: finite elements, beams, Hermitian interpolation, shape functions
Generation of shape functions for straight beam elements Charles E. Augarde Originally published in Computers and Structures, 68 (1998) 555-560 Corrections highlighted. Abstract Straight beam finite elements
More informationME751 Advanced Computational Multibody Dynamics
ME751 Advanced Computational Multibody Dynamics October 24, 2016 Antonio Recuero University of Wisconsin-Madison Quote of the Day If a cluttered desk is a sign of a cluttered mind, of what, then, is an
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 informationPRACTICE 2 PROYECTO Y CONSTRUCCIÓN DE PUENTES. 1º Máster Ingeniería de Caminos. E.T.S.I. Caminos, canales y puertos (Ciudad Real) 01/06/2016
PRACTICE 2 PROYECTO Y CONSTRUCCIÓN DE PUENTES 1º Máster Ingeniería de Caminos E.T.S.I. Caminos, canales y puertos (Ciudad Real) 01/06/2016 AUTHOR: CONTENT 1. INTRODUCTION... 3 2. BRIDGE GEOMETRY AND MATERIAL...
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 informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
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 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 information7. Hierarchical modeling examples
7. Hierarchical modeling examples The objective of this chapter is to apply the hierarchical modeling approach discussed in Chapter 1 to three selected problems using the mathematical models studied in
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 informationSTRESS 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 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 informationCOORDINATE TRANSFORMATIONS
COORDINAE RANSFORMAIONS Members of a structural system are typically oriented in differing directions, e.g., Fig. 17.1. In order to perform an analysis, the element stiffness equations need to be expressed
More informationChapter 2: Deflections of Structures
Chapter 2: Deflections of Structures Fig. 4.1. (Fig. 2.1.) ASTU, Dept. of C Eng., Prepared by: Melkamu E. Page 1 (2.1) (4.1) (2.2) Fig.4.2 Fig.2.2 ASTU, Dept. of C Eng., Prepared by: Melkamu E. Page 2
More informationCreating Axisymmetric Models in FEMAP
Creating Axisymmetric Models in FEMAP 1. Introduction NE/Nastran does not support 2-d axisymmetric elements. 3-d axisymmetric models are supported, and can be generated with a few additional steps. The
More informationTable of Contents. Preface...xvii. Part 1. Level
Preface...xvii Part 1. Level 1... 1 Chapter 1. The Basics of Linear Elastic Behavior... 3 1.1. Cohesion forces... 4 1.2. The notion of stress... 6 1.2.1. Definition... 6 1.2.2. Graphical representation...
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 informationCivil Engineering Design (1) Analysis and Design of Slabs 2006/7
Civil Engineering Design (1) Analysis and Design of Slabs 006/7 Dr. Colin Caprani, Chartered Engineer 1 Contents 1. Elastic Methods... 3 1.1 Introduction... 3 1. Grillage Analysis... 4 1.3 Finite Element
More informationIf the number of unknown reaction components are equal to the number of equations, the structure is known as statically determinate.
1 of 6 EQUILIBRIUM OF A RIGID BODY AND ANALYSIS OF ETRUCTURAS II 9.1 reactions in supports and joints of a two-dimensional structure and statically indeterminate reactions: Statically indeterminate structures
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 informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 06
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 06 In the last lecture, we have seen a boundary value problem, using the formal
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 informationHIGHER-ORDER THEORIES
HIGHER-ORDER THEORIES THIRD-ORDER SHEAR DEFORMATION PLATE THEORY LAYERWISE LAMINATE THEORY J.N. Reddy 1 Third-Order Shear Deformation Plate Theory Assumed Displacement Field µ u(x y z t) u 0 (x y t) +
More informationHIGHER-ORDER THEORIES
HIGHER-ORDER THEORIES Third-order Shear Deformation Plate Theory Displacement and strain fields Equations of motion Navier s solution for bending Layerwise Laminate Theory Interlaminar stress and strain
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 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 informationD : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each.
GTE 2016 Q. 1 Q. 9 carry one mark each. D : SOLID MECHNICS Q.1 single degree of freedom vibrating system has mass of 5 kg, stiffness of 500 N/m and damping coefficient of 100 N-s/m. To make the system
More informationChapter 3 Variational Formulation & the Galerkin Method
Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation
More informationApplications in Solid Mechanics
Companies, 4 CHAPTER 9 Applications in Solid 9. INTRODUCTION The bar and beam elements discussed in Chapters 4 are line elements, as only a single coordinate axis is required to define the element reference
More information46th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference April 2005 Austin, Texas
th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics & Materials Conference - April, Austin, Texas AIAA - AIAA - Bi-stable Cylindrical Space Frames H Ye and S Pellegrino University of Cambridge, Cambridge,
More informationMechanics Research Communications
Mechanics Research Communications 85 (27) 5 Contents lists available at ScienceDirect Mechanics Research Communications journal h om epa ge: www.elsevier.com/locate/mechrescom Load s temporal characteristics
More information1.105 Solid Mechanics Laboratory Fall 2003
1.105 Solid Mechanics Laboratory Fall 200 Experiment 7 Elastic Buckling. The objectives of this experiment are To study the failure of a truss structure due to local buckling of a compression member. To
More informationMechanical Design in Optical Engineering
OPTI Buckling Buckling and Stability: As we learned in the previous lectures, structures may fail in a variety of ways, depending on the materials, load and support conditions. We had two primary concerns:
More informationAERSYS KNOWLEDGE UNIT
-7016 1. INTRODUCTION The scope of this document is to provide a clarification and a deeper understanding of the two different ways to move the mid plane of the element out of the nodal plane. Although
More informationTwo Tier projects for students in ME 160 class
ME 160 Introduction to Finite Element Method Spring 2016 Topics for Term Projects by Teams of 2 Students Instructor: Tai Ran Hsu, Professor, Dept. of Mechanical engineering, San Jose State University,
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 informationTriangular Plate Displacement Elements
Triangular Plate Displacement Elements Chapter : TRIANGULAR PLATE DISPLACEMENT ELEMENTS TABLE OF CONTENTS Page. Introduction...................... Triangular Element Properties................ Triangle
More informationHeat Transfer Analysis
Heat Transfer 2011 Alex Grishin MAE 323 Chapter 8: Grishin 1 In engineering applications, heat is generally transferred from one location to another and between bodies. This transfer is driven by differences
More informationELASTICITY AND FRACTURE MECHANICS. Vijay G. Ukadgaonker
THEORY OF ELASTICITY AND FRACTURE MECHANICS y x Vijay G. Ukadgaonker Theory of Elasticity and Fracture Mechanics VIJAY G. UKADGAONKER Former Professor Indian Institute of Technology Bombay Delhi-110092
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 information3D Semiloof Thin Beam Elements
3D Semiloof Thin Beam Elements General Element Name Z,w,θz Y,v,θy Element Group X,u,θx Element Subgroup Element Description Number Of Nodes BSL3, BSL4 y 1 4 x z 2 Semiloof 3 Curved beam elements in 3D
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 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 informationVerification Examples. FEM-Design. version
FEM-Design 6.0 FEM-Design version. 06 FEM-Design 6.0 StruSoft AB Visit the StruSoft website for company and FEM-Design information at www.strusoft.com Copyright 06 by StruSoft, all rights reserved. Trademarks
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 informationSolution Manual A First Course in the Finite Element Method 5th Edition Logan
Solution Manual A First Course in the Finite Element Method 5th Edition Logan Instant download and all chapters Solution Manual A First Course in the Finite Element Method 5th Edition Logan https://testbandata.com/download/solution-manual-first-course-finite-elementmethod-5th-edition-logan/
More informationLecture Slides. Chapter 4. Deflection and Stiffness. The McGraw-Hill Companies 2012
Lecture Slides Chapter 4 Deflection and Stiffness The McGraw-Hill Companies 2012 Chapter Outline Force vs Deflection Elasticity property of a material that enables it to regain its original configuration
More informationMarch 24, Chapter 4. Deflection and Stiffness. Dr. Mohammad Suliman Abuhaiba, PE
Chapter 4 Deflection and Stiffness 1 2 Chapter Outline Spring Rates Tension, Compression, and Torsion Deflection Due to Bending Beam Deflection Methods Beam Deflections by Superposition Strain Energy Castigliano
More information