Indentation responses of crystalline materials have been found to be radically different at micrometer and nanometer scales. The latter is usually thought to be controlled by the nucleation of dislocations. To explore this physical process, a dislocation mechanics study is performed to determine the conditions for the nucleation of a finite number of dislocations under a two-dimensional wedge indenter, using the Rice-Thomson nucleation criterion. The configurational force on the dislocation consists of the applied force, the image force, and the interaction force between dislocations. Dislocations reach equilibrium positions when the total driving force equals the effective Peierls stress, giving a set of nonlinear equations that can be solved using the Newton-Raphson method. When the apex angle of the wedge indenter increases, the critical contact size for dislocation nucleation increases rapidly, indicating that dislocation multiplication near a blunt wedge tip is extremely difficult. This geometric dependence agrees well with experimental findings.
1. Introduction
The ability to quantitatively model mesoscopic material behavior is a critical step in understanding reliability analysis and material design for applications ranging from aerospace components to nanoelectronics. Many experimental studies in the past decade have shown that the contact plasticity is size-dependent (i.e., indentation size effects, ISEs) and thus is essentially a multiscale problem [1–7]. This phenomenon is believed to be due to the collective behavior of dislocation nucleation and storage under the contact. Statistical theories have been developed to understand the dislocation microstructures that form and evolve during contact events, but most of these theories are only applicable when the contact size is larger than micrometers. A number of recent experiments demonstrate the importance of individual dislocation nucleation events, such as pop-in excursions on the load-displacement
curves [7, 8], which are thought to govern the indentation size effects at nanoscale, while the size effects at microscale are associated with strain gradient effects and geometrically necessary dislocations [1, 2, 7].
A variety of analyses have been carried out using discrete dislocation and phenomenological strain gradient plasticity to model ISE across the length scales [9, 10]. At larger length scale (submicron to micrometer scales), both strain gradient plasticity theories that are based on geometrically necessary dislocations and discrete dislocation simulations have provided important insights in understanding the dependence of microhardness on
contact size. Many recent theories aim to modify these theories for
nanohardness measurements [4, 10], but the connection to the dislocation microstructure still remains elusive. At sufficiently small-length scale, the size of the contact region under high stress may be comparable to the
dislocation spacing, so that the indentation behavior may enter the dislocation-nucleation-controlled regime. This work attempts to explore this line from the study of the dislocation nucleation and pileup under a two-dimensional wedge contact tip.
2. Problem Formulation and Solution Method
Figure 1 gives the problem definition where the apex half angle is α,
the applied load is P (per unit length out of the plane), and the contact size is a. Suppose that a dislocation b=b(cosθ,sinθ,0)T can be nucleated from the indenter tip, then a negative dislocation −b will be left. This dislocation pair will lead to tilting of the indenter, but such a modification is usually very small [11] and will be neglected in this work. The
micromechanical model consists of three parts: (1) the determination of the configurational force on the dislocation, (2) the use of the Rice-Thomson criterion for dislocation nucleation, and (3) the calculation of the equilibrium dislocation positions. In the continuum dislocation modeling, we use the Rice-Thomson dislocation nucleation criterion [12, 13]. If the driving force on a dislocation at a distance η from a stress concentration is greater than the Peierls force, the dislocation can be emitted from the step. The dislocation glides only if the driving force overcomes the Peierls force. In the Rice-Thomson model, parameter η characterizes the size of the dislocation emission process zone. We can either regard it as a material parameter or determine it by atomistic simulations. With the information of η and the effective Peierls stress τp, we can easily decide whether a dislocation nucleates and where the equilibrium position of the nucleated dislocation would be, as long as the driving force on the dislocation is determined.
Schematic illustration and geometric conventions of the contact problem.
The driving force on the dislocation b has three parts:
J=Jappl+Jimage+Jint,
where Jappl is the applied driving force, Jimage is the image force, and Jint is the interaction force due to the
dislocation −b and/or other pairs of dislocations. The
applied driving force can be computed from the stress field [11]:
σxx=2yπa(μ1−ν)cotαπ∫−11cosh−1(1/t)(x/a−t)2dt[(x/a−t)2+(y/a)2]2,σyy=2y3πa3(μ1−ν)cotαπ∫−11cosh−1(1/t)dt[(x/a−t)2+(y/a)2]2,σxy=2y2πa2(μ1−ν)cotαπ∫−11cosh−1(1/t)(x/a−t)dt[(x/a−t)2+(y/a)2]2,
and the Peach-Koehler formula:
Jappl=[(σ⋅b)×(0,0,1)T]⋅(cosθ,sinθ,0)T.
Those curves of Jappl(1−ν)/μb against r/b will collapse onto each other by replotting Jappl(1−ν)/μb against r/a, clearly because the indenter is self-similar. The state of stress at the apex comprises a finite shear stress superposed on an infinite hydrostatic pressure. The applied driving force at the apex can be written as
(1−ν)μbJappl|r→0=cotαπf(θ),
where f(θ) is a dimensionless function that only depends on the slip-plane angle of the order of unity. This stress state is radically different from the stress singularity caused by a crack tip, or a flat-ended punch under normal contact, or an arbitrary indenter under tangential contact [14, 15]. The
implications on the dislocation nucleation behavior will be discussed shortly.
The image force can be determined by considering two opposing semi-infinite cracks in an undislocated plane, subjected to tractions that will cancel the tractions induced by the dislocation in an uncracked plane. Using the complex function analysis, the stress field at z=x+iy, caused by a dislocation b at z0=x0+iy0, in an uncracked plane is given by [16]
σ˜xx+σ˜yy=2[Φ(z)+Φ(z)¯],σ˜yy+iσ˜xy=Φ(z)¯+Ω(z)+(z¯−z)Φ′(z),
where the holomorphic complex functions Φ and Ω are
Φ(z)=Bz−z0,Ω(z)=Bz¯0−z0(z−z0)2+B¯z−z0,B=μπi(1+κ)(bx+iby)=(μ1−ν)14π(−ibx+by),
and κ=3−4ν for plane strain.
Now consider an undislocated plane with two opposing semi-infinite cracks, subjected to the corrected tractions −σ˜(|x|>a,y=0). Following Tada’s solution 4.5 [17], the correction stress field can be determined by
ZI(z)=(∫−∞−a+∫a∞)dt−σ˜yy(t,0)πt2−a2a2−z2(1t−z+2za2),ZII(z)=(∫−∞−a+∫a∞)dtσ˜xy(t,0)πt2−a2a2−z2(1t−z),
where ZI(z) and ZII(z) are the mode-I and mode-II Westergaard stress functions, respectively. The stress components inside the substrate are
σ^xx=(ReZI−yImZI′)+(2ImZII+yReZII′),σ^yy=(ReZI+yImZI′)+(−yReZII′),σ^xy=(−yReZI′)+(ReZII−yImZII′).
In (7), we have assumed that the contact region is full bonded, so that the crack analogy can be used. This is equivalent to the infinite friction condition, or to the contact between two identical solids (regardless of the friction condition). For two different materials under finite frictional contact, the shear stress inside the contact should be considered in the calculation of image force. Finite friction usually does not change the load-displacement curve noticeably [18], but may change the stress distribution and thus lead to a large contribution to the image force.
The image force Jimage can be computed from the stress fields in (8) and the Peach-Koehler formula. The
interaction force Jint can be easily computed from (5) and (8) by substituting −b into b,
and 0 into z0. The image force becomes less important when the contact size increases, since its magnitude decays rapidly with respect to the distance from the free surface. The interaction force decays with respect to the location of the leading dislocation.
We are concerned with the critical load or contact size to nucleate a dislocation. Using the Rice-Thomson criterion, a straightforward dimensional analysis gives that
acrtb=Πa(ηb,τpμ(1−ν),α,θ).
After the dislocation is nucleated, the equilibrium position of the dislocation is determined by J=τpb, leading to
reqa=Πr(ab,ηb,τpμ(1−ν),α,θ),
with a>acrt. It is anticipated that when a≫b, the driving force is dominated by the applied stress field, so that req/a becomes size-independent since the applied stress field is self-similar. The relationship between req/a and a/b gives a qualitative measure of the dislocation density under the indenter.
When there are many dislocations, the driving force on the kth dislocation consists of (1) and the interaction force from other pair of dislocations. The latter can be easily evaluated from (5) and (8). To determine the equilibrium positions of N dislocations, we use the Newton-Raphson method to solve N force-balance equations. The root finding process is very slow because of the N-dimensional functions and the shallow slope of driving force curves near equilibria.
3. Results and Discussion
Figure 2 shows the total driving
force with respect to the dislocation location. The zone where the interaction force dominates does not proportionally increase with the contact size. Consequently, the driving force on a fictitious dislocation at r=η will eventually be nonzero and larger than the Peierls stress, assuming that the Peierls stress is smaller than Jappl(r=0).
The configurational force on the dislocation plotted against the dislocation position with varying contact size.
Figure 3(a) shows the critical
contact size to nucleate one and many dislocations. The Rice-Thomson model
introduces a length scale η, so that the first dislocation nucleation occurs at a finite contact size for this otherwise self-similar contact problem. It is found that the larger the wedge half-angle, the more difficult to nucleate a dislocation. The image force and interaction force are independent of the wedge angle. The applied driving force, however, is proportional to cotα. In addition, the shear stress near the wedge tip is bounded, so that there is a critical angle above which the maximum applied driving force is lower than the Peierls stress. The second and subsequent dislocations occur at higher contact size, since the existing dislocation exerts strong back stress to prevent further dislocation nucleation. The diverging trend with respect to α is more evident, suggesting that dislocation multiplication for blunt wedges is extremely difficult. Consequently, when we are concerned with the effective dislocation number and density under a wedge contact, we must incorporate the bulk dislocation nucleation sources. The analysis presented in this paper is therefore only applicable for nanocontact experiments. This indenter geometry effects were also observed in a series of experiments [6, 19] where the sharper cube corner indenter tip (with half angle of ~35º) shows higher hardness value as compared to the Berkovich indenter tip (with half angle of ~70º) at the same indentation depth in the nanometer length scale, resulting in more dramatic indentation size effects.
(a) The critical contact size to nucleate dislocations with respect to the wedge half-angle. (b) The equilibrium positions of nucleated dislocations with respect to the contact size.
Figure 3(b) shows the normalized dislocation equilibrium positions, req/a, with respect to the contact size. When a/b→∞, (10) will be independent of the contact size, and the effective plasticity zone will
be proportional to the contact size. Consequently, the dislocation density is approximately linear with b/a at large contact size. (Of course, when the contact size is large, it is the bulk nucleation sources and the statistically stored dislocations that govern the plastic zone and contact behavior.) However, this relationship will break down for small-contact size, because of the discrete nature of the dislocation nucleation and motion, and the length scale η. It thus becomes inappropriate to estimate plastic zone size in terms of contact size and relate the dislocation density to hardness in the nanohardness measurements. A more rational treatment would require a closer look at the connection between the mean contact pressure and the dislocation microstructure evolution at nanometer length scale.
4. Summary
A dislocation mechanics study is performed here to examine the conditions for dislocation nucleation from, and pileup under, a wedge indenter. The configurational force on the dislocation is evaluated by linear elastic analysis, and the Rice-Thomson criterion is used for the dislocation nucleation. The contact strength map in Figure 3(a) shows that dislocation nucleation and multiplication are extremely difficult for blunt indenters, mainly because the stress state near the wedge tip consists of a finite shear stress and an infinite hydrostatic component. The ratio of the effective plastic zone size to the contact size is found to approach a constant when there are sufficient numbers of dislocations at large contact size. However, at nanoscale we should directly investigate the relationship between contact
pressure and dislocation microstructure, instead of using dislocation-density-based hardening law.
Acknowledgments
Y. F. Gao would like to acknowledge the support from the Joint Institute of Advanced Materials at the University of Tennessee. J. Lou gratefully acknowledges the startup fund from Rice University.
NixW. D.GaoH.Indentation size effects in crystalline materials: a law for strain gradient plasticity199846341142510.1016/S0022-5096(97)00086-0SwadenerJ. G.GeorgeE. P.PharrG. M.The correlation of the indentation size effect measured with indenters of various shapes200250468169410.1016/S0022-5096(01)00103-XLouJ.ShrotriyaP.BuchheitT.YangD.SoboyejoW. O.Nanoindentation study of plasticity length scale effects in LIGA Ni microelectromechanical systems structures2003183719728FengG.NixW. D.Indentation size effect in MgO200451659960310.1016/j.scriptamat.2004.05.034ZongZ.LouJ.AdewoyeO. O.ElmustafaA. A.HammadF.SoboyejoW. O.Indentation size effects in the nano- and micro-hardness of fcc single crystal metals20064341-217818710.1016/j.msea.2006.06.137LouJ.ShrotriyaP.AllamehS.BuchheitT.SoboyejoW. O.Strain gradient plasticity length scale parameters for LIGA Ni MEMs thin films20064411-229930710.1016/j.msea.2006.08.048NixW. D.GreerJ. R.FengG.LilleoddenE. T.Deformation at the nanometer and micrometer length scales: effects of strain gradients and dislocation starvation200751563152315710.1016/j.tsf.2006.01.030BeiH.GeorgeE. P.HayJ. L.PharrG. M.Influence of indenter tip geometry on elastic deformation during nanoindentation20059541404550110.1103/PhysRevLett.95.045501WidjajaA.NeedlemanA.van der GiessenE.The effect of indenter shape on sub-micron indentation according to discrete dislocation plasticity2007151S121S13110.1088/0965-0393/15/1/S11HuangY.ZhangF.HwangK. C.NixW. D.PharrG. M.FengG.A model of size effects in nano-indentation20065481668168610.1016/j.jmps.2006.02.002JohnsonK. L.1985Cambridge, UKCambridge University PressRiceJ. R.ThomsonR.Ductile versus brittle behavior of crystals1973291739710.1080/14786437408213555YuH. H.ShrotriyaP.GaoY. F.KimK.-S.Micro-plasticity of surface steps under adhesive contact: part I-surface yielding controlled by single-dislocation nucleation200755348951610.1016/j.jmps.2006.09.003GaoY. F.LucasB. N.HayJ. C.OliverW. C.PharrG. M.Nanoscale incipient asperity sliding and interface micro-slip assessed by the measurement of tangential contact stiffness200655765365610.1016/j.scriptamat.2006.05.006GaoY. F.XuH. T.OliverW. C.PharrG. M.A comparison of Coulomb friction and friction stress models based on multidimensional nanocontact experimentsto appear in Journal of Applied MechanicsSuoZ.Singularities interacting with interfaces and cracks198925101133114210.1016/0020-7683(89)90072-3TadaH.ParisP. C.IrwinG. R.2000New York, NY, USAASME PressGaoY. F.XuH. T.OliverW. C.PharrG. M.Effective elastic modulus of film-on-substrate systems under normal and tangential contact2008562402416MirshamsR. A.ParakalaP.Nanoindentation of nanocrystalline Ni with geometrically different indenters20043721-225226010.1016/j.msea.2004.01.010