The. Consortium. Continuum Mechanics. Original notes by Professor Mike Gunn, South Bank University, London, UK Produced by the CRISP Consortium Ltd

Similar documents
8.7 Associated and Non-associated Flow Rules

Geo-E2010 Advanced Soil Mechanics L Wojciech Sołowski. 07 March 2017

FE FORMULATIONS FOR PLASTICITY

Computational Solid Mechanics Computational Plasticity

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

Constitutive models: Incremental plasticity Drücker s postulate

ME111 Instructor: Peter Pinsky Class #21 November 13, 2000

Lecture Triaxial Stress and Yield Criteria. When does yielding occurs in multi-axial stress states?

GEO E1050 Finite Element Method Mohr-Coulomb and other constitutive models. Wojciech Sołowski

Analysis of cold rolling a more accurate method

A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE

Pressure-sensitivity Effects on Toughness Measurements of Compact Tension Specimens for Strain-hardening Solids

INTRODUCING THE SHEAR-CAP MATERIAL CRITERION TO AN ICE RUBBLE LOAD MODEL

Planar Transformations and Displacements

Lower bound solutions for bearing capacity of jointed rock

Adam Paweł Zaborski. 8 Plasticity. reloading. 1. Bauschinger s effect. 2. unchanged yielding limit. 3. isotropic hardening

PLASTIC ANALYSIS 1.0 INTRODUCTION

SAMARIS. Draft report

Classical fracture and failure hypotheses

Micro I. Lesson 5 : Consumer Equilibrium

ELASTO-PLASTIC ANALYSIS OF CORRUGATED SANDWICH STEEL PANELS

University of Pretoria Department of Mechanical & Aeronautical Engineering MOW 227, 2 nd Semester 2014

Churilova Maria Saint-Petersburg State Polytechnical University Department of Applied Mathematics

Elio Sacco. Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e LM

0.1 Practical Guide - Surface Integrals. C (0,0,c) A (0,b,0) A (a,0,0)

A Simple And Efficient FEM-Implementation Of The Modified Mohr-Coulomb Criterion Clausen, Johan Christian; Damkilde, Lars

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

Conference Paper A Finite Element Approach for the Elastic-Plastic Behavior of a Steel Pipe Used to Transport Natural Gas

REFINED STRAIN ENERGY OF THE SHELL

Plasticity R. Chandramouli Associate Dean-Research SASTRA University, Thanjavur

MECHANICS OF MATERIALS

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

GBT (GENERALISED BEAM THEORY)-BASED ELASTIC - PLASTIC POST-BUCKLING ANALYSIS OF STAINLESS STEEL THIN-WALLED MEMBERS

Implementation of an Isotropic Elastic-Viscoplastic Model for Soft Soils using COMSOL Multiphysics

Problem " Â F y = 0. ) R A + 2R B + R C = 200 kn ) 2R A + 2R B = 200 kn [using symmetry R A = R C ] ) R A + R B = 100 kn

The effect of dynamic bending moments on the ratchetting behavior of stainless steel pressurized piping elbows

HSNV122 - Thermoplasticity and metallurgy in large deformations in simple tension

Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity

FEM for elastic-plastic problems

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

Concept Question Comment on the general features of the stress-strain response under this loading condition for both types of materials

THE GENERAL ELASTICITY PROBLEM IN SOLIDS

STRESS-STRAIN-DILATANCY RELATIONSHIPS OF NORMALLY CONSOLIDATED DHAKA CLAY

A General Damage Initiation and Evolution Model (DIEM) in LS-DYNA

Tutorial #1 - CivE. 205 Name: I.D:

σ = 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.

Lecture Thermodynamics 9. Entropy form of the 1 st law. Let us start with the differential form of the 1 st law: du = d Q + d W

KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS

Lesson 16 Electric Potential

V-0. Review of Probability

Rational Expressions. Definition: P, is an algebraic expression that can be written as the quotient of two polynomials, P. A rational expression, Q

PLATE GIRDERS - I 1.0 INTRODUCTION

CONTENTS. Lecture 1 Introduction. Lecture 2 Physical Testing. Lecture 3 Constitutive Models

An elasto-plastic model to describe the undrained cyclic behavior of saturated sand with initial static shear

Analysis of Masonry Vaulted Systems: The Barrel Vaults

BENCHMARK 2 Cup Drawing of Anisotropic Thick Steel Sheet

Algorithmic Formulations of Evolutionary Anisotropic Plasticity Models Based on Non-Associated Flow Rule

A NUMERICAL MODEL OF FINITE DIFFERENCE (F.D) FOR DYNAMIC PILE DRIVING

Numerical Analysis Differential Methods

STRUCTURAL ENGINEERING AND GEOMECHANICS - Modeling For Geotechnical Engineering Applications - Pradeep U. Kurup

UNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 20, 2011 Professor A. Dolovich

