MEG 741 Energy and Variational Methods in Mechanics I

Size: px
Start display at page:

Download "MEG 741 Energy and Variational Methods in Mechanics I"

Transcription

1 MEG 74 Energ and Variational Methods in Mechanics I Brendan J. O Toole, Ph.D. Associate Professor of Mechanical Engineering Howard R. Hghes College of Engineering Universit of Nevada Las Vegas TBE B- (7) bj@me.nlv.ed Chapter : Differential Form of Basic Eqations -

2 Part Vectors, Tensors, and Matrices -

3 Scalar, Vector & Tensor Eamples Speed (Scalar or th order Tensor) It has magnitde onl and is represented b a single vale. Velocit (Vector or st order tensor) It has magnitde and direction, 3 components are needed to fll define a vector, e.g.: v (v, v, v ). Stress ( nd order tensor) Stress at a point has different magnitdes in different directions. Nine components are needed to define stress. σ σ σ σ σ σ σ σ σ σ -3

4 Sbscript Notation, σ, σ First Sbscript Signifies the ais to which the otward normal of the face is parallel Second Sbscript Indicates the ais to which the stress is parallel σ is a stress acting on the -face in the -direction which implies that it is a shear stress. σ is a stress acting on the -face in the -direction which implies that it is a normal stress. The second sbscript is sometimes omitted for normal stresses. If onl one sbscript is given, assme that it is a normal stress. -4

5 Sign Convention for Stresses on a Control Element of an Elastic Bod σ τ τ τ τ σ τ τ σ σ τ τ τ τ σ τ τ σ Positive Stress Directions On Positive Faces Positive Stress Directions On Negative Faces -5

6 Sign Convention For Stresses Positive Faces of Control Element A positive face has an otward normal vector in the positive direction of an ais. Stresses are positive when the are in the positive direction of an ais Negative Faces of Control Element A negative face has an otward normal vector in the negative direction of an ais. Stresses are positive when the are in the negative direction of an ais -6

7 Sbscript Notation for Stresses Normal Stresses (σ or σ ) Shear Stresses (τ or σ ) The Greek smbol σ is sall sed to represent a normal stress and τ is sed to represent a shear stress. When working with vector or tensor eqations it is necessar to se one smbol to represent all tpes of stresses. In general, σ is sed to represent all stresses and the sbscripts designate whether it is a normal stress or shear stress. -7

8 Tensors & Indicial Notation Scalars and vectors are sbsets of the general definition of a tensor qantit. Higher order tensors have more components. Einstein developed a shorthand notation to simplif writing tensor eqations. With indicial notation, the sbscripts i and j are sed to represent all aes (,, and ). For eample, σ ij is a shorthand version of writing all nine components of the entire stress tensor. -8

9 Eamples Using Indicial Notation Vectors ( st order tensors) The vector (v, v, v ) can be written as v i or v j, the sbscript represents all three components. nd order tensors (Stress and Strain) The stress (or strain) at a point is represented b nine different vectors. σ ij is the shorthand method for writing all nine components where i,, and j,,. ij where, and are normal strains and the rest are shear strains; e.g: γ /. -9

10 Smmation Convention Component Form of Vector A A a e a e a 3 e 3, where e, e, and e 3 are base vectors Smmation Notation for Vector A 3 j j A a e j a e j Shorthand Notation for Vector A j k m A a e a e a e j k m j -

11 Second Order Tensor Component Form Inditial Notation Form Φ φ e e φ e e φ 3 e e 3 φ e e φ e e φ 3 e e 3 φ 3 e 3 e φ 3 e 3 e φ 33 e 3 e 3 Φ φ ij e i e j Unit Second Order Tensor, I I δ ij e i e j δ ij if if i i j j -

12 Dot Prodct Between Second Order Tensors (Doble Dot Prodct) Similar to single dot prodct, encontered periodicall in mechanics derivations Φ : Ψ (φ ij e i e j ) : (ψ mn e m e n ) φ ij ψ mn (e j. e m ) (e i. e n ) φ ij ψ mn δ jm δ in φ im ψ mn δ in φ nm ψ mn φ ij ψ ji -

13 Matrices and Matri Operations The components of a tensor can be written down in a rectanglar arra (or matri). Sqare matri: n n (has order n) Diagonal Elements: elements with the same row and colmn nmber Diagonal Matri: Onl non-ero elements are along the diagonal Identit Matri [I]: Diagonal matri with all non-ero elements -3

14 Matrices and Matri Operations Trace of a Matri: Sm of the diagonal elements Smmetric Matri: [A] [A] T Skew Smmetric Matri: [A] T - [A] Review on or own: Matri addition, mltiplication Calclating determinants, Cramer s Rle Matri Transformations -4

15 Part Theor of Elasticit -5

16 Governing Eqations For Elastic Bodies Kinematics Strain-Displacement Eqations Thermodnamics st and nd Laws of Thermodnamics (Compatibilit Eqations) Constittive Eqations Stress-Strain Relations Kinetics Conservation of momenta -6

