Two new expanding cavity models for indentation deformations of elastic strain-hardening materials

Two expanding cavity models (ECMs) are developed for describing indentation deformations of elastic power-law hardening and elastic linear-hardening materials. The derivations are based on two elastic–plastic solutions for internally pressurized thick-walled spherical shells of strain-hardening materials. Closed-form formulas are provided for both conical and spherical indentations, which explicitly show that for a given indenter geometry indentation hardness depends on Young’s modulus, yield stress and strain-hardening index of the indented material. The two new models reduce to Johnson’s ECM for elastic-perfectly plastic materials when the strain-hardening effect is not considered. The sample numerical results obtained using the two newly developed models reveal that the indentation hardness increases with the Young’s modulus and strain-hardening level of the indented material. For conical indentations the values of the indentation hardness are found to depend on the sharpness of the indenter: the sharper the indenter, the larger the hardness. For spherical indentations it is shown that the hardness is significantly affected by the strain-hardening level when the indented material is stiff (i.e., with a large ratio of Young’s modulus to yield stress) and/or the indentation depth is large. When the indentation depth is small such that little or no plastic deformation is induced by the spherical indenter, the hardness appears to be independent of the strain-hardening level. These predicted trends for spherical indentations are in fairly good agreement with the recent finite element results of Park and Pharr.     

Influence of penetration depth and mechanical properties on contact radius determination for spherical indentation

Knowledge of the relationship between the penetration depth and the contact radius is required in order to determine the mechanical properties of a material starting from an instrumented indentation test. The aim of this work is to propose a new penetration depth–contact radius relationship valid for most metals which are deformed plastically by parabolic and spherical indenters. Numerical simulation results of the indentation of an elastic–plastic half-space by a frictionless rigid paraboloïd of revolution show that the contact radius–indentation depth relationship can be represented by a power law, which depends on the reduced Young’s modulus of the contact, on the strain hardening exponent and on the yield stress of the indented material. In order to use the proposed formulation for experimental spherical indentations, adaptation of the model is performed in the case of a rigid spherical indenter. Compared to the previous formulations, the model proposed in the present study for spherical indentation has the advantage of being accurate in the plastic regime for a large range of contact radii and for materials of well-developed yield stress. Lastly, a simple criterion, depending on the material mechanical properties, is proposed in order to know when piling-up appears for the spherical indentation.     

Comparison between Berkovich,Vickers and conical indentation tests: A three-dimensional numerical simulation study

Three-dimensional numerical simulations of Berkovich, Vickers and conical indenter hardness tests were carried out to investigate the influence of indenter geometry on indentation test results of bulk and composite film/substrate materials. The strain distributions obtained from the three indenters tested were studied, in order to clarify the differences in the load–indentation depth curves and hardness values of both types of materials. For bulk materials, the differentiation between the results obtained with the three indenters is material sensitive. The indenter geometry shape factor, β, for evaluating Young’s modulus for each indenter, was also estimated. Depending on the indenter geometry, distinct mechanical behaviours are observed for composite materials, which are related to the size of the indentation region in the film. The indentation depth at which the substrate starts to deform plastically is sensitive to indenter geometry.     

The indenter tip radius effect in micro- and nanoindentation hardness experiments.

Nix and Gao established an important relation between the microindentation hardness and indentation depth. Such a relation has been verified by many microindentation experiments (indentation depths in the micrometer range), but it does not always hold in nanoindentation experiments (indentation depths approaching the nanometer range). Indenter tip radius effect has been proposed by Qu et al. and others as possibly the main factor that causes the deviation from Nix and Gao's relationship. We have developed an indentation model for micro- and nanoindentation, which accounts for two indenter shapes, a sharp, conical indenter and a conical indenter with a spherical tip. The analysis is based on the conventional theory of mechanism-based strain gradient plasticity established from the Taylor dislocation model to account for the effect of geometrically necessary dislocations. The comparison between numerical result and Feng and Nix's experimental data shows that the indenter tip radius effect indeed causes the deviation from Nix-Gao relation, but it seems not be the main factor. The project supported by the National Natural Science Foundation of China (10121202) and the Ministry of Education of China (20020003023)     

