Chapter 11. Mechanical Behavior of Materials

Similar documents
Chapter 11. Mechanical Behavior of Materials

Chapter Outline: Failure

Chapter Outline Dislocations and Strengthening Mechanisms. Introduction

Chapter 8 Strain Hardening and Annealing

Fracture. Brittle vs. Ductile Fracture Ductile materials more plastic deformation and energy absorption (toughness) before fracture.

Mechanical Properties

Chapter 7 Dislocations and Strengthening Mechanisms. Dr. Feras Fraige

MECHANICAL PROPERTIES AND TESTS. Materials Science

STRENGTHENING MECHANISM IN METALS

26. Irradiation Induced Mechanical Property Changes: Hardening and Embrittlement

Chapter Outline: Failure

Engineering materials

1) Fracture, ductile and brittle fracture 2) Fracture mechanics

ME -215 ENGINEERING MATERIALS AND PROCESES

Chapter Outline Mechanical Properties of Metals How do metals respond to external loads?

CHAPTER 4 1/1/2016. Mechanical Properties of Metals - I. Processing of Metals - Casting. Hot Rolling of Steel. Casting (Cont..)

26. Irradiation Induced Mechanical Property Changes: Hardening and Embrittlement

5. A round rod is subjected to an axial force of 10 kn. The diameter of the rod is 1 inch. The engineering stress is (a) MPa (b) 3.

ME 254 MATERIALS ENGINEERING 1 st Semester 1431/ rd Mid-Term Exam (1 hr)

IMPERFECTIONSFOR BENEFIT. Sub-topics. Point defects Linear defects dislocations Plastic deformation through dislocations motion Surface

Chapter 2: Mechanical Behavior of Materials

Chapter 7. Mechanical properties 7.1. Introduction 7.2. Stress-strain concepts and behaviour 7.3. Mechanical behaviour of metals 7.4.

Mechanical Properties of Metals. Goals of this unit

ENGR 151: Materials of Engineering LECTURE #12-13: DISLOCATIONS AND STRENGTHENING MECHANISMS

Introduction to Materials Science, Chapter 8, Failure. Failure. Ship-cyclic loading from waves.

CHAPTER 8 DEFORMATION AND STRENGTHENING MECHANISMS PROBLEM SOLUTIONS

Strengthening Mechanisms

Problems to the lecture Physical Metallurgy ( Materialkunde ) Chapter 6: Mechanical Properties

Fundamental Course in Mechanical Processing of Materials. Exercises

Creep and High Temperature Failure. Creep and High Temperature Failure. Creep Curve. Outline

Chapter 8. Deformation and Strengthening Mechanisms

Mechanical behavior of crystalline materials - Stress Types and Tensile Behaviour

Mechanical behavior of crystalline materials- Comprehensive Behaviour

Tensile/Tension Test Advanced Topics

21 Fracture and Fatigue Revision

MSE 3143 Ceramic Materials

Deformation, plastic instability

Engineering Materials

Movement of edge and screw dislocations

Introduction to Engineering Materials ENGR2000 Chapter 7: Dislocations and Strengthening Mechanisms. Dr. Coates

Mechanical Behaviour of Materials Chapter 10 Fracture morpholgy

Materials Engineering 272-C Fall 2001, Lectures 9 & 10. Introduction to Mechanical Properties of Metals

Assume that the growth of fatigue cracks in the plate is governed by a Paris type of law, i.e. da

CREEP CREEP. Mechanical Metallurgy George E Dieter McGraw-Hill Book Company, London (1988)

a. 50% fine pearlite, 12.5% bainite, 37.5% martensite. 590 C for 5 seconds, 350 C for 50 seconds, cool to room temperature.

Metal working: Deformation processing II. Metal working: Deformation processing II

3. Mechanical Properties of Materials

Introduction to Materials Science

Lecture 08 Fracture Toughness and Toughening Mechanisms Ref: Richerson, Modern Ceramic Engineering, Ch17, Marcel Dekker, 1992

MT 348 Outline No MECHANICAL PROPERTIES

Problem Set 2 Solutions

ENGINEERING SOLID MECHANICS Fundamentals and Applications

Fundamentals of Plastic Deformation of Metals

Quiz 1 - Mechanical Properties and Testing Chapters 6 and 8 Callister

Mechanical Properties of Materials

Introduction to Engineering Materials ENGR2000 Chapter 8: Failure. Dr. Coates

INGE Engineering Materials. Chapter 7: Part 2. Mechanical Failure. INGE Engineering Materials. Failure Analysis

High Temperature Materials. By Docent. N. Menad. Luleå University of Technology ( Sweden )

Ductile and Brittle Fracture of 1018 Steel and. 304 Stainless Steel Using Charpy Impact Test. Eman Mousa Alhajji. North Carolina State University

11/2/2018 7:58 PM. Chapter 6. Mechanical Properties of Metals. Mohammad Suliman Abuhaiba, Ph.D., PE

Structural design criteria

Types of Strain. Engineering Strain: e = l l o. Shear Strain: γ = a b

Chapter 14 Fracture Mechanics

PLASTIC DEFORMATION AND THE ONSET OF TENSILE INSTABILITY

Strengthening Mechanisms

Materials and their structures

Single vs Polycrystals

Homework 4 on Dislocations, Yield Stress, Hardness, Creep, Grain Size

Module-6. Dislocations and Strengthening Mechanisms

Subject Index. STP 1207-EB/Dec. 1994

REVISED PAGES IMPORTANT TERMS AND CONCEPTS REFERENCES QUESTIONS AND PROBLEMS. 166 Chapter 6 / Mechanical Properties of Metals

CE 221: MECHANICS OF SOLIDS I CHAPTER 3: MECHANICAL PROPERTIES OF MATERIALS

Engineering 45: Properties of Materials Final Exam May 9, 2012 Name: Student ID number:

FME201 Solid & Structural Mechanics I Dr.Hussein Jama Office 414

Chapter 4 MECHANICAL PROPERTIES OF MATERIAL. By: Ardiyansyah Syahrom

Part IA Paper 2: Structures and Materials MATERIALS Examples Paper 3 Stiffness-limited Design; Plastic Deformation and Properties

A THERMOMECHANICAL FATIGUE CRACK INITIATION MODEL FOR DIRECTIONALLY-SOLIDIFIED NI-BASE SUPERALLOYS

1-6.4 THE CRACK TIP: THE INGLIS EQUATION

