Plate Theory. Section 11: PLATE BENDING ELEMENTS

Similar documents
Plate Theory. Section 13: PLATE BENDING ELEMENTS

Kirchhoff and Mindlin Plates

Plates on elastic foundation

KINEMATICS OF RIGID BODIES

6.5 Plate Problems in Rectangular Coordinates

99/105 Comparison of OrcaFlex with standard theoretical results

Equations of Motion. Figure 1.1.1: a differential element under the action of surface and body forces

KINEMATICS OF RIGID BODIES

1 Bending of a beam with a rectangular section

The Algebra (al-jabr) of Matrices

Chapter 5 Bending Moments and Shear Force Diagrams for Beams

1 Part II: Numerical Integration

Shear and torsion interaction of hollow core slabs

Module 1. Energy Methods in Structural Analysis

ROTATION IN 3D WORLD RIGID BODY MOTION

Multiple Integrals. Review of Single Integrals. Planar Area. Volume of Solid of Revolution

Counting intersections of spirals on a torus

Job No. Sheet 1 of 8 Rev B. Made by IR Date Aug Checked by FH/NB Date Oct Revised by MEB Date April 2006

Electromagnetics P5-1. 1) Physical quantities in EM could be scalar (charge, current, energy) or vector (EM fields).

Explain shortly the meaning of the following eight words in relation to shells structures.

Introduction to Finite Elements in Engineering. Tirupathi R. Chandrupatla Rowan University Glassboro, New Jersey

Approximate Large Deflection Analysis of Thin Rectangular Plates under Distributed Lateral Line Load

Generalizations of the Basic Functional

Multiple Integrals. Review of Single Integrals. Planar Area. Volume of Solid of Revolution

SUPPLEMENTARY INFORMATION

Columns and Stability

Classical Mechanics. From Molecular to Con/nuum Physics I WS 11/12 Emiliano Ippoli/ October, 2011

Section - 2 MORE PROPERTIES

Partial Derivatives. Limits. For a single variable function f (x), the limit lim

First midterm topics Second midterm topics End of quarter topics. Math 3B Review. Steve. 18 March 2009

Differentiability. Alex Nita

13: Diffusion in 2 Energy Groups

Heat flux and total heat

The margin is too narrow to contain a truly remarkable proof.

Summary: Method of Separation of Variables

Vorticity. curvature: shear: fluid elements moving in a straight line but at different speeds. t 1 t 2. ATM60, Shu-Hua Chen

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

and that at t = 0 the object is at position 5. Find the position of the object at t = 2.

INTERFACE DESIGN OF CORD-RUBBER COMPOSITES

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method. Version 2 CE IIT, Kharagpur

ME 141. Lecture 10: Kinetics of particles: Newton s 2 nd Law

A little harder example. A block sits at rest on a flat surface. The block is held down by its weight. What is the interaction pair for the weight?

Simulation of Eclipsing Binary Star Systems. Abstract

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

Reference. Vector Analysis Chapter 2

Chapter 1 VECTOR ALGEBRA

MASKING OF FERROMAGNETIC ELLIPTICAL SHELL IN TRANSVERSE MAGNETIC FIELD

8 factors of x. For our second example, let s raise a power to a power:

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

13.4 Work done by Constant Forces

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.

CAPACITORS AND DIELECTRICS

9-1 (a) A weak electrolyte only partially ionizes when dissolved in water. NaHCO 3 is an

AQA Further Pure 2. Hyperbolic Functions. Section 2: The inverse hyperbolic functions

The Wave Equation I. MA 436 Kurt Bryan

Plane curvilinear motion is the motion of a particle along a curved path which lies in a single plane.

Some Methods in the Calculus of Variations

M344 - ADVANCED ENGINEERING MATHEMATICS

Analysis for Transverse Sensitivity of the Microaccelerometer

Lecture 13 - Linking E, ϕ, and ρ

Recitation 3: More Applications of the Derivative

ES.182A Topic 30 Notes Jeremy Orloff

Eigen Values and Eigen Vectors of a given matrix

Math 113 Exam 1-Review

Plane curvilinear motion is the motion of a particle along a curved path which lies in a single plane.

LINEAR ALGEBRA APPLIED