17 Kinematics Kinematic eqations are sed to describe the motion and deformation of deformable bodies Lagrangian (or referential) Description The motion of the bod, B, is described relative to a reference configration C R (sall the nstressed state of the bod). The crrent coordinates (,, ) are epressed in terms of the reference coordinates (,, ), which are sometimes called the material coordinates, and time, t. Spatial (or Elerian) Description The motion is referred to the crrent configration of the bod. The spatial description focses on a region of space, instead of on a given bod. Elerian descriptions are sed most often in flid mechanics and is sometimes referred to as a control volme approach. -7

18 Lagrangian Strain, C O r Q r P Q ds P r Q P r P Q,, C Q ds P A bod ndergoes deformation, P and Q represent particles which occp the initial positions defined b r P & r Q. The final position of the two particles are represented b r P & r Q. The particles are separated b the infinitesimal distance, ds, in the ndeformed state and ds, in the deformed state. -8

19 Strains The strain tensor, є, is defined in terms of the change in distance between the two particles while moving from the ndeformed state to the deformed state. It can be shown that these distances are related to the displacements, P and Q. є is known as the Green-Lagrange strain tensor ik ( ) i, k k, i l, i l, k -9

20 - Cartesian Components of Strain (Green Lagrange Strain Tensor) w w v v w v w w v v w w w v v v w v w w v v w v displacement in the direction v displacement in the direction w displacement in the direction

21 Infinitesimal Strain Green-Lagrange strain tensor: ik ( ) i, k k, i l, i l, k For small strains, the higher order terms are negligible The redced form of the strain tensor is called the infinitesimal or Cach strain tensor: ik ( ) i, k k, i -

22 - Cartesian Components of Strain (Infinitesimal or Cach Strain Tensor) w w w v w v w v,,,,,,,,, γ γ γ displacement in the direction v displacement in the direction w displacement in the direction

23 -3 Matri Form of the Cach Strain Tensor D γ γ γ D is the differential operator matri

24 Compatibilit Eqations Assming we have epressions for the displacements and these epressions are differentiable, we can determine the strains sing the previosl defined eqations. However, if the si strains are known and we wish to determine the three displacements, o mst be sre that a niqe soltion eists. The compatibilit conditions ensre that a niqe set of displacements are fond and that there are no kinks, tears, gaps, or overlaps. Some of the strain components mst be related to ensre compatibilit. -4

25 -5 Compatibilit Relations The compatibilit relations are defined b considering derivatives of individal strain components. For eample, take the partial derivative of γ with respect to and : v v γ γ γ v γ Since Rearranging order of partial derivatives Sbstitting definitions of strain This is strain compatibilit eqation. There are 6 of these for a general 3-D problem.

26 Compatibilit Relations The si compatibilit relations can be written sing indicial notation as: ij mn im m n i For -D problems, this redces to a single eqation: j j n i jn m See Eq..8 for the si eqations. -6

27 Constittive Relations (Stress - Strain Eqations) The stress-strain eqations or constittive relations are a good eample of the seflness of the indicial notation. The most basic constittive relation that we first learn in Mechanics of Materials is the -D Hookes Law eqation: σ E, Stress eqals Yong's modls times strain. Later, we learned that if there is strain in more than one direction, the stress will be a fnction of all strains. -7

28 General 3-D Constittive Relations The first two (of nine) stress eqations can be written as shown below. Each of the nine stress components is a fnction of all nine strain components. Each Q variable is a different fnction of material properties. σ Q Q Q Q Q Q Q Q Q σ Q Q Q Q Q Q Q Q Q -8

29 General 3-D Constittive Relations Both the shorthand and fll matri versions of the 3-D constittive eqations are shown below. σ ij Q ijkl kl σ σ σ σ σ σ σ σ σ where i, j, k and l represent, and. Sometimes the nmbers, and 3 are sbstitted for, and. Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q -9

30 Simplifing the Constittive Relations The general constittive eqations have 8 elastic constants. Lckil, the nmber of constants is redced for most practical materials. Both the stress and strain tensors are smmetric. σ ij σ ji and ij ji This redces the nmber of independent stresses and strains to 6. -3

31 Constittive Relations Since there are onl si independent stresses and strains, another shorthand notation is introdced sing nmeric sbscripts. Let σ σ and σ σ and σ 3 σ and 3 σ 4 σ σ and 4 σ 5 σ σ and 5 σ 6 σ σ and 6-3

32 Constittive Relations This redces the nmber of independent material constants to 36. σ σ σ3 σ4 σ5 σ6 Q Q Q3 Q4 Q5 Q6 Q Q Q3 Q4 Q5 Q6 Q3 Q3 Q33 Q34 Q35 Q36 Q4 Q4 Q43 Q44 Q45 Q46 Q5 Q5 Q53 Q54 Q55 Q56 Q6 Q6 Q63 Q64 Q65 Q66 A thermodnamics proof can be sed to show that the Q matri itself is smmetric, or Q ij Q ji. This redces the nmber of independent constants to

