A Further Study on Knoop Indentation Plastic Deformation for Evaluating Residual Stress

Article information

Korean J. Met. Mater.. 2020;58(8):515-521
Publication date (electronic) : 2020 July 27
doi : https://doi.org/10.3365/KJMM.2020.58.8.515
1Department of Materials Science and Engineering, Seoul National University, Seoul 08826, Republic of Korea
2Research Center for Energy and Clean Technology, School of Materials Science and Engineering, Andong National University, Andong 36729, Republic of Korea
*Corresponding Author: Young-Cheon Kim Tel: +82-52-217-2344, E-mail: kimyc@anu.ac.kr

Equal contribution

- 김우주·이경열·김종형: 박사과정, 김영천·권동일: 교수

Received 2020 May 20; Accepted 2020 June 15.

Abstract

A method for evaluating residual stress using an instrumented indentation test was developed some decades ago. More recently, another method was developed, using a Knoop indenter. The conversion factor ratio, which is one of the key factors in the evaluation algorithm, has been taken to be 0.34, although this value comes from an experimental result and its physical meaning has not been examined. Here we examine the physical meaning of this conversion factor from the previous residual stress model, and calculate its ratio using analytical model of the stress field beneath the indenter. In this process, we assumed that the conversion factor ratio was the ratio of the projected area of the plastic zone generated during the Knoop indentation test. An analysis of the stress field beneath the indenter was performed by FE simulation. Actual nanoindentation was conducted after Knoop indentation testing, using the interface-bonding technique, to identify the plastic zone. In addition, the conversion factor ratio was also calculated for the case where residual stress was present, and the geometric ratio of the Knoop indenter was different. A comparison of our results with those from previous studies showed that the conversion factor ratio obtained using our assumption was in good agreement with previous studies.

1. Introduction

Residual stress is generated in materials during heat treatment machining, and this stress has emerged as an important factor in assessing a product. Numerous studies have investigated ways of measuring residual stress, including with XRD [1], hole-drilling [2], and other methods. Instrumented indentation testing (IIT) is an increasingly popular method because it is non-destructive and the test procedure is simple. In addition, it can be used in the field and can be used to assess various other mechanical properties as well [3-5,25-26].

Methods of evaluating residual stress using IIT have expanded over the past decade [6-9]. Lee and Kwon [6,7] recently suggested a technique for evaluating non-equibiaxial stress states using a Vickers indenter. After that, Choi et al. [10] proposed a way to evaluate residual stress using a Knoop indenter that permits evaluation even when the ratio of principal stress directions of the residual stress is unknown. In that study, each load difference was expressed by introducing a conversion factor (α) in two directions (the horizontal and vertical directions) as shown in Fig 1. They expressed the relationship between the principal stresses and the load difference as

Fig. 1.

Schematic indentation loading curves for Knoop indentation

(1) Lx = α//σresx + ασresy
(2) Ly = ασresx + α//σresy
(3) p = σresyσresx = LxLy - α//α1 - α//αLxLy

where σresx and σresy are the residual stresses in the x- and y-directions (assuming the x- and y-directions are the principal directions) respectively and ∆Lx and ∆Ly are the load differences obtained when the long Knoop axis is in the x- and y-direction. α// and α are conversion factors. Their ratio, α to α//, is called the conversion factor ratio and is essential to determining the stress ratio p. As the conversion factor ratio in their Knoop indentation studies, Han et al. [9] and Choi et al. [10] used 0.34, obtained from experimental results. In addition, they expressed the change in conversion factor ratio according to the change in the geometry of the Knoop indenter as a fitting equation using FE simulation. This research was extended by Kim et al. [11,12], but the physical meaning of this conversion factor ratio is not yet clear.

The present study investigates the physical meaning of the conversion factor ratio, α to α//. To identify it, we studied it theoretically, by analyzing the stress field beneath the indenter and also performed indentation tests. Then, from FE analyses of various situations where residual stress was present, we defined the physical meaning of the conversion factor ratio by comparison with previous experimental results.

2. Materials and Methods

2.1 Finite element analysis

Finite element analysis was used to study the size of the plastic zone beneath the indenter. Using Abaqus 6.14, a Knoop indenter was designed with a 7.11:1 ratio of major to minor axis and ∠172· and ∠130· respectively. The indenter body was a 3D discrete rigid type and consisted of 2114, 118 elements of the R3D4 type. The specimen had a cuboid shape, a 3D deformable type, and a C3D8 8-node brick element type of 67,156 elements. For more precise results, the mesh was set to fine at the indent center and was increased moving away from the center. The specimen was fixed for the simulation by applying boundary conditions at the bottom. The mechanical properties of the specimen were based on STS304L material, with an elastic modulus of 200 GPa, Poissons’s ratio 0.3, yield strength 360MPa and tensile strength 750MPa. The Knoop indentation simulation allowed the indenter body to penetrate 100 μm into the specimen surface and return to the original position considering nonlinear geometry.

