《纺织复合材料》课程参考文献（Micromechanisms of Fracture and Fatigue）Appendix（Ab initio Computational Methods、Mixed-mode Criteria of Crack Stability、Plastic Flow Rate Inside the Neck）

Appendix A Ab initio Computational Methods The objective of ab initio approaches based on applied quantum theory is to calculate stationary states for electrons in the electrostatic field of nuclei,i.e., the electronic structure (ES).The energy of this ground state can then serve as a potential energy for displacements of nuclei.From the point of view of ideal-strength (IS)calculations,the total energy of the system is the most important output of ab initio methods. The main advantage of ab initio methods is their independence of exper- imental data.Unlike in the case of empirical and semi-empirical methods, there is no need for calibration parameters.Thus,they can also be used for calculation of some structural and mechanical characteristics of hypothetical systems(prediction of properties of materials that have not yet been devel- oped)or study of materials behaviour under large deformations(far from the equilibrium state)that can give a better understanding of micromechanisms of materials failure. First attempts to develop applicable theories were made in the late 1920s [421,422,a few years after the foundations of modern quantum the- ory were laid(derivation of the Schrodinger equation).A very successful step forward was the Hartree-Fock method [422,423].This method yields very accurate bond lengths in molecules.On the other hand,the binding energies are generally not in good agreement with experimentally obtained energies. Moreover,for solids,the Hartree-Fock method has problems with a descrip- tion of band structures.The density functional theory (DFT)37,38]was invented to include correlation effects without using the very costly wave function methods.All the methods used within this book are based on the DFT. In DFT the energy is not obtained as eigenvalues of a wave function,but rather as a functional of the electron density.The complex problem of many interacting electrons is transformed into a much simpler study of single elec- tron interactions with an effective potential Uef created by other electrons and all nuclei.This is expressed by the Kohn-Sham equation (one-electron Schrodinger equation) 249

Appendix A Ab initio Computational Methods The objective of ab initio approaches based on applied quantum theory is to calculate stationary states for electrons in the electrostatic field of nuclei, i.e., the electronic structure (ES). The energy of this ground state can then serve as a potential energy for displacements of nuclei. From the point of view of ideal-strength (IS) calculations, the total energy of the system is the most important output of ab initio methods. The main advantage of ab initio methods is their independence of exper￾imental data. Unlike in the case of empirical and semi-empirical methods, there is no need for calibration parameters. Thus, they can also be used for calculation of some structural and mechanical characteristics of hypothetical systems (prediction of properties of materials that have not yet been devel￾oped) or study of materials behaviour under large deformations (far from the equilibrium state) that can give a better understanding of micromechanisms of materials failure. First attempts to develop applicable theories were made in the late 1920s [421, 422], a few years after the foundations of modern quantum the￾ory were laid (derivation of the Schr¨odinger equation). A very successful step forward was the Hartree–Fock method [422, 423]. This method yields very accurate bond lengths in molecules. On the other hand, the binding energies are generally not in good agreement with experimentally obtained energies. Moreover, for solids, the Hartree–Fock method has problems with a descrip￾tion of band structures. The density functional theory (DFT) [37, 38] was invented to include correlation effects without using the very costly wave function methods. All the methods used within this book are based on the DFT. In DFT the energy is not obtained as eigenvalues of a wave function, but rather as a functional of the electron density. The complex problem of many interacting electrons is transformed into a much simpler study of single elec￾tron interactions with an effective potential Ueff created by other electrons and all nuclei. This is expressed by the Kohn–Sham equation (one-electron Schr¨odinger equation) 249