33 Anisotrop and Material Smmetr General Anisotrop A material with different material properties in all directions ehibits general anisotrop. independent elastic constants are reqired to define the stress-strain relationship for this tpe of material. Orthotropic Material Has 3 mtall orthogonal planes of elastic smmetr e.g.: A material with the same properties in the and - directions has elastic smmetr abot the - plane. Most composite materials ehibit elastic smmetr abot three planes Reqires 9 independent elastic constants -33

34 Orthotropic Stress-Strain Eqations σ σ σ3 σ4 σ 5 σ6 or σ i Q ij j, where Q ij is know as the stiffness matri. Inverting the matri eqation ields: Q Q Q3 Q Q Q3 Q3 Q3 Q33 Q44 Q 55 Q66 S S S3 S S S3 S3 S3 S33 S44 S 55 S66 or i S ij σ j, where S ij is know as the compliance matri. [S] [Q] - σ σ σ3 σ4 σ 5 σ

35 Compliance Matri The compliance matri is easier to define in terms of elastic properties S S S3 S S S3 S3 S3 S33 S44 S 55 S66 σ σ σ3 σ4 σ 5 σ6 S E S 44 G 3 S E S 55 G 3 S 33 E 33 S 66 G S S ν E ν E S 3 S 3 ν 3 E ν 3 E 33 S 3 S 3 ν 3 E ν 3 E 33-35

36 Stiffness Matri TheQ ij components are fond b inverting S ij Q E ( ν 3 ν 3 ) Q 44 G 3 Q E ( ν 3 ν 3 ) Q 55 G 3 Q 33 E 33 ( ν ν ) Q 66 G Q E ( ν ν 3 ν 3 ) E ( ν ν 3 ν 3 ) Q 3 E ( ν 3 ν ν 3 ) E 33 ( ν 3 ν ν 3 ) Q 3 E ( ν 3 ν ν 3 ) E 33 ( ν 3 ν ν 3 ) ν ν ν 3 ν 3 ν 3 ν 3 ν ν 3 ν 3-36

37 -37 Isotropic Stress-Strain Relations ( )( ) E γ γ γ ν ν ν ν ν ν ν ν ν ν ν ν τ τ τ σ σ σ ν ν

Chapter 1: Differential Form of Basic Equations

Chapter 1: Differential Form of Basic Equations MEG 74 Energ and Variational Methods in Mechanics I Brendan J. O Toole, Ph.D. Associate Professor of Mechanical Engineering Howard R. Hghes College of Engineering Universit of Nevada Las Vegas TBE B- (7)

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

Concept of Stress at a Point

Concept of Stress at a Point Washkeic College of Engineering Section : STRONG FORMULATION Concept of Stress at a Point Consider a point ithin an arbitraril loaded deformable bod Define Normal Stress Shear Stress lim A Fn A lim A FS

More information

3 2D Elastostatic Problems in Cartesian Coordinates

3 2D Elastostatic Problems in Cartesian Coordinates D lastostatic Problems in Cartesian Coordinates Two dimensional elastostatic problems are discssed in this Chapter, that is, static problems of either plane stress or plane strain. Cartesian coordinates

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

FEA Solution Procedure

FEA Solution Procedure EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations

More information

8 Properties of Lamina

8 Properties of Lamina 8 Properties of Lamina 8- ORTHOTROPIC LAMINA An orthotropic lamina is a sheet with unique and predictable properties and consists of an assemblage of fibers ling in the plane of the sheet and held in place

More information

EE2 Mathematics : Functions of Multiple Variables

EE2 Mathematics : Functions of Multiple Variables EE2 Mathematics : Fnctions of Mltiple Variables http://www2.imperial.ac.k/ nsjones These notes are not identical word-for-word with m lectres which will be gien on the blackboard. Some of these notes ma

More information

Formal Methods for Deriving Element Equations

Formal Methods for Deriving Element Equations Formal Methods for Deriving Element Eqations And the importance of Shape Fnctions Formal Methods In previos lectres we obtained a bar element s stiffness eqations sing the Direct Method to obtain eact

More information

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University 9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)

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

And similarly in the other directions, so the overall result is expressed compactly as,

And similarly in the other directions, so the overall result is expressed compactly as, SQEP Tutorial Session 5: T7S0 Relates to Knowledge & Skills.5,.8 Last Update: //3 Force on an element of area; Definition of principal stresses and strains; Definition of Tresca and Mises equivalent stresses;

More information

4.4 Moment of a Force About a Line

4.4 Moment of a Force About a Line 4.4 Moment of a orce bot a Line 4.4 Moment of a orce bot a Line Eample 1, page 1 of 3 1. orce is applied to the end of gearshift lever DE. Determine the moment of abot shaft. State which wa the lever will

More information

Geometry of Span (continued) The Plane Spanned by u and v

Geometry of Span (continued) The Plane Spanned by u and v Geometric Description of Span Geometr of Span (contined) 2 Geometr of Span (contined) 2 Span {} Span {, } 2 Span {} 2 Geometr of Span (contined) 2 b + 2 The Plane Spanned b and If a plane is spanned b