Find the equation of a plane perpendicular to the line x = 2t + 1, y = 3t + 4, z = t 1 and passing through the point (2, 1, 3).

Stability Analysis of a Geometrically Imperfect Structure using a Random Field Model

SEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by

CMSC 425: Lecture 4 Geometry and Geometric Programming

MODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Uniaxial Model: Strain-Driven Format of Elastoplasticity

4 Strain true strain engineering strain plane strain strain transformation formulae

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

HÅLLFASTHETSLÄRA, LTH Examination in computational materials modeling

10.2 Polar Equations and Graphs

3D and 2D Formulations of Incremental Stress-Strain Relations for Granular Soils

Stable ion beam transport through periodic electrostatic structures: linear and non-linear effects

A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials

Instability of seabed under wave induced loading

DEFORMATION THEORY OF PLASTICITY

SHAPE OPTOMIZATION OF H-BEAM FLANGE FOR MAXIMUM PLASTIC ENERGY DISSIPATION

ANSYS Mechanical Basic Structural Nonlinearities

ME 323 Examination #2 April 11, 2018

Characterization of Material Parameters

Module 2 Stresses in machine elements. Version 2 ME, IIT Kharagpur

Elio Sacco. Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e LM

Casimir Force Between the Two Moving Conductive Plates.

6.1 The Linear Elastic Model

Elasto-Plastic Behaviour of Thick Plates with a Higher-Order Shear Deformation Theory

Numerical Modeling of Direct Shear Tests on Sandy Clay

Finite Element Analysis of V-Bending of Polypropylene Using Hydrostatic-Pressure-Dependent Plastic Constitutive Equation*

7.6 Stress in symmetrical elastic beam transmitting both shear force and bending moment

Analysis of tunnel and super structures for excavation

MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design

Chapter 6: Plastic Theory

Tensor Transformations and the Maximum Shear Stress. (Draft 1, 1/28/07)

torsion equations for lateral BucKling ns trahair research report r964 July 2016 issn school of civil engineering

CMSC 425: Lecture 7 Geometric Programming: Sample Solutions

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

Failure mechanisms and corresponding shape factors of shallow foundations

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

Modelling the behaviour of plastics for design under impact

Figure 1: General Plane Motion (Translation and Rotation)

Chapter 15: Visualizing Stress and Strain

1.1 The Equations of Motion

Transcription:

The C R I S P Consortium Continuum Mechanics Original notes b Professor Mike Gunn, South Bank Universit, London, UK Produced b the CRISP Consortium Ltd

THOR OF STRSSS In a three dimensional loaded bod, there are si indeendent comonents of stress at a oint Shear stress notation: ab face of cube direction B moment euilibrium:,, In an stressed material there are alwas three mutuall erendicular lanes on which the shear stresses are ero. These are called the rincial lanes. The direct stresses acting on these lanes are called rincial stresses (,, ). c In a stressed bod the comonents of stress are,,,, and. We can also describe the stress state with resect to another set of aes (e.g. a, b, c). b a

Stress invariants are functions of the stress comonents which are indeendent of the ais sstem chose. For eamle, c b a is a stress invariant. Another stress invariant is where: ( ) ( ) ( ) ( ) ( ) ( ) 6 6 6 6 6 6 ca bc ab a c c b b a Note that these general definitions of and reduce to those given later for triaial stress conditions (see lecture on Cam-cla). We will encounter these invariants more freuentl eressed in terms of the rincial stresses, i.e. ( ) ( ) ( ) Note,, are also stress invariants. There is a simle geometric interretation of and in rincial stress sace. The line is called the sace diagonal or the hdrostatic ais OA ( ) AB AB B A O (,, )

LASTICIT Generalised Hooke s Law: ε ν ν ε ν ν ε ν ν γ γ γ ( ν) ( ν ) ( ν ) PLASTICIT INTRODUCTION The theor of elasticit allows the calculation of stresses and strains in a loaded bod when the bod is linear and elastic. The theor of lasticit allows the calculation of stresses and strains in a loaded bod when lastic ielding takes lace. Histor (lasticit) 8 Cauch invents concets of stress and strain 885 Boussines gives his solutions for elastic stress distributions (still used in geotechnical engineering toda) (For 00 ears mathematicians slave awa obtaining elasticit solutions...) 956 Finite element method invented (mathematicians are now redundant) 4

Histor (Plasticit) 77 Coulomb identifies two comonents in the strength of soil - cohesion and friction 864 Criteria are ut forward for (Tresca) limits to the elastic behaviour of metals 9 (Von Mises) 940 s tes of statement are now seen to be necessar to comletel describe lastic stress-strain relations: a) ield criteria b) flow rule c) hardening law 950 R. Hill s book Mathematical Theor of Plasticit is ublished. PLASTIC BHAVIOUR FOR ON- DIMNSIONAL LOADING Direct strain in direction is ε. B D Tical stress -strain relation for an elastic, work hardening lastic material (e.g. metal allo) A O C OA is elastic. A is a ield oint. is the uniaial ield stress. BC and D are elastic unloading and reloading (arallel to OA). On reloading to B the ield stress has increased to. The material is harder and has strain (or work) hardened. 5