250 A Ab initio Computational Methods (-△+Uef(r)-e)(r)=0, where si represents the one-electron energies and i are the one-electron wave functions. The wave functions are then occupied in accordance with the Pauli prin- ciple and a new field is obtained by solving the Poisson equation for point nuclei shielded by the electronic density p(r)=>lv:(r)2. i.occ In the case of periodic crystalline materials,the one-electron wave func- tions are expanded into appropriate basis sets and satisfy the Bloch theorem. The Kohn-Sham equation is solved iteratively until the solution becomes self-consistent.This means that the electron density,determined from the effective one-electron potential,must generate the same effective potential (which is again a functional of the electron density).The self-consistent cycle usually starts with a guess of the effective potential(superposition of atomic- like potentials)and then the input and output potentials are appropriately mixed before starting a new iteration.The quality and speed of the conver- gence of such calculations is related not only to the choice of a suitable basis but also to the sophistication of the iterative process.The necessary correc- tions for exchange and correlation are also included in the effective potential Uef.This seems to be the crucial point of ab initio calculations because the exchange-correlation(XC)functional is not known exactly and must be ap- proximated.The first (and the simplest)attempt to build an approximation of the XC energy functional in the DFT is the local density approximation (LDA)[424].The LDA is local in the sense that the electron exchange and correlation energy at any point in space is a functional of the electron density at that point only.As a consequence of this,LDA fails in situations where the density undergoes rapid changes (molecules).An improvement to this can be made by considering the gradient of the electron density.The density gradient corrections are implemented in the so-called Generalized Gradient Approximation(GGA).While there is only one LDA there are several differ- ent parameterizations of the GGA [425-427]. Nevertheless,it is the choice of the basis wave functions that makes the main difference between various methods used for the ES calculations.The better we choose them (according to the character of the problem),the smaller number of them is needed for the description of one-electron wave functions.Commonly used bases are augmented (APW)and orthogonalized (OPW)plane waves,linear muffin-tin orbitals (LMTO),linear combination of atomic orbitals (LCAO),of Gaussian orbitals (LCGO)and linear aug- mented Slater-type orbitals (LASTO),augmented spherical waves (ASW), etc.The Korringa-Kohn-Rostoker(KKR)method proceeds by the use of the Green function of the Kohn-Sham equation and is also called Green's func-

250 A Ab initio Computational Methods (−Δ + Ueff (r) − εi)ψi(r)=0, where εi represents the one-electron energies and ψi are the one-electron wave functions. The wave functions are then occupied in accordance with the Pauli prin￾ciple and a new field is obtained by solving the Poisson equation for point nuclei shielded by the electronic density ρ(r) =  i,occ |ψi(r)| 2 . In the case of periodic crystalline materials, the one-electron wave func￾tions are expanded into appropriate basis sets and satisfy the Bloch theorem. The Kohn–Sham equation is solved iteratively until the solution becomes self-consistent. This means that the electron density, determined from the effective one-electron potential, must generate the same effective potential (which is again a functional of the electron density). The self-consistent cycle usually starts with a guess of the effective potential (superposition of atomic￾like potentials) and then the input and output potentials are appropriately mixed before starting a new iteration. The quality and speed of the conver￾gence of such calculations is related not only to the choice of a suitable basis, but also to the sophistication of the iterative process. The necessary correc￾tions for exchange and correlation are also included in the effective potential Ueff . This seems to be the crucial point of ab initio calculations because the exchange-correlation (XC) functional is not known exactly and must be ap￾proximated. The first (and the simplest) attempt to build an approximation of the XC energy functional in the DFT is the local density approximation (LDA) [424]. The LDA is local in the sense that the electron exchange and correlation energy at any point in space is a functional of the electron density at that point only. As a consequence of this, LDA fails in situations where the density undergoes rapid changes (molecules). An improvement to this can be made by considering the gradient of the electron density. The density gradient corrections are implemented in the so-called Generalized Gradient Approximation (GGA). While there is only one LDA there are several differ￾ent parameterizations of the GGA [425–427]. Nevertheless, it is the choice of the basis wave functions that makes the main difference between various methods used for the ES calculations. The better we choose them (according to the character of the problem), the smaller number of them is needed for the description of one-electron wave functions. Commonly used bases are augmented (APW) and orthogonalized (OPW) plane waves, linear muffin-tin orbitals (LMTO), linear combination of atomic orbitals (LCAO), of Gaussian orbitals (LCGO) and linear aug￾mented Slater-type orbitals (LASTO), augmented spherical waves (ASW), etc. The Korringa–Kohn–Rostoker (KKR) method proceeds by the use of the Green function of the Kohn–Sham equation and is also called Green’s func-