More information

Professor Terje Haukaas University of British Columbia, Vancouver The M4 Element. Figure 1: Bilinear Mindlin element.

Professor Terje Haukaas University of British Columbia, Vancouver   The M4 Element. Figure 1: Bilinear Mindlin element. Professor Terje Hakaas University of British Colmbia, ancover www.inrisk.bc.ca The M Element variety of plate elements exist, some being characterized as Kirchhoff elements, i.e., for thin plates, and

More information

Differentiation of Exponential Functions

Differentiation of Exponential Functions Differentiation of Eponential Fnctions The net derivative rles that o will learn involve eponential fnctions. An eponential fnction is a fnction in the form of a constant raised to a variable power. The

More information

A Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane

A Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane Adv. Theor. Appl. Mech., Vol. 1, 8, no. 1, 9 A Decomposition Method for Volme Flx and Average Velocit of Thin Film Flow of a Third Grade Flid Down an Inclined Plane A. Sadighi, D.D. Ganji,. Sabzehmeidani

More information

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2 MATH 307 Vectors in Rn Dr. Neal, WKU Matrices of dimension 1 n can be thoght of as coordinates, or ectors, in n- dimensional space R n. We can perform special calclations on these ectors. In particlar,

More information

1. Solve Problem 1.3-3(c) 2. Solve Problem 2.2-2(b)

1. Solve Problem 1.3-3(c) 2. Solve Problem 2.2-2(b) . Sole Problem.-(c). Sole Problem.-(b). A two dimensional trss shown in the figre is made of alminm with Yong s modls E = 8 GPa and failre stress Y = 5 MPa. Determine the minimm cross-sectional area of

More information

MATHEMATICAL FUNDAMENTALS I. Michele Fitzpatrick

