Plastic Anisotropy: Relaxed Constraints, Theoretical Textures

Similar documents
Chapter 6: Plastic Theory

STRESS UPDATE ALGORITHM FOR NON-ASSOCIATED FLOW METAL PLASTICITY

Reference material Reference books: Y.C. Fung, "Foundations of Solid Mechanics", Prentice Hall R. Hill, "The mathematical theory of plasticity",

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

MECHANICS OF MATERIALS. EQUATIONS AND THEOREMS

Characterizations of Aluminum Alloy Sheet Materials Numisheet 2005

Lecture #7: Basic Notions of Fracture Mechanics Ductile Fracture

Lecture 8. Stress Strain in Multi-dimension

Bone Tissue Mechanics

Constitutive Equations

Chapter 2: Elasticity

An Algorithm for Determining the Active Slip Systems in Crystal Plasticity Based on the Principle of Maximum Dissipation

28 GEOMETRY. + III A tr(a 1 Ã B) 1, (1.4.7a) with

Analytical formulation of Modified Upper Bound theorem

1. Background. is usually significantly lower than it is in uniaxial tension

Constitutive models: Incremental plasticity Drücker s postulate

Texture, Microstructure & Anisotropy A.D. (Tony) Rollett

* Many components have multiaxial loads, and some of those have multiaxial loading in critical locations

ANSYS Mechanical Basic Structural Nonlinearities

Lecture #8: Ductile Fracture (Theory & Experiments)

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

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

1 Stress and Strain. Introduction

Continuum Mechanics. Continuum Mechanics and Constitutive Equations

Study of Axes Rotation during Simple Shear Tests on Aluminum Sheets

INCREASING RUPTURE PREDICTABILITY FOR ALUMINUM

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

arxiv: v1 [cond-mat.mtrl-sci] 23 May 2017

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading

Lecture #10: Anisotropic plasticity Crashworthiness Basics of shell elements

DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS

Fatigue and Fracture

6.37 Determine the modulus of resilience for each of the following alloys:

Rheology. What is rheology? From the root work rheo- Current: flow. Greek: rhein, to flow (river) Like rheostat flow of current

ELASTOPLASTICITY THEORY by V. A. Lubarda

Module-4. Mechanical Properties of Metals

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

Classical fracture and failure hypotheses

Inverse identification of plastic material behavior using. multi-scale virtual experiments

Lecture contents. Stress and strain Deformation potential. NNSE 618 Lecture #23

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

INTRODUCTION TO STRAIN

CHAPTER 6 MECHANICAL PROPERTIES OF METALS PROBLEM SOLUTIONS

6.4 A cylindrical specimen of a titanium alloy having an elastic modulus of 107 GPa ( psi) and

Multiaxial Fatigue. Professor Darrell F. Socie. Department of Mechanical Science and Engineering University of Illinois at Urbana-Champaign

Multiaxial Fatigue Life Prediction Methods for Notched Bars of Ti-6Al-4V ABSTRACT

MAE 322 Machine Design Lecture 2. Dr. Hodge Jenkins Mercer University

A rate-dependent Hosford-Coulomb model for predicting ductile fracture at high strain rates

Ádx 1 + Áu 1 + u 1 ˆ

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

3D Stress Tensors. 3D Stress Tensors, Eigenvalues and Rotations

Mechanical Properties

ME 2570 MECHANICS OF MATERIALS

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

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

Rock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth

N = N A Pb A Pb. = ln N Q v kt. = kt ln v N

MHA042 - Material mechanics: Duggafrågor

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

Outline. Tensile-Test Specimen and Machine. Stress-Strain Curve. Review of Mechanical Properties. Mechanical Behaviour

A Simple and Accurate Elastoplastic Model Dependent on the Third Invariant and Applied to a Wide Range of Stress Triaxiality

Finite Element Method in Geotechnical Engineering