The correlation of the indentation size effect measured with indenters of various shapes

Experimental results are presented which show that the indentation size effect for pyramidal and spherical indenters can be correlated. For a pyramidal indenter, the hardness measured in crystalline materials usually increases with decreasing depth of penetration, which is known as the indentation size effect. Spherical indentation also shows an indentation size effect. However, for a spherical indenter, hardness is not affected by depth, but increases with decreasing sphere radius. The correlation for pyramidal and spherical indenter shapes is based on geometrically necessary dislocations and work-hardening. The Nix and Gao indentation size effect model (J. Mech. Phys. Solids 46 (1998) 411) for conical indenters is extended to indenters of various shapes and compared to the experimental results.     

A study on robust indentation techniques to evaluate elastic–plastic properties of metals

Spherical indentation is studied based on numerical analysis and experiment, to develop robust testing techniques to evaluate isotropic elastic–plastic material properties of metals. The representative stress and plastic strain concept is critically investigated via finite element analysis, and some conditions for the representative values are suggested. The representative values should also be a function of material properties, not only indenter angle for sharp indenter and indentation depth for spherical indenter. The pros and cons of shallow and deep spherical indentation techniques are also discussed. For an indentation depth of 20% of an indenter diameter, the relationships between normalized indentation parameters and load–depth data are characterized, and then numerical algorithm to estimate material elastic–plastic curve is presented. From the indentation load–depth curve, the new approach provides stress–strain curve and the values of elastic modulus, yield strength, and strain-hardening exponent with an average error of less than 5%. The method is confirmed to be valid for various elastic properties of indenter. Experimental validation of the approach then is performed by using developed micro-indentation system. For the material severely disobeying power law hardening, a modified method to reduce errors of predicted material properties is contrived. It is found that our method is robust enough to get ideal power law properties, and applicable to input of more complex physics.     

金属材料的强度与应力-应变关系的球压入测试方法

压入法获取材料单轴应力-应变关系和抗拉强度对服役结构完整性评价有重要的基础意义.假定材料均匀连续、各向同性、应力应变关系符合Hollomon律,基于能量等效假定,即代表性体积单元(representativevolume element, RVE)的vonMises等效和有效变形域内能量中值等效假定,本文提出了关联材料载荷、深度、球压头直径和Hollomon律的四参数半解析球压入(semi-analyticalspherical indentation,SSI)模型.通过球压入载荷-深度试验关系获得材料的应力-应变关系和抗拉强度.考虑压入过程中的损伤效应,针对金属材料提出了用于球压入测试的材料弹性模量修正模型.对11种延性金属材料完成了球压入试验,采用本文提出的球压入试验方法测到的弹性模量、应力-应变关系和抗拉强度与单轴拉伸试验结果吻合良好.      

Size effects in the conical indentation of an elasto-plastic solid

The size effect in conical indentation of an elasto-plastic solid is predicted via the Fleck and Willis formulation of strain gradient plasticity (Fleck, N.A. and Willis, J.R., 2009, A mathematical basis for strain gradient plasticity theory. Part II: tensorial plastic multiplier, J. Mech. Phys. Solids, 57, 1045–1057). The rate-dependent formulation is implemented numerically and the full-field indentation problem is analyzed via finite element calculations, for both ideally plastic behavior and dissipative hardening. The isotropic strain-gradient theory involves three material length scales, and the relative significance of these length scales upon the degree of size effect is assessed. Indentation maps are generated to summarize the sensitivity of indentation hardness to indent size, indenter geometry and material properties (such as yield strain and strain hardening index). The finite element model is also used to evaluate the pertinence of the Johnson cavity expansion model and of the Nix–Gao model, which have been extensively used to predict size effects in indentation hardness.     