MATHEMATICAL FUNDAMENTALS I. Michele Fitzpatrick MTHEMTICL FUNDMENTLS I Michele Fitpatrick OVERVIEW Vectors and arras Matrices Linear algebra Del( operator Tensors DEFINITIONS vector is a single row or column of numbers. n arra is a collection of vectors

More information

Incompressible Viscoelastic Flow of a Generalised Oldroyed-B Fluid through Porous Medium between Two Infinite Parallel Plates in a Rotating System

Incompressible Viscoelastic Flow of a Generalised Oldroyed-B Fluid through Porous Medium between Two Infinite Parallel Plates in a Rotating System International Jornal of Compter Applications (97 8887) Volme 79 No., October Incompressible Viscoelastic Flow of a Generalised Oldroed-B Flid throgh Poros Medim between Two Infinite Parallel Plates in

More information

Chapter 2 Introduction to the Stiffness (Displacement) Method. The Stiffness (Displacement) Method

Chapter 2 Introduction to the Stiffness (Displacement) Method. The Stiffness (Displacement) Method CIVL 7/87 Chater - The Stiffness Method / Chater Introdction to the Stiffness (Dislacement) Method Learning Objectives To define the stiffness matrix To derive the stiffness matrix for a sring element

More information

CS 450: COMPUTER GRAPHICS VECTORS SPRING 2016 DR. MICHAEL J. REALE

CS 450: COMPUTER GRAPHICS VECTORS SPRING 2016 DR. MICHAEL J. REALE CS 45: COMPUTER GRPHICS VECTORS SPRING 216 DR. MICHEL J. RELE INTRODUCTION In graphics, we are going to represent objects and shapes in some form or other. First, thogh, we need to figre ot how to represent

More information

Complex Variables. For ECON 397 Macroeconometrics Steve Cunningham

Complex Variables. For ECON 397 Macroeconometrics Steve Cunningham Comple Variables For ECON 397 Macroeconometrics Steve Cnningham Open Disks or Neighborhoods Deinition. The set o all points which satis the ineqalit

More information

Motion in an Undulator

Motion in an Undulator WIR SCHAFFEN WISSEN HEUTE FÜR MORGEN Sven Reiche :: SwissFEL Beam Dnamics Grop :: Pal Scherrer Institte Motion in an Undlator CERN Accelerator School FELs and ERLs On-ais Field of Planar Undlator For planar

More information

Section 7.4: Integration of Rational Functions by Partial Fractions

Section 7.4: Integration of Rational Functions by Partial Fractions Section 7.4: Integration of Rational Fnctions by Partial Fractions This is abot as complicated as it gets. The Method of Partial Fractions Ecept for a few very special cases, crrently we have no way to

More information

Chapter 2 Difficulties associated with corners

Chapter 2 Difficulties associated with corners Chapter Difficlties associated with corners This chapter is aimed at resolving the problems revealed in Chapter, which are cased b corners and/or discontinos bondar conditions. The first section introdces

More information

10.4 Solving Equations in Quadratic Form, Equations Reducible to Quadratics

10.4 Solving Equations in Quadratic Form, Equations Reducible to Quadratics . Solving Eqations in Qadratic Form, Eqations Redcible to Qadratics Now that we can solve all qadratic eqations we want to solve eqations that are not eactl qadratic bt can either be made to look qadratic

More information

Chapter 6 Momentum Transfer in an External Laminar Boundary Layer

Chapter 6 Momentum Transfer in an External Laminar Boundary Layer 6. Similarit Soltions Chapter 6 Momentm Transfer in an Eternal Laminar Bondar Laer Consider a laminar incompressible bondar laer with constant properties. Assme the flow is stead and two-dimensional aligned

More information

Period #5: Strain. A. Context. Structural mechanics deals with the forces placed upon mechanical systems and the resulting deformations of the system.

Period #5: Strain. A. Context. Structural mechanics deals with the forces placed upon mechanical systems and the resulting deformations of the system. Period #5: Strain A. Contet Strctral mechanics deals with the forces placed pon mechanical sstems and the reslting deformations of the sstem. Solid mechanics on the smaller scale relates stresses in 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

Classify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled.

Classify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled. Jnction elements in network models. Classify by nmber of ports and examine the possible strctres that reslt. Using only one-port elements, no more than two elements can be assembled. Combining two two-ports

More information

Symmetry and Properties of Crystals (MSE638) Stress and Strain Tensor

Symmetry and Properties of Crystals (MSE638) Stress and Strain Tensor Symmetry and Properties of Crystals (MSE638) Stress and Strain Tensor Somnath Bhowmick Materials Science and Engineering, IIT Kanpur April 6, 2018 Tensile test and Hooke s Law Upto certain strain (0.75),

More information

2.2 Relation Between Mathematical & Engineering Constants Isotropic Materials Orthotropic Materials

2.2 Relation Between Mathematical & Engineering Constants Isotropic Materials Orthotropic Materials Chapter : lastic Constitutive quations of a Laminate.0 Introduction quations of Motion Symmetric of Stresses Tensorial and ngineering Strains Symmetry of Constitutive quations. Three-Dimensional Constitutive

More information

Section 1: Review of Elasticity

Section 1: Review of Elasticity Finite Elements in Elasticit Section : Review of Elasticit Stress & Strain 2 Constitutive Theor 3 Energ Methods Section : Stress and Strain Stress at a point Q : F = lim A 0 A F = lim A 0 A F = lim A 0

More information

σ = F/A. (1.2) σ xy σ yy σ zy , (1.3) σ xz σ yz σ zz The use of the opposite convention should cause no problem because σ ij = σ ji.

σ = F/A. (1.2) σ xy σ yy σ zy , (1.3) σ xz σ yz σ zz The use of the opposite convention should cause no problem because σ ij = σ ji. Cambridge Universit Press 978-0-521-88121-0 - Metal Forming: Mechanics Metallurg, Third Edition Ecerpt 1 Stress Strain An understing of stress strain is essential for analzing metal forming operations.

More information

σ = F/A. (1.2) σ xy σ yy σ zx σ xz σ yz σ, (1.3) The use of the opposite convention should cause no problem because σ ij = σ ji.

σ = F/A. (1.2) σ xy σ yy σ zx σ xz σ yz σ, (1.3) The use of the opposite convention should cause no problem because σ ij = σ ji. Cambridge Universit Press 978-1-107-00452-8 - Metal Forming: Mechanics Metallurg, Fourth Edition Ecerpt 1 Stress Strain An understing of stress strain is essential for the analsis of metal forming operations.

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

1. State-Space Linear Systems 2. Block Diagrams 3. Exercises

1. State-Space Linear Systems 2. Block Diagrams 3. Exercises LECTURE 1 State-Space Linear Sstems This lectre introdces state-space linear sstems, which are the main focs of this book. Contents 1. State-Space Linear Sstems 2. Block Diagrams 3. Exercises 1.1 State-Space

More information

16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations

16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations 6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =

More information

Digital Image Processing. Lecture 8 (Enhancement in the Frequency domain) Bu-Ali Sina University Computer Engineering Dep.

Digital Image Processing. Lecture 8 (Enhancement in the Frequency domain) Bu-Ali Sina University Computer Engineering Dep. Digital Image Processing Lectre 8 Enhancement in the Freqenc domain B-Ali Sina Uniersit Compter Engineering Dep. Fall 009 Image Enhancement In The Freqenc Domain Otline Jean Baptiste Joseph Forier The

More information

Primary dependent variable is fluid velocity vector V = V ( r ); where r is the position vector

Primary dependent variable is fluid velocity vector V = V ( r ); where r is the position vector Chapter 4: Flids Kinematics 4. Velocit and Description Methods Primar dependent ariable is flid elocit ector V V ( r ); where r is the position ector If V is known then pressre and forces can be determined

More information

A FOUR-NODED PLANE ELASTICITY ELEMENT BASED ON THE SEPARATION OF THE DEFORMATION MODES

A FOUR-NODED PLANE ELASTICITY ELEMENT BASED ON THE SEPARATION OF THE DEFORMATION MODES A FOUR-NODD LAN LASICIY LMN BASD ON H SARAION OF H DFORMAION MODS A. DÓSA D. RADU Abstract: his paper presents a for-noded qadrilateral finite element with translational degrees of freedom for plane elasticit

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

Homework 5 Solutions

Homework 5 Solutions Q Homework Soltions We know that the colmn space is the same as span{a & a ( a * } bt we want the basis Ths we need to make a & a ( a * linearly independent So in each of the following problems we row

More information

Constitutive Equations

Constitutive Equations Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module

More information

TENSOR TRANSFORMATION OF STRESSES

TENSOR TRANSFORMATION OF STRESSES GG303 Lecture 18 9/4/01 1 TENSOR TRANSFORMATION OF STRESSES Transformation of stresses between planes of arbitrar orientation In the 2-D eample of lecture 16, the normal and shear stresses (tractions)

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

Numerical verification of the existence of localization of the elastic energy for closely spaced rigid disks

Numerical verification of the existence of localization of the elastic energy for closely spaced rigid disks Nmerical verification of the existence of localization of the elastic energy for closely spaced rigid disks S. I. Rakin Siberian State University of transport Rssia, 6349, Novosibirsk, Dsy Kovalchk street,

More information

Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS

Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS 3. System Modeling Mathematical Modeling In designing control systems we mst be able to model engineered system dynamics. The model of a dynamic system

More information

Partial Differential Equations with Applications

Partial Differential Equations with Applications Universit of Leeds MATH 33 Partial Differential Eqations with Applications Eamples to spplement Chapter on First Order PDEs Eample (Simple linear eqation, k + = 0, (, 0) = ϕ(), k a constant.) The characteristic

More information

Geometric Image Manipulation. Lecture #4 Wednesday, January 24, 2018

Geometric Image Manipulation. Lecture #4 Wednesday, January 24, 2018 Geometric Image Maniplation Lectre 4 Wednesda, Janar 4, 08 Programming Assignment Image Maniplation: Contet To start with the obvios, an image is a D arra of piels Piel locations represent points on the

More information

Time-adaptive non-linear finite-element analysis of contact problems

Time-adaptive non-linear finite-element analysis of contact problems Proceedings of the 7th GACM Colloqim on Comptational Mechanics for Yong Scientists from Academia and Indstr October -, 7 in Stttgart, German Time-adaptive non-linear finite-element analsis of contact problems

More information

Crystal Micro-Mechanics

Crystal Micro-Mechanics Crystal Micro-Mechanics Lectre Classical definition of stress and strain Heng Nam Han Associate Professor School of Materials Science & Engineering College of Engineering Seol National University Seol

More information

Reduction of over-determined systems of differential equations

Reduction of over-determined systems of differential equations Redction of oer-determined systems of differential eqations Maim Zaytse 1) 1, ) and Vyachesla Akkerman 1) Nclear Safety Institte, Rssian Academy of Sciences, Moscow, 115191 Rssia ) Department of Mechanical