E S dition event Vector Mechanics for Engineers: Dynamics h Due, next Wednesday, 07/19/2006! 1-30

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

Improper Integrals, and Differential Equations

T 1 T 2 T 3 T 4 They may be illustrated by triangular patterns of numbers (hence their name) as shown:

Magnetic forces on a moving charge. EE Lecture 26. Lorentz Force Law and forces on currents. Laws of magnetostatics

Section 6.1 INTRO to LAPLACE TRANSFORMS

CE 160 Lab 2 Notes: Shear and Moment Diagrams for Beams

MECHANICS OF MATERIALS

Solutions to Physics: Principles with Applications, 5/E, Giancoli Chapter 16 CHAPTER 16

The contact stress problem for a piecewisely defined punch indenting an elastic half space. Jacques Woirgard

Solution Manual. for. Fracture Mechanics. C.T. Sun and Z.-H. Jin

Problems for HW X. C. Gwinn. November 30, 2009

PHYSICS 211 MIDTERM I 21 April 2004

Week 10: Line Integrals

Bending and Free Vibration Analysis of Isotropic Plate Using Refined Plate Theory

Torsional buckling of thin-walled cylinders under circumferentially varying thermal loads

Math 1B, lecture 4: Error bounds for numerical methods

Consequently, the temperature must be the same at each point in the cross section at x. Let:

PDE Notes. Paul Carnig. January ODE s vs PDE s 1

21.6 Green Functions for First Order Equations

The Complete Part Design Handbook

Introduction to the Calculus of Variations

Math 8 Winter 2015 Applications of Integration

A finite thin circular beam element for out-of-plane vibration analysis of curved beams

Interactive Simulation of Elasto-Plastic Materials using the Finite Element Method

16 Newton s Laws #3: Components, Friction, Ramps, Pulleys, and Strings

l 2 p2 n 4n 2, the total surface area of the

Things to Memorize: A Partial List. January 27, 2017

20 MATHEMATICS POLYNOMIALS

Polynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

Calculus - Activity 1 Rate of change of a function at a point.

Maths in Motion. Theo de Haan. Order now: 29,95 euro

Transcription:

Section : PLATE BENDING ELEMENTS Plte Theor A plte is structurl element hose mid surfce lies in flt plne. The dimension in the direction norml to the plne is referred to s the thickness of the plte. A plte is chrcteried b smll thickness reltive to the other dimensions of the plte. Pltes m be clssified b the reltive thickness of the structurl element. We cn spek of:. Reltivel thick pltes ith smll displcements. Reltivel thin pltes ith smll displcements. Ver thin pltes ith lrge deflections. Etremel thin pltes (membrnes) ith lrge or smll deflections If plte is curved s opposed to flt, e refer to these structurl elements s shells. We strt the discussion ith clssicl thin plte theor ttributed to Kirchoff

Section : PLATE BENDING ELEMENTS Plte Geometr nd Deformtion The bsic derivtion of the field equtions for thin pltes strts ith the geometr shon in the figure belo: The plte surfces re t nd dthe mid-surfce is t ± t

Section : PLATE BENDING ELEMENTS Note tht t << t << b c nd if t is more thn bout one-tenth the spn of the plte, then trnsverse sher deformtions must be ccounted for nd the plte is sid to be thick thick. Also note tht the trnsverse deflection is much smller thn the thickness t, i.e., t << These re things tht should be checked hen running finite element nlsis ith thin plte elements. Anlogous to bems there re severl tpes of plte theories ssocited ith geometric ti ssumptions.

Section : PLATE BENDING ELEMENTS Kirchoff Assumptions nd Theoreticl Formultion Consider differentil slice from plte ith cutting plnes perpendiculr to the -is s shon belo. The lod q from the previous figure cuses the plte to deform lterll in the positive direction. The deflection t n point P in the plte is functionll dependent on its position reltive e to the - plne, i.e.,, ( ) nd the plte does not stretch in the -direction. The line -b drn perpendiculr to the mid-surfce before loding remins perpendiculr to the mid-surfce fter loding.

