Next Article in Journal
Liquid Copper and Iron Production from Chalcopyrite, in the Absence of Oxygen
Next Article in Special Issue
Control of Welding Residual Stress in Large Storage Tank by Finite Element Method
Previous Article in Journal
High-Temperature Tensile Properties of Hastelloy X Produced by Laser Powder Bed Fusion with Different Heat Treatments
Previous Article in Special Issue
Micromechanics-Based Low Cycle Fatigue Life Prediction Model of ECAPed Aluminum Alloy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Binary Newton Calculation Method of Residual Stress Based on the Indentation Energy Difference Theory

College of New Energy, China University of Petroleum (East China), Qingdao 266580, China
*
Author to whom correspondence should be addressed.
Metals 2022, 12(9), 1439; https://doi.org/10.3390/met12091439
Submission received: 4 August 2022 / Revised: 26 August 2022 / Accepted: 26 August 2022 / Published: 29 August 2022
(This article belongs to the Special Issue Residual Stress and Fatigue of Metals)

Abstract

:
Residual stress is a key parameter to evaluate the structural reliability of energy equipment. The indentation method has the characteristics of being nondestructive and easy to operate to calculate the residual stress of test materials, which has a broad application prospect in the field testing of energy equipment. However, because of the effect of preloading and data acquisition delay, the problem of indentation data fluctuation is prominent, and the indentation energy coefficient fitted by the traditional least square method is not consistent with the theoretical law, making it difficult to carry out stable calculations. In this paper, the Newton iteration formula of a binary nonlinear formula is derived based on the univariate Newton iteration formula, which is introduced into the data processing of residual stress, which increases the weight of the data in the stability stage and reduces the influence of the fluctuation data on the fitting results, so that the indentation energy coefficient is accurately calculated. Combined with the basic principle of indentation energy difference theory, the precise and efficient measurement of residual stress is realized.

1. Introduction

Residual stress [1,2] refers to the internal stress that remains in an object to maintain self-phase equilibrium after the disappearance of external force or the action of an inhomogeneous temperature field [3]. The existence of residual stress will not only reduce the structural strength but also cause stress corrosion cracking and fatigue damage, which destroy the integrity of the structure. Therefore, in order to ensure the safety and reliability of the structure, it is necessary to accurately test the magnitude and distribution of residual stress and formulate a pertinent control strategy. At present, the testing methods for residual stress mainly include the mechanical test method and the physical test method [4]. The mechanical testing methods mainly include the drilling method and slotting method, which cause local damage to materials during the measurement process and are difficult to use in service or waiting-for-service equipment; however, physical test methods, such as X-ray diffraction and ultrasound, fluctuate greatly when testing nonuniform tissue and strong stress gradient structures (such as welded structures), and their application in engineering sites is limited [1,5,6]. In contrast, the indentation method has the advantage of convenience without sampling and testing and has a good engineering application potential in the field test of service equipment.
In the early stages, it was found that the hardness value decreased when residual tensile stress existed; when residual compressive stress was present, the hardness value of the test increased. Therefore, by comparing the variation pattern of hardness values, the magnitude and direction of the residual stress can be estimated qualitatively. Subsequently, Tsui and Bolshakov [7] found that the residual stress did not change the true hardness of the material but affected the real indentation area. In 1998, Suresh et al. [8] put forward the indentation method for quantitatively testing residual stress based on the assumption of hardness invariance. However, the test results of this method depend on the indentation topography, and only the equivalent average residual stress can be measured. Aiming at this problem, Peng et al. [9] proposed a new method to calculate the biaxial residual stress by means of the method of indentation energy difference. In this method, the Knoop [10,11] indenter was indented in the direction of biaxial stress, the indentation energy obtained by the integration of the load–depth curve was used as the calculation parameter, and the biaxial stress component was obtained without observing the indentation topography. However, in actual measurement, because of the delay of the control system of the indentation equipment and the large fluctuation of the data points, the original data fitting effect was poor, the residual stress calculation result deviated greatly from the actual result, and sometimes cannot even be calculated.
In view of the above problems, based on the univariate Newton iteration method [12,13], the Newton iteration formula [14] of binary nonlinear formulas was deduced, which is introduced into the fitting process of the indentation data. By increasing the weight of the data in the stability stage, the influence of the fluctuation data on the fitting results is reduced, and the accurate and efficient solution of calculation of the biaxial residual stress based on energy difference theory is realized.

