In the digital image correlation research of fatigue crack growth rate, the accuracy of the crack tip position determines the accuracy of the calculation of the stress intensity factor, thereby affecting the life prediction. This paper proposes a Gauss-Newton iteration method for solving the crack tip position. The conventional linear fitting method provides an iterative initial solution for this method, and the preconditioned conjugate gradient method is used to solve the ill-conditioned matrix. A noise-added artificial displacement field is used to verify the feasibility of the method, which shows that all parameters can be solved with satisfactory results. The actual stress intensity factor solution case shows that the stress intensity factor value obtained by the method in this paper is very close to the finite element result, and the relative error between the two is only ??0.621%; The Williams coefficient obtained by this method can also better define the contour of the plastic zone at the crack tip, and the maximum relative error with the test plastic zone area is ??11.29%. The relative error between the contour of the plastic zone defined by the conventional method and the area of the experimental plastic zone reached a maximum of 26.05%. The crack tip coordinates, stress intensity factors, and plastic zone contour changes in the loading and unloading phases are explored. The results show that the crack tip change during the loading process is faster than the change during the unloading process; the stress intensity factor during the unloading process under the same load condition is larger than that during the loading process; under the same load, the theoretical plastic zone during the unloading process is higher than that during the loading process.
Bing Yang
,
Zhanjiang Wei
,
Zhen Liao
,
Shuwei Zhou
,
Shoune Xiao
,
Tao Zhu
,
Guangwu Yang
,
Mingmeng Wang
. Optimisation Method for Determination of Crack Tip Position Based on Gauss-Newton Iterative Technique[J]. Chinese Journal of Mechanical Engineering, 2021
, 34(4)
: 70
-70
.
DOI: 10.1186/s10033-021-00585-0
In the digital image correlation research of fatigue crack growth rate, the accuracy of the crack tip position determines the accuracy of the calculation of the stress intensity factor, thereby affecting the life prediction. This paper proposes a Gauss-Newton iteration method for solving the crack tip position. The conventional linear fitting method provides an iterative initial solution for this method, and the preconditioned conjugate gradient method is used to solve the ill-conditioned matrix. A noise-added artificial displacement field is used to verify the feasibility of the method, which shows that all parameters can be solved with satisfactory results. The actual stress intensity factor solution case shows that the stress intensity factor value obtained by the method in this paper is very close to the finite element result, and the relative error between the two is only ??0.621%; The Williams coefficient obtained by this method can also better define the contour of the plastic zone at the crack tip, and the maximum relative error with the test plastic zone area is ??11.29%. The relative error between the contour of the plastic zone defined by the conventional method and the area of the experimental plastic zone reached a maximum of 26.05%. The crack tip coordinates, stress intensity factors, and plastic zone contour changes in the loading and unloading phases are explored. The results show that the crack tip change during the loading process is faster than the change during the unloading process; the stress intensity factor during the unloading process under the same load condition is larger than that during the loading process; under the same load, the theoretical plastic zone during the unloading process is higher than that during the loading process.
[1] D R J Owen, A J Fawkes. Engineering fracture mechanics: Numerical methods and applications. Swansea: Pineridge Press, 1983.
[2] P Paul, K Murthy, D Chakraborty. A strain gage technique for mode I notch stress intensity factor of sharp V-notched configurations. Theoretical and Applied Fracture Mechanics, 2018, 94: 57-70.
[3] D Chakraborty, K Murthy. A strain gage technique for the determination of mixed mode stress intensity factors of orthotropic materials. Composite Structures, 2017, 160: 185-194.
[4] L V Stepanova, O N Belova. The digital photoelasticity method and finite element analysis for determination of the multi-point crack-tip field series expansions for notched semi-disks. Journal of Physics: Conference Series, 2021, 1745(1): 1-10.
[5] V V Avtaev, N O Yakovlev, V S Erasov, et al. Determination of the crack tip position by digital image correlation during static and cyclic fracture toughness tests of aluminum alloys. Russian Metallurgy (Metally), 2019, https://doi.org/10.1134/S0036029519040025.
[6] V V Avtaev, N O Yakovlev. Digital image correlation method study of the static fracture toughness and the fracture strength of a thin-sheet aluminum alloy. Russian Metallurgy (Metally), 2020. https://doi.org/10.1134/S0036029520100043.
[7] B Yang, M N James. Fatigue crack growth rate curve based on the CJP model and its application method. Journal of Mechanical Engineering, 2018, 54(18): 76-84. (in Chinese)
[8] J M Vasco-Olmo, F A Díaz, F V Antunes, et al. Experimental evaluation of CTOD in constant amplitude fatigue crack growth from crack tip displacement fields. Frattura ed Integrità Strutturale, 2017, 11(41): 157-165.
[9] J M Vasco-Olmo, F A Díaz, F V Antunes, et al. Characterisation of fatigue crack growth using digital image correlation measurements of plastic CTOD. Theoretical and Applied Fracture Mechanics, 2019, 101: 332-341.
[10] J M Vasco-Olmo, F Garrido, F V Antunes, et al. Plastic CTOD as fatigue crack growth characterising parameter in 2024-T3 and 7050-T6 aluminium alloys using DIC. Fatigue & Fracture of Engineering Materials & Structures, 2020, 43(8): 1-12.
[11] H Gao, Z Zhang, W Jiang, et al. Deformation fields measurement of crack tip under high-frequency resonant loading using a novel hybrid image processing method. Shock & Vibration, 2018. https://doi.org/10.1155/2018/1928926.
[12] M Mokhtarishirazabad, P Lopez-Crespo, B Moreno, et al. Evaluation of crack-tip fields from DIC data: A parametric study. International Journal of Fatigue, 2016, 89: 11-19.
[13] Y Cao, Q Ma, W Guo. Digital image correlation based high-speed crack tip locating method and its application. Journal of Beijing Institute of Technology, 2016, 25: 7-16.
[14] Z Yue, Y Song, P Li, et al. Applications of digital image correlation (DIC) and the strain gage method for measuring dynamic mode I fracture parameters of the white marble specimen. Rock Mechanics and Rock Engineering, 2019, 52(11): 4203-4216.
[15] J Tong. Full-field characterisation of crack tip deformation and fatigue crack growth using digital image correlation-A review. Fatigue & Fracture of Engineering Materials & Structures, 2018, 41(9): 1855-1869.
[16] M Zhu, Y W Lu, C Lupton, et al. In situ near-tip normal strain evolution of a growing fatigue crack. Fatigue & Fracture of Engineering Materials & Structures, 2016, 39(8): 950-955.
[17] S Yoneyama, T Ogawa, Y Kobayashi. Evaluating mixed-mode stress intensity factors from full-field displacement fields obtained by optical methods. Engineering Fracture Mechanics, 2006, 74(9): 1399-1412.
[18] M Zanganeh, P Lopez-Crespo, Y H Tai, et al. Locating the crack tip using displacement field data: A comparative study. Strain, 2013, 49(2): 102-115.
[19] W Zhang, S Li, D Dong, et, al. Locating the crack tip algorithm based on displacement field. Journal of Aerospace Power, 2018, 9(33): 2198-2204. (in Chinese)
[20] M S Fernández. Uncertainty evaluation of crack-tip parameters from displacement field data. Engineering Fracture Mechanics, 2018, 205: 333-346.
[21] Williams. On the stress distribution at the base of a stationary crack. Journal of Applied Mechanics, 1957, 24(1): 109-114.
[22] T Fan. Foundation of fracture theory. Beijing: Science Press, 2003. (in Chinese)
[23] M Mokhtarishirazabad, P Lopez-Crespo, M Zanganeh. Stress intensity factor monitoring under cyclic loading by digital image correlation. Fatigue & Fracture of Engineering Materials & Structures, 2018, 41(10): 2162-2171.
[24] L Parnas, G B ?mer, E Tezcan. Strain gage methods for measurement of opening mode stress intensity factor. Engineering Fracture Mechanics, 1996, 55(3): 485-492.
[25] M Y Dehnavi, S Khaleghian, A Emami, et al. Utilizing digital image correlation to determine stress intensity factors. Polymer Testing, 2014, 37: 28-35.
[26] B Pan, H Xie, Z Guo, et al. Full-field strain measurement using a two-dimensional Savitzky-Golay digital differentiator in digital image correlation. Optical Engineering, 2007, 46(3): 033601.
[27] J M Vasco-Olmo, F A Díaz, M N James, et al. Crack tip plastic zone evolution during an overload cycle and the contribution of plasticity-induced shielding to crack growth rate changes. Fatigue & Fracture of Engineering Materials & Structures, 2018, 41(10): 2172-2186.
[28] D Wang. Nonline equations solution and optimal method. Beijing: People's Education Press, 1979. (in Chinese)
[29] P J Lagace. Power flow methods for improving convergence. Conference of the IEEE Industrial Electronics Society, 2012: 1387-1392.
[30] S H Ju, S H Liu, K W Liu. Measurement of stress intensity factors by digital camera. International Journal of Solids & Structures, 2006, 43(5): 1009-1022.
[31] X Jia, X Hu, Y Song. Calculation formula for stress intensity factors of CT specimens based on three dimensional finite element solutions. Mechanical Engineering Materials, 2015, 9(12): 30-34. (in Chinese)