More information

Lecture 5. Differential Analysis of Fluid Flow Navier-Stockes equation

Lecture 5. Differential Analysis of Fluid Flow Navier-Stockes equation Lectre 5 Differential Analsis of Flid Flo Naier-Stockes eqation Differential analsis of Flid Flo The aim: to rodce differential eqation describing the motion of flid in detail Flid Element Kinematics An

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

The Dual of the Maximum Likelihood Method

The Dual of the Maximum Likelihood Method Department of Agricltral and Resorce Economics University of California, Davis The Dal of the Maximm Likelihood Method by Qirino Paris Working Paper No. 12-002 2012 Copyright @ 2012 by Qirino Paris All

More information

Chem 4501 Introduction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics. Fall Semester Homework Problem Set Number 10 Solutions

Chem 4501 Introduction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics. Fall Semester Homework Problem Set Number 10 Solutions Chem 4501 Introdction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics Fall Semester 2017 Homework Problem Set Nmber 10 Soltions 1. McQarrie and Simon, 10-4. Paraphrase: Apply Eler s theorem

More information

Pulses on a Struck String

Pulses on a Struck String 8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a

More information

Elements of Coordinate System Transformations

Elements of Coordinate System Transformations B Elements of Coordinate System Transformations Coordinate system transformation is a powerfl tool for solving many geometrical and kinematic problems that pertain to the design of gear ctting tools and

More information

LINEAR COMBINATIONS AND SUBSPACES

LINEAR COMBINATIONS AND SUBSPACES CS131 Part II, Linear Algebra and Matrices CS131 Mathematics for Compter Scientists II Note 5 LINEAR COMBINATIONS AND SUBSPACES Linear combinations. In R 2 the vector (5, 3) can be written in the form

More information

The Heat Equation and the Li-Yau Harnack Inequality

The Heat Equation and the Li-Yau Harnack Inequality The Heat Eqation and the Li-Ya Harnack Ineqality Blake Hartley VIGRE Research Paper Abstract In this paper, we develop the necessary mathematics for nderstanding the Li-Ya Harnack ineqality. We begin with

More information

PLATE AND PANEL STRUCTURES OF ISOTROPIC, COMPOSITE AND PIEZOELECTRIC MATERIALS, INCLUDING SANDWICH CONSTRUCTION

