Chapter 6 2D Elements Plate Elements

Size: px
Start display at page:

Download "Chapter 6 2D Elements Plate Elements"

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

MMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM

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

Chapter 11 Three-Dimensional Stress Analysis. Chapter 11 Three-Dimensional Stress Analysis

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

Kirchhoff Plates: Field Equations

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

Chapter 5 Structural Elements: The truss & beam elements

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

More information

1.1 The Equations of Motion

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

THE GENERAL ELASTICITY PROBLEM IN SOLIDS

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

CH.7. PLANE LINEAR ELASTICITY. Multimedia Course on Continuum Mechanics

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

KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS

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

The Plane Stress Problem

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

The Plane Stress Problem

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

MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design

MAE4700/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 information

Aircraft Structures Kirchhoff-Love Plates

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

ME 7502 Lecture 2 Effective Properties of Particulate and Unidirectional Composites

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

Aircraft Structures Structural & Loading Discontinuities

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

Lecture 15 Strain and stress in beams

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

More information

Two-Dimensional Formulation.

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

Bilinear Quadrilateral (Q4): CQUAD4 in GENESIS

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

More information

Chapter 3: BASIC ELEMENTS. solid mechanics)

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

Survey of Wave Types and Characteristics

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

UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES

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

More information

Finite element modelling of structural mechanics problems

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

Finite Element Method in Geotechnical Engineering

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

Variational formulation of the B.V.P. in terms of displacements is. This principle is particularly useful in analyzing large deflection

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

The Finite Element Method for Solid and Structural Mechanics

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

18 Equilibrium equation of a laminated plate (a laminate)

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

COPYRIGHTED MATERIAL. Index

COPYRIGHTED MATERIAL. Index Index A Admissible function, 163 Amplification factor, 36 Amplitude, 1, 22 Amplitude-modulated carrier, 630 Amplitude ratio, 36 Antinodes, 612 Approximate analytical methods, 647 Assumed modes method,

More information

MITOCW MITRES2_002S10linear_lec07_300k-mp4

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 information

Consider a slender rod, fixed at one end and stretched, as illustrated in Fig ; the original position of the rod is shown dotted.

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

INTERFACE CRACK IN ORTHOTROPIC KIRCHHOFF PLATES

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

ME FINITE ELEMENT ANALYSIS FORMULAS

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

REVIEW FOR EXAM II. Dr. Ibrahim A. Assakkaf SPRING 2002

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

Stability Analysis of Laminated Composite Thin-Walled Beam Structures

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

A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS

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

More information

Outline. Organization. Stresses in Beams

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

4 Strain true strain engineering strain plane strain strain transformation formulae

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

JEPPIAAR ENGINEERING COLLEGE

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

x 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

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

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

MAE 323: Chapter 4. Plane Stress and Plane Strain. The Stress Equilibrium Equation

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

MAE 323: Chapter 6. Structural Models

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

More information

Chapter 12 Plate Bending Elements. Chapter 12 Plate Bending Elements

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

More information

The Plane Stress Problem

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

Note on Mathematical Development of Plate Theories

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

BOUNDARY EFFECTS IN STEEL MOMENT CONNECTIONS

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

Lecture notes Models of Mechanics

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

Theories of Straight Beams

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

Bending of Simply Supported Isotropic and Composite Laminate Plates

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

GEOMETRIC NONLINEAR ANALYSIS

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

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

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

More information

Rigid and Braced Frames

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

Calculus of the Elastic Properties of a Beam Cross-Section

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

Strain Transformation and Rosette Gage Theory

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

Using MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,

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

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

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

The Finite Element Method for the Analysis of Linear Systems

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

7. Hierarchical modeling examples

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

Free Vibration Analysis of Rectangular Plates Using a New Triangular Shear Flexible Finite Element

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

15 INTERLAMINAR STRESSES

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

Linear Strain Triangle and other types of 2D elements. By S. Ziaei Rad

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

Example 3.7 Consider the undeformed configuration of a solid as shown in Figure 3.60.

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

Table of Contents. Preface... 13

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

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

Vibrational Power Flow Considerations Arising From Multi-Dimensional Isolators. Abstract

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

Lecture 21: Isoparametric Formulation of Plane Elements.

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

CRITERIA FOR SELECTION OF FEM MODELS.

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

More information

EML4507 Finite Element Analysis and Design EXAM 1

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

More information

1 Bending of beams Mindlin theory

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

Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian

Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:

More information

Geometric Stiffness Effects in 2D and 3D Frames

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

More information

1 Lecture 1 Fundamental Concepts. 1.1 Introduction

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

C:\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.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 information

Some Remarks on the Reissner Mindlin Plate Model

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

1 Differential Equations for Solid Mechanics

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

Bending Analysis of Isotropic Rectangular Plate with All Edges Clamped: Variational Symbolic Solution

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

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

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

Exercise solutions: concepts from chapter 5

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

UNIVERSITY OF HAWAII COLLEGE OF ENGINEERING DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING

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

Accepted Manuscript. R.C. Batra, J. Xiao S (12) Reference: COST Composite Structures. To appear in:

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

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 4 Pure Bending

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

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

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

More information

Stress Functions. First Semester, Academic Year 2012 Department of Mechanical Engineering Chulalongkorn University

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

NONLINEAR CONTINUUM FORMULATIONS CONTENTS

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

On the Missing Modes When Using the Exact Frequency Relationship between Kirchhoff and Mindlin Plates

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

Development of Truss Equations

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

LECTURE 13 Strength of a Bar in Pure Bending

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

Hierarchic Isogeometric Large Rotation Shell Elements Including Linearized Transverse Shear Parametrization

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

Applications of Gauss-Radau and Gauss-Lobatto Numerical Integrations Over a Four Node Quadrilateral Finite Element

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

Structures. Carol Livermore Massachusetts Institute of Technology

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

Bending of Shear Deformable Plates Resting on Winkler Foundations According to Trigonometric Plate Theory

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

Chapter 2 Overview of the Simulation Methodology

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

Mechanics of Inflatable Fabric Beams

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

More information

MAE 323: Lecture 1. Review

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

Stress and Strain ( , 3.14) MAE 316 Strength of Mechanical Components NC State University Department of Mechanical & Aerospace Engineering

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

ACCURATE MODELLING OF STRAIN DISCONTINUITIES IN BEAMS USING AN XFEM APPROACH

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

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

4.5 The framework element stiffness matrix

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

Vibration Analysis of Isotropic and Orthotropic Plates with Mixed Boundary Conditions

Vibration 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