Chapter 6 2D Elements Plate Elements
|
|
- Darleen Atkinson
- 5 years ago
- Views:
Transcription
1 Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I
2 Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda s Lecture Etending (discretization) to higher Dimensions Structural Elements Plate Elements Method of Finite Elements I
3 Institute of Structural Engineering Page 3 Toda s Lecture (in more detail) In the previous lectures: The Galerkin method was presented The isoparametric concept was introduced The bar and beam element were presented Some numerical integration methods were presented In toda s lecture: Concepts from previous lectures are combined to formulate elements for the above cases Linear 3D elasticit equations are used as a starting point Plane stress/plane strain and plates are presented as special cases Method of Finite Elements I
4 Institute of Structural Engineering Page 4 3D elasticit problem Problem variables: n t Γ t Displacements: u u = u u z Stresses and strains: Γ u Γ 0 Ω ε ε ε ε z = ε ε εz εz εz ε z Method of Finite Elements I σ σ τ τ z = τ σ τ z τz τ z σ z
5 Institute of Structural Engineering Page 5 3D elasticit problem In FE analsis Voigt notation is tpicall used for strains and stresses: ε ε ε ε z = ε ε εz ε ε ε z z zz ε ε ε ε ε ε z ε z = = 2ε γ 2ε z γ z 2 εz γ z engineering strain σ Method of Finite Elements I σ τ τ z = τ σ τ z τ τ σ z z zz σ σ σ σ z = τ τ z τ z
6 Institute of Structural Engineering Page 6 3D elasticit problem n t Γ t Boundar conditions: Ω σn = 0 at Γ σn = t at Γ u= u at Γ u t 0 Γ u Γ 0 Method of Finite Elements I
7 Institute of Structural Engineering Page 7 3D elasticit equations Equilibrium equations: In vector form σ+ F=0 Component wise σ τ τz F z = 0 τ σ τ z F z = 0 τ τ z z σ z Fz z = 0 Where FF is the applied bod force: F F = F F z Method of Finite Elements I
8 Institute of Structural Engineering Page 8 Strain definition: 3D elasticit equations In vector form s 1 ε= u= u+ u 2 Method of Finite Elements I T ( ) Component wise u 1 u u 1 u uz z 1 u u u 1 u u z ε = z 1 uz u 1 u u z uz z 2 z z Using Voigt notation u u uz z ε = u u + u uz + z u u z + z
9 Institute of Structural Engineering Page 9 3D elasticit equations Constitutive equation, using Voigt notation can be convenientl written as: σ = Eε 1 ν ν ν ν 1 ν ν σ ε σ ν ν 1 ν ε σ 1 2 z E ν ε z = τ (1 )(1 2 ) 2 ν ν + γ 1 2ν τ z γ z 2 τ z γ z 1 2ν Method of Finite Elements I
10 Institute of Structural Engineering Page 10 3D elasticit equations The weak form of the problem can be obtained using the Galerkin method as: T T T ε ( w) Eε( u) dω= w FdΩ+ w t dγ Ω Ω Γ where w w = w w z is the weight function Method of Finite Elements I
11 Institute of Structural Engineering Page 11 Plane stress/strain In several cases of practical interest: The third dimension of the problem is: Either ver small Or ver large but includes no variation in the problem parameters The problem equations can be simplified resulting in a 2D problem Method of Finite Elements I
12 Institute of Structural Engineering Page 12 Plane stress assumptions First let s consider a structure where: The length in one dimension is much smaller than the other two Loads are applied onl within a plane The in-plane stresses, strains and displacements are constant in the third dimension Normal and shear stresses in the third dimension are negligible Such a structure is said to be in a state of plane stress
13 Institute of Structural Engineering Page 13 Plane stress equations Assuming that stresses along the third dimension are zero: Assumptions: Equilibrium equations: σ σ τ σ z = σ = σ = τ = τ z (, ) (, ) (, ) = τ z = 0 σ τ τ + σ F F = 0 = 0 Strain definition Constitutive equation u u ε z = ε z u = ε =, εz = 0 γ εz = 0 u u + ε 2 z σ 1 ν 0 ε E σ = ν 1 0 ε 1 ν σ 1 ν γ How is this derived? Method of Finite Elements I
14 Institute of Structural Engineering Page 14 Eamples of plane stress problems Thin Plate With Central Hole Circular Plate Under Edge Loadings
15 Institute of Structural Engineering Page 15 Plane strain assumptions Net we consider a structure where: The length in one dimension is much larger than the other two Loads are applied onl within a plane Loads are constant in the third dimension Displacements and strains along the third dimension are negligible Such a structure is said to be in a state of plane strain
16 Institute of Structural Engineering Page 16 Plane strain equations Net we assume that strains along the third dimension are zero: Assumptions: Equilibrium equations: u = u(, ) u = u(, ) u z = 0 σ τ τ + σ F F = 0 = 0 Strain definition Constitutive equation u ε εz = 0 u ε = ε =, γ z = 0 γ γ z = 0 u u + σ 1 ν ν 0 ε E σ = ν 1 ν 0 ε (1 ν)(1 2 ν) σ + 1 2ν γ Method of Finite Elements I
17 Institute of Structural Engineering Page 17 Eamples of plane strain P z z Long Clinders Under Uniform Loading Semi-Infinite Regions Under Uniform Loadings
18 Institute of Structural Engineering Page 18 Weak form of the problem The weak form of the problem is similar to the one used for the 3D problem: T T T ε ( w) Eε( u) dω= w FdΩ+ w t dγ Ω Ω Γ where: stresses, strains, displacements and constitutive matrices correspond to the plain stress/strain case integration is carried out over an area instead of a volume Method of Finite Elements I
19 Institute of Structural Engineering Page 19 Discretization Several options are available for discretizing the plane stress/strain problem One of the most common choices is isoparametric Lagrange elements In the following we will briefl review some basic properties of isoparametric elements Some specific, and widel used, elements will be presented for plane stress/strain analsis Method of Finite Elements I
20 Institute of Structural Engineering Page 20 Using the isoparametric concept, geometr can be discretized as: n Isoparametric formulation = h(,) rs; = h(,) rs i i i i i= 1 i= 1 n and displacements as: n u = h(,) r su ; u = h(,) r su i i i i i= 1 i= 1 n where i, i, ui, ui are nodal values of the spatial coordinates and displacement components
21 Institute of Structural Engineering Page 21 Isoparametric quadrilaterals For a 4 noded linear isoparametric quadrilateral (q4), coordinates r and s are defined based on the following transformation: 4(, ) 4 4 3(, ) ( 1,1) 3 (1,1) s r 2(, ) 2 2 1( 1, 1) Method of Finite Elements I 1 ( 1, 1) 2 (1, 1)
22 Institute of Structural Engineering Page 22 Isoparametric quadrilaterals Then the displacement and geometr approimations specialize to: 4 4 = h(,) rs; = h(,) rs i i i i i= 1 i= u = h(,) r su ; u = h(,) r su i i i i i= 1 i= 1 with 1 1 h1 = (1 r)(1 s); h2 = (1 + r)(1 s); h3 = (1 + r)(1 + s); h4 = (1 r)(1 + s); 4 4
23 Institute of Structural Engineering Page 23 Isoparametric quadrilaterals h1 h2 h3 h4
24 Institute of Structural Engineering Page 24 Isoparametric quadrilaterals The shape functions can also be written in matri form as: u1 u 1 u 2 u h1 0 h2 0 h3 0 h4 0 u2 = u 0 h1 0 h2 0 h3 0 h u 4 3 u N u 3 u4 u 4 u h1 0 h2 0 h3 0 h4 0 2 = 0 h1 0 h2 0 h3 0 h 4 3 N u= Nu = N
25 Institute of Structural Engineering Page 25 Shape function derivatives Derivatives with respect to the spatial coordinates can be obtained as: r r r = s s s r J = J = Γ r r 1 Γ= J In addition, the infinitesimal surface element can be transformed as: dω= d d = det( J) dr ds
26 Institute of Structural Engineering Page 26 Shape function derivatives Using the above transformations, strain components can be obtained, for instance: ε u hi(,) rs hi(,) rs hi(,) rs with hi, (,) rs = = Γ 11 + Γ 12 r s 4 = = hi, (,) rs ui i= 1 The above can be written in matri form as: ε u u u u 2 = h1, 0 h2, 0 h3, 0 h4, 0 u3 u u u
27 Institute of Structural Engineering Page 27 Shape function derivatives All strain components can be similarl computed and gathered in a matri using Voigt notation: u1 u 1 u 2 h1, 0 h2, 0 h3, 0 h4, 0 u2 ε = 0 h1, 0 h2, 0 h3, 0 h4, u3 h1, h1, h2, h2, h3, h3, h4, h 4, u 3 B u4 u 4 u ε= Bu
28 Institute of Structural Engineering Page 28 Stiffness matri and load vectors Net the strains can be substituted in the weak form to obtain the stiffness matri as: K T = B EBdΩ Ω In the above, the domain of integration is defined with respect to the sstem and the shape functions with respect to r and s. Therefore a change of variables has to be performed as follows: 1 1 K T T = B EB dω= B EB det( J ) drds Ω 1 1
29 Institute of Structural Engineering Page 29 Stiffness matri and load vectors In a similar manner load vectors can be obtained, for instance due to a bod force: 1 1 T f NFdet( J) drds F = 1 1 For surface tractions, integration has to be carried out along the sides of the quadrilateral, e.g.: 1 T ft = N tdet( Js1) dr where det(jj ssss ) is the Jacobean determinant 1 of the side where the load is applied
30 Institute of Structural Engineering Page 30 Stiffness matri and load vectors B reviewing the derived epressions: 1 1 T K B EBdet( J) drds = T f NFdet( J) drds F = T ft = N tdet( Js1) dr 1 It can be observed that: integration over a general quadrilateral was reduced to integration over a square the epressions to be integrated are now significantl more complicated due to the use of the chain rule analtical integration might not be possible
31 Institute of Structural Engineering Page 31 Stiffness matri and load vectors To overcome this difficult, numerical integration is tpicall emploed, thus reducing integrals to sums: 1 1 T T K = B EBdet( J) drds = ww i jb ( ri, sj) EB( ri, sj)det[ J( ri, sj)] 1 1 i j where wi, ri, si are the weights and coordinates of the Gauss points used Notice that the Jacobian is also evaluated at the different Gauss points since it is not constant in general!
32 Institute of Structural Engineering Page 32 Numerical integration Gauss point weights and coordinates are tpicall precomputed can be found in tables. For instance for a linear quadrilateral: s r In 1D: w w 1 = 1, r =, 3 1 = 1, r = The above can be combined to obtain the 2D coordinates and weights
33 Institute of Structural Engineering Page 33 Other isoparametric elements The procedure presented above can be performed for different kinds of elements just b changing the shape functions and Gauss points used, for instance: 8 noded quadratic quadrilateral (q8) Shape functions 1 h1 = (1 r )(1 s ) 1 r s 4 ( ) Method of Finite Elements I s r 1 h2 = (1 + r )(1 s ) r s 1 h3 = (1 + r )(1 + s ) r+ s 1 h4 = (1 r )(1 + s ) 4 1 r+ s h h h h = 8 = = = 2 ( 1 r )( 1 s) 2 2 ( 1+ r)( 1 s ) 2 2 ( 1 r )( 1+ s) 2 2 ( 1 r)( 1 s ) 2 ( ) ( ) ( )
34 Institute of Structural Engineering Page 34 Other isoparametric elements The procedure presented above can be performed for different kinds of elements just b changing the shape functions and Gauss points used, for instance: 8 noded quadratic quadrilateral (q8) Gauss points In 1D: s r w w w = 5 / 9, r = 0.6 = 8 / 9, r = = 5 / 9, r = 0.6, 3 3 The above can be combined to obtain the 2D coordinates and weights
35 Institute of Structural Engineering Page 35 Other isoparametric elements The procedure presented above can be performed for different kinds of elements just b changing the shape functions and Gauss points used, for instance: Linear or constant strain triangle (CST) Shape functions 3 3 h h 1 2 = r = s h = r s 3 1 s r 2
36 Institute of Structural Engineering Page 36 Other isoparametric elements The procedure presented above can be performed for different kinds of elements just b changing the shape functions and Gauss points used, for instance: Linear or constant strain triangle (CST) Gauss points w = 1, r = 1/3, s = 1/ s r 2
37 Institute of Structural Engineering Page 37 The plate problem We net consider structures where: One dimension (thickness) is sufficientl smaller that the other two dimensions The geometr is planar Loads are applied both in plane and normal to the plane The problem can again be reduced to a 2D problem
38 Institute of Structural Engineering Page 38 The plate problem Depending on the loading, the problem can be reduced to two separate problems: In plane loading Normal loading Plane stress (or membrane) problem Plate bending problem
39 Institute of Structural Engineering Page 39 Plate theories Two alternative theories eist for the bending of plates which are in complete analog to beam theories: Beam theor Bernoulli-Euler Kinematic assumptions Plane sections that are normal to the ais of the beam remain plane and normal to the ais. Valid for Beams with small height to length ratio Plate theor Kirchhoff-Love Kinematic assumptions Straight lines that are normal to the mid-surface remain straight and normal to the mid-surface. Valid for Thin plates
40 Institute of Structural Engineering Page 40 Plate theories Two alternative theories eist for the bending of plates which are in complete analog to beam theories: Beam theor Timoshenko Kinematic assumptions Plane sections that are normal to the ais of the beam remain plane but not necessaril normal to the ais. Valid for Beams with an height to length ratio Plate theor Mindlin-Reissner Kinematic assumptions Straight lines that are normal to the mid-surface remain straight but not necessaril normal to the mid-surface. Valid for Thin and thick plates
41 Institute of Structural Engineering Page 41 Mindlin-Reissner plate theor In the following we will focus on the Mindlin-Reissner theor, since most elements used in practice are based on that theor. The basic equations for this theor can be derived from the 3D elasticit equations b emploing the following assumptions: Deflections don t var along the thickness direction In plane displacements, and as a result strains and stresses, var linearl along the thickness Normal stresses along the thickness direction vanish Method of Finite Elements I
42 Institute of Structural Engineering Page 42 Mindlin-Reissner plate theor The kinematic assumptions can be written as: u = zθ (, ); u = zθ (, ); u = u(, ) z z where θ are the rotations of the normal to the mid-plane along, θ the z and z planes z It should be noted that in this theor rotations are not defined as the derivatives of the deflections Instead the are independent variables B A z A θ B γ z
43 Institute of Structural Engineering Page 43 Mindlin-Reissner plate theor The kinematic assumptions can be written as: u = zθ (, ); u = zθ (, ); u = u(, ) z z Substituting in the strain definition and separating the in-plane from the transverse shear components: ε u θ u θ ε = = z γ u u θ θ + + uz u uz θ γ + z z = = γ z u u z uz θ + z
44 Institute of Structural Engineering Page 44 Mindlin-Reissner plate theor Using plane stress constitutive relations, the stresses can be obtained: θ σ 1 ν 0 E θ σ = z ν ν τ 1 ν 0 0 θ 2 θ + uz θ τ z E = τ z 2(1 + ν ) uz θ In plane stresses var linearl along the thickness Transverse shear stresses are constant along the thickness
45 Institute of Structural Engineering Page 45 Mindlin-Reissner plate theor B separating integration along the thickness and on the area, the weak form can be written as: 1 ν 0 ε γ da dz da dz w p da t/2 t/2 E E z ε ε γ ν ε γ z γ z z z 1 ν + = 2(1 + ν) γ t/2 A t/2 A z A 1 ν γ where: t ε, ε, γ, γ z, γ z w z pz A is the plate thickness are strain components computed using the weight function is the weight function components in the z direction is the normal load applied to the plate is the area of the plate
46 Institute of Structural Engineering Page 46 Mindlin-Reissner plate theor Once integration along the thickness has been carried out, the weak form can be re-written as: T κ Eκ + T b γ Eγ s = A A A with: da da w p da z z κ θ θ θ θ + = E b 1 ν 0 3 Et = ν 12( 1 ν ) 1 ν γ uz θ = uz θ E s Ekt 1 0 = 2(1 + ν ) 0 1 k is a corrective factor to account for variation of shear stresses, usuall k = 5/6
47 Institute of Structural Engineering Page 47 Mindlin-Reissner plate theor Moments and shear forces are obtained from the curvatures and transverse shear stresses through the constitutive equations: M θ M 1 ν 0 3 Et θ = Ebκ, M = ν ( 1 ν ) M 1 ν 0 0 θ 2 θ + uz Q θ Et 1 0 Q= Esγ, = Q 2(1 ν ) uz θ
48 Institute of Structural Engineering Page 48 Isoparametric plate elements Since rotations and deflections are independent, the are also discretized independentl: n n n u = h(,) rs u; θ = h(,) rsθ ; θ = h(,) rsθ z i zi i i i i i= 1 i= 1 i= 1 One of the most commonl used elements for this kind of analsis is the isoparametric quadrilateral previousl presented The shape functions, derivative definitions, jacobian determinant and Gauss points used for this element are identical to the ones presented for the plane stress element
49 Institute of Structural Engineering Page 49 Isoparametric plate elements Element geometr, and degree of freedom definitions: u z4 θ u z3 θ 4 3 θ θ 4 3 z u z1 θ u 1 z2 θ 2 θ θ 1 2
50 Institute of Structural Engineering Page 50 Isoparametric plate elements Curvatures and shear strains can then be obtained as: uz1 θ 1 θ 1 uz2 θ 2 0 h1, 0 0 h2, 0 0 h3, 0 0 h4, 0 θ 2 κ = 0 0 h1, 0 0 h1, 0 0 h1, 0 0 h1, uz3 0 h1, h1, 0 h2, h2, 0 h3, h3, 0 h4, h 4, θ 3 B b θ 3 u z4 θ 4 θ 4 u κ B u = b
51 Institute of Structural Engineering Page 51 Isoparametric plate elements Curvatures and shear strains can then be obtained as: uz1 θ 1 θ 1 uz2 θ 2 h1, h1 0 h2, h2 0 h3, h3 0 h4, h4 0 θ 2 γ = h1, 0 h1 h2, 0 h2 h3, 0 h3 h4, 0 h 4 uz3 B s θ 3 θ 3 u z4 θ 4 θ 4 u γ B u = s
52 Institute of Structural Engineering Page 52 Isoparametric plate elements The element stiffness matri can be obtained as: T T b b b drds s s s K= BEB det( J) + BEBdet( J) drds 1 1 T f = N p det( J) drds 1 1 z Numerical integration can be carried out as for the plane stress/strain elements.
53 Institute of Structural Engineering Page 53 Isoparametric plate elements The element presented: Is simple to formulate and implement Offers good performance for relativel large thickness to length ratios (>1/10) Suffers from locking problems for smaller ratios Method of Finite Elements I
54 Institute of Structural Engineering Page 54 The shear locking problem Let s consider the presented element in a state of pure bending: Element deformation Actual deformation (for a thin plate) Transverse shear strain develops No transverse shear θ θ 1 2 θ θ 1 2 As the thickness decreases, the spurious energ produced b this shear deformation dominates the total energ leading to unrealisticall small deflections!
55 Institute of Structural Engineering Page 55 The shear locking problem Several remedies have been proposed for the shear locking problem such as: Alternative numerical integration schemes Alternative interpolation for the shear strain components Specialized variational formulations B emploing one of the above alternatives, the element can also be used for thin plates! Method of Finite Elements I
Finite elements for plates and shells. Advanced Design for Mechanical System LEC 2008/11/04
Finite elements for plates and shells Sstem LEC 008/11/04 1 Plates and shells A Midplane The mid-surface of a plate is plane, a curved geometr make it a shell The thickness is constant. Ma be isotropic,
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 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 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 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 information1.1 The Equations of Motion
1.1 The Equations of Motion In Book I, balance of forces and moments acting on an component was enforced in order to ensure that the component was in equilibrium. Here, allowance is made for stresses which
More informationTHE GENERAL ELASTICITY PROBLEM IN SOLIDS
Chapter 10 TH GNRAL LASTICITY PROBLM IN SOLIDS In Chapters 3-5 and 8-9, we have developed equilibrium, kinematic and constitutive equations for a general three-dimensional elastic deformable solid bod.
More informationCH.7. PLANE LINEAR ELASTICITY. Multimedia Course on Continuum Mechanics
CH.7. PLANE LINEAR ELASTICITY Multimedia Course on Continuum Mechanics Overview Plane Linear Elasticit Theor Plane Stress Simplifing Hpothesis Strain Field Constitutive Equation Displacement Field The
More informationKINEMATIC RELATIONS IN DEFORMATION OF SOLIDS
Chapter 8 KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS Figure 8.1: 195 196 CHAPTER 8. KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS 8.1 Motivation In Chapter 3, the conservation of linear momentum for a
More informationThe Plane Stress Problem
. 4 The Plane Stress Problem 4 Chapter 4: THE PLANE STRESS PROBLEM 4 TABLE OF CONTENTS Page 4.. INTRODUCTION 4 3 4... Plate in Plane Stress............... 4 3 4... Mathematical Model.............. 4 4
More informationThe Plane Stress Problem
14 The Plane Stress Problem IFEM Ch 14 Slide 1 Plate in Plane Stress Thickness dimension or transverse dimension z Top surface Inplane dimensions: in, plane IFEM Ch 14 Slide 2 Mathematical Idealization
More informationMAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design
MAE4700/5700 Finite Element Analsis for Mechanical and Aerospace Design Cornell Universit, Fall 2009 Nicholas Zabaras Materials Process Design and Control Laborator Sible School of Mechanical and Aerospace
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 informationME 7502 Lecture 2 Effective Properties of Particulate and Unidirectional Composites
ME 75 Lecture Effective Properties of Particulate and Unidirectional Composites Concepts from Elasticit Theor Statistical Homogeneit, Representative Volume Element, Composite Material Effective Stress-
More informationAircraft Structures Structural & Loading Discontinuities
Universit of Liège Aerospace & Mechanical Engineering Aircraft Structures Structural & Loading Discontinuities Ludovic Noels Computational & Multiscale Mechanics of Materials CM3 http://www.ltas-cm3.ulg.ac.be/
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 informationTwo-Dimensional Formulation.
Two-Dimensional Formulation mi@seu.edu.cn Outline Introduction( 概论 ) Two vs. Three-Dimensional Problems Plane Strain( 平面应变 ) Plane Stress( 平面应力 ) Boundar Conditions( 边界条件 ) Correspondence between Plane
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 informationChapter 3: BASIC ELEMENTS. solid mechanics)
Chapter 3: BASIC ELEMENTS Section 3.: Preliminaries (review of solid mechanics) Outline Most structural analsis FE codes are displacement based In this chapter we discuss interpolation methods and elements
More informationSurvey of Wave Types and Characteristics
Seminar: Vibrations and Structure-Borne Sound in Civil Engineering Theor and Applications Surve of Wave Tpes and Characteristics Xiuu Gao April 1 st, 2006 Abstract Mechanical waves are waves which propagate
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 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 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 informationVariational formulation of the B.V.P. in terms of displacements is. This principle is particularly useful in analyzing large deflection
4-8 CATIGLIANO' THEOREM OF LEAT WORK. Principle of stationar potential energ: Variational formulation of the B.V.P. in terms of displacements is achieved b the principle. This principle is particularl
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 information18 Equilibrium equation of a laminated plate (a laminate)
z ν = 1 ν = ν = 3 ν = N 1 ν = N z 0 z 1 z N 1 z N z α 1 < 0 α > 0 ν = 1 ν = ν = 3 Such laminates are often described b an orientation code [α 1 /α /α 3 /α 4 ] For eample [0/-45/90/45/0/0/45/90/-45/0] Short
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 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 informationConsider a slender rod, fixed at one end and stretched, as illustrated in Fig ; the original position of the rod is shown dotted.
4.1 Strain If an object is placed on a table and then the table is moved, each material particle moves in space. The particles undergo a displacement. The particles have moved in space as a rigid bod.
More informationINTERFACE CRACK IN ORTHOTROPIC KIRCHHOFF PLATES
Gépészet Budapest 4-5.Ma. G--Section-o ITERFACE CRACK I ORTHOTROPIC KIRCHHOFF PLATES András Szekrénes Budapest Universit of Technolog and Economics Department of Applied Mechanics Budapest Műegetem rkp.
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 informationREVIEW FOR EXAM II. Dr. Ibrahim A. Assakkaf SPRING 2002
REVIEW FOR EXM II. J. Clark School of Engineering Department of Civil and Environmental Engineering b Dr. Ibrahim. ssakkaf SPRING 00 ENES 0 Mechanics of Materials Department of Civil and Environmental
More informationStability Analysis of Laminated Composite Thin-Walled Beam Structures
Paper 224 Civil-Comp Press, 2012 Proceedings of the Eleventh International Conference on Computational Structures echnolog, B.H.V. opping, (Editor), Civil-Comp Press, Stirlingshire, Scotland Stabilit nalsis
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 informationOutline. Organization. Stresses in Beams
Stresses in Beams B the end of this lesson, ou should be able to: Calculate the maimum stress in a beam undergoing a bending moment 1 Outline Curvature Normal Strain Normal Stress Neutral is Moment of
More information4 Strain true strain engineering strain plane strain strain transformation formulae
4 Strain The concept of strain is introduced in this Chapter. The approimation to the true strain of the engineering strain is discussed. The practical case of two dimensional plane strain is discussed,
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 informationx y plane is the plane in which the stresses act, yy xy xy Figure 3.5.1: non-zero stress components acting in the x y plane
3.5 Plane Stress This section is concerned with a special two-dimensional state of stress called plane stress. It is important for two reasons: () it arises in real components (particularl in thin components
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 informationMAE 323: Chapter 4. Plane Stress and Plane Strain. The Stress Equilibrium Equation
The Stress Equilibrium Equation As we mentioned in Chapter 2, using the Galerkin formulation and a choice of shape functions, we can derive a discretized form of most differential equations. In Structural
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 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 informationThe Plane Stress Problem
. 14 The Plane Stress Problem 14 1 Chapter 14: THE PLANE STRESS PROBLEM 14 TABLE OF CONTENTS Page 14.1. Introduction 14 3 14.1.1. Plate in Plane Stress............... 14 3 14.1.. Mathematical Model..............
More informationNote on Mathematical Development of Plate Theories
Advanced Studies in Theoretical Phsics Vol. 9, 015, no. 1, 47-55 HIKARI Ltd,.m-hikari.com http://d.doi.org/10.1988/astp.015.411150 Note on athematical Development of Plate Theories Patiphan Chantaraichit
More informationBOUNDARY EFFECTS IN STEEL MOMENT CONNECTIONS
BOUNDARY EFFECTS IN STEEL MOMENT CONNECTIONS Koung-Heog LEE 1, Subhash C GOEL 2 And Bozidar STOJADINOVIC 3 SUMMARY Full restrained beam-to-column connections in steel moment resisting frames have been
More informationLecture notes Models of Mechanics
Lecture notes Models of Mechanics Anders Klarbring Division of Mechanics, Linköping University, Sweden Lecture 7: Small deformation theories Klarbring (Mechanics, LiU) Lecture notes Linköping 2012 1 /
More informationTheories of Straight Beams
EVPM3ed02 2016/6/10 7:20 page 71 #25 This is a part of the revised chapter in the new edition of the tetbook Energy Principles and Variational Methods in pplied Mechanics, which will appear in 2017. These
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 informationGEOMETRIC NONLINEAR ANALYSIS
GEOMETRIC NONLINEAR ANALYSIS The approach for solving problems with geometric nonlinearity is presented. The ESAComp solution relies on Elmer open-source computational tool [1] for multiphysics problems.
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 informationRigid and Braced Frames
RH 331 Note Set 12.1 F2014abn Rigid and raced Frames Notation: E = modulus of elasticit or Young s modulus F = force component in the direction F = force component in the direction FD = free bod diagram
More informationCalculus of the Elastic Properties of a Beam Cross-Section
Presented at the COMSOL Conference 2009 Milan Calculus of the Elastic Properties of a Beam Cross-Section Dipartimento di Modellistica per l Ingegneria Università degli Studi della Calabria (Ital) COMSOL
More informationStrain Transformation and Rosette Gage Theory
Strain Transformation and Rosette Gage Theor It is often desired to measure the full state of strain on the surface of a part, that is to measure not onl the two etensional strains, and, but also the 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 informationShear and torsion correction factors of Timoshenko beam model for generic cross sections
Shear and torsion correction factors of Timoshenko beam model for generic cross sections Jouni Freund*, Alp Karakoç Online Publication Date: 15 Oct 2015 URL: http://www.jresm.org/archive/resm2015.19me0827.html
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 informationThe Finite Element Method for the Analysis of Linear Systems
Swiss Federal Institute of Technolog Page The Finite Element Method for the Analsis of Linear Sstems Prof. Dr. Michael Havbro Faber Swiss Federal Institute of Technolog ETH Zurich, Switzerland Swiss Federal
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 informationFree Vibration Analysis of Rectangular Plates Using a New Triangular Shear Flexible Finite Element
International Journal of Emerging Engineering Research and echnolog Volume Issue August PP - ISSN - (Print) & ISSN - (Online) Free Vibration Analsis of Rectangular Plates Using a New riangular Shear Fleible
More information15 INTERLAMINAR STRESSES
15 INTERLAMINAR STRESSES 15-1 OUT-OF-PLANE STRESSES Classical laminate plate theor predicts the stresses in the plane of the lamina,, and τ but does not account for out-of-plane stresses, τ and τ. It assumes
More informationLinear Strain Triangle and other types of 2D elements. By S. Ziaei Rad
Linear Strain Triangle and other tpes o D elements B S. Ziaei Rad Linear Strain Triangle (LST or T6 This element is also called qadratic trianglar element. Qadratic Trianglar Element Linear Strain Triangle
More informationExample 3.7 Consider the undeformed configuration of a solid as shown in Figure 3.60.
162 3. The linear 3-D elasticity mathematical model The 3-D elasticity model is of great importance, since it is our highest order hierarchical model assuming linear elastic behavior. Therefore, it provides
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 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 informationVibrational Power Flow Considerations Arising From Multi-Dimensional Isolators. Abstract
Vibrational Power Flow Considerations Arising From Multi-Dimensional Isolators Rajendra Singh and Seungbo Kim The Ohio State Universit Columbus, OH 4321-117, USA Abstract Much of the vibration isolation
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 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 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 information1 Bending of beams Mindlin theory
1 BENDNG OF BEAMS MNDLN THEORY 1 1 Bending of beams Mindlin theory Cross-section kinematics assumptions Distributed load acts in the xz plane, which is also a plane of symmetry of a body Ω v(x = 0 m Vertical
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 informationGeometric Stiffness Effects in 2D and 3D Frames
Geometric Stiffness Effects in D and 3D Frames CEE 41. Matrix Structural Analsis Department of Civil and Environmental Engineering Duke Universit Henri Gavin Fall, 1 In situations in which deformations
More information1 Lecture 1 Fundamental Concepts. 1.1 Introduction
1.1 Introduction 1 ecture 1 Fundamental Concepts The finite element method is a numerical method for solving problems of engineering and physical science. Useful for problems with complicated geometries,
More informationC:\Users\whit\Desktop\Active\304_2012_ver_2\_Notes\4_Torsion\1_torsion.docx 6
C:\Users\whit\Desktop\Active\304_2012_ver_2\_Notes\4_Torsion\1_torsion.doc 6 p. 1 of Torsion of circular bar The cross-sections rotate without deformation. The deformation that does occur results from
More informationSome Remarks on the Reissner Mindlin Plate Model
Some Remarks on the Reissner Mindlin Plate Model Alexandre L. Madureira LNCC Brazil Talk at the Workshop of Numerical Analysis of PDEs LNCC February 10, 2003 1 Outline The 3D Problem and its Modeling Full
More information1 Differential Equations for Solid Mechanics
1 Differential Eqations for Solid Mechanics Simple problems involving homogeneos stress states have been considered so far, wherein the stress is the same throghot the component nder std. An eception to
More informationBending Analysis of Isotropic Rectangular Plate with All Edges Clamped: Variational Symbolic Solution
Journal of Emerging Trends in Engineering and pplied Sciences (JETES) (5): 86-85 Scholarlink Research Institute Journals, 0 (ISSN: -706) jeteas.scholarlinkresearch.org Journal of Emerging Trends in Engineering
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 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 informationExercise solutions: concepts from chapter 5
1) Stud the oöids depicted in Figure 1a and 1b. a) Assume that the thin sections of Figure 1 lie in a principal plane of the deformation. Measure and record the lengths and orientations of the principal
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 informationAccepted Manuscript. R.C. Batra, J. Xiao S (12) Reference: COST Composite Structures. To appear in:
Accepted Manuscript Finite deformations of curved laminated St. Venant-Kirchhoff beam using layerwise third order shear and normal deformable beam theory (TSNDT) R.C. Batra, J. Xiao PII: S0263-8223(12)00486-2
More informationEMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 4 Pure Bending
EA 3702 echanics & aterials Science (echanics of aterials) Chapter 4 Pure Bending Pure Bending Ch 2 Aial Loading & Parallel Loading: uniform normal stress and shearing stress distribution Ch 3 Torsion:
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 informationStress Functions. First Semester, Academic Year 2012 Department of Mechanical Engineering Chulalongkorn University
Lecture Note 3 Stress Functions First Semester, Academic Year 01 Department of Mechanical Engineering Chulalongkorn Universit 1 Objectives #1 Describe the equation of equilibrium using D/3D stress elements
More informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More informationOn the Missing Modes When Using the Exact Frequency Relationship between Kirchhoff and Mindlin Plates
On the Missing Modes When Using the Eact Frequenc Relationship between Kirchhoff and Mindlin Plates C.W. Lim 1,*, Z.R. Li 1, Y. Xiang, G.W. Wei 3 and C.M. Wang 4 1 Department of Building and Construction,
More informationDevelopment of Truss Equations
CIVL 7/87 Chapter 3 - Truss Equations - Part /53 Chapter 3a Development of Truss Equations Learning Objectives To derive the stiffness matri for a bar element. To illustrate how to solve a bar assemblage
More informationLECTURE 13 Strength of a Bar in Pure Bending
V. DEMENKO MECHNCS OF MTERLS 015 1 LECTURE 13 Strength of a Bar in Pure Bending Bending is a tpe of loading under which bending moments and also shear forces occur at cross sections of a rod. f the bending
More informationHierarchic Isogeometric Large Rotation Shell Elements Including Linearized Transverse Shear Parametrization
Proceedings of the 7th GACM Colloquium on Computational Mechanics for Young Scientists from Academia and Industry October 11-13, 2017 in Stuttgart, Germany Hierarchic Isogeometric Large Rotation Shell
More informationApplications of Gauss-Radau and Gauss-Lobatto Numerical Integrations Over a Four Node Quadrilateral Finite Element
Avaiable online at www.banglaol.info angladesh J. Sci. Ind. Res. (), 77-86, 008 ANGLADESH JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH CSIR E-mail: bsir07gmail.com Abstract Applications of Gauss-Radau
More informationStructures. Carol Livermore Massachusetts Institute of Technology
C. Livermore: 6.777J/.7J Spring 007, Lecture 7 - Structures Carol Livermore Massachusetts Institute of Technolog * With thanks to Steve Senturia, from whose lecture notes some of these materials are adapted.
More informationBending of Shear Deformable Plates Resting on Winkler Foundations According to Trigonometric Plate Theory
J. Appl. Comput. ech., 4(3) (018) 187-01 DOI: 10.055/JAC.017.3057.1148 ISSN: 383-4536 jacm.scu.ac.ir Bending of Shear Deformable Plates Resting on Winkler Foundations According to Trigonometric Plate Theor
More informationChapter 2 Overview of the Simulation Methodology
Chapter Overview of the Simulation Methodolog. Introduction This chapter presents an overview of the simulation methodolog that comprises both the art and science involved in simulating phsical phenomena.
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 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 informationStress and Strain ( , 3.14) MAE 316 Strength of Mechanical Components NC State University Department of Mechanical & Aerospace Engineering
(3.8-3.1, 3.14) MAE 316 Strength of Mechanical Components NC State Universit Department of Mechanical & Aerospace Engineering 1 Introduction MAE 316 is a continuation of MAE 314 (solid mechanics) Review
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 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 information4.5 The framework element stiffness matrix
45 The framework element stiffness matri Consider a 1 degree-of-freedom element that is straight prismatic and symmetric about both principal cross-sectional aes For such a section the shear center coincides
More informationVibration Analysis of Isotropic and Orthotropic Plates with Mixed Boundary Conditions
Tamkang Journal of Science and Engineering, Vol. 6, No. 4, pp. 7-6 (003) 7 Vibration Analsis of Isotropic and Orthotropic Plates with Mied Boundar Conditions Ming-Hung Hsu epartment of Electronic Engineering
More information