Strain gradient plasticity solution for an internally pressurized thick-walled spherical shell of an elastic–plastic material

An analytical solution is presented for an internally pressurized thick-walled spherical shell of an elastic strain-hardening plastic material. A strain gradient plasticity theory is used to describe the constitutive behavior of the material undergoing plastic deformations, whereas the generalized Hooke’s law is invoked to represent the material response in the elastic region. The solution gives explicit expressions for the stress, strain and displacement components. The inner radius of the shell enters these expressions not only in non-dimensional forms but also with its own dimensional identity, unlike classical plasticity-based solutions. As a result, the current solution can capture the size effect. The classical plasticity-based solution of the same problem is shown to be a special case of the present solution. Numerical results for the maximum effective stress in the shell wall are also provided to illustrate applications of the newly derived solution. The new solution can be used to construct improved expanding cavity models in indentation mechanics that incorporate both the strain-hardening and indentation size effects.     

基于应变率阶跃测试的晶体铜压入应变率敏感性研究

纳米压入测试可以原位获取材料的诸多力学性能,包括弹性模量,硬度,屈服应力,应变率敏感指数等。本文利用应变率阶跃测试技术对多晶铜试样的应变率敏感性进行测试分析,硬度-位移曲线表明压头下方所存在的变形梯度对各阶跃应变率下的硬度值存在明显影响;采用基于晶体细观机制的塑性应变梯度理论对压入变形梯度效应予以修正,比较了修正与未修正数据所得的应变率敏感指数,在有效剔除压入变形梯度影响的基础上,应变率阶跃测试可实现单次压入下材料应变率敏感性的测试表征。     

Evaluation of the stress–strain curve of metallic materials by spherical indentation

A method for deducing the stress–strain uniaxial properties of metallic materials from instrumented spherical indentation is presented along with an experimental verification.An extensive finite element parametric analysis of the spherical indentation was performed in order to generate a database of load vs. depth of penetration curves for classes of materials selected in order to represent the metals commonly employed in structural applications. The stress–strain curves of the materials were represented with three parameters: the Young modulus for the elastic regime, the stress of proportionality limit and the strain-hardening coefficient for the elastic–plastic regime.The indentation curves simulated by the finite element analyses were fitted in order to obtain a continuous function which can produce accurate load vs. depth curves for any combination of the constitutive elastic–plastic parameters. On the basis of this continuous function, an optimization algorithm was then employed to deduce the material elastic–plastic parameters and the related stress–strain curve when the measured load vs. depth curve is available by an instrumented spherical indentation test.The proposed method was verified by comparing the predicted stress–strain curves with those directly measured for several metallic alloys having different mechanical properties.This result confirms the possibility to deduce the complete stress–strain curve of a metal alloy with good accuracy by a properly conducted instrumented spherical indentation test and a suitable interpretation technique of the measured quantities.     

A numerical approach to indentation technique to evaluate material properties of film-on-substrate systems

Spherical indentation approach (Lee et al., 2005, Lee et al., 2010) for the evaluation of bulk material properties is extended to that for elastic–plastic properties of film-on-substrate systems. Our interest focuses on single isotropic, metallic, and elastic–plastic film on a substrate, and we do not consider the size effects in plasticity behavior. We first determine the optimal data acquisition location, where the strain gradient is the least and the effect of friction is negligible. Dimensional analysis affords the mapping parameters as functions of normalized indentation variables. An efficient way is further introduced to reduce both the number of analyses and the regression order of mapping functions. The new numerical approach to the film indentation technique is then proposed by examining the finite element solutions at the optimal point. With the new approach, the values of elastic modulus, yield strength, and strain-hardening exponent of film materials are successfully obtained from the spherical indentation tests. We have shown that the effective property ranges such as indenter properties, substrate modulus, and E/E_s ratio can be extended without additional simulations and even loss of accuracy. For other ranges of variables or other properties, which are not dealt with in this study, this methodology is applicable through resetting FEA variables and finding proper normalized parameters.     

