0
2006 NADAI AWARD PAPER

A Comprehensive Theory of Yielding and Failure for Isotropic Materials OPEN ACCESS

[+] Author and Article Information
Richard M. Christensen

Department of Aeronautics and Astronautics, Stanford University, Stanford, CA 94305 and Chemistry, Materials and Life Sciences Department,Lawrence Livermore National Laboratorychristensen@stanford.edu

J. Eng. Mater. Technol 129(2), 173-181 (Apr 01, 2007) (9 pages) doi:10.1115/1.2712847 History:

A theory of yielding and failure for homogeneous and isotropic materials is given. The theory is calibrated by two independent, measurable properties and from those it predicts possible failure for any given state of stress. It also differentiates between ductile yielding and brittle failure. The explicit ductile-brittle criterion depends not only upon the material specification through the two properties, but also and equally importantly depends upon the type of imposed stress state. The Mises criterion is a special (limiting) case of the present theory. A close examination of this case shows that the Mises material idealization does not necessarily imply ductile behavior under all conditions, only under most conditions. When the first invariant of the yield/failure stress state is sufficiently large relative to the distortional part, brittle failure will be expected to occur. For general material types, it is shown that it is possible to have a state of spreading plastic flow, but as the elastic-plastic boundary advances, the conditions for yielding on it can change over to conditions for brittle failure because of the evolving stress state. The general theory is of a three-dimensional form and it applies to full density materials for which the yield/failure strength in uniaxial tension is less than or at most equal to the magnitude of that in uniaxial compression.

FIGURES IN THIS ARTICLE
<>

The failure of materials generates research at all length scales from the electronic state to the atomic scale to nano to micro and on to macro (macroscopic) scales. The resulting information is considerably enhanced when the effects at the various scales can be interrelated. Usually the intention is to produce a particular result at the macro scale with this behavior being controlled by the mechanisms operative at the smaller scales. However, it can be difficult to confidently and securely approach the intended macro scale when its characterization is so shrouded in doubt and uncertainty. It would be highly advantageous to have a complete and comprehensive account of failure at the macro scale, one which transcends the various materials classes. The alternative is to have descriptors that are unique to each class or subclass of materials but with great uncertainty as to the range and limits of validity of each characterization. This latter condition represents the current status.

The objective here is to present and then probe a reasonably complete macroscopic theory of yielding and failure for homogeneous and isotropic materials. Much of the formalism will be synthesized here from various publications that have recently appeared but of necessity have been given in somewhat fragmented and unrelated forms. Using this new formulation, new results will be found for the yielding, plastic flow and failure in several important problems or classes of problems. This communication completes the cycle of recent papers on the failure of materials mentioned above. In whatever direction future work in the field goes, the present work may help to stimulate further interest and related activity.

It is an evident irony of the history of mechanics that the many books written on the strength of materials basically had almost nothing to do with that subject. Such books were virtually confined to the linear range of elastic behavior with minimal or no attention to failure. This occurred because the understanding of failure as an organized discipline was nonexistent. This state continued until the advent of fracture mechanics, about which more will be said later. The sparse historical scene of successful research upon materials failure did have one major prominence and this account should begin by acknowledging the subtle but profound contribution of Coulomb (1). Mohr (2) put Coulomb’s failure result into a form allowing easy utility. The fact that the Coulomb-Mohr failure form does not successfully account for many of the physical effects does not detract from the efforts of either scientist. In the time frames of their separate works, their grasp of the problem was completely beyond compare. Later, the Mises criterion was given (3), but only as an adjunct to a special case of the Coulomb criterion, namely, the Tresca form. Both the Mises and Tresca criteria apply only to the yielding of very ductile metals. The Coulomb-Mohr form was intended to apply across the spectrum of materials types, as is the interest here. A history of strength and failure treatments has been given by Paul (4), and a brief historical summary by Christensen (5).

At the most elementary level it is sometimes said that ceramics are brittle, many but not all metals are ductile, and some types of polymers are ductile and others are brittle. While there is a degree of truth in this assertion, it can be extremely misleading, or even worse. An example will be given later wherein a material commonly considered as being completely ductile when placed in a particularly important special state of stress fails in a brittle manner. Relative to failure, all materials can behave either in a ductile or a brittle manner depending upon the state of stress that they are under and other environmental influences.