Section : PLATE BENDING ELEMENTS This is consistent ith the Kirchoff ssumptions hich re stted s follos:. Stright lines norml to the mid surfce remin norml. This implies tht Hoever γ γ γ tht is right ngles in the plne of the plte do not necessril remin right ngles fter the plte is loded.. Thickness chnges re neglected nd lines norml to the mid surfce do not undergo etension. This implies ε Free surfce Too thin to hve Norml stresses re considered negligible. σ resonble Vrition in σ. Membrne forces re neglected. Plne stress solutions re superimposed lter in the plne of the membrne. In-plne deformtions in the nd directions t the mid surfce re ero. u (,,) v(,,)

Section : PLATE BENDING ELEMENTS The ssumption regrding: γ γ cn be interpretted grphicll in the folloing to figures: γ b rigid bod rottion γ γ γ γ γ γ γ One cn esil see in the second figure tht normls do not remin norml in the presence of γ strin.

Section : PLATE BENDING ELEMENTS The sme grphicl interprettion cn be mde regrding γ γ γ γ γ γ

Section : PLATE BENDING ELEMENTS For Kirchoff pltes displcements in the -direction re ssumed to vr linerl from the mid surfce: ( ) u,, As result strins in the direction re chrcteried b the prtil differenctil eqution to the right: u ε

Section : PLATE BENDING ELEMENTS ( ) v,, Simrlrl for displcements in the direction: hich leds to: v ε Finll, the in plne sher strin is given b the epression u v v u γ

Section : PLATE BENDING ELEMENTS The curvture of the pltes re defined s Thus the strins cn be epressed s ε γ ε ε γ

Section : PLATE BENDING ELEMENTS For pltes the ssumption of plne stress is mde in the direction. Thus the constitutive reltionship is [ ] D ε ε σ σ here τ γ [ ] E D

Section : PLATE BENDING ELEMENTS Note tht the stresses vr linerl from the middle surfce, just like bending stresses in bems. Also note tht the sher stresses (τ ) produced b bending lso vr linerl from the middle surfce. The sher stresses τ nd τ re present nd required for equilibrium, lthough the corresponding strins re ssumed negligible. Prbolic vritions of the stresses re ssumed s in bem theor. The bending stresses cn be simplified to resultnt moments (M, M, M ). These moments re resultnts of the liner stress vritions through the thickness

Section : PLATE BENDING ELEMENTS The moments M,M,, M ( torsionl moment) s ell s M d M d the sher forces nd shon to the right re per unitl length quntities. M d d M d d

Section : PLATE BENDING ELEMENTS t d d d d M d M σ σ Here t d d M τ here ε ε γ E E ε ε σ σ γ τ net pge

Section : PLATE BENDING ELEMENTS t M σ Thus t d M M M τ σ σ t d E t ( ) Et

Section : PLATE BENDING ELEMENTS In similr fshion t d σ If e sum forces in the -direction nd sum moments bout the nd es ( nd M re lloed to vr cross the differentil element, hence the use of grdients in the figure belo t σ nd on the net pge) ssuming pressure q is pplied to the plte d d d d F qda d d d d d d qdd dd dd q

Section : PLATE BENDING ELEMENTS In similr fshion summing moments bout the -is leds to M d d d d dd d d ddd qdd M M M d M M d dd M M dd Similrl, summing moments bout the -is ields (sho for homeork) M M

Section : PLATE BENDING ELEMENTS Thus the equtions of equilibrium t point for plte re M M q M M These equilibrium equtions long ith the moment-curvture reltionships leds to the folloing non-homogenous fourth order, mied prtil differentil eqution in tht must be solved for pltes: solved for pltes: ( ) q Et If is knon, the strins re knon. If the strins re knon the stresses re knon.

Section : PLATE BENDING ELEMENTS Plte Element Formultion bsed on Kirchoff Theor In pper published in 98 Hrbok nd Hrudle cited 88 tpes of plte element formultions. In this section one is presented nd issues ssocited ith the formultion re discussed. Consider the bsic rectngulr plte element ith four nodes ech ith three degrees of freedom, i.e., one out of plne displcement () nd to in plne rottions ( nd ). This element is depicted belo. Shon in the figure re possible nodl forces (f i ) nd nodl moments (f i nd f i ): Keep in mind tht unknon deformtions (displcements nd rottions) re temed up ith pplied ctions (forces nd moments).