Question Grade Maximum Grade Total 100

MECHANICAL PROPERTIES PROPLEM SHEET

High temperature applications

THE USE OF FIELD INDENTATION MICROPROBE IN MEASURING MECHANICAL PROPERTIES OF WELDS

A discrete dislocation plasticity analysis of grain-size strengthening

3. MECHANICAL PROPERTIES OF STRUCTURAL MATERIALS

Chapter 7: Mechanical Properties 1- Load 2- Deformation 3- Stress 4- Strain 5- Elastic behavior

Chapter 6 Mechanical Properties

Imperfections, Defects and Diffusion

MECHANICAL PROPERTIES OF MATERIALS

Creep failure Strain-time curve Effect of temperature and applied stress Factors reducing creep rate High-temperature alloys

Defects and Diffusion

Module #0. Introduction. READING LIST DIETER: Ch. 1, pp. 1-6

Concepts of stress and strain

Single-crystal Modeling of Ni-based Superalloys for Gas Turbine Blades

بسم الله الرحمن الرحیم. Materials Science. Chapter 7 Mechanical Properties

where n is known as strain hardening exponent.

Mechanical Properties

THE MECHANICAL PROPERTIES OF STAINLESS STEEL

The strength of a material depends on its ability to sustain a load without undue deformation or failure.

BFF1113 Engineering Materials DR. NOOR MAZNI ISMAIL FACULTY OF MANUFACTURING ENGINEERING

Materials Issues in Fatigue and Fracture. 5.1 Fundamental Concepts 5.2 Ensuring Infinite Life 5.3 Failure 5.4 Summary

Transcription:

Chapter 11. Mechanical Behavior of Materials 11.1 Introduction... 2 11.2. Mechanical Testing... 2 11.2.1. Uniaxial Tensile Loading Test... 2 11.2.2. Biaxial Testing: Tube Burst... 6 11.2.3. The Von Mises Equivalent Stress... 7 11.2.4. Hardness Testing... 9 11.2.5. Impact Testing... 10 11.3 Microstructural Aspects of Deformation... 12 11.3.1. Critical Resolved Shear Stress and the Yield Stress... 12 11.3.2. Dislocations as cause of work hardening... 14 11.3.3. Yield Strength Increase Mechanisms... 14 11.4. Creep Deformation... 15 Larson-Miller Plots... 17 11.4.1 Creep Mechanisms... 18 Thermally enhanced glide... 18 Thermally induced dislocation glide (Climb and glide mechanism)... 19 Stress Induced Diffusional Flow Difference (Nabarro-Herring Creep):... 19 Coble Creep:... 19 11.5. Material Fracture... 20 11.5.1. Ductile Failure by Diffuse Necking: The Considère Criterion... 21 11.5.2. Ductile Fracture... 22 11.5.3. Fracture due to Crack Growth... 24 Griffith Fracture Theory... 24 Problems... 28 References... 31 Light Water Reactor Materials Donald Olander and Arthur Motta 1

11.1 Introduction When in service, materials may be subjected to loads of various intensities, types and duration. The response of the material to these applied loads is termed the mechanical behavior of the material, and it is one of the most important factors to be considered for materials design. The most important questions to be answered are: How and when does the material undergo permanent deformation? When does the material fail or fracture? These apparently simple questions encompass the range of questions addressed in the field of mechanical behavior of materials. The questions have different answers depending on whether the load is applied quickly or over a period of time, whether the material is pulled in one or in more than one direction at once (uniaxial or multi-axial loading), whether the load is cyclic, and whether there are pre-existing flaws on the material. The ability of the material to resist deformation and failure under such conditions is measured by various mechanical properties, such as strength, ductility and toughness. These macroscopic material properties are discussed in this chapter as well as their link to the material microstructure. 11.2. Mechanical Testing Our quantitative knowledge of materials behavior comes from subjecting such materials to various mechanical tests and deriving measurable material properties. To instruct the discussion, various mechanical tests and their mechanical properties are described in the following 11.2.1. Uniaxial Tensile Loading Test In a tensile test a specimen of uniform cross section (such as a cylinder) is uniaxially loaded and its deformation measured as a function of the applied load (Figure 11.1). Light Water Reactor Materials Donald Olander and Arthur Motta 2

Figure 11.1 A schematic depiction of specimen deformation during a uniaxial tensile loading test. In figure 11.1, a specimen of initial gauge length l o and cross-sectional area A o is subjected to an increasing tensile load P, applied along its axis and the resulting specimen deformation is measured. Upon application of the load the specimen length increases to l and the cross-sectional area is reduced to A. The engineering stress is defined as: P σ eng = (11.1) A o but since during deformation the specimen cross sectional area is reduced from A o to A, the true stress is actually P σ true = (11.2) A The engineering strain is ε eng l l l o = (11.3) o However, the true strain is the integral of the increments of strain along the specimen length: ε true l dl l = = ln l l (11.4) l o o which is always somewhat larger than the engineering strain as long as deformation is uniform along gauge length. From volume conservation Light Water Reactor Materials Donald Olander and Arthur Motta 3

Al = A l (11.5) o o Note that ε eng is the strain in the loading direction; by Poisson s ratio this implies a strain of ε t =νε eng in the two transverse directions; if the material is isotropic, volume conservation yields ε t =-0.5ε eng. Although volume is conserved during plastic deformation, it is ot so during elastic deformation, for which a theoretically sisotropic material undergoing linerar elastiuc deformation should show ν=0.25. In reality most metals show ν~0.33. From (11.5) and (11.4) ε true A da Ao = = ln A A (11.6) Ao We are now ready to examine the results of a tensile test. If the applied stress is plotted against the specimen strain we obtain the stress-strain curve, a schematic example of which is shown in Figure 11.2. σ true-stress/true-strain engineering-stress/engineering-strain σ UTS σ y elastic energy ε uniform ε f ε Figure 11.2: Schematic stress-strain behavior for an austenitic steel. Light Water Reactor Materials Donald Olander and Arthur Motta 4

As load is applied the material initially starts to deform such that the strain is linearly proportional to the load, according to Hooke s Law, i.e. σ = Eε (11.7) where E is Young s modulus. The region of validity of Hooke s law is the elastic region. This region is characterized by small reversible deformation: that is, once the load is removed, the strain disappears. In that region, the engineering stress/engineering strain curve is coincident with the true stress/true strain curve. Some typical values of the elastic constants are shown in Table 11.1. Material E (Gpa) G(GPa) ν (Poisson's ratio) Aluminum 70.3 26.1 0.345 Gold 78 27 0.44 Iron 211.4 81.6 0.293 Nickel 199.5 76 0.312 Tungsten 411 160.6 0.28 Zirconium 1120 dir 99 32.43 (ave.) Zirconium 0001 dir 125 Table 11.1 Typical elastic and shear constants and Poisson s ratios of different materials Eventually the material deformation starts to deviate from Hooke s law as a result of plastic, or irreversible, deformation. In order to establish a definite onset of plastic deformation,, a minimum value of deviation from Hooke s law (normally 0.2%) is taken to be the point at which plastic deformation starts to occur. The stress at which this happens is called the yield stress, σ y, and is an important property of engineered materials. Beyond the yield point, the material starts to deform plastically. The stress required for further deformation continues to increase because of work hardening, (also called strain hardening) often described by an equation of the type: σ = Kε (11.8) n T where n is the work hardening coefficient and K is a constant, called the strength coefficient (note that the equation (11.8) describes only the work hardening region, starting at yield.) At any given strain the stress required to continue plastic deformation is called the flow stress. Thus the flow stress is equal to the yield stress at the yield point but is higher than the yield stress in the work hardening region. The material deformation occurring in the elastic region and in the beginning of the work hardening region is uniform, that is, all the material within the gauge section participates equally in the deformation process. At some point diffuse necking sets in, causing Light Water Reactor Materials Donald Olander and Arthur Motta 5

deformation to become non-uniform. A small variation in cross sectional area in a section of the material can cause a slightly larger stress, which causes larger deformation in turn further diminishing the area and so on. Once specimen necking starts, the material fails, because the cross-sectional area in the necked region becomes progressively smaller, and thus no additional load is needed to cause further deformation. The stress at which necking occurs is called the ultimate tensile strength of the material or σ UTS. This is the maximum load-bearing capacity of the material. The strain at which necking occurs is the uniform strain, ε uniform, and is one measure of the ductility of the material (i.e. how much can the material deform before failing). As deformation proceeds, the neck becomes progressively thinner until the material fractures. The strain at which this occurs is the fracture strain, ε f, and represents another measure of material ductility. Another quantity that can be obtained from Figure 11.2 is the toughness of the material, represented by the area under the curve, associated with plastic deformation. This area is the full integrated area under the curve minus the shaded region on the right, which represents the elastic energy stored in the material and which is recovered upon fracture. The more the material can deform before fracturing, the tougher the material is. Notice that in a tensile test the true-stress/true-strain curve is always higher than the engineering stress- engineering strain curve. Although at small strains the two curves are similar, as the strains become larger they start to deviate from each other, with the truestress/true-strain curve not showing the apparent stress decrease shown in the engineering curve. 11.2.2. Biaxial Testing: Tube Burst Although the uniaxial tensile test described above is experimentally convenient and its theoretical interpretation is made straightforward by the fact that there is only one nonzero component of the stress tensor (the normal stress in the direction of the stress, σ xx ), the relevant states of stress for the loading of fuel cladding loading and other core materials are often different and exhibit greater degree of stress triaxiality. For example the stress-state of fuel-element cladding loaded internally by fission-gas pressure more closely resembles that in a long thin-walled cylindrical tube closed at both ends and pressurized by a gas. This state of stress is a biaxial state of stress, in which the cladding wall is simultaneously stressed in the axial and hoop directions. Such a state of stress is best studied by a burst test. Creep rupture testing of unirradiated and irradiated steel tubing has also been performed by pressurizing closed tubing with an inert gas. According to elasticity theory (see Chapter 6) the stresses in such a condition are pr σ θ = t (11.9) Light Water Reactor Materials Donald Olander and Arthur Motta 6

2 pr σ z = t (11.10) σ r = 0 (11.11) where p is the internal pressure, R is the tube radius and t is the tube thickness. To analyze the plastic deformation of the material and to derive a yield condition in such multiaxial stress configurations, it is necessary to define an equivalent stress. 11.2.3. The Von Mises Equivalent Stress By definition, a specimen subjected to uniaxial loading yields when the stress reaches σ y. When a specimen is subjected to multiaxial loading, a yield criterion based on an equivalent stress must be developed. Such a criterion has been developed by von Mises, and is derived here for the case of the tube burst test. The criterion is based on the principal stresses in a material. [1] Von Mises general yield criterion is based on the difference between the total energy density under a multiaxial state of stress E el and the energy density resulting in the material when subjected to the mean of the three principal stresses, E el. When this * difference reaches a critical value E, i.e., when el * Eel Eel > Eel (11.12) the material will undergo plastic deformation. The mean of the three principal stresses is given by 1 σ = ( σ1+ σ2 + σ3 ) (11.13) 3. The elastic strain energy E el as a result of applied stresses is given by equation 6.31: E el 1 2 2 2 ν 1 2 2 2 = ( σ1 + σ2 + σ3) ( σ1σ2 + σ1σ3+ σ2σ3) + ( σ12 + σ13+ σ23) (11.14) 2E E 2G which for a solid under hydrostatic stress (no shear stresses) simplifies to E el = 1 2 2 2 ν ( σ1 σ2 σ3) ( σσ 1 2 σσ 1 3 σσ 2 3 2E + + E + + ) (11.15) Grouping the elements together yields the value of the homogeneous (average) elastic energy density E el 1 2 ν ( σ σ σ 6E 2 = 1+ 2 + 3) (11.16) Light Water Reactor Materials Donald Olander and Arthur Motta 7

Subtracting equation (11.16) from (11.15) we obtain 1+ ν 2 2 2 Eel Eel = ( σ1 σ2) + ( σ2 σ3) + ( σ1 σ3) 6E (11.17) By setting equation (11.17) equal to the case where there is only one applied stress allows to obtain the equivalent stress σ eq 1+ ν 2 1+ ν 2 2 2 2 σeq = ( σ1 σ2) + ( σ2 σ3) + ( σ1 σ3) 6E 6E (11.18) from where we obtain the von Mises equivalent stress 1 σeq = σ σ + σ σ + σ σ 2 2 2 2 1/2 ( 1 2) ( 2 3) ( 1 3) (11.19) Since the equivalent stress derived above defines yielding in a multiaxial stress state then equation (11.19) defines a yield locus, which is a surface in three-dimensional space. The yield locus for the von Mises condition for the case in which σ 3 =0 (a plane stress condition) is illustrated in Figure 11.3 Figure 11.3: Yield locus generated by the von Mises equivalent stress criterion [2]. The dotted lines indicate the Tresca yield criterion. For an isotropic material, a uniaxial tensile test results in a strain along the principal stress axis equal to twice that of the transverse directions and thus in the elastic regime σ1 = 2σ2 = 2σ 3 (11.20) For that case when σ 1=σ y, it is easy to show that σ eq=σ y Light Water Reactor Materials Donald Olander and Arthur Motta 8

We note here that the yield locus shown in Fig.11.3 and it is symmetric with respect to the deformation direction, i.e. it is the same in tension or compression. This does not always have to be the case and materials such as zirconium are anisotropic and show higher yield stresses along particular directions of deformation. Clearly, the yield locus depends on the mechanisms of deformation at work and it is possible that the operation of other plastic deformation mechanisms could intervene to cut the yield locus (i.e. achieve plastic deformation at a lower stress than the von Mises criterion). Figure 11.4 shows the yield locus for textured hexagonal close packed sheet material, which is seen not to obey equation (11.19). Shown in Figure 11.4 is also a line showing that under compressive states of stress twinning is favored over slip. Figure 11.4 Yield Loci for slip and twinning in textured hexagonal close packed sheet, after [3] 11.2.4. Hardness Testing Hardness testing is a comparatively easy means to obtain information about the deformation behavior of a material. In hardness testing, an indenter which can be of various shapes, but typically a ball or a pyramid- is applied with a certain force to a material. The indenter is much harder than the material to be tested, so that only the material tested plastically deforms. For a given force, the indenter will penetrate a certain depth into the material. The depth of indenter penetration can be measured microscopically by the area of indentation. Naturally, the depth of penetration and indentation area will be significantly larger in softer materials. As the indenter penetrates the material, plastic flow takes place under the indenter. The stress required for such flow to occur is the flow stress (defined above as the stress required to sustain deformation at a given strain). The presence of the elastic region nearby the indent causes a constraint on the deformation with the result that a much larger pressure than the yield stress is needed to cause plastic flow. For a circular indenter the relationship between the two is [2]. Light Water Reactor Materials Donald Olander and Arthur Motta 9

4P p= (2.5 to3.0) σ 2 y (11.21) π d where p is the pressure on the indenter, P is the applied load, d is the indent diameter. Thus σ y can be estimated from the test above. Several types of hardness testing exist, depending on indenter shape and whether indent area or depth is measured. Typical examples includebrinell, Meyer, Vickers and Rockwell. It is also possible to perform this test in a microscale, for example by applying an indenter that is smaller than the dimensions of one grain. 11.2.5. Impact Testing Another type of test that is of great relevance to the nuclear industry is the Charpy Impact test. Impact tests can be used to measure the fracture resistance of the material and can be at least qualitatively related to the toughness of the material. The Charpy test has been widely used in the nuclear industry for assessing the degree of embrittlement suffered by the pressure vessel after exposure to neutron irradiation. In the Charpy test a large hammer swung from a pendulum is released from an initial height h 1 towards a specimen in which a v-notch groove has been fabricated (Figure 11.5). The high speed of the hammer and the presence of the notch cause the material to be deformed at a high strain rate and in a triaxial state of stress, both of which favor fracture. The kinetic energy of the hammer is sufficient to break the sample. Subsequent to specimen fracture the hammer continues along its arc, rising to a height h 2, lower than h 1. The energy absorbed by the specimen during fracture is given by Eabs = mg( h1 h 2) (11.22) where m is the mass of the hammer is g is the gravitational constant. A sample made from a brittle material (such as glass) absorbs little energy as it breaks, while a ductile sample absorbs more energy before fracturing. The test is repeated at several temperatures. For most materials the energy absorbed shows an abrupt decrease at a certain temperature such that a ductile-to-brittle transition temperature (DBTT) can be defined. Light Water Reactor Materials Donald Olander and Arthur Motta 10

h 2 Figure 11.5. Schematic geometry of a Charpy v-notch Impact test. This is illustrated in Figure 11.6, which shows schematically the energy absorbed during Charpy Impact testing for a steel as a function of testing temperature (in an actual Charpy test the specimen would be positioned horizontally, shown vertically for clarity). The absorbed energy is relatively high at elevated temperature, undergoes a decrease in a relatively narrow temperature regime as the temperature is decreased, then stays at a lower level thereafter. The energies absorbed during fracture at high and low temperature are called the upper-shelf energy and lower-shelf energy, respectively. The midpoint of the transition from the upper shelf to the lower shelf energy is the ductile- to-brittle transition temperature, defined specifically in the case of pressure vessel steels to be the temperature at which the energy absorbed is 41 J. Analyses of the fracture surfaces from testing at low temperature and at high temperature reveals differences in fracture behavior: at low temperature the fracture surface is flat and shiny, indicating brittle failure by cleavage, while at high temperature the fracture surface is dull indicating ductile fracture. The degree of ductile failure can be evaluated by the fraction of the cross-sectional area that shows ductile behavior. Especially for low-tomoderate strength steels the Charpy test is a good indicator of the material toughness, as the fracture toughness often exhibits the same behavior as the Charpy energy, and therefore this is a useful test for assessing material ductility [2]. Light Water Reactor Materials Donald Olander and Arthur Motta 11

E abs (J) upper-shelf energy Brittle 41 J lower-shelf energy Ductile (DBTT) Testing temperature Figure 11.6 Schematic representation of energy absorbed during Charpy test as a function of testing temperature 11.3 Microstructural Aspects of Deformation 11.3.1. Critical Resolved Shear Stress and the Yield Stress As discussed in Chapter 7, plastic deformation in real crystals is effected by the movement of dislocations. As discussed in section 7.3 a slip system is activated when the Schmid factor for this system causes the resolved shear stress to be bigger than the critical value for slip. For a single crystal this predicts the onset of slip and thus the yield point. However, real materials are polycrystalline and contain grains of different orientations, more or less favorably oriented to the loading direction for slip to occur. As a result, some grains will deform more than others so that the apparently homogeneous repartition of strain actually results from grain-to grain inhomogeneous deformation. Because the deformation of a given crystallite (grain) is actually constrained by other neighboring grains, grain-to-grain incompatibilities arise during deformation which need to be accommodated by further plastic deformation. A simple illustration of such a process is shown in two dimensions in Figure 11.10. Consider two neighboring grains, which are oriented differently to the loading axis such that their shear deformation occurs in different directions. Light Water Reactor Materials Donald Olander and Arthur Motta 12

Figure 11.10. Schematic of role of various slip systems in ensuring deformation compatibility In such a situation, if the grains were free, the slip system operational on the left-hand side grain deforms the grain to the left, while the opposite happens on the grain on the right. To maintain the grain boundary together, it is clear that another slip system needs to be operational such that the deformation can be compensated in the other direction, thus keeping the grain boundary together. Thus, two slip systems are needed to maintain compatibility of deformation. In the case of a three-dimensional solid, the process above repeats itself in three different orthogonal directions, which means six independent slip systems need to be operational so that grain-to-grain compatibility of deformation can be achieved. However, the imposition of constant volume means one of those slip systems is not independent and the result is the conclusion that five slip systems are needed to induce deformation in a polycrystalline aggregate. To put it another way, we can assume that each grain deforms in the same way as the overall sample (i.e. it follows the same constitutive law). The strain tensor is written ε ε ε ε ε ε ε ε ε ε 11 12 13 = 21 22 23 31 32 33 (11.23) But from symmetry the off-diagonal terms ε ij =ε ji b and so there are six independent variables, and with the assumption that volume is conserved, then five independent variables remain. This means again that five independent slip systems need to be active in a given crystal in order for plastic deformation to occur. Light Water Reactor Materials Donald Olander and Arthur Motta 13

For a given load direction, the resolved shear stress τ RSS acting on a particular slip direction inclined by λ with respect to the loading axis and contained in a plane inclined by φ from the loading axis is F cosλ cosφ = τrss (11.24) A where F is the applied force and A the cross sectional area. When this stress is enough to cause the slip system under consideration to be operative, this is equal to the critical resolved shear stress. The factor cosλ cosφ is the Schmid factor for the slip system in question. Thus the yield stress in a polycrystalline material is related to the fifth largest Schmid factor for the active slip systems in the material. Because the fifth highest Schmid factor is much lower than the first, the necessary stress for plastic deformation is much higher than the critical resolved shear stress for a single crystal. This means that the yield stress for a polycrystal is larger than that of a single crystal by the Taylor factor m σ y = mτ (11.25) CRS S where m is approximately equal to the fifth largest Schmid factor in the sample. Because several grains are involved in ensuring deformation compatibility, the average Taylor factor m is used to describe yield in polycrystalline materials. The exact value of m depends on many factors, including the availability of slip systems, the ease of switching from one slip system to another, the texture of the material and the avialbility of other deformation mechanisms such as twinning. simulations conducted in both fcc and bcc materials show that the Taylor factor is 3-3.1 [4]. As mentioned above, although cubic materials normally can have five independent active slip systems, due to the paucity of slip systems, this is not true for hcp crystals, and plastic deformation is achieved with the help of twinning (see chapter 17). 11.3.2. Dislocations as cause of work hardening Once the material starts to yield, dislocation tangling will end up generating more dislocations by processes analogous to the Frank-Read sources described in Chapter 7. The increase in dislocation density with increasing deformation causes the stress required to cause further deformation to increase according to equation 7.36. This is one the principal causes of work hardening discussed in section 11.2 above. 11.3.3. Yield Strength Increase Mechanisms Several mechanisms exist for increasing the yield strength of a material and all of those involve hindering dislocations motion by either creating obstacles that force dislocations to bend around it, requiring more energy or by altering the dislocation itself such that its motion is impeded. Light Water Reactor Materials Donald Olander and Arthur Motta 14

As shown in equation 7.36, the critical stress to induce motion of dislocations through a field of obstacles is the Orowan Stress 2Gb σ c = (11.26) l where G is the shear modulus, b the burgers vector, and l is the spacing between obstacles. Any mechanism that creates obstacles that impede dislocations motion will increase in the stress necessary for dislocation motion and thus increase the yield stress. Equation (11.26) is valid for a uniformly distributed array of obstacles. For a random array of obstacles, the correct value has been shown [5] to be where α is a constant between 0.8 and 0.9. 2Gb σc = α (11.27) l In the same fashion as work hardening cold working increases the dislocations density which then represent an addition of barriers to the motion of other dislocations. A similar effect is caused by grain refinement. As the grain size becomes smaller, the distance between grain boundaries decreases, and since these also impede dislocation motion, the stress required for deformation correspondingly increases. Alloy hardening is also possible by the addition of foreign solute atoms, either in solid solution or as precipitate particles. In both cases the mechanism is related to slowing down dislocation motion. In the case of solid solution hardening impurities present in the bulk of the crystal migrate to the dislocation core. If the impurity atom size is substantially different from that of the matrix atoms, such segregation to the dislocations core will cause the dislocations to be semi- anchored by the atoms segregated. In the case of precipitation hardening, the added elements come out of solid solution and become obstacles to dislocations glide. The dislocations either have to cut through a precipitate particle (often harder than the matrix), or bend around it in the manner shown in the derivation of the Frank-Read stress. In either case glide is only possible at higher stresses. A review of strengthening mechanisms is available in other references [4, 5]. 11.4. Creep Deformation The materials tests described in the previous sections aim at determining the instantaneous material response, that is, when the structure of the material does not appreciably change during the test, as is the case at low homologous temperatures (T<0.3T m ) where T m is the melting temperature. It is possible to observe time-dependent deformation if the load is applied over a long period of time, especially above 0.3 T m. One simple example is to think of the deformation of a spring. If it is deformed under load in the elastic regime (σ<σ y ) and quickly unloaded, it returns to its original shape. However if, instead, the loading is applied over a long period (e.g. many months), it is Light Water Reactor Materials Donald Olander and Arthur Motta 15

possible that the strain will not be recovered, i.e. there will be plastic deformation. This slow (time-dependent) plastic deformation at a stress lower than the yield stress is called creep. Whether creep or yielding is the limiting factor in engineering design depends on the temperature and on the stress. It is possible to generalize plastic deformation in a material by considering under which conditions of stress and temperature each deformation mechanisms is observed. This has been done graphically by Ashby and Frost [6]. At high enough stress, the theoretical shear strength of the material is reached and the materials fractures, without the intervention of lattice defects. This theoretical shear strength decreases slightly with temperature as the material becomes softer. However, as mentioned in Chapter 7, this theoretical shear strength is seldom reached in reality as dislocations can mediate plastic deformation effectively, at much lower stresses (on the order of 10-3 to 10-4 G). When the applied stress is smaller than that required for dislocation motion, slow strain rate deformation is possible by various creep mechanisms dfescribed in this section. Phenomenological Description of Creep Creep of a material subjected to a constant load is often observed to be divided into three stages ε I II III t Figure 11.7: Schematic description of creep strain versus time. Figure 11.7 shows schematically the creep strain as a function of time [7]. In Stage I the initially high strain rate decreases with further deformation. This is likely due to processes akin to strain hardening in which the creep deformation causes the material structure to change such that the strain rate is reduced. When the twin processes of thermally induced deformation recovery and deformation induced hardening match, the material reaches Stage II, or secondary creep, which is characterized by a constant strain rate. This stage may last the greater part of the material lifetime. Finally, Stage III or tertiary creep causes the strain rate to increase and leads to material failure. Light Water Reactor Materials Donald Olander and Arthur Motta 16

The life of the material under creep deformation is usually dominated by steady state creep, which is sensitive to both stess and temperature. As a result, the creep strain rate under steady state creep ε SS is often described empirically by a power law equation of the form Q k kt B SS = Aσ e ε (11.28) where A and k are material constants, σ is the applied stress and Q is the activation energy for the relevant creep process. Larson-Miller Plots Assuming the dominance of steady state creep, Larson-Miller plots are a convenient method of predicting and comparing the creep strengths of engineering alloys. Taking the logarithm of both sides in equation (11.28) we obtain Q ln ε = + f ( σ) (11.29) kt B k where f( σ ) = ln( Aσ ). If a linear relationship exists between ε and the time to failure t f (such as exists in secondary creep) then B ε = (11.30) t f where B is a constant and where C is a constant B Q f( σ ) ln = ln t f + C = t k T f B (11.31) From equation (11.31) we define the Larson-Miller parameter LM. Q LM = T[ln t f + C] = (11.32) k Here Q/k B is a material constant, and the terms on the left hand side of the equation are the parameters of the creep test. The LM parameter is a good indicator of a material s creep resistance. Different materials show different LM plots. Figure 11.8 shows a LM plot for various candidate materials for advanced reactors [8, 9]. B Light Water Reactor Materials Donald Olander and Arthur Motta 17

Figure 11.8: Larson-Miller plots for various ferritic-martensitic steels considered for advanced reactors.[9] Example: Use the LM plot to determine the maximum allowable stress for HT-9, when operated at a temperature of 700 K, so that failure does not occur before 10 4 h. The LM parameter under these conditions is LM=700(ln(10 4 )+30)x10-3 =27.8 so that the maximum stress allowed would be about 150 MPa. 11.4.1 Creep Mechanisms During creep, plastic deformation which was not possible during instantaneous loading is made possible by thermally activated processes that allow permanent deformation to occur over a period of time. Creep occurs by several mechanisms and because it depends on thermally activated processes normally dislocation climb and glide or other thermally activated dislocations motion process or diffusional flow of atoms -, creep is more prevalent at higher temperatures. These mechanisms normally involve the thermally assisted motion of dislocations or a flux of point defects as briefly reviewed here. Thermally enhanced glide Light Water Reactor Materials Donald Olander and Arthur Motta 18

Computer simulations have shown that a dislocation in a crystal is constantly vibrating around a lattice position. If the critical shear stress for the dislocation is almost reached, thermal energy can provide the additional amount that would get the dislocations over the hump and allow it to glide. Thermally induced dislocation glide (Climb and glide mechanism) If a dislocation is stopped by an obstacle in the microstructure which exerts a force on it greater than the force from the applied stress, it is possible for glide to continue if the dislocation can absorb a net flow of point defects so that it can climb to another plane where the force exerted by the obstacle is less and glide is possible. This is then normally controlled by the diffusion of point defects to the dislocations allowing it to climb. Stress Induced Diffusional Flow Difference (Nabarro-Herring Creep): When a stress is applied to a solid, the equilibrium concentration of point defects is altered (see Chapter 4), such that, relative to the unstressed value, the vacancy concentration is higher in the presence of a tensile stress and lower in the presence of a compressive stress. This causes the vacancy concentration to be higher in the tensile regions than in the compressive region. Therefore a flux of atoms occurs to the grain boundary under tension and preferential grain growth occurs causing a deformation in the direction of the applied load. Coble Creep: In the same manner as in Nabarro-Herring creep, Coble creep occurs because of a vacancy gradient induced by the sample loading. However in this case mass transport is achieved through the grain boundaries of the material. This causes also an overall grain extension and deformation in the direction of loading. Thermally Induced Grain Boundary Sliding The strength of the material near the grain boundaries decreases with increasing temperature so that grain boundary sliding and consequent formation of cracks can occur at high temperature. Deformation Mechanism Maps The regions of applicability of these various deformation mechanisms can be made visually explicit in regions of applicability in a single map. This has been done by Ashby and co-workers [10] for various materials. One example is shown in Figure 11.9. The region at σ> σ y, deformation occurs by general yielding and at even higher stresses we reach the ideal strength of the material. At σ < σ y plasticity can occur by several creep mechanisms, activated at different temperatures. Light Water Reactor Materials Donald Olander and Arthur Motta 19

1.E+00 Normalized Shear Stress 1.E-01 1.E-02 1.E-03 1.E-04 1.E-05 1.E-06 1.E-07 Diffusional Creep by Grain Bondary Diffusion (Coble Creep) Theoretical Shear Strength Dislocation Glide Dislocation Core diffusion Power Law (dislocation Climb and Glide) creep Lattice Diffusion Diffusional Creep by Bulk Diffusion (Nabarro- Herring Creep) 1.E-08 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 T/T m Figure 11.9: schematic deformation mechanism map based on pure nickel with a grain size of 1 mm [10] Other creep mechanisms are operational under irradiation as reviewed in Chapter 27. 11.5. Material Fracture At some point during material deformation the process undergoes a qualitative change which leads the material to fracture, which is discussed in this section.. Fracture can occur as a bulk deformation and rupture or as the result of the propagation of a flaw, such as a crack. In bulk fracture, after the onset of plastic deformation, at some point the material can undergo a loss of load carrying capacity as a result of necking in the gauge section (the ultimate tensile strength illustrated in Figure 11.2). Ductile fracture is usually a result of an accumulation of damage in the material in the form of voids which eventually lead to failure. Brittle failure can occur either by cleavage, in which a microcrack propagates along favorable crystallographic directions ior by intergranular fracture, normally resulting from grain boundary decohesion because of grain boundary segregation of tramp elements in the alloy or by cracking of particles such as carbides Light Water Reactor Materials Donald Olander and Arthur Motta 20

which precipitate at the grain boundary. These mechanisms are reviewed in the following 11.5.1. Ductile Failure by Diffuse Necking: The Considère Criterion As the material starts to deform plastically, the strain is uniform (see Figure 11.2), until the UTS where necking starts to appear. Once a small deviation from uniform strain appears in the gauge section, such as shown in Figure 11.1, a loss of load carrying capacity develops, as additional strain causes the stress in the neck to increase, which causes further deformation and so on. However, as the material deforms, it also work hardens, which raises the flow stress. Thus a criterion can be developed whereby the material becomes unstable when it starts to deform faster than it can strain harden [11]. This is called the Considère criterion, as derived below. The first thing to note is that the maximum load occurs at the instability and thus and so P = σ. A => dp = Adσ + σda (11.33) da d d d A = σ σ ε σ σ = dε = (11.34) Consider further that when a rod such as in Figure 11.1, subjected to a load P, initial length l and initial cross sectional area A, the material is deformed plastically to a length l+dl which results in a diminution of area to A. This causes the material to strain harden. If the material deforms according to equation (11.8), σ The strain hardening caused by an incremental strain dε is dσ dε n = Kε (11.35) Knε nσ ε n 1 = = (11.36) At the same time the area is reduced by an amount given by conservation of volume This implies that Al. = A'( l+ dl) (11.37) da dl = = dε (11.38) A l So from equations (11.36) and (11.34) Error! Reference source not found., we obtain a criterion for the development of a plastic instability ε = n (11.39) Light Water Reactor Materials Donald Olander and Arthur Motta 21

This equation, the Considère criterion[12], states that at the onset of a necking instability the strain is equal to the strain-hardening exponent n. 11.5.2. Ductile Fracture Once necking occurs, the actual fracture of, the material normally occurs by a process of damage accumulation. Due to the presence of grain boundaries, and specially hard second phase particles or inclusions (such as intermetallic precipitates, carbides, sulfides or hydrides), deformation incompatibilities develop within the material as it plastically deforms. During uniform deformation these deformation incompatibilities are accounted for by other slip systems that allow grain-to-grain interfaces to remain whole as shown schematically in Figure 11.10. Light Water Reactor Materials Donald Olander and Arthur Motta 22

(a) (b) (c) (d) Figure 11.11 Schematic representation of void formation and coalescence in ductile failure. As a result of dislocation pileups and interactions, or because enough damage has accumulated nearby the second phase particles and interfaces, this process of deformation accommodation may start to break down. This causes damage in the form of voids to nucleate at grain boundaries or at precipitate-matrix interfaces. It is also possible that the second phase itself may crack. In any case, the result is the formation of voids in the highly strained areas as illustrated in Fig.11.11. Upon further straining of the material such voids grow and eventually coalesce. To correctly predict failure in this case it is necessary then to have both a void growth law and a void coalescence condition. Many such models have been developed [13] but are beyond the scope of this book. As a result of the formation of voids before fracture (b) and their coalescence (c), the fracture surface will normally exhibit a cup-and-cone shape (d). Also the fracture surface will exhibit ductile dimples corresponding to the voids, as shown in Fig.11.12 picture of dimples Fig.11.12 Light Water Reactor Materials Donald Olander and Arthur Motta 23

11.5.3. Fracture due to Crack Growth Brittle fracture occurs without significant overall plasticity of the material. This can occur by the formation of a crack and its propagation. The presence of a crack or a flaw in a material can severely affect its strength and fracture. The problem then is to calculate what controls crack propagation, and this is covered in the following section. Griffith Fracture Theory Figure 11.12 shows the process of formation of a crack of unit length in a solid under load. If a section of this material is cut, introducing a crack of depth c, a certain amount of elastic energy is released. This represents a driving force for crack propagation, since elastic energy stored is thus released. In counterbalance to this driving force, the propagation of a crack, (thus causing the material to fracture) requires the formation of two surfaces. The energy required to form these new surfaces is a large contributor to the fracture energy. As a result, the condition for unstable crack propagation can be derived from a comparison of energy for new surface creation with the elastic energy released upon crack advance. This was the theory derived by Griffith in 1921. However, except in very brittle materials, this theory tends to underpredict the work required to advance a crack. This is because in materials which exhibit some ductility it is necessary to induce some plastic deformation prior to fracture. This occurs in a region ahead of the crack called the plastic zone. The Griffith theory of fracture with a correction for the plastic zone is derived in the following. Consider a solid under load in which a crack is introduced, as shown in Figure 11.12. As mentioned above this will release elastic energy and create surface energy. Light Water Reactor Materials Donald Olander and Arthur Motta 24

σ plastic zone σ c Figure 11.12 A crack advancing in a solid under load The energy cost to creating two new surfaces and a plastic zone ahead of a crack of depth c is given by E = (2 γ + γ ) c= γ * c (11.40) surf surf plastic where γ surf is the energy to create a new surface and γ plastic is the energy necessary to create a plastic zone at the tip of the crack. The elastic strain energy released per unit volume is E el 2 2 Eε σ = = (11.41) 2 2E We assume all elastic energy is contained in a half-cylinder whose cross section is the semi-circle radius c around the crack. The total elastic energy released by the creation of the crack length c is then: E el 2 2 π c σ = 2 2E (11.42) Light Water Reactor Materials Donald Olander and Arthur Motta 25

The presence of the crack relieves the strain energy in the shaded region. Thus a solid o with initial elastic energy would have and E el E el 2 2 o πσ c = Eel (11.43) 4E Etot = Esurf + Eel (11.44) E tot πσ c 4E 2 2 * = γ c (11.45) To find the condition for crack propagation we set the derivative of the energy relative to crack size equal to zero and rearrange de tot 0 dc = (11.46) detot dc πσ c 4E 2 * = γ 2 = 0 (11.47) * * 2Eγ 2Eγ σ = => ( σ c) crit = = K IC πc π (11.48) The right hand side of equation (11.48) contains mostly material properties (the size of the plastic zone in γ* is also a function of the stress intensity). The level of this parameter that causes failure under a specific set of conditions (plane-strain deformation in mode I) is called the fracture toughness of the material K IC. In actual tests, a geometric correction is needed to take into account the increased stress intensity in the vicinity of a crack. K = Yσ c (11.49) IC where Y is a constant between 1 and 2, depending on the geometry of the crack. [7] Various closed forms of Y for many different geometries are available [7]. If the maximum allowable flaw is substituted into equation (11.50), then the critical stress intensity factor is obtained when the design stress is obtained. Equations (11.48) and (11.49) indicate that the higher the elastic modulus, the surface energies and the larger the plastic zone, the higher the fracture toughness.. Light Water Reactor Materials Donald Olander and Arthur Motta 26

where Κ Ι is the stress at the crack tip, and ρ is the crack curvature [7]. The K I refers to the fact that this crack is opening under mode I loading. Mode I loading is the most common mode of loading for engineering materials. Thus in general K = Y σ c (11.50) I where Y is a geometrical factor depending on the type of crack and whose value is normally between 1 and 2. In an inert environment with KI < KIC, cracks are indefinitely stable. In an aggressive environment, (Chapter 25) failure can occur by slow crack growth, i.e. if KIC > KI > KISCC (11.51) K IC (critical stress intensity or fracture toughness) is measured in laboratory with carefully prepared cracks in fracture toughness specimens. Light Water Reactor Materials Donald Olander and Arthur Motta 27

Problems 11.1 A metal rod of initial length l 0 and cross sectional area A 0 is tested uniaxially by fixing one end of the rod to a stationary point and loading the other with a spring of force constant K. The change in length of the rod, l-l 0, is related to the change in length of the spring, Δx, by: Δx = Δx Y - a(l-l 0 ) 2 + b(l-l 0 ) where the first term on the right is the spring extension at the yield point of the rod and a and b are empirical constants. (a) Express the true stress σ x in terms of the true (logarithmic) strain ε x and the parameters K, A 0, l 0, a, b, and Δx Y. Do not assume a functional relation between these two variables(e.g., Eq(7.3)). Neglect elastic strain. (b) Use the criterion for plastic instability to determine the true strain at the UTS. 11.2 The true stress-true strain curve in the work-hardening region can be written in dimensionless form: σ σy σ Y ( ε ε ) n = A Y where the subscript Y denotes the yield point, A is a dimensionless constant, and n is the work-hardening coefficient. Define X = (σ UTS - σ Y )/σ Y, where UTS denotes the ultimate tensile stress. (a) What is the equation relating X to A and n? (b) solve this equation for A = 23 and n = 0.7. (c) show that X increases with n at constant A 11.3 For a particular metal, the true stress and the true strain obey the following plastic constitutive law: σ x = σ Y exp(ε x - ε Y ). Prove that this material has no work-hardening region. Light Water Reactor Materials Donald Olander and Arthur Motta 28

11.4 A closed cylindrical tube is pressurized to exceed yielding by 20%. The constitutive equation for the material of the tube in the plastic region is: σ x - σ Y = Cε x n (elastic strain at the yield point is negligible). (a) What are the stresses in the tube wall in terms of the uniaxial yield stress of the material? (b) The plastic strength of the material is defined by: S = (σ * - σ Y )/ε *, where σ * and ε * are the equivalent stress and strain, respectively. Determine S. (c) What are the radial, azimuthal, and axial plastic strain components in the tube wall? 11.5 A spherical shell of radius R and wall thickness δ is pressurized internally to a pressure p. (a) What stress components are nonzero and what are their values(use a force balance as in Sect. 5 for the pressurized cylinder) (b) If the uniaxial yield stress is σ Y, at what internal pressure will the shell yield? 11.6 The creep law of a particular steel is found in uniaxial tests to obey: ε = T 6.4 4 161σ exp( 3 10 / ) where the strain rate is in hr -1, the stress is in MPa, and the temperature in Kelvin. A component fabricated from this steel operates at 330 o C. The principal normal stress components at a particular location are σ x = 380 MPa, σ y = 220 MPa and σ z = 180 MPa. How long is required for the creep strain (in any direction) to reach 2%? 11.7 A weight of mass m is suspended from a closed thin-wall tube of radius R and wall thickness δ with internal pressure p. The yield strength of the tube material is σ Y. Allowing for a 25% safety factor, what is the maximum gas pressure in the tube? Light Water Reactor Materials Donald Olander and Arthur Motta 29

11.8 A tensile test specimen is initially of length l 0 and cross sectional area A 0. The specimen is subjected to varying loads P and the deformation Δl measured at each. The data are correlated by the empirical equation: P = c + bδl - a(δl) 2 where a, b and c are positive constants. (a) What are the engineering yield and ultimate stresses of this material? (b) Show how the true stress-true strain curve can be deduced from the above equation. Neglect elastic deformations and assume all deformation to be plastic. 7.9 When tested uniaxially at a particular temperature, the creep rupture behavior of a 1/m batch of metal tubing obeys the equation: σ x = C( tr), where C and m are material properties. This equation is Eq (7.7) of the notes with the temperature dependence included in the constant. (a) When this tubing is sealed at both ends and pressurized at constant pressure p, what is the rupture time? The tube diameter is D and its wall thickness is δ. (b) If the internal pressure varies linearly with time(p=kt), what is the time to failure? 11.10 A rod just fits between two rigid plates when its temperature is T o. The rod is heated quickly and uniformly to a temperature T but the restraining plates do not move. heated rod Rigid plates (a) What is the initial axial stress in the rod when the temperature reaches T? Hint: Eq (5.14c) can be viewed as the total strain equal to the sum of the elastic and thermal strains, or ε tot = ε el + ε th. (b) With time, the stress computed in part (a) is reduced because of creep (this is called creep relaxation). Since the thermal strain does not change with time, the total strain rate is the sum of the time rate of change of the elastic strain plus the creep strain rate, or ε tot = ε el + ε. At the final temperature T, the creep rate is cr n ε cr = Bσ, where the constant B includes the temperature effect. Derive the equation giving the variation of the axial stress in the rod with time. Hint: in both parts (a) and (b) keep in mind that the total length of the rod remains constant. Light Water Reactor Materials Donald Olander and Arthur Motta 30