2. Residual Stress Calculation Model Design Based on the Newton Method

According to the theory of indentation energy difference, obtaining the biaxial residual stress of the material requires three indentation tests on the specimen: the first is the test in the stress-free state, and the other two are the orthogonal indentation of the Knoop indenter on the surface of the specimen. The resulting load–depth curve is shown in Figure 1. From an energy point of view, the indentation work and the stress work are all converted into strain energy during each indentation process [15,16,17]. Considering that the stress work cannot be directly solved, the residual stress can be correlated by calculating the indentation energy difference between the stressed and the stress-free states [18].
The calculation of the indentation energy difference at the same indentation depth requires the indentation energy of the three indentation processes, which can be directly obtained by integrating the load–depth curve, as shown in Formula (1) and Figure 2. Since the indentation curve is a quadratic function, the indentation energy can be obtained by substituting the slope of the curve instead of the integral according to the integral solution method of the quadratic function, in which the curve slope is called the indentation energy coefficient, as shown in Formula (2).
W F = 0 h m a x F d h = 0 h m a x C h 2 d h
W F = 1 3 C h m a x 3
It can be seen from Formula (2) that the accurate fitting of the load–depth data is the key to the solution of the indentation energy in the study of welding residual stress measured by the energy difference method [19,20]. In his study, Peng Wei found that due to the deviation of the shape of the equipment indenter and the delay of the control system, the existence of unstable region data caused the curve fitted by Formula (3) to have a large coincidence error with the origin in the initial position, while the fitting with Formula (4) could restrict the original data, and the slope obtained by Formula (3) was compared with that of Formula (3). The fitting curves of the stress-free and the unidirectional stress state do not conform to the theoretical model.
F = a h 2 + b h + c
F = C h + u 2
In order to fit the indentation data accurately, the author designed an indentation data fitting algorithm based on Newton method. The author first calculated the sum of variance S2 of the experimental data and the fitting formula
S 2 = i = 1 n F i C h i + u 2 2
When S2 is the smallest, the fitting effect is best at this time, i.e., the function must have a unique extremum and a minimum. The basic calculation method based on the extremum of the function should meet the requirements of
S 2 C = 2 C γ 4 , 0 + 8 C u γ 3 , 0 + 12 C u 2 γ 2 , 0 + 8 C u 3 γ 1 , 0 + 2 n C u 4 2 γ 2 , 1 4 u γ 1 , 1 2 u 2 γ 0 , 1 = 0 S 2 u = 4 C 2 γ 3 , 0 + 12 C 2 u γ 2 , 0 + 12 C 2 u 2 γ 1 , 0 + 4 n C 2 u 3 4 C γ 1 , 1 4 C u γ 1 , 0 = 0
γ j , k = i = 1 n u i j F i k ;   k = 0 ,   γ j , 0 = i = 1 n u i j ;   j = 0 ,   γ 0 , k = i = 1 n F i k ;   j   a n d   k   a r e   n o t   b o t h   0 .
Thus, the fit of Formula (3) is transformed into solving a binary nonlinear system of Formula (6).
Based on the Newton iterative formula of nonlinear formulas of binary functions, the author rewrote the abovementioned binary nonlinear formulas into the following functions:
A C , u = C r 4 , 0 + 4 C u r 3 , 0 + 6 C u 2 r 2 , 0 + 4 C u 3 r 1 , 0 + n C u 4 r 2 , 1 2 u r 1 , 1 u 2 r 0 , 1 B C , u = C 2 r 3 , 0 + 3 C 2 u r 2 , 0 + 3 C 2 u 2 r 1 , 0 + n C 2 u 3 C r 1 , 0 C u r 0 , 1
This function satisfies the basic requirements of Newton method for solving binary nonlinear formulas, so it can be solved by substituting Newton iterative formula.
C k + 1 = C k + A C k , u k B h C k , u k B C k , u k A u C k , u k B c C k , u k A u C k , u k A c C k , u k B u C k , u k u k + 1 = u k + B C k , u k A C C k , u k A C k , u k B C C k , u k B C C k , u k A u C k , u k A C C k , u k B u C k , u k
A c C , u , A u C , u , B c C , u , B u C , u is the first-order partial differential of the function A C , u and B C , u for c and u.
The indentation energy coefficient can be accurately fitted by using the binary Newton method, and the three-indentation data are correlated according to the calculated indentation energy to obtain the formula of the biaxial indentation energy difference shown in Figure 3 (when the indentation energy difference is equal to the same depth, the stress-free indentation work minus the stressed indentation work).
Δ W F 1 = W F 0 W F 1 Δ W F 2 = W F 0 W F 2
Peng found that the energy difference of the Knoop indenter had superposition and directional characteristics after a lot of experimental studies [9,21,22]. The superposition characteristic [23] points out that the indentation energy difference under the biaxial stress is equal to the sum of the two stress components of the biaxial stress as the uniaxial stress, respectively, and the difference in the relative direction is the same as the original biaxial stress, as shown in Formula (10). The directional characteristics show that the difference of the indentation energy produced when the long axis of the Knoop indenter is indented perpendicular to the stress direction is about φ times that in parallel, and φ is the directional coefficient of the Knoop indenter [24], as shown in Formula (11).
Δ W F = Δ W 11 + Δ W 12
Δ W 21 = φ Δ W 11
The combination of the above two formulas immediately converts the indentation energy difference under any biaxial stress to the sum of the difference of the two equivalent uniaxial stresses [25,26] perpendicular to the long axis direction of the indenter to obtain Formula (12)
Δ W 1 e q = φ φ 2 1 φ Δ W F 1 Δ W F 2 Δ W 2 e q = φ φ 2 1 φ Δ W F 2 Δ W F 1
Biaxial residual stress components can be calculated according to the relationship between the difference of the indentation energy and residual stress.
Δ W i e q = k V σ i + j = 1 3 A j σ i σ y j , i = 1 , 2
At this point, the Knoop indentation method is used to measure biaxial residual stress model. The binary Newton iteration formula commonly used in engineering calculation is used in the above model to process the measured data to obtain the energy coefficient, which avoids the error caused by the error of the specific point in the test data to the calculation results. Combined with the indentation energy difference theory, the problem that the calculation error caused by the traditional fitting method is large or even does not conform to the practical law in the actual calculation process is solved and provides theoretical support for the algorithm design of the residual stress calculation flow automation shown in Figure 4.