PLATE AND PANEL STRUCTURES OF ISOTROPIC, COMPOSITE AND PIEZOELECTRIC MATERIALS, INCLUDING SANDWICH CONSTRUCTION PLATE AND PANEL STRUCTURES OF ISOTROPIC, COMPOSITE AND PIEZOELECTRIC MATERIALS, INCLUDING SANDWICH CONSTRUCTION SOLID MECHANICS AND ITS APPLICATIONS Volume 10 Series Editor: G.M.L. GLADWELL Department

More information

Constitutive Relations

Constitutive Relations Constitutive Relations Andri Andriyana, Ph.D. Centre de Mise en Forme des Matériaux, CEMEF UMR CNRS 7635 École des Mines de Paris, 06904 Sophia Antipolis, France Spring, 2008 Outline Outline 1 Review of

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

Research Article Equivalent Elastic Modulus of Asymmetrical Honeycomb

Research Article Equivalent Elastic Modulus of Asymmetrical Honeycomb International Scholarl Research Network ISRN Mechanical Engineering Volume, Article ID 57, pages doi:.5//57 Research Article Equivalent Elastic Modulus of Asmmetrical Honecomb Dai-Heng Chen and Kenichi

More information

Assignment Fall 2014

Assignment Fall 2014 Assignment 5.086 Fall 04 De: Wednesday, 0 December at 5 PM. Upload yor soltion to corse website as a zip file YOURNAME_ASSIGNMENT_5 which incldes the script for each qestion as well as all Matlab fnctions

More information

. D CR Nomenclature D 1

. D CR Nomenclature D 1 . D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the

More 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

3.4-Miscellaneous Equations

3.4-Miscellaneous Equations .-Miscellaneos Eqations Factoring Higher Degree Polynomials: Many higher degree polynomials can be solved by factoring. Of particlar vale is the method of factoring by groping, however all types of factoring

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

Chapter 6 2D Elements Plate Elements

Chapter 6 2D Elements Plate Elements Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda

More information

APPENDIX B MATRIX NOTATION. The Definition of Matrix Notation is the Definition of Matrix Multiplication B.1 INTRODUCTION

APPENDIX B MATRIX NOTATION. The Definition of Matrix Notation is the Definition of Matrix Multiplication B.1 INTRODUCTION APPENDIX B MAIX NOAION he Deinition o Matrix Notation is the Deinition o Matrix Mltiplication B. INODUCION { XE "Matrix Mltiplication" }{ XE "Matrix Notation" }he se o matrix notations is not necessary

More information

OPTI-502 Optical Design and Instrumentation I John E. Greivenkamp Final Exam In Class Page 1/16 Fall, 2013

OPTI-502 Optical Design and Instrumentation I John E. Greivenkamp Final Exam In Class Page 1/16 Fall, 2013 OPTI-502 Optical Design and Instrmentation I John E. Greivenkamp Final Exam In Class Page 1/16 Fall, 2013 Name Closed book; closed notes. Time limit: 120 mintes. An eqation sheet is attached and can be

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

Elements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004

Elements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004 Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic

More information

Second-Order Wave Equation

Second-Order Wave Equation Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order

More information

DIAGONALIZATION OF THE STRESS TENSOR

DIAGONALIZATION OF THE STRESS TENSOR DIAGONALIZATION OF THE STRESS TENSOR INTRODUCTION By the use of Cauchy s theorem we are able to reduce the number of stress components in the stress tensor to only nine values. An additional simplification

More information

Setting The K Value And Polarization Mode Of The Delta Undulator

Setting The K Value And Polarization Mode Of The Delta Undulator LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions

More information

Get Solution of These Packages & Learn by Video Tutorials on Matrices

Get Solution of These Packages & Learn by Video Tutorials on  Matrices FEE Download Stud Package from website: wwwtekoclassescom & wwwmathsbsuhagcom Get Solution of These Packages & Learn b Video Tutorials on wwwmathsbsuhagcom Matrices An rectangular arrangement of numbers

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

OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIELD OF A POLYHEDRAL BODY WITH LINEARLY INCREASING DENSITY 1

OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIELD OF A POLYHEDRAL BODY WITH LINEARLY INCREASING DENSITY 1 OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIEL OF A POLYHERAL BOY WITH LINEARLY INCREASING ENSITY 1 V. POHÁNKA2 Abstract The formla for the comptation of the gravity field of a polyhedral body

More information

The Cross Product of Two Vectors in Space DEFINITION. Cross Product. u * v = s ƒ u ƒƒv ƒ sin ud n

The Cross Product of Two Vectors in Space DEFINITION. Cross Product. u * v = s ƒ u ƒƒv ƒ sin ud n 12.4 The Cross Prodct 873 12.4 The Cross Prodct In stdying lines in the plane, when we needed to describe how a line was tilting, we sed the notions of slope and angle of inclination. In space, we want

More information

EDEXCEL NATIONAL CERTIFICATE/DIPLOMA. PRINCIPLES AND APPLICATIONS of FLUID MECHANICS UNIT 13 NQF LEVEL 3 OUTCOME 3 - HYDRODYNAMICS