2.2 Indentation testing

Knoop indentation and nanoindentation were performed to compare the results of the FE analyses with the plastic zone generated during indentation. The materials were 3 mm × 3 mm × 1 mm S20C and STS304L cuboids. To observe the plastic zone beneath the Knoop indenter, an interface-bonding technique [13,14] was used. This technique is well known as a way to observe the deformation beneath the indent in the indentation testing of metallic materials. Even though the plastic behavior during indentation in an interfaced specimen does not completely match that in bulk indentation, several studies [13-20] have reported that similar deformation mechanics can be observed. Specimen preparation is as follows. First, the specimen is cut exactly in half. Next, the cut surfaces are polished precisely and bonded with a high-strength adhesive. Then, Knoop indentation is performed along the long and short axes of the cut line. After that, the bonded specimen is immersed in acetone to dissolve the bond. Finally, nanoindentation is carried out on the crosssection of the de-bonded surface.

Knoop indentation was carried out on the top of the bonded specimen. AIS3000 (Frontics. Inc., Seoul, Korea) was used for indentation testing with 588 N (60 kgf) loading and unloading. For the nanoindentation, the cross-section of the specimen was polished once again with 0.3 μm alumina particles. Nanoindentations were made in the area beneath the Knoop indent using an Ultra Nanoindentation system (Anton Paar, Switzerland). In the experiments, a Berkovich indenter was used and the maximum loading was 50 mN.

3. Results and Discussion

In previous studies [9,10], the conversion factor (α) in the Knoop indentation test was defined as:

(4) σres α = ΔL

Since the conversion factor is a relationship between stress and load, it can be inferred that it has units of area. Since the Knoop indenter has two-fold symmetry and axes of different lengths (unlike the usual Vickers indenter), it can respond more sensitively to the residual stress direction. This implies that the conversion factor is related to the projected area where the residual stress arises.

For simplicity, imagine that the projected areas corresponding to the short-axis direction and long-axis direction of the Knoop indenter are taken as the conversion factors. Since the areas in the short-axis direction and the long-axis direction have the same height, their ratio is the length of the short-axis and the long-axis and the conversion factor ratio can be expressed as the ratio of these lengths. The difference between the value of this ratio (about 0.14 [1:7.11] from the indenter geometry) and that found in previous research [9,10] indicates that the projected area of the Knoop indenter geometry is not directly related to the conversion factor. However, considering the previous residual stress model [6,7], the indentation load difference (∆L) is related to the deviatoric stress and the deviatoric stress affects plastic deformation. From this, it can be inferred that the conversion factor, which introduces a relationship between the load difference (∆L) and the residual stress is related to the concept of plastic deformation rather than the Knoop indenter itself. Therefore, it is reasonable to look for the meaning of the conversion factor in plastic deformation during indentation.

The stress field beneath the indenter has been of interest to many researchers. The plastic zone can be analyzed simply using the hemisphere-shaped plastic zone, and the cavity model [21]. In this way, we can approach the plastic zone mathematically. In the expanded cavity model of Gao et al. [22], an equation for the plastic zone radius beneath the indenter in conical indentation is proposed:

(5) (ca)3 = 13·Eσy·cotα

where a is the contact radius, c is the radius of the plastic zone, E is the elastic modulus, σythe yield strength and α the semi-included angle of the conical indenter. From equation (5), the projected area R of each plastic zone arising in a conical indentation having an angle corresponding to the long and short axes of the Knoop indenter can be calculated, and yields a ratio value of 0.27 (Areashortplastic zone/Arealongplastic zone) as the conversion factor ratio. This result is closer to 0.34, the experimental value used in previous studies, than the 0.14 obtained from the Knoop indenter geometry. However, it still differs from the previous result because it uses a conical indenter, which differs geometrically from the Knoop indenter. In particular, a pile-up primarily occurs in the short-axis direction, and less in the long-axis direction during the Knoop indentation [12,23]. As a result, the plastic behaviors in each direction from a conical indenter are quite different from each other. In addition, the expanded cavity model basically assumes that the plastic zone beneath the indenter is hemispherical, so it cannot be directly applied to the Knoop indentation. Thus, the ratio must be examined more precisely. To do this we used the FE simulation to study the plastic zone beneath the indenter.