3. Derivation of the Newton Formula for Binary Equation

The Newton iteration method is one of the most powerful techniques to solve the problem of nonlinear equations. In this method, a series of approximate solutions are found to solve the nonlinear equation by selecting an initial estimator and approaching iteratively.
For the univariate equation: f(x) = 0, the Newton iteration formula is as follows
x n + 1 = x n f x n f x n
where n is the number of iterations, x n , x n + 1 is the approximate root of equation f x = 0 ,   and   f x n is the first-order differential of the function f x at x n .
After understanding the iterative solution process of univariate nonlinear equations, the Newton iteration method for binary nonlinear equation can be deduced.
If we have continuous partial derivatives of more than two orders in a certain neighborhood of function z = f x , y   a t   x 0 ,   y 0 , and x 0 + h , y 0 + k is at any point within this neighborhood, there is
f x 0 + h , y 0 + k f x 0 ,   y 0 +   h x f x , y | x = x 0 + k y f x , y | y = y 0
h = x x 0 ,   k = y y 0
If f x , y = 0 , the above equation can be converted to
f x k , y k + x x k f x x k , y k + y y k f y x k , y k = 0
Similarly, let z = g x , y be continuous in a certain neighborhood at point x 0 ,   y 0 , and there are continuous partial derivatives up to the 2nd order; then, x 0 + h , y 0 + k at any point within this neighborhood will also have
g x 0 + h , y 0 + k g x 0 ,   y 0 +   h x g x , y | x = x 0 + k y g x , y | y = y 0
h = x x 0 ,   k = y y 0
Then, equation g x , y = 0 can be reduced approximately to
g x 0 + h , y 0 + k g x 0 ,   y 0 +   h x g x , y | x = x 0 + k y g x , y | y = y 0
Therefore, the equation set can be obtained as follows:
f x k , y k + x x k f x x k , y k + y y k f y x k , y k = 0 g x k , y k + x x k g x x k , y k + y y k g y x k , y k = 0
For the above equation, when g x   x k , y k f y x k , y k f x x k , y k g y x k , y k 0 , it can be converted to
x = x k + f x k , y k g y x k , y k g x k , y k f y x k , y k g x x k , y k f y x k , y k f x x k , y k g y x k , y k y = y k + g x k , y k f x x k , y k f x k , y k g x x k , y k g x x k , y k f y x k , y k f x x k , y k g y x k , y k
The Newton iteration method for nonlinear equations of binary functions can be used to obtain approximate values by the iteration method when k is taken as a non-zero natural number x k , y k . δ are the error controls required for data processing; when x k + 1 , y k + 1 δ δ > 0 , the root of the original equation is x k , y k . Then, the Newton iteration formula of the nonlinear equation of the binary function is as follows:
x k + 1 = x k + f x k , y k g y x k , y k g x k , y k f y x k , y k g x x k , y k f y x k , y k f x x k , y k g y x k , y k y k + 1 = y k + g x k , y k f x x k , y k f x k , y k g x x k , y k g x x k , y k f y x k , y k f x x k , y k g y x k , y k