The two terms, yielding and failure, have imprecise definitions that usually allow a wide latitude of interpretation. This imprecision underlies an uncertain basis of operation. The obvious exception to this situation was the development of fracture mechanics. Fracture mechanics was one of the technical achievements of modern mechanics. A typical application of fracture mechanics involves determining the stability, under imposed stress, of imperfections and stress risers such as cracks and edge notches, holes, attachments, etc. In contrast, the means of applying fracture mechanics to the problem of the failure of homogeneous materials under uniform stress states has been far less clear. There are many opinions on this subject, but little substance beyond general statements. Thus, for the failure of homogeneous materials, the integration of fracture mechanics into a more general formalism has not been successfully accomplished in the past.

Failure criteria have usually been formulated in terms of stresses, but over the historical time span, failure criteria have occasionally been postulated in terms of strains. The view here is that trying to specify failure in terms of strains is inappropriate and internally inconsistent. Stress must be used in order to have compatibility with fracture mechanics in the brittle range and with dislocation mechanics in the ductile range. To not have union with these two anchor points of physical reality would be extremely serious. Stress, not strain, will be used here for these well grounded reasons.

In the modern era there have been many attempts to find criteria more general than just that for the perfectly ductile response or alternatively the fracture controlled response. The references cited above give many previous references to such works. A sampling of these efforts should include the following. Drucker and Prager (6) gave a two-parameter yield criterion of conical form in principle stress space. Paul (4) proposed a three-parameter pyramidal type yield surface. Wronski and Pick (7) applied Paul’s criteria to polymers. Raghava, Caddell, and Yeh (8) proposed a criterion similar to parts of the present forms, and applied it to polymers. Stassi (9) also discussed similar forms, but without application. Pae (10) applied a three-parameter criterion to polymers. Wilson (11) applied the Drucker-Prager criterion to metals. Jaeger and Cook (12) discussed many three or more parameter models for application to geological materials. None of these approaches possess the combined attributes of involving only a few adjustable parameters, preferably only two, along with the power and flexibility to recreate many different physical effects for many different types of materials. Also, none of these works deal with the essential problem of providing indicators (for any stress state of interest) that differentiate a materials capability for undergoing ductile flow as opposed to the undesirable outcome of brittle failure.

The present work will give specific meanings to the terms yielding and failure, ones that will offer a useful distinction between them. There are two technical keys to the following developments. These are: (i) the explicit integration of a fracture mechanism into the yielding versus failure formalism, and (ii) the derivation of an explicit criterion that determines whether a failure mode in any particular state of stress is expected to be of ductile or brittle nature. Not surprisingly, these developments (i) and (ii) are found to be interrelated, but still take separate forms in the final set of equations and conditions.

In the next section, the governing relations of this new theory will be given.

The references from which various aspects at this new theory of yielding and failure are collected are Christensen (5,13-17). The initial work in 1997 identified a nondimensional properties grouping that spanned the range from completely ductile to very brittle behaviors. The work in the year 2000 was a view of ductile versus brittle behavior, approached in a very mathematical formalism. A broader treatment in 2004 brought in an explicit fracture mechanism and involved comparison with experimental results for a wide variety of materials types. In 2005 an explicit criterion was derived for distinguishing expected brittle failure behaviors from those of plastic yielding, expressed in terms of the imposed stress state and a particular material characteristic. The first work in 2006 was an examination of plastic flow potentials, needed when the condition of ductile behavior exists. The second work in 2006 involved a detailed comparison with the Coulomb-Mohr and Drucker-Prager (6) theories. The collection of these individual works and particularly the present amalgam of all of them serve to form this comprehensive account of yielding and brittle failure.

The mathematical conditions to be given in this section for yielding and failure are primarily taken from the above references, but for the background details, the particular references should be consulted. A uniform notation and terminology will be adopted. In this section and the following sections some previously open issues will be clarified and closed, and some important special problems will be examined and some critical interpretations given.

The governing yield/failure function for isotropic materials was found by taking a polynomial expansion, through terms of second degree, of the invariants of the stress tensor. This procedure givesDisplay Formula

ασ̂ii+32(1+α)ŝijŝij1
(1)
where sij is the deviatoric stress tensor and σ̂ij is nondimensional stress with
σ̂ij=σijκ
where
α=CT1,0α<
κ=C
C and T are the uniaxial yield/failure stresses in compression and tension. Alpha, α, is the nondimensional properties grouping that affords special advantages.

When the yield/failure in uniaxial compression is not to be used, α and κ will be shown later to take other forms. Relation 1 as written out in terms of components is given in the last section.

Relation 1 is part of the governing criteria, the other essential part is an explicit fracture criterion given byDisplay Formula