In describing lastic behaviour, the following simlifications ( idealisations ) are often made. lastic, erfectl lastic (steel behaviour is uite like this) lastic, linear strain hardening lastic Rigid-lastic (often assumed in collase calculations) ILD FUNCTIONS - BASIC IDAS is held constant and is increased until ielding starts (or vice versa) The combinations of and that cause ielding are described b a ield function. ield functions can be reresented b lines in D stress sace and surfaces in D stress sace. 6

For elastic-erfectl lastic behaviour: LASTIC Plastic ielding on this ield curve or surface Stress states inside the ield surface are elastic. Stress states outside the ield surface are imossible, b definition. If strain hardening takes lace, there are two ossibilities: A B A B C O O C The ield surface eands uniforml - this is called isotroic hardening A B A B C ' O O C' The ield surface is dragged along - this is called kinematic hardening. The differences between these two assumtions become imortant if unloading takes lace. 7

ILD FUNCTIONS - XAMPLS a) TRSCA ielding takes lace when the maimum value of -, -, - euals a critical value (k). This can be interreted (b considering Mohr s circles) as being euivalent to limiting the maimum shear stress on an lane in the material to being less than or eual to k. k k Consider the ield condition in (, ) stress sace when 0. k k k k k For triaial stress conditions and we can lot the ield function in (, ) sace. k k k 8

b) VON MISS Plastic ielding takes lace when: ( ) ( ) ( ) where is the ield stress in uniaial tension. Note (from our earlier definition of ) that the above euation can be written as: Structural and mechanical engineers call the stress arameter either: i) the Von Mises stress, ii) the Von Mises euivalent stress, or iii) the euivalent stress and use the smbol rather than. Von Mises ield function in (, ) sace when 0 : This is an ellise with its major ais along. The Tresca ield criterion is shown thus ----- for comarison. Clearl k. In triaial stress conditions and the ield function in (, ) sace is as shown: STRAINS For ever stress comonent or invariant there is a corresonding strain comonent so that the work done er unit volume in elastic deformation is / stress strain. W.D ε 9

Stress Strain ε ε ε γ γ γ ε ε ε Stress invariant Strain invariant (defined reviousl) ε (volumetric strain) ε ε ε ε ε ε (defined reviousl) ε (deviatoric strain) ε ( ε ε ) ( ε ε ) ( ε ε ) ( γ γ γ ) ε ( ε ε ) ( ε ε ) ( ε ε ) 0

PLASTIC STRAINS e The lastic strain is the strain which remains on comletel unloading the alied stress. ε ε e ε T ε lastic strain ε e elastic strain ε T total strain CALCULATION OF LASTIC STRAINS Suose that the soil is in some stress state ( [,, ] T ). An increment of stress is now alied to the soil ( ). The resulting strains, if the behaviour is elastic, can be calculated as ε C where C is a suare matri containing elastic constants. CALCULATION OF PLASTIC STRAINS In contrast to elastic behaviour, it is found that lastic strains are strongl deendent on the current total stresses. C Thus the lastic strains generated a ield oint C are not deendent on the direction of the stress ath when ielding starts. A X O B X i.e. aths OAC and OBC lead to the same lastic strains.

PLASTIC POTNTIALS Mathematicall, the te of behaviour described above can convenientl be eressed in terms of a otential function such that the derivatives of the otential function define the ratios of the lastic strains. Plastic strain comonents are often lotted in stress sace:, Current stress state often called the strain increment vector (, ), Mathematicall we write:, constant otential G d ε m () where G is the lastic otential, For metals the lastic otential function is the same as the ield function i.e. G F and hence we can write: F m d ε () (m is a scalar number which deends on the hardening law and the details of a articular analsis. It is called the lastic multilier). when G F we sa there is normalit (the strain increment vector is normal to the ield surface) - also called associated flow. uations like () and () above are called flow rules. The govern the ratios of the lastic strain comonents.

Flow Rule - amle with Von Mises (a) in, sace when 0 F (, ) 0 dε dε F m m F m m ( ) ( ), dε X X, dε B insection the strain increment vector is erendicular to the ellise reviousl described. (b) in, sace, dε ( ) ielding takes lace with ero volumetric strains dε 0 lastic P, dε

() The Von Mises ield Surface c Hdrostatic ais b a () The Tresca ield Surface c Hdrostatic ais b a 4

() The Drucker-Prager ield Surface c Hdrostatic ais b a (4) The Mohr-Coulomb ield Surface c Hdrostatic ais b a 5