A.1 TB-LMTO-ASA Code 251 tion (GF)method.The pseudopotential approach applied mostly to solids containing no d-or f-electrons is also widely used.A detailed description of these methods may be found in many books and articles,e.g.,in 86,88,89. The atomic configurations corresponding to deformed structures usually have lower symmetries and,at the strength limit,they are very far from the lowest- energy equilibrium state.Therefore,to obtain reliable structural energy dif ferences,the full-potential methods (i.e.,without any shape approximation of the crystal potential and the electronic charge density)have to be utilized in such studies.At present,several codes are available,e.g.,WIEN,VASP, FHI,FLEUR,FPLO,FPLMTO,ABINIT,SIESTA,etc. All the ab initio calculations of IS presented in this book were performed by using the following three computational codes:LMTO-ASA,WIEN and VASP. A.1 TB-LMTO-ASA Code Linearized ab initio methods have been successfully utilized for solving many problems in materials science [428,429].One of the most effective approaches for early first principles calculations was the LMTO formalism which has been continuously developed since 1980 [428-430.This method is very appropriate for self-consistent calculations. In the LMTO-ASA code [431],the crystal potential U is approximated by a muffin-tin(MT)shape potential which is composed of a set of spherically symmetric potentials inside slightly overlapping spheres around individual nuclei and a constant potential in the interstitial region outside the spheres (Figure A.1).Atomic-like orbitals derived for the MT potential constitute a suitable basis set.The tight binding approximation,which is also imple- mented into the code [432],assumes that the full Hamiltonian may be ap- proximated by that of an isolated atom centred at each lattice point.The atomic orbitals (eigenfunctions of the single-atom Hamiltonian)are assumed to be negligible at distances exceeding the lattice constant. In all presented calculations,the LMTO method is used in the framework of an atomic sphere approximation(ASA)which is particularly suitable for closely packed structures like fcc,hcp or bcc [428].The size of a spherically symmetric potential is assumed to be equal to that of the Wigner-Seitz cell (Figure A.1(b)).This suppresses the interstitial region and neglects the ki- netic energy of related free electrons.Owing to the necessary space-filling condition,ASA represents a physically plausible model only for a description of an infinite system of atomic spheres

A.1 TB-LMTO-ASA Code 251 tion (GF) method. The pseudopotential approach applied mostly to solids containing no d- or f-electrons is also widely used. A detailed description of these methods may be found in many books and articles, e.g., in [86, 88, 89]. The atomic configurations corresponding to deformed structures usually have lower symmetries and, at the strength limit, they are very far from the lowest￾energy equilibrium state. Therefore, to obtain reliable structural energy dif￾ferences, the full-potential methods (i.e., without any shape approximation of the crystal potential and the electronic charge density) have to be utilized in such studies. At present, several codes are available, e.g., WIEN, VASP, FHI, FLEUR, FPLO, FPLMTO, ABINIT, SIESTA, etc. All the ab initio calculations of IS presented in this book were performed by using the following three computational codes: LMTO-ASA, WIEN and VASP. A.1 TB-LMTO-ASA Code Linearized ab initio methods have been successfully utilized for solving many problems in materials science [428,429]. One of the most effective approaches for early first principles calculations was the LMTO formalism which has been continuously developed since 1980 [428–430]. This method is very appropriate for self-consistent calculations. In the LMTO-ASA code [431], the crystal potential U is approximated by a muffin-tin (MT) shape potential which is composed of a set of spherically symmetric potentials inside slightly overlapping spheres around individual nuclei and a constant potential in the interstitial region outside the spheres (Figure A.1). Atomic-like orbitals derived for the MT potential constitute a suitable basis set. The tight binding approximation, which is also imple￾mented into the code [432], assumes that the full Hamiltonian may be ap￾proximated by that of an isolated atom centred at each lattice point. The atomic orbitals (eigenfunctions of the single-atom Hamiltonian) are assumed to be negligible at distances exceeding the lattice constant. In all presented calculations, the LMTO method is used in the framework of an atomic sphere approximation (ASA) which is particularly suitable for closely packed structures like fcc, hcp or bcc [428]. The size of a spherically symmetric potential is assumed to be equal to that of the Wigner–Seitz cell (Figure A.1(b)). This suppresses the interstitial region and neglects the ki￾netic energy of related free electrons. Owing to the necessary space-filling condition, ASA represents a physically plausible model only for a description of an infinite system of atomic spheres

252 A Ab initio Computational Methods 5 0 3 3 1012345 G- [001] [110] -7012 [001) 34 [110] (a) (b) Figure A.1 Muffin-tin shape potential in(110)plane of a general bcc crystal of the lattice constant a =3au with spherical radii:(a)rMT 1au,and (b)rMT =rws 1.48au A.2 Wien 95-w2k Codes The program package WIEN [433]does not use any shape approximation to the potential.The crystal environment is divided into a region of non- overlapping atomic spheres(centred at individual atomic sites)and an inter- stitial region as can be seen in Figure A.2.In order to describe ES reliably and effectively,two different basis sets are employed:a product of linear com- bination of radial functions and spherical harmonics is used inside the spheres whereas the wave functions in the interstitial region are expanded into a linear combination of plane waves.The solution in both regions must be continuous at the spherical boundaries.Each basis function is then defined as a plane- wave in the interstitial region connected smoothly to a linear combination of atomic-like functions in the spheres,thus providing an efficient represen- tation throughout the space.A similar representation is used for potentials and charge densities.The method is called the linear augmented plane wave (LAPW)method [428,434]. S Figure A.2 Illustration of a crystal model-three atomic spheres (S1-S3)with (embedded in the interstitial region I with

252 A Ab initio Computational Methods (a) (b) Figure A.1 Muffin-tin shape potential in (110) plane of a general bcc crystal of the lattice constant a = 3 au with spherical radii: (a) rMT = 1 au, and (b) rMT = rW S = 1.48 au A.2 Wien 95 – w2k Codes The program package WIEN [433] does not use any shape approximation to the potential. The crystal environment is divided into a region of non￾overlapping atomic spheres (centred at individual atomic sites) and an inter￾stitial region as can be seen in Figure A.2. In order to describe ES reliably and effectively, two different basis sets are employed: a product of linear com￾bination of radial functions and spherical harmonics is used inside the spheres whereas the wave functions in the interstitial region are expanded into a linear combination of plane waves. The solution in both regions must be continuous at the spherical boundaries. Each basis function is then defined as a plane￾wave in the interstitial region connected smoothly to a linear combination of atomic-like functions in the spheres, thus providing an efficient represen￾tation throughout the space. A similar representation is used for potentials and charge densities. The method is called the linear augmented plane wave (LAPW) method [428, 434]. S1 S2 S3 I Figure A.2 Illustration of a crystal model – three atomic spheres (S1 – S3) with potential US(r) = 2 lm Ulm(r)Ylm(ˆr) embedded in the interstitial region I with UI (r) = 2 K UK(r)eiKr

A.3 VASP Code 253 In order to increase the flexibility of the basis (to improve upon the lin- earization of wave functions)and to make possible a consistent treatment of semicore and valence states in one energy window(to ensure orthogonality) additional basis functions can be added.They are called local orbitals [435 and consist of a linear combination of two radial functions at two different energies (e.g.,at the 3s and 4s energy)and one energy derivative (at one of these energies).The local orbitals are normalized and have zero value and slope at the spherical boundaries. A.3 VASP Code Another way of avoiding a problem with plane wave basis set in the vicinity of atomic nuclei,where the number of plane waves would exceed any practical limits(perhaps except for H or Li),is to substitute the exact potential by a pseudopotential. LU cut Ψ al-electron cut pseudo U- Figure A.3 Comparison of a wavefunction in the Coulomb potential of the nucleus (dashed lines)to that in the pseudopotential (solid lines) Construction of the pseudo-wavefunctions is schematically described in Figure A.3.The Coulomb potential of the nucleus and the corresponding wavefunction are represented by dashed lines.The solid line displays the pseudopotential and the pseudo-wavefunction.The real and the pseudo-

A.3 VASP Code 253 In order to increase the flexibility of the basis (to improve upon the lin￾earization of wave functions) and to make possible a consistent treatment of semicore and valence states in one energy window (to ensure orthogonality) additional basis functions can be added. They are called local orbitals [435] and consist of a linear combination of two radial functions at two different energies (e.g., at the 3s and 4s energy) and one energy derivative (at one of these energies). The local orbitals are normalized and have zero value and slope at the spherical boundaries. A.3 VASP Code Another way of avoiding a problem with plane wave basis set in the vicinity of atomic nuclei, where the number of plane waves would exceed any practical limits (perhaps except for H or Li), is to substitute the exact potential by a pseudopotential. Upseudo r rcut pseudo r rcut  U 0 0 Figure A.3 Comparison of a wavefunction in the Coulomb potential of the nucleus (dashed lines) to that in the pseudopotential (solid lines) Construction of the pseudo-wavefunctions is schematically described in Figure A.3. The Coulomb potential of the nucleus and the corresponding wavefunction are represented by dashed lines. The solid line displays the pseudopotential and the pseudo-wavefunction. The real and the pseudo-

254 A Ab initio Computational Methods wavefunction (and also the potentials)match above a certain cutoff radius Tcut· The pseudopotential approach has been implemented in the VASP (Vi- enna ab initio Simulation Package)code.This code was developed at the Institut fiir Materialphysik,Universitat Wien.The VASP currently supports three types of pseudopotentials:norm-conserving pseudopotentials,ultrasoft pseudopotentials and projector augmented wave pseudopotentials [436-438. In all three cases,the core electrons (at lower energy levels than valence elec- trons)are precalculated in an atomic environment and kept "frozen"during the remaining calculations

254 A Ab initio Computational Methods wavefunction (and also the potentials) match above a certain cutoff radius rcut. The pseudopotential approach has been implemented in the VASP (Vi￾enna ab initio Simulation Package) code. This code was developed at the Institut f¨ur Materialphysik, Universit¨at Wien. The VASP currently supports three types of pseudopotentials: norm-conserving pseudopotentials, ultrasoft pseudopotentials and projector augmented wave pseudopotentials [436–438]. In all three cases, the core electrons (at lower energy levels than valence elec￾trons) are precalculated in an atomic environment and kept “frozen” during the remaining calculations

Appendix B Mixed-mode Criteria of Crack Stability In order to describe the crack stability under mixed-mode loading,various concepts within the framework of LEFM were proposed (e.g.,[162-164).In this brief overview,the stress intensity factors are denoted only by capital letters KI,KI and KI,as is common for remote quantities.However,the criteria are relevant also to local stress intensity factors usually denoted as k1,k2 and k3. B.1 Energy Criterion The criterion postulates that the total energy Gr+Gr+Gr,released by the system to create a unit of a new surface,equals the crack growth resistance Ge: Ge=GI+GI+GI1I. (B.1) When the crack grows in its original plane(or propagates in a self-similar manner),the energy criterion can be expressed in terms of the effective K- factor as K=VK+K+年 (B.2) where Ke is the critical stress intensity factor (fracture toughness),K=3-4v for plane strain,K=(3-v)/(1+v)for plane stress and v is the Poisson's ratio. 255

Appendix B Mixed-mode Criteria of Crack Stability In order to describe the crack stability under mixed-mode loading, various concepts within the framework of LEFM were proposed (e.g., [162–164]). In this brief overview, the stress intensity factors are denoted only by capital letters KI , KII and KIII , as is common for remote quantities. However, the criteria are relevant also to local stress intensity factors usually denoted as k1, k2 and k3. B.1 Energy Criterion The criterion postulates that the total energy GI+GII+GIII , released by the system to create a unit of a new surface, equals the crack growth resistance Gc: Gc = GI + GII + GIII . (B.1) When the crack grows in its original plane (or propagates in a self-similar manner), the energy criterion can be expressed in terms of the effective K￾factor as Kc = K2 I + K2 II + 4 κ + 1 K2 III , (B.2) where Kc is the critical stress intensity factor (fracture toughness), κ = 3−4ν for plane strain, κ = (3 − ν)/(1 + ν) for plane stress and ν is the Poisson’s ratio. 255

256 B Mixed-mode Criteria of Crack Stability B.2 Criterion of Linear Damage Accumulation Each of Equations B.1 and B.2 can be understood as a condition of subse- quently reaching the critical level of material damage during loading.If the partial damages accumulate independently,the mixed-mode criterion +m+G皿≤1 GI GIe GIle GIlle can be accepted.Clearly,for the special case GIc GIIc GIIIe Ge, the criterion of linear damage accumulation transfers to the energy criterion. Under the dominant mode I loading,the problem of shear friction in modes II and III on the rough fracture surfaces is not crucial due to sufficiently large opening displacements of crack flanks.Consequently,the equality of critical crack driving forces in all modes might be assumed and the energy criterion can be accepted. Since the assumption of self-similar crack propagation does not hold un- der a general remote mixed-mode loading,one must be careful when applying both the above-mentioned criteria.They appear to be useful only for a predic- tion of the moment of unstable fracture or,generally,in the case of a dominant external mode I loading of microscopically tortuous cracks.The description of long fatigue crack propagation under external mixed-mode loading condi- tions by means of these criteria is rather limited in the sense that the growth direction cannot be predicted. B.3 Criterion of Minimal Deformation Energy The deformation energy density S is minimal in those elements near the crack front,where the ratio of the hydrostatic tensile stress and the octahedral stress is maximal [439.In such elements,therefore,the first local fracture is expected to occur and the crack propagates towards such damaged sites in the matrix.In the plane mixed-mode I+II,the condition of minimal S yields the stability condition as 16G Ka=2R-aea)K+2aaam)Kki+a2amKa, where ai(m)are angular functions and G is the shear modulus.The criterion generally enables us to predict the moment of crack nucleation as well as the crack growth direction in the 3D homogeneous continuum.Its disadvantage is a rather low transparency

256 B Mixed-mode Criteria of Crack Stability B.2 Criterion of Linear Damage Accumulation Each of Equations B.1 and B.2 can be understood as a condition of subse￾quently reaching the critical level of material damage during loading. If the partial damages accumulate independently, the mixed-mode criterion GI GIc + GII GIIc + GIII GIIIc ≤ 1 can be accepted. Clearly, for the special case GIc = GIIc = GIIIc = Gc, the criterion of linear damage accumulation transfers to the energy criterion. Under the dominant mode I loading, the problem of shear friction in modes II and III on the rough fracture surfaces is not crucial due to sufficiently large opening displacements of crack flanks. Consequently, the equality of critical crack driving forces in all modes might be assumed and the energy criterion can be accepted. Since the assumption of self-similar crack propagation does not hold un￾der a general remote mixed-mode loading, one must be careful when applying both the above-mentioned criteria. They appear to be useful only for a predic￾tion of the moment of unstable fracture or, generally, in the case of a dominant external mode I loading of microscopically tortuous cracks. The description of long fatigue crack propagation under external mixed-mode loading condi￾tions by means of these criteria is rather limited in the sense that the growth direction cannot be predicted. B.3 Criterion of Minimal Deformation Energy The deformation energy density S is minimal in those elements near the crack front, where the ratio of the hydrostatic tensile stress and the octahedral stress is maximal [439]. In such elements, therefore, the first local fracture is expected to occur and the crack propagates towards such damaged sites in the matrix. In the plane mixed-mode I+II, the condition of minimal S yields the stability condition as K2 eff = 16G 2 (κ − 1)  a11(θm)K2 I + 2a12(θm)KIKII + a22(θm)K2 II , where aij (θm) are angular functions and G is the shear modulus. The criterion generally enables us to predict the moment of crack nucleation as well as the crack growth direction in the 3D homogeneous continuum. Its disadvantage is a rather low transparency

B.4 Criterion of Maximal Principal Stress 257 B.4 Criterion of Maximal Principal Stress This criterion can be derived from the energy criterion by expressing the effective crack driving force Gef as a function of growth angle 0 and by searching its maximal value.It can be shown that this maximum corresponds to a crack perpendicular to the maximum principal stress.Considering the plane model,stresses in close vicinity of the crack tip can be expressed as follows: 1 00=- COS V2Tr Tre= cos(KIsin0+Kn(3cos0-1)). 2V2Tr If Tro=0,then the hoop stress oe becomes the principal stress.This condition determines the growth angle 0m as well as the related effective stress intensity value: +8 Kef =Ki cos3Kncos 0m 2 2 sin 2 The criterion can be used for predictions of crack stability and growth in homogeneous materials.The criterion does not involve the antiplane mode III and,consequently,it can be applied only in the framework of 2D models. Nevertheless,a very good applicability of the maximum principal stress crite- rion to propagation of long fatigue cracks was already proven(e.g.,[399,440). This success can be understood from the micromechanical point of view.In- deed,the fatigue crack tries to maximize the opening as well as the friction between the flanks.Consequently,it inclines to the maximum opening mode I perpendicularly to the maximal main stress direction,thereby also minimiz- ing the shear modes (friction)

B.4 Criterion of Maximal Principal Stress 257 B.4 Criterion of Maximal Principal Stress This criterion can be derived from the energy criterion by expressing the effective crack driving force Geff as a function of growth angle θ and by searching its maximal value. It can be shown that this maximum corresponds to a crack perpendicular to the maximum principal stress. Considering the plane model, stresses in close vicinity of the crack tip can be expressed as follows: σθ = 1 √ 2πr cos θ 2  KI cos2 θ 2 − 3 2 KII sin θ  , τrθ = 1 2 √ 2πr cos θ 2 (KI sin θ + KII (3 cos θ − 1)). If τrθ = 0, then the hoop stress σθ becomes the principal stress. This condition determines the growth angle θm as well as the related effective stress intensity value: tan θm 2 = KI 4KII ± 1 4  KI KII 2 + 8, Keff = KI cos3 θm 2 − 3KII cos2 θm 2 sin θm 2 . The criterion can be used for predictions of crack stability and growth in homogeneous materials. The criterion does not involve the antiplane mode III and, consequently, it can be applied only in the framework of 2D models. Nevertheless, a very good applicability of the maximum principal stress crite￾rion to propagation of long fatigue cracks was already proven (e.g., [399,440]). This success can be understood from the micromechanical point of view. In￾deed, the fatigue crack tries to maximize the opening as well as the friction between the flanks. Consequently, it inclines to the maximum opening mode I perpendicularly to the maximal main stress direction, thereby also minimiz￾ing the shear modes (friction)

Appendix C Plastic Flow Rate Inside the Neck C.1 Ideal Model When taking Equations2.25,2.26 and rangesλ1∈(0.8,2)，λ2∈(0,1)into account,the growth of the semi-axes of elliptical voids at the onset of necking can be described as a aoe0.8ep, b≈bo三ao (C.1) while a=a0e2sp,b≈a0ep (C.2) holds in the final stages of necking just before the final fracture.In order to calculate the active volume V2,the surface area of the elliptical void Se has to be considered as Se4兰2r(b2+ab). (C.3) By inserting Equations C.1 and C.2 into Equation C.3 and transferring the surface of the ellipsoid to that of the equivalent sphere of the same surface area one obtains aSel =4mage.,Se =4ndpe2.e. A corresponding error of this transfer is less than 10%for 0<p<1.Thus, the variation of the void surface during the necking process can be expressed as Se1≈4 rageRep, (C.4) where 0.4<K<2.5.Let us denote 6 the mean distance from the void surface to the outer boundary of the volume V2(see Figure 2.30)and n the number 259

Appendix C Plastic Flow Rate Inside the Neck C.1 Ideal Model When taking Equations 2.25, 2.26 and ranges λ1 ∈ 0.8, 2, λ2 ∈ 0, 1 into account, the growth of the semi-axes of elliptical voids at the onset of necking can be described as a = a0e0.8εp , b ≈ b0 ≡ a0 (C.1) while a = a0e2εp , b ≈ a0eεp (C.2) holds in the final stages of necking just before the final fracture. In order to calculate the active volume V2, the surface area of the elliptical void Sel has to be considered as Sel ∼= 2π  b2 + ab . (C.3) By inserting Equations C.1 and C.2 into Equation C.3 and transferring the surface of the ellipsoid to that of the equivalent sphere of the same surface area one obtains aSel = 4πa2 0e0,4εp , b Sel = 4πa2 0e2,5εp . A corresponding error of this transfer is less than 10% for 0 < εp < 1. Thus, the variation of the void surface during the necking process can be expressed as Sel ≈ 4πa2 0eκεp , (C.4) where 0.4 <κ< 2.5 . Let us denote δ the mean distance from the void surface to the outer boundary of the volume V2 (see Figure 2.30) and n the number 259