Section : PLATE BENDING ELEMENTS The nodl deformtion vector t n node i is identified s here { } d i i i i The negtive sign on is due to the fct tht negtive displcement () is required to produce positive rottion bout the -is (recll the right hnd rule). The totl element deformtion vector is {} { } T d

Section : PLATE BENDING ELEMENTS Becuse there re degrees of freedom for the rectngulr element term polnomil in nd is selected. A complete cubic polnomil hs terms. (, ) 7 8 9 5 6 to other terms We need to more terms from the qurtic polnomil line of Pscl s polnomil tringle. To understnd ho to pick these terms e need to discuss Continuit Conditions. To tpes on continuit re considered. The re clssified in the folloing mnner: C Continuit - The element is C comptible if long n side of n element is completel specified b the degrees of freedom on long tht side. Thus no gps pper long the edge here to elements meet. C Continuit - The element is C comptible if it is C comptible, nd the slope the side is completel specified b the degrees of freedom on long tht side.

Section : PLATE BENDING ELEMENTS With C continuit the slopes from to elements ill mtch long n edge here the elements meet. This hppens in the folloing depiction of elements vieed on edge: Where in the follo figure slopes do not mtch long the edges of elements here the meet. With Kirchoff plte theor, the slopes must mtch, giving the smoother description of the deformed surfce in the first figure.

Section : PLATE BENDING ELEMENTS If e dd nd terms s the other terms in the polnomil for then e get to the required terms: (, ) 7 8 The element stisfies the trnsverse deflection continuit. Along n edge beteen to nodes, s the edge beteen nodes i nd j in the figure belo 9 5 6 the deflection ill vr in cubic mnner nd, derivtive of, ill vr qudrticll.

Section : PLATE BENDING ELEMENTS Tht is for reltive to the locl element coordinte is, then ( ) 7, ( ) 7, The constnts,, nd 7 cn be uniquel determined long this edge of the element b four degrees of freedom long tht edge: i, j, i nd j. Consequentl C continuit ill hold nd pprentl C is stisfied. But ( ), 8 5

Section : PLATE BENDING ELEMENTS is cubic in, hich requires more boundr conditions to find, 5, 8 nd. But onl to more deformtion conditions eist long this edge of the element, i.e., i nd j. This g g,, i j slope is not uniquel defined nd slope discontinuit cn occur long this edge. Hence this element does not stisf C continuit conditions. This does not men tht the constnts through cn not be determined. To see this g consider tht M

Section : PLATE BENDING ELEMENTS With 7 6 5 9 8 or { } [ ]{ } C d

Section : PLATE BENDING ELEMENTS Finding the inverse of this reltionship cn be trivill stted s { } [ ] { } d C With { } [ ] { } d C here { }{ } P here { } P

Section : PLATE BENDING ELEMENTS Then or { P}[ C] {} d here { N }{ d } { N} { P}[ C] is the mtri of shpe funtions.

Section : PLATE BENDING ELEMENTS Returning the ttention to the strins in plte element γ ε ε 8 7 6 6 γ 9 8 5 9 6 8 7 6 6 6 6 We cn rite { } { }{ } { } { }{ } { }[ ] { } d C B B ε

Section : PLATE BENDING ELEMENTS here 6 6 { } 6 6 6 6 B nd { } { }[ ] C B B The element stiffness mtri is c b [ ] [ ] [ ] [ ] d d B D B t k c b T T

Section : PLATE BENDING ELEMENTS There re mn different s of interpolting deformtions over plte, leding to mn different tpes of plte elements hving been proposed for Kirchoff s plte theor. The sme issue comes up for Mindlin-Reisner plte theor. Some proposed elements for Kirchoff s theor hve encountered problems ith locking - refusl of the element to deform in some circumstnces, s the set of equtions becomes over-constrined, due to inconsistent ssumptions being mde. This is referred to s sher locking. Other elements hve hd the opposite problem of lcking stiffness to resist certin eroenerg modes of deformtion. These spurious modes of deformtion then dominte the solution. Another nme for this is hourglss deformtion, s the unnted in-plne deformtion of qudrtic element resembles this shpe.

Section : PLATE BENDING ELEMENTS