# Examining the Effect of Pileup on the Accuracy of Sharp Indentation Testing

Authors: Zhongbin Tang; Yazhou Guo; Zhaohu Jia; Yulong Li; Qiuming Wei
Journal: Advances in Materials Science and Engineering
 (2015)

---

## Abstract

One essential problem related to instrumented indentation is the effect of pileup, which could introduce significant errors on the measured hardness and elastic modulus. In this work, we have assessed some critical issues associated with instrumented indentation by means of numerical simulation. Dimensional analysis is adopted to acquire the principal governing parameters of the process, such as the ratio of yield strength to elastic modulus(Y/E) and the ratio of indentation depth to maximum penetration depth (h/hmax). The dimensional analysis provides a general understanding on the relationships among these dimensionless variables and their influence on the test results. Three kinds of indenters, that is, conical, Berkovich, and Vickers indenter, are comparatively studied. The dependence of pileup on the dimensionless parameters and its effect on the simulation results are also given. These findings are instructive for measuring the mechanical properties by instrumented indentation, especially when the effect of pileup is involved.


---

## Body

## 1. Introduction

Instrumented indentation, particularly nanoindentation, has been widely used for characterizing the mechanical properties of materials at small scales [1–6]. One of the great advantages of this technique is that the mechanical properties, such as hardness and elastic modulus, can be directly determined from the indentation load-displacement curves, while the measurement of the actual hardness impression is not necessary. Recent developments of the nanoindentation technique make it capable of recording experimental data at load and displacement as small as micro-Newton and nanometer levels.The commonly used definition of hardness is as follows:(1)H=PA,where H is hardness, P the load applied to the test surface, and A the projected contact area. A was calculated by some researchers using the following equation [1, 4]:(2)A=f1hc=f2h,where hc is the contact depth and h is the resulting penetration depth. If pileup effect is negligible, hc could be given by(3)hc=h-εPS,where ε is a constant that depends on indenter geometry and S is the elastic stiffness of the contact [4]. One should notice that (3) cannot account for the plastic phenomenon of pileup, which may lead to significant change of the true contact area.The elastic unloading curve of nanoindentation is approximately expressed in the power-law relation [4, 7]:(4)PPmax⁡=αh-hfhmax⁡-hfn,where α and n are empirically determined fitting parameters, depending on the geometry of the indenter. hf is the final displacement of the indenter after complete unloading. Then S could be acquired by calculating the slope of the initial portion of the unloading curve:(5)S=dPdhh=hmax⁡=αPmax⁡nhmax⁡-hf-1.The relationship between the effective elastic modulus and the projected contact area and unloading stiffness is given by [1, 4](6)Eeff=π×S2βA,where β is a dimensionless factor that accounts for deviations in stiffness caused by the lack of axial symmetry for Berkovich or Vickers indenters. The effective elastic modulus is defined by (7)1Eeff=1-ν2E+1-νi2Ei,where E and ν are elastic modulus and Poisson’s ratio of the specimen and Ei and νi are those of the indenter. The elastic modulus of the specimen, E, could then be calculated by (7). In this work, the indenter is taken as rigid body for numerical simulation, and then E=Eeff(1-ν2).From (1) to (7), the hardness as well as the elastic modulus of the material could be determined in the following sequence:(a)
Fit the unloading curve of the indentation; find values ofβ and n, and determine S by (5).(b)
Find contact depthhc by (3).(c)
Determine the projected contact areaA by (2).(d)
Obtain hardness by (1).(e)
DeriveEeff by (6) and then elastic modulus E by (7).The difficulty related to the above method is the determination of the contact depth hc and the true projected contact area A, especially when the phenomenon of pileup occurs. The resulting contact area might be greatly underestimated by using these equations if pileup effect exists. Bolshakov and Pharr [7] found by numerical simulation that pileup was significant for conical indenter if the ratio of final displacement to the maximum depth of penetration is larger than 0.7 (i.e., hf/hmax>0.7). Y.-T Cheng and C.-M Cheng [2, 8–11] did a series of analysis on conical indentation by finite element method (FEM). They analyzed their data by the method of dimensional analysis, which provided us a general understanding of this process. The dependencies of H and E on some critical dimensionless variables, such as Y/E, ν, and θ (Y is the yield strength of the tested material and θ is the half angle of the indenter), were investigated. However, their numerical simulations were limited to 2D, which may represent the conical indentation fairly well but are not so representative for some 3D indenters such as Berkovich indenter and Vickers indenter. Dao et al. [12] conducted comprehensive analysis on the data derived from indentation by finite element simulation, where they also used 2D axial symmetric model instead of the real configuration. Giannakopoulos et al. [13] and Larsson et al. [14] investigated Vickers indenter and Berkovich indenter separately, but detailed comparison among these sharp indenters, such as Vickers, Berkovich, and conical, was not provided. Although the geometry functions of Berkovich indenter and Vickers indenter can be designed to be the same as that of the conical indenter, the true projected contact areas at the same indentation depth may still be different considering the effect of pileup.Apart from geometric difference of indenters, there are several other factors that may influence the measurement of instrumented indentation. The thickness and properties of the test materials and the substrate effect have been studied systematically by many researchers [16–21]. The stress state within the material also affects the testing results of indentation [22, 23]. In fact, Larsson et al. successfully related stress state to indentation testing results and found a way of determining residual stress within the material [24, 25] by sharp indentation.In this paper, we will focus on the differences of indenter geometries and their effect on the indentation results by means of numerical simulation. Dimensional analysis will be adopted to acquire the principal governing parameters. The effect of pileup on the simulation results will also be given.Three kinds of indenters, namely, conical indenter, Berkovich indenter, and Vickers indenter, will be discussed separately. The half angle of the conical indenter is set to be 70.3°, while the half angles of Berkovich indenter and Vickers indenter are chosen to be 65.3° and 68°, respectively. The relationship between the projected contact areaA and the contact depth hc (i.e., f1(hc)) can be obtained through simple calculation and all of the three indenters have equivalent geometry functions:(8)A=f1hc≈24.52hc2.