4. Results and discussion

4.1. Analysis of Newton Fitting Effect of Indentation Data

In order to verify the accuracy of the calculation results, the indentation data were obtained by multiple tests using the indentation equipment in the adjacent area of Q345 weldless sample steel plate. After the unloading stage was removed by data processing, the indentation data were fitted according to the above algorithm and compared with the fitting results of the traditional one-dimensional binary equation. The fitting results are shown in Table 1. The indentation energy coefficient results of the traditional fitting method were 66,400 ± 9000 and the Newton fitting method were 75,200 ± 3000. This indicates that the repeatability of the energy coefficient C obtained by the Newton method was very high, while the numerical error obtained by the traditional one-dimensional binary equation fitting method was large. In contrast, the residual stress data fitting method designed by Newton method was accurate, fast, and stable, which can meet the requirements of data processing.
The standard deviation analysis of the indentation curve obtained by the binary Newton fitting method and the raw data is shown in Figure 5. The indentation data can be divided into a stable fitting area and an unstable fitting area according to the maximum error of data fitting in the middle stage of the acquisition process. Most of the indentation data can meet the requirements of residual stress measured by indentation method, so it is in the stable fitting area. Since the preloading process and the delay of the sensor during the measurement with the indentation instrument will cause errors in the data at the initial and end of loading phases of the indentation, these values must be removed during the calculation of the end of measurement as “noise”, as they are in the unstable fitting region because they do not conform to the measurement law.
Because of the large amount of data and the difficulty of analyzing the error data, it is almost impossible to remove the “noise” manually. While this issue was not considered in the traditional fitting process, as all data were fitted with the same weight, and with this method, although the degree of fit was high, the results deviated to a large extent, which was an overfitting. It can be seen from the graph that the Newton fitting algorithm had a good fitting effect in the middle stage of the curve, indicating that the middle stage data (stable fitting area) occupied a significant weight in the fitting process, while the standard deviation of the initial stage and the end stage (unstable fitting area) with more error data was large, and the fitting weight was small. Therefore, the binary Newton fitting method increased the weight of the steady-phase data in the fitting process and reduced the influence of fluctuating data on the fitting results.

4.2. Measurement and Analysis of Biaxial Residual Stress

Because X-ray diffraction (XRD) has the advantage of high measurement accuracy, the accuracy of the Newton method in calculating residual stress was verified by comparing X-ray diffraction and the indentation method.
The testing material was an Inconel718 alloy welded plate. The width of the fusion zone was 20 mm, and the residual stress was measured by the X-ray diffraction method and the indentation method in the weld metal, as shown in Figure 6. The nonuniform temperature field caused by the welding process was the essential reason for the residual stress. The higher the temperature, the greater the thermal expansion. However, the thermal expansion will be constrained by the surrounding material subjected to relatively low temperature and vice versa. The uneliminated stress is called the welding residual stress. The measurement results are shown in Figure 7, where σ1 is the residual stress parallel to the welding direction, and σ2 is the residual stress perpendicular to the welding direction.
It can be seen from the test results that the residual stress measured by the indentation method was more stable than that by X-ray diffraction. The average error between the indentation method and the X-ray diffraction method was within 10%, and the fluctuation of the calculation result was small, which proves that the test result was stable and reliable and has good engineering applicability.