Figure 2 shows the plastic zone beneath the indenter as observed in the Knoop indentation simulation. The plastic zone is identified in the FE analysis as being in the PEEQ 0.1 ~ 0.5 range; that zone is marked in Fig 2 and the other areas are marked in blue. The ratio of the projected areas of the plastic zone in the long-axis and short-axis directions in the simulation is 0.321. This result differs by about 5% from the ratio used in previous studies [9,10], which is in reasonably good agreement since the previous value was obtained experimentally. However, the experimental value of Han [9] was obtained by inverting the results in the presence of residual stresses, and therefore we need further simulation to determine the residual stressed state. Before the Knoop indentation simulation, 100 and 200 MPa of residual stress was applied in the x and y directions, respectively, and the ratio was calculated from the projected area of the plastic zone. The ratios of the projected plastic-zone area in each residual stress state, as shown in Table 1, have an error range of about 10% from the 0.34 obtained by Han [9]. This suggests that the real experiment and the simulation are not completely identical. In particular, Choi [10] and Kim [24] stressed the manufacturing limitations of indenters used for measuring residual stress. It is difficult to manufacture indenters with ideal ratios, unlike the indenters designed in simulations. In addition, the mechanical properties used in our FE simulation may differ from those from which the experimental relations in previous studies were derived. Bearing these factors in mind, we believe that the results obtained here show reasonable agreement with those of previous studies, and we believe that our result should be considered in agreement with Choi’s study: that the conversion factor ratio is independent of the residual stress [10].

Fig. 2.

Simulation of the plastic zone beneath a Knoop indentation

Ratio of the projected area of the plastic zone in FE simulation results and experimental result

Figure 3b)-c) shows the results of nanoindentation on the interface-bonded surface used to verify the plastic zone beneath the Knoop indent. As shown in Fig 3, nanoindentation was conducted in a region where the plastic zone was expected to occur and the hardness was measured on each cross-sectional region of the long and short axes. The average hardnesses of S20C and STS304L were 2.95 GPa (±0.34 GPa) and 5.95 GPa (±0.30 GPa), respectively. These values were obtained by performing nanoindentation on a region unaffected by Knoop indentation. The load-displacement curves of both the plastic zone and the unaffected region are shown as Fig 4a)-b). The hardness values from nanoindentation fluctuated slightly. The point with a value that was greater than that obtained from the side part is indicated in red in Fig 3b)-c) to distinguish it from points of equal or lower values than the average. Based on this, a virtual plastic semi-ellipsoid region (referring to the simulation results) was drawn. The depth of the ellipsoid plastic zone was the same in the cross-sectional area of the plastic zone, since the indentation depth was the same on both the long-axis and short-axis. In addition, since the area equation of the semi-ellipsoid is expressed as a function of the length of the indentation direction and the surface direction, if the depth is the same, then the ratio of the crosssectional areas is determined by the ratio of the length of the surface direction. The projected area ratios of the plastic zone in S20C and STS304L obtained from the above process were 0.32 and 0.35, respectively. These values are relatively consistent with the conversion factor ratio, 0.34, from, previous studies [9,10], and it is thus reasonable to assume that the conversion factor ratio is the ratio of the projected area of the plastic zone beneath the indenter.

Fig. 3.

a) Nanoindentation results on interface-bonded sample (b) S20C, c) STS304L)

Fig. 4.

Load-depth (displacement) curves of a) S20C and b) STS304L

The Knoop indenter used in this study had a ratio of long-axis to minor axis of 7.11:1. However, depending on difficulties in the experimental environment, such as a very narrow area or an elbow of a pipeline, the length in the x- and y-directions of the Knoop indenter may change, and that means the conversion factor ratio must also be prepared for each change. Choi et al. [10], in reporting this issue, expressed a fitting equation by observing the change of the conversion factor ratio according to the indenter geometry diagonal ratio through FE analysis. Kim et al. [24], also identified this change in their modified Berkovich work. Therefore, for our assumption to be valid, it is necessary to find the ratio of the plastic zone occurring during indentation relative to the diagonal change in the Knoop indenter.

Modified Knoop indenters with long to short-axis ratios of 3.5:1 and 14:1 were designed using FE simulation. Table 1 shows the ratio of the projected area of the plastic zone obtained using the analysis procedure described above. For diagonal ratios of 3.5:1 and 14:1, the ratios were 0.189 and 0.518, respectively; the conversion factor ratio increased when the long-axis was shorter compared to the long-axis of 7.11:1, and vice versa. These results are very similar to those in Choi and Kim [10,24] as shown in Fig 5, and support the conclusion that the conversion factor ratio depends on the indenter geometric shape and not on the magnitude of residual stress, as mentioned in their studies.

Fig. 5.

Relation between conversion factor ratio and diagonal ratio of Knoop indenter

4. Conclusion

We discussed the physical meaning of the conversion factor ratio used as an experimental result when estimating residual stress using indentation. Based on the estimation algorithm, we expected the conversion factor ratio to be related to plastic deformation. FEA studies and hardness mapping were conducted for verification, assuming that the conversion factor ratio was the ratio of the projected area of the plastic zone beneath the indenter. Our results were compared with experimental results from previous studies, and were confirmed to agree well, even when residual stress was present and the indenter diagonal was changed. Finally, based on the results of our analyses, we defined the conversion factor ratio to be the ratio of the projected area of the plastic zone under the indenter.