(MPa) compute (a) The traction vector acting on an internal material plane with normal n ( e1 e

Theory of Plasticity. Lecture Notes

MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD

A Critical Plane-energy Model for Multiaxial Fatigue Life Prediction. of Homogeneous and Heterogeneous Materials. Haoyang Wei

Loading σ Stress. Strain

The objective of this experiment is to investigate the behavior of steel specimen under a tensile test and to determine it's properties.

ON THE USE OF HOMOGENEOUS POLYNOMIALS TO DEVELOP ANISOTROPIC YIELD FUNCTIONS WITH APPLICATIONS TO SHEET FORMING

Elastic Properties of Solid Materials. Notes based on those by James Irvine at

Continuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms

Modelling the behaviour of plastics for design under impact

Constitutive model for quasi-static deformation of metallic sandwich cores

IAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.

CONSIDERATIONS CONCERNING YIELD CRITERIA INSENSITIVE TO HYDROSTATIC PRESSURE

The Finite Element Method II

CHAPTER 3 THE EFFECTS OF FORCES ON MATERIALS

MODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Isotropic Elastic Models: Invariant vs Principal Formulations

06 - kinematic equations kinematic equations

Materials and Structures. Indian Institute of Technology Kanpur

Failure from static loading

Material parameter identification for the numerical simulation of deep-drawing drawing of aluminium alloys

MATERIAL MODELS FOR THE

FE FORMULATIONS FOR PLASTICITY

COMPARISON OF THE TESTS CHOSEN FOR MATERIAL PARAMETER IDENTIFICATION TO PREDICT SINGLE POINT INCREMENTAL FORMING FORCES

By drawing Mohr s circle, the stress transformation in 2-D can be done graphically. + σ x σ y. cos 2θ + τ xy sin 2θ, (1) sin 2θ + τ xy cos 2θ.

GG303 Lab 12 11/7/18 1

MECHANICS OF MATERIALS

Chapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature

Mechanics of Earthquakes and Faulting

Lecture Notes #10. The "paradox" of finite strain and preferred orientation pure vs. simple shear

Linear Elasticity ( ) Objectives. Equipment. Introduction. ε is then

Onedimensional ratchetting Twodimensional ratchetting Response of unified models Response of 2M1C model Response of polycrystalline models

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

Single Crystal Gradient Plasticity Part I

2 Constitutive Models of Creep

Non-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises

Failure surface according to maximum principal stress theory

Strain Transformation equations

MATERIAL MECHANICS, SE2126 COMPUTER LAB 2 PLASTICITY

Transcription:

1 Plastic Anisotropy: Relaxed Constraints, Theoretical Textures Texture, Microstructure & Anisotropy Last revised: 11 th Oct. 2016 A.D. Rollett

2 The objective of this lecture is to complete the description of plastic anisotropy. Vectorization of stress and strain Definition of Taylor factor Comparison of single with multiple slip Description of Relaxed Constraints

References Kocks, Tomé & Wenk: Texture & Anisotropy (Cambridge); chapter 8, 1996. Detailed analysis of plastic deformation and texture development. Reid: Deformation Geometry for Materials Scientists, 1973. Older text with many nice worked examples. Be careful of his examples of calculation of Taylor factor because, like Bunge & others, he does not use von Mises equivalent stress/strain to obtain a scalar value from a multiaxial stress/strain state. Hosford: The Mechanics of Crystals and Textured Polycrystals, 1993 (Oxford). Written from the perspective of a mechanical metallurgist with decades of experimental and analytical experience in the area. Khan & Huang: Continuum Theory of Plasticity. Written from the perspective of continuum mechanics. De Souza Neto, Peric & Owen: Computational Methods for Plasticity, 2008 (Wiley). Written from the perspective of continuum mechanics. 3

4 Vectorization of stress, strain The lack of dependence on hydrostatic stress and volume constancy permits a different vectorization to be used (as compared with Voigt/matrix notation for anisotropic elasticity). The following set of basis tensors provides a systematic approach. The first tensor, b (1) represents deviatoric tension on Z, second plane strain compression in the Z plane, and the second row are the three simple shears in the {XYZ} system. The first two tensors provide a basis for the π-plane.

5 Vectorization: 2 " b (1) 1 1 0 0% $ ' $ 0 1 0' 6 $ # 0 0 2 ' & b (3) 1 " 0 0 0% $ ' $ 0 0 1' 2 $ ' # 0 1 0& b (4) 1 " 0 0 1% $ ' $ 0 0 0' 2 $ ' # 1 0 0& " b (2) 1 1 0 0% $ ' $ 0 1 0' 2 $ # 0 0 0 ' & b (5) 1 " 0 1 0% $ ' $ 1 0 0' 2 $ ' # 0 0 0& The components of the vectorized stresses and strains are then contractions (projections) with the basis tensors: σ λ = σ:b (λ), ε λ = ε:b (λ) There is, in fact, a 6 th eigentensor that separates out the hydrostatic components of stress and strain

6 Vectorization: 3 This preserves work conjugacy, i.e. σ:ε = Σ λ σ λ ε λ. The Lequeu vectorization scheme was almost the same as this but with the first two components interchanged. It is also useful to see the vector forms in terms of the regular tensor components. $ S λ = 2σ σ σ 33 11 22 & % 6 $ D λ = 2D D D 33 11 22 & % 6, σ 22 σ 11 2, D 22 D 11 2 ', 2σ 23, 2σ 31, 2σ 12 ) ( ', 2D 23, 2D 31, 2D 12 ) (

7 Maximum work, summary The sum of the shears for the actual set of active systems is less than any hypothetical set (Taylor s hypothesis). This also shows that we have obtained an upper bound on the stress required to deform the crystal because we have approached the solution for the work rate from above: any hypothetical solution results in a larger work rate than the actual solution.

8 Taylor factor, M We can take ratios of stresses, or strains to define the Taylor factor, M; if a simple choice of deformation axes is made, such as uniaxial tension, then the indices can be dropped to obtain the typical form of the equation, σ = <M>τ crss. For the polycrystal, the arithmetic mean of the Taylor factors is typically used to represent the ratio between the macroscopic flow stress and the critical resolved shear stress. This relies on the assumption that weak (low Taylor factor) grains cannot deform much until the harder ones are also deforming plastically, and that the harder (high Taylor factor) grains are made to deform by a combination of stress concentration and work hardening around them.

9 Taylor factor, multiaxial stress For multiaxial stress states, one may use the effective stress, e.g. the von Mises stress (defined in terms of the stress deviator tensor, S=σ-(σ ii /3), and also known as effective stress). Note that the equation below provides the most self-consistent approach for calculating the Taylor factor for multiaxial deformation because it ties everything back to the yield stress in uniaxial tension. M = σ vm τ σ vonmises σ vm = = s Δγ (s) 3 2 S : S dε vm = dw τ c dε vm = σ : dε τ c dε vm

10 Taylor factor, multiaxial strain Similarly for the von Mises equivalent strain increment (where dε p is the plastic strain increment which has zero trace, i.e. dε ii =0). dε vonmises dε vm = 2 3 dε p : dε p = 2 3 1 2 dε ij : dε ij = " 2 % $ '{( dε 11 dε 22 ) 2 + ( dε 22 dε 33 ) 2 + ( dε 33 dε 11 ) 2 }+ 4 # 9 & 3 dε 2 23 + dε 2 + dε 2 31 12 M = σ vm τ = s Δγ (s) dε vm = dw τ c dε vm { } = σ : dε τ c dε vm Compare with single slip: Schmid factor = cosφcosλ = τ/σ Formula for von Mises strain corrected 11 x 16

11 Uniaxial compression/tension So, for axisymmetric straining paths, the orientation dependence within a Standard Stereographic Triangle (SST) is as follows. The velocity gradient has the form (upper +/- means compression, lower +/- means tension):! # L = # # " ±δ 0 0 0 ±δ 0 0 0 2δ $ & & & % The von Mises equivalent strain for such a tensile strain is always 2δ.

12 Taylor factor(orientation) Hosford: mechanics of xtals...

13 Texture hardening Note that the Taylor factor is largest at both the 111 and 110 positions and a minimum at the 100 position. Thus a cubic material with a perfect <111> or <110> fiber will be 1.5 times as strong in tension or compression as the same material with a <100> fiber texture. This is not as dramatic a strengthening as can be achieved by other means, e.g. precipitation hardening, but it is significant. Also, it can be achieved without sacrificing other properties.

14 Uniaxial deformation: single slip Recall the standard picture of the single crystal under tensile load. In this case, we can define angles between the tensile direction and the slip plane normal, φ, and also between the tensile direction and the slip direction, λ. Given an applied tensile stress (force over area) on the crystal, we can calculate the shear stress resolved onto the particular slip system as τ=σ cosφ cosλ. This simple formula (think of using only the direction cosine for the slip plane and direction that corresponds to the tensile axis) then allows us to rationalize the variation of stress with testing angle (crystal orientation) with Schmid's Law concerning the existence of a critical resolved shear stress (CRSS).

15 Single slip: Schmid factors Re-write the relationship in terms of (unit) vectors that describe the slip plane, n, and slip direction, b: index notation and tensor notation are used interchangeably. τ = ˆ b σ ˆ n τ = σ ij b i n j

16 Single: multiple comparison It is interesting to consider the difference between multiple slip and single slip stress levels because of its relevance to deviations from the Taylor model. Hosford presents an analysis of the ratio between the stress required for multiple slip and the stress for single slip under axisymmetric deformation conditions (in {111}<110> slip).

17 Single: Multiple comparison Max. difference (=1.65) between single & multiple slip Figure 6.2 from Ch. 6 of Hosford, showing the ratio of the Taylor factor to the reciprocal Schmid factor, M/(1/m), for axisymmetric flow with {111}<110> slip. Orientations near 110 exhibit the largest ratios and might therefore be expected to deviate most readily from the Taylor model.

18 Single: Multiple comparison The result expressed as a ratio of the Taylor factor to the reciprocal Schmid factor is that orientations near 100 or 111 exhibit negligible differences whereas orientations near 110 exhibit the largest differences. This suggests that the latter orientations are the ones that would be expected to deviate most readily from the Taylor model. In wire drawing of bcc metals, this is observed*: the grains tend to deform in plane strain into flat ribbons. The flat ribbons then curl around each other in order to maintain compatibility. *W. F. Hosford, Trans. Met. Soc. AIME, 230 (1964), pp. 12.

19 Plane Strain Compression So, for plane strain straining paths, the velocity gradient has the form: L = # ±δ 0 0 & % ( 0 0 0 % ( $ 0 0 δ' The von Mises equivalent strain for such a tensile strain is always 2/ 3δ = 1.155δ. Note that this value leads to different results for the Taylor factor, compared to examples in, e.g., Reid and in Bunge (but is consistent with the definitions in the LApp and VPSC codes).

20 Relaxed Constraints An important modification of the Taylor model is the Relaxed Constraints (RC) model. It is important because it improves the agreement between experimental and calculated textures. The model is based on the development of large aspect ratios in grain shape with increasing strain. It makes the assumption that certain components of shear strain generate displacements in volumes that are small enough that they can be neglected. Thus any shear strain developed parallel to the short direction in an elongated grain (e.g. in rolling) produces a negligible volume of overlapping material with a neighboring grain.

21 Relaxed Constraints, contd. In the case of rolling, both the ε 13 and the ε 23 shears produce negligible compatibility problems at the periphery of the grain. Thus the is a relaxation of the strict enforcement of compatibility inherent in the Taylor model.

22 Relaxed Constraints, contd. Full constraints vs. Relaxed constraints D C 11 = s s D C 22 = D C 23 = D C 31 = D C 12 = s s s γ (s) (s) m 11 γ (s) (s) m 22 γ (s) m 23 (s) γ (s ) (s ) m 31 γ (s ) (s) m 12 D C 11 = s γ (s) (s) m 11 D C 22 = γ (s) (s) m 22 σ C 23 = 0 σ C 31 = 0 D C 12 = s s γ (s ) (s) m 12 5 strain BCs mixed strain & stress BCs

23 (rate sensitive) With only 3 boundary conditions on the strain rate, the same equation must be satisfied but over fewer components. D C = ε 0 s m (s) :σ c τ (s) n (s ) m (s) sgn m (s) :σ c ( ) Note: the Bishop-Hill maximum work method can still be applied to find the operative stress state: one uses the 3-fold vertices instead of the 6- and 8-fold vertices of the single crystal yield surface.

24 Relaxed Constraints Despite the crude nature of the RC model, experience shows that it results in superior prediction of texture development in both rolling and torsion. It is only an approximation! Better models, such as the LAMEL model, the self-consistent model (and finiteelement models) account for grain shape more accurately.

25 Effect of RC on texture development In (fcc) rolling, the stable orientation approaches the Copper instead of the Taylor/Dillamore position.

26 Y.S. for textured polycrystal FC RC Kocks: Ch.10 Note sharp vertices for strong textures at large strains. Yield surfaces based on highly elongated grains and the

27 Summary Vectorization of stress, strain tensors. Definition, explanation of the Taylor factor. Comparison of single and multiple slip. Relaxed Constraints model with mixed boundary conditions, allowing for the effect of grain shape on anisotropy.