σ̂111+αifα1
(3)
or in dimensional formDisplay Formula
σ1TifTC2
(3′)
where σ1 is the largest principal stress. The relationships of 1,3 to some historical criteria as well as their physical interpretations will be given at the end of this section and in the following section.

The two criteria 1,3 in the range of α1 are both employed to find whichever of the two is the more limiting. Also, the two relations 1,3 are in some sense coupled through the appearance of α in both of them. For α1 relation 1 stands alone. It is best to view the variation of α on a log scale having the range from to . Then the value α=1 in 3 designates the center of the log scale range, logα=0. This center point is where the fracture criterion 3 begins to take effect. Compared with the yield/failure criterion 1, the fracture criterion 3 can be shown to only have an infinitesimal difference from 1 in the immediate vicinity of α=1, logα=0, but it becomes of increasing significance as α increases beyond this value.

The fracture restriction 3 naturally corresponds to brittle behavior, but under certain conditions relation 1 also can embody brittle failure characteristics. This is why the form in 1 is termed as the yield/failure function. This condition is shown schematically in Fig. 1 for biaxial stress states with α=2. The intersection of the yield function in 1 and the fracture condition in 3 establishes the division into ductile and brittle regions. It is seen in Fig. 1 that the mean normal stress part of the failure stress state locates the position of the ductile/brittle dividing line. This ductile-brittle delineation generalizes directly to triaxial stress conditions and gives the D-B criterion as

σ̂ii<2α1+αDuctile
σ̂ii>2α1+αBrittle
where σ̂ii is found from the stresses satisfying relation 1 as an equality and α is taken over its full range. That is, σ̂ii is the first stress invariant (mean normal stress) of the yield/failure stress state. The right-hand side of 4 specifies the material type through the value of α, while the left-hand side specifies the part of the yield/failure stress state that controls whether the behavior will be brittle or ductile.

The region near the dividing line between ductile and brittle behavior specified by 4 is in reality likely to be transitional rather than discontinuous. Furthermore, this region is more likely to be dominated by the brittle inclination than by the plastic flow tendency. As a separate matter, any time the fracture criterion 3 is controlling, that by definition is a brittle behavior. These general behaviors are similar to those modeled by Harlin and Willis (18), wherein they used macroscopic criteria with the involvement of triaxial stress conditions to distinguish the ductile from the brittle response of growth activated cracks.

Relations 1,2,3,4 comprise the entire yield/failure criteria. Only two properties are involved, the uniaxial tensile and compressive yield/failure values. Strain hardening could be involved in this methodology, but it is not included here in order to focus upon the essential but idealized aspects of the ductile versus brittle behavior. The yield stress in a ductile behavior is taken to be that at the point of major deviation between the previous elastic region and the following plastic flow region. The idealized behavior is as shown in Fig. 2. It is clear that in the context of the brittle and ductile behaviors of Fig. 2 stress alone cannot distinguish between the three conditions of: (i) brittle failure, (ii) ductile yield, or (iii) ductile failure. Thus, it is necessary to have the separately determined criterion to distinguish brittle from ductile response. The complete theory does not predict the strain at failure after ductile flow, only that such a flow state exists. In some cases the amount of plastic flow before failure may be very small and it may be somewhat difficult to distinguish these cases from brittle failure.

It is to be expected that the ductile-brittle criterion 4 would be controlled by the mean normal stress part of the total stress tensor at yield or failure. This is similar to and in fact related to the effect of temperature on the same ductile, brittle behaviors, see also Harlin and Willis (18). Temperature and pressure are the two most fundamental control variables for such effects. Again it is emphasized that in physical reality there would not be an abrupt dividing line between the two states, but rather a region or band of most rapid variation from ductile to brittle.

To complete this theory it is necessary to specify the governing characteristics of plastic flow when criteria 4 specifies that the ductile regime controls. Following Hill (19) the plastic flow potential is necessarily taken to be different from the yield/failure function, in 1, thereby leading to a nonassociative form. The decomposition of strain components into elastic and plastic parts then has the plastic strain increments given byDisplay Formula

ε̇ijp=λGσij
(5)
where G is the independent plastic flow potential. It was found (in the previous references) that for homogeneous materials G takes the formDisplay Formula
G=αA1+ασ̂kk+32ŝijŝij
(6)
where A is a constant, an order of magnitude less than one, and with considerable support for a value of about A115, from the data of Richmond and co-workers: Spitzig, Sober, and Richmond (20), Spitzig, Sober, and Richmond (21), and Spitzig and Richmond (22). The hypothetical value, A=1, would correspond to the associative form. The dilatational term in 6 is due to the mechanism by which materials in plastic flow create a small amount of voids and vacancies.