EDEXCEL NATIONAL CERTIFICATE/DIPLOMA. PRINCIPLES AND APPLICATIONS of FLUID MECHANICS UNIT 13 NQF LEVEL 3 OUTCOME 3 - HYDRODYNAMICS EDEXCEL NATIONAL CERTIFICATE/DIPLOMA PRINCIPLES AND APPLICATIONS of FLUID MECHANICS UNIT 3 NQF LEVEL 3 OUTCOME 3 - HYDRODYNAMICS TUTORIAL - PIPE FLOW CONTENT Be able to determine the parameters of pipeline

More information

AN ISOGEOMETRIC SOLID-SHELL FORMULATION OF THE KOITER METHOD FOR BUCKLING AND INITIAL POST-BUCKLING ANALYSIS OF COMPOSITE SHELLS

AN ISOGEOMETRIC SOLID-SHELL FORMULATION OF THE KOITER METHOD FOR BUCKLING AND INITIAL POST-BUCKLING ANALYSIS OF COMPOSITE SHELLS th Eropean Conference on Comptational Mechanics (ECCM ) 7th Eropean Conference on Comptational Flid Dynamics (ECFD 7) 5 Jne 28, Glasgow, UK AN ISOGEOMETRIC SOLID-SHELL FORMULATION OF THE KOITER METHOD

More information

Discontinuous Fluctuation Distribution for Time-Dependent Problems

Discontinuous Fluctuation Distribution for Time-Dependent Problems Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation

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

Homotopy Perturbation Method for Solving Linear Boundary Value Problems

Homotopy Perturbation Method for Solving Linear Boundary Value Problems International Jornal of Crrent Engineering and Technolog E-ISSN 2277 4106, P-ISSN 2347 5161 2016 INPRESSCO, All Rights Reserved Available at http://inpressco.com/categor/ijcet Research Article Homotop

More information

A State Space Based Implicit Integration Algorithm. for Differential Algebraic Equations of Multibody. Dynamics

A State Space Based Implicit Integration Algorithm. for Differential Algebraic Equations of Multibody. Dynamics A State Space Based Implicit Integration Algorithm for Differential Algebraic Eqations of Mltibody Dynamics E. J. Hag, D. Negrt, M. Ianc Janary 28, 1997 To Appear Mechanics of Strctres and Machines Abstract.

More information

Module #4. Fundamentals of strain The strain deviator Mohr s circle for strain READING LIST. DIETER: Ch. 2, Pages 38-46

Module #4. Fundamentals of strain The strain deviator Mohr s circle for strain READING LIST. DIETER: Ch. 2, Pages 38-46 HOMEWORK From Dieter 2-7 Module #4 Fundamentals of strain The strain deviator Mohr s circle for strain READING LIST DIETER: Ch. 2, Pages 38-46 Pages 11-12 in Hosford Ch. 6 in Ne Strain When a solid is

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

Lesson 81: The Cross Product of Vectors

Lesson 81: The Cross Product of Vectors Lesson 8: The Cross Prodct of Vectors IBHL - SANTOWSKI In this lesson yo will learn how to find the cross prodct of two ectors how to find an orthogonal ector to a plane defined by two ectors how to find

More information

9. Tensor product and Hom

9. Tensor product and Hom 9. Tensor prodct and Hom Starting from two R-modles we can define two other R-modles, namely M R N and Hom R (M, N), that are very mch related. The defining properties of these modles are simple, bt those

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

Bertrand s Theorem. October 8, µr 2 + V (r) 0 = dv eff dr. 3 + dv. f (r 0 )

Bertrand s Theorem. October 8, µr 2 + V (r) 0 = dv eff dr. 3 + dv. f (r 0 ) Bertrand s Theorem October 8, Circlar orbits The eective potential, V e = has a minimm or maximm at r i and only i so we mst have = dv e L µr + V r = L µ 3 + dv = L µ 3 r r = L µ 3 At this radis, there

More information

Concepts Introduced. Digital Electronics. Logic Blocks. Truth Tables

Concepts Introduced. Digital Electronics. Logic Blocks. Truth Tables Concepts Introdced Digital Electronics trth tables, logic eqations, and gates combinational logic seqential logic Digital electronics operate at either high or low voltage. Compters se a binary representation

More information

ECON3120/4120 Mathematics 2, spring 2009

ECON3120/4120 Mathematics 2, spring 2009 University of Oslo Department of Economics Arne Strøm ECON3/4 Mathematics, spring 9 Problem soltions for Seminar 4, 6 Febrary 9 (For practical reasons some of the soltions may inclde problem parts that

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

Localization in Undrained Deformation 1

Localization in Undrained Deformation 1 Localization in Undrained Deformation 1 J. W. Rdnicki Dept. of Civil and Env. Engn. and Dept. of Mech. Engn. Northwestern University Evanston, IL 6001-3109 John.Rdnicki@gmail.com Janary 7, 009 1 To appear

More information