## 2. Dimensional Analysis

### 2.1. Dimensional Analysis for Loading

For a simple instrumented indentation test, the procedure is composed of two steps, that is, loading and unloading. During the loading process, the produced loadP and projected contact area A are functions of all the independent governing parameters, which include the material parameters (Young’s modulus E, yield strength Y, and Poisson’s ratio ν), the indenter geometry parameter (those in (2)), and the indentation depth h. For a given indenter, P and A can be written as(9)P=fpLE,Y,ν,h,A=fAE,Y,ν,h.Among these governing parameters, we may choose two basic dimensions, that is, [E] and [h]. The dimensions of the rest of the variables can then be described as(10)P=Eh2,A=h2,Y=E,ν=E0h0.Since the influence of ν on the measurement of hardness and elastic modulus is relatively small [1] (e.g., ν=0.25±0.1 produces only about 5% uncertainty in the calculated value of E), we will not incorporate the effect of ν for simplicity in this work. Conducting dimensional analysis by applying the Π theorem [26], the following equations can be obtained:(11)ψΠΠp1,Πp2=0,φΠΠA1,ΠA2=0,where ΠP1=PE-1h-2, ΠP2=ΠA2=YE-1, and ΠA1=Ah-2. Then the above equations can be rewritten as(12)P=Eh2ϕPLYE,(13)A=h2ϕAYE.Note that ϕPL and ϕA are different functions from ψΠ and φΠ, respectively. It is easy to see from (12) and (13) that P and A are all proportional to h2. The hardness can then be calculated by(14)H=PmaxA=EϕHYE.Clearly, the hardness depends only on the material parameters E and Y, and it is independent from the indentation depth. Moreover, the dimensionless hardness, that is, H/E, is a function of only one parameter Y/E from (14).


### 2.2. Dimensional Analysis for Unloading

During unloading, the imposed loadP depends on the following independent variables: E, Y, h, and maximum indentation depth hmax, while the final indentation depth hf depends only on E, Y, and hmax (the effect of ν is not considered here). We have(15)P=fPUE,Y,h,hmax,hf=fhfE,Y,hmax.Again, the dimensional analysis gives(16)P=Eh2ϕPUYE,hhmax,(17)hf=hmaxϕhfYE.Equation (17) indicates that the final indentation depth is proportional to the maximum indentation depth and is dependent on Y/E. Since the slope of the unloading curve is essential for calculating E and H, we take derivative with respect to h on both sides of (16) and obtain the values at h=hmax:(18)dPdhh=hmax=Ehmax2ϕPUYE,1+ϕPU˙YE,1.Or equivalently, (19)dPdhh=hmax=S=2ϕPUYE,1+ϕPU˙YE,1=EhmaxϕSYE.The dimensionless form of the elastic contact stiffness S/(Ehmax) depends on Y/E.


## 2.1. Dimensional Analysis for Loading

For a simple instrumented indentation test, the procedure is composed of two steps, that is, loading and unloading. During the loading process, the produced loadP and projected contact area A are functions of all the independent governing parameters, which include the material parameters (Young’s modulus E, yield strength Y, and Poisson’s ratio ν), the indenter geometry parameter (those in (2)), and the indentation depth h. For a given indenter, P and A can be written as(9)P=fpLE,Y,ν,h,A=fAE,Y,ν,h.Among these governing parameters, we may choose two basic dimensions, that is, [E] and [h]. The dimensions of the rest of the variables can then be described as(10)P=Eh2,A=h2,Y=E,ν=E0h0.Since the influence of ν on the measurement of hardness and elastic modulus is relatively small [1] (e.g., ν=0.25±0.1 produces only about 5% uncertainty in the calculated value of E), we will not incorporate the effect of ν for simplicity in this work. Conducting dimensional analysis by applying the Π theorem [26], the following equations can be obtained:(11)ψΠΠp1,Πp2=0,φΠΠA1,ΠA2=0,where ΠP1=PE-1h-2, ΠP2=ΠA2=YE-1, and ΠA1=Ah-2. Then the above equations can be rewritten as(12)P=Eh2ϕPLYE,(13)A=h2ϕAYE.Note that ϕPL and ϕA are different functions from ψΠ and φΠ, respectively. It is easy to see from (12) and (13) that P and A are all proportional to h2. The hardness can then be calculated by(14)H=PmaxA=EϕHYE.Clearly, the hardness depends only on the material parameters E and Y, and it is independent from the indentation depth. Moreover, the dimensionless hardness, that is, H/E, is a function of only one parameter Y/E from (14).


## 2.2. Dimensional Analysis for Unloading

During unloading, the imposed loadP depends on the following independent variables: E, Y, h, and maximum indentation depth hmax, while the final indentation depth hf depends only on E, Y, and hmax (the effect of ν is not considered here). We have(15)P=fPUE,Y,h,hmax,hf=fhfE,Y,hmax.Again, the dimensional analysis gives(16)P=Eh2ϕPUYE,hhmax,(17)hf=hmaxϕhfYE.Equation (17) indicates that the final indentation depth is proportional to the maximum indentation depth and is dependent on Y/E. Since the slope of the unloading curve is essential for calculating E and H, we take derivative with respect to h on both sides of (16) and obtain the values at h=hmax:(18)dPdhh=hmax=Ehmax2ϕPUYE,1+ϕPU˙YE,1.Or equivalently, (19)dPdhh=hmax=S=2ϕPUYE,1+ϕPU˙YE,1=EhmaxϕSYE.The dimensionless form of the elastic contact stiffness S/(Ehmax) depends on Y/E.


## 3. Numerical Analysis

Finite element analysis (FEA) of the indentation process has been carried out by using the commercial software Abaqus. The above-mentioned scaling relationships are carefully examined. These dimensionless functions, such asϕPL, ϕA, ϕhf, and ϕS, are evaluated. The effect of pileup on the calculated hardness and elastic modulus will also be discussed.As mentioned before, three kinds of indention, namely, conical, Vickers, and Berkovich, are simulated. All of the indenters have the same projected area function (8). The dimension of the tested specimen was set to be larger than ten times of the contact length in every direction to minimize the influence of the boundary. Similar to experimental test, the bottom surface of the specimen is set to be fixed and all of the remaining surfaces are set to be free. The indenter is free along the indentation direction only and the loads are given through applying displacements on the indenter. The contact between the indenter and specimen is set to be frictionless and the whole process is considered to be time-independent.The indenters were taken as rigid body, while the tested materials were regarded as elastic-perfect-plastic. The yield strengths of the materials are set to range from 80 MPa to 8 GPa, which covers most of the engineering materials. Two elastic modulus values, that is, 100 GPa and 200 GPa, are considered. Poisson’s ratio is set to be 0.38, corresponding to a material of niobium that will be used for comparison with the experimental observations. Actually, the values of Poisson’s ratio taken from 0.2 to 0.4 will not lead to any significant difference [8].A total number of 60,000 brick elements are meshed in the material and about 3,000 quadrilateral elements are meshed in the indenter. An example of the mesh of half of the conical indentation model is given in Figure1. Intense mesh is given in the region of contact and sparse mesh in the far field to compromise between accuracy and efficiency. Mesh dependence of the simulation results was conducted first and further increase of element number led to little change of theP-h curves.Figure 1
Finite element meshes of the numerical model of conical indentation (only half model is shown for better view). The image in the red box provides the details of the mesh in the region of contact.Similar to experiments, theP-h curves are extracted directly from the calculated results of each simulation. Other properties such as hardness and elastic modulus are derived from theP-h curves following the procedure given in Section 1 of this paper. The errors related to FEA include modeling error, discretization error, and numerical error. Since all of the simulation models are identical in the sense of modeling, loading, meshing, contact condition, and solver and the only differences are the geometry of indenters and the properties of testing materials, the simulated results can accurately reflect the role of these different inputs in instrumented indentation. Actually, one advantage of the “numerical experiment” is that it excludes all other factors that may influence the results and focuses on one or two concerned terms, which is very hard to achieve in real experiments.


## 4. Results and Discussion

Figure2 displays the displacement field in the loading direction at a penetration depth of 3 μm for the three indentations. The current Y and E are 800 MPa and 100 GPa, respectively. Only half model is shown for better view. Note from the legend that positive values of displacement (U3) represent penetration into the material, while negative values indicate pileup. Clearly, pileup exists in each of the situations, demonstrating that corrections are required when using traditional process equations (2) and (3). Moreover, although the geometry functions of the three indenters are the same, their pileup patterns are rather different. These differences indicate that axial symmetric conical indenter may not be representative for the 3-side Berkovich indenter or 4-side Vickers indenter, especially when pileup plays a role.Figure 2
Displacement fields at indentation depth of 3μm for conical indenter (a), Vickers indenter (b), and Berkovich indenter (c). Only half of each model is shown. Note that the negative values of displacement (U3) indicate pileup of the material.
(a)
(b)
(c)The evolution of forceP with respect to indentation depth h is plotted in Figure 3. The loading part of each indentation curve is fitted by the power law function. The fitting exponents, b, for all the simulations are between 1.93 and 2.02, corresponding well with the prediction by (12). The true contact area, Ac, between the indenter and the material is easy to obtain from the simulation. The true projected contact area, A, can be derived from Ac through simple geometric relationships. Figure 4 displays theA-h curves for the three kinds of indentations. Although they behave differently in Figure 2, the three kinds of indenters present identicalA-h curves given the same Y and E. The loading part of theA-h curves is also fitted by power law function and it gives the values of b ranging from 2.0 to 2.1, in accordance with the prediction by (13). A reference area Ar, defined as the projected area calculated from indentation depth h by (8) is given in Figure 4. Obviously, Ar is significantly smaller than the true projected areas when the indentation depth is larger than 1 μm, signifying essential error may be introduced if the effect of pileup is ignored.Figure 3
Force-displacement curves for different kinds of indenters. The solid lines are the power law fitting of the loading parts.E=100 GPa, Y=800 MPa.Figure 4
True projected area-displacement curves for three kinds of indenters.Ar is the reference area that is calculated from the indentation depth by using the area function of the indenters. The solid lines are the power law fitting of the loading parts. E=100 GPa, Y=800 MPa.The dimensionless functions, that is,ϕPL=P/(Eh2), ϕA=A/h2, ϕhf=hf/hmax, and ϕS=S/(Eh), were plotted with respect to Y/E in Figures 5, 6, 7, and 8, respectively. It is obvious from Figure 5 that these derived quantities lie in a single curve for each indenter. Thus ϕPL is a function of Y/E, as predicted by (12). Nevertheless, these functions are not exactly the same for each indenter. When Y/E is larger than 0.01 (corresponding to materials like soda-lime glass, high strength steels, titanium alloys, and so on), the differences among the data derived from the three indenters become evident. Similarly, the dimensionless function for projected contact area ϕA is a function of Y/E and is different for each indenter. But, different from ϕPL, the discrepancy of ϕA becomes significant when Y/E is smaller than 0.02 (corresponding to materials like aluminum alloys, low carbon steels, and so on). The dimensionless function ϕhf decreases monotonously with increasing Y/E and appears independent of the geometry of the indenter (see Figure 7). hf/hmax is a very useful quantity, because it is easy to obtain from theP-h curve and many researchers [7–9, 14] took it as an indicator for pileup or sink-in. The dimensionless contact stiffness ϕS becomes smaller as Y/E increases but the differences among different indenters become more significant. Since the contact stiffness is essential for the derivation of the elastic modulus (see (6)), its scattering among different indenters may lead to significant variation of the calculated elastic modulus.Figure 5
The variation of dimensionless functionϕPL=P/(Eh2) with respect to Y/E for different kinds of indenters.Figure 6
The variation of dimensionless functionϕA=A/h2 with respect to Y/E for different kinds of indenters.Figure 7
The variation of dimensionless functionϕhf=hf/hmax with respect to Y/E for different kinds of indenters.Figure 8
The variation of dimensionless functionϕS=S/(Eh) with respect to Y/E for different kinds of indenters.The projected contact areaA is an important quantity that directly influences the measured hardness. The dependence of the normalized projected contact area A/Ar on hf/hmax is given in Figure 9. The true projected contact areas (solid symbols) and those calculated by the Oliver and Pharr (O/P) model [1] (hollow symbols) are all presented. The red dashed line, where A/Ar is 1.0, denotes no pileup or sink-in. All of the hollow symbols are below the dashed line, showing that the O/P model considers sink-in only and could lead to significant error if hf/hmax is larger than 0.8. The solid symbols intersect with the dashed line at around hf/hmax=0.85 for all the three indenters, indicating that direct penetration depth and area function could be used to calculate projected contact area A at this point (8). Our simulation results also reveal that less than 10% error will be introduced if pileup or sink-in effect is ignored when 0.8<hf/hmax<0.9 (corresponding 0.008<Y/E<0.02) for Berkovich indentation. However, this error increases to more than 20% for conical or Vickers indentation. This observation demonstrates that Berkovich indenter is less sensitive to pileup or sink-in than Vickers and conical indenters within this region. If pileup and sink-in are ignored, the maximum error could be larger than 50% for the conical indenter.Figure 9
Influence ofhf/hmax on the projected contact area. The projected contact area is normalized by Ar, the area calculated by the indenter geometry function. The solid icons stand for the true projected contact area derived from finite element models, while the hollow icons represent the projected contact area calculated by Oliver and Pharr model [1].The derived elastic modulus of the material depends on the contact stiffnessS, the projected contact area A, and the correction factor β (6). The differences of S and A could lead to variation of E for different indenters. Figure 10 displays the dependence of calculated E on hf/hmax. The solid symbols represent E derived from the true projected contact area, while hollow symbols represent E derived from the O/P model. Obviously, all of the values of calculated E are larger than the input value and the O/P method could introduce more error if improper correction factor β is used. Figure 10 shows that the error of calculated E increases with increment of hf/hmax, demonstrating that the correction factor β should not be a constant but a function of hf/hmax; that is, β=a+b(hf/hmax). For the non-work-hardening materials in this work, a and b are 1.021 and 0.358, 0.832 and 0.526, and 0.759 and 0.512 for the Berkovich, Vickers, and conical indenters, respectively. Another noteworthy feature in Figure 10 is that Berkovich indentation will yield higher elastic modulus than Vickers and conical indentations, indicating again that replacing Berkovich and Vickers by conical indenter is not appropriate. An approximate relationship of E among different indenters can be given as EBerkovich=1.06EVickers=1.13EConical. This estimation leads to less than 5% error when hf/hmax<0.95 (corresponding to Y/E>0.002).Figure 10
Influence ofhf/hmax on the calculated elastic modulus. The modulus is normalized by the input E of the numerical model. The solid icons stand for the elastic modulus derived from true projected contact areas (6), while the hollow icons represent the modulus derived by Oliver and Pharr model [1]. The solid and dashed lines are linear fitting of solid icons and hollow icons, respectively.Figure11 illustrates the variation of the normalized hardness HO/P/Htrue with respect to hf/hmax. HO/P and Htrue are derived by using the O/P model and the true projected contact area, respectively. The O/P hardness, as expected, is larger than the true hardness for all cases. When hf/hmax<0.8, HO/P/Htrue≈1.1 and this approximation will introduce no more than 5% error for all the three indenters. However, this ratio increases sharply when hf/hmax>0.8. For conical indentation, HO/P/Htrue rises monotonously up to 1.54, but the maximum values for Berkovich and Vickers indentation are less than 1.3. Thus, conical indentation fails again to represent Berkovich and Vickers indentations.Figure 11
Effect ofhf/hmax on HO/P/Htrue, the ratio of hardness derived by Oliver and Pharr model to that calculated by using true projected contact area.The hardness is normalized by yield stress of the material and plotted againsthf/hmax in Figure 12. In the case that H is calculated by the O/P method, similar trends of H/Y are observed for all the three indenters in Figure 12, although there are still differences in the magnitudes. When hf/hmax approaches unity, the maximum values of H/Y are achieved and are larger than 4.2 for all indenters. When true hardness is applied, different patterns of H/Y are found for different indenters. The relationship between hardness and yield stress, C=H/Y, was firstly proposed by Tabor [27], where C is the constraint factor. Tabor analyzed the rigid/plastic deformation of Vickers indentation and claimed that the hardness of metals is about 2.9~3.0 times its flow stress at 8~10% strain; that is, C≈2.9~3.0. Since rigid/plastic deformation corresponds to hf/hmax=1, our numerical results could be compared with C. Our simulation shows that H/Y≈3.3 for Vickers indentation when hf/hmax approaches unity, slightly larger than Tabor’s value. Different from Vickers and Berkovich, the value of H/Y for conical indentation plateaus at about 2.7, which is close to that calculated by Bolshakov and Pharr [7]. Clearly, C is not a constant but varies with hf/hmax. The value of C also depends on the indenter geometry. For soft material, that is, hf/hmax close to 1 (or Y/E close to zero from Figure 7), C reaches 2.7, 3.3, and 3.6 for conical, Vickers, and Berkovich indentation, respectively. However, C might be lower than 2.0 for some hard materials.Figure 12
Influence ofhf/hmax on H/Y. The solid icons stand for the hardness derived from true projected contact areas, while the hollow icons represent those derived by Oliver and Pharr model [1].It should be noted thathf/hmax=0.8 (or Y/E=0.02) is a critical point, below which the true projected contact area, hardness, and Tabor factor are almost proportional to the values derived from Oliver and Pharr model. In that case, the results by traditional method could still be used if a correction factor is applied (hardness divided by, e.g., 1.1). However, this factor could not be used if hf/hmax>0.8 (or Y/E<0.02), where the relation between true values and Oliver-Pharr results becomes different for different indenters. As Y/E is less than 0.02 for most of the engineering materials, the corrections should be more careful.Our experimental results [15] on niobium could be used to compare with some of the simulation results. The commercially pure niobium disks with diameter of 10.0 mm and thickness of 1.0 mm were processed by high pressure torsion (HPT) at room temperature. The Von Mises equivalent strain at the edge of the disk reaches about 90, while it is theoretically zero at the disk center. The large deformation by HPT substantially enhances the yield strength of the material while reducing its ductility. One distinct characteristic of HPT-ed materials is their vanished work-hardening ability, which is in keeping with the assumption of elastic-perfect-plastic material in our simulation. In fact, the material undergoes deformation even at the center of the disk during HPT. The Berkovich hardness is 1.2 GPa at the niobium disk center [15], corresponding to about 400 MPa for yield stress Y from Figure 12. Considering that the strength of as received niobium is less than 200 MPa, the enhancement in strength and at the same time reduction in ductility at the disk center would be significant. Therefore, assuming elastic-perfect-plastic material property at the disk center will not lead to essential error. The HPT process does not change the elastic modulus remarkably. Thus, the value of Y/E is larger at the edge than at the center of the disk. Figure 13 presents the indentation geometry at the center (a) and edge (b) of the HPT processed niobium disk. It can be observed that significant pileup appears at the center (a2), while neither sink-in nor pileup is observed at the edge (b2). The values of hf/hmax are 0.93 and 0.81 for Figures 13(a) and 13(b), respectively. The corresponding values of Y/E are estimated to be 0.0057 and 0.019 from Figure 7. This experimental observation of pileup or sink-in is well predicted by Figure 9.Figure 13
Appearance of the indentation geometry by AFM at the center (a) and edge (b) of the HPTed niobium disk. It can be observed that, at the center of the disk, significant pileup effect is present (a2), while, at the edge, neither sink-in nor pileup effect is observed (b2) [15].
(a)
(b)


## 5. Summary and Conclusion

Three kinds of indentation, namely, conical, Berkovich, and Vickers, were comparatively studied by means of numerical simulations. The effect of pileup on the derived quantities was carefully examined. Dimensionless analysis was adopted to acquire the principal governing parameters of the process. The main findings are summarized in the following:(1)
The dimensionless functionsϕPL, ϕA, ϕhf, and ϕS, which represent loading, projected contact area, final indentation depth, and contact stiffness, respectively, are all functions of Y/E only for the elastic-perfect-plastic material considered in this work. However, these functions are not the same for different kinds of indenters.(2)
Whenhf/hmax>0.85, the effect of pileup becomes significant. Failure to account for the effect of pileup may lead to more than 50% error in calculating elastic modulus and hardness.(3)
The derived elastic modulus could be different if different indenters are adopted, even if the effect of pileup is considered. An approximate relationship ofE of the sample measured using different indenters can be given by EBerkovich=1.06EVickers=1.13EConical. This estimation leads to less than 5% error when hf/hmax<0.95 (corresponding to Y/E>0.002).(4)
The Tabor factor,C=H/Y, is not a constant but varies with hf/hmax. For soft material, that is, hf/hmax close to unity (or Y/E close to zero from Figure 7), C reaches 2.7, 3.3, and 3.6 for conical, Vickers, and Berkovich indentation, respectively.The above conclusions are acquired based on a group of specific materials whose strain hardening is ignored and the effect of Poisson’s ratio is not considered. Moreover, the time-dependent (or strain-rate dependent) characteristics of the indentation process for some materials are not reckoned in for present work. Therefore, application of these results should be conditional. The role of strain hardening and strain rate hardening might be significant in determining the properties of the testing materials and should be further investigated.



---

*Source: 528729-2015-09-06.xml*