5. Conclusions

Indentation equipment is poorly fitted to the original data due to the delay of its control system and the large fluctuation of data points. The traditional fitting method does not conform to the theoretical law due to its overfitting of the loading slope (also called the indentation energy coefficient), which leads to large deviation of the measurement results from the actual results, and subsequent calculations cannot be performed. Therefore, in this paper, starting from the Newtonian solution of the nonlinear binary system of equation and combining it with the fitting of indentation data, the solution formula of the residual stress indentation energy coefficient was derived, and the interference of the non-stable zone data was excluded by the binary Newton iteration. The main work was as follows.
(1) The Newton iterative formula for the root of the binary nonlinear equation was deduced through the relationship between the univariate Taylor formula and the Newton iteration method for solving the univariate nonlinear equation. The indentation energy coefficient problem was transformed into the root problem of the binary nonlinear equation. The Newton iterative formula was introduced, and the fast calculation formula of the indentation energy coefficient was established.
(2) The fitting effects of the traditional least squares method and the binary Newton method were compared and analyzed. It was found that the traditional least squares method fitted all data with the same weight, the influence of the data fluctuation could not be ruled out, and the resulting deviation was too large due to overfitting. The binary Newton method improved the fitting weight of the indentation curve in the stable fitting region, and the fitting result was more stable and accurate.
(3) The stability of the binary Newton method for calculating the energy coefficient was verified by using the Q345 weldless sample steel plate test, and the accuracy of the binary Newton calculation method of residual stress based on the theory of indentation energy difference was verified by comparison with the X-ray method.

Author Contributions

Conceptualization, D.K. and B.Y.; methodology D.K.; software, D.K.; validation, D.K., B.Y., and P.H.; formal analysis, D.K. and B.Y.; investigation, B.Y.; resources, B.Y.; data curation, D.K.; writing—original draft preparation, D.K. and B.Y.; writing—review and editing, D.K., B.Y. and P.H.; visualization, D.K.; supervision, B.Y.; project administration, B.Y.; funding acquisition, B.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grants No.: 51805546), the Nature Science Foundation of Shandong Province (Grants No.: ZR2019BEE050).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in this article.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

C Indentation energy coefficient k Elasticity coefficient
u Fit Deviation V Volume of Knoop indenter intruded into the specimen
F LoadφKnoop indenter direction coefficient
h Depth Δ W i e q Indentation energy difference under uniaxial stress perpendicular to the long axis of the indenter
W F Indentation work A j Material-related plasticity coefficient
S Sum of all data fit standard deviations Δ W i j Indentation energy difference. i, j are the indenter direction and stress direction, respectively
Δ W F i Indentation energy difference in one direction W F i Indentation work in a certain direction
a , b , c Three coefficients of a quadratic function of one variable γ i , j Parameters associated with parameters u and F. i and j are the powers of u and F, respectively
A C , u The nonlinear function constructed in the process of solving C B C , u The nonlinear function constructed in the process of solving C

