Machine-Precise Evaluation of Stress Intensity Factors with the Consistent Boundary Element Method

Machine-Precise Evaluation of Stress Intensity Factors with the Consistent Boundary Element Method

Ney Augusto Dumont Osmar Alexandre Do Amaral Neto

Department of Civil and Environmental Engineering, Pontifical Catholic University of Rio de Janeiro, Brazil

Available online: 
| Citation



As classically proposed in the technical literature, the boundary element modeling of cracks is best carried out by resorting to a hypersingular fundamental solution – in the frame of the so-called dual formulation – since with the singular fundamental solution alone, the ensuing topological issues would not be adequately tackled. A more natural approach might rely on the direct representation of the crack tip singularity, as already proposed in the frame of the hybrid boundary element method, with implementation of generalized Westergaard stress functions. On the other hand, recent mathematical assessments indicate that the conventional boundary element formulation – based on Kelvin’s fundamental solution – is, in fact, able to precisely represent high stress gradients and deal with extremely convoluted topologies provided only that the numerical integrations be properly resolved. We propose in this paper that inde- pendent of the configuration, a cracked structure is geometrically represented as it would appear in real-world laboratory experiments, with crack openings in the range of micrometers. (The nanometer range is actually mathematically feasible, but not realistic in terms of continuum mechanics.) Owing to the newly developed numerical integration scheme, machine precision evaluation of all quantities may be achieved and stress results consistently evaluated at interior points arbitrarily close to crack tips. Importantly, no artificial topological issues are introduced, linear algebra conditioning is kept well under control, and arbitrarily high convergence of results is always attainable. The present develop- ments apply to two-dimensional problems. Some numerical illustrations show that highly accurate results are obtained for cracks represented with just a few quadratic, generally curved, boundary ele- ments – and a few Gauss–Legendre integration points per element – and that the numerical evaluation of the J-integral turns out to be straightforward and actually the most reliable means of obtaining stress intensity factors. Higher-order boundary elements lead to still better results.


boundary elements, fracture mechanics, machine precision integration, stress intensity factor


[1] Inglis, C.E., Stress in a plate due to the presence of cracks and sharp corners. Transactions of the institute of Naval Archtects, 55, pp. 219–241, 1913.

[2] Griffth, A.A., The phenomena of rupture and flow in solids. Philosophical Transactions, 221, pp. 163–198, 1920.

[3] Irwin, G.R., Onset of fast crack propagation in high strength steel and aluminum alloys. Sagamore Research Conference Proceedings, 2, pp. 289–305, 1956.

[4] Westergaard, H.M., Bearing pressures and cracks. Journal of Applied Mechanics, 6, pp. 49–53, 1939.

[5] Irwin, G.R., Analysis of stresses and strains near the end of a crack traversing a plate. Journal of Applied Mechanics, 24, pp. 361–364, 1957.

[6] Dugdale, D.S., Yielding of steel sheets containing slit. Journal of Mechanical Physics Solids, 8(2), p. 100–104, 1960.

[7] Rice, J.R., A path independent integral and the approximate analysis of strain con- centration by notches and cracks. Journal of Applied Mechanics, 35(2), pp. 379–386, 1968.

[8] Jiangbo, S., Jun, S., Pin, Z., Zengjie, D. & Huijui, Z., Study of the relationship between J-integral and COD parameters under mixed mode I + II loading in aluminum alloy  Ly 12. International Journal of Fracture, 104(4), p. 409–423, 2000. https://doi. org/10.1023/a:1007647230722

[9] Cruse, T.A. & Van Buren, W., Three dimensional elastic stress analysis of a fracture specimen with an edge crack. International Journal for Fracture Mechanics, 7(1), pp. 1–15, 1971.

[10] Dumont, N.A. & Lopes, A.A.O., On the explicit evaluation of stress intensity factors in the hybrid boundary element method. Fatigue & Fracture of Engineering Materials & Structures, 26(2), pp. 151–165, 2002. 2695.2003.00591.x

[11] Dumont, N.A. & Mamani, E.Y., Generalized Westergaard stress functions as fundamen- tal solutions. CMES – Computer Modeling in Engineering & Sciences, 78, pp. 109–150, 2011. doi:10.3970/cmes.2011.078.109

[12] Dumont, N.A., Mamani, E.Y. & Cardoso, M.L., A boundary element implementation for fracture mechanics problems using generalized Westergaard stress functions. European Journal of Computational Mechanics, 27(5–6), pp. 401–424, 2018. doi:10.1080/ 17797179.2018.1499188

[13] Dumont, N.A., The boundary element method revisited. Boundary Elements and Other Mesh Reduction Methods XXXII, 50, pp. 227–238, 2010. doi:10.2495/BE100201

[14] Dumont, N.A., The collocation boundary element method revisited: Perfect code for 2D problems. International Journal of Computational Methods and Experimental Measurements, 6(6), pp. 965–975, 2018.

[15] Dumont, N.A. & Kurz, T.G., Analytical 3D boundary element implementation of flat triangle and quadrilateral elements for potential and linear elasticity problems. Boundary Elements and Other Mesh Reduction Methods 42, WIT Transactions on Engineering Sciences, volume 126, pp. 1–11, 2019. doi:10.2495/BE420011

[16] Amaral Neto, O.A., Consistent application of the boundary element method to fracture mechanics problems. Master’s thesis, PUC-Rio, Rio de Janeiro, Brazil, 2020.

[17] Amaral Neto, O.A. & Dumont, N.A., Machine-precision fracture mechanics evaluations with the consistent boundary element method. Proceedings of CILAMCE – XLI Iberian Latin-American Congress on Computational Methods in Engineering, pp. 1–07, 2020.

[18] Dumont, N.A., On the eficient numerical evaluation of integrals with complex singularity poles. Engineering Analysis with Boundary Elements, 13(2), pp. 155–168, 1994.

[19] Newman Jr., J.C., An improved method of collocation for the stress analysis of cracked plates with various shaped boundaries. Technical report, NASA technical note, D-6376, 1971.

[20] Namakian, R., Shodja, H.M. & Mashayekhi, M., Fully enriched weight functions in mesh-free methods for the analysis of linear elastic fracture mechanics problems. Engineering Analysis with Boundary Elements, 43, pp. 1–18, 2014. enganabound.2014.02.006

[21] Zhuang, X., Augarde, C. & Bordas, S., Accurate fracture modeling using meshless- methods, the visibility criterion and level sets: formulation and 2d modeling. International Journal for Numerical Methods in Engineering, 86(2), p. 249–268, 2011. https://

[22] Murakami, Y., Stress Intensity Factors Handbook. Pergamon Press, 1987.