Acknowledgements

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. 2020R1A5A6017701), and the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (No. NRF-2018R1D1A3B07048712).

References

1. Withers P. J., Schajer G. S.. Practical Residual Stress Measurement Methods p. 163–194. Wiley. Chichester, UK: 2013.
2. Schajer G. S.. Exp. Mech 50:159. 2010;
3. Ahn J., Kwon D.. J. Mater. Res 16:3170. 2001;
4. Kim Y.-C., Kim K.-H., Kim J.-Y., Kwon D.. Int. J. Plast 49:1. 2013;
5. Jeon S.-W., Lee K.-W, Kim J. Y., Kim W. J., Park C. P., Kwon D.. Exp. Mech 57:1013. 2017;
6. Lee Y.-H., Kwon D.. Scr. Mater 49:459. 2003;
7. Lee Y.-H., Kwon D.. Acta. Mater 52:155. 2004;
8. Zeng K., Giannakopoulos A. E., Rowcliffe D., Meier P.. J. Am. Cera. Soc 81:689. 1998;
9. Han J. H., Lee J. S., Lee Y. H., Choi M. J., Lee G. J., Kim K. H., Kwon D. I.. Key Eng. Mater 345-346:1125. 2007;
10. Choi M.-J., Kang S.-K., Kang I., Kwon D.. J. Mater. Res 27:121. 2012;
11. Kim Y.-C., Choi M. J., Kwon D., Kim J. Y.. Met. Mater. Int 21:850. 2015;
12. Kim Y.-C., Ahn H., Kwon D., Kim J. Y.. Met. Mater. Int 22:12. 2016;
13. Yoo B.-G., Jang J.-I.. J. Phys. D: Appl. Phys 41:074017. 2008;
14. Yoo B.-G., Oh J.-H., Kim Y.-J., Jang J.-I.. Intermetallics 18:1872. 2010;
15. Mulhearn T. O.. J. Mech. Phys. Solids 7:85. 1959;
16. Jana S., Ramaurty U., Chattopadhyay K., Kawamura Y.. Mater. Sci. Eng. A. 375-377:1191. 2004;
17. Ramaurty U., Jana S., Bhowmick R., Kawamura Y., Chattopadhyay K.. Acta Mater 53:705. 2004;
18. Bhowmick R., Raghavan R., Chattopadhyay K., Ramaurty U.. Acta Mater 54:4221. 2006;
19. Zhang H., Zing X., Subhash G., Kecskes L. J., Dwoding R. J.. Acta Mater 53:3849. 2005;
20. Li W. H., Zhang T. H., Xing D. M., Wei B. C., Wang Y. R., Dong Y. D.. J. Mater. Res 21:75. 2006;
21. Hill R.. The mathematical theory of plasticity p. 97–103. Claredon press. Oxford, UK: 1950.
22. Gao X. L., Jing X. N., Subhash G.. Int. J. Solids Struct 43:2193. 2006;
23. Giannakopoulos A. E., Zisis Th.. Int. J. Solids struct 48:175. 2011;
24. Kim J.-H., Xu H., Choi M. J., Heo E., Kim Y.-C., Kwon D.. J. Mater. Res 33:3849. 2018;
25. Kim J.-h., Lee J., Kim W., Kim J., Kang S.-K., Kwon D.. Korean J. Met. Mater 57:289. 2019;
26. Kim W., Choi S., Kim J., Choi S., Choi M., Kwon D.. Met. Mater. Int 2020;https://doi.org/10.1007/s12540-020-00753-2.

Article information Continued

Fig. 1.

Schematic indentation loading curves for Knoop indentation

Fig. 2.

Simulation of the plastic zone beneath a Knoop indentation

Fig. 3.

a) Nanoindentation results on interface-bonded sample (b) S20C, c) STS304L)

Fig. 4.

Load-depth (displacement) curves of a) S20C and b) STS304L

Fig. 5.

Relation between conversion factor ratio and diagonal ratio of Knoop indenter

Table 1.

Ratio of the projected area of the plastic zone in FE simulation results and experimental result

Diagonal ratio of indenter FE simulation & experiment Residual stress state
Ave.
Stress-free X-direction 100 MPa X-direction 200 MPa Y-direction 100 MPa Y-direction 200 MPa
7.11:1 FEA 0.321 0.310 0.309 0.280 0.296 0.303
7.11:1 Experiment S20C 0.319 - - - - -
STS304L 0.349
14:1 FEA 0.199 0.189 0.191 0.184 0.185 0.189
3.5:1 FEA 0.504 0.527 0.521 0.534 0.507 0.518