References

  1. Huchings, M.T.; Withers, P.J.; Holden, T.M.; Lorentzen, T. Introduction to the characterization of residual stress by neutron diffraction. Exp. Dermatol. 2005, 13, 505–511. [Google Scholar] [CrossRef]
  2. Peng, W.; Jiang, W.; Jin, Q.; Wan, Y.; Luo, Y.; Ren, L.; Zhang, K.; Tu, S.-T. Reduction of welding residual stress in the head-cylinder joint of a large rectifying tower by finite element method and experimental study. Int. J. Press. Vessel. Pip. 2021, 191, 104311. [Google Scholar] [CrossRef]
  3. Yan, S.; Meng, Z.; Chen, B.; Tan, C.; Song, X.; Wang, G. Prediction of temperature field and residual stress of oscillation laser welding of 316LN stainless steel. Opt. Laser Technol. 2022, 145, 107493. [Google Scholar] [CrossRef]
  4. Huang, X.-D.; Zhang, X.-M.; Ding, H. A novel relaxation-free analytical method for prediction of residual stress induced by mechanical load during orthogonal machining. Int. J. Mech. Sci. 2016, 115–116, 299–309. [Google Scholar] [CrossRef]
  5. Liu, J.; Yuan, H. Prediction of residual stress relaxations in shot-peened specimens and its application for the rotor disc assessment. Mater. Sci. Eng. A 2010, 527, 6690–6698. [Google Scholar] [CrossRef]
  6. Marshall, D.B.; Lawn, B.R. Residual stress effects in sharp contact cracking. J. Mater. Sci. 1979, 14, 2001–2012. [Google Scholar] [CrossRef]
  7. Bolshakov, A.; Oliver, W.C.; Pharr, G.M. Influences of stress on the measurement of mechanical properties using nanoindentation: Part II. Finite element simulations. J. Mater. Res. 1996, 11, 760–768. [Google Scholar] [CrossRef]
  8. Suresh, S.; Giannakopoulos, A.E. A new method for estimating residual stresses by instrumented sharp indentation. Acta Mater. 1998, 46, 5755–5767. [Google Scholar] [CrossRef]
  9. Peng, W.; Jiang, W.; Sun, G.; Yang, B.; Shao, X.; Tu, S.-T. Biaxial residual stress measurement by indentation energy difference method: Theoretical and experimental study. Int. J. Press. Vessel. Pip. 2022, 195, 104573. [Google Scholar] [CrossRef]
  10. Phani, P.S.; Oliver, W.C.; Pharr, G.M. Influences of elasticity on the measurement of power law creep parameters by nanoindentation. J. Mech. Phys. Solids 2021, 154, 104527. [Google Scholar] [CrossRef]
  11. Li, W.; Huang, C.; Yu, M.; Liao, H. Investigation on mechanical property of annealed copper particles and cold sprayed copper coating by a micro-indentation testing. Mater. Des. 2013, 46, 219–226. [Google Scholar] [CrossRef]
  12. Crisfield, M.A. A faster modified newton-raphson iteration. Comput. Methods Appl. Mech. Eng. 1979, 20, 267–278. [Google Scholar] [CrossRef]
  13. Marshall, D.B.; Noma, T.; Evans, A.G. A simple method for determining elastic-modulus–to-hardness ratios using Knoop indentation measurements. J. Am. Ceram. Soc. 1982, 65, c175–c176. [Google Scholar] [CrossRef]
  14. Momani, S.; Abuasad, S. Application of He’s variational iteration method to Helmholtz equation. Chaos Solitons Fractals 2006, 27, 1119–1123. [Google Scholar] [CrossRef]
  15. Kang, S.K.; Kim, J.Y.; Park, C.P.; Kim, H.U.; Kwon, D. Conventional Vickers and true instrumented indentation hardness determined by instrumented indentation tests. J. Mater. Res. 2010, 25, 337–343. [Google Scholar] [CrossRef]
  16. He, A.; Liang, Y.; Zhao, O. Flexural buckling behaviour and resistances of circular high strength concrete-filled stainless steel tube columns. Eng. Struct. 2020, 219, 110893. [Google Scholar] [CrossRef]
  17. Saadatfard, H.; Niknejad, A.; Liaghat, G.; Hatami, S. A novel general theory for bending and plastic hinge line phenomena in indentation and flattening processes. Thin-Walled Struct. 2019, 136, 150–161. [Google Scholar] [CrossRef]
  18. Moharrami, R.; Sanayei, M. Improvement of indentation technique for measuring general biaxial residual stresses in austenitic steels. Precis. Eng. 2020, 64, 220–227. [Google Scholar] [CrossRef]
  19. Deng, D. FEM prediction of welding residual stress and distortion in carbon steel considering phase transformation effects. Mater. Des. 2009, 30, 359–366. [Google Scholar] [CrossRef]
  20. Dean Deng, H.M. Prediction of welding residual stress in multi-pass butt-welded modified 9Cr–1Mo steel pipe considering phase transformation effects. Comput. Mater. Sci. 2006, 37, 209–219. [Google Scholar] [CrossRef]
  21. Lee, J.H.; Lee, H.; Hyun, H.C.; Kim, M. Numerical approaches and experimental verification of the conical indentation techniques for residual stress evaluation. J. Mater. Res. 2010, 25, 2212–2223. [Google Scholar] [CrossRef]
  22. Larsson, P.-L. On the determination of biaxial residual stress fields from global indentation quantities. Tribol. Lett. 2014, 54, 89–97. [Google Scholar] [CrossRef]
  23. Lee, Y.-H.; Kwon, D. Estimation of biaxial surface stress by instrumented indentation with sharp indenters. Acta Mater. 2004, 52, 1555–1563. [Google Scholar] [CrossRef]
  24. Kim, Y.-C.; Choi, M.-J.; Kwon, D.; Kim, J.-Y. Estimation of principal directions of Bi-axial residual stress using instrumented Knoop indentation testing. Met. Mater. Int. 2015, 21, 850–856. [Google Scholar] [CrossRef]
  25. Peng, G.; Xu, F.; Chen, J.; Wang, H.; Hu, J.; Zhang, T. Evaluation of non-equibiaxial residual stresses in metallic materials via instrumented spherical indentation. Metals 2020, 10, 440. [Google Scholar] [CrossRef]
  26. Zhang, T.; Cheng, W.; Peng, G.; Ma, Y.; Jiang, W.; Hu, J.; Chen, H. Numerical investigation of spherical indentation on elastic-power-law strain-hardening solids with non-equibiaxial residual stresses. MRS Commun. 2019, 9, 360–369. [Google Scholar] [CrossRef]