On establishing elastic–plastic properties of a sphere by indentation testing

Instrumented indentation is a popular technique for determining mechanical properties of materials. Currently, the evaluation techniques of instrumented indentation are mostly limited to a flat substrate being indented by various shaped indenters (e.g., conical or spherical). This work investigates the possibility of extending instrumented indentation to non-flat surfaces. To this end, conical indentation of a sphere is investigated where two methodologies for establishing mechanical properties are explored. In the first approach, a semi-analytical approach is employed to determine the elastic modulus of the sphere utilizing the elastic unloading response (the “unloading slope”). In the second method, reverse analysis based on finite element analysis is used, where non-dimensional characteristic functions derived from the force–displacement response are utilized to determine the elastic modulus and yield strength. To investigate the accuracies of the proposed methodologies, selected numerical experiments have been performed and excellent agreement was obtained.     

Mechanics of indentation of plastically graded materials—I: Analysis

The introduction of controlled gradients in plastic properties is known to influence the resistance to damage and cracking at contact surfaces in many tribological applications. In order to assess potentially beneficial effects of plastic property gradients in tribological applications, it is essential first to develop a comprehensive and quantitative understanding of the effects of yield strength and strain hardening exponent on contact deformation under the most fundamental contact condition: normal indentation. To date, however, systematic and quantitative studies of plasticity gradient effects on indentation response have not been completed. A comprehensive parametric study of the mechanics of normal indentation of plastically graded materials was therefore undertaken in this work by recourse to finite element method (FEM) computations. On the basis of a large number of computational simulations, a general methodology for assessing instrumented indentation response of plastically graded materials is formulated so that quantitative interpretations of depth-sensing indentation experiments could be performed. The specific case of linear variation in yield strength with depth below the indented surface is explored in detail. Universal dimensionless functions are extracted from FEM simulations so as to predict the indentation load versus depth of penetration curves for a wide variety of plastically graded engineering metals and alloys for interpretation of, and comparisons with, experimental results. Furthermore, the effect of plasticity gradient on the residual indentation pile-up profile is systematically studied. The computations reveal that pile-up of the graded alloy around the indenter, for indentation with increasing yield strength beneath the surface, is noticeably higher than that for the two homogeneous reference materials that constitute the bounding conditions for the graded material. Pile-up is also found to be an increasing function of yield strength gradient and a decreasing function of frictional coefficient. The stress and plastic strain distributions under the indenter tip with and without plasticity gradient are also examined to rationalize the predicted trends. In Part II of this paper, we compare the predictions of depth-sensing indentation and pile-up response with experiments on a specially made, graded model Ni-W alloy with controlled gradients in nanocrystalline grain size.     

Polymer indentation: Numerical analysis and comparison with a spherical cavity model

Three dimensional analyses of indentation of a polymer by a rigid indenter are carried out. The polymer is characterized by a finite strain elastic-viscoplastic constitutive relation and the calculations are carried out with a dynamic finite element program used to simulate quasi-static behavior. Two types of indenter are considered; a conical indenter and a pyramidal indenter. For each indenter type, different values of the sharpness of the indenter are considered and two rates of indentation are analyzed. Significant sink-in is found to occur in all cases considered. The amount of sink-in is found to be smaller for sharper indenters. The calculated values of both the nominal and true hardness do not differ significantly for the two indenter shapes. An expanding spherical cavity model is also considered and the predictions of this model are compared with the finite element results for various indenter shapes and indentation rates. The spherical cavity model is found to give fairly good agreement with the predictions of the finite element analyses for the nominal polymer hardness for both indenter shapes.     

Simulation of micro-indentation hardness of FCC single crystals by mechanism-based strain gradient crystal plasticity

The size effect observed in the micro-indentation of FCC single crystal copper is modelled by the employment of mechanism-based strain gradient crystal plasticity (MSG-CP). The total slip resistance in each active system is assumed to be due to a mixed population of forest obstacles arising from both statistically stored and geometrically necessary dislocations. The MSG-CP constitutive model is implemented into the Abaqus/Standard FE platform by developing the User MATerial subroutine UMAT. The simulation of micro-indentation hardness on (0 0 1) and (1 1 1) single crystal copper, with a conical indenter having a sharp tip, a conical indenter with a spherical tip and a three-sided Berkovich indenter, is undertaken. The phenomena of pile-up and sink-in have been observed in the simulation and dealt with by appropriate use of the contact analysis function in Abaqus. These phenomena have been taken into account in the determination of the contact areas and hence the average indentation depth for anisotropic single crystals. The depth dependence of the micro-indentation hardness, the size effect, is calculated. The micro-hardness results from the simulation are compared with those of the published experimental ones in the literature and a good agreement is found.     

A computational study on the instrumented sharp indentations with dual indenters

Finite element analysis was performed to investigate the indentation response of elasto-plastic solids for conical indenters of half included angles of 60° and 70.3°. The interdependence indentation parameters resulting from a single indentation load–depth curve is considered. Regarding dimensional analysis, several dimensionless relationships are constructed as functions of the reduced elastic modulus-loading curvature ratio E¹/C and the strain hardening exponent n. Further, the duality between corresponding parameters with dual indenters is explored. Finally, a new method based on dual indenters is proposed to extract the strain hardening exponent and the reduced elastic modulus of an indented material. The accuracy of this method is verified and discussed with experimental data from the literature and representative materials.     

考虑压头曲率半径和应变梯度的微压痕分析

在压头尖端曲率半径取100nm的前提下，采用Chen和Wang的应变梯度理论，对微压痕实验进行了系统的数值分析. 首先通过拟合载荷-位移实验曲线的后半段来确定材料的屈服应力和幂硬化指数值,然后用有限元方法数值模拟压痕实验，并将计算得到的整段载荷-位移曲线及硬度-位移曲线和实验结果进行了比较. 结果表明应变梯度理论所预测的计算结果和实验结果很好地符合,包括压痕深度在亚微米和微米范围内的整段曲线.     

考虑压头曲率半径的J2形变理论压痕有限元分析

采用经典J2形变理论,在考虑压头曲率半径的前提下,针对几组微压痕实验进行了有限元数值模拟,给出理论计算结果并与实验进行了系统比较．结果发现考虑了压头曲率半径后的经典J2形变理论得到的整段计算曲线不能和实验曲线吻合,可以推理出即使考虑压头曲率半径的的影响,经典J2形变理论也不能很好地解释实验现象．     

Mechanics of indentation for piezoelectric thin films on elastic substrate

Frictionless normal indentation problem of rigid flat-ended cylindrical, conical and spherical indenters on piezoelectric film, which is either in frictionless contact with or perfectly bonded to an elastic half-space (substrate), is investigated. Both conducting and insulating indenters are considered. With Hankel transform, the general solutions of the homogeneous governing equations for the piezoelectric layer and the elastic half-space are presented. Using the boundary conditions for a vertical point force or a point electric charge, and the boundary conditions on the film/substrate interface, the Green’s functions can be obtained by solving sets of simultaneous linear algebraic equations. The solution of the indentation problem is obtained by integrating these Green’s functions over the contact area with unknown surface tractions or electric charge distribution, which will be determined from the boundary conditions on the contact surface between the indenter and the film. The solution is expressed in terms of dual integral equations that are converted to a Fredholm integral equation of the second kind and solved numerically. Numerical examples are also presented. The comparison between two film/substrate bonding conditions is made. It shows that the indentation rigidity of the film/substrate system is lower when the film is in frictionless contact with the substrate. The effects of the Young’s modulus and Poisson’s ratio of the elastic substrate, indenter electrical condition and indenter prescribed electric potential on the indentation responses are presented.