In most (but not all) situations and for most homogeneous and isotropic materials the first term in 6 is small compared to the second term because of the smallness of constant A and can be neglected leavingDisplay Formula

Gŝijŝij
(7)
which gives the plastic flow potential as being completely distortional in character. In 7 the 32 factor in 6 is absorbed into λ in 5.

This is a two-property theory calibrated by the uniaxial tensile and compressive yield/failure values. It should be noted that the uniaxial compressive test can be difficult to perform. See Lassila et al. (23) for a detailed account of the measures and care needed to produce a reliable test for uniaxial compression. An example of an alternative set of two independent yield/failure prescriptions is that of uniaxial tension and uniaxial tension under pressure. In this case it can be shown that α and κ are given byDisplay Formula

α=2Θ3Θα1
(8)
where
Θ=dTdpp=0
and then
κ=(1+α)Tp=0
where p is the pressure and T is the uniaxial tensile yield/failure stress superimposed upon p. The restriction in 8 to the range α1 can be removed, but with some additional involvement of fracture relation 3 in the range α1.

Finally, it should be noted that if the fracture criterion 3 were used as a stand-alone criterion it would correspond to the maximum normal stress criterion favored by Rankine and Lamé. That criterion did not succeed as a single, all encompassing criterion, but it certainly has a special and vital role to play here as a restriction tightly coordinated with the yield/failure function 1. With regard to macroscopically homogeneous materials, the fracture criterion 3 is here seen as effectively a mode I fracture event activated by whatever inhomogeneities control the behavior on the micro scale. Since this theory is only for macroscopically homogeneous and isotropic materials, the micro scale inhomogeneities possibly causing fracture must not disturb the macro scale conditions of homogeneity and isotropy. It also follows that any initial micro scale porosity must be small.

Cases represented by the limiting values of α will be discussed as well as some intermediate values of α. The case of α=0 leaves the yield/failure function in 1 as reducing to the Mises criterion. This α=0 case will be referred to as that of a Mises material. It is commonly accepted that a Mises material is the example of a perfectly ductile solid, this implying that there exists ductility under all conditions. However, the present results show that this notion is not correct. The ductile-brittle criterion 4 is not always satisfied as being ductile for α=0. In fact, it can be shown that criterion 1,4 reveal that in principle stress space the infinite cylindrical form of the Mises criterion has ductile and brittle regions as shown in Fig. 3. The case of uniaxial tension is on the ductile side of the D-B dividing plane in Fig. 3 as are most of the common stress states such as unequal biaxial stresses. But the special case of equal biaxial tension is right on the dividing line between the ductile and brittle regions. More will be said about this important case later. As seen from Fig. 3 any stress state at yield/failure with a sufficiently large tensile hydrostatic component will be brittle. An example is as follows: When a high intensity compressive dilatational wave reflects from a free surface it creates a large triaxial tensile stress state condition. This can cause a micro scale void nucleation as a precursor to the macro scale brittle failure condition of spallation, even for what are normally thought of as ductile materials.

For values of α>0 the yield/failure function 1 is a paraboloid in principal stress space, Fig. 4. Its axis makes equal angles with the three principal stress axes. For α>1 the fracture criterion 3 comes into effect. Right at α=1 the three planes prescribed by the fracture criterion 3 are just tangent to the yield/failure paraboloid. As α increases beyond α=1 the fracture criterion takes three slices out of the paraboloid. These are here called fracture cut offs. An illustration of these fracture cutoffs at α=2 is shown in Fig. 1 for biaxial stress, with this value of α being typical of cast iron. The paraboloid is divided into ductile and brittle regions by its intersection with the plane normal to its axis prescribed by 4. This ductile-brittle division must also deviate around the fracture cutoffs prescribed by 3.

The limiting case of α is the brittle limit with the tensile strength T being negligible compared with the compressive strength C. This limiting case is perhaps best viewed through the shear stress at yield/failure. From 1 the shear stress τ at α is very small compared with C. Now consider a state of shear stress τ superimposed upon pressure p. From 1 at α the yield/failure shear stress under pressure is given byDisplay Formula

τ=pC
(9)
The associated ductile-brittle criterion 4 gives the result
p>C3Ductile
p<C3Brittle
The fracture criterion 3 givesDisplay Formula
τ=p
(11)
which is by definition brittle.

Conditions 9,10,11 give the complete behavior as

τ=pforpCBrittle
andDisplay Formula
τ=pCforpCDuctile
(12)
Relations 12 are shown in Fig. 5. Again, at p=0 the material cannot sustain any shear stress. Sufficiently large pressure supports a ductile behavior in shear. The case just considered corresponds to the condition σii0. For the condition of positive mean normal stress, σii>0, as α, relation 1 can be used to show that it becomes justDisplay Formula
σiiT
(13)
This relation is more restrictive than the fracture condition 3 and from 4 it is of a brittle nature. The physical interpretation of 13 is that approaching this brittle limit of α any stress state with positive σii satisfying 13 as an equality causes the material to disassociate and disintegrate. Although it is important to verify reasonable and rational behavior in the limiting case of α, probably the practical maximum range for applicability of this theory would be for ceramics and glasses. These materials would have values of α from about 5 to 25.

Now, two general examples will be given, first an equal biaxial stress state and then unequal biaxial stresses in a 2:1 ratio. For equal biaxial stresses letDisplay Formula

σ11=σ22=σ
(14)
The yield/failure criterion 1 givesDisplay Formula
σ̂=1(1+α)(α±1+α+α2)
(15)
For the case of tension and at α=0 then 15 gives σ̂=1. With this result then it follows that the left-hand side of 4 is given by
σ̂ii=2
The right-hand side of the ductile-brittle criterion 4 at α=0 is simply 2. Thus at α=0 and in equal biaxial tension the yield/failure state is borderline between being ductile or brittle. For α>0 it is always brittle thus for all values of α
EqualBiaxialTensionBrittle

Using 15 in the case of compression with criterion 4 then gives for all values of α

EqualBiaxialCompressionDuctile
The significance of this result that equal biaxial tension is brittle is that this is the stress state existing in thin spherical pressure vessels. The failure is of the brittle-type, even though the material, such as aluminum, for example, is nominally considered as being ductile. Thus thin spherical pressure vessels can represent a considerable safety hazard no matter how ductile the composing material may be, in accordance with practical experience.

The example just considered suggests looking at the case of cylindrical pressure vessels. In the thin cylindrical portion of the pressure vessel the stress state is given by biaxial stresses with

σ11=σ
σ22=σ2
The yield failure criterion 1 givesDisplay Formula
σ̂=1(1+α)(α±α2+43α+43)
(17)
Using the stress states for tension and compression from 17 in the ductile-brittle criterion 4 gives
Tension,α<12Ductile
Tension,α>12Brittle
and
CompressionDuctile
for all values of α. In the tensile case, the D-B dividing line at α=12 corresponds to TC=23, which is right in the range of many polymers.

In both of these latter two examples symbolic of thin shell pressure vessels, the fracture criterion 3 is not as limiting as the results shown above from the yield/failure criterion 1. These examples are for the two different biaxial stress states. When interpreted for application to thin shell pressure vessels, only the tensile conditions would be meaningful. Instability would supercede the compressive yield results.

Now consider more complex three-dimensional problems with nonuniform stresses. The Lamé problems for thick cylinders and spheroids are the classical problems of this type. In particular a thick spherical annulus with internal or external pressure will be considered.

The internal pressurization of the spherical annulus is as shown in Fig. 6. In the case of an ideal Mises material the complete solution is given by Hill (19). The solution for the more general material model involving yield and or failure governed by 1,2,3,4,5,6,7 will be developed here. For sufficiently small pressures the entire region will be elastic, but at some pressure the material at the inner radius will commence yielding, and as pressure further increases the elastic plastic boundary at r=ρ, Fig. 6, will move outward.

The stress solution for the entire elastic deformation case is given by

σr=p[(br)31(ba)31]and
σθ=σϕ=p[12(br)3+1(ba)31]
where p is the internal pressure.

Now find the pressure at which yielding first commences at r=a. The yield/failure function 1 becomesDisplay Formula

α(σ̂r+2σ̂θ)+(1+α)(σ̂rσ̂θ)2=1
(19)
Combining 18,19 at r=a gives the pressure for initial yielding asDisplay Formula
p̂=2(1λ)3(1+α)(λα+1+α+λ2α2)
(20)
where
λ=(ab)3
and
p̂=pκ
The symbol λ used here should not be confused with that in 5.

Next, the question of whether this initial yielding is actually of ductile yielding or brittle failure must be answered. Using the stresses at r=a in the ductile-brittle criterion 4 gives

3p̂(ba)31<2α1+αDuctile
3p̂(ba)31>2α1+αBrittle
where p̂ is given by 20. Interest here will be confined to cases of α1, thus the fracture criterion 3 does not enter the considerations. Materials with α>1 would almost never be used in applications of internal pressurization since they would be too brittle. Combining 20,21 gives
(ab)3<1α21α+α2Ductile
(ab)3>1α21α+α2Brittle
For example for α=122 became
ab<0.794Ductile
ab>0.794Brittle
The ductile case will be considered here. If cases with α2 were to be considered they would be found to always be of the brittle failure type. For values of α between 1 and 2 the previous criterion can be used to find the values of ab for which the ductile or brittle condition applies.

For pressures greater than that in 20 the medium deforms into a plastic interior region and an elastic exterior region, with the division at r=ρ. For the elastic region the stresses are given by

σr=A[(br)31]
σθ=σϕ=A[12(br)3+1]
where A is a constant to be determined. At r=ρ the elastic stresses must satisfy the yield/failure criterion 19, which determines A asDisplay Formula
Aκ=2(ρb)33(1+α)[α(ρb)3+1+α+(ρb)6α2]
(25)
For later use it is necessary to have σ̂r at r=ρ which is found asDisplay Formula
σrr=ρ=2[1(ρb)3]3(1+α)[α(ρb)31+α+(ρb)6α2]
(26)

With the elastic region determined, now the stresses in the plastic region will be found, being careful to only do this under conditions of ductile behavior. The governing equilibrium equation isDisplay Formula

dσrdr+2(σrσθ)r=0
(27)
Use the yield function 19 to solve for σ̂θ and substitute that into 27 givingDisplay Formula
(1+α)dσ̂rdr+2r(αη)=0
(28)
whereDisplay Formula
η=1+α+α23α(1+α)σ̂r
(29)
This nonlinear differential equation can be integrated to give σ̂r as a function of r. Carrying out this process, and evaluating the constant of integration such that at r=ρ, σ̂r from 28,29 is continuous with the elastic region stress 26 then givesDisplay Formula
3αlnrρ=[η+αln(α+η)]+[η1+αln(α+η1)]
(30)
whereDisplay Formula
η1=ησ̂r=σ̂rr=ρ
(31)
The symbol σ̂rr=ρ in 31 is given by 26.

The boundary conditions at r=a isDisplay Formula

σ̂r=p̂atr=a
(32)
which when substituted into 30 givesDisplay Formula
[η2+αln(α+η2)]=3αlnρa+[η1+αln(α+η1)]
(33)
whereDisplay Formula
η2=ησ̂r=p̂
(34)

Now the full stress solution in the plastic region has been found, subject to the condition that ductile plastic flow occurs, rather than brittle failure. This will be illustrated in a particular example. Equation 33 provides the relationship between the applied internal pressure p̂ and the location of the elastic-plastic boundary, ρ. Once this is known, then the stress σ̂r follows from 29,30 as a function of r. The other stresses follow from 24,25.

Take a particular example of ab=13 and α=1 which is typical of some metals and polymers. The ductile-brittle criterion 22 applied to the elastic region at r=ρ gives

(ρb)3<12Ductileor
(ρb)<0.794Ductile
Thus plastic flow ceases when sufficient pressure has caused the elastic plastic boundary to reach the value in 35. Brittle failure follows thereafter.

For this particular example the solution given here determines that for first yield

p̂=0.442Firstyield
and the pressure at which ρ becomes that shown in 35 is
p̂=1.70Terminationofplasticflowinceptionofbrittlefailure
For comparison the solution for a Mises material, α=0, gives
p̂=0.642Firstyield
and
p̂=2.20Fullyplasticyield
The stress components in this solution example for α=1 follow directly from the previous results. The plastic strain increments can similarly be found by using the plastic flow potential 6 or the simplified form 7.

An important characteristic has emerged from this solution and example. It is that in elastic, plastic flow problems, as the flow progresses and the boundary between the two regions changes, so can the conditions existing at the elastic plastic boundary which switch over from that of plastic flow to that of brittle failure. Thus the conditions for plastic flow, at a particular configuration, can cease to exist and plastic flow is interrupted by brittle failure. Although this behavior was found from the present idealized theoretical formulation, physical reality can be expected to at least show some features or aspects of this behavior. This possibility of switching the failure mode and type follows from the dependence of the ductile-brittle criterion upon the prevailing state of stress as well as the material type. This complex aspect of physical behavior with regard to the moving elastic plastic boundary appears to have been unrecognized before now.

The corresponding and companion problem of a thick spherical annulus subjected to external pressure can be solved by the same general method as that just discussed. Yielding first occurs on the inner surface at r=a and the external pressure to cause this is given byDisplay Formula

p̂=23(1λ)(1+α)(α+1+α+α2)
(36)
where λ=(ab)3 as before. An examination of the ductile-brittle criterion 4 reveals that the further deformation is always that of ductile plastic flow for all values of α, in contrast to the complex situation just discussed for the internal pressurization problem, where plastic flow can be interrupted by brittle failure under certain conditions. An example with ab=13 and α=1, 36 gives the external pressure at first yielding as p̂=0.877, about twice the value required by the internal pressurization case.

Although the present solutions show that the resulting stresses take more complicated forms than in the simplified Mises material case, it is still completely practical to obtain such results. The formulation is well posed and amenable to numerical solution in complex problems.

The orientations of failure surfaces reveal specific patterns and signs of behavior. Certainly brittle material orientations are expected to be much different from ductile material orientations. This topic will now be examined with the objective to determine the failure surface orientations for uniaxial tension and compression as a function of the material type, specified by the value of α in 2. Failure surface orientations in this context have been studied before (15), but there was some uncertainty as to whether the yield/failure function or the plastic flow potential should be used in certain operations. This question can now be formally and finally answered.

The usual plastic increment under ductile flow conditions is already stated as 5, where G is the plastic flow potential. The failure that occurs at the end of a ductile flow process is not simply a continuation of the previous plastic flow. It is an inherently unstable process, and it is taken here to be controlled by the yield/failure function, rather than the plastic flow potential. Take the increment of the failure strain, ε̇ijf, to be given asDisplay Formula

ε̇ijf=Λf()σij
(37)
where f() is the yield/failure function on the left-hand side of 1. This should not be confused with using the associated flow rule because that terminology only applies to the stable ductile plastic flow. The form 37 is here being applied to the final failure process.

For uniaxial tension, 1,37 then give

ε̇11f=Λ(2+α)
ε̇22f=ε̇33f=Λ(1+α)
For uniaxial compression 1,37 give
ε̇11f=Λ(2+α)
ε̇22f=ε̇33f=Λ(1+2α)
Now take a rotated coordinate system xi such that 1 is in the plane of the failure surface and 2 is normal to it. Let angle ϕ be the angle between the axial direction, 1, and 1 in the failure surface. Then the failure strain increments in the rotated coordinates are given by
(ε̇11f)Λ=(2+α)cosϕ(1α)sin2ϕ
(ε̇22f)Λ=(2+α)sin2ϕ(1α)cos2ϕ
(ε̇12f)Λ=3sinϕcosϕ

Take the failure strain increment (ε̇11f) in the plane of the failure surface as being nonactive and vanishing while the other two increments in 40 give strains that undergo unlimited change in the failure process. Setting (ε̇11f)=0 then gives the orientation of the failure surface asDisplay Formula

tanϕ=2+α1α,α1,Tension
(41)
At α=1, relation 41 gives ϕ=90deg. For α>1 the fracture cutoff 3 comes into effect giving brittle failure which results inDisplay Formula
ϕ=90deg,α1,Tension
(42)

For the case of uniaxial compression the rotation of coordinates givesDisplay Formula

(ε̇11f)Λ=(2+α)cos2ϕ+(12α)sin2ϕ
(43)
Similarly to the tensile case, for this failure strain increment in the plane of the failure surface taken as inactive and vanishing while the other 2 strain components grow to unbounded values then gives the orientation of the compressive failure surface asDisplay Formula
tanϕ=2+α1+2αCompression
(44)
The failure surface orientations are shown in Fig. 7 as a function of α. At α=0 (a Mises material) the failure surface orientation is at the octahedral angle for both tension and compression.

The uniaxial tensile case shown in Fig. 7 has the ductile versus brittle regions as shown. This behavior coordinates perfectly with the ductile-brittle criterion 4 which shows the change over at α=1 for the case of uniaxial tension.

In the case of uniaxial compression the ductile-brittle criterion 4 is in the ductile range of behavior for all values of α except that it becomes borderline between ductile and brittle as α. It is interesting to determine the stress components on the failure surface at the orientations given by 44. Let σ be the normal stress and τ be the shear stress on the failure surface. Using 44 and the appropriate coordinate rotations for stress it is found that

σ̂=(2+α)3(1+α)and
τ̂=(1+2α)(2+α)3(1+α)

The two failure surface stresses in 45 for values at α=0, α=1, and α are given in Table 1. It must be remembered that σ̂=σκ and that κ would vary greatly over the range of the α’s. These nondimensional stresses have a special variation with α. The shear stress shown in Table 1 hardly varies over the full range of α while the magnitude of the normal stress diminishes with increasing α. The latter variation can be interpreted as the indication that failure surface orientation is such that the compressive normal stress on the failure surface tends toward small magnitudes for large values of α, but has larger magnitudes of σ̂ at the smaller values of α. In other words, the magnitude of the normal stress is more controlling at large α’s and less controlling at small α’s, which is intuitively expected. Also shown in Table 1 is the ratio τσ. Note that τσ has a type of antisymmetry relative to the central point at logα=0.

Finally it is perhaps obvious that the yield/failure function rather than the plastic flow potential was the correct quantity to use in finding the failure orientations. If the plastic flow potential 7 were used, the failure surface orientations would have been unvarying with respect to α, always at the octahedral angle and the same in compression as in tension.

The previous results will not be summarized here other than to broadly say that this is an approach to the difficult problem of characterizing yielding and failure for general materials. The central focus has been to distinguish brittle failure from ductile yielding type response. It has always been known that this discrimination must depend upon a specification of the material type but the present work also shows that this distinction crucially depends upon the type of stress state under consideration. For example, shear stress gives a different result than does uniaxial tension. When one brings in the full stress tensor, the number of possible combinations becomes boundless and a generalized approach becomes necessary. These resulting special cases display a great variety of different ductile-brittle transition circumstances, many of which have been examined here and shown to be of importance.

It is of some relevance to consider which single stress state, if any, best characterizes a simplified but realistic form of the ductile-brittle characteristics of materials. Uniaxial tension certainly is the logical choice. From the yield/failure criterion 1, for uniaxial tension it follows thatDisplay Formula

σ̂11T=11+α
(46)
Substituting this into the ductile-brittle criterion 4 givesDisplay Formula
α<1Ductile
(47)
orDisplay Formula
TC>12Ductile
(48)
For α>1, the response in this stress state is brittle. This gives the transition in failure type for uniaxial tension as being at the value α=1. This is the same value as for the inception of the fracture criterion in 3. The value of α=1 is at the center of the logα scale. For simple shear stress, criterion 4 directly shows that its transition is at α=2. Other stress states show widely varying values of α at which their ductile to brittle transitions occur.

At an application level, often no attempt is made to distinguish ductile and brittle responses, other than by intuitive guidelines. In this situation failure would be considered in an inclusive sense as encompassing both yielding and complete rupture. Then the previous results could simply be designated as generalized failure criteria. In this case everything needed is directly specified by conditions 1,2,3, rewritten here in terms of components asDisplay Formula

(1T1C)(σ11+σ22+σ33)+1TC{12[(σ11σ22)2+(σ22σ33)2+(σ33σ11)2]+3(σ122+σ232+σ312)}1
(49)
andDisplay Formula
σ1TifTC2
(50)
where σ1 is the largest principal stress. Whichever of 49 or 50 is the more limiting determines the failure condition. Relations 49,50 are as easy to use as are the Mises or Tresca criteria, but 49,50 have far greater generality. The Mises criterion corresponds to 49 with T=C while 50 is seen to be inapplicable at T=C.

This work was performed under the auspices of the U.S. Department of Energy by the University of California, Lawrence Livermore National Laboratory under Contract No. W-7405-Eng-48. This work was partially supported by the Office of Naval Research, Dr. Y. D. S. Rajapakse, Program Manager.

Copyright © 2007 by American Society of Mechanical Engineers
Topics: Brittleness , Stress , Failure
View article in PDF format.

References

Figures

Grahic Jump Location
Figure 1

Biaxial stress state

Grahic Jump Location
Figure 2

Brittle failure and ductile yield/failure

Grahic Jump Location
Figure 3

Mises material, alpha=0, ductile yield and brittle failure regions

Grahic Jump Location
Figure 4

Yield/failure paraboloid

Grahic Jump Location
Figure 5

Shear stress superimposed upon pressure, alpha=infinity limit, Eq. 12

Grahic Jump Location
Figure 6

Spherical annulus with internal pressure

Grahic Jump Location
Figure 7

Failure angles, uniaxial tension and compression

Tables

Table Grahic Jump Location
Table 1
Stresses on failure surfaces, Eqs. 45

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In