Figure 1. Load-Displacement Curve Obtained by Indentation equipment.
Figure 1. Load-Displacement Curve Obtained by Indentation equipment.
Metals 12 01439 g001
Figure 2. Graphical Representation of Indentation energy.
Figure 2. Graphical Representation of Indentation energy.
Metals 12 01439 g002
Figure 3. Graphical representation of biaxial indentation energy difference: (a) the indentation energy difference when the indenter is pressed crosswise; (b) the indentation energy difference when the indenter is pressed longitudinally.
Figure 3. Graphical representation of biaxial indentation energy difference: (a) the indentation energy difference when the indenter is pressed crosswise; (b) the indentation energy difference when the indenter is pressed longitudinally.
Metals 12 01439 g003
Figure 4. Flow chart of biaxial residual stress algorithm.
Figure 4. Flow chart of biaxial residual stress algorithm.
Metals 12 01439 g004
Figure 5. Analysis of standard deviation between fit curve and raw data: (a) fitting standard deviation of indentation measurement points; (b) indentation curve fitting stable zone division.
Figure 5. Analysis of standard deviation between fit curve and raw data: (a) fitting standard deviation of indentation measurement points; (b) indentation curve fitting stable zone division.
Metals 12 01439 g005
Figure 6. Inconel 718 weld zone residual stress measurement diagram.
Figure 6. Inconel 718 weld zone residual stress measurement diagram.
Metals 12 01439 g006
Figure 7. Measurement of Residual Stress in Inconel718 Weld Zone by X-ray Diffraction and the Indentation Method.
Figure 7. Measurement of Residual Stress in Inconel718 Weld Zone by X-ray Diffraction and the Indentation Method.
Metals 12 01439 g007
Table 1. Comparison of The Compression Energy Coefficient C between the Newton Method and the Conventional Method.
Table 1. Comparison of The Compression Energy Coefficient C between the Newton Method and the Conventional Method.
Standard Value12345678910
Newton method75,000 ± 500073,76573,03778,77172,17578,04977,50779,44675,19873,72772,405
Conventional method68,19565,26870,13965,39863,50272,91875,66865,64062,37962,375
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kong, D.; Yang, B.; Hu, P. Binary Newton Calculation Method of Residual Stress Based on the Indentation Energy Difference Theory. Metals 2022, 12, 1439. https://doi.org/10.3390/met12091439

AMA Style

Kong D, Yang B, Hu P. Binary Newton Calculation Method of Residual Stress Based on the Indentation Energy Difference Theory. Metals. 2022; 12(9):1439. https://doi.org/10.3390/met12091439

Chicago/Turabian Style

Kong, Desheng, Bin Yang, and Peiyu Hu. 2022. "Binary Newton Calculation Method of Residual Stress Based on the Indentation Energy Difference Theory" Metals 12, no. 9: 1439. https://doi.org/10.3390/met12091439

APA Style

Kong, D., Yang, B., & Hu, P. (2022). Binary Newton Calculation Method of Residual Stress Based on the Indentation Energy Difference Theory. Metals, 12(9), 1439. https://doi.org/10.3390/met12091439

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop